ModularGapConductanceConstraint.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "ADMortarConstraint.h"

class GapFluxModelBase;

class ModularGapConductanceConstraint : public ADMortarConstraint
{
public:
  static InputParameters validParams();

  ModularGapConductanceConstraint(const InputParameters & parameters);

  virtual void initialSetup() override;

protected:
  virtual ADReal computeQpResidual(Moose::MortarType mortar_type) override;

  std::vector<UserObjectName> _gap_flux_model_names;

  std::vector<const GapFluxModelBase *> _gap_flux_models;

  const std::vector<std::string> _disp_name;
  const unsigned int _n_disp;
  std::vector<const ADVariableValue *> _disp_secondary;
  std::vector<const ADVariableValue *> _disp_primary;

private:
  virtual ADReal computeSurfaceIntegrationFactor() const;

  virtual ADReal computeGapLength() const;

  void computeGapRadii(const ADReal & gap_length);

  enum class GapGeometry
  {
    AUTO,
    PLATE,
    CYLINDER,
    SPHERE,
  } _gap_geometry_type;

  virtual void setGapGeometryParameters(const InputParameters & params,
                                        const Moose::CoordinateSystemType coord_sys,
                                        unsigned int axisymmetric_radial_coord,
                                        GapGeometry & gap_geometry_type);

  void deduceGeometryParameters();

  ADReal _gap_width;

  ADReal _surface_integration_factor;

  Point & _p1;
  Point & _p2;

  ADReal _r1;
  ADReal _r2;
  ADReal _radius;

  const Real _max_gap;

  ADReal _adjusted_length;

  const MooseVariable * const _disp_x_var;
  const MooseVariable * const _disp_y_var;
  const MooseVariable * const _disp_z_var;

  ADReal _normal_pressure;

  friend class GapFluxModelBase;
};