INSProjection Class Reference

This class computes the "projection" part of the "split" method for solving incompressible Navier-Stokes. More...

#include <INSProjection.h>

Inheritance diagram for INSProjection:

Public Member Functions
	INSProjection (const InputParameters &parameters)
virtual	~INSProjection ()

Protected Member Functions
virtual Real	computeQpResidual ()
virtual Real	computeQpJacobian ()
virtual Real	computeQpOffDiagJacobian (unsigned jvar)

Protected Attributes
const VariableValue &	_a1
const VariableValue &	_a2
const VariableValue &	_a3
const VariableGradient &	_grad_p
unsigned	_a1_var_number
unsigned	_a2_var_number
unsigned	_a3_var_number
unsigned	_p_var_number
unsigned	_component
const MaterialProperty< Real > &	_rho

Detailed Description

This class computes the "projection" part of the "split" method for solving incompressible Navier-Stokes.

This is a time-varying equation for u that is coupled to both the acceleration "a" and the pressue.

Do not use, USE INSChorinCorrector and related classes instead.

Definition at line 28 of file INSProjection.h.

Constructor & Destructor Documentation

INSProjection()

INSProjection::INSProjection

(

const InputParameters &

parameters

)

Definition at line 40 of file INSProjection.C.

  : Kernel(parameters),

    // Coupled variables
    _a1(coupledValue("a1")),
    _a2(_mesh.dimension() >= 2 ? coupledValue("a2") : _zero),
    _a3(_mesh.dimension() == 3 ? coupledValue("a3") : _zero),

    // Gradients
    _grad_p(coupledGradient("p")),

    // Variable numberings
    _a1_var_number(coupled("a1")),
    _a2_var_number(_mesh.dimension() >= 2 ? coupled("a2") : libMesh::invalid_uint),
    _a3_var_number(_mesh.dimension() == 3 ? coupled("a3") : libMesh::invalid_uint),
    _p_var_number(coupled("p")),

    // Required parameters
    _component(getParam<unsigned>("component")),

    // Material properties
    _rho(getMaterialProperty<Real>("rho_name"))
{
}

~INSProjection()

virtual INSProjection::~INSProjection ( )

inlinevirtual

Definition at line 33 of file INSProjection.h.

{}

Member Function Documentation

computeQpJacobian()

Real INSProjection::computeQpJacobian ( )

protectedvirtual

Definition at line 86 of file INSProjection.C.

{
  // There will be a diagonal component from the time derivative term...
  return 0.;
}

computeQpOffDiagJacobian()

Real INSProjection::computeQpOffDiagJacobian ( unsigned  jvar )

protectedvirtual

Definition at line 93 of file INSProjection.C.

{
  if (((jvar == _a1_var_number) && (_component == 0)) ||
      ((jvar == _a2_var_number) && (_component == 1)) ||
      ((jvar == _a3_var_number) && (_component == 2)))
  {
    // The symmetric term's Jacobian is only non-zero when the
    // component of 'a' being differentiated is the same as _component.
    return -_phi[_j][_qp] * _test[_i][_qp];
  }

  else if (jvar == _p_var_number)
  {
    return (1. / _rho[_qp]) * (_grad_phi[_j][_qp](_component) * _test[_i][_qp]);
  }

  else
    return 0;
}

computeQpResidual()

Real INSProjection::computeQpResidual ( )

protectedvirtual

Definition at line 66 of file INSProjection.C.

{
  // Vector object for a
  RealVectorValue a(_a1[_qp], _a2[_qp], _a3[_qp]);

  // Vector object for test function (only the component'th entry is non-zero)
  RealVectorValue test;
  test(_component) = _test[_i][_qp];

  // "Symmetric" part, -a.test
  Real symmetric_part = -a(_component) * _test[_i][_qp];

  // The pressure part, (1/_rho[_qp]) * (grad(p).v)
  Real pressure_part = (1. / _rho[_qp]) * (_grad_p[_qp] * test);

  // Return the result
  return symmetric_part + pressure_part;
}

Member Data Documentation

_a1

const VariableValue& INSProjection::_a1

protected

Definition at line 41 of file INSProjection.h.

Referenced by computeQpResidual().

_a1_var_number

unsigned INSProjection::_a1_var_number

protected

Definition at line 49 of file INSProjection.h.

Referenced by computeQpOffDiagJacobian().

_a2

const VariableValue& INSProjection::_a2

protected

Definition at line 42 of file INSProjection.h.

Referenced by computeQpResidual().

_a2_var_number

unsigned INSProjection::_a2_var_number

protected

Definition at line 50 of file INSProjection.h.

Referenced by computeQpOffDiagJacobian().

_a3

const VariableValue& INSProjection::_a3

protected

Definition at line 43 of file INSProjection.h.

Referenced by computeQpResidual().

_a3_var_number

unsigned INSProjection::_a3_var_number

protected

Definition at line 51 of file INSProjection.h.

Referenced by computeQpOffDiagJacobian().

_component

unsigned INSProjection::_component

protected

Definition at line 55 of file INSProjection.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

_grad_p

const VariableGradient& INSProjection::_grad_p

protected

Definition at line 46 of file INSProjection.h.

Referenced by computeQpResidual().

_p_var_number

unsigned INSProjection::_p_var_number

protected

Definition at line 52 of file INSProjection.h.

Referenced by computeQpOffDiagJacobian().

_rho

const MaterialProperty<Real>& INSProjection::_rho

protected

Definition at line 58 of file INSProjection.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

---

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

• navier_stokes/include/kernels/INSProjection.h
• navier_stokes/src/kernels/INSProjection.C