\[ (k(u) (u-u_\infty), v) \]

#include <NLBoundaryConvectiveHeatFluxIntegrator.h>

Inheritance diagram for Moose::MFEM::NLBoundaryConvectiveHeatFluxIntegrator:

Public Member Functions
	NLBoundaryConvectiveHeatFluxIntegrator (mfem::Coefficient &k, mfem::Coefficient &dk_du, mfem::Coefficient &duinf_du, mfem::Coefficient &gf_offset, mfem::Coefficient &gf)
	Construct the nonlinear convective heat flux residual and Jacobian action. More...

virtual void	AssembleElementVector (const mfem::FiniteElement &el, mfem::ElementTransformation &Tr, const mfem::Vector &elfun, mfem::Vector &elvect) override

virtual void	AssembleElementGrad (const mfem::FiniteElement &el, mfem::ElementTransformation &Tr, const mfem::Vector &elfun, mfem::DenseMatrix &elmat) override

Protected Attributes
mfem::SumCoefficient	_shifted_gf_coef
	Coefficient for the shifted state, u - u_inf. More...

mfem::ProductCoefficient	_k_uinf_coef
	Coefficient for the inwards flux contribution, k(u) * u_inf. More...

mfem::ProductCoefficient	_net_flux_du_coef
	Coefficient for the dk/du contribution, dk/du * (u - u_inf). More...

mfem::ProductCoefficient	_k_duinf_du_coef
	Coefficient for the d u_inf / du contribution, k(u) * d u_inf / du. More...

mfem::ProductCoefficient	_duinf_du_flux_coef
	Negated coefficient for the d u_inf / du Jacobian term, -k(u) * d u_inf / du. More...

mfem::BoundaryLFIntegrator	_inwards_flux
	Boundary load integrator for the inwards flux term. More...

mfem::BoundaryMassIntegrator	_outwards_flux
	Boundary mass integrator for the outwards flux term. More...

mfem::BoundaryMassIntegrator	_jacobian_k_component
	Jacobian contribution from the h(T) term. More...

mfem::BoundaryMassIntegrator	_jacobian_dk_du_component
	Jacobian contribution from the dh/dT term. More...

mfem::BoundaryMassIntegrator	_jacobian_duinf_du_component
	Jacobian contribution from the dT_inf/dT term. More...

mfem::SumIntegrator	_jacobian_action {0}
	Sum of the boundary mass integrators forming the Jacobian action. More...

Detailed Description

\[ (k(u) (u-u_\infty), v) \]

Definition at line 26 of file NLBoundaryConvectiveHeatFluxIntegrator.h.

Constructor & Destructor Documentation

◆ NLBoundaryConvectiveHeatFluxIntegrator()

Moose::MFEM::NLBoundaryConvectiveHeatFluxIntegrator::NLBoundaryConvectiveHeatFluxIntegrator	(	mfem::Coefficient &	k,
		mfem::Coefficient &	dk_du,
		mfem::Coefficient &	duinf_du,
		mfem::Coefficient &	gf_offset,
		mfem::Coefficient &	gf
	)

Construct the nonlinear convective heat flux residual and Jacobian action.

Definition at line 16 of file NLBoundaryConvectiveHeatFluxIntegrator.C.

   : _shifted_gf_coef(gf, gf_offset, 1.0, -1.0),  // u-u_inf
     _k_uinf_coef(k, gf_offset),                  // k(u)*u_inf
     _net_flux_du_coef(dk_du, _shifted_gf_coef),  // dk/du*(u-u_inf)
     _k_duinf_du_coef(k, duinf_du),               // k(u) * du_inf/du
     _duinf_du_flux_coef(-1.0, _k_duinf_du_coef), // -k(u) * du_inf/du
     _inwards_flux(_k_uinf_coef),
     _outwards_flux(k),
     _jacobian_k_component(k),
     _jacobian_dk_du_component(_net_flux_du_coef),
     _jacobian_duinf_du_component(_duinf_du_flux_coef)
 {
   _jacobian_action.AddIntegrator(&_jacobian_k_component);
   _jacobian_action.AddIntegrator(&_jacobian_dk_du_component);
   _jacobian_action.AddIntegrator(&_jacobian_duinf_du_component);
 }

Member Function Documentation

◆ AssembleElementGrad()

void Moose::MFEM::NLBoundaryConvectiveHeatFluxIntegrator::AssembleElementGrad	(	const mfem::FiniteElement &	el,
		mfem::ElementTransformation &	Tr,
		const mfem::Vector &	elfun,
		mfem::DenseMatrix &	elmat
	)

overridevirtual

Definition at line 51 of file NLBoundaryConvectiveHeatFluxIntegrator.C.

 {
   // (k (1 - du_inf/du) + dk_du (u-u_inf)) phi, v);
   _jacobian_action.AssembleElementGrad(el, Tr, elfun, elmat);
 }

◆ AssembleElementVector()

void Moose::MFEM::NLBoundaryConvectiveHeatFluxIntegrator::AssembleElementVector	(	const mfem::FiniteElement &	el,
		mfem::ElementTransformation &	Tr,
		const mfem::Vector &	elfun,
		mfem::Vector &	elvect
	)

overridevirtual

Definition at line 39 of file NLBoundaryConvectiveHeatFluxIntegrator.C.

 {
   mfem::Vector ext_flux_elvect; // vector to store flux coming from outside of domain
   _outwards_flux.AssembleElementVector(el, Tr, elfun, elvect);   // (k(u)*u, v)
   _inwards_flux.AssembleRHSElementVect(el, Tr, ext_flux_elvect); // (k(u)*u_inf, v)
   elvect -= ext_flux_elvect;                                     // (k(u)*(u-u_inf),v)
 }