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 
197 }
198 
199 }
200 
201 }
202 
203 #endif // GEOS_PHYSICSSOLVERS_CONTACT_CONTACTFIELDS_HPP_
#define DECLARE_FIELD(NAME, KEY, TYPE, DEFAULT, PLOTLEVEL, RESTARTFLAG, DESCRIPTION)
Generates a traits struct.
Definition: MeshFields.hpp:39
@ Stick
element is closed: no jump across the discontinuity.