|
[Sponsors] |
Velocity curve in channel flow by pisoFoam and LES |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 19, 2013, 09:49 |
Velocity curve in channel flow by pisoFoam and LES
|
#1 |
Member
X Meng
Join Date: Jun 2012
Location: Scotland
Posts: 89
Rep Power: 14 |
Dear all,
I am struggling with a trouble in pisoFoam using LES for channel flow. Basically, the velocity curve from experiment has the maximum value of about 1.15 times of inlet uniform velocity, yeah? But I got much larger max value like 1.4 times of inlet uniform velocity, and it changes sharply from wall to centre. My question is, what parameters I may make wrong, or anything I can do to develop that? P.S. I am using oneEqEddy model. OF version 2.2.0. Cheers |
|
September 23, 2014, 10:14 |
|
#2 |
Member
Florian Ries
Join Date: Feb 2014
Location: Darmstadt, Germany
Posts: 88
Rep Power: 12 |
Hi,
that sounds like to coarse mesh. Which LES Model do you use??? Is your case really turbulent in the LES simulation. Is your case statistically stationary. Have you checked this?? kind regards Florian |
|
May 18, 2015, 09:19 |
|
#3 |
New Member
marzieh goodarzi
Join Date: Feb 2014
Location: Tehran-Iran
Posts: 11
Rep Power: 12 |
hi
I have the same problem but with pimple solver. I used oneEqEddy of LES model, and used mapped inlet Method, but the maximum velocity is 1.4 that exceed the experimental data. yplus of mesh is under 1 and about 162,000 node please help me ,I am confused Last edited by m.goudarzi; May 18, 2015 at 09:22. Reason: n |
|
May 19, 2015, 01:41 |
|
#4 |
Member
Florian Ries
Join Date: Feb 2014
Location: Darmstadt, Germany
Posts: 88
Rep Power: 12 |
hi m.goudarzi,
do you use Van Driest damping??? The OneEqModel behaves wrong at the wall. You have to use van Driest damping!! Or you can use dynOneEq (without Van Driest damping!), but there is a bug in this model Is your domain long enough??? Your domain has to be twice the longest turbulent structures (Nyquist criteria) i.e. in a pipe the largest turbulent structures is about 6-7 times the pipe diameter, this leads to a length of your domain of 16 times the diameter. Is your averaging time long enough?? How have you checked that?? Is your flow statistically statinoary before you start averaging?? How have you checked that? kind regards Florian |
|
May 23, 2015, 14:05 |
|
#5 | |||
New Member
marzieh goodarzi
Join Date: Feb 2014
Location: Tehran-Iran
Posts: 11
Rep Power: 12 |
Hi dear Florian
Thank you for the quick reply Quote:
Quote:
Quote:
Thank you marzieh |
||||
May 23, 2015, 14:21 |
|
#6 | |
New Member
marzieh goodarzi
Join Date: Feb 2014
Location: Tehran-Iran
Posts: 11
Rep Power: 12 |
This is my LES properties
Quote:
|
||
May 24, 2015, 06:13 |
|
#7 |
Member
Florian Ries
Join Date: Feb 2014
Location: Darmstadt, Germany
Posts: 88
Rep Power: 12 |
Hi,
how do you simulate the channel??? I would do it like this: 1. Use pimpleFoam (no relaxation in final iteration, do residualControl with very low values!) 2. use fvOptions with pressureGradientExplicitSource (no mapping method!) 3. initialize your flow with isotropic turbulence (boxTurb) and superimpose a powerlaw profile 4. For LES use Gauss linear (no Upwind!!), corrected (or orthogonal for channel) and backward (ddt) 6. Lets say flowdirection is z and the heigth of the channel is h. As a result of the turbulence structures, the domain has to be: Lz=12h, Lx=4h, Ly=h 7. Run your simulation for 100 flow through domain and then start averaging. (How to average see tutorial channel95) 8. Include a model for your subgrid-stress tensor Tsgs in the oneEq. You need an averaged Tsgs for postprocessing. How to compare LES data with DNS data: With LES you simulate a filtered velocity field. The mean value of a filtered signal is equal with the mean value of the unfiltered signal. So you can compare a mean value from a LES with a mean value of a DNS. This is not true for the RMS value. The RMS value of a LES is smaller then the RMS value of the DNS. To compare this, you have to add the mean subgrid stress tensor TSgs to the covariance matrix of the LES. Averaging time: The random error Err of the average depends on the integral time scale Tint. In example, for the mean your random errror is proprtinal: Err ~ sqrt(2*Tint/tav) where tav is your averaging time. Tint is a function of space and is much higher near the wall. As a result of this, you have to average much longer at the wall. You should do a reference simulation with Smagorinsky model. For channel use Smagorinsky parameter Cs=0.05. kind regards Florian |
|
June 4, 2015, 15:45 |
LES simulation of curved channel
|
#8 |
New Member
marjan
Join Date: Aug 2013
Posts: 11
Rep Power: 13 |
Hi everyone, I model a curved channel with N-S equation in curvilinear coordinates.
I imposed a fix pressure gradient in x- direction which decomposed along the channel in curvilinear coordinate( stream wise and span wise direction) to drive the flow. So the pressure gradient alter in steam wise and span wise direction along the channel . But my solution is not stable. so i Think it s because of unbalance between the body forces. Does some one has any idea? Should I use a fixed pressure gradient in stream wise while imposing the centrifugal force as well? |
|
June 15, 2015, 06:12 |
pisoFoam problems
|
#9 |
New Member
Gizela
Join Date: May 2015
Posts: 11
Rep Power: 11 |
Hi everybody, i´ḿ new to openFoam and i´m trying to run pisoFoam for a channel with a pier inside. I received an error (see image) and from that i tried to change somethings in my boundary conditions but still giving the same error.
So, please I need your help to find out what is going wrong, where is my failure and how can I do it right. it is attached some files (blockMeshDict, U, p, epsilon). help please. Kind regards Gizela Last edited by Gizazu; June 24, 2015 at 12:02. |
|
June 21, 2015, 08:28 |
|
#10 |
New Member
marzieh goodarzi
Join Date: Feb 2014
Location: Tehran-Iran
Posts: 11
Rep Power: 12 |
Hi,
Thank you very much Florian for your answer I have done all these steps, but the result did not change I do not know where is the problem |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Friction velocity in Channel flow | impecca | OpenFOAM Post-Processing | 2 | April 11, 2014 02:33 |
setup problems - LES pipe flow with cyclic BC (1) and direct mapped inlet (2) | florian_krause | OpenFOAM | 22 | June 13, 2013 22:25 |
channelFoam for LES of Flow over a Wavy Wall | sunliming | OpenFOAM | 4 | May 1, 2013 16:59 |
unsteady LES pisoFoam questions | garrison | OpenFOAM Running, Solving & CFD | 13 | August 17, 2012 21:12 |
calculating mean streamwise velocity in LES | maysmech | OpenFOAM | 4 | September 6, 2010 12:38 |