|
[Sponsors] |
Can not use very high pressure with sonicFoam and rhoCentralFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 20, 2016, 09:56 |
Can not use very high pressure with sonicFoam and rhoCentralFoam
|
#1 |
Senior Member
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 23 |
I am trying to simulate shock wave propagation using sonicFoam and rhoCentralFoam. The idea is to patch a very high pressure for a small spherical region and let the shock wave propagate through the cubical domain.
The simulation goes on nicely without any errors or divergence issues. So, no other setup issues. The problem is with patching high pressure using setField. I would like to patch as high as 1e11 Pa, but no matter what I do the value is always limited to 105686784 Pa. For very few trials I managed to patch the high pressure values but in the very next time-step the highest values is switched to 105686784 Pa. Is sonicFoam or rhoCentralFoam limiting pressure values some how? I understand that the equation of state for Ideal gas plays role in updating vector fields inside the PISO loop. But I dont know for sure if that is the cause for this discrepancy. I am very confused and dont know what to do about it. Any hint, suggestion or explanation? |
|
September 21, 2016, 04:44 |
|
#2 |
Senior Member
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 23 |
If I try to patch pressure of order 10^8 (e.g. 9E+8), everything works fine. But any higher order, it simply does not work.
|
|
September 22, 2016, 03:44 |
Solution!!
|
#3 |
Senior Member
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 23 |
It turns out that its just matter of precision. I mean,
Code:
writePrecision 8;//12;// timePrecision 8;//12;// I would still appreciate some comments from other users. |
|
September 22, 2016, 03:49 |
|
#4 |
Senior Member
Join Date: Oct 2013
Posts: 397
Rep Power: 19 |
This is strange, if anything I would have expected a limit in one of the thermophysical models. Using a lower precision should also not cause this as far as I can tell.
|
|
October 22, 2024, 14:05 |
|
#5 |
New Member
Yas
Join Date: Oct 2024
Posts: 1
Rep Power: 0 |
I am trying to do the same thing but I keep getting negative temperatures, and can't get it to solve past 1MPa, would you kindly be able to send me one of your simulations that you've done at a high pressure so I can see where I'm going wrong
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Unable to validate results with sonicFoam and rhoCentralFoam | lordvon | OpenFOAM Running, Solving & CFD | 8 | December 5, 2016 15:25 |
Which equations are solved by rhoCentralFoam, rhoSimplecFoam and sonicFoam? | gamma_user | OpenFOAM Verification & Validation | 1 | April 19, 2015 13:13 |
sonicFoam vs rhoCentralFoam forwardStep tutorial | hk318i | OpenFOAM Running, Solving & CFD | 2 | March 12, 2015 12:15 |
rhoCentralFoam / sonicFoam | VSass | OpenFOAM Running, Solving & CFD | 1 | June 26, 2013 13:39 |
is sonicFoam better in convergence than rhoPimplefoam or rhoCentralFoam? | immortality | OpenFOAM Running, Solving & CFD | 1 | February 11, 2013 23:26 |