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

Ahmed body simulation

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

Like Tree1Likes
  • 1 Post By Far

Reply
 
LinkBack Thread Tools Display Modes
Old   May 30, 2012, 12:18
Default Ahmed body simulation
  #1
New Member
 
alessandro
Join Date: Oct 2010
Posts: 28
Rep Power: 6
alenglaro is on a distinguished road
Here the continuation of the thread about meshing of Ahmed body
alenglaro is offline   Reply With Quote

Old   May 30, 2012, 12:23
Default
  #2
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,916
Blog Entries: 6
Rep Power: 39
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
Mesh Quality


Do you think that the turbulence intensity has such a huge impact on drag?
Far is offline   Reply With Quote

Old   May 30, 2012, 12:28
Default
  #3
New Member
 
alessandro
Join Date: Oct 2010
Posts: 28
Rep Power: 6
alenglaro is on a distinguished road
I don't know how much is relevant but with default 10% cd was about 0.35. With 0.25% 0.289.
alenglaro is offline   Reply With Quote

Old   May 30, 2012, 12:33
Default
  #4
New Member
 
alessandro
Join Date: Oct 2010
Posts: 28
Rep Power: 6
alenglaro is on a distinguished road
Quote:
Originally Posted by scipy View Post
Did you read any of the Ahmed body papers? Usual turbulence intensity of the free stream airflow in a windtunnel is around 0.25 %, so me using 1 % is already a bit higher and more "real world"-like, but 2.5 is unnecessary.

You (and most others) seem to confuse a no slip wall and a no shear wall. No slip wall is what I've used for the ground and the ahmed body (since that's what they are, stationary walls with shear stress or 0 slip). However, for the side and top wall - if you are going after ultimate accuracy and recreating exact wind tunnel conditions, then those walls should also have a boundary layer (same as the ground and the Ahmed body itself), but since most of the time they're far enough away to only affect the solution a little bit (if left as no slip), they can be given no shear stress boundary conditions so their viscosity effects will be disregarded. Since this is mathematically the same as a symmetry BC for Fluent, I used that for simplicity.

/edit, seems you have read about turbulence intensity
[/QUOTE] Sorry, 2.5% was a typing mistake

About no-slip wall and no-shear wall, I knew the the difference, I was just wondering if using one instead of the other could affect the solution. But if it's for simplicity now I try with symmetry
alenglaro is offline   Reply With Quote

Old   May 30, 2012, 13:50
Default
  #5
New Member
 
alessandro
Join Date: Oct 2010
Posts: 28
Rep Power: 6
alenglaro is on a distinguished road
Simulation completed
Here settings used:
-Coupled with pseudo-transient, first order on momentum and turbulence kinetic energy and dissipation rate for 25 iterations and then switch to second order
- explicit relaxation factors to default except for Pressure and momentum at 0.4 and turbulent viscosity to 0.95.
-Automatic time-step method with Timescale Factor=10
-bcs with 1% turbulence intensity and symmetry also on top and lateral surfaces
-fmg initialization with default settings

Convergence in 168 iterations and 55 minutes ( about 20 seconds per iteration)
Cd=0.27444
Cl=4.2129e-02
(unfortunatly i forgot to hard copy cl and cd histories)
maximum residuals at 10^-5
Attached Images
File Type: jpg ahmed_far_fmg_pseudo.jpg (47.1 KB, 64 views)
alenglaro is offline   Reply With Quote

Old   May 30, 2012, 14:07
Default
  #6
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,916
Blog Entries: 6
Rep Power: 39
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
did you compare to experimental values? How much error is there?

What is the slant angle in current simulation and are you going to try the URANS?

These results are taken with mesh3? Try to refine mesh to 4 million and see the results.
Far is offline   Reply With Quote

Old   May 30, 2012, 14:53
Default
  #7
New Member
 
alessandro
Join Date: Oct 2010
Posts: 28
Rep Power: 6
alenglaro is on a distinguished road
I've found different experimental values for cd, from 0.25 to 0.295. For the LSTM case i found two values: 0.279 and 0.285. The ercoftac case 9.4 on which I rely doesn't report cd but only velocity profiles and pressure coefficients in some locations along the body. I have to compare those values but it will take some time. For cl I found a value of 0.004.
I used 35 configuration. I don't think I will perform a URANS simulation, I'll rather try different turbulence models to match velocity profiles as much as possibile.
alenglaro is offline   Reply With Quote

Old   May 30, 2012, 14:58
Default
  #8
New Member
 
alessandro
Join Date: Oct 2010
Posts: 28
Rep Power: 6
alenglaro is on a distinguished road
Results are taken on mesh 3. Do you think wake zone should be refined? And at the beginning of the slanted surface? That point should be critical for flow separation
alenglaro is offline   Reply With Quote

Old   May 30, 2012, 15:01
Default
  #9
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,916
Blog Entries: 6
Rep Power: 39
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
http://css.engineering.uiowa.edu/~me...medcarrodi.pdf

I think mesh is similar to used by Rodi in above paper. However it is always necessary to check the mesh independence.
Far is offline   Reply With Quote

Old   May 30, 2012, 16:16
Default
  #10
New Member
 
alessandro
Join Date: Oct 2010
Posts: 28
Rep Power: 6
alenglaro is on a distinguished road
I checked yplus on veicle surface and it's between 0.05 and 1.05. Is is necessary to have such low values? Shouldn't it be sufficient between 1 and 5?
alenglaro is offline   Reply With Quote

Old   May 30, 2012, 18:27
Default
  #11
Senior Member
 
Join Date: Aug 2011
Posts: 114
Rep Power: 6
scipy is on a distinguished road
Send a message via Skype™ to scipy
The most recent paper I could find that listed Cd and Cl for the Ahmed body was "Experiments and numerical simulations on the aerodynamics of the Ahmed body" by W. Meile, G. Brenn, A. Reppenhagen, B. Lechner, A. Fuchs in which they reported the experimental Cd of 0.279 and Cl of 0.004 for the 35 slant.

This means the agreement of Cd with the experimental value was within 0.95 % and Cl was unfortunately an order of magnitude higher at 0.044 (11x higher). However, the same paper lists their numerical figures at 0.276 for Cd and 0.013 for Cl which is again 3.25x higher than experiment.
scipy is offline   Reply With Quote

Old   May 30, 2012, 22:21
Default
  #12
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,916
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:
I checked yplus on veicle surface and it's between 0.05 and 1.05. Is is necessary to have such low values? Shouldn't it be sufficient between 1 and 5?
Yes, you are correct. Y+ between 1 and 5 is sufficient.
Far is offline   Reply With Quote

Old   May 31, 2012, 10:48
Default
  #13
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,916
Blog Entries: 6
Rep Power: 39
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
few updates are available here

NEW : Tutorial on road vehicle
alenglaro likes this.
Far is offline   Reply With Quote

Old   May 31, 2012, 17:15
Default
  #14
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,916
Blog Entries: 6
Rep Power: 39
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
I have used "Higher order term relaxation" technique to accelerate the convergence.
Far is offline   Reply With Quote

Old   December 14, 2012, 12:44
Default
  #15
Member
 
Li Huang
Join Date: Jan 2011
Posts: 30
Rep Power: 6
lihuang is on a distinguished road
Quote:
Originally Posted by Far View Post
few updates are available here

NEW : Tutorial on road vehicle
Hi Far,

Thanks for the tutorial. But i am having trouble to create a full car body mesh when i mirror your blocking. Basically what i did was, first, mirror the geometry and blocking, and set associations for some edges and vertices. Everywhere else looks good, except the part showing in the attached pic. The problem is i have some mesh lines converged to some points. But i have checked the Pre-mesh params for that edge and parallel ones to make sure proper params. How can I fix this part of mesh?
Thanks for your help!
Attached Images
File Type: jpg full body.jpg (63.8 KB, 39 views)
lihuang is offline   Reply With Quote

Old   December 14, 2012, 13:04
Default
  #16
Member
 
Li Huang
Join Date: Jan 2011
Posts: 30
Rep Power: 6
lihuang is on a distinguished road
Quote:
Originally Posted by lihuang View Post
Hi Far,

Thanks for the tutorial. But i am having trouble to create a full car body mesh when i mirror your blocking. Basically what i did was, first, mirror the geometry and blocking, and set associations for some edges and vertices. Everywhere else looks good, except the part showing in the attached pic. The problem is i have some mesh lines converged to some points. But i have checked the Pre-mesh params for that edge and parallel ones to make sure proper params. How can I fix this part of mesh?
Thanks for your help!
Hi Far,
The problem was just solved by remove association of that edge in problem. Don't know why it happened tho.
Thanks again!
lihuang is offline   Reply With Quote

Old   December 14, 2012, 13:06
Default
  #17
Far
Super Moderator
 
Far's Avatar
 
Sijal Ahmed Memon (turboenginner@gmail.com)
Join Date: Mar 2009
Location: Islamabad Pakistan
Posts: 3,916
Blog Entries: 6
Rep Power: 39
Far will become famous soon enoughFar will become famous soon enough
Send a message via Skype™ to Far
.................
Far 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
Ahmed Body Simulation nick FLUENT 4 May 10, 2010 10:40
Can someone please e-mail me geometry files for an Ahmed Body? ha85 ANSYS Meshing & Geometry 0 May 1, 2010 18:38
LES simulation of bluff body at Re(d)=1200 arash_nl FLUENT 0 June 1, 2009 13:16
Submarged body simulation in Fluent 6.3.26 Behnam FLUENT 0 April 22, 2009 06:10
Can't get good results on ahmed body research Michael Main CFD Forum 8 September 17, 2008 13:55


All times are GMT -4. The time now is 09:23.