[y99]

Originally Posted by bigphil

Hi,

There seems to be a compatibility issue between the MPI used by your PETSc installation and the MPI used by your OpenFOAM installation.

One solution is to install solids4foam with PETSc support. To do this, unset the PETSC_DIR variable, i.e.

> unset PETSC_DIR
> ./Allwclean
> ./Allwmake

Or, if you do want to use solids4foam with PETSc then I suggest installing OpenFOAM using the system MPI (set in foam-extend-4.1/etc/bashrc), and similarly ensure PETSc was installed using the system MPI.

Philip

Thank you Phillip for the answer, now I have compiled solids4Foam. I'd like to ask your advice again. I have another question in recent weeks. I added a modified fluid solver to fluidModels and it worked fine on its own. During the operation of the coupled solid solver, it runs well when U and rho are set to 0, but when U and rho are not set to 0, the solid solver will diverge in the first one or two steps of iteration. The output of running solids4Foan is left at the end. Do you have any suggestions?
In addition, I want to know how to transfer the force at the interface. Where should I find the source code of the force transfer at the fluid-solid interface?

Best regards
Yan




/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : _7bdb509494-20201222 OPENFOAM=2012
Arch : "LSB;label=32;scalar=64"
Exec : solids4Foam
Date : Apr 26 2023
Time : 17:46:17
Host : LAPTOP-KN3DB15U
PID : 4869
I/O : uncollated
Case : /mnt/d/Desktop/FSIiswtest
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

--> FOAM Warning :
From static Foam::IOstreamOption::compressionType Foam::IOstreamOption::compressionEnum(const Foam::word&, Foam::IOstreamOption::compressionType)
in file db/IOstreams/IOstreams/IOstreamOption.C at line 115
Unknown compression specifier 'uncompressed', using compression off
/*---------------------------------------------------------------------------*\
| 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 iswFluid
Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: velocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Selecting patchDistMethod meshWave
Reading g from constant directory
No finite volume options present

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSSTdensity
RAS
{
RASModel kOmegaSSTdensity;
turbulence on;
printCoeffs on;
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;
decayControl false;
kInf 0;
omegaInf 0;
}


PIMPLE: Operating solver in PISO mode

Resetting the dimensions of p
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.62458e-07, No Iterations 305
time step continuity errors : sum local = 6.91646e-12, global = -1.53114e-22, cumulative = -1.53114e-22
Courant Number mean: 0.000205662 max: 0.00409088
Selecting solidModel linearGeometryTotalDisplacement
Selecting dynamicFvMesh staticFvMesh
Creating solidTraction boundary condition
Creating fixedDisplacement boundary condition
under-relaxation method: fixed
Creating the mechanicalModel
Selecting mechanical law linearElastic
additionalMeshCorrection: false
Selecting interfaceToInterfaceMapping AMI
Time = 0.001

Setting traction on solid interfaces
Interpolating face values using AMI
Create AMI zone-to-zone interpolator
AMI: Creating addressing and weights between 170 source faces and 170 target faces
AMI: Patch source sum(weights) min:1 max:1 average:1
AMI: Patch target sum(weights) min:1 max:1 average:1
interface-to-interface face error: 8.89911e-16
Total force on fluid interface 0: (-1.38556e-13 0 -395.214)
Total force on solid interface 0: (1.38556e-13 0 395.214)

Evolving solid solver
Solving the momentum equation for D
setCellDisplacements: reading cellDisplacements
Corr, res, relRes, matRes, iters
Both residuals have converged
34, 7.86847e-07, 6.01779e-07, 0, 1

Interpolating point values using AMI
zoneA point orientation (< 0), max: -0.707107, min: -1, nIncorrectPoints: 0/342
Interpolating point values using AMI
FSI relative residual1 norm for interface 0: 1
FSI residual2 norm for interface 0: 1

Time = 0.001, iteration: 1
Current fsi under-relaxation factor (fixed): 0.4
Maximal accumulated displacement of interface 0: 0.0142375
GAMG: Solving for cellMotionUx, Initial residual = 1, Final residual = 2.27741e-17, No Iterations 1
GAMG: Solving for cellMotionUz, Initial residual = 1, Final residual = 6.59782e-17, No Iterations 1
Evolving fluid model: iswFluid
PIMPLE: iteration 1
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 255.079, No Iterations 1000
time step continuity errors : sum local = 0.012809, global = -2.75858e-06, cumulative = -2.75858e-06
DICPCG: Solving for p_rgh, Initial residual = 0.992327, Final residual = 0.0893359, No Iterations 118
time step continuity errors : sum local = 0.00151754, global = -2.75862e-06, cumulative = -5.5172e-06
GAMG: Solving for p_rgh, Initial residual = 0.932017, Final residual = 3.056e+124, No Iterations 1000
time step continuity errors : sum local = 4.97476e+121, global = 3.0744e+101, cumulative = 3.0744e+101
smoothSolver: Solving for omega, Initial residual = 0.00928351, Final residual = 7.20246e-09, No Iterations 14
bounding omega, min: -3.5626 max: 12.9472 average: 0.36269
smoothSolver: Solving for k, Initial residual = 1, Final residual = 8.92327e-09, No Iterations 20
bounding k, min: -1.19153e-06 max: 5.34667e-06 average: 9.90071e-07
smoothSolver: Solving for rho, Initial residual = 0.769647, Final residual = 6.66153e-09, No Iterations 31
Setting traction on solid interfaces
Interpolating face values using AMI
Total force on fluid interface 0: (-6.0442e+129 -1.68682e+100 7.27115e+136)
Total force on solid interface 0: (6.0442e+129 1.68682e+100 -7.27115e+136)

Evolving solid solver
Solving the momentum equation for D
Corr, res, relRes, matRes, iters
Both residuals have converged
39, 2.61483e-07, 7.98765e-07, 0, 1

Interpolating point values using AMI
Interpolating point values using AMI
FSI relative residual1 norm for interface 0: 1
FSI residual2 norm for interface 0: 1

Time = 0.001, iteration: 2
Current fsi under-relaxation factor (fixed): 0.4
Maximal accumulated displacement of interface 0: 2.61942e+132
GAMG: Solving for cellMotionUx, Initial residual = 1, Final residual = 2.84138e-17, No Iterations 1
GAMG: Solving for cellMotionUz, Initial residual = 1, Final residual = 1.31025e-17, No Iterations 1
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libpthread.so.0
#3 Foam::face::centre(Foam::UList<Foam::Vector<double > > const&) const at ??:?
#4 Foam::face::sweptVol(Foam::UList<Foam::Vector<doub le> > const&, Foam::UList<Foam::Vector<double> > const&) const at ??:?
#5 Foam:rimitiveMesh::movePoints(Foam::Field<Foam:: Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&) at ??:?
#6 Foam:olyMesh::movePoints(Foam::Field<Foam::Vecto r<double> > const&) at ??:?
#7 Foam::fvMesh::movePoints(Foam::Field<Foam::Vector< double> > const&) at ??:?
#8 Foam::dynamicMotionSolverFvMesh::update() at ??:?
#9 Foam::fluidSolidInterface::moveFluidMesh() at ??:?
#10 Foam::fluidSolidInterfaces::AitkenCouplingInterfac e::evolve() at ??:?
#11 ? at ??:?
#12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#13 ? at ??:?
Floating point exception

[bigphil]

Glad to hear you got it installed. How did you solve the installation problem?

If you are using the Aitken's coupling method, you can find the overall FSI coupling procedure in the evolve function in solids4foam/src/solids4FoamModels/fluidSolidInterfaces/AitkenCouplingInterface/AitkenCouplingInterface.C. Within this function, it calls a function updateForce, defined in the base class, i.e. solids4foam/src/solids4FoamModels/fluidSolidInterfaces/fluidSolidInterface/fluidSolidInterface.C

[y99]

Originally Posted by bigphil

Glad to hear you got it installed. How did you solve the installation problem?

If you are using the Aitken's coupling method, you can find the overall FSI coupling procedure in the evolve function in solids4foam/src/solids4FoamModels/fluidSolidInterfaces/AitkenCouplingInterface/AitkenCouplingInterface.C. Within this function, it calls a function updateForce, defined in the base class, i.e. solids4foam/src/solids4FoamModels/fluidSolidInterfaces/fluidSolidInterface/fluidSolidInterface.C

Thank you very much for clearing up my confusion.
Regarding the installation issue, I used the second method you mentioned and resolved it by using the system MPI.