doxygen/modules/NavierStokesMaterial_8h_source.html

//* This file is part of the MOOSE framework

//* https://www.mooseframework.org

//*

//* All rights reserved, see COPYRIGHT for full restrictions

//* https://github.com/idaholab/moose/blob/master/COPYRIGHT

//*

//* Licensed under LGPL 2.1, please see LICENSE for details

//* https://www.gnu.org/licenses/lgpl-2.1.html


#pragma once


#include "Material.h"


// Forward Declarations

class NavierStokesMaterial;

class IdealGasFluidProperties;


template <>

InputParameters validParams<NavierStokesMaterial>();


class NavierStokesMaterial : public Material

{

public:

  NavierStokesMaterial(const InputParameters & parameters);


protected:

  virtual void computeProperties();


  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;


  // Also store the "cal A" matrices at each quadrature point as material

  // poperties.  These are defined in the compns notes and are needed for

  // computing strong and weak residual values and Jacobian entries, so it's

  // good if we reuse them...

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


  // The "velocity column" matrices.  _calC[_qp][k] is a tensor with the

  // velocity vector in the k'th column.  See notes for additional details.

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


  // The energy equation inviscid flux matrix components.  There are n_vars of

  // these for each dimension, so in 3D, 3*5=15 different matrices.

  // See notes for additional details.

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


  // Convenient storage for all of the velocity gradient components so

  // we can refer to them in a loop.

  std::vector<const VariableGradient *> _vel_grads;


  // Coupled values needed to compute strong form residuals

  // for SUPG stabilization...

  const VariableValue & _u_vel;

  const VariableValue & _v_vel;

  const VariableValue & _w_vel;


  // Temperature is needed to compute speed of sound

  const VariableValue & _temperature;


  // Enthalpy is needed in computing energy equation strong residuals

  const VariableValue & _enthalpy;


  // Main solution variables are all needed for computing strong residuals

  const VariableValue & _rho;

  const VariableValue & _rho_u;

  const VariableValue & _rho_v;

  const VariableValue & _rho_w;

  const VariableValue & _rho_E;


  // Time derivative values for dependent variables

  const VariableValue & _drho_dt;

  const VariableValue & _drhou_dt;

  const VariableValue & _drhov_dt;

  const VariableValue & _drhow_dt;

  const VariableValue & _drhoE_dt;


  // Gradients

  const VariableGradient & _grad_rho;

  const VariableGradient & _grad_rho_u;

  const VariableGradient & _grad_rho_v;

  const VariableGradient & _grad_rho_w;

  const VariableGradient & _grad_rho_E;


  // The real-valued material properties representing the element stabilization

  // parameters for each of the equations.

  MaterialProperty<Real> & _hsupg;

  MaterialProperty<Real> & _tauc;

  MaterialProperty<Real> & _taum;

  MaterialProperty<Real> & _taue;


  // The (vector-valued) material property which is the strong-form

  // residual at each quadrature point.

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


  // Fluid properties

  const IdealGasFluidProperties & _fp;


private:

  // To be called from computeProperties() function to compute _hsupg

  void computeHSUPG(unsigned int qp);


  // To be called from computeProperties() function to compute _tauc, _taum, _taue

  void computeTau(unsigned int qp);


  // To be called from computeProperties() function to compute the strong residual of each equation.

  void computeStrongResiduals(unsigned int qp);

};