ERF_AdvanceFitch.cpp File Reference

#include <ERF_Fitch.H>
#include <ERF_IndexDefines.H>
#include <ERF_Constants.H>
#include <ERF_Interpolation_1D.H>

Include dependency graph for ERF_AdvanceFitch.cpp:

Functions
AMREX_FORCE_INLINE AMREX_GPU_DEVICE Real	compute_A (const Real z, const Real hub_height, const Real rotor_rad)
AMREX_FORCE_INLINE AMREX_GPU_DEVICE Real	compute_Aijk (const Real z_k, const Real z_kp1, const Real hub_height, const Real rotor_rad)

Function Documentation

compute_A()

AMREX_FORCE_INLINE AMREX_GPU_DEVICE Real compute_A	(	const Real	z,
		const Real	hub_height,
		const Real	rotor_rad
	)

{
 
    Real d  = std::min(std::fabs(z - hub_height), rotor_rad);
    Real theta = std::acos(d/rotor_rad);
    Real A_s = rotor_rad*rotor_rad*theta - d*std::pow(std::max(rotor_rad*rotor_rad - d*d,0.0), 0.5);
    Real A = PI*rotor_rad*rotor_rad/2.0 - A_s;
 
    return A;
}

Referenced by compute_Aijk().

Here is the caller graph for this function:

compute_Aijk()

AMREX_FORCE_INLINE AMREX_GPU_DEVICE Real compute_Aijk	(	const Real	z_k,
		const Real	z_kp1,
		const Real	hub_height,
		const Real	rotor_rad
	)

{
 
    Real A_k   = compute_A(z_k, hub_height, rotor_rad);
    Real A_kp1 = compute_A(z_kp1, hub_height, rotor_rad);
 
    Real check = (z_k - hub_height)*(z_kp1 - hub_height);
    Real A_ijk;
    if(check > 0){
        A_ijk = std::fabs(A_k -A_kp1);
    }
    else{
        A_ijk = A_k + A_kp1;
    }
 
    return A_ijk;
}

Referenced by Fitch::source_terms_cellcentered().

Here is the call graph for this function:

Here is the caller graph for this function: