erf_cloud_chamber_wall_stress Namespace Reference

Functions
bool	wall_stress_requires_most_state (const erf_wall_thermodynamics::Boundary &walls) noexcept
template<int DIR, bool HIGH>
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real	stored_tau_from_physical_traction (Real traction) noexcept
template<int DIR, bool HIGH>
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE GpuArray< Real, AMREX_SPACEDIM >	wall_traction_at_cell (int i, int j, int k, const Array4< const Real > &state, const Array4< const Real > &base_state, const Array4< const Real > &u, const Array4< const Real > &v, const Array4< const Real > &w, const erf_wall_thermodynamics::FaceWall &wall, const GpuArray< Real, AMREX_SPACEDIM > &dx_inv, const Real rdOcp, const bool cloudy) noexcept
template<int DIR, int TANGENT, bool HIGH>
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real	averaged_physical_traction_component (int i, int j, int k, const Box &domain, const Array4< const Real > &state, const Array4< const Real > &base_state, const Array4< const Real > &u, const Array4< const Real > &v, const Array4< const Real > &w, const erf_wall_thermodynamics::FaceWall &wall, const GpuArray< Real, AMREX_SPACEDIM > &dx_inv, const Real rdOcp, const bool cloudy) noexcept
template<int DIR, int TANGENT, bool HIGH>
AMREX_FORCE_INLINE void	apply_tangent_stress (const Box &bx, const Box &domain, const Array4< const Real > &state, const Array4< const Real > &base_state, const Array4< const Real > &u, const Array4< const Real > &v, const Array4< const Real > &w, const Array4< Real > &tau, const erf_wall_thermodynamics::FaceWall &wall, const GpuArray< Real, AMREX_SPACEDIM > &dx_inv, const Real rdOcp, const bool cloudy)
template<int DIR, bool HIGH>
AMREX_FORCE_INLINE void	apply_face (const Box &bx, const Box &domain, const Array4< const Real > &state, const Array4< const Real > &base_state, const Array4< const Real > &u, const Array4< const Real > &v, const Array4< const Real > &w, const Array4< Real > &tau12, const Array4< Real > &tau13, const Array4< Real > &tau23, const GpuArray< Real, AMREX_SPACEDIM > &dx_inv, const erf_wall_thermodynamics::Boundary &walls, const Real rdOcp, const bool cloudy)
void	apply (const Box &bx, const Box &domain, const Array4< const Real > &state, const Array4< const Real > &base_state, const Array4< const Real > &u, const Array4< const Real > &v, const Array4< const Real > &w, const Array4< Real > &tau12, const Array4< Real > &tau13, const Array4< Real > &tau23, const GpuArray< Real, AMREX_SPACEDIM > &dx_inv, const erf_wall_thermodynamics::Boundary &walls, const Real rdOcp, const bool cloudy)
void	apply (const Box &bx, const Box &domain, const Array4< const Real > &state, const Array4< const Real > &u, const Array4< const Real > &v, const Array4< const Real > &w, const Array4< Real > &tau12, const Array4< Real > &tau13, const Array4< Real > &tau23, const GpuArray< Real, AMREX_SPACEDIM > &dx_inv, const erf_wall_thermodynamics::Boundary &walls)
bool	has_momentum_wall (const erf_wall_thermodynamics::Boundary &walls) noexcept

Function Documentation

apply() [1/2]

void erf_cloud_chamber_wall_stress::apply ( const Box &  bx, const Box &  domain, const Array4< const Real > &  state, const Array4< const Real > &  base_state, const Array4< const Real > &  u, const Array4< const Real > &  v, const Array4< const Real > &  w, const Array4< Real > &  tau12, const Array4< Real > &  tau13, const Array4< Real > &  tau23, const GpuArray< Real, AMREX_SPACEDIM > &  dx_inv, const erf_wall_thermodynamics::Boundary &  walls, const Real  rdOcp, const bool  cloudy  )

inline

{
    AMREX_ALWAYS_ASSERT_WITH_MESSAGE(
        !wall_stress_requires_most_state(walls) ||
        (base_state.nComp() > BaseState::p0_comp &&
         state.nComp() > RhoTheta_comp &&
         (!cloudy || state.nComp() > RhoQ1_comp)),
        "Cloud Chamber MOST momentum stress requires valid thermodynamic and base-state arrays");
    apply_face<0,false>(bx, domain, state, base_state, u, v, w, tau12, tau13, tau23, dx_inv, walls, rdOcp, cloudy);
    apply_face<0,true >(bx, domain, state, base_state, u, v, w, tau12, tau13, tau23, dx_inv, walls, rdOcp, cloudy);
    apply_face<1,false>(bx, domain, state, base_state, u, v, w, tau12, tau13, tau23, dx_inv, walls, rdOcp, cloudy);
    apply_face<1,true >(bx, domain, state, base_state, u, v, w, tau12, tau13, tau23, dx_inv, walls, rdOcp, cloudy);
    apply_face<2,false>(bx, domain, state, base_state, u, v, w, tau12, tau13, tau23, dx_inv, walls, rdOcp, cloudy);
    apply_face<2,true >(bx, domain, state, base_state, u, v, w, tau12, tau13, tau23, dx_inv, walls, rdOcp, cloudy);
}

Referenced by apply(), and erf_slow_rhs_pre().

Here is the call graph for this function:

Here is the caller graph for this function:

apply() [2/2]

void erf_cloud_chamber_wall_stress::apply ( const Box &  bx, const Box &  domain, const Array4< const Real > &  state, const Array4< const Real > &  u, const Array4< const Real > &  v, const Array4< const Real > &  w, const Array4< Real > &  tau12, const Array4< Real > &  tau13, const Array4< Real > &  tau23, const GpuArray< Real, AMREX_SPACEDIM > &  dx_inv, const erf_wall_thermodynamics::Boundary &  walls  )

inline

{
    AMREX_ALWAYS_ASSERT_WITH_MESSAGE(
        !wall_stress_requires_most_state(walls),
        "Cloud Chamber MOST momentum requires the base-state-aware wall-stress overload");
    apply(bx, domain, state, Array4<const Real>{}, u, v, w,
          tau12, tau13, tau23, dx_inv, walls, Real(0.0), false);
}

Here is the call graph for this function:

apply_face()

template<int DIR, bool HIGH>

AMREX_FORCE_INLINE void erf_cloud_chamber_wall_stress::apply_face	(	const Box &	bx,
		const Box &	domain,
		const Array4< const Real > &	state,
		const Array4< const Real > &	base_state,
		const Array4< const Real > &	u,
		const Array4< const Real > &	v,
		const Array4< const Real > &	w,
		const Array4< Real > &	tau12,
		const Array4< Real > &	tau13,
		const Array4< Real > &	tau23,
		const GpuArray< Real, AMREX_SPACEDIM > &	dx_inv,
		const erf_wall_thermodynamics::Boundary &	walls,
		const Real	rdOcp,
		const bool	cloudy
	)

{
    if constexpr (HIGH) {
        if (bx.bigEnd(DIR) != domain.bigEnd(DIR)) { return; }
    } else {
        if (bx.smallEnd(DIR) != domain.smallEnd(DIR)) { return; }
    }
 
    const auto wall = walls[2*DIR + (HIGH ? 1 : 0)];
    if (!erf_cloud_chamber_wall_flux::requires_tangential_momentum(wall)) {
        return;
    }
 
    if constexpr (DIR == 0) {
        apply_tangent_stress<DIR,1,HIGH>(bx, domain, state, base_state, u, v, w,
                                         tau12, wall, dx_inv, rdOcp, cloudy);
        apply_tangent_stress<DIR,2,HIGH>(bx, domain, state, base_state, u, v, w,
                                         tau13, wall, dx_inv, rdOcp, cloudy);
    } else if constexpr (DIR == 1) {
        apply_tangent_stress<DIR,0,HIGH>(bx, domain, state, base_state, u, v, w,
                                         tau12, wall, dx_inv, rdOcp, cloudy);
        apply_tangent_stress<DIR,2,HIGH>(bx, domain, state, base_state, u, v, w,
                                         tau23, wall, dx_inv, rdOcp, cloudy);
    } else {
        apply_tangent_stress<DIR,0,HIGH>(bx, domain, state, base_state, u, v, w,
                                         tau13, wall, dx_inv, rdOcp, cloudy);
        apply_tangent_stress<DIR,1,HIGH>(bx, domain, state, base_state, u, v, w,
                                         tau23, wall, dx_inv, rdOcp, cloudy);
    }
}

Here is the call graph for this function:

apply_tangent_stress()

template<int DIR, int TANGENT, bool HIGH>

AMREX_FORCE_INLINE void erf_cloud_chamber_wall_stress::apply_tangent_stress	(	const Box &	bx,
		const Box &	domain,
		const Array4< const Real > &	state,
		const Array4< const Real > &	base_state,
		const Array4< const Real > &	u,
		const Array4< const Real > &	v,
		const Array4< const Real > &	w,
		const Array4< Real > &	tau,
		const erf_wall_thermodynamics::FaceWall &	wall,
		const GpuArray< Real, AMREX_SPACEDIM > &	dx_inv,
		const Real	rdOcp,
		const bool	cloudy
	)

{
    Box target = amrex::surroundingNodes(amrex::surroundingNodes(bx, DIR), TANGENT);
    const int face = HIGH ? domain.bigEnd(DIR) + 1 : domain.smallEnd(DIR);
    target.setSmall(DIR, face);
    target.setBig(DIR, face);
    target.setSmall(TANGENT,
                    amrex::max(target.smallEnd(TANGENT), domain.smallEnd(TANGENT)));
    target.setBig(TANGENT,
                  amrex::min(target.bigEnd(TANGENT), domain.bigEnd(TANGENT) + 1));
    if (target.isEmpty()) { return; }
 
    amrex::ParallelFor(target,
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
    {
        const Real traction = averaged_physical_traction_component<DIR,TANGENT,HIGH>(
            i, j, k, domain, state, base_state, u, v, w, wall, dx_inv,
            rdOcp, cloudy);
        tau(i,j,k) = stored_tau_from_physical_traction<DIR,HIGH>(traction);
    });
}

Here is the call graph for this function:

averaged_physical_traction_component()

template<int DIR, int TANGENT, bool HIGH>

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real erf_cloud_chamber_wall_stress::averaged_physical_traction_component ( int  i, int  j, int  k, const Box &  domain, const Array4< const Real > &  state, const Array4< const Real > &  base_state, const Array4< const Real > &  u, const Array4< const Real > &  v, const Array4< const Real > &  w, const erf_wall_thermodynamics::FaceWall &  wall, const GpuArray< Real, AMREX_SPACEDIM > &  dx_inv, const Real  rdOcp, const bool  cloudy  )

noexcept

{
    int cell_lo[AMREX_SPACEDIM] = {i, j, k};
    int cell_hi[AMREX_SPACEDIM] = {i, j, k};
    const int node = cell_lo[TANGENT];
    cell_lo[TANGENT] = amrex::max(domain.smallEnd(TANGENT), node - 1);
    cell_hi[TANGENT] = amrex::min(domain.bigEnd(TANGENT), node);
    cell_lo[DIR] = HIGH ? domain.bigEnd(DIR) : domain.smallEnd(DIR);
    cell_hi[DIR] = cell_lo[DIR];
 
    const auto traction_lo = wall_traction_at_cell<DIR,HIGH>(
        cell_lo[0], cell_lo[1], cell_lo[2], state, base_state, u, v, w, wall,
        dx_inv, rdOcp, cloudy);
    const auto traction_hi = wall_traction_at_cell<DIR,HIGH>(
        cell_hi[0], cell_hi[1], cell_hi[2], state, base_state, u, v, w, wall,
        dx_inv, rdOcp, cloudy);
    return Real(0.5) * (traction_lo[TANGENT] + traction_hi[TANGENT]);
}

has_momentum_wall()

bool erf_cloud_chamber_wall_stress::has_momentum_wall ( const erf_wall_thermodynamics::Boundary &  walls )

inlinenoexcept

{
    for (const auto& wall : walls) {
        if (erf_cloud_chamber_wall_flux::requires_tangential_momentum(wall)) {
            return true;
        }
    }
    return false;
}

Referenced by erf_slow_rhs_pre().

Here is the call graph for this function:

Here is the caller graph for this function:

stored_tau_from_physical_traction()

template<int DIR, bool HIGH>

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real erf_cloud_chamber_wall_stress::stored_tau_from_physical_traction ( Real  traction )

noexcept

{
    amrex::ignore_unused(DIR);
    return HIGH ? -traction : traction;
}

wall_stress_requires_most_state()

bool erf_cloud_chamber_wall_stress::wall_stress_requires_most_state ( const erf_wall_thermodynamics::Boundary &  walls )

inlinenoexcept

{
    for (const auto& wall : walls) {
        if (wall.momentum.model == erf_wall_thermodynamics::MomentumModel::BulkAero &&
            wall.momentum.provider == erf_wall_thermodynamics::CoefficientProvider::MOST) {
            return true;
        }
    }
    return false;
}

Referenced by apply().

Here is the caller graph for this function:

wall_traction_at_cell()

template<int DIR, bool HIGH>

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE GpuArray<Real, AMREX_SPACEDIM> erf_cloud_chamber_wall_stress::wall_traction_at_cell ( int  i, int  j, int  k, const Array4< const Real > &  state, const Array4< const Real > &  base_state, const Array4< const Real > &  u, const Array4< const Real > &  v, const Array4< const Real > &  w, const erf_wall_thermodynamics::FaceWall &  wall, const GpuArray< Real, AMREX_SPACEDIM > &  dx_inv, const Real  rdOcp, const bool  cloudy  )

noexcept

{
    const auto u_t = erf_cloud_chamber_wall_flux::
        relative_tangential_velocity_cell_centered(DIR, i, j, k, u, v, w, wall);
    const Real U_t = erf_cloud_chamber_wall_flux::tangential_magnitude(u_t);
    const Real rho = state(i,j,k,Rho_comp);
    const bool most_momentum =
        wall.momentum.model == erf_wall_thermodynamics::MomentumModel::BulkAero &&
        wall.momentum.provider == erf_wall_thermodynamics::CoefficientProvider::MOST;
    const erf_cloud_chamber_wall_flux::MomentumWallSample sample{
        rho, u_t, U_t, Real(0.5) / dx_inv[DIR],
        most_momentum ? state(i,j,k,RhoTheta_comp) / rho : Real(0.0),
        most_momentum && cloudy ? state(i,j,k,RhoQ1_comp) / rho : Real(0.0),
        most_momentum ? base_state(i,j,k,BaseState::p0_comp) : Real(0.0),
        most_momentum ? rdOcp : Real(0.0),
        DIR == 2 ? (HIGH ? -1 : 1) : 0};
    return erf_cloud_chamber_wall_flux::evaluate_momentum_traction(wall, sample).
        traction_on_fluid;
}

Here is the call graph for this function: