GEOS
ContactFields.hpp
Go to the documentation of this file.
1 /*
2  * ------------------------------------------------------------------------------------------------------------
3  * SPDX-License-Identifier: LGPL-2.1-only
4  *
5  * Copyright (c) 2016-2024 Lawrence Livermore National Security LLC
6  * Copyright (c) 2018-2024 TotalEnergies
7  * Copyright (c) 2018-2024 The Board of Trustees of the Leland Stanford Junior University
8  * Copyright (c) 2023-2024 Chevron
9  * Copyright (c) 2019- GEOS/GEOSX Contributors
10  * All rights reserved
11  *
12  * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
13  * ------------------------------------------------------------------------------------------------------------
14  */
15 
20 #ifndef GEOS_PHYSICSSOLVERS_CONTACT_CONTACTFIELDS_HPP_
21 #define GEOS_PHYSICSSOLVERS_CONTACT_CONTACTFIELDS_HPP_
22 
23 #include "mesh/MeshFields.hpp"
25 
26 namespace geos
27 {
31 namespace fields
32 {
33 
34 namespace contact
35 {
36 
37 DECLARE_FIELD( totalBubbleDisplacement,
38  "totalBubbleDisplacement",
39  array2d< real64 >,
40  0,
41  LEVEL_0,
42  WRITE_AND_READ,
43  "Total bubble displacements at the faces" );
44 
45 DECLARE_FIELD( incrementalBubbleDisplacement,
46  "incrementalBubbleDisplacement",
47  array2d< real64 >,
48  0,
49  LEVEL_3,
50  WRITE_AND_READ,
51  "Incremental bubble displacements for the current time step on the nodes" );
52 
53 DECLARE_FIELD( iterativePenalty,
54  "iterativePenalty",
55  array2d< real64 >,
56  1.e5,
57  LEVEL_0,
58  WRITE_AND_READ,
59  "Penalty coefficients used in the iterative procedure of the Augmented Lagrangian Method" );
60 
61 DECLARE_FIELD( rotationMatrix,
62  "rotationMatrix",
63  array3d< real64 >,
64  0,
65  LEVEL_0,
66  WRITE_AND_READ,
67  "An array that holds the rotation matrices on the fracture" );
68 
69 DECLARE_FIELD( dispJump,
70  "displacementJump",
71  array2d< real64 >,
72  0,
73  LEVEL_0,
74  WRITE_AND_READ,
75  "Displacement jump vector in the local reference system" );
76 
77 DECLARE_FIELD( dispJump_n,
78  "displacementJump_n",
79  array2d< real64 >,
80  0,
81  NOPLOT,
82  WRITE_AND_READ,
83  "Displacement jump vector in the local reference system at the current time-step" );
84 
85 DECLARE_FIELD( slip,
86  "slip",
87  array1d< real64 >,
88  0,
89  LEVEL_0,
90  WRITE_AND_READ,
91  "Slip" );
92 
93 DECLARE_FIELD( tangentialTraction,
94  "tangentialTraction",
95  array1d< real64 >,
96  0,
97  LEVEL_0,
98  WRITE_AND_READ,
99  "Tangential traction" );
100 
101 DECLARE_FIELD( deltaSlip,
102  "deltaSlip",
103  array2d< real64 >,
104  0.0,
105  LEVEL_0,
106  WRITE_AND_READ,
107  "Slip increment" );
108 
109 DECLARE_FIELD( deltaSlip_n,
110  "deltaSlip_n",
111  array2d< real64 >,
112  0.0,
113  NOPLOT,
114  WRITE_AND_READ,
115  "Initial slip increment at this time step" );
116 
117 DECLARE_FIELD( deltaDispJump,
118  "deltaDisplacementJump",
119  array2d< real64 >,
120  0,
121  NOPLOT,
122  WRITE_AND_READ,
123  "Delta displacement jump vector" );
124 
125 DECLARE_FIELD( oldDispJump,
126  "oldDisplacementJump",
127  array2d< real64 >,
128  0,
129  NOPLOT,
130  WRITE_AND_READ,
131  "Displacement jump vector at the previous time-step" );
132 
133 DECLARE_FIELD( traction,
134  "traction",
135  array2d< real64 >,
136  0,
137  LEVEL_0,
138  WRITE_AND_READ,
139  "Fracture traction vector in the local reference system." );
140 
141 DECLARE_FIELD( traction_n,
142  "traction_n",
143  array2d< real64 >,
144  0,
145  NOPLOT,
146  WRITE_AND_READ,
147  "Initial fracture traction vector in the local reference system at this time-step." );
148 
149 DECLARE_FIELD( deltaTraction,
150  "deltaTraction",
151  array2d< real64 >,
152  0,
153  NOPLOT,
154  NO_WRITE,
155  "An array that holds the traction increments on the fracture." );
156 
157 DECLARE_FIELD( dTraction_dJump,
158  "dTraction_dJump",
159  array3d< real64 >,
160  0,
161  NOPLOT,
162  NO_WRITE,
163  "Derivative of the traction w.r.t. the displacement jump." );
164 
165 DECLARE_FIELD( dTraction_dPressure,
166  "dTraction_dPressure",
167  array1d< real64 >,
168  0,
169  NOPLOT,
170  NO_WRITE,
171  "Derivative of the traction w.r.t. to the fluid pressure." );
172 
173 DECLARE_FIELD( fractureState,
174  "fractureState",
175  array1d< integer >,
177  LEVEL_0,
178  WRITE_AND_READ,
179  "Fracture state." );
180 
181 DECLARE_FIELD( oldFractureState,
182  "oldFractureState",
183  array1d< integer >,
185  NOPLOT,
186  NO_WRITE,
187  "Fracture state at the previous timestep." );
188 
189 DECLARE_FIELD( targetIncrementalJump,
190  "targetIncrementalJump",
191  array2d< real64 >,
192  0,
193  NOPLOT,
194  WRITE_AND_READ,
195  "It's the target incremental jump in a timestep (e.g., slip coming from RS)." );
196 
198 
200  "yield0",
202  0,
203  NOPLOT,
204  WRITE_AND_READ,
205  "Accelerated variable (yield) at outer iteration v (two iterations ago)" );
206 
207 DECLARE_FIELD( yield1,
208  "yield1",
210  0,
211  NOPLOT,
212  WRITE_AND_READ,
213  "Accelerated variable (yield) at outer iteration v + 1 (previous iteration)" );
214 
215 DECLARE_FIELD( yield1_tilde,
216  "yield1_tilde",
218  0,
219  NOPLOT,
220  WRITE_AND_READ,
221  "Unaccelerated variable (yield) at outer iteration v + 1 (previous iteration)" );
222 
223 DECLARE_FIELD( yield2,
224  "yield2",
226  0,
227  NOPLOT,
228  WRITE_AND_READ,
229  "Accelerated variable (yield) at outer iteration v + 2 (current iteration)" );
230 
231 DECLARE_FIELD( yield2_tilde,
232  "yield2_tilde",
234  0,
235  NOPLOT,
236  WRITE_AND_READ,
237  "Unaccelerated variable (yield) at outer iteration v + 2 (current iteration)" );
238 
239 DECLARE_FIELD( relaxationFactor0,
240  "relaxationFactor0",
242  0,
243  NOPLOT,
244  WRITE_AND_READ,
245  "Old relaxation factor" );
246 
247 DECLARE_FIELD( relaxationFactor1,
248  "relaxationFactor1",
250  0,
251  NOPLOT,
252  WRITE_AND_READ,
253  "New relaxation factor" );
254 
255 }
256 
257 }
258 
259 }
260 
261 #endif // GEOS_PHYSICSSOLVERS_CONTACT_CONTACTFIELDS_HPP_
#define DECLARE_FIELD(NAME, KEY, TYPE, DEFAULT, PLOTLEVEL, RESTARTFLAG, DESCRIPTION)
Generates a traits struct.
Definition: MeshFields.hpp:39
Array< T, 1 > array1d
Alias for 1D array.
Definition: DataTypes.hpp:175
@ Stick
element is closed: no jump across the discontinuity.