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

Sphere at Re=1e+7

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 13, 2018, 08:06
Smile Sphere at Re=1e+7
  #1
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
Hi,
i'm trying to simulate a 3D-Sphere in water at Reynolds Number 1e+7 (i.e. 10'000'000) in Star-CCM+, but at the moment i'm not achieving the correct CD result of about 0.19. My physics options are:

Velocity Inlet V = 10.014 m/s
Density rho= 998.6 kg/m3
Sphere Diameter D = 1m
Turbulence Model = K-w
Time = Steady
Fluid Type = Liquid.

The Mesh options on the body are:
Number of Prism Layer = 10
Stretching Factor = 1.5


Can someone please help me?
I really appreciate much.

-Danny-
danny9019 is offline   Reply With Quote

Old   June 13, 2018, 08:25
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by danny9019 View Post
Hi,
i'm trying to simulate a 3D-Sphere in water at Reynolds Number 1e+7 (i.e. 10'000'000) in Star-CCM+, but at the moment i'm not achieving the correct CD result of about 0.19. My physics options are:

Velocity Inlet V = 10.014 m/s
Density rho= 998.6 kg/m3
Sphere Diameter D = 1m
Turbulence Model = K-w
Time = Steady
Fluid Type = Liquid.

The Mesh options on the body are:
Number of Prism Layer = 10
Stretching Factor = 1.5


Can someone please help me?
I really appreciate much.

-Danny-



Computing viscous drag at such Re number require a very fine grid. You cannot model the wall BC if your goal is computing the viscous forces
FMDenaro is offline   Reply With Quote

Old   June 13, 2018, 08:30
Default
  #3
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
So? What i can do?
Thanks very much.
danny9019 is offline   Reply With Quote

Old   June 13, 2018, 08:42
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by danny9019 View Post
So? What i can do?
Thanks very much.
What about your present grid resolution at the wall?
FMDenaro is offline   Reply With Quote

Old   June 13, 2018, 09:03
Default
  #5
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
The first element near to the body is 0.0004m lenght. The Total prism layer thickness is 0.05m.
Thanks again
danny9019 is offline   Reply With Quote

Old   June 13, 2018, 09:08
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by danny9019 View Post
The first element near to the body is 0.0004m lenght. The Total prism layer thickness is 0.05m.
Thanks again



you need to compare this height to the BL thickness, a fact expressed by the y+ value. To compute accurately the viscous stress at the wall you need at least 3-4 nodes at y+<1.
FMDenaro is offline   Reply With Quote

Old   June 13, 2018, 09:11
Default
  #7
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
So i have to obtain a y+ value less than 1?
danny9019 is offline   Reply With Quote

Old   June 13, 2018, 09:15
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by danny9019 View Post
So i have to obtain a y+ value less than 1?

at least 3-4 nodes at y+<1
FMDenaro is offline   Reply With Quote

Old   June 14, 2018, 06:52
Default
  #9
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
Hi FMDenaro, i've setted the grid in this way:

Number of prism layer = 22
Total prism layer lenght = 5 cm

Achieving an y+ maximum of 0.92 like in the attached picture.

Now my Cd is 0.12
What can i improve to rise the Cd value to the right 0.19?
Thanks so much.

danny9019 is offline   Reply With Quote

Old   June 14, 2018, 06:58
Default
  #10
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by danny9019 View Post
Hi FMDenaro, i've setted the grid in this way:

Number of prism layer = 22
Total prism layer lenght = 5 cm

Achieving an y+ maximum of 0.92 like in the attached picture.

Now my Cd is 0.12
What can i improve to rise the Cd value to the right 0.19?
Thanks so much.


Still refining the grid close to the wall. Try also different model
FMDenaro is offline   Reply With Quote

Old   June 14, 2018, 07:02
Default
  #11
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
I'm running the simulation with the "All y+ treatment". Do you think it is better use "Low y+ model?

How can i refine the grid close to the wall? Increasing the number of layer? Because with 40 the Cd is still 0.12.

Many thanks.
danny9019 is offline   Reply With Quote

Old   June 14, 2018, 08:51
Default
  #12
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by danny9019 View Post
I'm running the simulation with the "All y+ treatment". Do you think it is better use "Low y+ model?

How can i refine the grid close to the wall? Increasing the number of layer? Because with 40 the Cd is still 0.12.

Many thanks.



More layers in the range y+<1 so that you can better describe the BL. Try to compare the results obtained with other models.
FMDenaro is offline   Reply With Quote

Old   June 14, 2018, 09:20
Default
  #13
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
In order to do this, i have to split the sphere in the regions where y+ is less than 1 and more than 1, right? In order to construct two different boundary layer?

Thanks again
danny9019 is offline   Reply With Quote

Old   June 14, 2018, 09:25
Default
  #14
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
That depends on the grid generation code you are using. Using RANS you can exploit the symmetry of the geometry to solve only a half sphere, using structured grid in the BL and unstructured grid in the outer region.
FMDenaro is offline   Reply With Quote

Old   June 14, 2018, 09:48
Default
  #15
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
I'm currently using RANSE and analyzing a quarter sphere exploiting the double symmetry. If i have get you i've to create two different boundary layer where: one more refined with more layer And the other less refined with less layers.
I have to obtain a value of y+ costant along the body right?
Many thanks FMDenaro
danny9019 is offline   Reply With Quote

Old   June 14, 2018, 09:52
Default
  #16
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by danny9019 View Post
I'm currently using RANSE and analyzing a quarter sphere exploiting the double symmetry. If i have get you i've to create two different boundary layer where: one more refined with more layer And the other less refined with less layers.
I have to obtain a value of y+ costant along the body right?
Many thanks FMDenaro



yes, you need to retain 3-4 grid points below y+=1 along the whole body
FMDenaro is offline   Reply With Quote

Old   June 14, 2018, 10:01
Default
  #17
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Still refining the grid close to the wall. Try also different model
So in the regions where is Yellow/ red i have to create a refined mesh and where is blue a coarse mesh?
Many thanks.
danny9019 is offline   Reply With Quote

Old   June 14, 2018, 10:04
Default
  #18
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Unfortunately, I don't see any figure...
FMDenaro is offline   Reply With Quote

Old   June 14, 2018, 10:07
Default
  #19
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
In The picture attached before this morning
danny9019 is offline   Reply With Quote

Old   June 14, 2018, 11:32
Default
  #20
Member
 
Danny
Join Date: Jun 2018
Location: Rome, Italy
Posts: 36
Rep Power: 8
danny9019 is on a distinguished road
Quote:
Originally Posted by danny9019 View Post
Hi FMDenaro, i've setted the grid in this way:

Number of prism layer = 22
Total prism layer lenght = 5 cm

Achieving an y+ maximum of 0.92 like in the attached picture.

Now my Cd is 0.12
What can i improve to rise the Cd value to the right 0.19?
Thanks so much.

This picture FMDenaro
danny9019 is offline   Reply With Quote

Reply

Tags
drag coefficient, mesh 3d, starccm+


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Possible to select Cells next to a SPHERE surface? keepfit OpenFOAM 2 November 30, 2016 16:47
[snappyHexMesh] Sphere in a channel by snappyHexMesh arsalan.dryi OpenFOAM Meshing & Mesh Conversion 0 November 14, 2016 17:01
[ICEM] meshing a sphere - large deviation from perfect sphere murx ANSYS Meshing & Geometry 25 August 15, 2012 13:37
CFX analysis of a sphere - drag too low Rob Findlay CFX 6 March 26, 2007 11:11
meshing F1 front wing Steve FLUENT 0 April 17, 2003 13:37


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