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

DNS with openfoam

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   October 7, 2008, 15:28
Default Has anyone been able to succes
  #1
Senior Member
 
N/A
Join Date: Mar 2009
Posts: 188
Rep Power: 8
harishg is on a distinguished road
Has anyone been able to successfully match the DNS results (especially Reynolds stresses) of Moser, Kim and Mansour at higher friction Re (395, 590) for channel flow? Is the numerical error of the cubic discretization available for the convection term comparable to spectral methods and high-order compact schemes, commonly employed for DNS computations ?

Thanks,
Harish
harishg is offline   Reply With Quote

Old   October 8, 2008, 03:45
Default Hi harish, I did some stuff
  #2
Senior Member
 
Cedric DUPRAT
Join Date: Mar 2009
Location: Belgium
Posts: 179
Rep Power: 8
cedric_duprat is on a distinguished road
Hi harish,

I did some stuff on that point.
you can have a look there : http://www.cfd-online.com/OpenFOAM_D...es/1/8544.html
If you have more questions, don't hesitate but ask them in the previous thread, to keep channel flow questions together :-)

Regards,

Cedric
cedric_duprat is offline   Reply With Quote

Old   October 8, 2008, 04:29
Default Hello Everybody My name is Su
  #3
Senior Member
 
Suresh kumar Kannan
Join Date: Mar 2009
Location: Luxembourg, Luxembourg, Luxembourg
Posts: 129
Rep Power: 8
kumar is on a distinguished road
Hello Everybody
My name is Suresh kumar, i am a P.hd student and i have started my study a month back at the University of Luxembourg. I have decided to use OpenFOAM for peforming some Multiphase calculations of a jet.
Since i have just installed OpenFOAM on my personal machine which is UBUNTU 8.04. I am facing some problems in installing the software.
I am posting the message that i get after i run the script OpenFOAM installation test, please give me some suggestions on resolving it


Checking basic setup...
-------------------------------------------------------------------------------
Shell: bash
Host: lareole
OS: Linux version 2.6.24-19-generic
-------------------------------------------------------------------------------


Checking main OpenFOAM env variables...
-------------------------------------------------------------------------------
Environment_variable Set_to_file_or_directory Valid Crit
-------------------------------------------------------------------------------
$WM_PROJECT_INST_DIR /home/kumar/OpenFOAM yes yes
$WM_PROJECT_USER_DIR /home/kumar/OpenFOAM/kumar-1.5 yes no
$WM_THIRD_PARTY_DIR /home/kumar/OpenFOAM/ThirdParty yes yes
-------------------------------------------------------------------------------


Checking the OpenFOAM env variables set on the PATH...
-------------------------------------------------------------------------------
Environment_variable Set_to_file_or_directory Valid Path Crit
-------------------------------------------------------------------------------
$WM_PROJECT_DIR /home/kumar/OpenFOAM/OpenFOAM-1.5 yes yes yes

$FOAM_APPBIN ...M-1.5/applications/bin/linuxGccDPOpt yes yes yes
$FOAM_USER_APPBIN ...r-1.5/applications/bin/linuxGccDPOpt yes yes no
$WM_DIR /home/kumar/OpenFOAM/OpenFOAM-1.5/wmake yes yes yes
-------------------------------------------------------------------------------


Checking the OpenFOAM env variables set on the LD_LIBRARY_PATH...
-------------------------------------------------------------------------------
Environment_variable Set_to_file_or_directory Valid Path Crit
-------------------------------------------------------------------------------
$FOAM_LIBBIN ...nFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt yes yes yes
$FOAM_USER_LIBBIN ...OpenFOAM/kumar-1.5/lib/linuxGccDPOpt yes yes no
$MPI_ARCH_PATH ...penmpi-1.2.6/platforms/linuxGccDPOpt yes yes yes
-------------------------------------------------------------------------------


Third party software
-------------------------------------------------------------------------------
Software Version Location
-------------------------------------------------------------------------------
WARNING: gcc version does not match gcc supplied with this release of OpenFOAM
Supplied version: 4.3.1
User version : 4.2.3
Minimum required: 4.2.0

gcc 4.2.3
WARNING: Conflicting installations:
OpenFOAM settings : /bin/gcc
current path : /usr/bin/gcc
CRITICAL ERROR

gzip 1.3.12 /bin/gzip
tar 1.19 /bin/tar
icoFoam
WARNING: Conflicting installations:
OpenFOAM settings : /home/kumar/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/icoFoam
current path :
CRITICAL ERROR

-------------------------------------------------------------------------------


Checking networking...
-------------------------------------------------------------------------------
Action Result Crit
-------------------------------------------------------------------------------
Pinging_lareole Successful yes
Pinging_localHost Successful yes
Test_rsh: Unsuccessful_connection_refused* yes
Test_ssh: Unsuccessful_connection_refused* yes
FATAL ERROR: No remote shell available.
OpenFOAM 1.5 enviroment requires either ssh and/or rsh.
Contact your system administrator.


-------------------------------------------------------------------------------

The system test has evoked 1 fatal error(s).

The foam installation contains 2 critical error(s).

Review the output for warning messages and consult
the installation guide for trouble shooting.

done.
if anybody could give me some suggestions it would be good. To be noted i have changed the settings of the compiler to SYSTEM in the file settings.sh.
bye
kumar is offline   Reply With Quote

Old   October 8, 2008, 06:38
Default Hi Suresh! At first: How is
  #4
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,914
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi Suresh!

At first: How is your question related to this thread? (see nr.5 on http://openfoamwiki.net/index.php/Ho..._Message_Board)

At your question: try running some tutorials. If they work, then the script is oversensitive.
@the remote shell: you won't need that unless you are doing parallel calculations

Bernhard
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   October 10, 2008, 08:03
Default Hi Harish, I would like to
  #5
Senior Member
 
John Deas
Join Date: Mar 2009
Posts: 160
Rep Power: 8
johndeas is on a distinguished road
Hi Harish,

I would like to discuss with you about DNS in OpenFOAM, but your email is not available. Could you contact me at john.deas@gmail.com ?

If anybody is interested by a direct contact on DNS to discuss the matter and to work out solutions together, do not hesitate.

JD
johndeas is offline   Reply With Quote

Old   October 10, 2008, 22:11
Default This should be an open communi
  #6
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
This should be an open community using an open platform for CFD, so, if there are no reasons to keep the discussion private as it seems, why don't making it here? Maybe readers and new users can benefit of that in the future!

Regards,
Alberto
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   September 3, 2009, 10:01
Default
  #7
Member
 
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 8
zebu83 is on a distinguished road
I have some question about channel DNS simulation with OpenFOAM.

I found in another thread [Channel flow DNS channelDNS program (that Alberto post long time ago), I compiled it and it work perfectly!
Then I created a new directory (channelDNS) and in this directory I run my simulation of a channel with this parameter:
L*H*W: 3.5*2*1
nx*ny*nz: 80*130*60
nu=4*10-5 (to have a Re_tau=180 more or less)

I use:
Quote:
ddtSchemes
{
default CrankNicholson 0.5;
}

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

divSchemes
{
default none;
div(phi,U) Gauss cubic;
}

laplacianSchemes
{
default none;
laplacian(nu,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
}

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

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p;
}
that is like fvSchemes of boxTurb tutorial, except of temporal schemes (CrankNicholson 0.5 instead of euler).

When I visualized the results I notice that in the middle of the channel there are some strange bubble, but they don't come from the wall streaks!

This cause a bed Umean and u_rms profile.

Where/what are the errors?
Probably interpolation schemes or resolution?

Thank you very much in advance.

Best regard
MT
zebu83 is offline   Reply With Quote

Old   November 13, 2009, 23:13
Default
  #8
Member
 
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 8
feijooos is on a distinguished road
Marco,

you mentioned that you used the channelDNS code and that it compiled, correct? I just downloaded this code from the forum, but I have some trouble compiling it. Would you mind sharing how you got it to work? I am using OF 1.5

Thanks,

Eelco
feijooos is offline   Reply With Quote

Old   November 14, 2009, 18:15
Default
  #9
Member
 
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 8
zebu83 is on a distinguished road
Good night Eelco,

what are your trouble? Can you give me more informations?

I am using OF 1.5.x, that has some corrections from 1.5 version. To download and install OF 1.5.x there are many posts that explain it.

You must read also all posts of the discussion that i quote before, there some interesting things.
Remember: channelDNS was written for OF 1.4 or 1.3 I don't remember, so some header file are now different. For example, in OF 1.5 there is a function that calculate the field average, so you don't need "createAverage.C".

I hope this instructions are useful for you.

MT
zebu83 is offline   Reply With Quote

Old   November 15, 2009, 13:36
Default
  #10
Member
 
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 8
feijooos is on a distinguished road
Hey Marco,

Yeah, I realize that it is not necessary to compile the createAverages.C and calculateAverages.C files, because those already get calculated in the 1.5 version. Alberto told me, that in 1.5 I can just select a 'laminar' model instead of LES and that will take care of it.

Btw, how are your comparisons against the DNS results from Kim (1987)?

Thanks,

Eelco
feijooos is offline   Reply With Quote

Old   November 15, 2009, 14:27
Default
  #11
Member
 
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 8
zebu83 is on a distinguished road
Very bad!

Not really very bad, but there is a little problem with "postChannel" utility: it makes a wrong statistic.
I notice that, if you want to calculate "urms" or "vrms" etc., postChannel uses first time average, then field average. I don't know if I explain it well...
I can't use it, postChannel, if I want correct statistics because the first fields are not good for statistics, I must use only stationary fields! You can see it if you look, for example, at Retau, if it is more or less constant it's ok, instead you must wait!

If it is not right, please correct me!
I loose some day for this things.

Then when I finish my simulation I would write a post about this.

Best regards

MT
zebu83 is offline   Reply With Quote

Old   November 15, 2009, 14:48
Default
  #12
Member
 
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 8
feijooos is on a distinguished road
Marco,

What was your resolution for these calculations?
I am using a domain size of: 2pi x 1 x pi, with 196*129*160 points.

Thanks,

Eelco
feijooos is offline   Reply With Quote

Old   November 15, 2009, 15:11
Default
  #13
Member
 
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 8
zebu83 is on a distinguished road
I modify "channelOodles" tutorial like this: 3.5 x 2 x 2 with 80*180*120, and I have Retau = 180 (more or less).
Your resolution is better then my resolution, but what about your Retau?

Are your results good? Is it like Moser or Marusic DNS?

The most evident difference are urms, my urms are 2.5 instead 2.8! But other quantities are not so bad. I thik that my channel is too small, and also finite volume are less conservative than spectral method, for example.

What do you thik about this?
I am interesting on it.
What I want to do it's a good channel DNS, but it is not so simple that I thought!!!

Bye

MT
zebu83 is offline   Reply With Quote

Old   November 15, 2009, 20:31
Default
  #14
Member
 
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 8
feijooos is on a distinguished road
My cases are still running. They should be done, I'll let you know what I get.

I was playing around with the postChannel utility, but it doesn't do anything. It doesn't write out anything for me. I checked it out in the utilities folder and it seems like everything is commented out. Have you checked this with 1.5?

Thanks, Eelco
feijooos is offline   Reply With Quote

Old   November 15, 2009, 20:47
Default
  #15
Member
 
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 8
feijooos is on a distinguished road
hmm I think I understand what is happening with postChannel. I see what you mean by using the first directory.

How do you calculate Retau, with utau equal to sqrt(gradP*delta/rho)?

Thanks, Eelco
feijooos is offline   Reply With Quote

Old   November 16, 2009, 07:27
Default
  #16
Member
 
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 8
zebu83 is on a distinguished road
I don't have my data, because I'm not in my office these day, but when I come back I'll post some result to explain better what I mean.

You can use "gradP" to calculate Retau, but only if you reach the stability. Instead you can calculate it by dU/dy (= wall grad U). There some utilities in OF for it, but I calculate it writing a little program myself.

I'll post all, when I come back in the office (tomorrow).

Have a nice day!
zebu83 is offline   Reply With Quote

Old   November 16, 2009, 16:23
Default
  #17
Member
 
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 8
feijooos is on a distinguished road
Thanks Marco, I'll try to get my things put some things together as well.

Eelco
feijooos is offline   Reply With Quote

Old   November 18, 2009, 06:01
Default
  #18
Member
 
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 8
zebu83 is on a distinguished road
Good morning Eelco.

I'm going to post some result, but first let specify something:
- if you want to calculate Retau with gradP, you have to do that: sqrt(gradP)*delta/ni, not
Quote:
sqrt(gradP*delta/rho)
. Because if you are in stationary condition you know that tau_wall = gradP, u_tau = sqrt(tau_wall/ro) [friction velocity] {pay attention here, OF gives you a gradP normalized with ro, so to calculate the friction velocity you don't have to divide by ro!}, then Retau = u_tau*delta/ni.
- if you use wall grad U (= dU/dy_wall) you can calculate Retau in any condition, because, as S. B. Pope [Turbulent flows, Cornell University, Cambridge, 2000] say, total stress tau = viscous stress - Reynolds stress: tau = ro*ni*dU/dy - ro*u'v'. Near the wall the Reynolds stress are zero, so you know that tau_wall = ro*ni*dU/dy_wall. You have only to calculate dU/dy_wall, for example (if you don't use OF post processing) with central differing (of course you must calculate yourself, using U field that OF saved at each time step).


These are two way to calculate Retau, of course the second it's always correct.

I'm going to post some pictures:


This is Umean log profile. You can notice that in the "log region" 30<y+<110 it don't reach the log law or not so much.


This is the stress picture: viscous stress, Reynolds stress and total stress. As you can notice, the result is not the straight line pass through y=-1.


These are the rms, as I say in a previous post, they are underestimated (if you look for example at the Moser or Marusic papers).


Turbulent structures.

As you can see from the pictures the result are quite good, but not at all! This is because I use "postChannel" in this case.
The flux is totally turbulent, as you can see from the last picture, and it also reach the stability, but when I try to calculate statistics postChannel uses all the data, because it makes first time average, then spatial average.
The correct way to do it is to calculate first spatial average at every field itself, then make the time average from the every spatial average. If you do like that you receive better results, but in my case I have a too small channel and my "rms" are not so good. It is due also to finite volume method, because it introduces a numerical diffusivity.

I hope it's all clear.

Ask me if you don't understand and, please, correct me if I am wrong.

Best regard

MT
zebu83 is offline   Reply With Quote

Old   November 18, 2009, 12:09
Default
  #19
Member
 
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 8
feijooos is on a distinguished road
Hey Marco,

Thanks for your post, this is really helpful. I can see what you mean, you are capturing the overall structures, but not the exact values. I agree that some of this is caused by numerical diffusion, not using a much more accurate spectral method.

I am currently running some cases with different mesh sizes. I will try to post-process these in the next few days. I will post my results with some images after I'm done with this.

Thanks for your help! Eelco
feijooos is offline   Reply With Quote

Old   November 18, 2009, 17:05
Default
  #20
Member
 
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 8
zebu83 is on a distinguished road
It's going to be very useful if you will be able to post some of your result!

Best regards

MT
zebu83 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
OpenFoam vs CFX5 mass balance in OpenFoam tangd OpenFOAM Running, Solving & CFD 33 May 23, 2010 16:36
CheckMesh error using a tutorial from OpenFOAM 114 with openFOAM 13 martapajon OpenFOAM Native Meshers: blockMesh 7 January 21, 2008 13:52
OpenFOAM users in Munich OpenFOAM benutzer in M%c3%bcnchen jaswi OpenFOAM 0 August 3, 2007 13:11
A new Howto on the OpenFOAM Wiki Compiling OpenFOAM under Unix mbeaudoin OpenFOAM Installation 2 April 28, 2006 08:54
Is GUI there for OpenFoam navaladi OpenFOAM 3 March 28, 2006 09:08


All times are GMT -4. The time now is 17:24.