|
[Sponsors] |
August 17, 2021, 01:27 |
Supercritical flow
|
#1 |
New Member
Piedao-Chuan
Join Date: Jul 2021
Posts: 3
Rep Power: 4 |
This is a supercritical flow (O2 、H2), but it always diverges in the calculation process:疯狂的::疯狂的:. I hope you can give me some suggestions. Thank you very much.
: Courant Number mean: 2.51125 max: 18.699 Time = 0.0003855 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 min/max(U)) = 0, 149.098 update D with binary diffusion model... species transprotation equation DILUPBiCG: Solving for O2, Initial residual = 0.00116254, Final residual = 2.02936e-06, No Iterations 3 DILUPBiCG: Solving for H2, Initial residual = 0.000823762, Final residual = 5.40367e-07, No Iterations 3 Calculating inert species N2 Solving energy equation DILUPBiCG: Solving for h, Initial residual = 0.00355332, Final residual = 4.344e-09, No Iterations 6 Starting correcting themo #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 in "/lib/x86_64-linux-gnu/libm.so.6" #4 pow in "/lib/x86_64-linux-gnu/libm.so.6" #5 Foam::PengRobinson<Foam::realGasSpecie>::Z(double, double) const at ??:? #6 Foam::realGasJanafThermo<Foam::PengRobinson<Foam:: realGasSpecie> >::cp(double, double) const at ??:? #7 Foam::realGasHePsiThermo<Foam:siReactionThermo, Foam::SpecieMixture<Foam::realGasReactingMixture<F oam::ChungHighPressureTransport<Foam::species::rea lGasThermo<Foam::realGasJanafThermo<Foam::PengRobi nson<Foam::realGasSpecie> >, Foam::sensibleEnthalpy> > > > >::calculate() at ??:? #8 Foam::realGasHePsiThermo<Foam:siReactionThermo, Foam::SpecieMixture<Foam::realGasReactingMixture<F oam::ChungHighPressureTransport<Foam::species::rea lGasThermo<Foam::realGasJanafThermo<Foam::PengRobi nson<Foam::realGasSpecie> >, Foam::sensibleEnthalpy> > > > >::correct() at ??:? #9 at ??:? #10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #11 at ??:? |
|
August 17, 2021, 12:28 |
|
#2 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
> Courant Number mean: 2.51125 max: 18.699
You should not have Courant numbers beyond 1. I recommend even lower values. Use the automatic time step control which ensures that.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
August 18, 2021, 05:34 |
|
#3 | |
New Member
Piedao-Chuan
Join Date: Jul 2021
Posts: 3
Rep Power: 4 |
Quote:
Thank you very much. I have made changes based on your suggestion. Unfortunately, the result is still divergent. Courant Number mean: 0.0499972 max: 0.11485 Time = 0.00164182 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 min/max(U)) = 0, 43.6587 update D with binary diffusion model... species transprotation equation DILUPBiCG: Solving for O2, Initial residual = 4.01292e-06, Final residual = 6.74252e-19, No Iterations 3 DILUPBiCG: Solving for H2, Initial residual = 3.93064e-06, Final residual = 2.12152e-18, No Iterations 3 Calculating inert species N2 Solving energy equation DILUPBiCG: Solving for h, Initial residual = 2.2811e-05, Final residual = 1.54394e-10, No Iterations 3 Starting correcting themo #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 in "/lib/x86_64-linux-gnu/libm.so.6" #4 pow in "/lib/x86_64-linux-gnu/libm.so.6" #5 Foam::PengRobinson<Foam::realGasSpecie>::Z(double, double) const at ??:? #6 Foam::realGasJanafThermo<Foam::PengRobinson<Foam:: realGasSpecie> >::cp(double, double) const at ??:? #7 Foam::realGasHePsiThermo<Foam:siReactionThermo, Foam::SpecieMixture<Foam::realGasReactingMixture<F oam::ChungHighPressureTransport<Foam::species::rea lGasThermo<Foam::realGasJanafThermo<Foam::PengRobi nson<Foam::realGasSpecie> >, Foam::sensibleEnthalpy> > > > >::calculate() at ??:? #8 Foam::realGasHePsiThermo<Foam:siReactionThermo, Foam::SpecieMixture<Foam::realGasReactingMixture<F oam::ChungHighPressureTransport<Foam::species::rea lGasThermo<Foam::realGasJanafThermo<Foam::PengRobi nson<Foam::realGasSpecie> >, Foam::sensibleEnthalpy> > > > >::correct() at ??:? #9 at ??:? #10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #11 at ??:? 浮点数例外 (核心已转储) |
||
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Will the results of steady state solver and transient solver be same? | carye | OpenFOAM Running, Solving & CFD | 9 | December 28, 2019 05:21 |
mass flow in is not equal to mass flow out | saii | CFX | 12 | March 19, 2018 05:21 |
transient, impregnating flow problem | fgommer | FLUENT | 0 | February 29, 2012 16:10 |
Flow meter Design | CD adapco Group Marketing | Siemens | 3 | June 21, 2011 08:33 |
Modelling transonic 2D flow over supercritical air | Teddy | FLUENT | 0 | November 20, 2007 17:03 |