[kimdh928]

Hello,

foam-extend 4.1
solids4Foam 2.1

Some of my FSI simulations crashed, and I thought it was caused by overestimation of fluid pressure at the fluid-solid interface.

I searched for solution to fix my problem and looked at the solids4Foam source code, and I found a BC named elasticWallPressure and relevant paper as below.

Tuković, Ž., Bukač, M., Cardiff, P., Jasak, H., & Ivanković, A. (2019, January). Added mass partitioned fluid–structure interaction solver based on a robin boundary condition for pressure. In OpenFOAM®: Selected Papers of the 11th Workshop (pp. 1-22). Cham: Springer International Publishing.

And I also found that elasticWallPressure BC was tested on 3dTube tutorials, so I expected it would work for my case.

So I imposed elasticWallPressure BC on elastically moving boundary in fluid domain, and the serial run was successful for short times.

But I usually run simulations in parallel because serial run is very time-consuming, so I tried parallel run and it crashed with error messages below.

[dohun@fpil02 elastic_pressure_wall_test]$ mpirun -np 4 solids4Foam -parallel
--------------------------------------------------------------------------
WARNING: No preset parameters were found for the device that Open MPI
detected:

  Local host:            fpil02
  Device name:           mlx5_0
  Device vendor ID:      0x02c9
  Device vendor part ID: 4119

Default device parameters will be used, which may result in lower
performance.  You can edit any of the files specified by the
btl_openib_device_param_files MCA parameter to set values for your
device.

NOTE: You can turn off this warning by setting the MCA parameter
      btl_openib_warn_no_device_params_found to 0.
--------------------------------------------------------------------------
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.1                                |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  | For copyright notice see file Copyright         |
\*---------------------------------------------------------------------------*/
Build  : 4.1-70b064d0f326
Exec   : solids4Foam -parallel
Date   : Sep 30 2025
Time   : 11:55:13
Host   : "fpil02"
PID    : 1374960
CtrlDict : "/home/dohun/case/fpil_ready/thesis/flexible/elastic_pressure_wall_test/system/controlDict"
Case   : /home/dohun/case/fpil_ready/thesis/flexible/elastic_pressure_wall_test
nProcs : 4
Slaves :
3
(
"fpil02.1374961"
"fpil02.1374962"
"fpil02.1374963"
)

Pstream initialized with:
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

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

/*---------------------------------------------------------------------------*\
|    For further information on the solids4Foam toolbox implementations,      |
|    please see the following publications:                                   |
|                                                                             |
|    P. Cardiff, A Karac, P. De Jaeger, H. Jasak, J. Nagy, A. Ivankovic,      |
|    Z. Tukovic: An open-source finite volume toolbox for solid mechanics and |
|    fluid-solid interaction simulations. arXiv:1808.10736v2, 2018, available |
|    at https://arxiv.org/abs/1808.10736.                                     |
|                                                                             |
|    Z. Tukovic, A. Karac, P. Cardiff, H. Jasak, A. Ivankovic: OpenFOAM       |
|    finite volume solver for fluid-solid interaction.  Transactions of       |
|    Famena, 42 (3), pp. 1-31, 2018, 10.21278/TOF.42301.                      |
\*---------------------------------------------------------------------------*/

Selecting physicsModel fluidSolidInterface
Selecting fluidSolidInterface Aitken

Selecting fluidModel pimpleFluid
--> FOAM Warning :
    From function void* Foam::dlOpen(const Foam::fileName&, bool)
    in file POSIX.C at line 1204
    dlopen error : libmesquiteMotionSolver.so: cannot open shared object file: No such file or directory
--> FOAM Warning :
    From function bool Foam::dlLibraryTable::open(const Foam::fileName&, bool)
    in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 97
    could not load "libmesquiteMotionSolver.so"
--> FOAM Warning :
    From function dynamicFvMesh::New(const IOobject& io)
    in file dynamicFvMesh/newDynamicFvMesh.C at line 46
    Loading of dynamic mesh library libmesquiteMotionSolver.so unsuccesful. Some dynamic mesh  methods may not be  available
Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: velocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
g field not found in constant directory: initialising to zero
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}


PIMPLE: max iterations = 50
    field U     : relTol 0, tolerance 1e-05
    field p     : relTol 0, tolerance 0.0005

Selecting solidModel nonLinearGeometryTotalLagrangianTotalDisplacement
--> FOAM Warning :
    From function void* Foam::dlOpen(const Foam::fileName&, bool)
    in file POSIX.C at line 1204
    dlopen error : libmesquiteMotionSolver.so: cannot open shared object file: No such file or directory
--> FOAM Warning :
    From function bool Foam::dlLibraryTable::open(const Foam::fileName&, bool)
    in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 97
    could not load "libmesquiteMotionSolver.so"
--> FOAM Warning :
    From function dynamicFvMesh::New(const IOobject& io)
    in file dynamicFvMesh/newDynamicFvMesh.C at line 46
    Loading of dynamic mesh library libmesquiteMotionSolver.so unsuccesful. Some dynamic mesh  methods may not be  available
Selecting dynamicFvMesh staticFvMesh
Creating solidTraction boundary condition
Creating fixedDisplacement boundary condition
    under-relaxation method: fixed
Creating the mechanicalModel
Selecting mechanical law neoHookeanElastic
additionalMeshCorrection: yes
Selecting interfaceToInterfaceMapping directMap
Time = 0.002

Setting traction on solid interfaces
Mapping face values using directMap
Total force on fluid interface 0: (-3.11837e-21 -0.00164846 1.49642e-21)
Total force on solid interface 0: (3.11837e-21 0.00164846 -2.44228e-21)

Evolving solid solver
Solving the total Lagrangian form of the momentum equation for D
setCellDisplacements: reading cellDisplacements
    Corr, res, relRes, matRes, iters
    Both residuals have converged
    80, 9.17692e-07, 3.0739e-07, 0, 0

Mapping point values using directMap
Mapping point values using directMap
Mapping point values using directMap
FSI relative residual1 norm for interface 0: 1
FSI residual2 norm for interface 0: 1

Time = 0.002, iteration: 1
Current fsi under-relaxation factor (fixed): 0.1
Setting acceleration at fluid side of the interface
Mapping face values using directMap
Maximal accumulated displacement of interface 0: 3.81251e-06
GAMG:  Solving for cellMotionUx, Initial residual = 1, Final residual = 0.000563794, No Iterations 2
GAMG:  Solving for cellMotionUy, Initial residual = 1, Final residual = 0.000391021, No Iterations 2
GAMG:  Solving for cellMotionUz, Initial residual = 1, Final residual = 0.000549124, No Iterations 2
Evolving fluid model: pimpleFluid
Courant Number mean: 0.0358118 max: 0.151555 velocity magnitude: 0.3
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 1.86162e-05, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 2.89152e-05, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.000915224, No Iterations 1
[fpil02:1374960] *** Process received signal ***
[fpil02:1374960] Signal: Segmentation fault (11)
[fpil02:1374960] Signal code:  (-6)
[fpil02:1374960] Failing at address: 0x3ed0014faf0
[fpil02:1374963] *** Process received signal ***
[fpil02:1374963] Signal: Segmentation fault (11)
[fpil02:1374963] Signal code:  (-6)
[fpil02:1374963] Failing at address: 0x3ed0014faf3
[fpil02:1374963] [ 0] /lib64/libc.so.6(+0x4eb50)[0x14a143f5fb50]
[fpil02:1374963] [ 1] /lib64/libc.so.6(gsignal+0x10f)[0x14a143f5facf]
[fpil02:1374963] [ 2] [fpil02:1374960] [ 0] /lib64/libc.so.6(+0x4eb50)[0x14a143f5fb50]
[fpil02:1374963] [ 3] /lib64/libc.so.6(+0x4eb50)[0x14bbc9cb9b50]
[fpil02:1374960] [ 1] /lib64/libc.so.6(gsignal+0x10f)[0x14bbc9cb9acf]
[fpil02:1374960] [ 2] /lib64/libc.so.6(+0x4eb50)[0x14bbc9cb9b50]
[fpil02:1374960] [ 3] /home/dohun/foam/foam-extend-4.1/lib/linux64GccDPInt32Opt/libfoam.so(_ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6_+0xb8)[0x14a14577a978]
[fpil02:1374963] [ 4] /home/dohun/foam/foam-extend-4.1/lib/linux64GccDPInt32Opt/libfoam.so(_ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6_+0xb8)[0x14bbcb4d4978]
[fpil02:1374960] [ 4] /home/dohun/foam/foam-extend-4.1/lib/linux64GccDPInt32Opt/libfoam.so(_ZN4FoammlERKNS_3tmpINS_5FieldIdEEEERKNS_5UListIdEE+0x52)[0x14a1457863c2]
[fpil02:1374963] [ 5] /home/dohun/foam/foam-extend-4.1/lib/linux64GccDPInt32Opt/libfoam.so(_ZN4FoammlERKNS_3tmpINS_5FieldIdEEEERKNS_5UListIdEE+0x52)[0x14bbcb4e03c2]
[fpil02:1374960] [ 5] /home/dohun/foam/dohun-4.1/lib/linux64GccDPInt32Opt/libsolids4FoamModels.so(_ZN4Foam37elasticWallPressureFvPatchScalarField12updateCoeffsEv+0x93d)[0x14a1466fc5ed]
[fpil02:1374963] [ 6] /home/dohun/foam/dohun-4.1/lib/linux64GccDPInt32Opt/libsolids4FoamModels.so(_ZN4Foam37elasticWallPressureFvPatchScalarField12updateCoeffsEv+0x93d)[0x14bbcc4565ed]
[fpil02:1374960] [ 6] /home/dohun/foam/dohun-4.1/lib/linux64GccDPInt32Opt/libsolids4FoamModels.so(_ZN4Foam8fvMatrixIdEC2ERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS_12dimensionSetE+0x3e2)[0x14a1463ad412]
[fpil02:1374963] [ 7] /home/dohun/foam/dohun-4.1/lib/linux64GccDPInt32Opt/libsolids4FoamModels.so(_ZN4Foam8fvMatrixIdEC2ERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS_12dimensionSetE+0x3e2)[0x14bbcc107412]
[fpil02:1374960] [ 7] /home/dohun/foam/foam-extend-4.1/lib/linux64GccDPInt32Opt/libfiniteVolume.so(_ZN4Foam2fv20gaussLaplacianSchemeIddE23fvmLaplacianUncorrectedERKNS_14GeometricFieldIdNS_13fvsPatchFieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldENS_7volMeshEEE+0xbc)[0x14a140a1419c]
[fpil02:1374963] [ 8] /home/dohun/foam/foam-extend-4.1/lib/linux64GccDPInt32Opt/libfiniteVolume.so(_ZN4Foam2fv20gaussLaplacianSchemeIddE23fvmLaplacianUncorrectedERKNS_14GeometricFieldIdNS_13fvsPatchFieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldENS_7volMeshEEE+0xbc)[0x14bbc676e19c]
[fpil02:1374960] [ 8] /home/dohun/foam/foam-extend-4.1/lib/linux64GccDPInt32Opt/libfiniteVolume.so(_ZN4Foam2fv20gaussLaplacianSchemeIddE12fvmLaplacianERKNS_14GeometricFieldIdNS_13fvsPatchFieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldENS_7volMeshEEE+0x6b)[0x14a140a10f9b]
[fpil02:1374963] [ 9] /home/dohun/foam/foam-extend-4.1/lib/linux64GccDPInt32Opt/libfiniteVolume.so(_ZN4Foam2fv20gaussLaplacianSchemeIddE12fvmLaplacianERKNS_14GeometricFieldIdNS_13fvsPatchFieldENS_11surfaceMeshEEERKNS3_IdNS_12fvPatchFieldENS_7volMeshEEE+0x6b)[0x14bbc676af9b]
[fpil02:1374960] [ 9] /home/dohun/foam/foam-extend-4.1/lib/linux64GccDPInt32Opt/libfiniteVolume.so(_ZN4Foam2fv15laplacianSchemeIddE12fvmLaplacianERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEES8_+0x56)[0x14a140a123c6]
[fpil02:1374963] [10] /home/dohun/foam/foam-extend-4.1/lib/linux64GccDPInt32Opt/libfiniteVolume.so(_ZN4Foam2fv15laplacianSchemeIddE12fvmLaplacianERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEES8_+0x56)[0x14bbc676c3c6]
[fpil02:1374960] [10] /home/dohun/foam/dohun-4.1/lib/linux64GccDPInt32Opt/libsolids4FoamModels.so(_ZN4Foam3fvm9laplacianIddEENS_3tmpINS_8fvMatrixIT_EEEERKNS_14GeometricFieldIT0_NS_12fvPatchFieldENS_7volMeshEEERKNS7_IS4_S9_SA_EERKNS_4wordE+0x5f)[0x14a146493aef]
[fpil02:1374963] [11] /home/dohun/foam/dohun-4.1/lib/linux64GccDPInt32Opt/libsolids4FoamModels.so(_ZN4Foam11fluidModels11pimpleFluid6evolveEv+0xef3)[0x14a1464dfde3]
[fpil02:1374963] [12] /home/dohun/foam/dohun-4.1/lib/linux64GccDPInt32Opt/libsolids4FoamModels.so(_ZN4Foam3fvm9laplacianIddEENS_3tmpINS_8fvMatrixIT_EEEERKNS_14GeometricFieldIT0_NS_12fvPatchFieldENS_7volMeshEEERKNS7_IS4_S9_SA_EERKNS_4wordE+0x5f)[0x14bbcc1edaef]
[fpil02:1374960] [11] /home/dohun/foam/dohun-4.1/lib/linux64GccDPInt32Opt/libsolids4FoamModels.so(_ZN4Foam20fluidSolidInterfaces23AitkenCouplingInterface6evolveEv+0x46)[0x14a14650a3a6]
[fpil02:1374963] [13] solids4Foam[0x401e20]
[fpil02:1374963] [14] /home/dohun/foam/dohun-4.1/lib/linux64GccDPInt32Opt/libsolids4FoamModels.so(_ZN4Foam11fluidModels11pimpleFluid6evolveEv+0xef3)[0x14bbcc239de3]
[fpil02:1374960] [12] /lib64/libc.so.6(__libc_start_main+0xe5)[0x14a143f4bd85]
[fpil02:1374963] [15] solids4Foam[0x40224e]
[fpil02:1374963] *** End of error message ***
/home/dohun/foam/dohun-4.1/lib/linux64GccDPInt32Opt/libsolids4FoamModels.so(_ZN4Foam20fluidSolidInterfaces23AitkenCouplingInterface6evolveEv+0x46)[0x14bbcc2643a6]
[fpil02:1374960] [13] solids4Foam[0x401e20]
[fpil02:1374960] [14] /lib64/libc.so.6(__libc_start_main+0xe5)[0x14bbc9ca5d85]
[fpil02:1374960] [15] solids4Foam[0x40224e]
[fpil02:1374960] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 1374960 on node fpil02 exited on signal 11 (Segmentation fault).
--------------------------------------------------------------------------
[fpil02:1374958] 3 more processes have sent help message help-mpi-btl-openib.txt / no device params found
[fpil02:1374958] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

ChatGPT proposed me to change decomposition method, so I also tried simple and hierarchical decomposition methods as well as scotch, but all crashed.

Is there any way to run parallel FSI simulations with the use of elasticWallPressure boundary condition in solids4Foam?

Any help would be appreciated.