ERF_AdvectionSrcForMom_EB.cpp File Reference

#include "AMReX_BCRec.H"
#include <ERF_Advection.H>
#include <ERF_AdvectionSrcForMom_N.H>
#include <ERF_AdvectionSrcForMom_T.H>

Include dependency graph for ERF_AdvectionSrcForMom_EB.cpp:

Functions
void	AdvectionSrcForMom_EB (const Box &bxx, const Box &bxy, const Box &bxz, const Array4< Real > &rho_u_rhs, const Array4< Real > &rho_v_rhs, const Array4< Real > &rho_w_rhs, const Array4< const Real > &, const Array4< const Real > &, const Array4< const Real > &, const Array4< const Real > &, const Array4< const Real > &, const Array4< const Real > &, const Array4< const Real > &, const Array4< const Real > &, const Array4< const Real > &, const Array4< const Real > &, const GpuArray< Real, AMREX_SPACEDIM > &, const Array4< const Real > &, const Array4< const Real > &, const Array4< const Real > &, const AdvType, const AdvType, const Real, const Real, const int, const int, const Box &)

Function Documentation

AdvectionSrcForMom_EB()

void AdvectionSrcForMom_EB	(	const Box &	bxx,
		const Box &	bxy,
		const Box &	bxz,
		const Array4< Real > &	rho_u_rhs,
		const Array4< Real > &	rho_v_rhs,
		const Array4< Real > &	rho_w_rhs,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const GpuArray< Real, AMREX_SPACEDIM > &	,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const Array4< const Real > &	,
		const	AdvType,
		const	AdvType,
		const	Real,
		const	Real,
		const int	,
		const int	,
		const Box &
	)

Function for computing the advective tendency for the momentum equations when using EB

Parameters

[in]	bxx	box over which the x-momentum is updated
[in]	bxy	box over which the y-momentum is updated
[in]	bxz	box over which the z-momentum is updated
[out]	rho_u_rhs	tendency for the x-momentum equation
[out]	rho_v_rhs	tendency for the y-momentum equation
[out]	rho_w_rhs	tendency for the z-momentum equation
[in]	u	x-component of the velocity
[in]	v	y-component of the velocity
[in]	w	z-component of the velocity
[in]	rho_u	x-component of the momentum
[in]	rho_v	y-component of the momentum
[in]	Omega	component of the momentum normal to the z-coordinate surface
[in]	z_nd	height coordinate at nodes
[in]	ax	Area fraction of x-faces
[in]	ay	Area fraction of y-faces
[in]	az	Area fraction of z-faces
[in]	vf	Volume fraction
[in]	cellSizeInv	inverse of the mesh spacing
[in]	mf_m	map factor at cell centers
[in]	mf_u	map factor at x-faces
[in]	mf_v	map factor at y-faces
[in]	horiz_adv_type	sets the spatial order to be used for lateral derivatives
[in]	vert_adv_type	sets the spatial order to be used for vertical derivatives

{
    BL_PROFILE_VAR("AdvectionSrcForMom_EB", AdvectionSrcForMom_EB);
 
//     auto dxInv = cellSizeInv[0], dyInv = cellSizeInv[1], dzInv = cellSizeInv[2];
 
    AMREX_ALWAYS_ASSERT(bxz.smallEnd(2) > 0);
 
#if 0
    // compute mapfactor inverses
    Box box2d_u(bxx);   box2d_u.setRange(2,0);   box2d_u.grow({3,3,0});
    Box box2d_v(bxy);   box2d_v.setRange(2,0);   box2d_v.grow({3,3,0});
 
    FArrayBox mf_u_invFAB(box2d_u,1,The_Async_Arena());
    FArrayBox mf_v_invFAB(box2d_v,1,The_Async_Arena());
 
    const Array4<Real>& mf_u_inv = mf_u_invFAB.array();
    const Array4<Real>& mf_v_inv = mf_v_invFAB.array();
 
    ParallelFor(box2d_u, box2d_v,
    [=] AMREX_GPU_DEVICE (int i, int j, int) noexcept
    {
        mf_u_inv(i,j,0) = 1. / mf_u(i,j,0);
    },
    [=] AMREX_GPU_DEVICE (int i, int j, int) noexcept
    {
        mf_v_inv(i,j,0) = 1. / mf_v(i,j,0);
    });
#endif
 
    ParallelFor(bxx, bxy, bxz,
    [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
    {
        rho_u_rhs(i, j, k) = 0.0;
    },
    [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
    {
        rho_v_rhs(i, j, k) = 0.0;
    },
    [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
    {
        rho_w_rhs(i, j, k) = 0.0;
    });
}