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
);
}

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?
thanks.

dear Omkar I have a foolish question!:D
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.

can anyone guide me through two questions i propounded at two posts above?

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.

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;

No. What I meant was that the pressure values stabilized quickly.

how can avoid high velocities at first iterations?maybe limiting of velocity to mach velocity?i tried but without useful result.
Efforts are thanked

has anyone done a progress on the issue of high velocities at first few iterations in total pressure model of BC?