-
-
Notifications
You must be signed in to change notification settings - Fork 83
/
Adapter.H
406 lines (280 loc) · 12.1 KB
/
Adapter.H
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
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
#ifndef PRECICEADAPTER_H
#define PRECICEADAPTER_H
#include "Interface.H"
// Conjugate Heat Transfer module
#include "CHT/CHT.H"
// Fluid-Structure Interaction module
#include "FSI/FSI.H"
// Fluid-Fluid module
#include "FF/FF.H"
// NOTE: If you want to couple a new variable, include your module's header here.
// You also need to include it in the Make/files file.
// In case you use additional OpenFOAM symbols, you may also need to specify
// the respective libraries in the Make/options.
// OpenFOAM header files
#include "Time.H"
#include "fvMesh.H"
// preCICE Solver Interface
#include <precice/precice.hpp>
namespace preciceAdapter
{
class Adapter
{
private:
//- Structure of the configuration of each coupling interface.
// Every interface needs to know the coupling mesh, the OpenFOAM
// patches that consist the coupling surface and the kinds
// of data that are exchanged.
struct InterfaceConfig
{
std::string meshName;
std::string locationsType;
bool meshConnectivity;
std::vector<std::string> patchNames;
std::vector<std::string> cellSetNames;
std::vector<std::string> writeData;
std::vector<std::string> readData;
};
//- Configuration interfaces
std::vector<struct InterfaceConfig> interfacesConfig_;
#ifdef ADAPTER_ENABLE_TIMINGS
//- Timers
clockValue timeInConfigRead_;
clockValue timeInPreciceConstruct_;
clockValue timeInMeshSetup_;
clockValue timeInInitialize_;
clockValue timeInCheckpointingSetup_;
clockValue timeInWrite_;
clockValue timeInAdvance_;
clockValue timeInRead_;
clockValue timeInCheckpointingRead_;
clockValue timeInCheckpointingWrite_;
clockValue timeInWriteResults_;
clockValue timeInFinalize_;
#endif
//- OpenFOAM runTime object
const Foam::Time& runTime_;
//- OpenFOAM fvMesh object
const Foam::fvMesh& mesh_;
// Configuration parameters used in the Adapter
//- Remember if there were errors in the read() method
bool errorsInConfigure = false;
//- preCICE configuration file name
Foam::word preciceConfigFilename_;
//- preCICE participant name
Foam::word participantName_;
//- Switch for checkpointing
bool checkpointing_ = false;
//- Type of timestep (fixed, adjustable)
bool adjustableTimestep_;
//- Should the (fixed) timestep be stored and used?
bool useStoredTimestep_ = false;
//- Switch to enable the ConjugateHeatTransfer module
bool CHTenabled_ = false;
//- Switch to enable the FluidStructureInteraction module
bool FSIenabled_ = false;
//- Switch to enable the FluidFluid module
bool FFenabled_ = false;
// NOTE: Add a switch for your new module here
//- Interfaces
std::vector<Interface*> interfaces_;
//- preCICE solver interface
precice::Participant* precice_ = nullptr;
//- preCICE solver interface initialized
bool preciceInitialized_ = false;
//- Conjugate Heat Transfer module object
CHT::ConjugateHeatTransfer* CHT_ = NULL;
//- Fluid-Structure Interaction module object
FSI::FluidStructureInteraction* FSI_ = NULL;
//- Fluid-Fluid module object
FF::FluidFluid* FF_ = NULL;
// NOTE: Add here a pointer for your new module object
// Timesteps
//- Timestep used by the solver
double timestepSolver_;
//- Stored (fixed) timestep
double timestepStored_;
// Checkpointing
//- Checkpointed time (value)
Foam::scalar couplingIterationTimeValue_;
//- Checkpointed time (index)
Foam::label couplingIterationTimeIndex_;
//- Checkpointed mesh points
Foam::pointField meshPoints_;
Foam::pointField oldMeshPoints_;
bool meshCheckPointed = false;
// TODO: Currently unused, see storeMeshPoints().
//- Checkpointed mesh volume
// bool oldVolsStored = false;
// Foam::volScalarField::Internal * oldVols_;
// Foam::volScalarField::Internal * oldOldVols_;
// int curTimeIndex_ = 0;
// Vectors of pointers to the checkpointed mesh fields and their copies
//- Checkpointed surfaceScalarField mesh fields
std::vector<Foam::surfaceScalarField*> meshSurfaceScalarFields_;
//- Checkpointed surfaceScalarField mesh fields (copies)
std::vector<Foam::surfaceScalarField*> meshSurfaceScalarFieldCopies_;
//- Checkpointed surfaceVectorField mesh fields
std::vector<Foam::surfaceVectorField*> meshSurfaceVectorFields_;
//- Checkpointed surfaceVectorField mesh fields (copies)
std::vector<Foam::surfaceVectorField*> meshSurfaceVectorFieldCopies_;
//- Checkpointed volVectorField mesh fields
std::vector<Foam::volVectorField*> meshVolVectorFields_;
//- Checkpointed volVectorField mesh fields (copies)
std::vector<Foam::volVectorField*> meshVolVectorFieldCopies_;
// TODO checkpoint for the V0 (Old volume) and V00 (Old-Old volume) fields.
//- Checkpointed volScalarField mesh fields
std::vector<Foam::volScalarField::Internal*> volScalarInternalFields_;
//- Checkpointed volScalarField mesh fields (copies)
std::vector<Foam::volScalarField::Internal*> volScalarInternalFieldCopies_;
// Vectors of pointers to the checkpointed fields and their copies
//- Checkpointed volScalarField fields
std::vector<Foam::volScalarField*> volScalarFields_;
//- Checkpointed volScalarField fields (copies)
std::vector<Foam::volScalarField*> volScalarFieldCopies_;
//- Checkpointed volVectorField fields
std::vector<Foam::volVectorField*> volVectorFields_;
//- Checkpointed volVectorField fields (copies)
std::vector<Foam::volVectorField*> volVectorFieldCopies_;
//- Checkpointed surfaceScalarField fields
std::vector<Foam::surfaceScalarField*> surfaceScalarFields_;
//- Checkpointed surfaceScalarField fields (copies)
std::vector<Foam::surfaceScalarField*> surfaceScalarFieldCopies_;
//- Checkpointed surfaceVectorField fields
std::vector<Foam::surfaceVectorField*> surfaceVectorFields_;
//- Checkpointed surfaceVectorField fields (copies)
std::vector<Foam::surfaceVectorField*> surfaceVectorFieldCopies_;
//- Checkpointed pointScalarField fields
std::vector<Foam::pointScalarField*> pointScalarFields_;
//- Checkpointed pointScalarField fields (copies)
std::vector<Foam::pointScalarField*> pointScalarFieldCopies_;
//- Checkpointed pointVectorField fields
std::vector<Foam::pointVectorField*> pointVectorFields_;
//- Checkpointed pointVectorField fields (copies)
std::vector<Foam::pointVectorField*> pointVectorFieldCopies_;
//- Checkpointed volTensorField fields
std::vector<Foam::volTensorField*> volTensorFields_;
//- Checkpointed volTensorField fields (copies)
std::vector<Foam::volTensorField*> volTensorFieldCopies_;
//- Checkpointed surfaceTensorField fields
std::vector<Foam::surfaceTensorField*> surfaceTensorFields_;
//- Checkpointed surfaceTensorField fields (copies)
std::vector<Foam::surfaceTensorField*> surfaceTensorFieldCopies_;
//- Checkpointed pointTensorField fields
std::vector<Foam::pointTensorField*> pointTensorFields_;
//- Checkpointed pointTensorField fields (copies)
std::vector<Foam::pointTensorField*> pointTensorFieldCopies_;
//- Checkpointed volSymmTensorField fields
std::vector<Foam::volSymmTensorField*> volSymmTensorFields_;
//- Checkpointed volSymmTensorField fields (copies)
std::vector<Foam::volSymmTensorField*> volSymmTensorFieldCopies_;
// NOTE: Declare additional vectors for any other types required.
// Configuration
//- Read the adapter's configuration file
bool configFileRead();
//- Check the adapter's configuration file
bool configFileCheck(const std::string adapterConfigFileName);
// Methods communicating with preCICE
//- Initialize preCICE and exchange the first data
void initialize();
//- Finalize and destroy preCICE
void finalize();
//- Advance preCICE
void advance();
//- Read the coupling data at each interface
void readCouplingData(double relativeReadTime);
//- Write the coupling data at each interface
void writeCouplingData();
//- Adjust the timestep of the solver according to preCICE
// and read data associated to the calculated time step length
void adjustSolverTimeStepAndReadData();
//- Determine if the coupling is still happening
bool isCouplingOngoing();
//- Determine if the coupling timestep has been completed
bool isCouplingTimeWindowComplete();
//- Determine if a checkpoint must be read
bool requiresReadingCheckpoint();
//- Determine if a checkpoint must be written
bool requiresWritingCheckpoint();
// Methods for checkpointing
//- Configure the mesh checkpointing
void setupMeshCheckpointing();
//- Configure the mesh checkpointing
void setupMeshVolCheckpointing();
//- Configure the checkpointing
void setupCheckpointing();
//- Make a copy of the runTime object
void storeCheckpointTime();
//- Restore the copy of the runTime object
void reloadCheckpointTime();
//- Store the locations of the mesh points
void storeMeshPoints();
//- Restore the locations of the mesh points
void reloadMeshPoints();
// Add mesh checkpoint fields, depending on the type
//- Add a surfaceScalarField mesh field
void addMeshCheckpointField(surfaceScalarField& field);
//- Add a surfaceVectorField mesh field
void addMeshCheckpointField(surfaceVectorField& field);
//- Add a volVectorField mesh field
void addMeshCheckpointField(volVectorField& field);
// TODO V0 and V00 checkpointed field.
//- Add the V0 and V00 checkpoint fields
void addVolCheckpointField(volScalarField::Internal& field);
// void addVolCheckpointFieldBuffer(volScalarField::Internal & field);
// Add checkpoint fields, depending on the type
//- Add a volScalarField to checkpoint
void addCheckpointField(volScalarField* field);
//- Add a volVectorField to checkpoint
void addCheckpointField(volVectorField* field);
//- Add a surfaceScalarField to checkpoint
void addCheckpointField(surfaceScalarField* field);
//- Add a surfaceVectorField to checkpoint
void addCheckpointField(surfaceVectorField* field);
//- Add a pointScalarField to checkpoint
void addCheckpointField(pointScalarField* field);
//- Add a pointVectorField to checkpoint
void addCheckpointField(pointVectorField* field);
// NOTE: Add here methods to add other object types to checkpoint,
// if needed.
//- Add a volTensorField to checkpoint
void addCheckpointField(volTensorField* field);
//- Add a surfaceTensorField to checkpoint
void addCheckpointField(surfaceTensorField* field);
//- Add a pointTensorField to checkpoint
void addCheckpointField(pointTensorField* field);
//- Add a volSymmTensorField to checkpoint
void addCheckpointField(volSymmTensorField* field);
//- Read the checkpoint - restore the mesh fields and time
void readMeshCheckpoint();
//- Read the checkpoint - restore the fields and time
void readCheckpoint();
//- Write the checkpoint - store the mesh fields and time
void writeMeshCheckpoint();
//- Write the checkpoint - store the fields and time
void writeCheckpoint();
// TODO Probably these can be included to the mesh checkpoints.
//- Read the volume checkpoint - restore the mesh volume fields
void readVolCheckpoint();
//- Write the volume checkpoint to a buffer - restore the mesh volume fields
void writeVolCheckpoint();
//- Destroy the preCICE interface and delete the allocated
// memory in a proper way. Called by the destructor.
void teardown();
public:
// Methods called by the functionObject
//- Constructor
Adapter(const Foam::Time& runTime, const Foam::fvMesh& mesh);
//- Setup the adapter's configuration
void configure();
//- Called by the functionObject's execute()
void execute();
//- Called by the functionObject's adjustTimeStep()
void adjustTimeStep();
//- Called by the functionObject's end()
void end();
//- Destructor
~Adapter();
};
}
#endif