|
[Sponsors] |
November 18, 2009, 16:07 |
|
#21 |
Member
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 17 |
Marco, one of my cases has 4 million mesh point. This one is going to take another 2 days. But I should have some results up from some other cases soon.
btw, just to be sure....what are your values for Ubar and nu? Thanks, Eelco |
|
November 18, 2009, 16:25 |
|
#22 |
Member
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 16 |
My Ubar is the same of "channelOodles" tutorial: Ubar = 0.1335, but I modify nu, I chose nu = 2*10^-5 to have Retau more or less = 180.
The previous pictures are adimensionalized with Retau or Utau, etc. so the results are generic. Best regards MT |
|
November 19, 2009, 09:31 |
|
#23 |
Member
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 16 |
Sorry, Eelco, could you post the set of your test case?
In my simulation I use: dimension (x*y*z): 3.5*2*2 number of cells (nz*ny*nz): 80*180*120 that means, for Retau = 200 (I use Ubar=0.1335, nu=2*10^-5), these adimensional units (x+, y+, z+): 8.75, *, 3.33 [*: means that I have a simple grading in y direction that is 10, so the dimension of the cell in the middle of the channel is 10 time the dimension of the boundary cell, so the range of y+={0.44 - 4.44} more or less]. If I well understand the set of your test case, your situation is more or less this: dimension (x*y*z): 2pi*1*pi number of cell (nz*ny*nz): 196*129*160 and if you have a Retau = 400 (that is more or less the same of channelOodles tutorial), you have these adimensional units (x+, y+, z+): 12.81, ?, 7.85 (?: means that I don't know if you put some grading). I can't use more cell because I have not more computational resource. bye MT |
|
November 19, 2009, 11:00 |
|
#24 |
Member
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 17 |
Marco, one of my cases just finished. They are taking longer than I thought. I checked my Retau and it is around 195, a little too high. I want to compare against the paper from Kim (1987) and he uses Retau equal to 180. I have to do some postprocessing now to see how it turned out. I will let you know.
|
|
November 19, 2009, 11:12 |
|
#25 |
Member
Eelco Gehring
Join Date: Mar 2009
Posts: 70
Rep Power: 17 |
Marco,
I have the velocity profile, u+ vs. y+ and the results are not good. I think I am missing something..... This profile is not with the fine mesh yet, this is with 96x129x80. |
|
January 20, 2011, 18:42 |
DNS and 1D spectra
|
#26 |
Member
Aldo Iannetti
Join Date: Feb 2010
Posts: 48
Rep Power: 16 |
Hi All,
Here attached a 1D assial spectra coming from a DNS symulation of a pipe using OF 1.6. Focus your attention on the higher wavenumbers. Is there anyone who knows the meaning of that strange behaviour? It seems there are some components caracterized by the high wavenumber that do not dissipate. Please post whatever you think in order to let the investigations go on. Let me know if something else is needed to understand |
|
January 27, 2011, 09:53 |
|
#27 |
Senior Member
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 15 |
Hi Aldo,
Is this spectrum in space or time? And what discretization do u use for the convective term? Regards, Steven |
|
January 27, 2011, 15:21 |
|
#28 |
Member
Aldo Iannetti
Join Date: Feb 2010
Posts: 48
Rep Power: 16 |
||
January 27, 2011, 15:31 |
|
#29 |
Senior Member
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 15 |
ok,
Did you use nyquist's criterion to determine the highest wavelength that you can represent on your mesh? How many cells do you have in the direction you took the spectrum? If your answer is to the first question is yes and to the second question is at least 200 you can start looking at your schemes. Do you use a full hex mesh? Do you use stretching? How much? If you use stretching and/or unstructered grids you can try either midPoint CDS or blended. Steven |
|
January 27, 2011, 15:41 |
|
#30 | |
Member
Aldo Iannetti
Join Date: Feb 2010
Posts: 48
Rep Power: 16 |
Yes, 300, full hex unstructured,
how to modify the fv schemes (for exemple)? Quote:
|
||
January 27, 2011, 15:58 |
|
#31 |
Senior Member
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 15 |
Just change:
Code:
div(phi,U) Gauss linear; Code:
div(phi,U) Gauss midPoint; Code:
div(phi,U) Gauss blended 0.95; blended supresses instabilities by using a small bit of upwind A question about physics: What is the smallest wavelength you expect based on your Reynold number and Kolmogorov theory? If your Reynolds number is really low there just might be no smaller eddies. |
|
January 27, 2011, 16:10 |
|
#32 | |
Member
Aldo Iannetti
Join Date: Feb 2010
Posts: 48
Rep Power: 16 |
Re_tao=200, Re=5500 I think is not due to the cell shape because I see the same behaviour in a similar simulation of a channel with perfect cells.
is midpoint useful for a perfect channel cells? (I've never used it) Blending with Upwind will introduce diffusivity.... Aldo Quote:
|
||
January 28, 2011, 02:48 |
|
#33 |
Senior Member
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 15 |
midPoint will be exactly the same as linear if you use cells without stretching
blended will indeed introduce diffusion I use sometimes either of these for DNS but that is on unstructured grids (polyhedral cells). Based on your Reynolds number you would expect physics at these wavelengths. What I find strange is that you have already dropped 6 orders of magnitude before the strange behavior in the at the high wavelengths start. Imagine where you would end up if the energy would keep dropping after say k = 45? About your mesh, how did you mesh a pipe with perfect hex cells? And where did you extract the spectrum? In the center of at the wall? |
|
January 28, 2011, 02:56 |
|
#34 |
Senior Member
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 15 |
This is one of my spectra for a DNS of a channel flow at Re_tau = 180.
Spectra extracted over the plane in the center of the channel. As you can see my mesh is 128 cells in streamwise direction so I have less information than you. Drop in magnitude almost the same. |
|
June 8, 2011, 03:29 |
|
#35 | |
New Member
Join Date: Oct 2010
Posts: 24
Rep Power: 15 |
Quote:
|
||
September 28, 2011, 18:23 |
|
#36 |
Member
Yong Wang
Join Date: Apr 2009
Posts: 34
Rep Power: 16 |
||
September 29, 2011, 08:14 |
|
#37 |
Member
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 16 |
Lijie here this is the link to the file that I wrote about:
http://www.cfd-online.com/OpenFOAM_D...S_tar-3861.unk I'm sorry for the late. MT |
|
October 6, 2011, 18:47 |
|
#38 | |
Member
Yong Wang
Join Date: Apr 2009
Posts: 34
Rep Power: 16 |
Quote:
Marco, have you ever tried to use some other div schemes? I am using my own solver in the frame of OpenFOAM. The size of the channel is 4pi*2*(4pi/3), and the mesh is 128*128*128 (uniform in the x and z directions, simpleGrading=10 in the y direction). My solver will crash when I use linear and midpoint as the div scheme. So, I am using Gamma0.1 and the solver works. Based on Jasak's paper about Gamma scheme, it should be similar with the central(linear) scheme. However, the velocity in the sublayer is smaller than Moser's DNS results. That means more dissipation has been introduced..... |
||
October 6, 2011, 21:53 |
|
#39 | |
Member
Yong Wang
Join Date: Apr 2009
Posts: 34
Rep Power: 16 |
Quote:
Steven, I am simulating the turb channel flow with Re_tau=180. I found your master thesis My solver will crash with linear/midpoint.....but works with blended0.995. hope I can get good result with it. |
||
October 11, 2011, 03:16 |
|
#40 | |
Member
Marco Turchi
Join Date: Jul 2009
Location: Pesaro, Italy
Posts: 31
Rep Power: 16 |
Quote:
my set case was: 3,5h*2h*2h; mesh: 80*180*120 (like you uniform the in x and the z direction, simpleGrading = 10 in the y direction); Re_tau=200 (more or less). I changed "fvSchemes" like this: ddtShemes: Euler; gradSchemes: Gauss linear (for all) 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 I hope this is useful. Marco |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Is GUI there for OpenFoam | navaladi | OpenFOAM | 4 | September 11, 2023 13:35 |
OpenFoam vs CFX5 mass balance in OpenFoam | tangd | OpenFOAM Running, Solving & CFD | 33 | May 23, 2010 16:36 |
[blockMesh] CheckMesh error using a tutorial from OpenFOAM 114 with openFOAM 13 | martapajon | OpenFOAM Meshing & Mesh Conversion | 7 | January 21, 2008 12: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 |