face_props.hpp File Reference

flow — face/cell material-property accessors for the variable-coefficient momentum operator. More...

#include <Kokkos_Core.hpp>
#include "mac_cutcell.hpp"

Go to the source code of this file.

Classes
struct	peclet::flow::UniformFaceProps
struct	peclet::flow::FieldFaceProps
struct	peclet::flow::VarFaceProps

Namespaces
namespace	peclet
namespace	peclet::flow

Detailed Description

flow — face/cell material-property accessors for the variable-coefficient momentum operator.

The variable-viscosity diffusion stencil (ibmBuildDiffusionVar) is templated on a small accessor supplying the per-cell time diagonal idiag(i) = rho/dt and the per-face viscosity beta between two adjacent cells. Two models:

• UniformFaceProps: two constants — reproduces the constant-mu operator (used to cross-check the Var kernel against the scalar ibmBuildDiffusion at the solution level).
• FieldFaceProps: idiag from a (constant here) rho, beta from a per-cell mu field, averaged at the face either arithmetically (default) or harmonically (correct for a viscosity jump — the shear stress mu*du/dy is continuous across a material interface).

Averaging is done in double; the caller casts the assembled band to float once (matching the constant path's (float)(idiag + 6*beta)), so no per-cell float rounding creeps into the sum.

Definition in file face_props.hpp.