Computes the force per unit length due to form loss, provided a form loss coefficient per unit length function. More...

#include <OneD3EqnMomentumFormLoss.h>

Inheritance diagram for OneD3EqnMomentumFormLoss:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	OneD3EqnMomentumFormLoss (const InputParameters &parameters)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)

void	validateDerivativeMaterialPropertyBase (const std::string &base)

const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

MaterialProperty< U > &	declarePropertyDerivativeTHM (const std::string &base, const std::string &var_name, const unsigned int i=0)
	Method for declaring derivative material properties. More...

const MaterialProperty< U > &	getMaterialPropertyDerivativeTHM (const std::string &base, const std::string &var_name, const unsigned int i=0)
	Method for retrieving derivative material properties. More...

const MaterialProperty< U > &	getMaterialPropertyDerivativeTHMPhase (const std::string &base, bool property_is_liquid, const std::string &var_name, bool var_is_liquid, const unsigned int i=0)
	Method for retrieving derivative material properties corresponding to phase-dependent derivatives. More...

Protected Attributes
const VariableValue &	_A
	area More...

const MaterialProperty< Real > &	_rho
	density More...

const MaterialProperty< Real > &	_drho_darhoA

const MaterialProperty< Real > &	_vel
	velocity More...

const MaterialProperty< Real > &	_dvel_darhoA

const MaterialProperty< Real > &	_dvel_darhouA

const MaterialProperty< Real > &	_K_prime
	form loss coefficient per unit length function More...

unsigned int	_arhoA_var_number

unsigned int	_arhouA_var_number

unsigned int	_arhoEA_var_number

Detailed Description

Computes the force per unit length due to form loss, provided a form loss coefficient per unit length function.

See RELAP-7 Theory Manual, pg. 72, Equation (239) {eq:form_loss_force_2phase}

Definition at line 23 of file OneD3EqnMomentumFormLoss.h.

Constructor & Destructor Documentation

◆ OneD3EqnMomentumFormLoss()

OneD3EqnMomentumFormLoss::OneD3EqnMomentumFormLoss ( const InputParameters & parameters )

Definition at line 33 of file OneD3EqnMomentumFormLoss.C.

   : DerivativeMaterialInterfaceTHM<Kernel>(parameters),
 
     _A(coupledValue("A")),
 
     _rho(getMaterialProperty<Real>("rho")),
     _drho_darhoA(getMaterialPropertyDerivativeTHM<Real>("rho", "arhoA")),
 
     _vel(getMaterialProperty<Real>("vel")),
     _dvel_darhoA(getMaterialPropertyDerivativeTHM<Real>("vel", "arhoA")),
     _dvel_darhouA(getMaterialPropertyDerivativeTHM<Real>("vel", "arhouA")),
 
     _K_prime(getMaterialProperty<Real>("K_prime")),
 
     _arhoA_var_number(coupled("arhoA")),
     _arhouA_var_number(coupled("arhouA")),
     _arhoEA_var_number(coupled("arhoEA"))
 {
 }

Member Function Documentation

◆ computeQpJacobian()

Real OneD3EqnMomentumFormLoss::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 61 of file OneD3EqnMomentumFormLoss.C.

 {
   return computeQpOffDiagJacobian(_var.number());
 }

◆ computeQpOffDiagJacobian()

Real OneD3EqnMomentumFormLoss::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 67 of file OneD3EqnMomentumFormLoss.C.

Referenced by computeQpJacobian().

 {
   if (jvar == _arhoA_var_number)
   {
     const Real drhou2_darhoA = _drho_darhoA[_qp] * _vel[_qp] * std::abs(_vel[_qp]) +
                                2.0 * _rho[_qp] * std::abs(_vel[_qp]) * _dvel_darhoA[_qp];
 
     return _K_prime[_qp] * 0.5 * drhou2_darhoA * _A[_qp] * _phi[_j][_qp] * _test[_i][_qp];
   }
   else if (jvar == _arhouA_var_number)
   {
     const Real drhou2_darhouA = 2.0 * _rho[_qp] * std::abs(_vel[_qp]) * _dvel_darhouA[_qp];
 
     return _K_prime[_qp] * 0.5 * drhou2_darhouA * _A[_qp] * _phi[_j][_qp] * _test[_i][_qp];
   }
   else if (jvar == _arhoEA_var_number)
     return 0.;
   else
     return 0;
 }

◆ computeQpResidual()

Real OneD3EqnMomentumFormLoss::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 54 of file OneD3EqnMomentumFormLoss.C.

 {
   return _K_prime[_qp] * 0.5 * _rho[_qp] * _vel[_qp] * std::abs(_vel[_qp]) * _A[_qp] *
          _test[_i][_qp];
 }

◆ declarePropertyDerivativeTHM()

MaterialProperty< U > & DerivativeMaterialInterfaceTHM< Kernel >::declarePropertyDerivativeTHM	(	const std::string &	base,
		const std::string &	var_name,
		const unsigned int	i = `0`
	)

protectedinherited

Method for declaring derivative material properties.

Template Parameters

U	The material property type

Parameters

base	The name of the property to take the derivative of
var_name	The name of the coupled variable vector containing the variable to take the derivative with respect to
i	The index of the variable to take the derivative with respect to within the coupled variable vector

Definition at line 88 of file DerivativeMaterialInterfaceTHM.h.

 {
   return this->template declarePropertyDerivative<U>(base, this->coupledName(var_name, i));
 }

◆ getMaterialPropertyDerivativeTHM()

const MaterialProperty< U > & DerivativeMaterialInterfaceTHM< Kernel >::getMaterialPropertyDerivativeTHM	(	const std::string &	base,
		const std::string &	var_name,
		const unsigned int	i = `0`
	)

protectedinherited

Method for retrieving derivative material properties.

Template Parameters

U	The material property type

Parameters

base	The name of the property to take the derivative of
var_name	The name of the coupled variable vector containing the variable to take the derivative with respect to
i	The index of the variable to take the derivative with respect to within the coupled variable vector

Definition at line 98 of file DerivativeMaterialInterfaceTHM.h.

 {
   // get the base property name
   const std::string prop_name = this->getMaterialPropertyName(base);
 
   // get the name of the variable which derivative is respect to
   const std::string der_var_name = this->coupledName(var_name, i);
 
   return this->template getMaterialPropertyByName<U>(
       this->derivativePropertyNameFirst(prop_name, der_var_name));
 }

◆ getMaterialPropertyDerivativeTHMPhase()

const MaterialProperty< U > & DerivativeMaterialInterfaceTHM< Kernel >::getMaterialPropertyDerivativeTHMPhase	(	const std::string &	base,
		bool	property_is_liquid,
		const std::string &	var_name,
		bool	var_is_liquid,
		const unsigned int	i = `0`
	)

protectedinherited

Method for retrieving derivative material properties corresponding to phase-dependent derivatives.

The intended behavior of this function is as follows. If the material property and derivative variable correspond to the same phase, then the material property is retrieved as usual. Else, a zero-valued material property is returned, as it is assumed that the material property depends only upon the variables corresponding to its phase.

Template Parameters

U	The material property type

Parameters

base	The name of the property to take the derivative of
property_is_liquid	Does the property correspond to the liquid phase?
var_name	The name of the coupled variable vector containing the variable to take the derivative with respect to
var_is_liquid	Does the derivative variable correspond to the liquid phase?
i	The index of the variable to take the derivative with respect to within the coupled variable vector

Definition at line 115 of file DerivativeMaterialInterfaceTHM.h.

 {
   if (property_is_liquid == var_is_liquid)
   {
     return getMaterialPropertyDerivativeTHM<U>(base, var_name, i);
   }
   else
   {
     const std::string prop_name = this->getMaterialPropertyName(base);
     const std::string der_var_name = this->coupledName(var_name, i);
     const std::string der_prop_name = this->derivativePropertyNameFirst(prop_name, der_var_name);
     return this->template getZeroMaterialProperty<U>(der_prop_name);
   }
 }

◆ validParams()

InputParameters OneD3EqnMomentumFormLoss::validParams ( )

static

Definition at line 16 of file OneD3EqnMomentumFormLoss.C.

 {
   InputParameters params = Kernel::validParams();
 
   params.addRequiredCoupledVar("arhoA", "Solution variable alpha*rho*A");
   params.addRequiredCoupledVar("arhouA", "Solution variable alpha*rho*u*A");
   params.addRequiredCoupledVar("arhoEA", "Solution variable alpha*rho*E*A");
 
   params.addRequiredCoupledVar("A", "Cross-sectional area");
   params.addRequiredParam<MaterialPropertyName>("rho", "Density property");
   params.addRequiredParam<MaterialPropertyName>("vel", "Velocity property");
   params.addClassDescription(
       "Computes a form loss term for the momentum equation for 1-phase flow");
 
   return params;
 }

Member Data Documentation

◆ _A

const VariableValue& OneD3EqnMomentumFormLoss::_A

protected

area

Definition at line 34 of file OneD3EqnMomentumFormLoss.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _arhoA_var_number

unsigned int OneD3EqnMomentumFormLoss::_arhoA_var_number

protected

Definition at line 48 of file OneD3EqnMomentumFormLoss.h.

Referenced by computeQpOffDiagJacobian().

◆ _arhoEA_var_number

unsigned int OneD3EqnMomentumFormLoss::_arhoEA_var_number

protected

Definition at line 50 of file OneD3EqnMomentumFormLoss.h.

Referenced by computeQpOffDiagJacobian().

◆ _arhouA_var_number

unsigned int OneD3EqnMomentumFormLoss::_arhouA_var_number

protected

Definition at line 49 of file OneD3EqnMomentumFormLoss.h.

Referenced by computeQpOffDiagJacobian().

◆ _drho_darhoA

const MaterialProperty<Real>& OneD3EqnMomentumFormLoss::_drho_darhoA

protected

Definition at line 38 of file OneD3EqnMomentumFormLoss.h.

Referenced by computeQpOffDiagJacobian().

◆ _dvel_darhoA

const MaterialProperty<Real>& OneD3EqnMomentumFormLoss::_dvel_darhoA

protected

Definition at line 42 of file OneD3EqnMomentumFormLoss.h.

Referenced by computeQpOffDiagJacobian().

◆ _dvel_darhouA

const MaterialProperty<Real>& OneD3EqnMomentumFormLoss::_dvel_darhouA

protected

Definition at line 43 of file OneD3EqnMomentumFormLoss.h.

Referenced by computeQpOffDiagJacobian().

◆ _K_prime

const MaterialProperty<Real>& OneD3EqnMomentumFormLoss::_K_prime

protected

form loss coefficient per unit length function

Definition at line 46 of file OneD3EqnMomentumFormLoss.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _rho

const MaterialProperty<Real>& OneD3EqnMomentumFormLoss::_rho

protected

density

Definition at line 37 of file OneD3EqnMomentumFormLoss.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _vel

const MaterialProperty<Real>& OneD3EqnMomentumFormLoss::_vel

protected

velocity

Definition at line 41 of file OneD3EqnMomentumFormLoss.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

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

thermal_hydraulics/include/kernels/OneD3EqnMomentumFormLoss.h
thermal_hydraulics/src/kernels/OneD3EqnMomentumFormLoss.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ OneD3EqnMomentumFormLoss()

Member Function Documentation

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ declarePropertyDerivativeTHM()

◆ getMaterialPropertyDerivativeTHM()

◆ getMaterialPropertyDerivativeTHMPhase()

◆ validParams()

Member Data Documentation

◆ _A

◆ _arhoA_var_number

◆ _arhoEA_var_number

◆ _arhouA_var_number

◆ _drho_darhoA

◆ _dvel_darhoA

◆ _dvel_darhouA

◆ _K_prime

◆ _rho

◆ _vel