CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Converging Diverging Nozzle with dbnsTurbFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 4, 2019, 01:30
Default Converging Diverging Nozzle with dbnsTurbFoam
  #1
Member
 
Saleh Abuhanieh
Join Date: Nov 2017
Posts: 68
Rep Power: 5
Saleh Abuhanieh is on a distinguished road
Hi Foamers,


I used different solvers to estimate the location of the shock wave for the subject case, namely: rhoSimpleFoam, sonicFoam, and rhoCentralFoam using OpenFoam 5.0. However the estimations/results for these solvers were not so accurate comparing to the experimental results.



I am trying to solve the same case using dbnsTurboFoam (density based solver available only in foam-extend 4.0) with no success till now, the value of h is going high very rapidly and floating point exception occurs at the end.


You may find the complete case folder via the below link, you can see that I tried different BC and thermophysicalProperties.


https://1drv.ms/f/s!Amq2JoBgpHbgmRHbcHXEXrg98zaG


Any kind of support will be appreciated.
Thank you in advance.



Saleh
Saleh Abuhanieh is offline   Reply With Quote

Old   January 5, 2019, 09:58
Default
  #2
Member
 
Saleh Abuhanieh
Join Date: Nov 2017
Posts: 68
Rep Power: 5
Saleh Abuhanieh is on a distinguished road
Well .. I know this solver is not so common.. but I am sure many Foamers did it..
Saleh Abuhanieh is offline   Reply With Quote

Old   December 12, 2019, 14:57
Default
  #3
Member
 
Thomas
Join Date: Nov 2017
Posts: 37
Rep Power: 5
hand90 is on a distinguished road
Hi

Were you able to solve the issue? I am using the same solver and encounter similar problems.

Many thanks
hand90 is offline   Reply With Quote

Old   December 13, 2019, 09:54
Default
  #4
Member
 
Saleh Abuhanieh
Join Date: Nov 2017
Posts: 68
Rep Power: 5
Saleh Abuhanieh is on a distinguished road
Hi Thomas,


Unfortunately no.


Regards,
Saleh Abuhanieh is offline   Reply With Quote

Old   December 13, 2019, 10:26
Default
  #5
Member
 
Thomas
Join Date: Nov 2017
Posts: 37
Rep Power: 5
hand90 is on a distinguished road
I assume I have a similar issue...The geometry I am modeling is a stator blade with a fixed inlet velocity, fixed inlet temperature and a waveTransmissive outlet. Maybe someone has an idea where the issue could be found. The Mesh seems alright with a non-orthogonality of Max: 45.89 and an average: 16.36.
I first suspected the thermophysical properties but they seem alright to me:

PHP Code:
thermoType
hPsiThermo
<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>;

// name, nMoles, mol weight, CP, Hf, mu, Pr;
mixture         air 1 28.9 1005 0 1.82e-05 0.71
This is the output I get when I run it. This clearly indicates that something goes wrong as the velocity also shoots up to unrealistic numbers.

PHP Code:
Create time

Create mesh 
for time 0

Initializing the GGI interpolator between master
/shadow patchesPERIODIC_SS/periodic_sh
Reading thermophysical properties

Selecting thermodynamics package hPsiThermo
<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
Reading field rho

Reading field U

Selecting numericFlux rusanovFluxBarthJespersenLimiter
Creating turbulence model

Selecting turbulence model type laminar

Starting time loop

Courant Number mean
0.15523467 max9.02943 velocity magnitude2

 Time 
0.0001
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -5.4748556e+08 max1.0516812e+09 average322283.12
new min120600
new max804000
bounding rho
min: -0.64592875 max2.5437391 average1.1398913
new min0.10000296
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -7.4175145e+11 max1627402.2 average: -1390944.7
new min120600
new max792543.54
bounding rho
min: -1448.9384 max2550.4832 average1.2205487
new min0.10193838
new max10
bounding p
min5152.4784 max1792217.6 average160869.08
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -7.7758231e+15 max7.7388903e+10 average: -1.820092e+10
new min120600
new max804000
bounding rho
min: -5902750.1 max3308708.2 average: -10.920043
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average230618.34
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -1.0034402e+18 max1.2001047e+15 average: -1.7356107e+12
new min120600
new max804000
bounding rho
min: -1.4206979e+09 max1.2240156e+09 average: -593.21117
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average320947.77
new max1000000
    ExecutionTime 
18.41 s

Courant Number mean
22.124283 max95043370 velocity magnitude38063491

 Time 
0.0002
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -1.4400817e+21 max7.3840623e+16 average: -2.918162e+15
new min120600
new max804000
bounding rho
min: -5.145617e+08 max3.9814147e+08 average: -729.44621
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average438814.94
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -2.7560433e+27 max1.1947501e+20 average: -2.7284239e+21
new min120600
new max804000
bounding rho
min: -9.6480279e+10 max7.897857e+10 average: -24990.598
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average717589.58
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -8.9519248e+29 max4.9622006e+26 average: -1.199368e+24
new min120600
new max804000
bounding rho
min: -1.77945e+14 max8.2541671e+13 average: -23022481
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average980527.11
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -1.8965922e+34 max1.8479372e+29 average: -3.1049571e+28
new min120600
new max804000
bounding rho
min: -1.053413e+16 max1.3694968e+16 average5.7596284e+09
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average1228861.4
new max1000000
    ExecutionTime 
30.67 s

.
.
.
.


Courant Number mean2.9573684e+58 max3.4137433e+64 velocity magnitude1.1977841e+65

 Time 
0.0008
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -1.7261537e+136 max3.5040164e+132 average: -1.9902149e+130
new min120600
new max804000
bounding rho
min: -1.6043521e+65 max9.3030893e+64 average: -2.2327806e+58
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average1617137.5
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -5.3658521e+140 max3.0051695e+135 average: -9.0106518e+134
new min120600
new max804000
bounding rho
min: -4.8938712e+77 max3.2369183e+77 average: -1.6393966e+71
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average1790213
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -4.5691667e+146 max7.9388043e+139 average: -4.5615589e+140
new min120600
new max804000
bounding rho
min: -2.4609042e+70 max1.6069814e+70 average7.3210329e+61
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average1786632.7
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -1.8277029e+150 max4.6411647e+145 average: -1.8115129e+144
new min120600
new max804000
bounding rho
min: -7.9445015e+75 max7.8572014e+75 average: -6.2886868e+67
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average1786791.7
new max1000000
    ExecutionTime 
105.26 s

Courant Number mean
4.0530947e+66 max5.2670735e+74 velocity magnitude3.6185028e+74

 Time 
0.0009
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -2.28913e+151 max2.635923e+149 average: -3.2903556e+145
new min120600
new max804000
bounding rho
min: -3.6605447e+76 max5.2267317e+76 average1.5500135e+70
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average1626119.9
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -5.3106867e+152 max1.2825681e+150 average: -1.4676302e+147
new min120600
new max804000
bounding rho
min: -1.9967951e+75 max1.1146483e+75 average2.161689e+68
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average1787851.8
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -8.9872287e+157 max2.4745995e+152 average: -9.0225504e+151
new min120600
new max804000
bounding rho
min: -2.4311779e+79 max2.5835559e+79 average1.4901175e+72
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average1785011.6
new max1000000
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -4.7798214e+160 max1.3470694e+157 average: -4.8954909e+154
new min120600
new max804000
bounding rho
min: -3.5739792e+79 max3.350678e+79 average: -7.7996527e+70
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average1787867
new max1000000
    ExecutionTime 
117.78 s

Courant Number mean
1.9003026e+72 max2.4619947e+78 velocity magnitude1.2660207e+79

 Time 
0.001
diagonal
:  Solving for rhoInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUxInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUyInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoUzInitial residual 0, Final residual 0No Iterations 0
diagonal
:  Solving for rhoEInitial residual 0, Final residual 0No Iterations 0
bounding h
min: -1.7871059e+243 max9.4626262e+159 average: -1.7681706e+237
new min120600
new max804000
bounding rho
min: -9.6538203e+80 max1.2541716e+81 average2.9290033e+74
new min0.1
new max10
bounding p
min3452.3917 max2301594.5 average1645121.1
new max1000000
Floating point exception 
(core dumped
hand90 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
Fluid Dissociation in Converging Diverging Nozzle Yuqi FLUENT 2 May 29, 2018 11:53
How to determine the distance between points on Converging diverging nozzle geometry Sorwar22 FLUENT 0 December 4, 2017 20:07
Help needed with my Converging diverging nozzle Sorwar22 FLUENT 6 November 8, 2017 17:24
Converging Diverging Nozzle in OpenFOAM danishdude OpenFOAM Running, Solving & CFD 1 September 15, 2012 00:12
help - Running Converging Diverging Nozzle with a shock wave Daco FLUENT 0 September 3, 2012 08:23


All times are GMT -4. The time now is 23:07.