This class implements the "No BC" boundary condition based on the "traction" form of the viscous stress tensor. More...

#include <INSMomentumNoBCBCTractionForm.h>

Inheritance diagram for INSMomentumNoBCBCTractionForm:

Public Member Functions
	INSMomentumNoBCBCTractionForm (const InputParameters &parameters)

virtual	~INSMomentumNoBCBCTractionForm ()

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

virtual Real	computeQpOffDiagJacobian (unsigned jvar)

Protected Attributes
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

unsigned	_u_vel_var_number

unsigned	_v_vel_var_number

unsigned	_w_vel_var_number

unsigned	_p_var_number

RealVectorValue	_gravity

unsigned	_component

bool	_integrate_p_by_parts

const MaterialProperty< Real > &	_mu

const MaterialProperty< Real > &	_rho

Detailed Description

This class implements the "No BC" boundary condition based on the "traction" form of the viscous stress tensor.

Definition at line 24 of file INSMomentumNoBCBCTractionForm.h.

Constructor & Destructor Documentation

◆ INSMomentumNoBCBCTractionForm()

INSMomentumNoBCBCTractionForm::INSMomentumNoBCBCTractionForm ( const InputParameters & parameters )

Definition at line 26 of file INSMomentumNoBCBCTractionForm.C.

   : INSMomentumNoBCBCBase(parameters)
 {
 }

◆ ~INSMomentumNoBCBCTractionForm()

virtual INSMomentumNoBCBCTractionForm::~INSMomentumNoBCBCTractionForm ( )

inlinevirtual

Definition at line 29 of file INSMomentumNoBCBCTractionForm.h.

29 {}

Member Function Documentation

◆ computeQpJacobian()

Real INSMomentumNoBCBCTractionForm::computeQpJacobian ( )

protectedvirtual

Definition at line 69 of file INSMomentumNoBCBCTractionForm.C.

 {
   // The extra contribution comes from the "2" on the diagonal of the viscous stress tensor
   return -_mu[_qp] * (_grad_phi[_j][_qp] * _normals[_qp] +
                       _grad_phi[_j][_qp](_component) * _normals[_qp](_component)) *
          _test[_i][_qp];
 }

◆ computeQpOffDiagJacobian()

Real INSMomentumNoBCBCTractionForm::computeQpOffDiagJacobian ( unsigned jvar )

protectedvirtual

Definition at line 78 of file INSMomentumNoBCBCTractionForm.C.

 {
   if (jvar == _u_vel_var_number)
     return -_mu[_qp] * _grad_phi[_j][_qp](_component) * _normals[_qp](0) * _test[_i][_qp];
  
   else if (jvar == _v_vel_var_number)
     return -_mu[_qp] * _grad_phi[_j][_qp](_component) * _normals[_qp](1) * _test[_i][_qp];
  
   else if (jvar == _w_vel_var_number)
     return -_mu[_qp] * _grad_phi[_j][_qp](_component) * _normals[_qp](2) * _test[_i][_qp];
  
   else if (jvar == _p_var_number)
   {
     if (_integrate_p_by_parts)
       return _phi[_j][_qp] * _normals[_qp](_component) * _test[_i][_qp];
     else
       return 0.;
   }
  
   else
     return 0.;
 }

◆ computeQpResidual()

Real INSMomentumNoBCBCTractionForm::computeQpResidual ( )

protectedvirtual

Definition at line 32 of file INSMomentumNoBCBCTractionForm.C.

 {
   // Compute n . sigma . v, where n is unit normal and v is the test function.
   RealTensorValue sigma;
  
   // First row
   sigma(0, 0) = 2. * _mu[_qp] * _grad_u_vel[_qp](0);
   sigma(0, 1) = _mu[_qp] * (_grad_u_vel[_qp](1) + _grad_v_vel[_qp](0));
   sigma(0, 2) = _mu[_qp] * (_grad_u_vel[_qp](2) + _grad_w_vel[_qp](0));
  
   // Second row
   sigma(1, 0) = _mu[_qp] * (_grad_v_vel[_qp](0) + _grad_u_vel[_qp](1));
   sigma(1, 1) = 2. * _mu[_qp] * _grad_v_vel[_qp](1);
   sigma(1, 2) = _mu[_qp] * (_grad_v_vel[_qp](2) + _grad_w_vel[_qp](1));
  
   // Third row
   sigma(2, 0) = _mu[_qp] * (_grad_w_vel[_qp](0) + _grad_u_vel[_qp](2));
   sigma(2, 1) = _mu[_qp] * (_grad_w_vel[_qp](1) + _grad_v_vel[_qp](2));
   sigma(2, 2) = 2. * _mu[_qp] * _grad_w_vel[_qp](2);
  
   // If the pressure term is integrated by parts, it is part of the
   // no-BC-BC, otherwise, it is not.
   if (_integrate_p_by_parts)
   {
     sigma(0, 0) -= _p[_qp];
     sigma(1, 1) -= _p[_qp];
     sigma(2, 2) -= _p[_qp];
   }
  
   // Set up test function
   RealVectorValue test;
   test(_component) = _test[_i][_qp];
  
   return -_normals[_qp] * (sigma * test);
 }

Member Data Documentation

◆ _component

unsigned INSMomentumNoBCBCBase::_component

protectedinherited

Definition at line 56 of file INSMomentumNoBCBCBase.h.

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

◆ _grad_u_vel

const VariableGradient& INSMomentumNoBCBCBase::_grad_u_vel

protectedinherited

Definition at line 45 of file INSMomentumNoBCBCBase.h.

Referenced by computeQpResidual().

◆ _grad_v_vel

const VariableGradient& INSMomentumNoBCBCBase::_grad_v_vel

protectedinherited

Definition at line 46 of file INSMomentumNoBCBCBase.h.

Referenced by computeQpResidual().

◆ _grad_w_vel

const VariableGradient& INSMomentumNoBCBCBase::_grad_w_vel

protectedinherited

Definition at line 47 of file INSMomentumNoBCBCBase.h.

Referenced by computeQpResidual().

◆ _gravity

RealVectorValue INSMomentumNoBCBCBase::_gravity

protectedinherited

Definition at line 55 of file INSMomentumNoBCBCBase.h.

◆ _integrate_p_by_parts

bool INSMomentumNoBCBCBase::_integrate_p_by_parts

protectedinherited

Definition at line 57 of file INSMomentumNoBCBCBase.h.

Referenced by INSMomentumNoBCBCLaplaceForm::computeQpOffDiagJacobian(), computeQpOffDiagJacobian(), computeQpResidual(), and INSMomentumNoBCBCLaplaceForm::computeQpResidual().

◆ _mu

const MaterialProperty<Real>& INSMomentumNoBCBCBase::_mu

protectedinherited

Definition at line 59 of file INSMomentumNoBCBCBase.h.

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

◆ _p

const VariableValue& INSMomentumNoBCBCBase::_p

protectedinherited

Definition at line 42 of file INSMomentumNoBCBCBase.h.

Referenced by INSMomentumNoBCBCLaplaceForm::computeQpResidual(), and computeQpResidual().

◆ _p_var_number

unsigned INSMomentumNoBCBCBase::_p_var_number

protectedinherited

Definition at line 53 of file INSMomentumNoBCBCBase.h.

Referenced by INSMomentumNoBCBCLaplaceForm::computeQpOffDiagJacobian(), and computeQpOffDiagJacobian().

◆ _rho

const MaterialProperty<Real>& INSMomentumNoBCBCBase::_rho

protectedinherited

Definition at line 60 of file INSMomentumNoBCBCBase.h.

◆ _u_vel

const VariableValue& INSMomentumNoBCBCBase::_u_vel

protectedinherited

Definition at line 39 of file INSMomentumNoBCBCBase.h.

◆ _u_vel_var_number

unsigned INSMomentumNoBCBCBase::_u_vel_var_number

protectedinherited

Definition at line 50 of file INSMomentumNoBCBCBase.h.

Referenced by computeQpOffDiagJacobian().

◆ _v_vel

const VariableValue& INSMomentumNoBCBCBase::_v_vel

protectedinherited

Definition at line 40 of file INSMomentumNoBCBCBase.h.

◆ _v_vel_var_number

unsigned INSMomentumNoBCBCBase::_v_vel_var_number

protectedinherited

Definition at line 51 of file INSMomentumNoBCBCBase.h.

Referenced by computeQpOffDiagJacobian().

◆ _w_vel

const VariableValue& INSMomentumNoBCBCBase::_w_vel

protectedinherited

Definition at line 41 of file INSMomentumNoBCBCBase.h.

◆ _w_vel_var_number

unsigned INSMomentumNoBCBCBase::_w_vel_var_number

protectedinherited

Definition at line 52 of file INSMomentumNoBCBCBase.h.

Referenced by computeQpOffDiagJacobian().

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

navier_stokes/include/bcs/INSMomentumNoBCBCTractionForm.h
navier_stokes/src/bcs/INSMomentumNoBCBCTractionForm.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ INSMomentumNoBCBCTractionForm()

◆ ~INSMomentumNoBCBCTractionForm()

Member Function Documentation

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

Member Data Documentation

◆ _component

◆ _grad_u_vel

◆ _grad_v_vel

◆ _grad_w_vel

◆ _gravity

◆ _integrate_p_by_parts

◆ _mu

◆ _p

◆ _p_var_number

◆ _rho

◆ _u_vel

◆ _u_vel_var_number

◆ _v_vel

◆ _v_vel_var_number

◆ _w_vel

◆ _w_vel_var_number