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

A STARTING POINT FOR DNS

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

Like Tree1Likes
  • 1 Post By grtabor

Reply
 
LinkBack Thread Tools Display Modes
Old   September 13, 2006, 23:11
Default The DNS solver requires inputs
  #1
New Member
 
Andy
Join Date: Mar 2009
Location: Fukuoka, Japan
Posts: 5
Rep Power: 8
andy_andy is on a distinguished road
The DNS solver requires inputs for Ea and K0 in "BOXTURBDICT" and for U0sigma, U0alpha, U0kupper and U0lower in "turbulence properties".

I want to simulate a turbulent flow for a known integral length scale and turbulent rms velocity and then apply it to a combustion event. I assume I can estimate these inputs from this and the thermodynamic state of the fluid. How do I go about that?

I have found some bits around and could have a guess at their meaning, however its always best to ask especially as it takes me weeks to find a solution. I would like to find the exact definitions of these parameters from literature or a website. Can anybody help?

I am happy to update the WIKI once I understand it all. I am sure others having the same problems.
andy_andy is offline   Reply With Quote

Old   September 24, 2006, 21:48
Default I`d welcome any input from oth
  #2
New Member
 
Andy
Join Date: Mar 2009
Location: Fukuoka, Japan
Posts: 5
Rep Power: 8
andy_andy is on a distinguished road
I`d welcome any input from others who have been through the same process.

Here is my progress to date.

The turbulence spectrum appears to be based on the following equation.

E(k)= (Ea*(k/k0)**4)exp(-2.0*(k/k0)**0.5)

That was easy enough.

However, the inputs for "turbulenceProperties" appear to be more difficult to define. They appear to be linked to the forcing of turbulence.

I have traced the terms U0sigma, U0alpha, U0kupper and U0lower back to UOprocess.C but cannot find the definition of each of these terms or any literature describing them.

My best guess at the moment is that U0kupper and U0lower are the maximum and minimum wavelengths used to excite the turbulence. I still cannot work out what U0sigma and U0alpha are exactly. Any help would be greatly appreciated.
andy_andy is offline   Reply With Quote

Old   September 25, 2006, 03:22
Default UO process uses an Uhlenbeck-O
  #3
Senior Member
 
Gavin Tabor
Join Date: Mar 2009
Posts: 181
Rep Power: 8
grtabor is on a distinguished road
UO process uses an Uhlenbeck-Ornstein stochastic process to generate a stream of random numbers - or rather not quite random as they are stochastic, and so each number is partly related to the one preceeding it. Sigma and alpha are the two parameters in the UO process, controlling the variance and mean (I think) respectively. Have a look at

http://en.wikipedia.org/wiki/Ornstein-Uhlenbeck_process

for more details. Unfortunately there are various different definitions of the process with different scalings for alpha and sigma, and I'm not clear which reference was used for the original algorithm; I always have to sort through the code to work out the precise definition.

Gavin
Z.Q. Niu likes this.
grtabor is offline   Reply With Quote

Old   November 9, 2006, 22:06
Default Thank you for your help. I
  #4
New Member
 
Andy
Join Date: Mar 2009
Location: Fukuoka, Japan
Posts: 5
Rep Power: 8
andy_andy is on a distinguished road
Thank you for your help.

I am still having problems with the DNS though. My solution always becomes unstable. Below is a summary of my problem.

Geometry 5mm x 5mm square 2 "cyclic" patches and an "empty" for the 3rd dimension

256x256 grid size

dT= 0.0001s

nu= 1.59e-5

Ea= 0.00056
K0= 10

UOsigma = 0.1
UOalpha = 0.8
U0Kupper = 300
U0Klower = 7

My K and Epsilon slowly grow until there is no solution. I have tried reducing the timestep and tuning each input to no use. Any ideas?

Andy
andy_andy is offline   Reply With Quote

Old   November 9, 2006, 22:15
Default Are you by any chance doing a
  #5
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Are you by any chance doing a 2-D DNS simulation? This will (of course) not work, as the physics is all wrong. If you are doing 3-D as you should, you should not be using empty patch type.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   February 21, 2008, 13:15
Default boxTurbDict seems to never be
  #6
Senior Member
 
John Deas
Join Date: Mar 2009
Posts: 160
Rep Power: 8
johndeas is on a distinguished road
boxTurbDict seems to never be accessed. I am a bit confused about the magnitude of the forcing in the current version.
johndeas is offline   Reply With Quote

Old   March 10, 2008, 08:23
Default Hi, I am trying to create a bo
  #7
ali
Member
 
Ali Heidari
Join Date: Mar 2009
Location: Surrey, London, United Kingdom
Posts: 39
Rep Power: 8
ali is on a distinguished road
Hi, I am trying to create a box with initial turbulent. I have the box dimention and Vrms, I don't know how to relate them to Ea and K0. the formula which is presented here for turbulence energy spectrum [E(k)= (Ea*(k/k0)**4)exp(-2.0*(k/k0)**0.5)] is not clear for me, I will be very grateful if some one can clue me in on this problem.
ali is offline   Reply With Quote

Old   February 13, 2012, 05:01
Default
  #8
New Member
 
Terrence Nguyen
Join Date: Jan 2012
Posts: 11
Rep Power: 5
hismother is on a distinguished road
Hi Andy,

are you still working on this. Since now I just begin with DNS and am facing the same problem as yours here. Would be nice, if you can show me how to overcome this.

Thanks in advance,

Terrence
hismother is offline   Reply With Quote

Old   July 30, 2013, 09:51
Default
  #9
Member
 
Jack
Join Date: Aug 2012
Posts: 42
Rep Power: 4
JackW is on a distinguished road
Hi all,

I know this is old, but I would love to know if anyone got anywhere,

I am guessing that sigma and alpha are the variance and mean, but I do not know how they relate to the taylor Re. This is what I would LOVE to know!

johndeas,
boxTurbDict controls the application boxTurb which generates the initial conditions on velocity.
This I don't think is too important as the forcing will eventually send the simulation to the correct Re.

Andy,
The upper limit on K seems very large, as in you are forcing too much of the spectrum and this could well be the simulation is blowing up.

Jack
JackW 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
Point coordinates via orinial point address after decomposing diana OpenFOAM Running, Solving & CFD 0 August 6, 2008 10:27
hex mesh by starting from a group of point Giulia Valerio FLUENT 0 November 8, 2007 05:47
floating point error when starting to iterate Ralf Schmidt FLUENT 1 June 2, 2006 06:00
Granular starting point not at a boundary Michael FLUENT 0 February 22, 2006 08:30
Starting Point to study CFD and code writing Lynn Main CFD Forum 3 November 23, 2005 05:31


All times are GMT -4. The time now is 10:44.