https://mooseframework.inl.gov
MFEMTransient.C
Go to the documentation of this file.
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 "MFEMTransient.h"
13 #include "MFEMProblem.h"
15 #include "TimeStepper.h"
16 
18 
21 {
23  params += TransientBase::validParams();
24  params.addClassDescription("Executioner for transient MFEM problems.");
25  return params;
26 }
27 
29  : TransientBase(params),
30  _mfem_problem(dynamic_cast<MFEMProblem &>(feProblem())),
31  _mfem_problem_data(_mfem_problem.getProblemData()),
32  _mfem_problem_solve(*this, getProblemOperators())
33 {
34  // If no ProblemOperators have been added by the user, add a default
35  if (getProblemOperators().empty())
36  {
37  _mfem_problem_data.eqn_system = std::make_shared<Moose::MFEM::TimeDependentEquationSystem>();
38  auto problem_operator =
39  std::make_shared<Moose::MFEM::TimeDomainEquationSystemProblemOperator>(_mfem_problem);
40  addProblemOperator(std::move(problem_operator));
41  }
42 }
43 
44 void
46 {
48 
49  // Set up initial conditions
53  getParam<MooseEnum>("assembly_level").getEnum<mfem::AssemblyLevel>());
54 
55  for (const auto & problem_operator : getProblemOperators())
56  {
57  problem_operator->SetGridFunctions();
58  problem_operator->Init(_mfem_problem_data.f);
59  }
60 }
61 
62 void
64 {
65  _dt_old = _dt;
66 
67  if (input_dt == -1.0)
69  else
70  _dt = input_dt;
71 
75 
76  // Advance time step of the MFEM problem. Time is also updated here, and
77  // _problem_operator->SetTime is called inside the ode_solver->Step method to
78  // update the time used by time dependent (function) coefficients.
80 
81  // Continue with usual TransientBase::takeStep() finalisation
83 
84  if (!lastSolveConverged())
85  {
86  _console << "Aborting as solve did not converge" << std::endl;
87  return;
88  }
89 
90  if (lastSolveConverged())
92  else
94 
95  // Set time to time old, since final time is updated in TransientBase::endStep()
96  _time = _time_old;
97 
99 }
100 
101 #endif
void timestepSetup() override
static InputParameters validParams()
Definition: TransientBase.C:69
bool & _last_solve_converged
Whether or not the last solve converged.
virtual std::vector< std::shared_ptr< ProblemOperatorBase > > & getProblemOperators()
Returns a pointer to the operator&#39;s equation system.
The main MOOSE class responsible for handling user-defined parameters in almost every MOOSE system...
virtual void preSolve()
Definition: TimeStepper.h:36
Real & _time
Current time.
virtual void takeStep(Real input_dt=-1.0) override
Do whatever is necessary to advance one step.
Definition: MFEMTransient.C:63
virtual bool converged() const
If the time step converged.
Definition: TimeStepper.C:197
Moose::MFEM::FESpaces fespaces
TimeStepper * _time_stepper
virtual void init() override
Initialize the executioner.
Definition: MFEMTransient.C:45
virtual Real computeConstrainedDT()
FEProblemBase & _problem
Here for backward compatibility.
mfem::BlockVector f
virtual void acceptStep()
This gets called when time step is accepted.
Definition: TimeStepper.C:175
Real & _time_old
Previous time.
MFEMTransient(const InputParameters &params)
Definition: MFEMTransient.C:28
Real & _dt
Current delta t... or timestep size.
Base class for transient executioners that use a FixedPointSolve solve object for multiapp-main app i...
Definition: TransientBase.h:27
static InputParameters validParams()
std::shared_ptr< Moose::MFEM::EquationSystem > eqn_system
virtual void init() override
Initialize the executioner.
virtual void onTimestepBegin() override
virtual void rejectStep()
This gets called when time step is rejected.
Definition: TimeStepper.C:185
virtual void postSolve()
Definition: TimeStepper.h:37
virtual void step()
Take a time step.
Definition: TimeStepper.C:166
virtual void addProblemOperator(std::shared_ptr< ProblemOperatorBase > problem_operator)
Add an MFEM problem operator. Takes ownership.
MFEMProblem & _mfem_problem
Definition: MFEMTransient.h:48
void addClassDescription(const std::string &doc_string)
This method adds a description of the class that will be displayed in the input file syntax dump...
const ConsoleStream _console
An instance of helper class to write streams to the Console objects.
virtual bool lastSolveConverged() const override
Whether or not the last solve converged.
registerMooseObject("MooseApp", MFEMTransient)
static InputParameters validParams()
Definition: MFEMTransient.C:20
Moose::MFEM::GridFunctions gridfunctions
MFEMProblemData & _mfem_problem_data
Definition: MFEMTransient.h:49