Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "ElementDamper.h"
      11             : 
      12             : // MOOSE includes
      13             : #include "Assembly.h"
      14             : #include "FEProblem.h"
      15             : #include "MooseVariableFE.h"
      16             : #include "SubProblem.h"
      17             : #include "SystemBase.h"
      18             : 
      19             : #include "libmesh/quadrature.h"
      20             : 
      21             : InputParameters
      22        6301 : ElementDamper::validParams()
      23             : {
      24        6301 :   InputParameters params = Damper::validParams();
      25        6301 :   params += MaterialPropertyInterface::validParams();
      26       18903 :   params.addRequiredParam<NonlinearVariableName>(
      27             :       "variable", "The name of the variable that this damper operates on");
      28        6301 :   return params;
      29           0 : }
      30             : 
      31         104 : ElementDamper::ElementDamper(const InputParameters & parameters)
      32             :   : Damper(parameters),
      33             :     MaterialPropertyInterface(this, Moose::EMPTY_BLOCK_IDS, Moose::EMPTY_BOUNDARY_IDS),
      34         104 :     _tid(parameters.get<THREAD_ID>("_tid")),
      35         104 :     _assembly(_subproblem.assembly(_tid, _sys.number())),
      36         104 :     _coord_sys(_assembly.coordSystem()),
      37         104 :     _var(_sys.getFieldVariable<Real>(_tid, parameters.get<NonlinearVariableName>("variable"))),
      38             : 
      39         104 :     _current_elem(_var.currentElem()),
      40         104 :     _q_point(_assembly.qPoints()),
      41         104 :     _qrule(_assembly.qRule()),
      42         104 :     _JxW(_assembly.JxW()),
      43             : 
      44         104 :     _u_increment(_var.increment()),
      45             : 
      46         104 :     _u(_var.sln()),
      47         208 :     _grad_u(_var.gradSln())
      48             : {
      49             :   mooseAssert(_var.count() == 1,
      50             :               "ElementDamper only supports scalar variables. Variable '" + _var.name() +
      51             :                   "' has multiple components.");
      52         104 : }
      53             : 
      54             : Real
      55        1554 : ElementDamper::computeDamping()
      56             : {
      57        1554 :   Real damping = 1.0;
      58        1554 :   Real cur_damping = 1.0;
      59             : 
      60       12929 :   for (_qp = 0; _qp < _qrule->n_points(); _qp++)
      61             :   {
      62       11375 :     cur_damping = computeQpDamping();
      63       11375 :     if (cur_damping < damping)
      64        1862 :       damping = cur_damping;
      65             :   }
      66             : 
      67        1554 :   return damping;
      68             : }
      69             : 
      70             : bool
      71        3749 : ElementDamper::variableDefinedOnElement(const Elem * elem) const
      72             : {
      73        3749 :   return _var.hasBlocks(elem->subdomain_id());
      74             : }