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