[treem22]

After observing the forum for a few days I decided I was going to try and implement solids4Foam with OpenFOAM-v1812. I am wary of foam-extend-4.1 due to the reported pimpleFoam and pisoFoam instabilities (https://sourceforge.net/p/foam-extend/tickets/39/). I was able to install OpenFOAM-v1812 and compile the solids4Foam release. I can successfully run the linearGeometryElasticBeam.openfoam tutorial case.

So, I'm branching out from there. I am attempting to run flow (forward flow only for now) through an aortic valve geometry using the two-way FSI features that solids4Foam affords. I'm copying a lot from the linearGeometryElasticBeam.openfoam tutorial case, but obviously providing a different geometry and boundary conditions. It seems that no matter what I do, the simulation crashes in the first FSI iteration just after solving for the fluid velocity. I'm assuming the fluid pressure is what's to be solved next, but the error messages aren't that helpful:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1812                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v1812 OPENFOAM=1812
Arch   : "LSB;label=32;scalar=64"
Exec   : solids4Foam
Date   : Sep 29 2020
Time   : 17:36:16
Host   : corvidpost5.corvidtec.com
PID    : 49239
I/O    : uncollated
Case   : /beegfs/users/mtree/OpenFOAM/mtree-v1812/run/AorticValve
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing 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 fluidSolidInteraction

Selecting fluidSolidInterface method fixedRelaxation

Selecting fluidModel icoFluid
Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: velocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Selecting patchDistMethod meshWave
g field not found in constant directory: initialising to zero
Selecting solidModel unsNonLinearGeometryTotalLagrangian
Selecting dynamicFvMesh staticFvMesh
Creating solidTraction boundary condition
    limiter coefficient: 1
Creating fixedDisplacement boundary condition
    under-relaxation method: fixed
Creating the mechanicalModel
Selecting mechanical law neoHookeanElastic
Selecting interfaceToInterfaceMapping directMap
Time = 0.0005

Setting traction on solid interfaces
Mapping face values using directMap
Total force on fluid interface 0: (-5.84492e-05 -0.00332125 -0.229065)
Total force on solid interface 0: (5.84492e-05 0.00332125 0.229065)

Evolving solid solver
Corr 0, relative residual = 0
FDICPCG: Solving for D, Initial residual = 0, Final residual = (0 0 0), No outer iterations = 0
 Max relative residual = 0, Relative residual = 0, enforceLinear = false
Mapping point values using directMap
Mapping point values using directMap
FSI relative residual1 norm for interface 0: 0
FSI residual2 norm for interface 0: 0

Time = 0.0005, iteration: 1
Current fsi under-relaxation factor: 0.4
Maximal accumulated displacement of interface 0: 0
GAMG:  Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG:  Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG:  Solving for cellMotionUz, Initial residual = 0, Final residual = 0, No Iterations 1
Evolving fluid model: icoFluid
volume continuity errors : sum local = 0, global = 0
Courant Number mean: 0.00123411 max: 0.154609 velocity magnitude: 0.0999959

PISO: Operating solver in PISO mode

DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 7.94578e-07, No Iterations 30
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 9.55747e-07, No Iterations 31
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 9.46624e-07, No Iterations 31
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib64/libc.so.6
#3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
#8  Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
#9  Foam::fluidModels::icoFluid::evolve() at ??:?
#10  Foam::fluidSolidInterfaces::fixedRelaxationCouplingInterface::evolve() at ??:?
#11  ? at ??:?
#12  __libc_start_main in /lib64/libc.so.6
#13  ? at ??:?

If anyone has any ideas, it will be much appreciated. My case directory (without a mesh) is attached.

[treem22]

I've changed changed pressure solvers without success. I've even backed away from FSI and tried to just solve with potentialFoam without success. There is something more fundamentally wrong with my case, but I can't seem to see it. Anyone see this glaringly obvious mistake I'm making?

I'm going to remove the leaflets from the aortic valve to see if I can just get flow through an anatomically-shaped pipe. If this doesn't work, it's got to be some compilation issue, right?

[treem22]

I can confirm this was a compilation issue. I cleaned my install and re-tried compiling and there are no errors. I'm still working on the aortic valve FSI sim, but the solvers are working correctly. All that's left is my ineptitude to operate them