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. More...
Real	_jump
	Vector normal to the internal interface. More...
Real	_jump_flux
	Vector normal to the internal 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 25 of file XFEMSingleVariableConstraint.h.

Constructor & Destructor Documentation

XFEMSingleVariableConstraint()

XFEMSingleVariableConstraint::XFEMSingleVariableConstraint

(

const InputParameters &

parameters

)

Definition at line 41 of file XFEMSingleVariableConstraint.C.

  : ElemElemConstraint(parameters),
    _alpha(getParam<Real>("alpha")),
    _jump(getParam<Real>("jump")),
    _jump_flux(getParam<Real>("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 61 of file XFEMSingleVariableConstraint.C.

{}

Member Function Documentation

computeQpJacobian()

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

overrideprotectedvirtual

Definition at line 107 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 71 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;
      }
      r += 0.5 * _test[_i][_qp] * _jump_flux;
      r += _alpha * (_u[_qp] - _u_neighbor[_qp] - _jump) * _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;
      }
      r += 0.5 * _test_neighbor[_i][_qp] * _jump_flux;
      r -= _alpha * (_u[_qp] - _u_neighbor[_qp] - _jump) * _test_neighbor[_i][_qp];
      break;
  }
  return r;
}

reinitConstraintQuadrature()

void XFEMSingleVariableConstraint::reinitConstraintQuadrature ( const ElementPairInfo &  element_pair_info )

overrideprotectedvirtual

Definition at line 64 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.

Definition at line 42 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 39 of file XFEMSingleVariableConstraint.h.

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

_jump

Real XFEMSingleVariableConstraint::_jump

protected

Vector normal to the internal interface.

Definition at line 45 of file XFEMSingleVariableConstraint.h.

Referenced by computeQpResidual().

_jump_flux

Real XFEMSingleVariableConstraint::_jump_flux

protected

Vector normal to the internal interface.

Definition at line 48 of file XFEMSingleVariableConstraint.h.

Referenced by computeQpResidual().

_use_penalty

bool XFEMSingleVariableConstraint::_use_penalty

protected

Use penalty formulation.

Definition at line 51 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 54 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