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)

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 =
    (qty_comp >= RhoScalar_comp && qty_comp < RhoScalar_comp + NSCALARS)
      ? BCVars::RhoScalar_bc_comp
      : qty_comp;
  if (policy.bc_comp > BCVars::RhoScalar_bc_comp) {
    policy.bc_comp -= (NSCALARS - 1);
  }
  policy.is_theta = qty_comp == RhoTheta_comp;
  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 || policy.is_q1;
 
  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 only through Q6.
    policy.coefficients.eddy_h_comp = EddyDiff::Q_h;
    policy.coefficients.eddy_v_comp = EddyDiff::Q_v;
    if (qty_comp >= RhoQ1_comp && qty_comp <= RhoQ6_comp) {
      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(), and ImplicitDiffForStateLU_T().

Here is the caller graph for this function: