CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Bugs

Supercritical flow

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 17, 2021, 02:27
Default Supercritical flow
  #1
New Member
 
Piedao-Chuan
Join Date: Jul 2021
Posts: 3
Rep Power: 3
chang_zc is on a distinguished road
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 ??:?
chang_zc is offline   Reply With Quote

Old   August 17, 2021, 13:28
Default
  #2
Senior Member
 
piu58's Avatar
 
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 727
Rep Power: 13
piu58 is on a distinguished road
> 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)
piu58 is offline   Reply With Quote

Old   August 18, 2021, 06:34
Default
  #3
New Member
 
Piedao-Chuan
Join Date: Jul 2021
Posts: 3
Rep Power: 3
chang_zc is on a distinguished road
Quote:
Originally Posted by piu58 View Post
> 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.

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 ??:?
浮点数例外 (核心已转储)
chang_zc 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
Will the results of steady state solver and transient solver be same? carye OpenFOAM Running, Solving & CFD 9 December 28, 2019 06:21
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 06:21
transient, impregnating flow problem fgommer FLUENT 0 February 29, 2012 17:10
Flow meter Design CD adapco Group Marketing Siemens 3 June 21, 2011 09:33
Modelling transonic 2D flow over supercritical air Teddy FLUENT 0 November 20, 2007 18:03


All times are GMT -4. The time now is 13:53.