CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > FLUENT

Compressible flow over two-dimensional bump geometry

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 15, 2013, 05:18
Question Compressible flow over two-dimensional bump geometry
  #1
Senior Member
 
Crank-Shaft's Avatar
 
Ovi
Join Date: Oct 2012
Location: Sydney, Australia
Posts: 141
Rep Power: 4
Crank-Shaft is on a distinguished road
Hello Forum regulars and guests,

I have recently started working with a 2D geometry consisting of a bump on a wall, a defined inlet and outlet and symmetry conditions on the rest of the exterior boundaries. This is described in detailed within the images attached. This is a benchmark study from NASA and I would like to use it as a reference for code verification. (For Details - http://turbmodels.larc.nasa.gov/bump.html)

I am looking for suggestions to implement the boundary conditions required, to match the reference results. I am trying to obtain steady state solutions using an ideal-gas model, with gravitational acceleration in -y axis. The current boundary conditions including the pressure outlet and inlet have an operating pressure of 101 325 Pa with a gauge pressure set to (101 325 x 0.02828) Pa. This pressure gradient results in a velocity of 69.7 m/s which satisfies the M=0.2 condition. I currently have the energy equation turned off, since I am working through the preliminary stages of the simulation setup. I have received several warnings and recommendations when I clicked Check Setup before commencing with the solver iterations, as shown in the image.

All the convergence criteria for the residuals were defined as 1e-6 however, the continuity simply remains at 1e-2 despite the large number of iterations performed. I tried changing the momentum relaxation factor from 0.8 to 0.5 however, this only had a minimal effect.

Please help me implement the boundary conditions outlined in the first two images and also provide some guidance to improve the absolute continuity residuals. I look forward to your comments.
Attached Images
File Type: jpg bumpBCpic.jpg (47.6 KB, 30 views)
File Type: jpg bumpBCpic2.jpg (39.6 KB, 26 views)
File Type: png PressureWarnings.PNG (11.2 KB, 17 views)
__________________
--
Mechanical Engineering
Sydney, Australia


Crank-Shaft is offline   Reply With Quote

Old   January 16, 2013, 04:46
Default
  #2
Senior Member
 
Crank-Shaft's Avatar
 
Ovi
Join Date: Oct 2012
Location: Sydney, Australia
Posts: 141
Rep Power: 4
Crank-Shaft is on a distinguished road
Here are some images showing the Pressure inlet settings and also the warnings received when using the Check Case feature prior to solving.

Please help me understand the messages and share some suggestions to correctly model this.
Attached Images
File Type: png Inlet.PNG (24.8 KB, 13 views)
File Type: png PressureWarnings.PNG (11.2 KB, 13 views)
__________________
--
Mechanical Engineering
Sydney, Australia


Crank-Shaft is offline   Reply With Quote

Old   February 12, 2013, 03:33
Default
  #3
Senior Member
 
Crank-Shaft's Avatar
 
Ovi
Join Date: Oct 2012
Location: Sydney, Australia
Posts: 141
Rep Power: 4
Crank-Shaft is on a distinguished road
Hello everyone,

I am surprised by the lack of any comments on this thread however, for those who are still interested I will post the results here. I have a comparison PDF file and some screenshots of the contour plots. Please provide any guidance or suggestions for improvements. I have not been able to accurately model the temperature boundary conditions and incorporate the energy equations yet however, I am not sure that is causing the difference in the Cf dataset.

All educational comments are welcome.
Attached Images
File Type: jpg TotalPressure.jpg (15.1 KB, 13 views)
File Type: jpg VelU.jpg (14.3 KB, 18 views)
File Type: jpg VelV.jpg (14.9 KB, 16 views)
Attached Files
File Type: pdf lowReSSTComparison.pdf (16.3 KB, 7 views)
__________________
--
Mechanical Engineering
Sydney, Australia


Crank-Shaft is offline   Reply With Quote

Old   March 8, 2013, 10:31
Unhappy
  #4
Senior Member
 
Crank-Shaft's Avatar
 
Ovi
Join Date: Oct 2012
Location: Sydney, Australia
Posts: 141
Rep Power: 4
Crank-Shaft is on a distinguished road
Hey CFD Online Members,

Hope you are all going well with your individual goals and projects. I am still working with the above geometry and the Cf values are significantly different to the benchmark dataset. The Cp values have much lower errors and this has been a consistent trend for more than 4 different simulation cases.

I am trying to find root causes for this discrepancy and thought about many different sources of errors. Is it possibly due to the 'bump' being a no-slip shear wall and the wall immediately ahead of it specified as zero-shear? It may be having some problems at the transition as the free-stream velocity strikes the edge of the bump profile.

I have also observed that the continuity residuals are not falling below 10^-2 level and this remains unchanged for both transient and steady-state flow simulations. Does this mean that there is an imbalance of mass continuity across the domain and could this be causing the Cf values to vary greatly?

I have attached an image and I look forward to your comments regarding the source of the large Cf errors and some fundamental concepts related to the continuity residuals.
Attached Images
File Type: jpg Capture2.jpg (13.9 KB, 10 views)
File Type: jpg Capture.jpg (20.8 KB, 11 views)
__________________
--
Mechanical Engineering
Sydney, Australia


Crank-Shaft is offline   Reply With Quote

Old   March 8, 2013, 16:19
Default
  #5
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,104
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
Quote:
Originally Posted by Crank-Shaft View Post
I am trying to find root causes for this discrepancy and thought about many different sources of errors. Is it possibly due to the 'bump' being a no-slip shear wall and the wall immediately ahead of it specified as zero-shear? It may be having some problems at the transition as the free-stream velocity strikes the edge of the bump profile.
Did you try a "real" symmetry boundary condition instead of a zero shear wall in front of the bump?
They are not equivalent when it comes to turbulent quantities.

And wouldnt it make sense to use a free stream boundary condition at the inlet?
flotus1 is offline   Reply With Quote

Old   March 8, 2013, 18:27
Default
  #6
Senior Member
 
Join Date: Nov 2010
Posts: 128
Rep Power: 6
taxalian is on a distinguished road
Hi Crank-Shaft,
I would do the following:

--> first of all check your yplus value that should be less than one if you don't use any wall function, this should influence your cf values.
--> then try to do this computation case with density based solver i.e. compressible mode as you are using currently the incompressible solver. This might help too, as it is stated on the main page of this test case that the computations have been carried out using compressible solver.
--> please try to use appropriate boundary conditions.

hope this helps.

best regards and good luck.
taxalian is offline   Reply With Quote

Old   March 10, 2013, 14:39
Default
  #7
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,944
Blog Entries: 6
Rep Power: 39
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Quote:
Originally Posted by Crank-Shaft View Post
The current boundary conditions including the pressure outlet and inlet have an operating pressure of 101 325 Pa with a gauge pressure set to (101 325 x 0.02828) Pa. This pressure gradient results in a velocity of 69.7 m/s which satisfies the M=0.2 condition.
You have no pressure gradient from inlet to outlet !

Set inlet pressure to 104190.471 pa , outlet to 101325 pa and operating pressure to 0 pa.

You can also try with velocity inlet.


Quote:
Did you try a "real" symmetry boundary condition instead of a zero shear wall in front of the bump?
They are not equivalent when it comes to turbulent quantities.
How slip wall differs from symmetry when calculating turbulent quantities?

PS. Are you sure about geometry and mesh? Did you try the solution on meshes available on website...
Far is offline   Reply With Quote

Old   March 11, 2013, 05:19
Default
  #8
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,104
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
Quote:
Originally Posted by Far View Post
How slip wall differs from symmetry when calculating turbulent quantities?

At walls with zero shear stress, wall-functions for the turbulent quantities are applied anyway.
This is not the case for symmetry boundary conditions.

The effect can be seen e.g. in a 2D channel, where you set one of the walls to symmetry and the other one to zero shear stress wall.
flotus1 is offline   Reply With Quote

Reply

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
About compressible flow at low mach hit Main CFD Forum 2 October 26, 2009 22:21
compressible flow maria teresa FLUENT 1 September 7, 2007 16:58
urgent help needed with 2d compressible flow James FLUENT 2 June 20, 2007 04:22
Is an uniform flow possible for this geometry? Michael Hu FLUENT 1 April 13, 2006 22:11
compressible channel flow.. R.D.Prabhu Main CFD Forum 0 July 17, 1998 17:23


All times are GMT -4. The time now is 15:27.