INSADMaterial.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"
#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;
};