PressureBase Class Reference

#include <PressureBase.h>

Inheritance diagram for PressureBase:

Public Member Functions
	PressureBase (const InputParameters &parameters)
virtual void	initialSetup () override

Static Public Member Functions
static InputParameters	validParams ()
static InputParameters	actionParams ()

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

Protected Attributes
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
std::vector< unsigned int >	_disp_var
	Variable numbers of coupled displacement variables. More...
const unsigned int	_ndisp
	Number of displacement variables. More...
const unsigned int	_component
	displacement component to apply the bc to More...
Moose::CoordinateSystemType	_coord_type
	Coordinate system type. More...
	usingTransientInterfaceMembers

Detailed Description

Definition at line 70 of file PressureBase.h.

Constructor & Destructor Documentation

PressureBase()

PressureBase::PressureBase

(

const InputParameters &

parameters

)

Definition at line 104 of file PressureBase.C.

  : PressureBaseTempl<false>(parameters),
    _q_dxi(nullptr),
    _q_deta(nullptr),
    _phi_dxi(nullptr),
    _phi_deta(nullptr),
    _use_displaced_mesh(this->template getParam<bool>("use_displaced_mesh")),
    _fe(libMesh::n_threads())
{
}

Member Function Documentation

actionParams()

InputParameters PressureBaseTempl< is_ad >::actionParams ( )

staticinherited

Definition at line 30 of file PressureBase.C.

{
  auto params = emptyInputParameters();
  params.addRequiredCoupledVar("displacements",
                               "The string of displacements suitable for the problem statement");
  return params;
}

computeFaceStiffness()

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

protected

Definition at line 116 of file PressureBase.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 computePressure() * _test[_i][_qp] * (inv_length * variation_b + rz_term);
}

computePressure()

virtual GenericReal<is_ad> PressureBaseTempl< is_ad >::computePressure ( ) const

protectedpure virtualinherited

Referenced by computeFaceStiffness(), and computeStiffness().

computeQpJacobian()

Real PressureBase::computeQpJacobian ( )

finaloverrideprotectedvirtual

Definition at line 184 of file PressureBase.C.

{
  if (_use_displaced_mesh)
    return computeStiffness(_component);

  return 0;
}

computeQpOffDiagJacobian()

Real PressureBase::computeQpOffDiagJacobian ( const unsigned int  jvar_num )

finaloverrideprotectedvirtual

Definition at line 193 of file PressureBase.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 > PressureBaseTempl< is_ad >::computeQpResidual ( )

finaloverrideprotectedvirtualinherited

Definition at line 99 of file PressureBase.C.

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

computeStiffness()

Real PressureBase::computeStiffness ( const unsigned int  coupled_component )

protected

Definition at line 157 of file PressureBase.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 computePressure() * _normals[_qp](_component) * _test[_i][_qp] * _phi[_j][_qp] *
           (2 / _q_point[_qp](0));
  }

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

  return 0;
}

initialSetup()

void PressureBaseTempl< is_ad >::initialSetup ( )

overridevirtualinherited

Definition at line 72 of file PressureBase.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());
  }

  if (block_ids.size())
    _coord_type = _fe_problem.getCoordSystem(*block_ids.begin());
  else
  {
    mooseInfo(
        "No connected blocks were found, the coordinate system type is obtained from the mesh.");
    _coord_type = _fe_problem.mesh().getUniqueCoordSystem();
  }
  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 PressureBase::precalculateQpJacobian ( )

finaloverrideprotectedvirtual

Definition at line 204 of file PressureBase.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 PressureBase::precalculateQpOffDiagJacobian ( const MooseVariableFEBase &  jvar )

finaloverrideprotectedvirtual

Definition at line 257 of file PressureBase.C.

{
  precalculateQpJacobian();
}

validParams()

InputParameters PressureBaseTempl< is_ad >::validParams ( )

staticinherited

Definition at line 19 of file PressureBase.C.

{
  InputParameters params = PressureBaseParent<is_ad>::validParams();
  params.addDeprecatedParam<unsigned int>(
      "component", "The component for the pressure", "This parameter is no longer necessary");
  params.addParam<bool>("use_displaced_mesh", true, "Whether to use the displaced mesh.");
  return params;
}

Member Data Documentation

_component

const unsigned int PressureBaseTempl< is_ad >::_component

protectedinherited

displacement component to apply the bc to

Definition at line 51 of file PressureBase.h.

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

_coord_type

Moose::CoordinateSystemType PressureBaseTempl< is_ad >::_coord_type

protectedinherited

Coordinate system type.

Definition at line 54 of file PressureBase.h.

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

_disp_var

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

protectedinherited

Variable numbers of coupled displacement variables.

Definition at line 45 of file PressureBase.h.

Referenced by computeQpOffDiagJacobian().

_fe

std::vector<std::unique_ptr<FEBase> > PressureBase::_fe

protected

Definition at line 91 of file PressureBase.h.

Referenced by precalculateQpJacobian().

_ndisp

const unsigned int PressureBaseTempl< is_ad >::_ndisp

protectedinherited

Number of displacement variables.

Definition at line 48 of file PressureBase.h.

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

_node_map

std::map<unsigned int, unsigned int> PressureBase::_node_map

protected

Definition at line 93 of file PressureBase.h.

Referenced by computeStiffness(), and precalculateQpJacobian().

_phi_deta

const std::vector<std::vector<Real> >* PressureBase::_phi_deta

protected

Definition at line 87 of file PressureBase.h.

Referenced by computeFaceStiffness(), and precalculateQpJacobian().

_phi_dxi

const std::vector<std::vector<Real> >* PressureBase::_phi_dxi

protected

Definition at line 86 of file PressureBase.h.

Referenced by precalculateQpJacobian().

_q_deta

const std::vector<RealGradient>* PressureBase::_q_deta

protected

Definition at line 85 of file PressureBase.h.

Referenced by computeFaceStiffness(), and precalculateQpJacobian().

_q_dxi

const std::vector<RealGradient>* PressureBase::_q_dxi

protected

Definition at line 84 of file PressureBase.h.

Referenced by precalculateQpJacobian().

_use_displaced_mesh

const bool PressureBase::_use_displaced_mesh

protected

Definition at line 88 of file PressureBase.h.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

usingTransientInterfaceMembers

PressureBaseTempl< is_ad >::usingTransientInterfaceMembers

protectedinherited

Definition at line 56 of file PressureBase.h.

---

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

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