INSBase.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 "Kernel.h"
 
class INSBase : public Kernel
{
public:
  static InputParameters validParams();
 
  INSBase(const InputParameters & parameters);
 
  virtual ~INSBase() {}
 
protected:
  virtual Real computeQpResidual() = 0;
  virtual Real computeQpJacobian() = 0;
  virtual Real computeQpOffDiagJacobian(unsigned jvar) = 0;
 
  virtual RealVectorValue convectiveTerm();
  virtual RealVectorValue dConvecDUComp(unsigned comp);
 
  virtual RealVectorValue strongViscousTermLaplace();
  virtual RealVectorValue strongViscousTermTraction();
  virtual RealVectorValue dStrongViscDUCompLaplace(unsigned comp);
  virtual RealVectorValue dStrongViscDUCompTraction(unsigned comp);
 
  virtual RealVectorValue weakViscousTermLaplace(unsigned comp);
  virtual RealVectorValue weakViscousTermTraction(unsigned comp);
  virtual RealVectorValue dWeakViscDUCompLaplace();
  virtual RealVectorValue dWeakViscDUCompTraction();
 
  virtual RealVectorValue strongPressureTerm();
  virtual Real weakPressureTerm();
  virtual RealVectorValue dStrongPressureDPressure();
  virtual Real dWeakPressureDPressure();
 
  virtual RealVectorValue gravityTerm();
 
  virtual RealVectorValue timeDerivativeTerm();
  virtual RealVectorValue dTimeDerivativeDUComp(unsigned comp);
 
  virtual Real tau();
  virtual Real dTauDUComp(unsigned comp);
 
  virtual Real tauNodal();
 
  RealVectorValue relativeVelocity() const;
 
  RealVectorValue strongViscousTermLaplaceRZ() const;
 
  RealVectorValue dStrongViscDUCompLaplaceRZ(const unsigned int comp) const;
 
  RealVectorValue strongViscousTermTractionRZ() const;
 
  RealVectorValue dStrongViscDUCompTractionRZ(const unsigned int comp) const;
 
  const VariablePhiSecond & _second_phi;
 
  // Coupled variables
  const VariableValue & _u_vel;
  const VariableValue & _v_vel;
  const VariableValue & _w_vel;
  const VariableValue & _p;
 
  const bool _picard;
  const VariableValue * const _u_vel_previous_nl;
  const VariableValue * const _v_vel_previous_nl;
  const VariableValue * const _w_vel_previous_nl;
 
  // Gradients
  const VariableGradient & _grad_u_vel;
  const VariableGradient & _grad_v_vel;
  const VariableGradient & _grad_w_vel;
  const VariableGradient & _grad_p;
 
  // Seconds
  const VariableSecond & _second_u_vel;
  const VariableSecond & _second_v_vel;
  const VariableSecond & _second_w_vel;
 
  // Time derivatives
  const VariableValue & _u_vel_dot;
  const VariableValue & _v_vel_dot;
  const VariableValue & _w_vel_dot;
 
  // Derivatives of time derivatives
  const VariableValue & _d_u_vel_dot_du;
  const VariableValue & _d_v_vel_dot_dv;
  const VariableValue & _d_w_vel_dot_dw;
 
  // Variable numberings
  unsigned _u_vel_var_number;
  unsigned _v_vel_var_number;
  unsigned _w_vel_var_number;
  unsigned _p_var_number;
 
  RealVectorValue _gravity;
 
  // Material properties
  const MaterialProperty<Real> & _mu;
  const MaterialProperty<Real> & _rho;
 
  const Real & _alpha;
  bool _laplace;
  bool _convective_term;
  bool _transient_term;
 
  const bool _disps_provided;
  const VariableValue & _disp_x_dot;
  const VariableValue & _disp_y_dot;
  const VariableValue & _disp_z_dot;
 
  const unsigned int _rz_radial_coord;
};