Build for OpenFOAM v11

[efirvida]

The source is updated to be compiled for OF-11,

• Fixing: OpenFOAM 11 and higher does not ship an fvCFD.h #306
• Recovering the primitivePatchInterpolation.h which was removed [this may need some other changes ]
• and update the namespaces and some build options dependencies

[efirvida]

Still not working something fails with the tutorials but I don't know how to fix it: this is the error log for perpendicular-flap tutorial using CalculiX as solid solver. Maybe related with the primitivePatchInterpolation.h

*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 11
Exec   : foamRun
Date   : Jan 01 2024
Time   : 12:09:28
Host   : "localhost.localdomain"
PID    : 12013
I/O    : uncollated
Case   : /home/efirvida/Desktop/dev/Doctorado/testcases/tutorials/perpendicular-flap/fluid-openfoam
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Selecting fvMeshMover motionSolver
Selecting motion solver: displacementLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Selecting solver incompressibleFluid
Selecting viscosity model Newtonian
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
No MRF models present

Courant Number mean: 0 max: 0

PIMPLE: No convergence criteria found


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode



Starting time loop

---[preciceAdapter] Loaded the OpenFOAM-preCICE adapter - v1.2.3.
---[preciceAdapter] Reading preciceDict...
---[preciceAdapter] [DEBUG]   precice-config-file : ../precice-config.xml
---[preciceAdapter] [DEBUG]   participant name    : Fluid
---[preciceAdapter] [DEBUG]   modules requested   : 
---[preciceAdapter] [DEBUG]   - FSI

---[preciceAdapter] [DEBUG]   interfaces : 
---[preciceAdapter] [DEBUG]   - mesh         : Fluid-Mesh
---[preciceAdapter] [DEBUG]     locations    : faceCenters
---[preciceAdapter] [DEBUG]     connectivity : 0
---[preciceAdapter] [DEBUG]     patches      : 
---[preciceAdapter] [DEBUG]       - flap
---[preciceAdapter] [DEBUG]     writeData    : 
---[preciceAdapter] [DEBUG]       - Force
---[preciceAdapter] [DEBUG]     readData     : 
---[preciceAdapter] [DEBUG]       - Displacement
---[preciceAdapter] [DEBUG] Configuring the FSI module...
---[preciceAdapter] [DEBUG]     user-defined solver type : 
---[preciceAdapter] [DEBUG]     restart from deformed : 1
---[preciceAdapter] [DEBUG]     pointDisplacement field name : pointDisplacement
---[preciceAdapter] [DEBUG]     cellDisplacement field name : cellDisplacement
---[preciceAdapter] [DEBUG]     force field name : Force
---[preciceAdapter] [DEBUG] Determining the solver type for the FSI module... (override by setting solverType to one of {compressible, incompressible, solid})
---[preciceAdapter] [DEBUG] Automatically determined solver type : incompressible
---[preciceAdapter] [DEBUG] Checking the timestep type (fixed vs adjustable)...
---[preciceAdapter] [DEBUG]   Timestep type: fixed.
---[preciceAdapter] [DEBUG] Creating the preCICE solver interface...
---[preciceAdapter] [DEBUG]   Number of processes: 1
---[preciceAdapter] [DEBUG]   MPI rank: 0
---[precice]  This is preCICE version 2.5.0
---[precice]  Revision info: no-info [git failed to run]
---[precice]  Build type: Release (without debug log)
---[precice]  Configuring preCICE with configuration "../precice-config.xml"
---[precice]  I am participant "Fluid"
---[preciceAdapter] [DEBUG]   preCICE solver interface was created.
---[preciceAdapter] [DEBUG] Creating interfaces...
---[preciceAdapter] [DEBUG] Number of face centres: 33
---[preciceAdapter] [DEBUG] Interface created on mesh Fluid-Mesh
---[preciceAdapter] [DEBUG] Adding coupling data writers...
---[preciceAdapter] [DEBUG] Added writer: Force.
---[preciceAdapter] [DEBUG] Adding coupling data readers...
---[preciceAdapter] [DEBUG] Added reader: Displacement.
---[preciceAdapter] [DEBUG] Initializing the preCICE solver interface...
---[precice]  Setting up primary communication to coupling partner/s
---[precice]  Primary ranks are connected
---[precice]  Setting up preliminary secondary communication to coupling partner/s
---[precice]  Prepare partition for mesh Fluid-Mesh
---[precice]  Receive global mesh Solid-Mesh
---[precice]  Setting up secondary communication to coupling partner/s
---[precice]  Secondary ranks are connected
---[precice]  iteration: 1 of 50, time-window: 1, time: 0 of 5, time-window-size: 0.01, max-timestep-length: 0.01, ongoing: yes, time-window-complete: no, write-iteration-checkpoint 
---[preciceAdapter] [DEBUG] Initializing preCICE data...
---[precice]  Compute "write" mapping from mesh "Fluid-Mesh" to mesh "Solid-Mesh".
---[precice]  Using tree-based preallocation for matrix C
---[precice]  Using tree-based preallocation for matrix A
---[precice]  Using rescaling. Solver converged after 363 of 10000 iterations due to sufficient relative convergence. Last residual norm: 1.07156e-08, limits: relative 1.57162e-08 (rtol 1e-09), absolute 1e-50, divergence 1.57162e+31(dtol 1e+30)
---[precice]  Mapping "Force" conservative from "Fluid-Mesh" (ID 0) to "Solid-Mesh" (ID 1) for dimension 0 with polynomial set to separate
---[precice]  Mapping "Force" conservative from "Fluid-Mesh" (ID 0) to "Solid-Mesh" (ID 1) for dimension 1 with polynomial set to separate
---[precice]  initializeData is skipped since no data has to be initialized.
---[preciceAdapter] preCICE was configured and initialized
---[preciceAdapter] [DEBUG] Reading coupling data...
---[preciceAdapter] [DEBUG] Adding in checkpointed fields...
---[preciceAdapter] [DEBUG] Checkpoint nu : volScalarField
---[preciceAdapter] [DEBUG] Checkpoint p : volScalarField
---[preciceAdapter] [DEBUG] Checkpoint Cc : volVectorField
---[preciceAdapter] [DEBUG] Checkpoint Force : volVectorField
---[preciceAdapter] [DEBUG] Checkpoint U : volVectorField
---[preciceAdapter] [DEBUG] Checkpoint cellDisplacement : volVectorField
---[preciceAdapter] [DEBUG] Checkpoint magSf : surfaceScalarField
---[preciceAdapter] [DEBUG] Checkpoint phi : surfaceScalarField
---[preciceAdapter] [DEBUG] Checkpoint Cf : surfaceVectorField
---[preciceAdapter] [DEBUG] Checkpoint Sf : surfaceVectorField
---[preciceAdapter] [DEBUG] Checkpoint pointDisplacement : pointVectorField
---[preciceAdapter] [DEBUG] Writing a checkpoint...
---[preciceAdapter] [DEBUG] Stored time value t = 0.000000
---[preciceAdapter] [DEBUG] Storing mesh points...
---[preciceAdapter] [DEBUG] Stored mesh points.
---[preciceAdapter] [DEBUG] Checkpoint for time t = 0.000000 was stored.
---[preciceAdapter] [DEBUG] Adjusting the solver's timestep...
---[preciceAdapter] [DEBUG] The solver's timestep is the same as the coupling timestep.
---[preciceAdapter] Setting the solver's endTime to infinity to prevent early exits. Only preCICE will control the simulation's endTime. Any functionObject's end() method will be triggered by the adapter. You may disable this behavior in the adapter's configuration.
Constructing face momentum Uf
Courant Number mean: 0 max: 0
Time = 0.01s

smoothSolver:  Solving for cellDisplacementx, Initial residual = 0, Final residual = 0, No Iterations 2
smoothSolver:  Solving for cellDisplacementy, Initial residual = 0, Final residual = 0, No Iterations 2
DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 0.000723746, No Iterations 68
DICPCG:  Solving for pcorr, Initial residual = 4.47129e-05, Final residual = 3.88864e-08, No Iterations 31
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 9.60943e-05, No Iterations 16
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 9.02537e-05, No Iterations 17
DICPCG:  Solving for p, Initial residual = 1, Final residual = 0.000833061, No Iterations 68
DICPCG:  Solving for p, Initial residual = 2.87207e-05, Final residual = 2.84167e-08, No Iterations 28
time step continuity errors : sum local = 7.82466e-09, global = -3.12757e-10
DICPCG:  Solving for p, Initial residual = 0.0125934, Final residual = 9.46941e-06, No Iterations 59
DICPCG:  Solving for p, Initial residual = 9.51191e-06, Final residual = 7.10002e-09, No Iterations 60
time step continuity errors : sum local = 1.95718e-09, global = 3.48354e-11
DICPCG:  Solving for p, Initial residual = 0.00839243, Final residual = 8.29905e-06, No Iterations 46
DICPCG:  Solving for p, Initial residual = 8.29111e-06, Final residual = 8.50421e-09, No Iterations 59
time step continuity errors : sum local = 2.35048e-09, global = 6.10765e-11
DICPCG:  Solving for p, Initial residual = 0.00691712, Final residual = 5.70749e-06, No Iterations 59
DICPCG:  Solving for p, Initial residual = 5.71007e-06, Final residual = 8.97883e-09, No Iterations 59
time step continuity errors : sum local = 2.48866e-09, global = -7.63994e-11, cumulative = -7.63994e-11
ExecutionTime = 0.596885 s  ClockTime = 1 s

---[preciceAdapter] [DEBUG] Writing coupling data...
---[preciceAdapter] [DEBUG] Advancing preCICE...
---[precice]  Mapping "Force" conservative from "Fluid-Mesh" (ID 0) to "Solid-Mesh" (ID 1) for dimension 0 with polynomial set to separate
---[precice]  Mapping "Force" conservative from "Fluid-Mesh" (ID 0) to "Solid-Mesh" (ID 1) for dimension 1 with polynomial set to separate
---[precice]  Compute "read" mapping from mesh "Solid-Mesh" to mesh "Fluid-Mesh".
---[precice]  Using tree-based preallocation for matrix C
---[precice]  Using tree-based preallocation for matrix A
---[precice]  Using rescaling. Solver converged after 363 of 10000 iterations due to sufficient relative convergence. Last residual norm: 1.07156e-08, limits: relative 1.57162e-08 (rtol 1e-09), absolute 1e-50, divergence 1.57162e+31(dtol 1e+30)
---[precice]  Mapping "Displacement" consistent from "Solid-Mesh" (ID 1) to "Fluid-Mesh" (ID 0) for dimension 0 with polynomial set to separate
---[precice]  Mapping "Displacement" consistent from "Solid-Mesh" (ID 1) to "Fluid-Mesh" (ID 0) for dimension 1 with polynomial set to separate
---[precice]  iteration: 2 of 50, time-window: 1, time: 0 of 5, time-window-size: 0.01, max-timestep-length: 0.01, ongoing: yes, time-window-complete: no, read-iteration-checkpoint 
---[preciceAdapter] [DEBUG] Reading a checkpoint...
---[preciceAdapter] [DEBUG] Reloaded time value t = 0.000000
---[preciceAdapter] [DEBUG] Moving mesh points to their previous locations...
smoothSolver:  Solving for cellDisplacementx, Initial residual = 0, Final residual = 0, No Iterations 2
smoothSolver:  Solving for cellDisplacementy, Initial residual = 0, Final residual = 0, No Iterations 2
---[preciceAdapter] [DEBUG] Moved mesh points to their previous locations.
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see https://petsc.org/release/faq/#valgrind and https://petsc.org/release/faq/
[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run 
[0]PETSC ERROR: to get more information on the crash.
[0]PETSC ERROR: Run with -malloc_debug to check if memory corruption is causing the crash.