INSADMaterial.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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"
#include "NSEnums.h"
 
class INSADObjectTracker;
 
class INSADMaterial : public Material
{
public:
  static InputParameters validParams();
 
  INSADMaterial(const InputParameters & parameters);
 
  virtual void subdomainSetup() override;
  virtual void resolveOptionalProperties() override;
 
protected:
  virtual void computeQpProperties() override;
 
  const ADVectorVariableValue & _velocity;
 
  const ADVectorVariableGradient & _grad_velocity;
 
  const ADVariableGradient & _grad_p;
 
  const ADMaterialProperty<Real> & _mu;
 
  const ADMaterialProperty<Real> & _rho;
 
  const ADVectorVariableValue * _velocity_dot;
 
  ADMaterialProperty<Real> & _mass_strong_residual;
 
  ADMaterialProperty<RealVectorValue> & _advective_strong_residual;
 
  ADMaterialProperty<RealVectorValue> & _td_strong_residual;
 
  ADMaterialProperty<RealVectorValue> & _gravity_strong_residual;
 
  ADMaterialProperty<RealVectorValue> & _boussinesq_strong_residual;
 
  ADMaterialProperty<RealVectorValue> & _coupled_force_strong_residual;
 
  ADMaterialProperty<RealVectorValue> & _mesh_velocity;
 
  ADMaterialProperty<RealVectorValue> & _relative_velocity;
 
  ADMaterialProperty<RealVectorValue> & _advected_mesh_strong_residual;
 
  // /// Future addition pending addition of INSADMMSKernel.
  // /// Strong residual corresponding to the mms function term
  // MaterialProperty<RealVectorValue> & _mms_function_strong_residual;
 
  const bool _use_displaced_mesh;
 
  const MooseArray<ADPoint> & _ad_q_point;
 
  const unsigned int _rz_radial_coord;
 
  const unsigned int _rz_axial_coord;
 
  const INSADObjectTracker * _object_tracker;
 
  bool _has_transient;
 
  bool _has_gravity;
 
  bool _has_boussinesq;
 
  bool _has_coupled_force;
 
  std::unordered_map<SubdomainID, const ADMaterialProperty<Real> *> _boussinesq_alphas;
  const ADMaterialProperty<Real> * _boussinesq_alpha;
 
  const ADVariableValue * _temperature;
 
  std::unordered_map<SubdomainID, const MaterialProperty<Real> *> _ref_temps;
  const MaterialProperty<Real> * _ref_temp;
 
  NS::ViscousForm _viscous_form;
 
  RealVectorValue _gravity_vector;
 
  std::vector<const ADVectorVariableValue *> _coupled_force_var;
 
  std::vector<const Function *> _coupled_force_vector_function;
 
  bool _has_advected_mesh;
 
  const ADVariableValue * _disp_x_dot;
 
  const ADVariableValue * _disp_y_dot;
 
  const ADVariableValue * _disp_z_dot;
 
  unsigned int _disp_x_num = libMesh::invalid_uint;
  unsigned int _disp_y_num = libMesh::invalid_uint;
  unsigned int _disp_z_num = libMesh::invalid_uint;
  unsigned int _disp_x_sys_num = libMesh::invalid_uint;
  unsigned int _disp_y_sys_num = libMesh::invalid_uint;
  unsigned int _disp_z_sys_num = libMesh::invalid_uint;
};