ERF_NativeScalarDiffusion.H File Reference

#include "ERF_ScalarDiffusion.H"
#include "ERF_IndexDefines.H"

Include dependency graph for ERF_NativeScalarDiffusion.H:

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

Go to the source code of this file.

Classes
struct	NativeScalarDiffusionPolicy

Functions
NativeScalarDiffusionPolicy	ResolveNativeScalarDiffusionPolicy (const int qty_comp, const DiffChoice &choice)
amrex::Vector< amrex::Real >	VerticalMolecularAlphaByPrim (const DiffChoice &choice)

Function Documentation

ResolveNativeScalarDiffusionPolicy()

NativeScalarDiffusionPolicy ResolveNativeScalarDiffusionPolicy ( const int  qty_comp, const DiffChoice &  choice  )

inline

Resolve conserved-state identity to ERF's native coefficient category.

{
    NativeScalarDiffusionPolicy policy;
    policy.scalar_comp = qty_comp - 1;
    policy.rhs_comp    = qty_comp;
    policy.flux_comp   = 0;
    policy.bc_comp     = cons_bc_comp(qty_comp);
    policy.is_theta = (qty_comp == RhoTheta_comp);
    // Water identity determines molecular diffusion, independently of ownership
    // and eligibility for the implicit turbulent vertical solve.
    policy.is_moist = (qty_comp >= RhoQ1_comp && qty_comp < (RhoQ1_comp + choice.n_qstate_moist));
    policy.is_q1    = (qty_comp == RhoQ1_comp);
    policy.is_q2    = (qty_comp == RhoQ2_comp);
    policy.scale_raw_vertical_flux =
        (policy.is_theta || qty_comp == RhoKE_comp || choice.implicit_q(qty_comp));
 
    if (policy.is_theta) {
        policy.coefficients.eddy_h_comp = EddyDiff::Theta_h;
        policy.coefficients.eddy_v_comp = EddyDiff::Theta_v;
        if (choice.molec_diff_type == MolecDiffType::ConstantAlpha) {
            policy.coefficients.molecular_coeff = choice.alpha_T;
        } else if (choice.molec_diff_type == MolecDiffType::Constant) {
            policy.coefficients.molecular_coeff = choice.rhoAlpha_T;
        }
    } else if (qty_comp == RhoKE_comp) {
        policy.coefficients.eddy_h_comp = EddyDiff::KE_h;
        policy.coefficients.eddy_v_comp = EddyDiff::KE_v;
    } else if (qty_comp >= RhoScalar_comp && qty_comp < RhoScalar_comp + NSCALARS) {
        policy.coefficients.eddy_h_comp = EddyDiff::Scalar_h;
        policy.coefficients.eddy_v_comp = EddyDiff::Scalar_v;
        if (choice.molec_diff_type == MolecDiffType::ConstantAlpha) {
            policy.coefficients.molecular_coeff = choice.alpha_C;
        } else if (choice.molec_diff_type == MolecDiffType::Constant) {
            policy.coefficients.molecular_coeff = choice.rhoAlpha_C;
        }
    } else {
        // Every q-state, including future appended species, uses the Q eddy
        // category. Molecular diffusion is defined for the moisture species
        // (number concentrations included), not for appended non-moist species.
        policy.coefficients.eddy_h_comp = EddyDiff::Q_h;
        policy.coefficients.eddy_v_comp = EddyDiff::Q_v;
        if (policy.is_moist) {
            if (choice.molec_diff_type == MolecDiffType::ConstantAlpha) {
                policy.coefficients.molecular_coeff = choice.alpha_C;
            } else if (choice.molec_diff_type == MolecDiffType::Constant) {
                policy.coefficients.molecular_coeff = choice.rhoAlpha_C;
            }
        }
    }
    return policy;
}

Referenced by DiffusionSrcForState_N(), DiffusionSrcForState_S(), DiffusionSrcForState_T(), if(), ImplicitDiffForStateLU_N(), ImplicitDiffForStateLU_S(), ImplicitDiffForStateLU_T(), and ImplicitVerticalFluxForReflux_T().

Here is the call graph for this function:

Here is the caller graph for this function:

VerticalMolecularAlphaByPrim()

amrex::Vector<amrex::Real> VerticalMolecularAlphaByPrim ( const DiffChoice &  choice )

inline

Molecular diffusivity by primitive component for the implicit vertical solve and the EB diffusion (ERF_SetupVertDiff.H): theta, every passive scalar (all NSCALARS of them) and the moisture species (number concentrations included); everything else – KE, and any non-moist species a scheme appends after its moisture species – is left at zero. With erf.molec_diff_type = Constant or ConstantAlpha these are the coefficients ResolveNativeScalarDiffusionPolicy() gives the explicit paths.

{
    amrex::Vector<amrex::Real> alpha_eff(NPRIMVAR_max, amrex::Real(0));
    const bool consA = (choice.molec_diff_type == MolecDiffType::ConstantAlpha);
    const amrex::Real alpha_T_eff = consA ? choice.alpha_T : choice.rhoAlpha_T;
    const amrex::Real alpha_C_eff = consA ? choice.alpha_C : choice.rhoAlpha_C;
    alpha_eff[PrimTheta_comp] = alpha_T_eff;
    for (int n = 0; n < NSCALARS; ++n) {
        alpha_eff[PrimScalar_comp + n] = alpha_C_eff;
    }
    for (int i = PrimQ1_comp; i < PrimQ1_comp + choice.n_qstate_moist; ++i) {
        alpha_eff[i] = alpha_C_eff;
    }
    return alpha_eff;
}