This class imposes a velocity direction component as a Dirichlet condition on the appropriate momentum equation. More...

#include <NSImposedVelocityDirectionBC.h>

Inheritance diagram for NSImposedVelocityDirectionBC:

[legend]

Public Member Functions
	NSImposedVelocityDirectionBC (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual ()

Protected Attributes
const VariableValue &	_rho

const VariableValue &	_u_vel

const VariableValue &	_v_vel

const VariableValue &	_w_vel

Real	_desired_unit_velocity_component

Detailed Description

This class imposes a velocity direction component as a Dirichlet condition on the appropriate momentum equation.

For example, in the x-direction, the residual equation becomes:

u1/|u| - u1_hat_desired = 0

or

u1 - u1_hat_desired*|u| = 0

or

rho*u1 - rho*u1_hat_desired*|u| = 0

where: u1 = the x-momentum component |u| = velocity magnitude u1_hat_desired = The desired velocity direction, \( \in (0,1) \)

Definition at line 41 of file NSImposedVelocityDirectionBC.h.

Constructor & Destructor Documentation

◆ NSImposedVelocityDirectionBC()

NSImposedVelocityDirectionBC::NSImposedVelocityDirectionBC ( const InputParameters & parameters )

Definition at line 41 of file NSImposedVelocityDirectionBC.C.

   : NodalBC(parameters),
     _rho(coupledValue(NS::density)),
     _u_vel(coupledValue(NS::velocity_x)),
     _v_vel(_mesh.dimension() == 2 ? coupledValue(NS::velocity_y) : _zero),
     _w_vel(_mesh.dimension() == 3 ? coupledValue(NS::velocity_z) : _zero),
     _desired_unit_velocity_component(getParam<Real>("desired_unit_velocity_component"))
 {
 }

Member Function Documentation

◆ computeQpResidual()

Real NSImposedVelocityDirectionBC::computeQpResidual ( )

protectedvirtual

Definition at line 52 of file NSImposedVelocityDirectionBC.C.

 {
   // The velocity vector
   RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
  
   // Specify desired velocity component
   return _u[_qp] - _rho[_qp] * _desired_unit_velocity_component * vel.norm();
 }