forked from idaholab/moose
/
LevelSetMeshRefinementTransfer.C
82 lines (72 loc) · 2.8 KB
/
LevelSetMeshRefinementTransfer.C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "LevelSetMeshRefinementTransfer.h"
// MOOSE includes
#include "Adaptivity.h"
#include "FEProblem.h"
#include "MooseVariable.h"
#include "MultiApp.h"
#include "LevelSetTypes.h"
registerMooseObject("LevelSetApp", LevelSetMeshRefinementTransfer);
InputParameters
LevelSetMeshRefinementTransfer::validParams()
{
InputParameters params = MultiAppCopyTransfer::validParams();
params.addClassDescription("Transfers the mesh from the master application to the sub "
"application for the purposes of level set reinitialization problems "
"with mesh adaptivity.");
params.set<MultiMooseEnum>("direction") = "TO_MULTIAPP";
params.suppressParameter<MultiMooseEnum>("direction");
ExecFlagEnum & exec = params.set<ExecFlagEnum>("execute_on");
exec.addAvailableFlags(LevelSet::EXEC_ADAPT_MESH, LevelSet::EXEC_COMPUTE_MARKERS);
exec = {LevelSet::EXEC_COMPUTE_MARKERS, LevelSet::EXEC_ADAPT_MESH};
params.set<bool>("check_multiapp_execute_on") = false;
params.suppressParameter<ExecFlagEnum>("execute_on");
return params;
}
LevelSetMeshRefinementTransfer::LevelSetMeshRefinementTransfer(const InputParameters & parameters)
: MultiAppCopyTransfer(parameters)
{
}
void
LevelSetMeshRefinementTransfer::initialSetup()
{
FEProblemBase & from_problem = _to_multi_app->problemBase();
for (unsigned int i = 0; i < _to_multi_app->numGlobalApps(); i++)
if (_to_multi_app->hasLocalApp(i))
{
FEProblemBase & to_problem = _to_multi_app->appProblemBase(i);
MooseVariable & to_var = to_problem.getStandardVariable(0, _to_var_name);
Adaptivity & adapt = to_problem.adaptivity();
adapt.setMarkerVariableName(to_var.name());
adapt.setCyclesPerStep(from_problem.adaptivity().getCyclesPerStep());
adapt.init(1, 0);
adapt.setUseNewSystem();
adapt.setMaxHLevel(from_problem.adaptivity().getMaxHLevel());
adapt.setAdaptivityOn(false);
}
}
void
LevelSetMeshRefinementTransfer::execute()
{
if (_current_execute_flag == LevelSet::EXEC_COMPUTE_MARKERS)
MultiAppCopyTransfer::execute();
else if (_current_execute_flag == LevelSet::EXEC_ADAPT_MESH)
{
for (unsigned int i = 0; i < _to_multi_app->numGlobalApps(); i++)
if (_to_multi_app->hasLocalApp(i))
{
FEProblemBase & to_problem = _to_multi_app->appProblemBase(i);
Adaptivity & adapt = to_problem.adaptivity();
adapt.setAdaptivityOn(true);
to_problem.adaptMesh();
adapt.setAdaptivityOn(false);
}
}
}