|
[Sponsors] |
Negative absolute pressure in rhoPorousMRFSimpleFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 14, 2013, 10:22 |
|
#21 |
Senior Member
Join Date: Nov 2009
Location: Michigan
Posts: 135
Rep Power: 16 |
Thank you, Niklas and maHein.
I monitored the outlet pressure using this code: pOutlet { type patchAverage; functionObjectLibs ("libsimpleFunctionObjects.so"); verbose true; patches (outlet); log true; fields ( p ); } |
|
March 5, 2013, 13:10 |
|
#22 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
then its incorrect to set totalPressure at outlet?
I have 2 totalPressure BC at inlet and outlet.at inlet it behaves somewhat good(except supersonic velocity at 3-4 begining time steps that should be subsonic in real but it modifies at later time steps and becomes subsonic) but when outlet opens pressure goes too little and velocity increases a lot and problem diverges. could anyone argue about the formula I have written at other thread also? Code:
p0_2/pow(1+(gamma-1)/2*magSqr(internalField(U))/(gamma*R*internalField(T)),3.5) |
|
March 5, 2013, 13:18 |
|
#23 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
dear Omkar I have a foolish question!
how do you plot residuals and othe variables like you did above? I want to plot residuals by iteration numbers like you did but pyFoam draws residuals( and not variables at all) by time step in x-axis. |
|
March 8, 2013, 04:56 |
|
#24 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
can anyone guide me through two questions i propounded at two posts above?
|
|
March 8, 2013, 09:56 |
|
#25 |
Senior Member
Join Date: Nov 2009
Location: Michigan
Posts: 135
Rep Power: 16 |
You can plot residuals using gnuplot as mentioned in this post:
http://www.cfd-online.com/Forums/ope...tml#post410537 I am not sure about using total pressure at both ends. When you specify total pressure at both ends, it is like fixing the loss across the system rather than letting the system decide. This is what I think, I am not sure about it though. Plus, as per the discussions above, if the flow is leaving the domain, OF will considered the input value as static pressure because the term 1.0 - pos(phip) will go to zero. |
|
March 12, 2013, 04:13 |
|
#26 |
New Member
Alessandro
Join Date: Feb 2013
Location: Italia
Posts: 15
Rep Power: 13 |
I was looking at the files posted by you and those modified by Niklas to understand the reason of the improvement. When you say "The damping for pressure is amazing" are you referring to the decrease in tolerance from 1e-6; to in 1e-9; the fvSolution?
... p { solver GAMG; tolerance 1e-9; relTol 0.1; |
|
March 12, 2013, 08:04 |
|
#27 |
Senior Member
Join Date: Nov 2009
Location: Michigan
Posts: 135
Rep Power: 16 |
No. What I meant was that the pressure values stabilized quickly.
|
|
March 14, 2013, 17:49 |
|
#28 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
how can avoid high velocities at first iterations?maybe limiting of velocity to mach velocity?i tried but without useful result.
Efforts are thanked |
|
March 15, 2013, 05:42 |
|
#29 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
has anyone done a progress on the issue of high velocities at first few iterations in total pressure model of BC?
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Negative Absolute Pressure | mike | FLUENT | 8 | December 18, 2014 04:50 |
Manifold Absolute Pressure in Engine Simulations | abminternet | OpenFOAM | 4 | July 4, 2011 12:45 |
seeking for help about a room with negative pressure | mengyue1 | FLUENT | 0 | November 26, 2009 07:40 |
negative static pressure | Selo melo | FLUENT | 5 | December 18, 2006 10:30 |
avoiding negative absolute pressure!!! | Tania | FLUENT | 8 | March 30, 2004 19:41 |