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 : // MOOSE Includes
11 : #include "MortarConsumerInterface.h"
12 : #include "InputParameters.h"
13 : #include "MooseObject.h"
14 : #include "FEProblemBase.h"
15 : #include "MooseMesh.h"
16 : #include "MortarData.h"
17 : #include "Assembly.h"
18 : #include "libmesh/quadrature.h"
19 :
20 : #include <algorithm>
21 :
22 : InputParameters
23 259468 : MortarConsumerInterface::validParams()
24 : {
25 : // Create InputParameters object that will be appended to the parameters for the inheriting object
26 259468 : InputParameters params = emptyInputParameters();
27 : // On a displaced mesh this will geometrically and algebraically ghost the entire interface
28 259468 : params.addRelationshipManager(
29 : "AugmentSparsityOnInterface",
30 : Moose::RelationshipManagerType::GEOMETRIC | Moose::RelationshipManagerType::ALGEBRAIC,
31 0 : [](const InputParameters & obj_params, InputParameters & rm_params)
32 : {
33 4048 : rm_params.set<bool>("use_displaced_mesh") = obj_params.get<bool>("use_displaced_mesh");
34 8096 : rm_params.set<BoundaryName>("secondary_boundary") =
35 8096 : obj_params.get<BoundaryName>("secondary_boundary");
36 8096 : rm_params.set<BoundaryName>("primary_boundary") =
37 8096 : obj_params.get<BoundaryName>("primary_boundary");
38 8096 : rm_params.set<SubdomainName>("secondary_subdomain") =
39 8096 : obj_params.get<SubdomainName>("secondary_subdomain");
40 8096 : rm_params.set<SubdomainName>("primary_subdomain") =
41 8096 : obj_params.get<SubdomainName>("primary_subdomain");
42 4048 : rm_params.set<bool>("ghost_point_neighbors") =
43 4048 : obj_params.get<bool>("ghost_point_neighbors");
44 4048 : rm_params.set<bool>("ghost_higher_d_neighbors") =
45 4048 : obj_params.get<bool>("ghost_higher_d_neighbors");
46 4048 : });
47 :
48 259468 : params.addRequiredParam<BoundaryName>("primary_boundary",
49 : "The name of the primary boundary sideset.");
50 259468 : params.addRequiredParam<BoundaryName>("secondary_boundary",
51 : "The name of the secondary boundary sideset.");
52 259468 : params.addRequiredParam<SubdomainName>("primary_subdomain", "The name of the primary subdomain.");
53 259468 : params.addRequiredParam<SubdomainName>("secondary_subdomain",
54 : "The name of the secondary subdomain.");
55 778404 : params.addParam<bool>(
56 : "periodic",
57 518936 : false,
58 : "Whether this constraint is going to be used to enforce a periodic condition. This has the "
59 : "effect of changing the normals vector for projection from outward to inward facing");
60 :
61 778404 : params.addParam<bool>(
62 : "debug_mesh",
63 518936 : false,
64 : "Whether this constraint is going to enable mortar segment mesh debug information. An exodus"
65 : "file will be generated if the user sets this flag to true");
66 :
67 778404 : params.addParam<bool>(
68 : "correct_edge_dropping",
69 518936 : false,
70 : "Whether to enable correct edge dropping treatment for mortar constraints. When disabled "
71 : "any Lagrange Multiplier degree of freedom on a secondary element without full primary "
72 : "contributions will be set (strongly) to 0.");
73 :
74 778404 : params.addParam<bool>(
75 : "interpolate_normals",
76 518936 : true,
77 : "Whether to interpolate the nodal normals (e.g. classic idea of evaluating field at "
78 : "quadrature points). If this is set to false, then non-interpolated nodal normals will be "
79 : "used, and then the _normals member should be indexed with _i instead of _qp");
80 :
81 778404 : params.addParam<bool>("ghost_point_neighbors",
82 518936 : false,
83 : "Whether we should ghost point neighbors of secondary face elements, and "
84 : "consequently also their mortar interface couples.");
85 778404 : params.addParam<Real>(
86 : "minimum_projection_angle",
87 518936 : 40.0,
88 : "Parameter to control which angle (in degrees) is admissible for the creation of mortar "
89 : "segments. If set to a value close to zero, very oblique projections are allowed, which "
90 : "can result in mortar segments solving physics not meaningfully, and overprojection of "
91 : "primary nodes onto the mortar segment mesh in extreme cases. This parameter is mostly "
92 : "intended for mortar mesh debugging purposes in two dimensions.");
93 :
94 778404 : params.addParam<bool>(
95 : "ghost_higher_d_neighbors",
96 518936 : false,
97 : "Whether we should ghost higher-dimensional neighbors. This is necessary when we are doing "
98 : "second order mortar with finite volume primal variables, because in order for the method to "
99 : "be second order we must use cell gradients, which couples in the neighbor cells.");
100 :
101 259468 : return params;
102 0 : }
103 :
104 : // Standard constructor
105 1334 : MortarConsumerInterface::MortarConsumerInterface(const MooseObject * moose_object)
106 1334 : : _mci_fe_problem(*moose_object->getCheckedPointerParam<FEProblemBase *>("_fe_problem_base")),
107 1334 : _mci_subproblem(*moose_object->getCheckedPointerParam<SubProblem *>("_subproblem")),
108 1334 : _mci_tid(moose_object->getParam<THREAD_ID>("_tid")),
109 1334 : _mci_mesh(_mci_subproblem.mesh()),
110 : // all geometric assembly information should be correct for nl system number 0
111 1334 : _mci_assembly(_mci_subproblem.assembly(_mci_tid, 0)),
112 1334 : _mortar_data(_mci_fe_problem.mortarData()),
113 1334 : _secondary_id(
114 1334 : _mci_mesh.getBoundaryID(moose_object->getParam<BoundaryName>("secondary_boundary"))),
115 1334 : _primary_id(_mci_mesh.getBoundaryID(moose_object->getParam<BoundaryName>("primary_boundary"))),
116 1334 : _secondary_subdomain_id(
117 1334 : _mci_mesh.getSubdomainID(moose_object->getParam<SubdomainName>("secondary_subdomain"))),
118 1334 : _primary_subdomain_id(
119 1334 : _mci_mesh.getSubdomainID(moose_object->getParam<SubdomainName>("primary_subdomain"))),
120 1334 : _secondary_set({_secondary_id}),
121 1334 : _interpolate_normals(moose_object->getParam<bool>("interpolate_normals")),
122 1334 : _phys_points_secondary(_mci_assembly.qPointsFace()),
123 1334 : _phys_points_primary(_mci_assembly.qPointsFaceNeighbor()),
124 1334 : _qrule_msm(_mci_assembly.qRuleMortar()),
125 1334 : _qrule_face(_mci_assembly.qRuleFace()),
126 1334 : _lower_secondary_elem(_mci_assembly.lowerDElem()),
127 1334 : _lower_primary_elem(_mci_assembly.neighborLowerDElem()),
128 1334 : _JxW_msm(_mci_assembly.jxWMortar()),
129 2668 : _msm_elem(_mci_assembly.msmElem())
130 : {
131 2668 : const bool displaced = moose_object->isParamValid("use_displaced_mesh")
132 2668 : ? moose_object->getParam<bool>("use_displaced_mesh")
133 1334 : : false;
134 :
135 : // Create the mortar interface if it hasn't already been created
136 1334 : _mci_fe_problem.createMortarInterface(
137 1334 : std::make_pair(_primary_id, _secondary_id),
138 0 : std::make_pair(_primary_subdomain_id, _secondary_subdomain_id),
139 : displaced,
140 2668 : moose_object->getParam<bool>("periodic"),
141 2668 : moose_object->getParam<bool>("debug_mesh"),
142 2668 : moose_object->getParam<bool>("correct_edge_dropping"),
143 2668 : moose_object->getParam<Real>("minimum_projection_angle"));
144 :
145 1334 : _amg = &_mci_fe_problem.getMortarInterface(
146 1334 : std::make_pair(_primary_id, _secondary_id),
147 1334 : std::make_pair(_primary_subdomain_id, _secondary_subdomain_id),
148 : displaced);
149 :
150 1334 : const auto & secondary_set = _mortar_data.getHigherDimSubdomainIDs(_secondary_subdomain_id);
151 1334 : const auto & primary_set = _mortar_data.getHigherDimSubdomainIDs(_primary_subdomain_id);
152 :
153 2668 : std::set_union(secondary_set.begin(),
154 : secondary_set.end(),
155 : primary_set.begin(),
156 : primary_set.end(),
157 1334 : std::inserter(_higher_dim_subdomain_ids, _higher_dim_subdomain_ids.begin()));
158 1334 : _boundary_ids = {_secondary_id, _primary_id};
159 1334 : }
160 :
161 : void
162 279422 : MortarConsumerInterface::setNormals()
163 : {
164 279422 : if (interpolateNormals())
165 278226 : _normals = amg().getNormals(*_lower_secondary_elem, _qrule_face->get_points());
166 : else
167 1196 : _normals = amg().getNodalNormals(*_lower_secondary_elem);
168 279422 : }
169 :
170 : void
171 210240 : MortarConsumerInterface::trimDerivative(const dof_id_type remove_derivative_index,
172 : ADReal & dual_number)
173 : {
174 210240 : auto md_it = dual_number.derivatives().nude_data().begin();
175 210240 : auto mi_it = dual_number.derivatives().nude_indices().begin();
176 :
177 210240 : auto d_it = dual_number.derivatives().nude_data().begin();
178 :
179 210240 : for (auto i_it = dual_number.derivatives().nude_indices().begin();
180 592000 : i_it != dual_number.derivatives().nude_indices().end();
181 381760 : ++i_it, ++d_it)
182 381760 : if (*i_it != remove_derivative_index)
183 : {
184 356240 : *mi_it = *i_it;
185 356240 : *md_it = *d_it;
186 356240 : ++mi_it;
187 356240 : ++md_it;
188 : }
189 :
190 210240 : std::size_t n_indices = md_it - dual_number.derivatives().nude_data().begin();
191 210240 : dual_number.derivatives().nude_indices().resize(n_indices);
192 210240 : dual_number.derivatives().nude_data().resize(n_indices);
193 210240 : }
|
