CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions > OpenFOAM CC Toolkits for Fluid-Structure Interaction

[solids4Foam] Crash Before GAMG Pressure Solver using OpenFOAM-v1812

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 29, 2020, 22:06
Default Crash Before GAMG Pressure Solver using OpenFOAM-v1812
  #1
New Member
 
Mike Tree
Join Date: Feb 2016
Location: Charlotte, NC
Posts: 23
Rep Power: 7
treem22 is on a distinguished road
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:

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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.
Attached Files
File Type: gz aortic_valve.tar.gz (8.6 KB, 3 views)
treem22 is offline   Reply With Quote

Old   October 1, 2020, 11:53
Default
  #2
New Member
 
Mike Tree
Join Date: Feb 2016
Location: Charlotte, NC
Posts: 23
Rep Power: 7
treem22 is on a distinguished road
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 is offline   Reply With Quote

Old   October 13, 2020, 09:59
Default
  #3
New Member
 
Mike Tree
Join Date: Feb 2016
Location: Charlotte, NC
Posts: 23
Rep Power: 7
treem22 is on a distinguished road
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
treem22 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
outerCorrectorResidualControl OpenFOAM v1812 vs. OpenFOAM 6 Kobi_Faust OpenFOAM Running, Solving & CFD 5 February 28, 2019 07:51
outerCorrectorResidualControl OpenFOAM v1812 vs. OpenFOAM 6 Kobi_Faust OpenFOAM Running, Solving & CFD 0 February 27, 2019 04:30
Map of the OpenFOAM Forum - Understanding where to post your questions! wyldckat OpenFOAM 9 March 30, 2017 05:19
Wrong fluctuation of pressure in transient simulation caitao OpenFOAM Running, Solving & CFD 2 March 5, 2015 21:33
Help for the small implementation in turbulence model shipman OpenFOAM Programming & Development 25 March 19, 2014 10:08


All times are GMT -4. The time now is 04:52.