This is the base class all materials should use if you are trying to use the Navier-Stokes Kernels. More...

#include <NavierStokesMaterial.h>

Inheritance diagram for NavierStokesMaterial:

Public Member Functions
	NavierStokesMaterial (const InputParameters &parameters)

Protected Member Functions
virtual void	computeProperties ()
	Must be called after the child class computes dynamic_viscocity. More...

Protected Attributes
const unsigned int	_mesh_dimension

const VariableGradient &	_grad_u

const VariableGradient &	_grad_v

const VariableGradient &	_grad_w

MaterialProperty< RealTensorValue > &	_viscous_stress_tensor

MaterialProperty< Real > &	_thermal_conductivity

MaterialProperty< Real > &	_dynamic_viscosity

MaterialProperty< std::vector< RealTensorValue > > &	_calA

MaterialProperty< std::vector< RealTensorValue > > &	_calC

MaterialProperty< std::vector< std::vector< RealTensorValue > > > &	_calE

std::vector< const VariableGradient * >	_vel_grads

const VariableValue &	_u_vel

const VariableValue &	_v_vel

const VariableValue &	_w_vel

const VariableValue &	_temperature

const VariableValue &	_enthalpy

const VariableValue &	_rho

const VariableValue &	_rho_u

const VariableValue &	_rho_v

const VariableValue &	_rho_w

const VariableValue &	_rho_E

const VariableValue &	_drho_dt

const VariableValue &	_drhou_dt

const VariableValue &	_drhov_dt

const VariableValue &	_drhow_dt

const VariableValue &	_drhoE_dt

const VariableGradient &	_grad_rho

const VariableGradient &	_grad_rho_u

const VariableGradient &	_grad_rho_v

const VariableGradient &	_grad_rho_w

const VariableGradient &	_grad_rho_E

MaterialProperty< Real > &	_hsupg

MaterialProperty< Real > &	_tauc

MaterialProperty< Real > &	_taum

MaterialProperty< Real > &	_taue

MaterialProperty< std::vector< Real > > &	_strong_residuals

const IdealGasFluidProperties &	_fp

Private Member Functions
void	computeHSUPG (unsigned int qp)

void	computeTau (unsigned int qp)

void	computeStrongResiduals (unsigned int qp)

Detailed Description

This is the base class all materials should use if you are trying to use the Navier-Stokes Kernels.

Note that the derived class just needs to compute dynamic_viscocity then call this class's computeProperties() function.

Also make sure that the derived class's validParams function just adds to this class's validParams.

Finally, note that this Material isn't registered with the MaterialFactory. The reason is that by itself this material doesn't work! You must derive from this material and compute dynamic_viscocity!

Definition at line 36 of file NavierStokesMaterial.h.

Constructor & Destructor Documentation

◆ NavierStokesMaterial()

NavierStokesMaterial::NavierStokesMaterial ( const InputParameters & parameters )

Definition at line 49 of file NavierStokesMaterial.C.

   : Material(parameters),
     _mesh_dimension(_mesh.dimension()),
     _grad_u(coupledGradient(NS::velocity_x)),
     _grad_v(_mesh_dimension >= 2 ? coupledGradient(NS::velocity_y) : _grad_zero),
     _grad_w(_mesh_dimension == 3 ? coupledGradient(NS::velocity_z) : _grad_zero),
  
     _viscous_stress_tensor(declareProperty<RealTensorValue>("viscous_stress_tensor")),
     _thermal_conductivity(declareProperty<Real>("thermal_conductivity")),
  
     // Declared here but _not_ calculated here
     // (See e.g. derived class, bighorn/include/materials/FluidTC1.h)
     _dynamic_viscosity(declareProperty<Real>("dynamic_viscosity")),
  
     // The momentum components of the inviscid flux Jacobians.
     _calA(declareProperty<std::vector<RealTensorValue>>("calA")),
  
     // "Velocity column" matrices
     _calC(declareProperty<std::vector<RealTensorValue>>("calC")),
  
     // Energy equation inviscid flux matrices, "cal E_{kl}" in the notes.
     _calE(declareProperty<std::vector<std::vector<RealTensorValue>>>("calE")),
     _vel_grads({&_grad_u, &_grad_v, &_grad_w}),
  
     // Coupled solution values needed for computing SUPG stabilization terms
     _u_vel(coupledValue(NS::velocity_x)),
     _v_vel(_mesh.dimension() >= 2 ? coupledValue(NS::velocity_y) : _zero),
     _w_vel(_mesh.dimension() == 3 ? coupledValue(NS::velocity_z) : _zero),
  
     _temperature(coupledValue(NS::temperature)),
     _enthalpy(coupledValue(NS::enthalpy)),
  
     // Coupled solution values
     _rho(coupledValue(NS::density)),
     _rho_u(coupledValue(NS::momentum_x)),
     _rho_v(_mesh.dimension() >= 2 ? coupledValue(NS::momentum_y) : _zero),
     _rho_w(_mesh.dimension() == 3 ? coupledValue(NS::momentum_z) : _zero),
     _rho_E(coupledValue(NS::total_energy)),
  
     // Time derivative values
     _drho_dt(coupledDot(NS::density)),
     _drhou_dt(coupledDot(NS::momentum_x)),
     _drhov_dt(_mesh.dimension() >= 2 ? coupledDot(NS::momentum_y) : _zero),
     _drhow_dt(_mesh.dimension() == 3 ? coupledDot(NS::momentum_z) : _zero),
     _drhoE_dt(coupledDot(NS::total_energy)),
  
     // Gradients
     _grad_rho(coupledGradient(NS::density)),
     _grad_rho_u(coupledGradient(NS::momentum_x)),
     _grad_rho_v(_mesh.dimension() >= 2 ? coupledGradient(NS::momentum_y) : _grad_zero),
     _grad_rho_w(_mesh.dimension() == 3 ? coupledGradient(NS::momentum_z) : _grad_zero),
     _grad_rho_E(coupledGradient(NS::total_energy)),
  
     // Material properties for stabilization
     _hsupg(declareProperty<Real>("hsupg")),
     _tauc(declareProperty<Real>("tauc")),
     _taum(declareProperty<Real>("taum")),
     _taue(declareProperty<Real>("taue")),
     _strong_residuals(declareProperty<std::vector<Real>>("strong_residuals")),
     _fp(getUserObject<IdealGasFluidProperties>("fluid_properties"))
 {
 }

Member Function Documentation

◆ computeHSUPG()

void NavierStokesMaterial::computeHSUPG ( unsigned int qp )

private

Definition at line 178 of file NavierStokesMaterial.C.

Referenced by computeProperties().

◆ computeProperties()

void NavierStokesMaterial::computeProperties ( )

protectedvirtual

Must be called after the child class computes dynamic_viscocity.

Must be called after the child class computes dynamic_viscosity.

Reimplemented in Air.

Definition at line 116 of file NavierStokesMaterial.C.

 {
   for (unsigned int qp = 0; qp < _qrule->n_points(); ++qp)
   {
     /******* Viscous Stress Tensor *******/
     // Technically... this _is_ the transpose (since we are loading these by rows)
     // But it doesn't matter....
     RealTensorValue grad_outer_u(_grad_u[qp], _grad_v[qp], _grad_w[qp]);
  
     grad_outer_u += grad_outer_u.transpose();
  
     Real div_vel = 0.0;
     for (unsigned int i = 0; i < 3; ++i)
       div_vel += (*_vel_grads[i])[qp](i);
  
     // Add diagonal terms
     for (unsigned int i = 0; i < 3; ++i)
       grad_outer_u(i, i) -= 2.0 / 3.0 * div_vel;
  
     grad_outer_u *= _dynamic_viscosity[qp];
  
     _viscous_stress_tensor[qp] = grad_outer_u;
  
     // Tabulated values of thermal conductivity vs. Temperature for air (k increases slightly with
     // T):
     // T (K)    k (W/m-K)
     // 273      0.0243
     // 373      0.0314
     // 473      0.0386
     // 573      0.0454
     // 673      0.0515
  
     // Pr = (mu * cp) / k  ==>  k = (mu * cp) / Pr = (mu * gamma * cv) / Pr.
     // TODO: We are using a fixed value of the Prandtl number which is
     // valid for air, it may or may not depend on temperature?  Since
     // this is a property of the fluid, it could possibly be moved to
     // the FluidProperties module...
     const Real Pr = 0.71;
     _thermal_conductivity[qp] = (_dynamic_viscosity[qp] * _fp.cp()) / Pr;
  
     // Compute stabilization parameters:
  
     // .) Compute SUPG element length scale.
     computeHSUPG(qp);
     // Moose::out << "_hsupg[" << qp << "]=" << _hsupg[qp] << std::endl;
  
     // .) Compute SUPG parameter values.  (Must call this after computeHSUPG())
     computeTau(qp);
     // Moose::out << "_tauc[" << qp << "]=" << _tauc[qp] << ", ";
     // Moose::out << "_taum[" << qp << "]=" << _taum[qp] << ", ";
     // Moose::out << "_taue[" << qp << "]=" << _taue[qp] << std::endl;
  
     // .) Compute strong residual values.
     computeStrongResiduals(qp);
     // Moose::out << "_strong_residuals[" << qp << "]=";
     // for (unsigned i=0; i<_strong_residuals[qp].size(); ++i)
     //   Moose::out << _strong_residuals[qp][i] << " ";
     // Moose::out << std::endl;
   }
 }

Referenced by Air::computeProperties().

◆ computeStrongResiduals()

void NavierStokesMaterial::computeStrongResiduals ( unsigned int qp )

private

Referenced by computeProperties().

◆ computeTau()

void NavierStokesMaterial::computeTau ( unsigned int qp )

private

Referenced by computeProperties().

Member Data Documentation

◆ _calA

MaterialProperty<std::vector<RealTensorValue> >& NavierStokesMaterial::_calA

protected

Definition at line 61 of file NavierStokesMaterial.h.

◆ _calC

MaterialProperty<std::vector<RealTensorValue> >& NavierStokesMaterial::_calC

protected

Definition at line 65 of file NavierStokesMaterial.h.

◆ _calE

MaterialProperty<std::vector<std::vector<RealTensorValue> > >& NavierStokesMaterial::_calE

protected

Definition at line 70 of file NavierStokesMaterial.h.

◆ _drho_dt

const VariableValue& NavierStokesMaterial::_drho_dt

protected

Definition at line 96 of file NavierStokesMaterial.h.

◆ _drhoE_dt

const VariableValue& NavierStokesMaterial::_drhoE_dt

protected

Definition at line 100 of file NavierStokesMaterial.h.

◆ _drhou_dt

const VariableValue& NavierStokesMaterial::_drhou_dt

protected

Definition at line 97 of file NavierStokesMaterial.h.

◆ _drhov_dt

const VariableValue& NavierStokesMaterial::_drhov_dt

protected

Definition at line 98 of file NavierStokesMaterial.h.

◆ _drhow_dt

const VariableValue& NavierStokesMaterial::_drhow_dt

protected

Definition at line 99 of file NavierStokesMaterial.h.

◆ _dynamic_viscosity

MaterialProperty<Real>& NavierStokesMaterial::_dynamic_viscosity

protected

Definition at line 55 of file NavierStokesMaterial.h.

Referenced by Air::computeProperties(), and computeProperties().

◆ _enthalpy

const VariableValue& NavierStokesMaterial::_enthalpy

protected

Definition at line 86 of file NavierStokesMaterial.h.

◆ _fp

const IdealGasFluidProperties& NavierStokesMaterial::_fp

protected

Definition at line 121 of file NavierStokesMaterial.h.

Referenced by computeProperties().

◆ _grad_rho

const VariableGradient& NavierStokesMaterial::_grad_rho

protected

Definition at line 103 of file NavierStokesMaterial.h.

◆ _grad_rho_E

const VariableGradient& NavierStokesMaterial::_grad_rho_E

protected

Definition at line 107 of file NavierStokesMaterial.h.

◆ _grad_rho_u

const VariableGradient& NavierStokesMaterial::_grad_rho_u

protected

Definition at line 104 of file NavierStokesMaterial.h.

◆ _grad_rho_v

const VariableGradient& NavierStokesMaterial::_grad_rho_v

protected

Definition at line 105 of file NavierStokesMaterial.h.

◆ _grad_rho_w

const VariableGradient& NavierStokesMaterial::_grad_rho_w

protected

Definition at line 106 of file NavierStokesMaterial.h.

◆ _grad_u

const VariableGradient& NavierStokesMaterial::_grad_u

protected

Definition at line 49 of file NavierStokesMaterial.h.

Referenced by computeProperties(), and NavierStokesMaterial().

◆ _grad_v

const VariableGradient& NavierStokesMaterial::_grad_v

protected

Definition at line 50 of file NavierStokesMaterial.h.

Referenced by computeProperties(), and NavierStokesMaterial().

◆ _grad_w

const VariableGradient& NavierStokesMaterial::_grad_w

protected

Definition at line 51 of file NavierStokesMaterial.h.

Referenced by computeProperties(), and NavierStokesMaterial().

◆ _hsupg

MaterialProperty<Real>& NavierStokesMaterial::_hsupg

protected

Definition at line 111 of file NavierStokesMaterial.h.

Referenced by computeHSUPG().

◆ _mesh_dimension

const unsigned int NavierStokesMaterial::_mesh_dimension

protected

Definition at line 47 of file NavierStokesMaterial.h.

◆ _rho

const VariableValue& NavierStokesMaterial::_rho

protected

Definition at line 89 of file NavierStokesMaterial.h.

◆ _rho_E

const VariableValue& NavierStokesMaterial::_rho_E

protected

Definition at line 93 of file NavierStokesMaterial.h.

◆ _rho_u

const VariableValue& NavierStokesMaterial::_rho_u

protected

Definition at line 90 of file NavierStokesMaterial.h.

◆ _rho_v

const VariableValue& NavierStokesMaterial::_rho_v

protected

Definition at line 91 of file NavierStokesMaterial.h.

◆ _rho_w

const VariableValue& NavierStokesMaterial::_rho_w

protected

Definition at line 92 of file NavierStokesMaterial.h.

◆ _strong_residuals

MaterialProperty<std::vector<Real> >& NavierStokesMaterial::_strong_residuals

protected

Definition at line 118 of file NavierStokesMaterial.h.

◆ _tauc

MaterialProperty<Real>& NavierStokesMaterial::_tauc

protected

Definition at line 112 of file NavierStokesMaterial.h.

◆ _taue

MaterialProperty<Real>& NavierStokesMaterial::_taue

protected

Definition at line 114 of file NavierStokesMaterial.h.

◆ _taum

MaterialProperty<Real>& NavierStokesMaterial::_taum

protected

Definition at line 113 of file NavierStokesMaterial.h.

◆ _temperature

const VariableValue& NavierStokesMaterial::_temperature

protected

Definition at line 83 of file NavierStokesMaterial.h.

◆ _thermal_conductivity

MaterialProperty<Real>& NavierStokesMaterial::_thermal_conductivity

protected

Definition at line 54 of file NavierStokesMaterial.h.

Referenced by computeProperties().

◆ _u_vel

const VariableValue& NavierStokesMaterial::_u_vel

protected

Definition at line 78 of file NavierStokesMaterial.h.

◆ _v_vel

const VariableValue& NavierStokesMaterial::_v_vel

protected

Definition at line 79 of file NavierStokesMaterial.h.

◆ _vel_grads

std::vector<const VariableGradient *> NavierStokesMaterial::_vel_grads

protected

Definition at line 74 of file NavierStokesMaterial.h.

Referenced by computeProperties().

◆ _viscous_stress_tensor

MaterialProperty<RealTensorValue>& NavierStokesMaterial::_viscous_stress_tensor

protected

Definition at line 53 of file NavierStokesMaterial.h.

Referenced by computeProperties().

◆ _w_vel

const VariableValue& NavierStokesMaterial::_w_vel

protected

Definition at line 80 of file NavierStokesMaterial.h.

The documentation for this class was generated from the following files:

navier_stokes/include/materials/NavierStokesMaterial.h
navier_stokes/src/materials/NavierStokesMaterial.C

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ NavierStokesMaterial()

Member Function Documentation

◆ computeHSUPG()

◆ computeProperties()

◆ computeStrongResiduals()

◆ computeTau()

Member Data Documentation

◆ _calA

◆ _calC

◆ _calE

◆ _drho_dt

◆ _drhoE_dt

◆ _drhou_dt

◆ _drhov_dt

◆ _drhow_dt

◆ _dynamic_viscosity

◆ _enthalpy

◆ _fp

◆ _grad_rho

◆ _grad_rho_E

◆ _grad_rho_u

◆ _grad_rho_v

◆ _grad_rho_w

◆ _grad_u

◆ _grad_v

◆ _grad_w

◆ _hsupg

◆ _mesh_dimension

◆ _rho

◆ _rho_E

◆ _rho_u

◆ _rho_v

◆ _rho_w

◆ _strong_residuals

◆ _tauc

◆ _taue

◆ _taum

◆ _temperature

◆ _thermal_conductivity

◆ _u_vel

◆ _v_vel

◆ _vel_grads

◆ _viscous_stress_tensor

◆ _w_vel