Converging Diverging Nozzle with dbnsTurbFoam

Saleh Abuhanieh · January 4, 2019, 02:30

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 · January 5, 2019, 10:58

Well .. I know this solver is not so common.. but I am sure many Foamers did it..

hand90 · December 12, 2019, 15:57

Hi

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

Many thanks

Saleh Abuhanieh · December 13, 2019, 10:54

Hi Thomas,

Unfortunately no.

Regards,

hand90 · December 13, 2019, 11:26

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 patches: PERIODIC_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 max: 9.02943 velocity magnitude: 2

 Time = 0.0001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -5.4748556e+08 max: 1.0516812e+09 average: 322283.12
new min: 120600
new max: 804000
bounding rho, min: -0.64592875 max: 2.5437391 average: 1.1398913
new min: 0.10000296
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -7.4175145e+11 max: 1627402.2 average: -1390944.7
new min: 120600
new max: 792543.54
bounding rho, min: -1448.9384 max: 2550.4832 average: 1.2205487
new min: 0.10193838
new max: 10
bounding p, min: 5152.4784 max: 1792217.6 average: 160869.08
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -7.7758231e+15 max: 7.7388903e+10 average: -1.820092e+10
new min: 120600
new max: 804000
bounding rho, min: -5902750.1 max: 3308708.2 average: -10.920043
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 230618.34
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -1.0034402e+18 max: 1.2001047e+15 average: -1.7356107e+12
new min: 120600
new max: 804000
bounding rho, min: -1.4206979e+09 max: 1.2240156e+09 average: -593.21117
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 320947.77
new max: 1000000
    ExecutionTime = 18.41 s

Courant Number mean: 22.124283 max: 95043370 velocity magnitude: 38063491

 Time = 0.0002
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -1.4400817e+21 max: 7.3840623e+16 average: -2.918162e+15
new min: 120600
new max: 804000
bounding rho, min: -5.145617e+08 max: 3.9814147e+08 average: -729.44621
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 438814.94
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -2.7560433e+27 max: 1.1947501e+20 average: -2.7284239e+21
new min: 120600
new max: 804000
bounding rho, min: -9.6480279e+10 max: 7.897857e+10 average: -24990.598
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 717589.58
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -8.9519248e+29 max: 4.9622006e+26 average: -1.199368e+24
new min: 120600
new max: 804000
bounding rho, min: -1.77945e+14 max: 8.2541671e+13 average: -23022481
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 980527.11
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -1.8965922e+34 max: 1.8479372e+29 average: -3.1049571e+28
new min: 120600
new max: 804000
bounding rho, min: -1.053413e+16 max: 1.3694968e+16 average: 5.7596284e+09
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 1228861.4
new max: 1000000
    ExecutionTime = 30.67 s

.
.
.
.


Courant Number mean: 2.9573684e+58 max: 3.4137433e+64 velocity magnitude: 1.1977841e+65

 Time = 0.0008
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -1.7261537e+136 max: 3.5040164e+132 average: -1.9902149e+130
new min: 120600
new max: 804000
bounding rho, min: -1.6043521e+65 max: 9.3030893e+64 average: -2.2327806e+58
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 1617137.5
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -5.3658521e+140 max: 3.0051695e+135 average: -9.0106518e+134
new min: 120600
new max: 804000
bounding rho, min: -4.8938712e+77 max: 3.2369183e+77 average: -1.6393966e+71
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 1790213
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -4.5691667e+146 max: 7.9388043e+139 average: -4.5615589e+140
new min: 120600
new max: 804000
bounding rho, min: -2.4609042e+70 max: 1.6069814e+70 average: 7.3210329e+61
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 1786632.7
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -1.8277029e+150 max: 4.6411647e+145 average: -1.8115129e+144
new min: 120600
new max: 804000
bounding rho, min: -7.9445015e+75 max: 7.8572014e+75 average: -6.2886868e+67
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 1786791.7
new max: 1000000
    ExecutionTime = 105.26 s

Courant Number mean: 4.0530947e+66 max: 5.2670735e+74 velocity magnitude: 3.6185028e+74

 Time = 0.0009
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -2.28913e+151 max: 2.635923e+149 average: -3.2903556e+145
new min: 120600
new max: 804000
bounding rho, min: -3.6605447e+76 max: 5.2267317e+76 average: 1.5500135e+70
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 1626119.9
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -5.3106867e+152 max: 1.2825681e+150 average: -1.4676302e+147
new min: 120600
new max: 804000
bounding rho, min: -1.9967951e+75 max: 1.1146483e+75 average: 2.161689e+68
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 1787851.8
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -8.9872287e+157 max: 2.4745995e+152 average: -9.0225504e+151
new min: 120600
new max: 804000
bounding rho, min: -2.4311779e+79 max: 2.5835559e+79 average: 1.4901175e+72
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 1785011.6
new max: 1000000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -4.7798214e+160 max: 1.3470694e+157 average: -4.8954909e+154
new min: 120600
new max: 804000
bounding rho, min: -3.5739792e+79 max: 3.350678e+79 average: -7.7996527e+70
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 1787867
new max: 1000000
    ExecutionTime = 117.78 s

Courant Number mean: 1.9003026e+72 max: 2.4619947e+78 velocity magnitude: 1.2660207e+79

 Time = 0.001
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
bounding h, min: -1.7871059e+243 max: 9.4626262e+159 average: -1.7681706e+237
new min: 120600
new max: 804000
bounding rho, min: -9.6538203e+80 max: 1.2541716e+81 average: 2.9290033e+74
new min: 0.1
new max: 10
bounding p, min: 3452.3917 max: 2301594.5 average: 1645121.1
new max: 1000000
Floating point exception (core dumped)

January 4, 2019, 02:30	Converging Diverging Nozzle with dbnsTurbFoam	#1
Saleh Abuhanieh Member Saleh Abuhanieh Join Date: Nov 2017 Posts: 83 Rep Power: 9	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

December 13, 2019, 11:26		#5
hand90 Member Thomas Join Date: Nov 2017 Posts: 37 Rep Power: 9	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 patches: PERIODIC_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 max: 9.02943 velocity magnitude: 2 Time = 0.0001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -5.4748556e+08 max: 1.0516812e+09 average: 322283.12 new min: 120600 new max: 804000 bounding rho, min: -0.64592875 max: 2.5437391 average: 1.1398913 new min: 0.10000296 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -7.4175145e+11 max: 1627402.2 average: -1390944.7 new min: 120600 new max: 792543.54 bounding rho, min: -1448.9384 max: 2550.4832 average: 1.2205487 new min: 0.10193838 new max: 10 bounding p, min: 5152.4784 max: 1792217.6 average: 160869.08 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -7.7758231e+15 max: 7.7388903e+10 average: -1.820092e+10 new min: 120600 new max: 804000 bounding rho, min: -5902750.1 max: 3308708.2 average: -10.920043 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 230618.34 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -1.0034402e+18 max: 1.2001047e+15 average: -1.7356107e+12 new min: 120600 new max: 804000 bounding rho, min: -1.4206979e+09 max: 1.2240156e+09 average: -593.21117 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 320947.77 new max: 1000000 ExecutionTime = 18.41 s Courant Number mean: 22.124283 max: 95043370 velocity magnitude: 38063491 Time = 0.0002 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -1.4400817e+21 max: 7.3840623e+16 average: -2.918162e+15 new min: 120600 new max: 804000 bounding rho, min: -5.145617e+08 max: 3.9814147e+08 average: -729.44621 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 438814.94 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -2.7560433e+27 max: 1.1947501e+20 average: -2.7284239e+21 new min: 120600 new max: 804000 bounding rho, min: -9.6480279e+10 max: 7.897857e+10 average: -24990.598 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 717589.58 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -8.9519248e+29 max: 4.9622006e+26 average: -1.199368e+24 new min: 120600 new max: 804000 bounding rho, min: -1.77945e+14 max: 8.2541671e+13 average: -23022481 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 980527.11 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -1.8965922e+34 max: 1.8479372e+29 average: -3.1049571e+28 new min: 120600 new max: 804000 bounding rho, min: -1.053413e+16 max: 1.3694968e+16 average: 5.7596284e+09 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 1228861.4 new max: 1000000 ExecutionTime = 30.67 s . . . . Courant Number mean: 2.9573684e+58 max: 3.4137433e+64 velocity magnitude: 1.1977841e+65 Time = 0.0008 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -1.7261537e+136 max: 3.5040164e+132 average: -1.9902149e+130 new min: 120600 new max: 804000 bounding rho, min: -1.6043521e+65 max: 9.3030893e+64 average: -2.2327806e+58 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 1617137.5 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -5.3658521e+140 max: 3.0051695e+135 average: -9.0106518e+134 new min: 120600 new max: 804000 bounding rho, min: -4.8938712e+77 max: 3.2369183e+77 average: -1.6393966e+71 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 1790213 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -4.5691667e+146 max: 7.9388043e+139 average: -4.5615589e+140 new min: 120600 new max: 804000 bounding rho, min: -2.4609042e+70 max: 1.6069814e+70 average: 7.3210329e+61 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 1786632.7 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -1.8277029e+150 max: 4.6411647e+145 average: -1.8115129e+144 new min: 120600 new max: 804000 bounding rho, min: -7.9445015e+75 max: 7.8572014e+75 average: -6.2886868e+67 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 1786791.7 new max: 1000000 ExecutionTime = 105.26 s Courant Number mean: 4.0530947e+66 max: 5.2670735e+74 velocity magnitude: 3.6185028e+74 Time = 0.0009 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -2.28913e+151 max: 2.635923e+149 average: -3.2903556e+145 new min: 120600 new max: 804000 bounding rho, min: -3.6605447e+76 max: 5.2267317e+76 average: 1.5500135e+70 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 1626119.9 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -5.3106867e+152 max: 1.2825681e+150 average: -1.4676302e+147 new min: 120600 new max: 804000 bounding rho, min: -1.9967951e+75 max: 1.1146483e+75 average: 2.161689e+68 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 1787851.8 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -8.9872287e+157 max: 2.4745995e+152 average: -9.0225504e+151 new min: 120600 new max: 804000 bounding rho, min: -2.4311779e+79 max: 2.5835559e+79 average: 1.4901175e+72 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 1785011.6 new max: 1000000 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -4.7798214e+160 max: 1.3470694e+157 average: -4.8954909e+154 new min: 120600 new max: 804000 bounding rho, min: -3.5739792e+79 max: 3.350678e+79 average: -7.7996527e+70 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 1787867 new max: 1000000 ExecutionTime = 117.78 s Courant Number mean: 1.9003026e+72 max: 2.4619947e+78 velocity magnitude: 1.2660207e+79 Time = 0.001 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 bounding h, min: -1.7871059e+243 max: 9.4626262e+159 average: -1.7681706e+237 new min: 120600 new max: 804000 bounding rho, min: -9.6538203e+80 max: 1.2541716e+81 average: 2.9290033e+74 new min: 0.1 new max: 10 bounding p, min: 3452.3917 max: 2301594.5 average: 1645121.1 new max: 1000000 Floating point exception (core dumped)

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Fluid Dissociation in Converging Diverging Nozzle	Yuqi	FLUENT	2	May 29, 2018 12:53
How to determine the distance between points on Converging diverging nozzle geometry	Sorwar22	FLUENT	0	December 4, 2017 21:07
Help needed with my Converging diverging nozzle	Sorwar22	FLUENT	6	November 8, 2017 18:24
Converging Diverging Nozzle in OpenFOAM	danishdude	OpenFOAM Running, Solving & CFD	1	September 15, 2012 01:12
help - Running Converging Diverging Nozzle with a shock wave	Daco	FLUENT	0	September 3, 2012 09:23

January 5, 2019, 10:58		#2
Saleh Abuhanieh Member Saleh Abuhanieh Join Date: Nov 2017 Posts: 83 Rep Power: 9	Well .. I know this solver is not so common.. but I am sure many Foamers did it..

December 12, 2019, 15:57		#3
hand90 Member Thomas Join Date: Nov 2017 Posts: 37 Rep Power: 9	Hi Were you able to solve the issue? I am using the same solver and encounter similar problems. Many thanks

December 13, 2019, 10:54		#4
Saleh Abuhanieh Member Saleh Abuhanieh Join Date: Nov 2017 Posts: 83 Rep Power: 9	Hi Thomas, Unfortunately no. Regards,