ERF_ComputeStress_N.cpp File Reference

#include <ERF_Diffusion.H>

Include dependency graph for ERF_ComputeStress_N.cpp:

Functions
void	ComputeStressConsVisc_N (Box bxcc, Box tbxxy, Box tbxxz, Box tbxyz, Real mu_eff, const Array4< const Real > &cell_data, Array4< Real > &tau11, Array4< Real > &tau22, Array4< Real > &tau33, Array4< Real > &tau12, Array4< Real > &tau13, Array4< Real > &tau23, const Array4< const Real > &er_arr, Array4< Real > &tau13i, Array4< Real > &tau23i)
void	ComputeStressVarVisc_N (Box bxcc, Box tbxxy, Box tbxxz, Box tbxyz, Real mu_eff, const Array4< const Real > &mu_turb, const Array4< const Real > &cell_data, Array4< Real > &tau11, Array4< Real > &tau22, Array4< Real > &tau33, Array4< Real > &tau12, Array4< Real > &tau13, Array4< Real > &tau23, const Array4< const Real > &er_arr, Array4< Real > &tau13i, Array4< Real > &tau23i)

Function Documentation

ComputeStressConsVisc_N()

void ComputeStressConsVisc_N	(	Box	bxcc,
		Box	tbxxy,
		Box	tbxxz,
		Box	tbxyz,
		Real	mu_eff,
		const Array4< const Real > &	cell_data,
		Array4< Real > &	tau11,
		Array4< Real > &	tau22,
		Array4< Real > &	tau33,
		Array4< Real > &	tau12,
		Array4< Real > &	tau13,
		Array4< Real > &	tau23,
		const Array4< const Real > &	er_arr,
		Array4< Real > &	tau13i,
		Array4< Real > &	tau23i
	)

Function for computing the stress with constant viscosity and without terrain.

Parameters

[in]	bxcc	cell center box for tau_ii
[in]	tbxxy	nodal xy box for tau_12
[in]	tbxxz	nodal xz box for tau_13
[in]	tbxyz	nodal yz box for tau_23
[in]	mu_eff	constant molecular viscosity
[in]	cell_data	to access rho if ConstantAlpha
[in,out]	tau11	11 strain -> stress
[in,out]	tau22	22 strain -> stress
[in,out]	tau33	33 strain -> stress
[in,out]	tau12	12 strain -> stress
[in,out]	tau13	13 strain -> stress
[in,out]	tau23	23 strain -> stress
[in]	er_arr	expansion rate
[in,out]	tau13i	contribution to stress from du/dz
[in,out]	tau23i	contribution to stress from dv/dz

{
    Real OneThird   = (one/three);
 
    // NOTE: mu_eff includes factor of 2
 
    if (cell_data)
    // constant alpha (stored in mu_eff)
    {
        // Cell centered strains
        ParallelFor(bxcc, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
        {
            Real rhoAlpha  = cell_data(i, j, k, Rho_comp) * mu_eff;
            tau11(i,j,k) = -rhoAlpha * ( tau11(i,j,k) - OneThird*er_arr(i,j,k) );
            tau22(i,j,k) = -rhoAlpha * ( tau22(i,j,k) - OneThird*er_arr(i,j,k) );
            tau33(i,j,k) = -rhoAlpha * ( tau33(i,j,k) - OneThird*er_arr(i,j,k) );
        });
 
        // Off-diagonal strains
        ParallelFor(tbxxy,tbxxz,tbxyz,
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real rho_bar = fourth*( cell_data(i-1, j  , k, Rho_comp) + cell_data(i, j  , k, Rho_comp)
                                + cell_data(i-1, j-1, k, Rho_comp) + cell_data(i, j-1, k, Rho_comp) );
            tau12(i,j,k) *= -rho_bar * mu_eff;
        },
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real rho_bar = fourth*( cell_data(i-1, j, k  , Rho_comp) + cell_data(i, j, k  , Rho_comp)
                                + cell_data(i-1, j, k-1, Rho_comp) + cell_data(i, j, k-1, Rho_comp) );
            tau13(i,j,k) *= -rho_bar * mu_eff;
 
            if (tau13i) tau13i(i,j,k) *= -rho_bar * mu_eff;
        },
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real rho_bar = fourth*( cell_data(i, j-1, k  , Rho_comp) + cell_data(i, j, k  , Rho_comp)
                                + cell_data(i, j-1, k-1, Rho_comp) + cell_data(i, j, k-1, Rho_comp) );
            tau23(i,j,k) *= -rho_bar * mu_eff;
 
            if (tau23i) tau23i(i,j,k) *= -rho_bar * mu_eff;
        });
    }
    else
    // constant mu_eff
    {
        // Cell centered strains
        ParallelFor(bxcc, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
        {
            tau11(i,j,k) = -mu_eff * ( tau11(i,j,k) - OneThird*er_arr(i,j,k) );
            tau22(i,j,k) = -mu_eff * ( tau22(i,j,k) - OneThird*er_arr(i,j,k) );
            tau33(i,j,k) = -mu_eff * ( tau33(i,j,k) - OneThird*er_arr(i,j,k) );
        });
 
        // Off-diagonal strains
        ParallelFor(tbxxy,tbxxz,tbxyz,
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            tau12(i,j,k) *= -mu_eff;
        },
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            tau13(i,j,k) *= -mu_eff;
 
            if (tau13i) tau13i(i,j,k) *= -mu_eff;
        },
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            tau23(i,j,k) *= -mu_eff;
 
            if (tau23i) tau23i(i,j,k) *= -mu_eff;
        });
    }
}

Referenced by erf_make_tau_terms().

Here is the call graph for this function:

Here is the caller graph for this function:

ComputeStressVarVisc_N()

void ComputeStressVarVisc_N	(	Box	bxcc,
		Box	tbxxy,
		Box	tbxxz,
		Box	tbxyz,
		Real	mu_eff,
		const Array4< const Real > &	mu_turb,
		const Array4< const Real > &	cell_data,
		Array4< Real > &	tau11,
		Array4< Real > &	tau22,
		Array4< Real > &	tau33,
		Array4< Real > &	tau12,
		Array4< Real > &	tau13,
		Array4< Real > &	tau23,
		const Array4< const Real > &	er_arr,
		Array4< Real > &	tau13i,
		Array4< Real > &	tau23i
	)

Function for computing the stress with variable viscosity and without terrain.

Parameters

[in]	bxcc	cell center box for tau_ii
[in]	tbxxy	nodal xy box for tau_12
[in]	tbxxz	nodal xz box for tau_13
[in]	tbxyz	nodal yz box for tau_23
[in]	mu_eff	constant molecular viscosity
[in]	mu_turb	variable turbulent viscosity
[in]	cell_data	to access rho if ConstantAlpha
[in,out]	tau11	11 strain -> stress
[in,out]	tau22	22 strain -> stress
[in,out]	tau33	33 strain -> stress
[in,out]	tau12	12 strain -> stress
[in,out]	tau13	13 strain -> stress
[in,out]	tau23	23 strain -> stress
[in]	er_arr	expansion rate
[in,out]	tau13i	contribution to stress from du/dz
[in,out]	tau23i	contribution to stress from dv/dz

{
    Real OneThird   = (one/three);
 
    // NOTE: mu_eff includes factor of 2
 
    if (cell_data)
    // constant alpha (stored in mu_eff)
    {
        // Cell centered strains
        ParallelFor(bxcc, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real rhoAlpha = cell_data(i, j, k, Rho_comp) * mu_eff;
            Real mu_11 = rhoAlpha + two * mu_turb(i, j, k, EddyDiff::Mom_h);
            Real mu_22 = mu_11;
            Real mu_33 = rhoAlpha + two * mu_turb(i, j, k, EddyDiff::Mom_v);
            tau11(i,j,k) = -mu_11 * ( tau11(i,j,k) - OneThird*er_arr(i,j,k) );
            tau22(i,j,k) = -mu_22 * ( tau22(i,j,k) - OneThird*er_arr(i,j,k) );
            tau33(i,j,k) = -mu_33 * ( tau33(i,j,k) - OneThird*er_arr(i,j,k) );
        });
 
        // Off-diagonal strains
        ParallelFor(tbxxy,tbxxz,tbxyz,
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real rho_bar = fourth*( cell_data(i-1, j  , k, Rho_comp) + cell_data(i, j  , k, Rho_comp)
                                + cell_data(i-1, j-1, k, Rho_comp) + cell_data(i, j-1, k, Rho_comp) );
            Real mu_bar = fourth*( mu_turb(i-1, j  , k, EddyDiff::Mom_h) + mu_turb(i, j  , k, EddyDiff::Mom_h)
                               + mu_turb(i-1, j-1, k, EddyDiff::Mom_h) + mu_turb(i, j-1, k, EddyDiff::Mom_h) );
            Real mu_12  = rho_bar*mu_eff + two*mu_bar;
            tau12(i,j,k) *= -mu_12;
        },
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real rho_bar = fourth*( cell_data(i-1, j, k  , Rho_comp) + cell_data(i, j, k  , Rho_comp)
                                + cell_data(i-1, j, k-1, Rho_comp) + cell_data(i, j, k-1, Rho_comp) );
            Real mu_bar = fourth*( mu_turb(i-1, j, k  , EddyDiff::Mom_v) + mu_turb(i, j, k  , EddyDiff::Mom_v)
                               + mu_turb(i-1, j, k-1, EddyDiff::Mom_v) + mu_turb(i, j, k-1, EddyDiff::Mom_v) );
            Real mu_13  = rho_bar*mu_eff + two*mu_bar;
            tau13(i,j,k) *= -mu_13;
 
            if (tau13i) tau13i(i,j,k) *= -mu_13;
        },
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real rho_bar = fourth*( cell_data(i, j-1, k  , Rho_comp) + cell_data(i, j, k  , Rho_comp)
                                + cell_data(i, j-1, k-1, Rho_comp) + cell_data(i, j, k-1, Rho_comp) );
            Real mu_bar = fourth*( mu_turb(i, j-1, k  , EddyDiff::Mom_v) + mu_turb(i, j, k  , EddyDiff::Mom_v)
                               + mu_turb(i, j-1, k-1, EddyDiff::Mom_v) + mu_turb(i, j, k-1, EddyDiff::Mom_v) );
            Real mu_23  = rho_bar*mu_eff + two*mu_bar;
            tau23(i,j,k) *= -mu_23;
 
            if (tau23i) tau23i(i,j,k) *= -mu_23;
        });
    }
    else
    // constant mu_eff
    {
        // Cell centered strains
        ParallelFor(bxcc, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real mu_11 = mu_eff + two * mu_turb(i, j, k, EddyDiff::Mom_h);
            Real mu_22 = mu_11;
            Real mu_33 = mu_eff + two * mu_turb(i, j, k, EddyDiff::Mom_v);
            tau11(i,j,k) = -mu_11 * ( tau11(i,j,k) - OneThird*er_arr(i,j,k) );
            tau22(i,j,k) = -mu_22 * ( tau22(i,j,k) - OneThird*er_arr(i,j,k) );
            tau33(i,j,k) = -mu_33 * ( tau33(i,j,k) - OneThird*er_arr(i,j,k) );
        });
 
        // Off-diagonal strains
        ParallelFor(tbxxy,tbxxz,tbxyz,
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real mu_bar = fourth*( mu_turb(i-1, j  , k, EddyDiff::Mom_h) + mu_turb(i, j  , k, EddyDiff::Mom_h)
                               + mu_turb(i-1, j-1, k, EddyDiff::Mom_h) + mu_turb(i, j-1, k, EddyDiff::Mom_h) );
            Real mu_12  = mu_eff + two*mu_bar;
            tau12(i,j,k) *= -mu_12;
        },
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real mu_bar = fourth*( mu_turb(i-1, j, k  , EddyDiff::Mom_v) + mu_turb(i, j, k  , EddyDiff::Mom_v)
                               + mu_turb(i-1, j, k-1, EddyDiff::Mom_v) + mu_turb(i, j, k-1, EddyDiff::Mom_v) );
            Real mu_13  = mu_eff + two*mu_bar;
            tau13(i,j,k) *= -mu_13;
 
            if (tau13i) tau13i(i,j,k) *= -mu_13;
        },
        [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept {
            Real mu_bar = fourth*( mu_turb(i, j-1, k  , EddyDiff::Mom_v) + mu_turb(i, j, k  , EddyDiff::Mom_v)
                               + mu_turb(i, j-1, k-1, EddyDiff::Mom_v) + mu_turb(i, j, k-1, EddyDiff::Mom_v) );
            Real mu_23  = mu_eff + two*mu_bar;
            tau23(i,j,k) *= -mu_23;
 
            if (tau23i) tau23i(i,j,k) *= -mu_23;
        });
    }
}

Referenced by erf_make_tau_terms().

Here is the call graph for this function:

Here is the caller graph for this function: