|
[Sponsors] |
influence of pMin in compressibleInterFoam - oscillatory pressure boundary |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 11, 2019, 09:05 |
influence of pMin in compressibleInterFoam - oscillatory pressure boundary
|
#1 |
Member
Join Date: Apr 2018
Location: UK
Posts: 78
Rep Power: 8 |
Hi there,
I am running some bubble dynamics simulations in compressibleInterFoam. I notice that pMin has a significant effect on the solution. The initial condition is atmospheric pressure. I set the pressure p to oscillate around atmospheric condition using uniformFixedValue BC. The issue is that when the pressure on the boundary goes to 0 or below 0, the solution diverges. So if p_oscillatory >= p_atm it causes problems with negative pressure/temperature or floating point exceptions! I understand that this is because of the EoS implementation in this solver, i.e. rho = p * psi which basically means that pressure can never be 0 otherwise density is also 0. So, I tried to bypass this by setting very small values of pMin and the simulation runs completely, however it seems that the values drastically affect the solution leading to unrealistic bubble shapes and in some cases even bubble breakup. I have tried the similar cases but using a simple fixedValue BC for pressure set at atmospheric and have good and stable results. I had also tried to implement fvOptions functionality but it seems it is not supported with compressibleInterFoam. I have been stuck on this issue for a while and I really need to use the oscillatory boundary, please help ! Thanks in advance, JM |
|
August 19, 2019, 15:09 |
can anyone help, please?
|
#2 |
Member
Join Date: Apr 2018
Location: UK
Posts: 78
Rep Power: 8 |
Hi, I am still having this problem, using OF v5.
Can anyone answer my question, please? |
|
August 19, 2019, 15:50 |
|
#3 |
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
Hello!
I dont use myself the solver compressibleInterFoam... Anyway, cause the name of the solver has comperssible in it, that means the pressure is principly by definition absolute not static!!! p_absolute = p_static + 1.013e5 In comperssible solver you dont define static but absolute pressure. i.e. It is not acceptable to have < 0 pressure! 0 absoulte pressure means absolutly vacuum! 0 static pressure means the pressure you define is 1.013e5! Any value between 0 and 1.013e5 means vacuum. And higher than 1.013e5 means higher than atmpspheric... That means the pressure in the p dictionary should be p + 1.013e5! The pMin (or pMax) do nothing else than to limit the pressure to the defined value by pMin value, if the calculated pressure went less than the pMin. i.e. if the calculated pressure went less than pMin, then the calculated value will be replaced by pMin! Please dont forget to replace also the internalField with: value + 1.013e5! Else the solver will suffer much to do the first calculation steps and divergence could happend... Regards Peter |
|
August 20, 2019, 06:25 |
incorrectly estimated pressure due to pMin value
|
#4 |
Member
Join Date: Apr 2018
Location: UK
Posts: 78
Rep Power: 8 |
Dear Peter,
Thank you for your reply! I am aware of all of this and why pMin needs to be specified in the dictionary, however in compressibleInterFoam one also specifies an Equation of state for the different phases. So an equation of state will be defined for air and water in this case. I do define the initial conditions based on absolute pressure, however, it seems that based on the value of pMin that I choose, a region of incorrectly estimated pressure. As mentioned in my previous post, I am using the solver for bubble dynamics problems. The bubble shape is highly influenced by the pressure field surrounding the bubble so this will result in an incorrect shape representation or artificial breakup. Looking at the pEqn for compressibleInterFoam, pMin seems to be used for the last non orthogonal correction. My mesh is cartesian so there should be no need for non orthogonal correctors so I set this to 0 in PIMPLE, however pMin still seems to be affecting the solution. I am therefore very confused by the outcome of my simulations and not sure what else I can try to solve this problem. Please help! Regards, JM |
|
August 20, 2019, 12:25 |
|
#5 |
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
Show me please the pressure dictionary?
Which values are inside? I still belive that you define the wrong pressures in the dictionary, If I understood your discribtion of the case. Why you have pMin at all? I belive you can delete it. The solver will not find it and will ignore it (I belive). I viewd the tutorials for the solver u r using. There is no pMin in any one! One more time I dont use this solver, but as a generall role it is used if it has been defined. Regards Peter |
|
August 21, 2019, 10:36 |
pMin specification
|
#6 | |
Member
Join Date: Apr 2018
Location: UK
Posts: 78
Rep Power: 8 |
Quote:
Hi Peter, Thank you for your reply. Not sure which tutorials you have looked at for compressibleInterFoam (using v5 here) but if you go on tutorials/multiphase/compressibleInterFoam/laminar/depthCharge2D and look in constant/thermophysicalProperties you will see that pMin needs to be defined. Code:
FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object thermophysicalProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // phases (water air); pMin 10000; sigma { type liquidProperties; phase water; } // ************************************************************************* // 0.orig/p_rgh reads: Code:
FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 1e5; boundaryField { farField { type uniformFixedValue; // oscillating pressure boundary uniformValue sine; uniformValueCoeffs { amplitude 95000; frequency 20000; scale 1; level 100000; // absolute pressure t0 0; } } defaultFaces { type empty; } } // ************************************************************************* // Hope this clarifies! |
||
August 23, 2019, 06:25 |
|
#7 |
Member
Join Date: Apr 2018
Location: UK
Posts: 78
Rep Power: 8 |
Hi Peter,
Any further advice with regards to my problem? Do you think my pressure specification is OK? Thanks, JM |
|
August 23, 2019, 09:25 |
|
#8 |
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
Hello and sorry for the late answer.
You are right pMin is needed. Usually pMin is to be found in fvSolution. That why I did not find it in the tutorial you mentioned. May I ask you to post the case please. Regards Peter |
|
Tags |
compressibleinterfoam, oscillatory pressure, pmin, uniformfixedvalue |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Centrifugal fan | j0hnny | CFX | 13 | October 1, 2019 14:55 |
CFD analaysis of Pelton turbine | amodpanthee | CFX | 31 | April 19, 2018 19:02 |
Centrifugal fan-reverse flow in outlet lesds to a mass in flow field | xiexing | CFX | 3 | March 29, 2017 11:00 |
Waterwheel shaped turbine inside a pipe simulation problem | mshahed91 | CFX | 3 | January 10, 2015 12:19 |
RPM in Wind Turbine | Pankaj | CFX | 9 | November 23, 2009 05:05 |