Pressure Class Reference

#include <Pressure.h>

Inheritance diagram for Pressure:

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeQpJacobian () override
virtual void	precalculateQpJacobian () override
virtual Real	computeQpOffDiagJacobian (const unsigned int jvar_num) override
virtual void	precalculateQpOffDiagJacobian (const MooseVariableFEBase &jvar) override
Real	computeStiffness (const unsigned int coupled_component)
Real	computeFaceStiffness (const unsigned int local_j, const unsigned int coupled_component)
virtual void	initialSetup () override
virtual GenericReal< is_ad >	computeQpResidual () override
GenericReal< is_ad >	computeFactor () const

Protected Attributes
unsigned int	_component
	displacement component to apply the bc to More...
const unsigned int	_ndisp
	Number of displacement variables. More...
const Real	_alpha
	_alpha Parameter for HHT time integration scheme More...
std::vector< unsigned int >	_disp_var
	Variable numbers of coupled displacement variables. More...
const FEBase *const &	_fe_side
const std::vector< RealGradient > *	_q_dxi
const std::vector< RealGradient > *	_q_deta
const std::vector< std::vector< Real > > *	_phi_dxi
const std::vector< std::vector< Real > > *	_phi_deta
const bool	_use_displaced_mesh
std::vector< std::unique_ptr< FEBase > >	_fe
std::map< unsigned int, unsigned int >	_node_map
Moose::CoordinateSystemType	_coord_type
	Coordinate system type. More...
	usingTransientInterfaceMembers
const Real	_factor
	Pressure value constant factor, function factor, and postprocessor factor. More...
const Function *const	_function
const PostprocessorValue *const	_postprocessor

Detailed Description

Definition at line 89 of file Pressure.h.

Member Function Documentation

computeFaceStiffness()

Real Pressure::computeFaceStiffness ( const unsigned int  local_j, const unsigned int  coupled_component  )

protected

Definition at line 122 of file Pressure.C.

Referenced by computeStiffness().

{
  //
  // Note that this approach will break down for shell elements, i.e.,
  //   topologically 2D elements in 3D space with pressure loads on
  //   the faces.
  //
  const Real phi_dxi = (*_phi_dxi)[local_j][_qp];
  const Real phi_deta = _phi_deta ? (*_phi_deta)[local_j][_qp] : 0;

  const RealGradient & dqdxi = (*_q_dxi)[_qp];
  const RealGradient out_of_plane(0, 0, 1);
  const RealGradient & dqdeta = _q_deta ? (*_q_deta)[_qp] : out_of_plane;
  // Here, b is dqdxi (cross) dqdeta
  // Then, normal is b/length(b)
  RealGradient b(dqdxi(1) * dqdeta(2) - dqdxi(2) * dqdeta(1),
                 dqdxi(2) * dqdeta(0) - dqdxi(0) * dqdeta(2),
                 dqdxi(0) * dqdeta(1) - dqdxi(1) * dqdeta(0));
  const Real inv_length = 1 / (b * _normals[_qp]);

  const unsigned int i = _component;
  const unsigned int j = coupled_component;

  // const int posneg[3][3] = {{0, -1, 1}, {1, 0, -1}, {-1, 1, 0}};
  const int posneg = 1 - (j + (2 - (i + 1) % 3)) % 3;

  // const unsigned int index[3][3] = {{0, 2, 1}, {2, 1, 0}, {1, 0, 2}};
  const unsigned int index = 2 - (j + (i + 2) % 3) % 3;

  const Real variation_b = posneg * (phi_deta * dqdxi(index) - phi_dxi * dqdeta(index));

  Real rz_term = 0;
  if (_coord_type == Moose::COORD_RZ && j == _subproblem.getAxisymmetricRadialCoord())
  {
    rz_term = _normals[_qp](i) * _phi[_j][_qp] / _q_point[_qp](0);
  }

  return computeFactor() * _test[_i][_qp] * (inv_length * variation_b + rz_term);
}

computeFactor()

GenericReal< is_ad > PressureTempl< is_ad >::computeFactor ( ) const

protectedinherited

Definition at line 108 of file Pressure.C.

Referenced by computeFaceStiffness(), and computeStiffness().

{
  GenericReal<is_ad> factor = _factor;

  if (_function)
    factor *= _function->value(_t + _alpha * _dt, _q_point[_qp]);

  if (_postprocessor)
    factor *= *_postprocessor;

  return factor;
}

computeQpJacobian()

Real Pressure::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 190 of file Pressure.C.

{
  if (_use_displaced_mesh)
    return computeStiffness(_component);

  return 0;
}

computeQpOffDiagJacobian()

Real Pressure::computeQpOffDiagJacobian ( const unsigned int  jvar_num )

overrideprotectedvirtual

Definition at line 199 of file Pressure.C.

{
  if (_use_displaced_mesh)
    for (unsigned int j = 0; j < _ndisp; ++j)
      if (jvar_num == _disp_var[j])
        return computeStiffness(j);

  return 0;
}

computeQpResidual()

GenericReal< is_ad > PressureTempl< is_ad >::computeQpResidual ( )

overrideprotectedvirtualinherited

Definition at line 101 of file Pressure.C.

{
  return computeFactor() * (_normals[_qp](_component) * _test[_i][_qp]);
}

computeStiffness()

Real Pressure::computeStiffness ( const unsigned int  coupled_component )

protected

Definition at line 163 of file Pressure.C.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

{
  if (_ndisp > 1)
  {
    const std::map<unsigned int, unsigned int>::iterator j_it = _node_map.find(_j);
    if (_test[_i][_qp] == 0 || j_it == _node_map.end())
      return 0;

    return computeFaceStiffness(j_it->second, coupled_component);
  }

  else if (_coord_type == Moose::COORD_RSPHERICAL)
  {
    return computeFactor() * _normals[_qp](_component) * _test[_i][_qp] * _phi[_j][_qp] *
           (2 / _q_point[_qp](0));
  }

  if (_coord_type == Moose::COORD_RZ)
  {
    return computeFactor() * _normals[_qp](_component) * _test[_i][_qp] * _phi[_j][_qp] /
           _q_point[_qp](0);
  }

  return 0;
}

initialSetup()

void PressureTempl< is_ad >::initialSetup ( )

overrideprotectedvirtualinherited

Definition at line 81 of file Pressure.C.

{
  auto boundary_ids = this->boundaryIDs();
  std::set<SubdomainID> block_ids;
  for (auto bndry_id : boundary_ids)
  {
    auto bids = _mesh.getBoundaryConnectedBlocks(bndry_id);
    block_ids.insert(bids.begin(), bids.end());
  }

  _coord_type = _fe_problem.getCoordSystem(*block_ids.begin());
  for (auto blk_id : block_ids)
  {
    if (_coord_type != _fe_problem.getCoordSystem(blk_id))
      mooseError("The Pressure BC requires subdomains to have the same coordinate system.");
  }
}

precalculateQpJacobian()

void Pressure::precalculateQpJacobian ( )

overrideprotectedvirtual

Definition at line 210 of file Pressure.C.

Referenced by precalculateQpOffDiagJacobian().

{
  if (_ndisp == 1)
    return;

  if (_fe[_tid] == nullptr)
  {
    const unsigned int dim = _sys.mesh().dimension() - 1;
    QBase * const & qrule = _assembly.writeableQRuleFace();
    _fe[_tid] = FEBase::build(dim, _var.feType());
    _fe[_tid]->attach_quadrature_rule(qrule);
  }

  _q_dxi = &_fe[_tid]->get_dxyzdxi();
  _phi_dxi = &_fe[_tid]->get_dphidxi();
  if (_coord_type == Moose::COORD_XYZ)
  {
    _q_deta = &_fe[_tid]->get_dxyzdeta();
    _phi_deta = &_fe[_tid]->get_dphideta();
  }

  _fe[_tid]->reinit(_current_side_elem);

  if (_coord_type == Moose::COORD_XYZ)
  {
    if (_q_deta->empty())
      _q_deta = nullptr;
    if (_phi_deta->empty())
      _phi_deta = nullptr;
  }

  // Compute node map (given elem node, supply face node)
  _node_map.clear();
  const unsigned int num_node_elem = _current_elem->n_nodes();
  const Node * const * elem_nodes = _current_elem->get_nodes();
  const unsigned int num_node_face = _current_side_elem->n_nodes();
  const Node * const * face_nodes = _current_side_elem->get_nodes();
  unsigned int num_found = 0;
  for (unsigned i = 0; i < num_node_elem; ++i)
  {
    for (unsigned j = 0; j < num_node_face; ++j)
      if (**(elem_nodes + i) == **(face_nodes + j))
      {
        _node_map[i] = j;
        ++num_found;
        break;
      }
    if (num_found == num_node_face)
      break;
  }
}

precalculateQpOffDiagJacobian()

void Pressure::precalculateQpOffDiagJacobian ( const MooseVariableFEBase &  jvar )

overrideprotectedvirtual

Definition at line 263 of file Pressure.C.

{
  precalculateQpJacobian();
}

validParams()

InputParameters PressureTempl< is_ad >::validParams ( )

staticinherited

Definition at line 21 of file Pressure.C.

{
  InputParameters params = PressureParent<is_ad>::validParams();
  params.addClassDescription("Applies a pressure on a given boundary in a given direction");
  params.addDeprecatedParam<unsigned int>(
      "component", "The component for the pressure", "This parameter is no longer necessary");
  params.addRequiredCoupledVar("displacements",
                               "The string of displacements suitable for the problem statement");
  params.addDeprecatedParam<Real>("constant",
                                  "The magnitude to use in computing the pressure",
                                  "Use 'factor' in place of 'constant'");
  params.addParam<Real>("factor", 1.0, "The magnitude to use in computing the pressure");
  params.addParam<FunctionName>("function", "The function that describes the pressure");
  params.addParam<PostprocessorName>("postprocessor",
                                     "Postprocessor that will supply the pressure value");
  params.addParam<Real>("alpha", 0.0, "alpha parameter required for HHT time integration scheme");
  params.set<bool>("use_displaced_mesh") = true;
  return params;
}

Member Data Documentation

_alpha

const Real PressureTempl< is_ad >::_alpha

protectedinherited

_alpha Parameter for HHT time integration scheme

Definition at line 55 of file Pressure.h.

_component

unsigned int PressureTempl< is_ad >::_component

protectedinherited

displacement component to apply the bc to

Definition at line 43 of file Pressure.h.

Referenced by computeFaceStiffness(), computeQpJacobian(), and computeStiffness().

_coord_type

Moose::CoordinateSystemType PressureTempl< is_ad >::_coord_type

protectedinherited

Coordinate system type.

Definition at line 73 of file Pressure.h.

Referenced by computeFaceStiffness(), computeStiffness(), and precalculateQpJacobian().

_disp_var

std::vector<unsigned int> PressureTempl< is_ad >::_disp_var

protectedinherited

Variable numbers of coupled displacement variables.

Definition at line 58 of file Pressure.h.

Referenced by computeQpOffDiagJacobian().

_factor

const Real PressureTempl< is_ad >::_factor

protectedinherited

Pressure value constant factor, function factor, and postprocessor factor.

Definition at line 49 of file Pressure.h.

_fe

std::vector<std::unique_ptr<FEBase> > PressureTempl< is_ad >::_fe

protectedinherited

Definition at line 68 of file Pressure.h.

Referenced by precalculateQpJacobian().

_fe_side

const FEBase* const& PressureTempl< is_ad >::_fe_side

protectedinherited

Definition at line 60 of file Pressure.h.

_function

const Function* const PressureTempl< is_ad >::_function

protectedinherited

Definition at line 50 of file Pressure.h.

_ndisp

const unsigned int PressureTempl< is_ad >::_ndisp

protectedinherited

Number of displacement variables.

Definition at line 46 of file Pressure.h.

Referenced by computeQpOffDiagJacobian(), computeStiffness(), and precalculateQpJacobian().

_node_map

std::map<unsigned int, unsigned int> PressureTempl< is_ad >::_node_map

protectedinherited

Definition at line 70 of file Pressure.h.

Referenced by computeStiffness(), and precalculateQpJacobian().

_phi_deta

const std::vector<std::vector<Real> >* PressureTempl< is_ad >::_phi_deta

protectedinherited

Definition at line 64 of file Pressure.h.

Referenced by computeFaceStiffness(), and precalculateQpJacobian().

_phi_dxi

const std::vector<std::vector<Real> >* PressureTempl< is_ad >::_phi_dxi

protectedinherited

Definition at line 63 of file Pressure.h.

Referenced by precalculateQpJacobian().

_postprocessor

const PostprocessorValue* const PressureTempl< is_ad >::_postprocessor

protectedinherited

Definition at line 51 of file Pressure.h.

_q_deta

const std::vector<RealGradient>* PressureTempl< is_ad >::_q_deta

protectedinherited

Definition at line 62 of file Pressure.h.

Referenced by computeFaceStiffness(), and precalculateQpJacobian().

_q_dxi

const std::vector<RealGradient>* PressureTempl< is_ad >::_q_dxi

protectedinherited

Definition at line 61 of file Pressure.h.

Referenced by precalculateQpJacobian().

_use_displaced_mesh

const bool PressureTempl< is_ad >::_use_displaced_mesh

protectedinherited

Definition at line 65 of file Pressure.h.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

usingTransientInterfaceMembers

PressureTempl< is_ad >::usingTransientInterfaceMembers

protectedinherited

Definition at line 75 of file Pressure.h.

---

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

• solid_mechanics/include/bcs/Pressure.h
• solid_mechanics/src/bcs/Pressure.C