CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Discrepancy in Pressure Driven Flow Behaviour in 211 vs 221

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   September 24, 2013, 12:15
Default Discrepancy in Pressure Driven Flow Behaviour in 211 vs 221, solver: rhoCentralFoam
  #1
cdm
Member
 
Join Date: May 2013
Location: Canada
Posts: 32
Rep Power: 4
cdm is on a distinguished road
Not sure if this should be posted under Bugs, but I'm hoping someone can verify/clarify an issue I'm having when transitioning from OF 211 to 221.

I spent some time investigating k-Omega SST wall function behaviour over a flat plate using the rhoCentralFoam solver and pressure driven flow. I wanted to ensure that the k-Omega SST wall functions would follow closely the standard log-law for a mesh with yPlusRAS of roughly 75 in the area of interest. However, in the process of investigating this, I discovered a more fundamental issue that hopefully someone might be able to clarify.

The case was set up as shown in the image below. The inlet uses totalPressure and totalTemperature boundary conditions, and the outlet is fixed pressure with all else being zeroGradient. The case was set up to approximate Ma = 0.8 subsonic, compressible flow.
flatplate_setup.jpg

The mesh was created in ICEM and converted for use in OF using fluent3DMeshToFoam. Pressure driven flow was chosen to match as closely as possible the simulation setup in ANSYS, as we sought to compare results of both utilities.

A number of issues came up, not least of which being that the wall shear stress calculated using the wallShearStress utility was significantly impacted by the chosen fvSchemes (mainly, the Interpolation Schemes). However, more importantly to this post, when plotting the inlet pressure/velocity distributions there is a significant difference between the versions tested of OpenFOAM, all else being equal between the cases. In version 2.2.1 inlet pressure and velocity was significantly impacted by changes to the Gamma differencing value specified in the interpolationSchemes of the fvSolution dictionary.

The plots below illustrate the differences in each of the tested cases, for pressure and velocity plotted over the inlet. Summarised, in version 2.1.1 the sensitivity to changes in the gamma differencing (or vanLeer) was negligible, while for 2.2.1 the sensitivity was severe.

For reference, Fluent yielded an inlet P = 101683 Pa, U = 273 m/s.

OpenFOAM v2.1.1:
U = 274 m/s (Gamma 1.0)
U = 274 m/s (Gamma 0.5)
U = 274 m/s (vanLeer)
P = 101665 Pa (Gamma 1.0)
P = 101589 Pa (Gamma 0.5)
P = 101583 Pa (vanLeer)

OpenFOAM v2.2.1:
U = 316 m/s (Gamma 1.0)
U = 292 m/s (Gamma 0.5)
U = 292 m/s (vanLeer)
P = 102845 Pa (Gamma 1.0)
P = 101771 Pa (Gamma 0.5)
P = 101737 Pa (vanLeer)

Inlet Pressure:
inlet_P.jpg

Inlet Velocity:
inlet_U.jpg

As is obvious from the plots, the pressure and velocity change at the inlet is negligible in 2.1.1 but is quite drastic in 2.2.1. The totalPressure boundary condition has not changed - the only change I can find is in the rhoCentralFoam energy equation: two previously "corrective" terms have been removed (appropriately). So why is there such a sensitivity suddenly to the interpolationSchemes? v2.1.1 matches very well with Fluent, but v.2.2.1 is all over the place.

Any ideas?

Last edited by cdm; September 26, 2013 at 19:23. Reason: Corrected title to add solver used; fixed typo
cdm is offline   Reply With Quote

Old   October 28, 2013, 19:03
Default
  #2
cdm
Member
 
Join Date: May 2013
Location: Canada
Posts: 32
Rep Power: 4
cdm is on a distinguished road
I hate to dredge up old posts, but I was hoping someone may have some insight into the above issue.
cdm is offline   Reply With Quote

Old   October 31, 2013, 16:40
Default
  #3
Member
 
Join Date: Jun 2012
Posts: 67
Rep Power: 5
maHein is on a distinguished road
Hello cdm,

have you tried to re-insert the two terms you mentioned into the solver? So one could see, if these terms are the source for the differences?

I used rhoCentralFoam for validation as well. I always got better results with the vanLeer scheme compared to results with the gamma scheme. I think, the vanLeer scheme is also used in the paper, where rhoCentralFoam was presented.

Regards,

maHein
maHein is offline   Reply With Quote

Old   November 1, 2013, 15:06
Default
  #4
cdm
Member
 
Join Date: May 2013
Location: Canada
Posts: 32
Rep Power: 4
cdm is on a distinguished road
Hi maHein,

Thanks for your reply. While I haven't added the two terms back into the energy equation, they should not be there anyway, so if their removal is what is causing this discrepancy, there is a worse underlying issue in rhoCentralFoam. That said, I haven't gotten into the world of compiling my own solvers nor have I had time to learn how to do so.

Which version of OF did you use for your validations with rhoCentralFoam? Did you encounter any discrepancies?
cdm is offline   Reply With Quote

Old   December 10, 2013, 17:39
Default
  #5
cdm
Member
 
Join Date: May 2013
Location: Canada
Posts: 32
Rep Power: 4
cdm is on a distinguished road
Just an update:

Upon further investigation, I found that the discrepancy is caused by the modification (and correction) of the energy equation in version 22x. See here and here for info on the modified term.
cdm is offline   Reply With Quote

Old   December 17, 2013, 08:48
Default
  #6
New Member
 
caduqued
Join Date: Apr 2009
Location: UK
Posts: 15
Rep Power: 8
caduqued is on a distinguished road
Quote:
Originally Posted by cdm View Post
Just an update:

Upon further investigation, I found that the discrepancy is caused by the modification (and correction) of the energy equation in version 22x. See here and here for info on the modified term.
Dear maHein,

Have you been able to check if this change corrects your initial case? Should we now assume that the latest version of OF correctly validates, for instance, in your case?

Regards,

caduqued
caduqued is offline   Reply With Quote

Reply

Tags
rhocentralfoam, totalpressure

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
internal flow BCs: pressure driven versus velocity driven mihaipruna OpenFOAM Running, Solving & CFD 22 March 6, 2014 11:06
Pressure driven flow in vertical pipe Hale OpenFOAM Pre-Processing 6 September 13, 2013 08:16
reference pressure and compressible flow bingo10 CFX 0 September 11, 2013 07:32
pressure BC in buoynacy driven flow Sasha FLUENT 3 October 11, 2006 10:08
Hydrostatic pressure in 2-phase flow modeling (CFX4.2) HB &DS CFX 0 January 9, 2000 14:19


All times are GMT -4. The time now is 16:26.