
[Sponsors] 
June 13, 2005, 17:08 
large number of grids, oscillation

#1 
Guest
Posts: n/a

I'm numerically solving continuity, momentum and energy equations for imcompressible flow. The scheme is one order accurate both in time and spacing. Theoretically, the decreasing grid size will give the more accurate results as long as the CFL condition is satified. But when I used large number of grids, the unexpected results of oscillation come out. Could anybody explain why this happen? I don't think the CFL is the reason. Is this common?


June 14, 2005, 04:14 
Re: large number of grids, oscillation

#2 
Guest
Posts: n/a

It is nonlinear instability, stability of your scheme was analyzed for linear equations where CFL is enough, when you increse number of cells, number of steps is increased and the influence in regions of big gradients is essential and it destroys ths solution


June 14, 2005, 05:41 
Re: large number of grids, oscillation

#3 
Guest
Posts: n/a

It sounds like your grid size is now minimizing numerical viscosity to the extent where you're resolving transient flow features which normally manifest themselves as periodic variation in residuals. A coarser grid would have sufficient numerical viscosity to damp out these transient features.
To check this, monitor the solution of a velocity component at the point of the max residual error, you might see this velocity oscillating periodically as well. 

June 14, 2005, 17:48 
Re: large number of grids, oscillation

#4 
Guest
Posts: n/a

Thanks a lot, it's quite enlighting. any reference concerning these topic?


June 15, 2005, 03:29 
Re: large number of grids, oscillation

#5 
Guest
Posts: n/a

No, just my experience!


June 15, 2005, 04:51 
Re: large number of grids, oscillation

#6 
Guest
Posts: n/a

find Quirk's article "Debate on great Riemann solver..." or try search "carbuncle phenomina" problems are similar m.a.


June 16, 2005, 02:01 
Re: large number of grids, oscillation

#7 
Guest
Posts: n/a

I'm trying to understand the logic in previous postings that with increased mesh density, that the 'diffusive effects' of the large mesh is removed & that the 'true' underlying flow filed begins to emerge  leading to oscillations.
Is there perhaps a chance that with increased mesh density, that numerical errors could begin increasing  leading to oscillations? Or, is there perhaps something inherent in the FV technique that could contribute to these issues? This last comment is based on conversations with colleagues involved in FEM work who noted that in comparason studies against St**CD, as the mesh desnity was increased, that the error began to also increase. They attributed this in part to the 'checkerboard scheme'. diaw... 

June 16, 2005, 04:41 
Re: large number of grids, oscillation

#8 
Guest
Posts: n/a

Hi diaw,
I'm quite interested in your last comment. Can you please explain a bit more of how the checkerboard scheme could cause such residual oscillation? I get similar problem in some cases where I try to use the more sophisticated turbulence model (RSM) and so I agree with you that the diffusive effect could have cause such a problem. 

June 16, 2005, 04:46 
Re: large number of grids, oscillation

#9 
Guest
Posts: n/a

hmm....It seems to deal with gas dynamics & high speed flow where shock wave may occur...Sounds like different stories for me..Haven't fully read through it anyway. The turbulent model used could have contributed some instabilities to the solution as well.


June 16, 2005, 05:50 
Re: large number of grids, oscillation

#10 
Guest
Posts: n/a

Pete,
The thing that has always concerned me with the use of staggered grids is that solutions are occuring which are essentially offset from each other. The u,v,w velocity & scalar fields are calculated & only interlink via interpolations between their nodes  say at the cell faces of the scalar cell. Often, I have seen that a simple arithmetic average is used to obtain face fluid properties  eg. a centraldifference approach as applied to the property itself. What happens when the property is itself varying? We are also dragging the pressure field up with us as we settle the velocity field. The techniques for doing this are, in themselves, subject to approximations & reliance on iteration to solution. With more iteration comes the risk of additional numerical error. I have always inherently felt that there was scope for problems at some point. With collocated schemes, with everything being calculated 'in syncronous' in one cell, I feel that there is less potential for longterm intrinsic instabilities to occur. One must also not forget, that the basic NavierStokes equations upon which we so stringly depend  are in themselves subject to approximations & simplifications. My inner feeling is that we have simplified out the inherent natural damping present in true fluids & now are becoming slaves to false sources of damping to settle the solutions. The experts in our community can possibly shoot holes in my arguments  but, these 'personal feelings' are based on deep investigations & numerous simulations.  The colleagues I refered to in the FEM community rely heavily on multigrid techniques & errorestimation. In the early days of the project they evaluated as many schemes as they could obtain data for. St**CD was burned in their memory for its increasing error with decreasing meshdimension (increasing mesh density). diaw... 

June 16, 2005, 09:32 
Re: large number of grids, oscillation

#11 
Guest
Posts: n/a

/********************************/ /*FLUENT PERSPECTIVE*/ /********************************/
Hi, I have a similar experience with my solution. I am trying to model bubbly flow in an upwardly pipeflow. When I have 10 nodes in radial direction I observe lesser number of oscillations in the solution where as when I increase nodes to 15 and 20 it becomes worse. Well what I infered with little reading is the pressure velocity coupling in FLUENT needs to undergo revision. One can refer paper by Kunz et al., (1997)where some similar behaviour is discussed. However for an end user it might be a bit hard to implement his reccommendation, as one has to tweak and modify the source code which might not be possible with UDF. CITATION: Kunz R. F., B. W. Siebert. W. K. Cope, N. F. Foster, S. P. Antal and S. M. Ettore, " A Coupled Phasic Exchange Algorithm for Three Dimensional Mutlifield Analysis of Heated Flows with Mass Transfer", Comp. and Fluids 27, 741768 (1998). LINK TO THE PAPER http://www.sciencedirect.com/science...&_coverDate=09%2F30%2F1998&_cdi=5694&_orig=search&_st=13&_sort=d &view=c&_acct=C000051270&_version=1&_urlVersion=0& _userid=1069263&md5=34214d96d9ea71a0f7ab620b036441 ab KP 

June 16, 2005, 10:12 
Re: large number of grids, oscillation

#12 
Guest
Posts: n/a

I think you have it the wrong way round. For convection dominated problems the usual FE arrangement or FD/FV scheme which stores velocity and pressure in the same location has too many unknowns. The staggered arrangement reduces the number of velocity components so that we have the same number of unknowns as equations.
In order to see the above consider the 1D mass and momentum equation for incompressible flow with 2 cells/elements and central differencing for the convection term and no diffusion term. The result of too many unknowns in collocated schemes is the checkboarding/uncontrolled modes/... or however one wishes to describe it. The solution is to add something unphysical such as pressure smoothing, momentum averaging, biased differencing,... or the like. 

June 17, 2005, 08:32 
Re: large number of grids, oscillation

#13 
Guest
Posts: n/a

The staggered arrangement was designed  as I understand it  to get around the 'checkerboard pressure field' problems...
Refer to Versteeg "An introduction to computational fluid dynamics  The Finite Volume Method"  Chapter 6  pg 136+ The description & workings are well laid out. Then move onto discussion of the Simple algorithm etc  pg 142. diaw... 

June 22, 2005, 03:19 
Re: large number of grids, oscillation

#14 
Guest
Posts: n/a

Hi,
Can you please pass me the full details of Kunz et al paper? Thanks. 

June 22, 2005, 10:20 
Re: large number of grids, oscillation

#15 
Guest
Posts: n/a

Pete,
Here is the citation; Kunz R. F., B. W. Siebert. W. K. Cope, N. F. Foster, S. P. Antal and S. M. Ettore, " A Coupled Phasic Exchange Algorithm for Three Dimensional Mutlifield Analysis of Heated Flows with Mass Transfer", Comp. and Fluids 27, 741768 (1998). Let me know if u do not have acess to it. KP 

June 22, 2005, 20:04 
Re: large number of grids, oscillation

#16 
Guest
Posts: n/a

I would concur that the pressurecoupling issue can be a problem in many FV schemes. This is what I was alluding to in a previous post.
Thanks for the paper link. diaw... 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Large Eddy Simulation with too high Mach number  Roland R  CFX  4  November 25, 2010 17:43 
DecomposePar unequal number of shared faces  maka  OpenFOAM PreProcessing  6  August 12, 2010 09:01 
air bubble is disappear increasing time using vof  xujjun  CFX  9  June 9, 2009 07:59 
BlockMeshmergePatchPairs  hjasak  OpenFOAM Native Meshers: blockMesh  11  August 15, 2008 07:36 
Unaligned accesses on IA64  andre  OpenFOAM  5  June 23, 2008 10:37 