
[Sponsors] 
September 12, 2015, 06:06 
Residuals of continuity not converging

#1 
New Member
Ren Zetian
Join Date: Sep 2015
Posts: 5
Rep Power: 8 
Hello ,
I'm modeling a flow through a 2D channel with a rectangular bluffbody of Channel Diameter = 25 mm Bluff body Diameter = 5 mm Velocity inlet of cold flow = 1.6 m/s Outlet is set outflow I intend to get a steady profile of velocity and pressure as starting conditions for other calculation. The "viscouslaminar" model is turned on, and used air with constant density. Pressurevelocity coupled scheme is selected. However , the residual of continuity does not converge after even 10000 iterations. And the velocity magnitude and the pressure profile seems asymmetric. the results change little after the first several hundreds iterations. convergence criteria for continuity, x,y,z velocity = 1e3. Any help or views on this would be greatly appreciated. Thanks! 

September 12, 2015, 06:44 

#2 
Senior Member
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 1,016
Rep Power: 24 
The problem here could be that your system is implicitly unsteady, so the continuity residual doesn't decrease.
The profile is not symmetric because it continues to change, this steady solution should be just a screenshot of an unsteady case. I suggest to switch from this result to unsteady simulation: I'm quite sure you will obtain convergence of continuity too.
__________________
Google is your friend and the same for the search button! 

September 12, 2015, 06:55 

#3 
New Member
Ren Zetian
Join Date: Sep 2015
Posts: 5
Rep Power: 8 
Thank you.
But even for an unsteady case, it is expected to get a timeaveraged profile with the steady solver, right? As can be seen from the pics, the pressure and velocity magnitude profile is some what timeaveraged, though not exactly symmetric yet. 

September 12, 2015, 06:56 

#4  
New Member
Ren Zetian
Join Date: Sep 2015
Posts: 5
Rep Power: 8 
Thank you.
But even for an unsteady case, it is expected to get a timeaveraged profile with the steady solver, right? As can be seen from the pics, the pressure and velocity magnitude profile is some what timeaveraged, though not exactly symmetric yet. Quote:


September 12, 2015, 06:58 

#5 
Senior Member
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 1,016
Rep Power: 24 
As I know, the steady solver solves for a steady solution: if the system physically doesn't reach a steady solution the solver can't converge.
If you want a "steady solution" you should perform unsteady simulation and then get averagetime profiles of what you want in post processing.
__________________
Google is your friend and the same for the search button! 

September 12, 2015, 15:16 

#6 
Member
amirhossein
Join Date: Jul 2014
Location: iran
Posts: 81
Rep Power: 9 
well first thing you have to do is that increase the diameter of channel , 25mm is small for it , increase up to 50mm
after that you should check your reynolds number , i think it is not laminar and also this problem is unsteady pay attention to time step to get correct converge 

September 12, 2015, 23:21 

#7  
New Member
Ren Zetian
Join Date: Sep 2015
Posts: 5
Rep Power: 8 
Thank you.
The geometry may not be changed. The Re based on the bluff body diameter is 500, and 2500 for the channel diameter. But upstream of the bluff body, the flow has fully developed. You also think that a transient solver should be applied? Quote:


September 13, 2015, 03:41 

#8  
Member
amirhossein
Join Date: Jul 2014
Location: iran
Posts: 81
Rep Power: 9 
Quote:
i am sure about geometry that should be increase , diameter 25 is not enough for this problem 

September 16, 2015, 00:26 

#9 
Senior Member
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 4,585
Rep Power: 53 
There's no guarantee that a steadystate solver will converge to the timeaveraged solution of an unsteady problem. If your problem is implicitly unsteady (or even worse explicitly unsteady) then your solution at different iterations of a steady state solver can still oscillate. In this case it's best to use an unsteady solver.


October 5, 2015, 10:27 

#10  
Member
wanghuo
Join Date: Aug 2014
Posts: 89
Rep Power: 9 
Quote:
What's the difference of implicitly unsteady and explicitly unsteady? 

September 9, 2020, 11:12 

#11  
Member
Sai Krishna
Join Date: May 2018
Posts: 33
Rep Power: 5 
Quote:
iam facing same issue with residuals not getting converged. iam doing conjugate heat transfer analysis for the external flow over my model using steady state formulation. i created monitors for temperature @ locations where re circulation happens and other important parts. those values are completely stable(constant). Net mass flow rate from fluxes also converged fully. But the residuals are oscillating about a particular value. continuity @ 10e1, k and w @ 10e4, all velocities and energy @10e6. u told about implicitly unsteady system for which we cant get residual convergence with steady simulation, what does it actually mean? i know its depends on physics of the problem, but how to identify it in the beginning? if its known in the beginning, we can go ahead with transient simulation right? can there be a system which doesnt reach physically steady state? I believe even in compressible high speed flows will have some unsteady behavior and perturbations in the beginning and reaches steady state in due course. Thanks for your answers 

Tags 
continuity, converge 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Help for the small implementation in turbulence model  shipman  OpenFOAM Programming & Development  25  March 19, 2014 11:08 
Micro Scale Pore, icoFoam  gooya_kabir  OpenFOAM Running, Solving & CFD  2  November 2, 2013 14:58 
AMI interDyMFoam for mixer nu problem  danny123  OpenFOAM Programming & Development  8  September 6, 2013 03:34 
dynamic Mesh is faster than MRF????  sharonyue  OpenFOAM Running, Solving & CFD  14  August 26, 2013 08:47 
Convergence moving mesh  lr103476  OpenFOAM Running, Solving & CFD  30  November 19, 2007 15:09 