XFEMSingleVariableConstraint Class Reference

#include <XFEMSingleVariableConstraint.h>

Inheritance diagram for XFEMSingleVariableConstraint:

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

Protected Member Functions
virtual void	reinitConstraintQuadrature (const ElementPairInfo &element_pair_info) override
virtual Real	computeQpResidual (Moose::DGResidualType type) override
virtual Real	computeQpJacobian (Moose::DGJacobianType type) override

Protected Attributes
Point	_interface_normal
	Vector normal to the internal interface. More...
Real	_alpha
	Stabilization parameter in Nitsche's formulation and penalty factor in the Penalty Method. More...
const Function &	_jump
	Change in variable value at the interface. More...
const Function &	_jump_flux
	Change in flux of variable value at the interface. More...
bool	_use_penalty
	Use penalty formulation. More...
std::shared_ptr< XFEM >	_xfem
	Pointer to the XFEM controller object. More...

Detailed Description

Definition at line 26 of file XFEMSingleVariableConstraint.h.

Constructor & Destructor Documentation

XFEMSingleVariableConstraint()

XFEMSingleVariableConstraint::XFEMSingleVariableConstraint

(

const InputParameters &

parameters

)

Definition at line 48 of file XFEMSingleVariableConstraint.C.

  : ElemElemConstraint(parameters),
    _alpha(getParam<Real>("alpha")),
    _jump(getFunction("jump")),
    _jump_flux(getFunction("jump_flux")),
    _use_penalty(getParam<bool>("use_penalty"))
{
  _xfem = std::dynamic_pointer_cast<XFEM>(_fe_problem.getXFEM());
  if (_xfem == nullptr)
    mooseError("Problem casting to XFEM in XFEMSingleVariableConstraint");
 
  const UserObject * uo =
      &(_fe_problem.getUserObjectBase(getParam<UserObjectName>("geometric_cut_userobject")));
 
  if (dynamic_cast<const GeometricCutUserObject *>(uo) == nullptr)
    mooseError("UserObject casting to GeometricCutUserObject in XFEMSingleVariableConstraint");
 
  _interface_id = _xfem->getGeometricCutID(dynamic_cast<const GeometricCutUserObject *>(uo));
}

~XFEMSingleVariableConstraint()

XFEMSingleVariableConstraint::~XFEMSingleVariableConstraint ( )

virtual

Definition at line 68 of file XFEMSingleVariableConstraint.C.

{}

Member Function Documentation

computeQpJacobian()

Real XFEMSingleVariableConstraint::computeQpJacobian ( Moose::DGJacobianType  type )

overrideprotectedvirtual

Definition at line 116 of file XFEMSingleVariableConstraint.C.

{
  Real r = 0;
 
  switch (type)
  {
    case Moose::ElementElement:
      if (!_use_penalty)
        r += -0.5 * _grad_phi[_j][_qp] * _interface_normal * _test[_i][_qp] -
             _phi[_j][_qp] * 0.5 * _grad_test[_i][_qp] * _interface_normal;
      r += _alpha * _phi[_j][_qp] * _test[_i][_qp];
      break;
 
    case Moose::ElementNeighbor:
      if (!_use_penalty)
        r += -0.5 * _grad_phi_neighbor[_j][_qp] * _interface_normal * _test[_i][_qp] +
             _phi_neighbor[_j][_qp] * 0.5 * _grad_test[_i][_qp] * _interface_normal;
      r -= _alpha * _phi_neighbor[_j][_qp] * _test[_i][_qp];
      break;
 
    case Moose::NeighborElement:
      if (!_use_penalty)
        r += 0.5 * _grad_phi[_j][_qp] * _interface_normal * _test_neighbor[_i][_qp] -
             _phi[_j][_qp] * 0.5 * _grad_test_neighbor[_i][_qp] * _interface_normal;
      r -= _alpha * _phi[_j][_qp] * _test_neighbor[_i][_qp];
      break;
 
    case Moose::NeighborNeighbor:
      if (!_use_penalty)
        r += 0.5 * _grad_phi_neighbor[_j][_qp] * _interface_normal * _test_neighbor[_i][_qp] +
             _phi_neighbor[_j][_qp] * 0.5 * _grad_test_neighbor[_i][_qp] * _interface_normal;
      r += _alpha * _phi_neighbor[_j][_qp] * _test_neighbor[_i][_qp];
      break;
  }
 
  return r;
}

computeQpResidual()

Real XFEMSingleVariableConstraint::computeQpResidual ( Moose::DGResidualType  type )

overrideprotectedvirtual

Definition at line 78 of file XFEMSingleVariableConstraint.C.

{
  Real r = 0;
 
  switch (type)
  {
    case Moose::Element:
      if (!_use_penalty)
      {
        r -= (0.5 * _grad_u[_qp] * _interface_normal +
              0.5 * _grad_u_neighbor[_qp] * _interface_normal) *
             _test[_i][_qp];
        r -= (_u[_qp] - _u_neighbor[_qp]) * 0.5 * _grad_test[_i][_qp] * _interface_normal;
        r += 0.5 * _grad_test[_i][_qp] * _interface_normal * _jump.value(_t, _u[_qp]);
      }
      r += 0.5 * _test[_i][_qp] * _jump_flux.value(_t, _u[_qp]);
      r += _alpha * (_u[_qp] - _u_neighbor[_qp] - _jump.value(_t, _u[_qp])) * _test[_i][_qp];
      break;
 
    case Moose::Neighbor:
      if (!_use_penalty)
      {
        r += (0.5 * _grad_u[_qp] * _interface_normal +
              0.5 * _grad_u_neighbor[_qp] * _interface_normal) *
             _test_neighbor[_i][_qp];
        r -= (_u[_qp] - _u_neighbor[_qp]) * 0.5 * _grad_test_neighbor[_i][_qp] * _interface_normal;
        r += 0.5 * _grad_test_neighbor[_i][_qp] * _interface_normal *
             _jump.value(_t, _u_neighbor[_qp]);
      }
      r += 0.5 * _test_neighbor[_i][_qp] * _jump_flux.value(_t, _u_neighbor[_qp]);
      r -= _alpha * (_u[_qp] - _u_neighbor[_qp] - _jump.value(_t, _u_neighbor[_qp])) *
           _test_neighbor[_i][_qp];
      break;
  }
  return r;
}

reinitConstraintQuadrature()

void XFEMSingleVariableConstraint::reinitConstraintQuadrature ( const ElementPairInfo &  element_pair_info )

overrideprotectedvirtual

Definition at line 71 of file XFEMSingleVariableConstraint.C.

{
  _interface_normal = element_pair_info._elem1_normal;
  ElemElemConstraint::reinitConstraintQuadrature(element_pair_info);
}

Member Data Documentation

_alpha

Real XFEMSingleVariableConstraint::_alpha

protected

Stabilization parameter in Nitsche's formulation and penalty factor in the Penalty Method.

Definition at line 44 of file XFEMSingleVariableConstraint.h.

Referenced by computeQpJacobian(), and computeQpResidual().

_interface_normal

Point XFEMSingleVariableConstraint::_interface_normal

protected

Vector normal to the internal interface.

Definition at line 40 of file XFEMSingleVariableConstraint.h.

Referenced by computeQpJacobian(), computeQpResidual(), and reinitConstraintQuadrature().

_jump

const Function& XFEMSingleVariableConstraint::_jump

protected

Change in variable value at the interface.

Definition at line 47 of file XFEMSingleVariableConstraint.h.

Referenced by computeQpResidual().

_jump_flux

const Function& XFEMSingleVariableConstraint::_jump_flux

protected

Change in flux of variable value at the interface.

Definition at line 50 of file XFEMSingleVariableConstraint.h.

Referenced by computeQpResidual().

_use_penalty

bool XFEMSingleVariableConstraint::_use_penalty

protected

Use penalty formulation.

Definition at line 53 of file XFEMSingleVariableConstraint.h.

Referenced by computeQpJacobian(), and computeQpResidual().

_xfem

std::shared_ptr<XFEM> XFEMSingleVariableConstraint::_xfem

protected

Pointer to the XFEM controller object.

Definition at line 56 of file XFEMSingleVariableConstraint.h.

Referenced by XFEMSingleVariableConstraint().

---

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

• xfem/include/constraints/XFEMSingleVariableConstraint.h
• xfem/src/constraints/XFEMSingleVariableConstraint.C