
[Sponsors] 
January 2, 2017, 08:42 
solver algorithm for LES

#1 
Member
raunak jung pandey
Join Date: Jun 2016
Posts: 64
Rep Power: 2 
Hello I am using Fluent 15.0 to carry out my LES simulations.
Which solver algorithm is appropriate for fast and accuracte results ? ITA/NITA/FSM/PISO/SIMPLEC My simulations run normally and converge when using SIMPLE/SIMPLEC and Coupled but takes long time My simulation shows follwoing error when using PISO and NITA Pressure Correction = 0.3 Moment Correction = 0.7 " divergence detected in AMG solver : pressure correction" 

January 2, 2017, 10:38 

#2 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,081
Rep Power: 35 
Try noniterative fractional step setting second order time and space accuracy


January 4, 2017, 09:29 

#3 
Member
raunak jung pandey
Join Date: Jun 2016
Posts: 64
Rep Power: 2 
[QUOTE=FMDenaro;631794]Try noniterative fractional step setting second order time and space accuracy[/QUOTE
Thank you for your answer. I am still facing convergence problem while using PISO and NITA. THis doesnt exist when I am using SIMPLE/SIMPLEC or Coupled and second order Implicit transeint formulation. WHat maybe the reasons for this ? 

January 4, 2017, 09:32 

#4 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,081
Rep Power: 35 
are you using the dynamic SGS model? Instability can happen for the eccessive value of the time step.
Could you plot the total kinetic energy versus the time? 

January 4, 2017, 10:34 

#5  
Member
raunak jung pandey
Join Date: Jun 2016
Posts: 64
Rep Power: 2 
Quote:
I am using LES WALE and Smagorinsky model. My time step size is 0.00001 sec giving courant no under 1. 

January 4, 2017, 12:24 

#6 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,081
Rep Power: 35 
why you can not? have you tried to reduce the CFL number of one order?


January 11, 2017, 07:59 

#7 
Member
Join Date: Jun 2016
Posts: 63
Rep Power: 3 
If you use PISO in combination with NITA, you should employ 23 neighbor corrections. Besides, a highquality mesh is desirable. NITA is naturally less stable than ITA since the coupling among equations is actually not enforced  there are only inner iterations, no outer iterations. Which implies that it is not applicable to all situations. The equation "faster computational method = the same stability & accuracy as the slower computational method" simply does not work
@FMDenaro  I am not aware of any simple method how to plot the turbulence kinetic energy in LES. Since it is defined via velocity fluctuations, one actually has to get there through instantaneous velocities and averaged velocities. To get averaged velocities means to save the instantaneous velocities at many time steps and to compute the mean and fluctuations after the calculation has finished (or have data from a previous calculation). Correct me if there is an easier way to get k. 

January 11, 2017, 10:35 

#8  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,081
Rep Power: 35 
Quote:
Just to be clear: I wrote total kinetic energy that is available from a LES solution. If the flow is incompressible then the solution v_bar allows to compute k_LES =v_bar . v_bar. Of course, that is not equivalent to k_bar. However, k_LES can be used to assess a statistical steady state or an energy pileup leading to numerical instability. what do you mean for: NITA is naturally less stable than ITA since the coupling among equations is actually not enforced  there are only inner iterations, no outer iterations. Which implies that it is not applicable to all situations. 

January 12, 2017, 06:14 

#9 
Member
Join Date: Jun 2016
Posts: 63
Rep Power: 3 
That is a pretty loose statement. I should have rather said that NITA is not suitable for every application. Obviously, strong coupling between equations would have adverse effects on convergence. I can imagine that to achieve a good convergence in multiphase flow modeling may be challenging using NITA. And just to quote the Fluent guide: the NITA solver is not recommended for highly viscous fluid flow.


January 12, 2017, 06:30 

#10 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,081
Rep Power: 35 
Despite of any high order accurate timeintegration, the fractional step method has an intrinsic error due to the timesplitting. I don't know much about multiphase flows but I do not understand the warning about viscous flows...


January 12, 2017, 08:17 

#11  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 627
Rep Power: 19 
Quote:
I created a fractional solver where I could include this correction and could demostrate that when including this nonlinear correction I could achieve results like iterative scheme for the test cases similar to what Milovan used (ie vortex shedding from cylinders). By switching this term off got results same as PISO. So yes, what you say has real effect and could be demostrated. This is why I created hybrid solver for multiphase flows where courant determines times step and using iterative scheme is expensive while NITA is unstable. I have been using this in FVUS for VOF calculations to cut the calculation cost by at least half. 

January 12, 2017, 09:31 

#12  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,081
Rep Power: 35 
Quote:
What type of correction have you implemented? Some years ago I deeply analysed the errors appearing in the fractional step method and the sources were highlighted in the splitting but also in the bc.s. Thus, some remedies were proposed to get full high order accuracy https://www.researchgate.net/publica...ary_conditions 

January 12, 2017, 10:00 

#13  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 627
Rep Power: 19 
Quote:
What I did was that I stored the momentum matrix and once the velocity correction was available I corrected flux and added another correction to source to momentum eqn. Then solved again (rather than setting up whole thing). I wrote down outline in a ppt, I try to find it and attach here. (it was long time ago so don't remember the steps exactly). Give me a bit of time to find that ppt. 

January 12, 2017, 10:11 

#14  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 627
Rep Power: 19 
Okay found the pdf.
slide 10 first chart is from Prof. Milovan's previous test. Second and third from my tests with on and off of this correction. This outlines the algo that I tried. Quote:


January 12, 2017, 10:48 

#15 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,081
Rep Power: 35 
No specific corrections in the BC.s? Have you checked the convergence slope on some analytical test?


January 12, 2017, 11:19 

#16  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 627
Rep Power: 19 
Quote:
Anyway did not work much on it but on a provided test case from industry it beat PISO of starCD by factor of 2 or so. (105 seconds vs 265 seconds of PISO in that test case). All in all I feel it is worth exploring but I have no time. May be in future I try it in FVUS, as for now there is lot to do there so no time. 

Tags 
les model 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Floating Point Exception Error  nyox  FLUENT  8  January 22, 2016 08:53 
Quarter Burner mesh with periosic condition  SamCanuck  FLUENT  2  August 31, 2011 11:34 
Working directory via command line  Luiz  CFX  4  March 6, 2011 21:02 
ns equation solver using simple algorithm  paresh.halder  FLUENT  0  December 3, 2010 17:00 
The correction on pressure equation of SIMPLE algorithm in MRFSimpleFOAM solver  renyun0511  OpenFOAM Running, Solving & CFD  0  November 10, 2010 02:47 