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

Car aerodynamics

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   February 25, 2008, 16:31
Default It's a grid I made some time a
  #41
New Member
 
Mattia
Join Date: Mar 2009
Posts: 26
Rep Power: 9
morfeus80 is on a distinguished road
It's a grid I made some time ago because at the moment I don't have a good software. Howerver the geometry is very simple.

Tonight I'm iterating with HRV's settings, tomorrow I'll check if the solution is more realistic of that I got with the other settings.
morfeus80 is offline   Reply With Quote

Old   February 27, 2008, 12:46
Default Would any of you have access t
  #42
Member
 
Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 82
Rep Power: 9
juho is on a distinguished road
Would any of you have access to 3d-models of real cars you could spare? Something I could find Cd and Cl figures on?

I'd like to compare results to experimental data.
juho is offline   Reply With Quote

Old   February 29, 2008, 18:34
Default Hallo, I made 2000 iteration
  #43
New Member
 
Mattia
Join Date: Mar 2009
Posts: 26
Rep Power: 9
morfeus80 is on a distinguished road
Hallo,
I made 2000 iteration with HRV's settings, starting from the previous solution. Now the flow field is realistic, completely different from what I had obteined with the other settings.
But the CD I calculated with liftDrag is 0.367: too high! Has anyone obtained a better result? What schemes and solvers could I try?
morfeus80 is offline   Reply With Quote

Old   March 1, 2008, 03:08
Default Hi, In my not quite converg
  #44
Member
 
Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 82
Rep Power: 9
juho is on a distinguished road
Hi,

In my not quite converged but Cd and Cl quite stable results realizableKE give a Cd of 0.45 and Cl 0.3. kOmegaSST gives Cd 0.36 and Cl 0.28. As the car is imaginary I don't have anything to compare it to but the realizableKE seems high. kOmegaSST might be in the ballpark, or close enough not to be able to tell if it's wrong.

I'm using prisms 1-20mm from the surface and tetras elsewhere.

I've tried to use V2F but I'm having great trouble with my y+ ~ 1 mesh. Not least due to memory limitations of my computer's memory but also lacking any convergence.

The pressure problem was solved by increasing the size of the domain. I'm now using a volume of about 25x5x5m which gives a blockage ratio of 4%.
juho is offline   Reply With Quote

Old   March 1, 2008, 07:52
Default Added turbulent drag to the so
  #45
Member
 
Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 82
Rep Power: 9
juho is on a distinguished road
Added turbulent drag to the solver and iterated some more with the kOmegaSST:

pressureDragCoefficient = 0.357712
viscDragCoefficient = 0.00803725
turbCoefficient = 0.021208
turbDragCoefficient = 0.386957
LiftCoefficient = 0.2582
juho is offline   Reply With Quote

Old   March 25, 2008, 06:27
Default Dear users, we are going to
  #46
Member
 
Tommaso Lucchini
Join Date: Mar 2009
Posts: 83
Rep Power: 9
lucchini is on a distinguished road
Dear users,

we are going to organize a session about OpenFOAM developments and applications in the Automotive field at the next OpenFOAM workshop (www.openfoamworkshop.org), which will be held in Milano, 10-11 July 2008.

Among the possible subjects, we encourage presentations in the area of car aerodynamics, since there is a lot of interest for it.

In the typical "workshop spirit", also presentations on open problems and/or with preliminary results are welcome.

People interested to make a presentation can send an abstract by e-mail to me or to the other session organizers of the Automotive session. The deadline is on the 31st March.

Best Regards,

Tommaso
lucchini is offline   Reply With Quote

Old   April 8, 2008, 05:26
Default Hi everybody, Still working
  #47
Senior Member
 
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 281
Rep Power: 10
vinz is on a distinguished road
Hi everybody,

Still working on the topic I'd like to know if other people are getting result with the ahmed body for instance. If it's the case, would it be possible to know the schemes and turbulence models you are using?
Tommaso, is it possible to have access to the abstract of people presenting at the workshop? Titles on the website look interesting and I would like to know a bit more about some of them.

Regards,

Vincent
vinz is offline   Reply With Quote

Old   February 18, 2009, 10:36
Default Dear all, Like many on the
  #48
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Delft, Netherlands
Posts: 326
Rep Power: 12
tomf is on a distinguished road
Send a message via MSN to tomf Send a message via Skype™ to tomf
Dear all,

Like many on the forum, I am investigating the flow for the Ahmed body. Currently investigating various turbulence models for the 35 degree slant angle case. With the k-Omega-SST I am getting results which look reasonable, however I am questioning the amount of iterations I need for convergence.

Using simpleFoam I start with 300 iterations using Gauss upwind as suggested on the forum for my convective terms. I then switch to second order (limitedUpwind for turbulence, Gamma 0.2 for velocity). After 6000 iterations, my forces are not yet converged, as can be seen from the image attached. For clarification, the first 500 iterations are truncated, so it runs to 5500. In the first 500 iterations, there are some oscillations and high values, which dissapear rather quickly, but from then convergence is very slow.



My fvSchemes for the latter part of the iterations looks like:

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) cellLimited Gauss linear 1.0;
}

divSchemes
{
default none;
div(phi,U) Gauss GammaV 0.2;
div(phi,k) Gauss linearUpwind cellLimited Gauss linear 1;
div(phi,omega) Gauss linearUpwind cellLimited Gauss linear 1;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear limited 0.7;
laplacian((1|A(U)),p) Gauss linear limited 0.7;
laplacian(DkEff,k) Gauss linear limited 0.7;
laplacian(DomegaEff,omega) Gauss linear limited 0.7;
laplacian(1,p) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}

snGradSchemes
{
default limited 0.7;
}

fluxRequired
{
default no;
p;
}

While my fvSolution is as:
solvers
{
p GAMG
{
tolerance 1e-06;
relTol 0.1;

smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;

cacheAgglomeration true;

nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
};
U smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-8;
relTol 0.1;
};
k smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-8;
relTol 0.1;
};
omega smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-8;
relTol 0.1;
};
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
}

relaxationFactors
{
p 0.3;
U 0.7;
k 0.7;
omega 0.7;
}

PISO
{
nCorrectors 2;
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
}

To be complete, checkMesh gives as output:
Mesh stats
points: 1726824
faces: 15463989
internal faces: 15147131
cells: 7343044
boundary patches: 12
point zones: 0
face zones: 0
cell zones: 0

Number of cells of each type:
hexahedra: 0
prisms: 1237404
wedges: 0
pyramids: 1540
tet wedges: 0
tetrahedra: 6104100
polyhedra: 0

Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
body_top 14666 7543 ok (non-closed singly connected)
body_front 9357 4807 ok (non-closed singly connected)
body_bottom 18400 9499 ok (non-closed singly connected)
body_slant 4413 2310 ok (non-closed singly connected)
body_base 5116 2685 ok (non-closed singly connected)
body_struts 1500 799 ok (non-closed singly connected)
body_side 26563 13581 ok (non-closed singly connected)
tunnel_road 101873 51568 ok (non-closed singly connected)
tunnel_inlet 1713 1042 ok (non-closed singly connected)
tunnel_infinity 10052 5508 ok (non-closed singly connected)
tunnel_outlet 3088 1730 ok (non-closed singly connected)
tunnel_symmetryplane120117 66275 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-4.176 -5.46623e-12 0) (7.308 3.3265 3.3265)
Mesh (non-empty) directions (1 1 1)
Mesh (non-empty, non-wedge) dimensions 3
Boundary openness (-1.43894e-17 5.86639e-15 1.04078e-15) OK.
Max cell openness = 5.96728e-16 OK.
Max aspect ratio = 12.1381 OK.
Minumum face area = 7.52672e-08. Maximum face area = 0.0367279. Face area magnitudes OK.
Min volume = 1.12097e-11. Max volume = 0.00212606. Total volume = 99.7353. Cell volumes OK.
Mesh non-orthogonality Max: 67.1824 average: 16.1836
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 2.47713 OK.

Mesh OK.


I was wondering if anyone has a suggestion to speed up convergence, or do others have similar amounts of iterations until convergence? Obviously there are also some oscillations in the forces, could this indicate too severe unsteady behavior for the steady state solver, or should these damp out as well?

Any experience and comments would be welcome and appreciated.

Kind regards,
Tom Fahner
tomf is offline   Reply With Quote

Old   February 18, 2009, 11:50
Default Dear Tom, in my experience
  #49
Member
 
Luca Gasparini
Join Date: Mar 2009
Location: Italy
Posts: 37
Rep Power: 9
luca_g is on a distinguished road
Dear Tom,

in my experience Gamma with a value much less than 1 is not very stable.
The "Gauss linearUpwind..." scheme is quite more expensive than Gamma X, so unless you really noticed an improvement I would use Gamma on the turbulence fields as well.
Although at convergence it should make no difference, you may try using 1 or 2 nNonOrthogonalCorrectors to add robustness and maybe speedup convergence (in terms of iteratin count), since although your grid quality is high you still have quite some non orthogonality.

Finally, given the large separation area do not necessarily expect a convergence to a steady state; a kind of oscillating wake is likely to occur.

Regards,

Luca
luca_g is offline   Reply With Quote

Old   February 19, 2009, 06:08
Default Dear Luca, Thank you for yo
  #50
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Delft, Netherlands
Posts: 326
Rep Power: 12
tomf is on a distinguished road
Send a message via MSN to tomf Send a message via Skype™ to tomf
Dear Luca,

Thank you for your reply. I have been experimenting in 2D at first and found that Gamma did not behave stable for the turbulence, while with linearUpwind I had a bit more succes. So therefore I adapted to that scheme. I will try out Gamma there. I will also try adding some non-orthogonal correctors. Thanks again.

Regards,
Tom
tomf is offline   Reply With Quote

Old   April 3, 2009, 16:26
Default
  #51
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 226
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Tom,

I think you could relax the tolerance on P and you would observe a quicker convergence.

Regards,

-Louis
louisgag is offline   Reply With Quote

Old   June 9, 2009, 00:15
Default CFD for car
  #52
New Member
 
Elingselasri
Join Date: Jun 2009
Location: Indonesia
Posts: 6
Rep Power: 9
linkse_edogawa is on a distinguished road
Send a message via Yahoo to linkse_edogawa Send a message via Skype™ to linkse_edogawa
Hi, I am newbie in CFD,. wanna ask about the best domain to simulate flows around the car. (length behind, above, beside).ALso, how small the grid close to the surfaces if i want to obtain the drag coefficient for low velocity car (is that related?).
very well, thank you.
linkse_edogawa is offline   Reply With Quote

Old   June 14, 2009, 11:18
Default
  #53
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 226
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Hi Elingselasri,
I suggest you use the same lengths as done in the wind-tunnel case you are trying to replicate. As for the grid size, that depends on the precision of the results you're looking for and the computer power available to you. In 3D a y+ of 60 and more is generally ok for a car travelling at highway speed... Also, if you are new to CFD you should start with RANS models and avoid LES for now.

Good luck,

-Louis
joy2000 likes this.
louisgag is offline   Reply With Quote

Old   June 15, 2009, 01:36
Default
  #54
New Member
 
Elingselasri
Join Date: Jun 2009
Location: Indonesia
Posts: 6
Rep Power: 9
linkse_edogawa is on a distinguished road
Send a message via Yahoo to linkse_edogawa Send a message via Skype™ to linkse_edogawa
Thank LOuis.. but i want to simulate with farfield condition.
linkse_edogawa is offline   Reply With Quote

Old   June 15, 2009, 10:33
Default
  #55
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 226
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
still.. find a wind tunnel test to replicate
louisgag is offline   Reply With Quote

Old   June 20, 2009, 14:00
Default
  #56
Member
 
Join Date: Jun 2009
Posts: 36
Rep Power: 9
Xwang is on a distinguished road
What kind of hardware is needed to do these computation? 4gb of ram are enough? Is it possible to use a cartesian grid? Xwang
Xwang is offline   Reply With Quote

Old   June 20, 2009, 20:44
Default
  #57
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 226
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
It's better than nothing.. I run 2D simulations with 2Gb ram on two processors. However for my 3D simulations I'm switching to 2-3 computers with each 2gb ram and I still expect my CFD runs to be long..... Serious vehicle CFD people use at least 20+ computers, but if you just want to learn you should be fine...
louisgag is offline   Reply With Quote

Old   June 21, 2009, 02:49
Default
  #58
Member
 
Join Date: Jun 2009
Posts: 36
Rep Power: 9
Xwang is on a distinguished road
Quote:
Originally Posted by louisgag View Post
It's better than nothing.. I run 2D simulations with 2Gb ram on two processors. However for my 3D simulations I'm switching to 2-3 computers with each 2gb ram and I still expect my CFD runs to be long..... Serious vehicle CFD people use at least 20+ computers, but if you just want to learn you should be fine...
I want to learn the program. Then I will buy a more powerful pc. Xwang
Xwang is offline   Reply With Quote

Old   November 3, 2009, 02:31
Default
  #59
New Member
 
Jarno
Join Date: Nov 2009
Posts: 1
Rep Power: 0
Jarno is on a distinguished road
Greetings!

First of all, sorry my bad english language but hope you get idea what i am trying!

I am looking for right kind of software for linux. I have a complete 3d-model from my racing car and i need to put it on somekind virtual wind tunnel that i can test different wings and under body air flows. What is best program for this idea? I hawe made 3d-model of my racing car with Blender, and it can also opened with Paraview program. Should i install OpenFOAM and run some test with that? Does Paraview program has any plugs in that can help me?
By the way, progam should be free...
--Jarno.
Jarno is offline   Reply With Quote

Old   November 8, 2009, 13:25
Default
  #60
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 226
Rep Power: 10
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Hi Jarno,

you should be able to convert your model into a mesh that will work with OpenFOAM. Do a search on the forum for "blender" and you'll likely find the information you're looking for.

I have not tried it personally so I can't really help you there.

good luck,


-Louis
louisgag 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
Aerodynamics vengi FLUENT 5 October 25, 2011 10:43
Aerodynamics Bonny Jacob Zachariah Phoenics 3 February 10, 2009 05:43
CFD in aerodynamics Ujjwal Bhaskar FLUENT 1 December 26, 2007 11:29
Use of Pro-Am in aerodynamics Javidan Ahmad CD-adapco 8 December 3, 2004 00:27
unsteady aerodynamics R.KRISHNAMURTHY Main CFD Forum 1 December 6, 2000 02:17


All times are GMT -4. The time now is 13:56.