CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   SU2 (https://www.cfd-online.com/Forums/su2/)
-   -   Pressure parameter in paraview (https://www.cfd-online.com/Forums/su2/133579-pressure-parameter-paraview.html)

Sachchit April 17, 2014 02:51

Pressure parameter in paraview
 
In SU2 v3.0.0 in the paraview output file how the parameter 'pressure' has been calculated? When I calculate pressure from pressure coefficient it is coming out different than the parameter 'pressure'. Is it p/(rho*U^2) or something else?

fpalacios April 18, 2014 09:28

Quote:

Originally Posted by Sachchit (Post 486637)
In SU2 v3.0.0 in the paraview output file how the parameter 'pressure' has been calculated? When I calculate pressure from pressure coefficient it is coming out different than the parameter 'pressure'. Is it p/(rho*U^2) or something else?

Please check

http://www.cfd-online.com/Forums/su2..._flow-csv.html

for more details,

Best,
Francisco

Sachchit April 19, 2014 14:25

Still some problem
 
2 Attachment(s)
Hi,

I have gone through the file but the formula seems to be the normal formula of Cp but the problem can be seen in the Images. Pressure is varying from 0.155 to 0.536 and Cp is varying from -1.23 to 1.29. Reference pressure is 101325 Pa. So ideally Cp should be -ve by that formula. By the way Cp is coming out accurate, verified by other other software suites. But I can't understand the value of pressure.
The formula of Pressure in that file is
Pressure = Gamma_Minus_One*Density*(node[iPoint]->GetSolution(nDim+1)/Density-0.5*Velocity2);
Can you please explain this formula?

Thanks in Advance.

Zef April 20, 2014 08:42

Sachchit, i don't understand your problem. The formula in your post is simply equation of state for perfect gas:
Pressure = (Gamma - 1)*Denstity*Internal_Energy.
node[iPoint]->GetSolution(nDim+1) is the last conservative variable, which equals to Density*Total_Energy, Total_Energy = Internal_Energy - 0.5*Velocity*Velocity, where Velocity is velocity magnitude (in your formula Velocity2 is squared velocity magnitude). So, node[iPoint]->GetSolution(nDim+1)/Density is total energy in specified point and (node[iPoint]->GetSolution(nDim+1)/Density-0.5*Velocity2) is the internal energy.

Sachchit April 20, 2014 09:04

Hi Zef,

Thanks for explaining formula. I am confused with Velocity2 being y-component of velocity. I have found out later that Velocity2 is square of velocity. But can you please explain the images that I have attached earlier? I am not able to understand the numbers which are coming up in Pressure. The Coefficient of Pressure is coming accurate and pressure calculated from that coefficient is also giving good results. But the problem lies in the Pressure shown by default from the SU2. I just want to know what it is?

Thanks Zef.

fpalacios April 20, 2014 19:08

Quote:

Originally Posted by Sachchit (Post 487134)
Hi Zef,

Thanks for explaining formula. I am confused with Velocity2 being y-component of velocity. I have found out later that Velocity2 is square of velocity. But can you please explain the images that I have attached earlier? I am not able to understand the numbers which are coming up in Pressure. The Coefficient of Pressure is coming accurate and pressure calculated from that coefficient is also giving good results. But the problem lies in the Pressure shown by default from the SU2. I just want to know what it is?

Thanks Zef.

I'm not sure if this is related (I need to check the config file) but, remember that SU2 can run dimensional or non-dimensional
for example if
REF_PRESSURE= 1.0
REF_TEMPERATURE= 1.0
REF_DENSITY= 1.0
then the code is running dimensional. You will find a description of the non-dimensionalization in http://su2.stanford.edu/documents/SU2_AIAA_ASM2013.pdf

Apart from that, the non-dimensionalization is implemented in
void CConfig::SetNondimensionalization(unsigned short val_nDim, unsigned short val_iZone)

in config_structure.cpp.

Best Regards,
Francisco


All times are GMT -4. The time now is 07:21.