XFEMPressure Class Reference

#include <XFEMPressure.h>

Inheritance diagram for XFEMPressure:

Public Member Functions
	XFEMPressure (const InputParameters &parameters)
virtual void	addPoints ()
virtual Real	computeQpResidual ()

Protected Attributes
const int	_component
const Real	_factor
const Function *const	_function
std::map< unsigned int, std::shared_ptr< ElementPairLocator > > *	_element_pair_locators
std::map< const Elem *, std::map< unsigned int, Point > >	_elem_qp_normal
std::map< const Elem *, std::map< unsigned int, Real > >	_elem_qp_JxW

Detailed Description

Definition at line 17 of file XFEMPressure.h.

Constructor & Destructor Documentation

XFEMPressure()

XFEMPressure::XFEMPressure

(

const InputParameters &

parameters

)

Definition at line 29 of file XFEMPressure.C.

  : DiracKernel(parameters),
    _component(getParam<unsigned int>("component")),
    _factor(getParam<Real>("factor")),
    _function(isParamValid("function") ? &getFunction("function") : NULL)
{
  GeometricSearchData & geo_search_data = _subproblem.geomSearchData();
  _element_pair_locators = &geo_search_data._element_pair_locators;
}

Member Function Documentation

addPoints()

void XFEMPressure::addPoints ( )

virtual

Definition at line 40 of file XFEMPressure.C.

{
  _elem_qp_normal.clear();
  _elem_qp_JxW.clear();
 
  for (std::map<unsigned int, std::shared_ptr<ElementPairLocator>>::iterator it_epl =
           _element_pair_locators->begin();
       it_epl != _element_pair_locators->end();
       ++it_epl)
  {
    ElementPairLocator & elem_pair_loc = *it_epl->second;
    // go over pair elements
    const std::list<std::pair<const Elem *, const Elem *>> & elem_pairs =
        elem_pair_loc.getElemPairs();
    for (std::list<std::pair<const Elem *, const Elem *>>::const_iterator it_ep =
             elem_pairs.begin();
         it_ep != elem_pairs.end();
         ++it_ep)
    {
      const Elem * elem1 = it_ep->first;
      const Elem * elem2 = it_ep->second;
      const ElementPairInfo & info = elem_pair_loc.getElemPairInfo(*it_ep);
 
      for (unsigned int i = 0; i < info._elem1_constraint_q_point.size(); ++i)
      {
        _elem_qp_normal[elem1][i] = -info._elem1_normal;
        _elem_qp_normal[elem2][i] = -info._elem2_normal;
        _elem_qp_JxW[elem1][i] = info._elem1_constraint_JxW[i];
        _elem_qp_JxW[elem2][i] = info._elem2_constraint_JxW[i];
        addPoint(elem1, info._elem1_constraint_q_point[i]);
        addPoint(elem2, info._elem2_constraint_q_point[i]);
      }
    }
  }
}

computeQpResidual()

Real XFEMPressure::computeQpResidual ( )

virtual

Definition at line 77 of file XFEMPressure.C.

{
  Real factor = _factor;
 
  if (_function)
    factor *= _function->value(_t, _current_point);
 
  Point normal = _elem_qp_normal[_current_elem][_qp];
  Real JxW = _elem_qp_JxW[_current_elem][_qp];
 
  return -factor * JxW * (normal(_component) * _test[_i][_qp]);
}

Member Data Documentation

_component

const int XFEMPressure::_component

protected

Definition at line 26 of file XFEMPressure.h.

Referenced by computeQpResidual().

_elem_qp_JxW

std::map<const Elem *, std::map<unsigned int, Real> > XFEMPressure::_elem_qp_JxW

protected

Definition at line 32 of file XFEMPressure.h.

Referenced by addPoints(), and computeQpResidual().

_elem_qp_normal

std::map<const Elem *, std::map<unsigned int, Point> > XFEMPressure::_elem_qp_normal

protected

Definition at line 31 of file XFEMPressure.h.

Referenced by addPoints(), and computeQpResidual().

_element_pair_locators

std::map<unsigned int, std::shared_ptr<ElementPairLocator> >* XFEMPressure::_element_pair_locators

protected

Definition at line 30 of file XFEMPressure.h.

Referenced by addPoints(), and XFEMPressure().

_factor

const Real XFEMPressure::_factor

protected

Definition at line 27 of file XFEMPressure.h.

Referenced by computeQpResidual().

_function

const Function* const XFEMPressure::_function

protected

Definition at line 28 of file XFEMPressure.h.

Referenced by computeQpResidual().

---

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

• xfem/include/dirackernels/XFEMPressure.h
• xfem/src/dirackernels/XFEMPressure.C