[Starcatcher]

Hi my solver is diverging. How can I understand which discretisation i can play around with to solve that issue?. Solver is buoyantPimpleFoam. Here is the log:


Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading thermophysical properties

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport polynomial;
thermo hPolynomial;
equationOfState icoTabulated;
specie specie;
energy sensibleEnthalpy;
}

Reading field U

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

Courant Number mean: 0 max: 0
deltaT = 0.00012000005
Time = 0.00012000005

PIMPLE: iteration 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 1.8730927e-12, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 3.7138159e-11, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 1.7817682e-12, No Iterations 1
DILUPBiCGStab: Solving for h, Initial residual = 5.7577156e-11, Final residual = 5.7577156e-11, No Iterations 0
DICPCG: Solving for p_rgh, Initial residual = 0.99798202, Final residual = 0.0099149812, No Iterations 311
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.7490191e-06, global = 2.7219716e-06, cumulative = 2.7219716e-06
DICPCG: Solving for p_rgh, Initial residual = 0.34455291, Final residual = 9.0605211e-09, No Iterations 473
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.7222694e-06, global = 2.7222691e-06, cumulative = 5.4442408e-06
ExecutionTime = 10.24 s ClockTime = 10 s

Courant Number mean: 0.00046187105 max: 71215.604
deltaT = 1.6850246e-12
Time = 0.00012000005

PIMPLE: iteration 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 2.531454e-10, Final residual = 2.531454e-10, No Iterations 0
DILUPBiCGStab: Solving for Uy, Initial residual = 8.75678e-10, Final residual = 8.75678e-10, No Iterations 0
DILUPBiCGStab: Solving for Uz, Initial residual = 3.2800995e-10, Final residual = 3.2800995e-10, No Iterations 0
DILUPBiCGStab: Solving for h, Initial residual = 0.099417735, Final residual = 1.8525106e-07, No Iterations 1
[0]
[0]
[0] --> FOAM FATAL ERROR: (openfoam-2212 patch=230612)
[0] Maximum number of iterations exceeded: 100 when starting from T0:281.04671 old T:74904.669 new T:56268.39 f:1.1808757e+10 p:100000 tol:0.028104671
[0]
[0] From Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hPolynomialThermo<Foam::icoTabulated<Foam::s pecie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hPolynomialThermo<Foam ::icoTabulated<Foam::specie> >, Foam::sensibleEnthalpy>]
[0] in file ./src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.
[0]
FOAM parallel run aborting
[0]
[0] #0 Foam::error:rintStack(Foam::Ostream&) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
[0] #1 Foam::error::simpleExit(int, bool) in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
[0] #2 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so
[0] #3 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so
[0] #4 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so
[0] #5 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/buoyantPimpleFoam
[0] #6 ? in /lib/x86_64-linux-gnu/libc.so.6
[0] #7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
[0] #8 ? in /usr/lib/openfoam/openfoam2212/platforms/linux64GccDPInt32Opt/bin/buoyantPimpleFoam
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.

[talop]

It can be anything for this unknown case. But the simulation crashes at predicting the temperature field so you can try reducing the relaxation factor for the energy equation.

https://www.openfoam.com/documentati...8H_source.html

Cheers!