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             : #ifdef MOOSE_MFEM_ENABLED
      11             : 
      12             : #include "MFEMScalarTimeAverageAux.h"
      13             : #include "MFEMProblem.h"
      14             : 
      15             : registerMooseObject("MooseApp", MFEMScalarTimeAverageAux);
      16             : 
      17             : InputParameters
      18        2110 : MFEMScalarTimeAverageAux::validParams()
      19             : {
      20        2110 :   InputParameters params = MFEMAuxKernel::validParams();
      21        4220 :   params.addClassDescription(
      22             :       "Calculates a running time average of a scalar coefficient projected onto an auxvariable");
      23        8440 :   params.addRequiredParam<MFEMScalarCoefficientName>("source", "Scalar coefficient to average");
      24        6330 :   params.addParam<mfem::real_t>("time_skip", 0.0, "Time to skip before beginning the average");
      25             : 
      26        2110 :   return params;
      27           0 : }
      28             : 
      29           6 : MFEMScalarTimeAverageAux::MFEMScalarTimeAverageAux(const InputParameters & parameters)
      30             :   : MFEMAuxKernel(parameters),
      31           6 :     _source_coefficient(getScalarCoefficient("source")),
      32           6 :     _result_coefficient(getScalarCoefficientByName(_result_var_name)),
      33           6 :     _average_var(_result_var.ParFESpace()),
      34          12 :     _skip(getParam<Real>("time_skip")),
      35           6 :     _time(getMFEMProblem().time()),
      36          12 :     _dt(getMFEMProblem().dt())
      37             : {
      38           6 : }
      39             : 
      40             : void
      41          18 : MFEMScalarTimeAverageAux::execute()
      42             : {
      43          18 :   if (_time <= _skip)
      44           6 :     return;
      45             : 
      46             :   // Linear blend: (1 - w) * avg_old + w * src
      47          12 :   const mfem::real_t w = _dt / (_time - _skip);
      48          12 :   mfem::SumCoefficient blend(_result_coefficient, _source_coefficient, 1.0 - w, w);
      49          12 :   _average_var.ProjectCoefficient(blend);
      50             : 
      51          12 :   _result_var = _average_var;
      52          12 : }
      53             : 
      54             : #endif