#include <XFEMPhaseTransitionMovingInterfaceVelocity.h>

Inheritance diagram for XFEMPhaseTransitionMovingInterfaceVelocity:

Public Member Functions
	XFEMPhaseTransitionMovingInterfaceVelocity (const InputParameters &parameters)

virtual	~XFEMPhaseTransitionMovingInterfaceVelocity ()

virtual Real	computeMovingInterfaceVelocity (unsigned int point_id) const override
	Compute the interface velocity for a point. More...

virtual void	initialize () override

unsigned int	numberPoints () const
	Compute total number of points that are used to define an interface. More...

Protected Attributes
Real	_diffusivity_at_positive_level_set
	Diffusivity in the positive level set region. More...

Real	_diffusivity_at_negative_level_set
	Diffusivity in the negative level set region. More...

Real	_equilibrium_concentration_jump
	Jump of the equilibrium concentrations at phase boundary. More...

const PointValueAtXFEMInterface *	_value_at_interface_uo
	Pointer to PointValueAtXFEMInterface object. More...

Detailed Description

Definition at line 19 of file XFEMPhaseTransitionMovingInterfaceVelocity.h.

Constructor & Destructor Documentation

◆ XFEMPhaseTransitionMovingInterfaceVelocity()

XFEMPhaseTransitionMovingInterfaceVelocity::XFEMPhaseTransitionMovingInterfaceVelocity ( const InputParameters & parameters )

Definition at line 30 of file XFEMPhaseTransitionMovingInterfaceVelocity.C.

   : XFEMMovingInterfaceVelocityBase(parameters),
     _diffusivity_at_positive_level_set(getParam<Real>("diffusivity_at_positive_level_set")),
     _diffusivity_at_negative_level_set(getParam<Real>("diffusivity_at_negative_level_set")),
     _equilibrium_concentration_jump(getParam<Real>("equilibrium_concentration_jump"))
 {
 }

◆ ~XFEMPhaseTransitionMovingInterfaceVelocity()

virtual XFEMPhaseTransitionMovingInterfaceVelocity::~XFEMPhaseTransitionMovingInterfaceVelocity ( )

inlinevirtual

Definition at line 23 of file XFEMPhaseTransitionMovingInterfaceVelocity.h.

23 {}

Member Function Documentation

◆ computeMovingInterfaceVelocity()

Real XFEMPhaseTransitionMovingInterfaceVelocity::computeMovingInterfaceVelocity ( unsigned int point_id ) const

overridevirtual

Compute the interface velocity for a point.

Parameters

point_id Point ID

Returns: Real Interface velocity

Implements XFEMMovingInterfaceVelocityBase.

Definition at line 40 of file XFEMPhaseTransitionMovingInterfaceVelocity.C.

 {
   Real value_positive = _value_at_interface_uo->getValueAtPositiveLevelSet()[point_id];
   Real value_negative = _value_at_interface_uo->getValueAtNegativeLevelSet()[point_id];
   RealVectorValue grad_positive = _value_at_interface_uo->getGradientAtPositiveLevelSet()[point_id];
   RealVectorValue grad_negative = _value_at_interface_uo->getGradientAtNegativeLevelSet()[point_id];
  
   // Current implementation only supports the case that the interface is moving horizontally
   return std::abs((_diffusivity_at_positive_level_set * grad_positive(0) -
                    _diffusivity_at_negative_level_set * grad_negative(0)) /
                   (value_positive - value_negative + _equilibrium_concentration_jump));
 }

◆ initialize()

void XFEMMovingInterfaceVelocityBase::initialize ( )

overridevirtualinherited

Definition at line 29 of file XFEMMovingInterfaceVelocityBase.C.

 {
   const UserObject * uo =
       &(_fe_problem.getUserObjectBase(getParam<UserObjectName>("value_at_interface_uo")));
  
   if (dynamic_cast<const PointValueAtXFEMInterface *>(uo) == nullptr)
     mooseError("UserObject casting to PointValueAtXFEMInterface in XFEMMovingInterfaceVelocity");
  
   _value_at_interface_uo = dynamic_cast<const PointValueAtXFEMInterface *>(uo);
 }