CriticalTimeStep.C

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

#include "CriticalTimeStep.h"

registerMooseObject("SolidMechanicsApp", CriticalTimeStep);

InputParameters
CriticalTimeStep::validParams()
{
  InputParameters params = ElementPostprocessor::validParams();
  params.addClassDescription(
      "Computes and reports the critical time step for the explicit solver.");
  params.addParam<MaterialPropertyName>(
      "density", "density", "Name of Material Property defining the density of the material.");
  params.addParam<MaterialPropertyName>(
      "density_scaling", "Name of material property to add mass scaling in explicit simulations.");
  params.addParam<Real>("factor", 1.0, "Factor to mulitply to the critical time step.");
  return params;
}

CriticalTimeStep::CriticalTimeStep(const InputParameters & parameters)
  : ElementPostprocessor(parameters),
    GuaranteeConsumer(this),
    _material_density(getMaterialProperty<Real>("density")),
    _density_scaling(parameters.isParamSetByUser("density_scaling")
                         ? &getMaterialPropertyByName<Real>("density_scaling")
                         : nullptr),
    _sqrt_effective_stiffness(getMaterialPropertyByName<Real>("effective_stiffness")),
    _factor(getParam<Real>("factor")),
    _critical_time(parameters.isParamValid("critical_time"))
{
}

void
CriticalTimeStep::initialSetup()
{
  if (!hasGuaranteedMaterialProperty("effective_stiffness", Guarantee::ISOTROPIC))
    paramError("CriticalTimeStep can only be used with elasticity tensor materials "
               "that guarantee isotropic tensors.");
}

void
CriticalTimeStep::initialize()
{
  _critical_time = std::numeric_limits<Real>::max();
}

void
CriticalTimeStep::execute()
{
  const Real dens = _material_density[0];
  const Real dens_scaling = _density_scaling ? (*_density_scaling)[0] : 0.0;

  // In the above line, density is inferred only at the first quadrature point
  // of each element. Since critical time step is computed across all elements and
  // a minimum is then taken, this is okay.
  _critical_time = std::min(_factor * _current_elem->hmin() * std::sqrt(dens + dens_scaling) /
                                (_sqrt_effective_stiffness[0]),
                            _critical_time);
}

void
CriticalTimeStep::finalize()
{
  gatherMin(_critical_time);
}

Real
CriticalTimeStep::getValue() const
{
  return _critical_time;
}

void
CriticalTimeStep::threadJoin(const UserObject & y)
{
  const auto & pps = static_cast<const CriticalTimeStep &>(y);
  _critical_time = std::min(pps._critical_time, _critical_time);
}