|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Andy
Join Date: Jul 2017
Posts: 62
Rep Power: 9 ![]() |
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: ![]() [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. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
talop
Join Date: Feb 2023
Posts: 6
Rep Power: 3 ![]() |
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! |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
change the under-relaxation factors | hotboy | FLUENT | 4 | May 28, 2015 22:25 |
Mass Transfer without Phase Change in InterPhaseChangeFoam | Parisa_Khiabani | OpenFOAM Running, Solving & CFD | 5 | May 26, 2015 09:40 |
No density change - rhoSimplecFoam | krapic | OpenFOAM | 1 | November 29, 2014 18:39 |
Help with writing a macro/UDF_file/script to change BC at regular time step intervals | Codophobia | Fluent Multiphase | 0 | April 24, 2014 10:07 |
Multicomponent fluid | Andrea | CFX | 2 | October 11, 2004 05:12 |