[quote=“fsimonis, post:4, topic:2138, full:true”]
Hi
Can you explain the solution in more detail, I still don’t understand it now, I can probably know that there is an error when the adapter passes data to preCICE, but at the moment I don’t have the idea of how to locate the specific error and solve it so that it can run normally under the new adapter version
Blockquote
Therefore, this is either a bug in the OpenFOAM adapter @Makis, or the solver you are using is performing some kind of mesh refinement/coarsening that leads to this difference, which is unsupported.
Output when running tutorials：
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 10 2024
Time : 23:09:22
Host : tlk-virtual-machine
PID : 23374
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<Foam::dictionary>&, 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 10 2024
Time : 23:09:23
Host : tlk-virtual-machine
PID : 23398
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 10 2024
Time : 23:09:25
Host : tlk-virtual-machine
PID : 23408
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 = 1.14 s ClockTime = 4 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)...