
[Sponsors] 
July 10, 2009, 04:19 
LES by Multigrid

#1 
Member
M. Nabi
Join Date: Jun 2009
Posts: 38
Rep Power: 7 
I am trying to solve a turbulent flow by LES, using multigrid. I am using multigrid for Pressure correction solver. But to increase the CFL number (larger time steps), I am solve the momentum equations implicitly, and I am trying to use Multigrid to increase their convergence rate. Thus what I do:
1. solve the momentum equations by multigrid 2. solve the pressure by Multigrid 3. correct the velocity field But I recognised that the flow sometimes, tends to become laminar. Is it so that the Multigrid kills the turbulence, if we use it for the momentum equations? It is clear that Multigrid removes the high frequency errors on the coarser grids. Does it also removes the high frequency eddies, and transfer the results to the finer grids ? I thank you in advance 

July 10, 2009, 05:31 

#2  
Member
Jed Brown
Join Date: Mar 2009
Posts: 56
Rep Power: 9 
Quote:
Quote:


July 10, 2009, 09:02 

#3 
Member
M. Nabi
Join Date: Jun 2009
Posts: 38
Rep Power: 7 
Thank you Jed for your kindly help.
Actually, I am using semi implicit method. I am using AdamsBashforthCrankNicolson. Second order AB for the convective term, and CN for the diffusion. To solve the diffusion (implicitly), I am useing the convection term and a part of diffusion (because of CN) as a source term. I don't fix any value, it is the case I update the turbulent viscosity and the source term after each MG cycle. During this update, the convection term stay unchanged, but the diffusion term changes because of turbulent viscosity and also participation of other velocity components in the diffusion part. I have in total three MGcycles for the momentum, which are coupled on eachothers. after the cycle for u, I have a cycle for v, and then w. update the viscosity and source term and then repeat the cycles, untill convergence. I am still useing a small CFL number. In anycase, I am trying to don't exceed the smallest turbulent time scale. Does this kind of splitting removes the second order accuracy in time (and also kills the eddies)? If wel, is it possible to adjust the accuracy by applying third or forth order AdamsBashforth for the convection? What is your advise to keep the second order accuracy? Do you have any idea? Thank you very much. 

July 10, 2009, 10:01 

#4  
Member
Jed Brown
Join Date: Mar 2009
Posts: 56
Rep Power: 9 
Using a different integration method for different terms in the equation amounts to fixing certain values while solving for others. This is basically what it means to be semiimplicit. Apparently you are iterating something to convergence, but I suspect it's not the whole system. To be fully implicit, you should be able to write the whole thing as
Code:
F(x',x) = 0 Code:
G(x) = F(x0 + ax,x) = 0 If you are not actually solving these fullycoupled nonlinear systems on each time step (e.g. with Newton iterations), then a modified equation analysis of the integration scheme will have splitting errors dependent on various derivatives involving the fastest time scales in the system. The point of highorder splittings are to raise the orders of these splitting errors, but in the presence of multiple time scales, splitting schemes rarely achieve their design order (or even first order) once they step significantly over the stiff time scale. For more on this, a good start is Code:
@article{knoll2003bat, title={{On balanced approximations for time integration of multiple time scale systems}}, author={Knoll, DA and Chacon, L. and Margolin, LG and Mousseau, VA}, journal={Journal of Computational Physics}, volume={185}, number={2}, pages={583611}, year={2003}, publisher={Elsevier} } Quote:
If you do manage to step significantly over advective CFL while maintaining accuracy, I'd certainly like to hear about it. 

July 10, 2009, 12:29 

#6  
Member
M. Nabi
Join Date: Jun 2009
Posts: 38
Rep Power: 7 
Quote:
In streamwise, and spanwise direction, I am using peiodic boundary condition. in the normal direction, I am using solid boundary, and plansymmetry. No slip for the bed, and freeslip for the top. I treat the sgs term dynamically. I use the method of Meneveau (1996). I tried also smagoronisky model, using loglaw model. 

July 10, 2009, 12:45 

#7 
Senior Member

Maybe it could be the frictional step that dissipate too much energy....
...no, i'm joking, i couldn't resist . The fractional procedure is ok and, if correctly implemented, it has a second order splitting error. Almost all the academic and research LES codes uses the fractional step with the AB for the convective terms and CN for the viscous ones. Four issues remains: 1) What spatial schemes are you using for the convective and diffusive terms? 2) Have you checked the overall accuracy of your solver? (i.e., with the Taylor Vortex Solution) 3) The symmetry condition is not feasible in LES (this is probably the cause of the laminarization). 4) What kind of grid spacing (in wall units) are you using? 

July 11, 2009, 13:25 

#8  
Member
M. Nabi
Join Date: Jun 2009
Posts: 38
Rep Power: 7 
Quote:
1) I am using central differencing for the convection and diffusion terms. I don't use upwind. 2) Still I didn't check for turbulence. I checked the solver for cubic cavity flow. It gives very accurate results. 3) It may can be the reason !!!!. But I saw some papers, in which, they use this condition to approximate the free surface. 4) till now, I didn't check the grid spacing beside the wall. But I think I am using very higher than yplus. Can this factor be the reason for damping the turbulence? Thanks 

July 11, 2009, 20:21 

#9 
Senior Member

1) Ok
2) You should be sure that your code gives you the the designed accuracy. The boundary consitions should also be carefully checked. For example, i was using fluent with the nodebased gradient computation and biperiodic boundary conditions...then i discovered that fluent has an inconsistency with this particular combination and that a costant error was created, regardless of the time or grid step. Thus you should be sure that with your setting your code gives you the required accuracy (the error goes down at the proper rate). 3) If you're performing a simulation of the turbulent channel flow then you can't apply the symmetry boundary condition. I don't know about other cases but, generally speaking, in LES you can't take advantage of any symmetry of the flow. 4) Of course it depends on the flow and the boundary coindition you're using. However, guidelines exists for LES of wallbounded flows. In any case, a too coarse grid is never a good choice. The coarsest grid i've ever used for the turbulent channel flow was a 64x33x64 completely uniform grid with wall functions at a turbulent Reynolds number Ret=2000. Even with this coarse grid and a slightly dissipative scheme i had "turbulence". 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
LES and combustion model  Margherita Cadorin  CFX  0  October 29, 2008 06:24 
multigrid  sureshkumar  Main CFD Forum  0  June 9, 2006 01:20 
Differences between a laminar code and a les one  ben  Main CFD Forum  9  February 16, 2005 23:40 
Some Questions about LES.  Bin Li  Main CFD Forum  2  February 20, 2004 09:58 
Multigrid applied to ke models  Paulo Zandonade  Main CFD Forum  9  May 24, 1999 08:10 