ERF_GetRhoAlphaForFaces.H File Reference

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void	getRhoAlphaForFaces (int i, int j, int k, int ioff, int joff, amrex::Real &rhoAlpha_lo, amrex::Real &rhoAlpha_hi, const amrex::Array4< const amrex::Real > &cell_data, const amrex::Array4< const amrex::Real > &mu_turb, const amrex::Real mu_eff, bool l_consA, bool l_turb, const int eddy_comp=EddyDiff::Mom_v)

Function Documentation

getRhoAlphaForFaces()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void getRhoAlphaForFaces	(	int	i,
		int	j,
		int	k,
		int	ioff,
		int	joff,
		amrex::Real &	rhoAlpha_lo,
		amrex::Real &	rhoAlpha_hi,
		const amrex::Array4< const amrex::Real > &	cell_data,
		const amrex::Array4< const amrex::Real > &	mu_turb,
		const amrex::Real	mu_eff,
		bool	l_consA,
		bool	l_turb,
		const int	eddy_comp = EddyDiff::Mom_v
	)

Average rhoAlpha from cell centers to lower and upper z-edges of staggered faces.

Parameters

[in]	i	face i-index
[in]	j	face j-index
[in]	k	face k-index
[in]	ioff	x-direction staggering offset
[in]	joff	y-direction staggering offset
[out]	rhoAlpha_lo	rhoAlpha averaged to the lower z-edge
[out]	rhoAlpha_hi	rhoAlpha averaged to the upper z-edge
[in]	cell_data	conserved cell-centered state
[in]	mu_turb	turbulent viscosity
[in]	mu_eff	effective molecular viscosity. NOTE: unlike the mu_eff in ERF_MakeTauTerms.cpp and ComputeStress*, this does not include the factor of 2, because the u/v correction strains carry a factor of 1/2.
[in]	l_consA	flag for conservative constant-alpha diffusion
[in]	l_turb	flag for turbulent diffusion
[in]	eddy_comp	eddy-viscosity component to average (EddyDiff::Mom_v, or EddyDiff::Mom_h for the terrain-metric coefficient)

{
    if (l_consA && l_turb) {
        rhoAlpha_lo = fourth * ( cell_data(i,j,k  ,Rho_comp       ) + cell_data(i-ioff,j-joff,k  ,Rho_comp       )
                               + cell_data(i,j,k-1,Rho_comp       ) + cell_data(i-ioff,j-joff,k-1,Rho_comp       ) ) * mu_eff
                    + fourth * (   mu_turb(i,j,k  ,eddy_comp) +   mu_turb(i-ioff,j-joff,k  ,eddy_comp)
                               +   mu_turb(i,j,k-1,eddy_comp) +   mu_turb(i-ioff,j-joff,k-1,eddy_comp) );
        rhoAlpha_hi = fourth * ( cell_data(i,j,k  ,Rho_comp       ) + cell_data(i-ioff,j-joff,k  ,Rho_comp       )
                               + cell_data(i,j,k+1,Rho_comp       ) + cell_data(i-ioff,j-joff,k+1,Rho_comp       ) ) * mu_eff
                    + fourth * (   mu_turb(i,j,k  ,eddy_comp) +   mu_turb(i-ioff,j-joff,k  ,eddy_comp)
                               +   mu_turb(i,j,k+1,eddy_comp) +   mu_turb(i-ioff,j-joff,k+1,eddy_comp) );
    }
    else if (l_turb) // with MolecDiffType::Constant or None
    {
        rhoAlpha_lo = mu_eff
                    + fourth * (   mu_turb(i,j,k  ,eddy_comp) +   mu_turb(i-ioff,j-joff,k  ,eddy_comp)
                               +   mu_turb(i,j,k-1,eddy_comp) +   mu_turb(i-ioff,j-joff,k-1,eddy_comp) );
        rhoAlpha_hi = mu_eff
                    + fourth * (   mu_turb(i,j,k  ,eddy_comp) +   mu_turb(i-ioff,j-joff,k  ,eddy_comp)
                               +   mu_turb(i,j,k+1,eddy_comp) +   mu_turb(i-ioff,j-joff,k+1,eddy_comp) );
    }
    else if (l_consA) // without an LES/PBL model
    {
        rhoAlpha_lo = fourth * ( cell_data(i,j,k  ,Rho_comp       ) + cell_data(i-ioff,j-joff,k  ,Rho_comp       )
                               + cell_data(i,j,k-1,Rho_comp       ) + cell_data(i-ioff,j-joff,k-1,Rho_comp       ) ) * mu_eff;
        rhoAlpha_hi = fourth * ( cell_data(i,j,k  ,Rho_comp       ) + cell_data(i-ioff,j-joff,k  ,Rho_comp       )
                               + cell_data(i,j,k+1,Rho_comp       ) + cell_data(i-ioff,j-joff,k+1,Rho_comp       ) ) * mu_eff;
    }
    else // with MolecDiffType::Constant or None - without an LES/PBL model
    {
        rhoAlpha_lo = mu_eff;
        rhoAlpha_hi = mu_eff;
    }
}

Referenced by ImplicitDiffForMomLU_N(), ImplicitDiffForMomLU_S(), and ImplicitDiffForMomLU_T().

Here is the caller graph for this function: