# beginner level - CG oscillations

 Register Blogs Members List Search Today's Posts Mark Forums Read

 January 11, 2006, 11:37 beginner level - CG oscillations #1 Peter Guest   Posts: n/a Sponsored Links Hi, I'm working on a CG solver for Navier-Stokes incompressible inviscous equations. Its purpose is purely visual. Advection - pressure cycle with some sources. If anyone has some idea about the following problem I would appreciate it. Jacobi relaxation works fine but CG solver produces strange effect, if speeds in scene exceed 4-6 voxels per frame. Like a frame to frame pulsation, in one frame the speeds are exagerrated, and another almost stopped. Sequence of images, played with every 2nd image, looks roughly continuous. Grid sizes range from 50^3 to 200^3. Basic unprecond. CG solver. Any number of iterations (10-200). Can it somehow converge to the wrong solution?

 January 11, 2006, 20:17 Re: beginner level - CG oscillations #2 diaw Guest   Posts: n/a >>>Can it somehow converge to the wrong solution? Yes, in my experience, the combination of scaling (dx, dy, dt) & flow velocity can push you in 2 directions. Sometimes, what you are observing is a 'mode' change & the solution wants to 'pull itself' towards one, or other mode. This can be as a result of the physics itself, or of the numeric scheme used. diaw...

 January 16, 2006, 06:55 Re: beginner level - CG oscillations #3 Peter Guest   Posts: n/a Is there anything that can be done, with preconditioning, or more precision (double vs float), or initial guess? Is it inevitable for things as such to happen if pressure/speed values go above some amount? If CG produces this, is it correct to guess that BiCG, or CGNR will behave the same way? Is there any paper on this subject of this alternate solution? _Sometimes_ the oscillations are smaller if number of iterations is set 10x more than usual. Thanks, Peter

 January 16, 2006, 09:42 Re: beginner level - CG oscillations #4 diaw Guest   Posts: n/a Hi Peter, The solutions are dictated by the physics, bc's & ic's of the problem you are simulating & the N-S themselves. Not so much by your choice of matrix solution technique, although a bad one will surely lead to odd results. The traditional approach has been to add convection-stabilisation &/ or artificial viscosity to prevent the solution becoming singular. In my experience, this constrains the flow into a viscous-dominated flow mode. It's ok for 'steady solvers' with very viscous, slow-moving flows. There appear to be other forms of flow modes - across the singularity - which require an acceleration term to be present - ie. unsteady flow solver. The mode-switch can happen very, very quickly, depending on the geometry & velocity. This next mode requires some cunning scaling to work through & understand. It isa lot of fun... to be sure. diaw...

 January 17, 2006, 11:43 Re: beginner level - CG oscillations #5 Peter Guest   Posts: n/a Thanks for your posts, diaw. Ona last thing, please Is there any book or article about this mode change you mention, and is there some other name for it? Peter PS I'll just post some more recent findings: * CG solution itself looks something like a standing wave, its pressure solution having positive and negative areas (30-50 voxes thick), which change in the direction of the wind source, and position of p=0 areas between doesn't depend much on wind speed, or grid topology. With every frame the sign changes on these areas: 1. +-+-+ 2. -+-+- 3. +-+-+ * if the grid is narrow (example 100*20*200) then everything is ok * if CG solver works on every 2nd frame, and plain Jacobi solver on frames in between, then everything is ok

 January 17, 2006, 18:20 Re: beginner level - CG oscillations #6 diaw Guest   Posts: n/a Books on 'mode change'... I doubt it... it's new research findings You may very well have standing wave phenomena - these can be a result of the physics, b/c's, i/c's, or of the numerics & discretisation used. The 'mode change' I refer to can be seen when certain flow structures intersect - there is a sudden change of flow regime, necessitating revised mesh & time steps, for example. This effect is a by-product of the physics, & the non-linearity of the N-S. In my experience, standing wave phenomena can occur in a physical situation, given the correct b/c's etc... even the old checkerboard pressure situation can represent a physical solution - a wave solution, not a bulk-flow solution. Perhaps you could look into your 'spatial sampling rate' a little more. Feel free to e-mail me some pics of your findings & we can discuss off-line. diaw...

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Clementhuon OpenFOAM Native Meshers: snappyHexMesh and Others 6 March 12, 2012 13:41 subhkirti OpenFOAM Native Meshers: snappyHexMesh and Others 0 October 13, 2010 04:10 andrea.pasquali OpenFOAM 5 February 8, 2010 04:08 steve Main CFD Forum 11 September 16, 2008 00:22 Vidya Raja FLUENT 7 September 26, 2006 20:00