# diagonal dominance

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

 October 29, 2005, 22:54 diagonal dominance #1 khai ching Guest   Posts: n/a Hi all: I am using GS/SOR to solve the Pressure Equation in SIMPLE algorithm. The code works well for a simple square duct & backward facing step problem. But, when I try to use the same code to simulate flow in a forward-facing step, my code can't take it and diverges... I found that diagonal dominance is difficult to be achieved particularly on cells near to the forward-step (stagnation region). Can this be the main problem in achieving convergence? What is the remedy? Shall I use a "different" relaxation factor in this region in order to ensure convergence?? Sometimes, backflow may happen at outlet... Any special BC treatment if backflow is detected at outlet? So far, I have set the normal velocity to zero should backflow occurs..(treat it as wall).. However, this does not solve my problem.. Can anyone share your experience? -khai ching-

 October 30, 2005, 23:27 Re: diagonal dominance #2 diaw Guest   Posts: n/a Khai Ching wrote: I found that diagonal dominance is difficult to be achieved particularly on cells near to the forward-step (stagnation region). Can this be the main problem in achieving convergence? What is the remedy? Shall I use a "different" relaxation factor in this region in order to ensure conve rgence?? ------- A few thoughts: If you are using a steady-solver, the problem may very well be that you do not have a steady solution to the flow field you have constructed - at least at the velocity you are using. Try testing this by lowering the inlet velocity & try for convergence. If you obtain a solution, then gradually raise the velocity until your 'lack of diagonal dominance' begins to show itself. Then move over to a 'transient solver', or a steady solver which is 'soft in time', to take you through the transition point. Take 'snapshots' of your solution at various run intervals & look at where the problem occurs physically. This excercise can provide amazing insights into what the solver is really having to cope with. You can also try testing your model mesh sensitivity by refining the mesh in the region where you are experiencing problems. Once the 'correct mesh limit' has been reached, then you will most likely have to move towards an unsteady solver. diaw...

 October 30, 2005, 23:40 Re: diagonal dominance #3 khaiching Guest   Posts: n/a Thanks diaw.. I have reduced the inlet velocities and it works!! Perhaps I should gradually raise the velocity to ensure convergence... Thanks for your advise..

 October 31, 2005, 00:24 Re: diagonal dominance #4 diaw Guest   Posts: n/a It is my pleasure... good luck with your work... 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 santiagomarquezd OpenFOAM Programming & Development 5 November 3, 2014 15:07 nileshjrane OpenFOAM Running, Solving & CFD 4 February 25, 2013 05:13 ankgupta8um OpenFOAM Running, Solving & CFD 7 January 15, 2011 14:38 nileshjrane OpenFOAM Running, Solving & CFD 8 August 26, 2010 12:50 zonexo Main CFD Forum 1 May 10, 2006 13:43

All times are GMT -4. The time now is 17:11.