When the volume coupling case of the old version is replicated in the new version adapter, the data does not match during the writing process

New statement:

The old version adapter added an AITS module and two Locationtype:volume and volsurf to implement volume coupling, and provided tutorials for verification. I successfully added the AITS module and two Locationtype:volume and volsurf in the new adapter version, and compiled successfully, but when I ran the original tutorials, I encountered the following error (all the data exchanged in the case I ran was similar to velocity, and I used the new version adapter's Locationtype:volumeCentre and the old version adapter's Locationtype:volume and volsurf, and all had the same error):

tlk@tlk-virtual-machine:~/OpenFOAM/tlk-v2006/run/openfoam-adapter-1.3.1_ooo/try/volCoup-oneWay/flowSolve$ ./run.sh wmake ../../solvers/flowSolveFoam g++ -std=c++11 -m64 -pthread -DOPENFOAM=2006 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -Wno-unknown-pragmas -O3 -DNoRepository -ftemplate-depth-100 -I/home/tlk/OpenFOAM/OpenFOAM-v2006/src/finiteVolume/lnInclude -I/home/tlk/OpenFOAM/OpenFOAM-v2006/src/meshTools/lnInclude -I/home/tlk/OpenFOAM/OpenFOAM-v2006/src/sampling/lnInclude -I/home/tlk/OpenFOAM/OpenFOAM-v2006/src/transportModels/compressible/lnInclude -I/home/tlk/OpenFOAM/OpenFOAM-v2006/src/thermophysicalModels/basic/lnInclude -I/home/tlk/OpenFOAM/OpenFOAM-v2006/src/thermophysicalModels/radiation/lnInclude -I/home/tlk/OpenFOAM/OpenFOAM-v2006/src/TurbulenceModels/turbulenceModels/lnInclude -I/home/tlk/OpenFOAM/OpenFOAM-v2006/src/TurbulenceModels/compressible/lnInclude -iquote. -IlnInclude -I/home/tlk/OpenFOAM/OpenFOAM-v2006/src/OpenFOAM/lnInclude -I/home/tlk/OpenFOAM/OpenFOAM-v2006/src/OSspecific/POSIX/lnInclude -fPIC -Xlinker --add-needed -Xlinker --no-as-needed Make/linux64GccDPInt32Opt/flowSolveFoam.o -L/home/tlk/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/lib \ -lfiniteVolume -lfvOptions -lmeshTools -lsampling -lcompressibleTransportModels -lfluidThermophysicalModels -lradiationModels -lspecie -lturbulenceModels -lcompressibleTurbulenceModels -latmosphericModels -lOpenFOAM -ldl \ -lm -o /home/tlk/OpenFOAM/tlk-v2006/platforms/linux64GccDPInt32Opt/bin/flowSolveFoam Preparing and running the flowSolve participant... /---------------------------------------------------------------------------\ | ========= | | | \ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \ / O peration | Version: v2006 | | \ / A nd | Website: www.openfoam.com | | \/ M anipulation | | *---------------------------------------------------------------------------*/ Build : v2006 OPENFOAM=2006 Arch : "LSB;label=32;scalar=64" Exec : blockMesh Date : Oct 09 2024 Time : 22:42:58 Host : tlk-virtual-machine PID : 20745 I/O : uncollated Case : /home/tlk/OpenFOAM/tlk-v2006/run/openfoam-adapter-1.3.1_ooo/try/volCoup-oneWay/flowSolve nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations

// * // Create time

Creating block mesh from "system/blockMeshDict" Creating block edges No non-planar block faces defined Creating topology blocks Creating topology patches

Creating block mesh topology --> FOAM Warning : From Foam::polyMesh::polyMesh(const Foam::IOobject&, Foam::pointField&&, const cellShapeList&, const faceListList&, const wordList&, const Foam::PtrList&, const Foam::word&, const Foam::word&, bool) in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 915 Found 6 undefined faces in mesh; adding to default patch defaultFaces

Check topology

Basic statistics
    Number of internal faces : 2
    Number of boundary faces : 14
    Number of defined boundary faces : 14
    Number of undefined boundary faces : 0
Checking patch -> block consistency

Creating block offsets Creating merge list (topological search)...

Creating polyMesh from blockMesh Creating patches Creating cells Creating points with scale 1 Block 0 cell size : i : 0.0123864236173 .. 0.00247728472346 j : 0.00232699027775 .. 0.0349048541662 k : 0.4 .. 0.4

Block 1 cell size :
    i : 0.0024958960536 .. 0.012479480268
    j : 0.00232699027775 .. 0.0349048541662
    k : 0.4 .. 0.4

Block 2 cell size :
    i : 0.0392156862745 .. 0.0392156862745
    j : 0.00232699027775 .. 0.0349048541662
    k : 0.4 .. 0.4

There are no merge patch pairs

Writing polyMesh with 0 cellZones

Mesh Information

boundingBox: (-0.5 0 0) (3 0.5 0.4) nPoints: 24696 nCells: 12013 nFaces: 48386 nInternalFaces: 23692

Patches

patch 0 (start: 23692 size: 41) name: inlet patch 1 (start: 23733 size: 41) name: outlet patch 2 (start: 23774 size: 293) name: top patch 3 (start: 24067 size: 81) name: slip-bottom patch 4 (start: 24148 size: 51) name: bottom patch 5 (start: 24199 size: 161) name: interface patch 6 (start: 24360 size: 24026) name: defaultFaces

End

/---------------------------------------------------------------------------\ | ========= | | | \ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \ / O peration | Version: v2006 | | \ / A nd | Website: www.openfoam.com | | \/ M anipulation | | *---------------------------------------------------------------------------*/ Build : v2006 OPENFOAM=2006 Arch : "LSB;label=32;scalar=64" Exec : checkMesh Date : Oct 09 2024 Time : 22:42:59 Host : tlk-virtual-machine PID : 20747 I/O : uncollated Case : /home/tlk/OpenFOAM/tlk-v2006/run/openfoam-adapter-1.3.1_ooo/try/volCoup-oneWay/flowSolve nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations

// * // Create time

Create mesh for time = 0

Time = 0

Mesh stats points: 24696 internal points: 0 faces: 48386 internal faces: 23692 cells: 12013 faces per cell: 6 boundary patches: 7 point zones: 0 face zones: 0 cell zones: 0

Overall number of cells of each type: hexahedra: 12013 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0

Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology
inlet 41 84 ok (non-closed singly connected)
outlet 41 84 ok (non-closed singly connected)
top 293 588 ok (non-closed singly connected)
slip-bottom 81 164 ok (non-closed singly connected)
bottom 51 104 ok (non-closed singly connected)
interface 161 324 ok (non-closed singly connected)
defaultFaces 24026 24696 ok (non-closed singly connected)

Checking faceZone topology for multiply connected surfaces... No faceZones found.

Checking basic cellZone addressing... No cellZones found.

Checking geometry... Overall domain bounding box (-0.5 0 0) (3 0.5 0.4) Mesh has 2 geometric (non-empty/wedge) directions (1 1 0) Mesh has 2 solution (non-empty) directions (1 1 0) All edges aligned with or perpendicular to non-empty directions. Boundary openness (-5.95502088768e-18 -1.17546934044e-16 -2.36226616626e-16) OK. Max cell openness = 2.13849956364e-16 OK. Max aspect ratio = 16.8525355069 OK. Minimum face area = 5.76461746671e-06. Maximum face area = 0.015686274512. Face area magnitudes OK. Min volume = 2.30584698668e-06. Max volume = 0.000547527124249. Total volume = 0.7. Cell volumes OK. Mesh non-orthogonality Max: 0 average: 0 Non-orthogonality check OK. Face pyramids OK. Max skewness = 5.82550106883e-14 OK. Coupled point location match (average 0) OK.

Mesh OK.

End

/---------------------------------------------------------------------------\ | ========= | | | \ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \ / O peration | Version: v2006 | | \ / A nd | Website: www.openfoam.com | | \/ M anipulation | | *---------------------------------------------------------------------------*/ Build : v2006 OPENFOAM=2006 Arch : "LSB;label=32;scalar=64" Exec : flowSolveFoam Date : Oct 09 2024 Time : 22:43:00 Host : tlk-virtual-machine PID : 20752 I/O : uncollated Case : /home/tlk/OpenFOAM/tlk-v2006/run/openfoam-adapter-1.3.1_ooo/try/volCoup-oneWay/flowSolve nProcs : 1 trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20) allowSystemOperations : Allowing user-supplied system call operations

// * // Create time

Create mesh for time = 0

PIMPLE: Operating solver in PISO mode

Reading thermophysical properties

Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState rhoConst; specie specie; energy sensibleEnthalpy; }

--> FOAM Warning : From bool Foam::GeometricField<Type, PatchField, GeoMesh>::readIfPresent() [with Type = double; PatchField = Foam::fvPatchField; GeoMesh = Foam::volMesh] in file /home/tlk/OpenFOAM/OpenFOAM-v2006/src/OpenFOAM/lnInclude/GeometricField.C at line 108 read option IOobject::MUST_READ or MUST_READ_IF_MODIFIED suggests that a read constructor for field rho would be more appropriate. Reading field U

Creating field Dpdt

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type laminar Selecting laminar stress model Stokes

Reading g

Reading hRef Calculating field g.h

Reading field p_rgh

Creating field dpdt

Creating field kinetic energy K

No MRF models present

Radiation model not active: radiationProperties not found Selecting radiationModel none No finite volume options present Courant Number mean: 0 max: 0

Starting time loop

---[preciceAdapter] Loaded the OpenFOAM-preCICE adapter - v1.3.1. ---[preciceAdapter] Reading preciceDict... SolverType: ||1---[precice] This is preCICE version 3.1.1 ---[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 "flowSolve" ---[preciceAdapter] Internal faces = 23692 ---[preciceAdapter] Patch size = 7 ---[preciceAdapter] Boundary faces = 24694 ---[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 flowSolve-Mesh-vol ---[precice] Gather mesh flowSolve-Mesh-vol ---[precice] Send global mesh flowSolve-Mesh-vol ---[precice] Prepare partition for mesh flowSolve-Mesh-volSurf ---[precice] Gather mesh flowSolve-Mesh-volSurf ---[precice] Send global mesh flowSolve-Mesh-volSurf ---[precice] Setting up secondary communication to coupling partner/s ---[precice] Secondary ranks are connected ---[precice] iteration: 1, time-window: 1, time: 0 of 10, time-window-size: 0.01, max-time-step-size: 0.01, ongoing: yes, time-window-complete: no, ---[preciceAdapter] preCICE was configured and initialized ---[preciceAdapter] [DEBUG] Writing data for interface 0 ---[preciceAdapter] [DEBUG] Writing data for interface 1 ---[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. Courant Number mean: 0 max: 0 Time = 0.01

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 PIMPLE: iteration 1 DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 7.78951375942e-08, No Iterations 2 DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 2.6959134457e-07, No Iterations 2 DICPCG: Solving for p_rgh, Initial residual = 0.815645845211, Final residual = 0.00751038250855, No Iterations 89 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 6.87802874908e-05, global = -1.0849523174e-05, cumulative = -1.0849523174e-05 DICPCG: Solving for p_rgh, Initial residual = 0.309794817979, Final residual = 8.7469094214e-09, No Iterations 215 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 2.43039149227e-10, global = 3.24487586703e-13, cumulative = -1.08495228495e-05 ExecutionTime = 0.77 s ClockTime = 3 s

---[preciceAdapter] [DEBUG] Write: first value of U at bufferIndex: 0 is 0.000000 ---[precice] ERROR: Input sizes are inconsistent attempting to write 3D data "Velocity" to mesh "flowSolve-Mesh-vol". You passed 36707 vertices and 38043 data components, but we expected 110121 data components (3 x 36707). Looking for any time directories without results (e.g. stray functionObjectProperties files, see openfoam-adapter issue #26 on GitHub)...

precice / openfoam-adapter