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