This class computes momentum equation residual and Jacobian viscous contributions for the "traction" form of the governing equations. More...

#include <INSMomentumTractionForm.h>

Inheritance diagram for INSMomentumTractionForm:

[legend]

Public Member Functions
	INSMomentumTractionForm (const InputParameters &parameters)

virtual	~INSMomentumTractionForm ()

Protected Member Functions
virtual Real	computeQpResidualViscousPart () override

virtual Real	computeQpJacobianViscousPart () override

virtual Real	computeQpOffDiagJacobianViscousPart (unsigned jvar) override

virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

virtual Real	computeQpOffDiagJacobian (unsigned jvar)

virtual Real	computeQpPGResidual ()

virtual Real	computeQpPGJacobian (unsigned comp)

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 ()
	Provides tau which yields superconvergence for 1D advection-diffusion. More...

Protected Attributes
unsigned	_component

bool	_integrate_p_by_parts

bool	_supg

const Function &	_ffn

const VariablePhiSecond &	_second_phi
	second derivatives of the shape function More...

const VariableValue &	_u_vel

const VariableValue &	_v_vel

const VariableValue &	_w_vel

const VariableValue &	_p

const VariableGradient &	_grad_u_vel

const VariableGradient &	_grad_v_vel

const VariableGradient &	_grad_w_vel

const VariableGradient &	_grad_p

const VariableSecond &	_second_u_vel

const VariableSecond &	_second_v_vel

const VariableSecond &	_second_w_vel

const VariableValue &	_u_vel_dot

const VariableValue &	_v_vel_dot

const VariableValue &	_w_vel_dot

const VariableValue &	_d_u_vel_dot_du

const VariableValue &	_d_v_vel_dot_dv

const VariableValue &	_d_w_vel_dot_dw

unsigned	_u_vel_var_number

unsigned	_v_vel_var_number

unsigned	_w_vel_var_number

unsigned	_p_var_number

RealVectorValue	_gravity

const MaterialProperty< Real > &	_mu

const MaterialProperty< Real > &	_rho

const Real &	_alpha

bool	_laplace

bool	_convective_term

bool	_transient_term

Detailed Description

This class computes momentum equation residual and Jacobian viscous contributions for the "traction" form of the governing equations.

Definition at line 24 of file INSMomentumTractionForm.h.

Constructor & Destructor Documentation

◆ INSMomentumTractionForm()

INSMomentumTractionForm::INSMomentumTractionForm ( const InputParameters & parameters )

Definition at line 24 of file INSMomentumTractionForm.C.

   : INSMomentumBase(parameters)
 {
 }

◆ ~INSMomentumTractionForm()

virtual INSMomentumTractionForm::~INSMomentumTractionForm ( )

inlinevirtual

Definition at line 29 of file INSMomentumTractionForm.h.

29 {}

Member Function Documentation

◆ computeQpJacobian()

Real INSMomentumBase::computeQpJacobian ( )

protectedvirtualinherited

Implements INSBase.

Reimplemented in INSMomentumLaplaceFormRZ, and INSMomentumTractionFormRZ.

Definition at line 93 of file INSMomentumBase.C.

 {
   Real jac = 0;
  
   // viscous term
   jac += computeQpJacobianViscousPart();
  
   // convective term
   if (_convective_term)
     jac += _test[_i][_qp] * dConvecDUComp(_component)(_component);
  
   if (_supg)
     jac += computeQpPGJacobian(_component);
  
   return jac;
 }

Referenced by INSMomentumTractionFormRZ::computeQpJacobian(), and INSMomentumLaplaceFormRZ::computeQpJacobian().

◆ computeQpJacobianViscousPart()

Real INSMomentumTractionForm::computeQpJacobianViscousPart ( )

overrideprotectedvirtual

Implements INSMomentumBase.

Definition at line 64 of file INSMomentumTractionForm.C.

 {
   // Viscous part, full stress tensor.  The extra contribution comes from the "2"
   // on the diagonal of the viscous stress tensor.
   return _mu[_qp] * (_grad_phi[_j][_qp] * _grad_test[_i][_qp] +
                      _grad_phi[_j][_qp](_component) * _grad_test[_i][_qp](_component));
 }

◆ computeQpOffDiagJacobian()

Real INSMomentumBase::computeQpOffDiagJacobian ( unsigned jvar )

protectedvirtualinherited

Implements INSBase.

Reimplemented in INSMomentumLaplaceFormRZ, and INSMomentumTractionFormRZ.

Definition at line 137 of file INSMomentumBase.C.

 {
   Real jac = 0;
   if (jvar == _u_vel_var_number)
   {
     Real convective_term = _convective_term ? _test[_i][_qp] * dConvecDUComp(0)(_component) : 0.;
     Real viscous_term = computeQpOffDiagJacobianViscousPart(jvar);
  
     jac += convective_term + viscous_term;
  
     if (_supg)
       jac += computeQpPGJacobian(0);
  
     return jac;
   }
   else if (jvar == _v_vel_var_number)
   {
     Real convective_term = _convective_term ? _test[_i][_qp] * dConvecDUComp(1)(_component) : 0.;
     Real viscous_term = computeQpOffDiagJacobianViscousPart(jvar);
  
     jac += convective_term + viscous_term;
  
     if (_supg)
       jac += computeQpPGJacobian(1);
  
     return jac;
   }
   else if (jvar == _w_vel_var_number)
   {
     Real convective_term = _convective_term ? _test[_i][_qp] * dConvecDUComp(2)(_component) : 0.;
     Real viscous_term = computeQpOffDiagJacobianViscousPart(jvar);
  
     jac += convective_term + viscous_term;
  
     if (_supg)
       jac += computeQpPGJacobian(2);
  
     return jac;
   }
  
   else if (jvar == _p_var_number)
   {
     if (_integrate_p_by_parts)
       jac += _grad_test[_i][_qp](_component) * dWeakPressureDPressure();
     else
       jac += _test[_i][_qp] * dStrongPressureDPressure()(_component);
  
     if (_supg)
     {
       RealVectorValue U(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
       jac += tau() * U * _grad_test[_i][_qp] * dStrongPressureDPressure()(_component);
     }
  
     return jac;
   }
  
   else
     return 0;
 }

Referenced by INSMomentumTractionFormRZ::computeQpOffDiagJacobian(), and INSMomentumLaplaceFormRZ::computeQpOffDiagJacobian().

◆ computeQpOffDiagJacobianViscousPart()

Real INSMomentumTractionForm::computeQpOffDiagJacobianViscousPart ( unsigned jvar )

overrideprotectedvirtual

Implements INSMomentumBase.

Definition at line 73 of file INSMomentumTractionForm.C.

 {
   // In Stokes/Laplacian version, off-diag Jacobian entries wrt u,v,w are zero
   if (jvar == _u_vel_var_number)
     return _mu[_qp] * _grad_phi[_j][_qp](_component) * _grad_test[_i][_qp](0);
  
   else if (jvar == _v_vel_var_number)
     return _mu[_qp] * _grad_phi[_j][_qp](_component) * _grad_test[_i][_qp](1);
  
   else if (jvar == _w_vel_var_number)
     return _mu[_qp] * _grad_phi[_j][_qp](_component) * _grad_test[_i][_qp](2);
  
   else
     return 0;
 }

◆ computeQpPGJacobian()

Real INSMomentumBase::computeQpPGJacobian ( unsigned comp )

protectedvirtualinherited

Definition at line 111 of file INSMomentumBase.C.

 {
   RealVectorValue U(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
   RealVectorValue d_U_d_U_comp(0, 0, 0);
   d_U_d_U_comp(comp) = _phi[_j][_qp];
  
   Real convective_term = _convective_term ? convectiveTerm()(_component) : 0;
   Real d_convective_term_d_u_comp = _convective_term ? dConvecDUComp(comp)(_component) : 0;
   Real viscous_term =
       _laplace ? strongViscousTermLaplace()(_component) : strongViscousTermTraction()(_component);
   Real d_viscous_term_d_u_comp = _laplace ? dStrongViscDUCompLaplace(comp)(_component)
                                           : dStrongViscDUCompTraction(comp)(_component);
   Real transient_term = _transient_term ? timeDerivativeTerm()(_component) : 0;
   Real d_transient_term_d_u_comp = _transient_term ? dTimeDerivativeDUComp(comp)(_component) : 0;
  
   return dTauDUComp(comp) * U * _grad_test[_i][_qp] *
              (convective_term + viscous_term + strongPressureTerm()(_component) +
               gravityTerm()(_component) + transient_term - _ffn.value(_t, _q_point[_qp])) +
          tau() * d_U_d_U_comp * _grad_test[_i][_qp] *
              (convective_term + viscous_term + strongPressureTerm()(_component) +
               gravityTerm()(_component) + transient_term - _ffn.value(_t, _q_point[_qp])) +
          tau() * U * _grad_test[_i][_qp] *
              (d_convective_term_d_u_comp + d_viscous_term_d_u_comp + d_transient_term_d_u_comp);
 }

Referenced by INSMomentumBase::computeQpJacobian(), and INSMomentumBase::computeQpOffDiagJacobian().

◆ computeQpPGResidual()

Real INSMomentumBase::computeQpPGResidual ( )

protectedvirtualinherited

Definition at line 67 of file INSMomentumBase.C.

 {
   RealVectorValue U(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
   RealVectorValue convective_term = _convective_term ? convectiveTerm() : RealVectorValue(0, 0, 0);
   RealVectorValue viscous_term =
       _laplace ? strongViscousTermLaplace() : strongViscousTermTraction();
   RealVectorValue transient_term =
       _transient_term ? timeDerivativeTerm() : RealVectorValue(0, 0, 0);
  
   return tau() * U * _grad_test[_i][_qp] *
          ((convective_term + viscous_term + transient_term + strongPressureTerm() +
            gravityTerm())(_component)-_ffn.value(_t, _q_point[_qp]));
  
   // For GLS as opposed to SUPG stabilization, one would need to modify the test function functional
   // space to include second derivatives of the Galerkin test functions corresponding to the viscous
   // term. This would look like:
   // Real lap_test =
   //     _second_test[_i][_qp](0, 0) + _second_test[_i][_qp](1, 1) + _second_test[_i][_qp](2, 2);
  
   // Real pg_viscous_r = -_mu[_qp] * lap_test * tau() *
   //                     (convective_term + viscous_term + strongPressureTerm()(_component) +
   //                      gravityTerm())(_component);
 }

Referenced by INSMomentumBase::computeQpResidual().

◆ computeQpResidual()

Real INSMomentumBase::computeQpResidual ( )

protectedvirtualinherited

Implements INSBase.

Reimplemented in INSMomentumLaplaceFormRZ, and INSMomentumTractionFormRZ.

Definition at line 40 of file INSMomentumBase.C.

 {
   Real r = 0;
  
   // viscous term
   r += computeQpResidualViscousPart();
  
   // pressure term
   if (_integrate_p_by_parts)
     r += _grad_test[_i][_qp](_component) * weakPressureTerm();
   else
     r += _test[_i][_qp] * strongPressureTerm()(_component);
  
   // body force term
   r += _test[_i][_qp] * (gravityTerm()(_component) - _ffn.value(_t, _q_point[_qp]));
  
   // convective term
   if (_convective_term)
     r += _test[_i][_qp] * convectiveTerm()(_component);
  
   if (_supg)
     r += computeQpPGResidual();
  
   return r;
 }

Referenced by INSMomentumTractionFormRZ::computeQpResidual(), and INSMomentumLaplaceFormRZ::computeQpResidual().

◆ computeQpResidualViscousPart()

Real INSMomentumTractionForm::computeQpResidualViscousPart ( )

overrideprotectedvirtual

Implements INSMomentumBase.

Definition at line 30 of file INSMomentumTractionForm.C.

 {
   // The component'th row (or col, it's symmetric) of the viscous stress tensor
   RealVectorValue tau_row;
  
   switch (_component)
   {
     case 0:
       tau_row(0) = 2. * _grad_u_vel[_qp](0);                  // 2*du/dx1
       tau_row(1) = _grad_u_vel[_qp](1) + _grad_v_vel[_qp](0); // du/dx2 + dv/dx1
       tau_row(2) = _grad_u_vel[_qp](2) + _grad_w_vel[_qp](0); // du/dx3 + dw/dx1
       break;
  
     case 1:
       tau_row(0) = _grad_v_vel[_qp](0) + _grad_u_vel[_qp](1); // dv/dx1 + du/dx2
       tau_row(1) = 2. * _grad_v_vel[_qp](1);                  // 2*dv/dx2
       tau_row(2) = _grad_v_vel[_qp](2) + _grad_w_vel[_qp](1); // dv/dx3 + dw/dx2
       break;
  
     case 2:
       tau_row(0) = _grad_w_vel[_qp](0) + _grad_u_vel[_qp](2); // dw/dx1 + du/dx3
       tau_row(1) = _grad_w_vel[_qp](1) + _grad_v_vel[_qp](2); // dw/dx2 + dv/dx3
       tau_row(2) = 2. * _grad_w_vel[_qp](2);                  // 2*dw/dx3
       break;
  
     default:
       mooseError("Unrecognized _component requested.");
   }
  
   // The viscous part, _mu[_qp] * tau : grad(v)
   return _mu[_qp] * (tau_row * _grad_test[_i][_qp]);
 }

◆ convectiveTerm()

RealVectorValue INSBase::convectiveTerm ( )

protectedvirtualinherited

Definition at line 96 of file INSBase.C.

 {
   RealVectorValue U(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
   return _rho[_qp] *
          RealVectorValue(U * _grad_u_vel[_qp], U * _grad_v_vel[_qp], U * _grad_w_vel[_qp]);
 }

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), and INSMomentumBase::computeQpResidual().

◆ dConvecDUComp()

RealVectorValue INSBase::dConvecDUComp ( unsigned comp )

protectedvirtualinherited

Definition at line 104 of file INSBase.C.

 {
   RealVectorValue U(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
   RealVectorValue d_U_d_comp(0, 0, 0);
   d_U_d_comp(comp) = _phi[_j][_qp];
  
   RealVectorValue convective_term = _rho[_qp] * RealVectorValue(d_U_d_comp * _grad_u_vel[_qp],
                                                                 d_U_d_comp * _grad_v_vel[_qp],
                                                                 d_U_d_comp * _grad_w_vel[_qp]);
   convective_term(comp) += _rho[_qp] * U * _grad_phi[_j][_qp];
  
   return convective_term;
 }

Referenced by INSMomentumBase::computeQpJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumBase::computeQpPGJacobian(), and INSMass::computeQpPGOffDiagJacobian().

◆ dStrongPressureDPressure()

RealVectorValue INSBase::dStrongPressureDPressure ( )

protectedvirtualinherited

Definition at line 226 of file INSBase.C.

 {
   return _grad_phi[_j][_qp];
 }

Referenced by INSMomentumBase::computeQpOffDiagJacobian(), and INSMass::computeQpPGJacobian().

◆ dStrongViscDUCompLaplace()

RealVectorValue INSBase::dStrongViscDUCompLaplace ( unsigned comp )

protectedvirtualinherited

Reimplemented in INSMassRZ, and INSMomentumLaplaceFormRZ.

Definition at line 134 of file INSBase.C.

 {
   RealVectorValue viscous_term(0, 0, 0);
   viscous_term(comp) = -_mu[_qp] * _second_phi[_j][_qp].tr();
  
   return viscous_term;
 }

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMomentumLaplaceFormRZ::dStrongViscDUCompLaplace(), and INSMassRZ::dStrongViscDUCompLaplace().

◆ dStrongViscDUCompTraction()

RealVectorValue INSBase::dStrongViscDUCompTraction ( unsigned comp )

protectedvirtualinherited

Reimplemented in INSMassRZ, and INSMomentumTractionFormRZ.

Definition at line 143 of file INSBase.C.

 {
   RealVectorValue viscous_term(0, 0, 0);
   viscous_term(comp) = -_mu[_qp] * (_second_phi[_j][_qp](0, 0) + _second_phi[_j][_qp](1, 1) +
                                     _second_phi[_j][_qp](2, 2));
   for (unsigned i = 0; i < 3; i++)
     viscous_term(i) += -_mu[_qp] * _second_phi[_j][_qp](i, comp);
  
   return viscous_term;
 }

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMomentumTractionFormRZ::dStrongViscDUCompTraction(), and INSMassRZ::dStrongViscDUCompTraction().

◆ dTauDUComp()

Real INSBase::dTauDUComp ( unsigned comp )

protectedvirtualinherited

Definition at line 298 of file INSBase.C.

 {
   Real nu = _mu[_qp] / _rho[_qp];
   RealVectorValue U(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
   Real h = _current_elem->hmax();
   Real transient_part = _transient_term ? 4. / (_dt * _dt) : 0.;
   return -_alpha / 2. * std::pow(transient_part + (2. * U.norm() / h) * (2. * U.norm() / h) +
                                      9. * (4. * nu / (h * h)) * (4. * nu / (h * h)),
                                  -1.5) *
          2. * (2. * U.norm() / h) * 2. / h * U(comp) * _phi[_j][_qp] /
          (U.norm() + std::numeric_limits<double>::epsilon());
 }

Referenced by INSMomentumBase::computeQpPGJacobian(), and INSMass::computeQpPGOffDiagJacobian().

◆ dTimeDerivativeDUComp()

RealVectorValue INSBase::dTimeDerivativeDUComp ( unsigned comp )

protectedvirtualinherited

Definition at line 250 of file INSBase.C.

 {
   Real base = _rho[_qp] * _phi[_j][_qp];
   switch (comp)
   {
     case 0:
       return RealVectorValue(base * _d_u_vel_dot_du[_qp], 0, 0);
  
     case 1:
       return RealVectorValue(0, base * _d_v_vel_dot_dv[_qp], 0);
  
     case 2:
       return RealVectorValue(0, 0, base * _d_w_vel_dot_dw[_qp]);
  
     default:
       mooseError("comp must be 0, 1, or 2");
   }
 }

Referenced by INSMomentumBase::computeQpPGJacobian(), and INSMass::computeQpPGOffDiagJacobian().

◆ dWeakPressureDPressure()

Real INSBase::dWeakPressureDPressure ( )

protectedvirtualinherited

Definition at line 232 of file INSBase.C.

 {
   return -_phi[_j][_qp];
 }

Referenced by INSMomentumBase::computeQpOffDiagJacobian().

◆ dWeakViscDUCompLaplace()

RealVectorValue INSBase::dWeakViscDUCompLaplace ( )

protectedvirtualinherited

Definition at line 202 of file INSBase.C.

 {
   return _mu[_qp] * _grad_phi[_j][_qp];
 }

◆ dWeakViscDUCompTraction()

RealVectorValue INSBase::dWeakViscDUCompTraction ( )

protectedvirtualinherited

Definition at line 208 of file INSBase.C.

 {
   return _mu[_qp] * _grad_phi[_j][_qp];
 }

◆ gravityTerm()

RealVectorValue INSBase::gravityTerm ( )

protectedvirtualinherited

Definition at line 238 of file INSBase.C.

 {
   return -_rho[_qp] * _gravity;
 }

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), and INSMomentumBase::computeQpResidual().

◆ strongPressureTerm()

RealVectorValue INSBase::strongPressureTerm ( )

protectedvirtualinherited

Definition at line 214 of file INSBase.C.

 {
   return _grad_p[_qp];
 }

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), and INSMomentumBase::computeQpResidual().

◆ strongViscousTermLaplace()

RealVectorValue INSBase::strongViscousTermLaplace ( )

protectedvirtualinherited

Reimplemented in INSMassRZ, and INSMomentumLaplaceFormRZ.

Definition at line 119 of file INSBase.C.

 {
   return -_mu[_qp] *
          RealVectorValue(_second_u_vel[_qp].tr(), _second_v_vel[_qp].tr(), _second_w_vel[_qp].tr());
 }

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), INSMomentumLaplaceFormRZ::strongViscousTermLaplace(), INSMassRZ::strongViscousTermLaplace(), and INSBase::strongViscousTermTraction().

◆ strongViscousTermTraction()

RealVectorValue INSBase::strongViscousTermTraction ( )

protectedvirtualinherited

Reimplemented in INSMassRZ, and INSMomentumTractionFormRZ.

Definition at line 126 of file INSBase.C.

 {
   return strongViscousTermLaplace() -
          _mu[_qp] *
              (_second_u_vel[_qp].row(0) + _second_v_vel[_qp].row(1) + _second_w_vel[_qp].row(2));
 }

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), INSMomentumTractionFormRZ::strongViscousTermTraction(), and INSMassRZ::strongViscousTermTraction().

◆ tau()

Real INSBase::tau ( )

protectedvirtualinherited

Definition at line 270 of file INSBase.C.

 {
   Real nu = _mu[_qp] / _rho[_qp];
   RealVectorValue U(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
   Real h = _current_elem->hmax();
   Real transient_part = _transient_term ? 4. / (_dt * _dt) : 0.;
   return _alpha / std::sqrt(transient_part + (2. * U.norm() / h) * (2. * U.norm() / h) +
                             9. * (4. * nu / (h * h)) * (4. * nu / (h * h)));
 }

Referenced by Advection::computeQpJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMass::computeQpPGJacobian(), INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), and Advection::computeQpResidual().

◆ tauNodal()

Real INSBase::tauNodal ( )

protectedvirtualinherited

Provides tau which yields superconvergence for 1D advection-diffusion.

Definition at line 281 of file INSBase.C.

 {
   Real nu = _mu[_qp] / _rho[_qp];
   RealVectorValue U(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
   Real h = _current_elem->hmax();
   Real xi;
   if (nu < std::numeric_limits<Real>::epsilon())
     xi = 1;
   else
   {
     Real alpha = U.norm() * h / (2 * nu);
     xi = 1. / std::tanh(alpha) - 1. / alpha;
   }
   return h / (2 * U.norm()) * xi;
 }

Referenced by Advection::computeQpJacobian(), and Advection::computeQpResidual().

◆ timeDerivativeTerm()

RealVectorValue INSBase::timeDerivativeTerm ( )

protectedvirtualinherited

Definition at line 244 of file INSBase.C.

 {
   return _rho[_qp] * RealVectorValue(_u_vel_dot[_qp], _v_vel_dot[_qp], _w_vel_dot[_qp]);
 }

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), and INSMomentumBase::computeQpPGResidual().

◆ weakPressureTerm()

Real INSBase::weakPressureTerm ( )

protectedvirtualinherited

Definition at line 220 of file INSBase.C.

 {
   return -_p[_qp];
 }

Referenced by INSMomentumBase::computeQpResidual().

◆ weakViscousTermLaplace()

RealVectorValue INSBase::weakViscousTermLaplace ( unsigned comp )

protectedvirtualinherited

Definition at line 155 of file INSBase.C.

 {
   switch (comp)
   {
     case 0:
       return _mu[_qp] * _grad_u_vel[_qp];
  
     case 1:
       return _mu[_qp] * _grad_v_vel[_qp];
  
     case 2:
       return _mu[_qp] * _grad_w_vel[_qp];
  
     default:
       return _zero[_qp];
   }
 }

◆ weakViscousTermTraction()

RealVectorValue INSBase::weakViscousTermTraction ( unsigned comp )

protectedvirtualinherited

Definition at line 174 of file INSBase.C.

 {
   switch (comp)
   {
     case 0:
     {
       RealVectorValue transpose(_grad_u_vel[_qp](0), _grad_v_vel[_qp](0), _grad_w_vel[_qp](0));
       return _mu[_qp] * _grad_u_vel[_qp] + _mu[_qp] * transpose;
     }
  
     case 1:
     {
       RealVectorValue transpose(_grad_u_vel[_qp](1), _grad_v_vel[_qp](1), _grad_w_vel[_qp](1));
       return _mu[_qp] * _grad_v_vel[_qp] + _mu[_qp] * transpose;
     }
  
     case 2:
     {
       RealVectorValue transpose(_grad_u_vel[_qp](2), _grad_v_vel[_qp](2), _grad_w_vel[_qp](2));
       return _mu[_qp] * _grad_w_vel[_qp] + _mu[_qp] * transpose;
     }
  
     default:
       return _zero[_qp];
   }
 }

Member Data Documentation

◆ _alpha

const Real& INSBase::_alpha

protectedinherited

Definition at line 107 of file INSBase.h.

Referenced by INSBase::dTauDUComp(), and INSBase::tau().

◆ _component

unsigned INSMomentumBase::_component

protectedinherited

Definition at line 43 of file INSMomentumBase.h.

Referenced by INSMomentumBase::computeQpJacobian(), INSMomentumTractionFormRZ::computeQpJacobian(), INSMomentumLaplaceFormRZ::computeQpJacobian(), computeQpJacobianViscousPart(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumTractionFormRZ::computeQpOffDiagJacobian(), INSMomentumLaplaceFormRZ::computeQpOffDiagJacobian(), computeQpOffDiagJacobianViscousPart(), INSMomentumBase::computeQpPGJacobian(), INSMomentumBase::computeQpPGResidual(), INSMomentumBase::computeQpResidual(), INSMomentumTractionFormRZ::computeQpResidual(), INSMomentumLaplaceFormRZ::computeQpResidual(), and computeQpResidualViscousPart().

◆ _convective_term

bool INSBase::_convective_term

protectedinherited

Definition at line 109 of file INSBase.h.

Referenced by INSMomentumBase::computeQpJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), INSMomentumBase::computeQpResidual(), and INSMomentumBase::INSMomentumBase().

◆ _d_u_vel_dot_du

const VariableValue& INSBase::_d_u_vel_dot_du

protectedinherited

Definition at line 91 of file INSBase.h.

Referenced by INSBase::dTimeDerivativeDUComp().

◆ _d_v_vel_dot_dv

const VariableValue& INSBase::_d_v_vel_dot_dv

protectedinherited

Definition at line 92 of file INSBase.h.

Referenced by INSBase::dTimeDerivativeDUComp().

◆ _d_w_vel_dot_dw

const VariableValue& INSBase::_d_w_vel_dot_dw

protectedinherited

Definition at line 93 of file INSBase.h.

Referenced by INSBase::dTimeDerivativeDUComp().

◆ _ffn

const Function& INSMomentumBase::_ffn

protectedinherited

Definition at line 46 of file INSMomentumBase.h.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMomentumBase::computeQpPGResidual(), and INSMomentumBase::computeQpResidual().

◆ _grad_p

const VariableGradient& INSBase::_grad_p

protectedinherited

Definition at line 78 of file INSBase.h.

Referenced by INSBase::strongPressureTerm().

◆ _grad_u_vel

const VariableGradient& INSBase::_grad_u_vel

protectedinherited

Definition at line 75 of file INSBase.h.

Referenced by INSMass::computeQpResidual(), computeQpResidualViscousPart(), INSBase::convectiveTerm(), INSBase::dConvecDUComp(), INSMomentumLaplaceFormRZ::strongViscousTermLaplace(), INSMassRZ::strongViscousTermLaplace(), INSMomentumTractionFormRZ::strongViscousTermTraction(), INSMassRZ::strongViscousTermTraction(), INSBase::weakViscousTermLaplace(), and INSBase::weakViscousTermTraction().

◆ _grad_v_vel

const VariableGradient& INSBase::_grad_v_vel

protectedinherited

Definition at line 76 of file INSBase.h.

Referenced by INSMass::computeQpResidual(), computeQpResidualViscousPart(), INSBase::convectiveTerm(), INSBase::dConvecDUComp(), INSMomentumLaplaceFormRZ::strongViscousTermLaplace(), INSMassRZ::strongViscousTermLaplace(), INSMomentumTractionFormRZ::strongViscousTermTraction(), INSMassRZ::strongViscousTermTraction(), INSBase::weakViscousTermLaplace(), and INSBase::weakViscousTermTraction().

◆ _grad_w_vel

const VariableGradient& INSBase::_grad_w_vel

protectedinherited

Definition at line 77 of file INSBase.h.

Referenced by INSMass::computeQpResidual(), computeQpResidualViscousPart(), INSBase::convectiveTerm(), INSBase::dConvecDUComp(), INSBase::weakViscousTermLaplace(), and INSBase::weakViscousTermTraction().

◆ _gravity

RealVectorValue INSBase::_gravity

protectedinherited

Definition at line 101 of file INSBase.h.

Referenced by INSBase::gravityTerm().

◆ _integrate_p_by_parts

bool INSMomentumBase::_integrate_p_by_parts

protectedinherited

Definition at line 44 of file INSMomentumBase.h.

Referenced by INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumTractionFormRZ::computeQpOffDiagJacobian(), INSMomentumLaplaceFormRZ::computeQpOffDiagJacobian(), INSMomentumBase::computeQpResidual(), INSMomentumTractionFormRZ::computeQpResidual(), and INSMomentumLaplaceFormRZ::computeQpResidual().

◆ _laplace

bool INSBase::_laplace

protectedinherited

Definition at line 108 of file INSBase.h.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), and INSMomentumBase::computeQpPGResidual().

◆ _mu

const MaterialProperty<Real>& INSBase::_mu

protectedinherited

Definition at line 104 of file INSBase.h.

Referenced by INSMomentumTractionFormRZ::computeQpJacobian(), INSMomentumLaplaceFormRZ::computeQpJacobian(), computeQpJacobianViscousPart(), INSMomentumLaplaceForm::computeQpJacobianViscousPart(), computeQpOffDiagJacobianViscousPart(), INSMomentumTractionFormRZ::computeQpResidual(), INSMomentumLaplaceFormRZ::computeQpResidual(), INSMomentumLaplaceForm::computeQpResidualViscousPart(), computeQpResidualViscousPart(), INSMomentumLaplaceFormRZ::dStrongViscDUCompLaplace(), INSMassRZ::dStrongViscDUCompLaplace(), INSBase::dStrongViscDUCompLaplace(), INSMomentumTractionFormRZ::dStrongViscDUCompTraction(), INSMassRZ::dStrongViscDUCompTraction(), INSBase::dStrongViscDUCompTraction(), INSBase::dTauDUComp(), INSBase::dWeakViscDUCompLaplace(), INSBase::dWeakViscDUCompTraction(), INSMomentumLaplaceFormRZ::strongViscousTermLaplace(), INSMassRZ::strongViscousTermLaplace(), INSBase::strongViscousTermLaplace(), INSMomentumTractionFormRZ::strongViscousTermTraction(), INSMassRZ::strongViscousTermTraction(), INSBase::strongViscousTermTraction(), INSBase::tau(), INSBase::tauNodal(), INSBase::weakViscousTermLaplace(), and INSBase::weakViscousTermTraction().

◆ _p

const VariableValue& INSBase::_p

protectedinherited

Definition at line 72 of file INSBase.h.

Referenced by INSMomentumTractionFormRZ::computeQpResidual(), INSMomentumLaplaceFormRZ::computeQpResidual(), and INSBase::weakPressureTerm().

◆ _p_var_number

unsigned INSBase::_p_var_number

protectedinherited

Definition at line 99 of file INSBase.h.

Referenced by INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumTractionFormRZ::computeQpOffDiagJacobian(), and INSMomentumLaplaceFormRZ::computeQpOffDiagJacobian().

◆ _rho

const MaterialProperty<Real>& INSBase::_rho

protectedinherited

Definition at line 105 of file INSBase.h.

Referenced by INSMass::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSBase::convectiveTerm(), INSBase::dConvecDUComp(), INSBase::dTauDUComp(), INSBase::dTimeDerivativeDUComp(), INSBase::gravityTerm(), INSBase::tau(), INSBase::tauNodal(), and INSBase::timeDerivativeTerm().

◆ _second_phi

const VariablePhiSecond& INSBase::_second_phi

protectedinherited

second derivatives of the shape function

Definition at line 66 of file INSBase.h.

Referenced by INSBase::dStrongViscDUCompLaplace(), and INSBase::dStrongViscDUCompTraction().

◆ _second_u_vel

const VariableSecond& INSBase::_second_u_vel

protectedinherited

Definition at line 81 of file INSBase.h.

Referenced by INSBase::strongViscousTermLaplace(), and INSBase::strongViscousTermTraction().

◆ _second_v_vel

const VariableSecond& INSBase::_second_v_vel

protectedinherited

Definition at line 82 of file INSBase.h.

Referenced by INSBase::strongViscousTermLaplace(), and INSBase::strongViscousTermTraction().

◆ _second_w_vel

const VariableSecond& INSBase::_second_w_vel

protectedinherited

Definition at line 83 of file INSBase.h.

Referenced by INSBase::strongViscousTermLaplace(), and INSBase::strongViscousTermTraction().

◆ _supg

bool INSMomentumBase::_supg

protectedinherited

Definition at line 45 of file INSMomentumBase.h.

Referenced by INSMomentumBase::computeQpJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumBase::computeQpResidual(), and INSMomentumBase::INSMomentumBase().

◆ _transient_term

bool INSBase::_transient_term

protectedinherited

Definition at line 110 of file INSBase.h.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), INSBase::dTauDUComp(), and INSBase::tau().

◆ _u_vel

const VariableValue& INSBase::_u_vel

protectedinherited

Definition at line 69 of file INSBase.h.

Referenced by Advection::computeQpJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumBase::computeQpPGJacobian(), INSMomentumBase::computeQpPGResidual(), Advection::computeQpResidual(), INSMomentumTractionFormRZ::computeQpResidual(), INSMomentumLaplaceFormRZ::computeQpResidual(), INSMassRZ::computeQpResidual(), INSBase::convectiveTerm(), INSBase::dConvecDUComp(), INSBase::dTauDUComp(), INSMomentumLaplaceFormRZ::strongViscousTermLaplace(), INSMassRZ::strongViscousTermLaplace(), INSMomentumTractionFormRZ::strongViscousTermTraction(), INSMassRZ::strongViscousTermTraction(), INSBase::tau(), and INSBase::tauNodal().

◆ _u_vel_dot

const VariableValue& INSBase::_u_vel_dot

protectedinherited

Definition at line 86 of file INSBase.h.

Referenced by INSBase::timeDerivativeTerm().

◆ _u_vel_var_number

unsigned INSBase::_u_vel_var_number

protectedinherited

Definition at line 96 of file INSBase.h.

Referenced by INSMass::computeQpOffDiagJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMassRZ::computeQpOffDiagJacobian(), and computeQpOffDiagJacobianViscousPart().

◆ _v_vel

const VariableValue& INSBase::_v_vel

protectedinherited

Definition at line 70 of file INSBase.h.

Referenced by Advection::computeQpJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumBase::computeQpPGJacobian(), INSMomentumBase::computeQpPGResidual(), Advection::computeQpResidual(), INSBase::convectiveTerm(), INSBase::dConvecDUComp(), INSBase::dTauDUComp(), INSBase::tau(), and INSBase::tauNodal().

◆ _v_vel_dot

const VariableValue& INSBase::_v_vel_dot

protectedinherited

Definition at line 87 of file INSBase.h.

Referenced by INSBase::timeDerivativeTerm().

◆ _v_vel_var_number

unsigned INSBase::_v_vel_var_number

protectedinherited

Definition at line 97 of file INSBase.h.

Referenced by INSMass::computeQpOffDiagJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), and computeQpOffDiagJacobianViscousPart().

◆ _w_vel

const VariableValue& INSBase::_w_vel

protectedinherited

Definition at line 71 of file INSBase.h.

Referenced by Advection::computeQpJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumBase::computeQpPGJacobian(), INSMomentumBase::computeQpPGResidual(), Advection::computeQpResidual(), INSBase::convectiveTerm(), INSBase::dConvecDUComp(), INSBase::dTauDUComp(), INSBase::tau(), and INSBase::tauNodal().

◆ _w_vel_dot

const VariableValue& INSBase::_w_vel_dot

protectedinherited

Definition at line 88 of file INSBase.h.

Referenced by INSBase::timeDerivativeTerm().

◆ _w_vel_var_number

unsigned INSBase::_w_vel_var_number

protectedinherited

Definition at line 98 of file INSBase.h.

Referenced by INSMass::computeQpOffDiagJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), and computeQpOffDiagJacobianViscousPart().

The documentation for this class was generated from the following files:

navier_stokes/include/kernels/INSMomentumTractionForm.h
navier_stokes/src/kernels/INSMomentumTractionForm.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ INSMomentumTractionForm()

◆ ~INSMomentumTractionForm()

Member Function Documentation

◆ computeQpJacobian()

◆ computeQpJacobianViscousPart()

◆ computeQpOffDiagJacobian()

◆ computeQpOffDiagJacobianViscousPart()

◆ computeQpPGJacobian()

◆ computeQpPGResidual()

◆ computeQpResidual()

◆ computeQpResidualViscousPart()

◆ convectiveTerm()

◆ dConvecDUComp()

◆ dStrongPressureDPressure()

◆ dStrongViscDUCompLaplace()

◆ dStrongViscDUCompTraction()

◆ dTauDUComp()

◆ dTimeDerivativeDUComp()

◆ dWeakPressureDPressure()

◆ dWeakViscDUCompLaplace()

◆ dWeakViscDUCompTraction()

◆ gravityTerm()

◆ strongPressureTerm()

◆ strongViscousTermLaplace()

◆ strongViscousTermTraction()

◆ tau()

◆ tauNodal()

◆ timeDerivativeTerm()

◆ weakPressureTerm()

◆ weakViscousTermLaplace()

◆ weakViscousTermTraction()

Member Data Documentation

◆ _alpha

◆ _component

◆ _convective_term

◆ _d_u_vel_dot_du

◆ _d_v_vel_dot_dv

◆ _d_w_vel_dot_dw

◆ _ffn

◆ _grad_p

◆ _grad_u_vel

◆ _grad_v_vel

◆ _grad_w_vel

◆ _gravity

◆ _integrate_p_by_parts

◆ _laplace

◆ _mu

◆ _p

◆ _p_var_number

◆ _rho

◆ _second_phi

◆ _second_u_vel

◆ _second_v_vel

◆ _second_w_vel

◆ _supg

◆ _transient_term

◆ _u_vel

◆ _u_vel_dot

◆ _u_vel_var_number

◆ _v_vel

◆ _v_vel_dot

◆ _v_vel_var_number

◆ _w_vel

◆ _w_vel_dot

◆ _w_vel_var_number