NavierStokesMaterial.h

Go to the documentation of this file.

//* 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 IdealGasFluidProperties;

class NavierStokesMaterial : public Material
{
public:
  static InputParameters validParams();

  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;

  // Specific total enthalpy is needed in computing energy equation strong residuals
  const VariableValue & _specific_total_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_et;

  // 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_et;

  // 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);
};