# Continuity convergence issue on a 3d wing

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

 March 21, 2012, 23:25 Continuity convergence issue on a 3d wing #1 Senior Member   Nick Join Date: Nov 2010 Posts: 126 Rep Power: 8 Sponsored Links Hi all, I'm modelling a section of a 3d wing with a c-grid topology using the k-kl-w transitional model in Fluent. I have velocity-inlet and pressure-outlet BCs and the side planes are periodic. I'm using Simple Quick as the solution method. Starting off with a steady solver, the residuals of continuity go down to 10^-3 and then they remain in this neighbourhood. I then switch to the unsteady solver with a time step of 10^-6 and residuals go down to 5*10^-5 however this is not a good convergence level for me as the lift and drag do not stabilize and the flow is not developed. My main issue is with the continuity equation, although the same mesh works well in terms of convergence in CFX, I struggle in Fluent. My mesh quality is decent and I have already tried extending the boundaries further away from the wing body. Also I've tried using the Coupled solution method and decreasing the timestep and even using symmetry side planes but continuity doesnt converge. I'd appreciate any suggestions to resolve this issue. Nick

 March 22, 2012, 04:50 #2 Senior Member   Join Date: Mar 2011 Location: Germany Posts: 431 Rep Power: 13 Hi Nick you need to explain your case in detail. First I would like to ask a question for myself...you are working with structured or unstructured grid?? in which package you have created the grid?? ICEM or GAMBIT or something else?? Firstly use outflow at the exit with velocity inlet or use pressure inlet with pressure outlet, secondly what's your Reynold No??what's your Mach No?? at what Angle of Attack you are performing your analysis?? do you really need to switch to unsteady solver?? BTW 5*10^-5 is not bad at all for the continuity equation..... for how many no. of time steps you have allowed unsteady solver to run before assuming that lift and drag are not converging?? explain these few points then i hope i will be able to help you Regards

 March 22, 2012, 08:46 #3 Senior Member   Nick Join Date: Nov 2010 Posts: 126 Rep Power: 8 Thanks for your response. My AOA is 4 degrees and the structured hex mesh was created in ICEM. Reynolds is 100,000 (chord-based). The steady solver only goes down to 10^-3. The unsteady solver's behavior doesn't change much with the number of time steps, it basically goes down to 5*10^-5. I am comparing two transitional models the SST transitional against K-Kl-w. The SST transitional does better in terms of continuity in the steady solver since it goes down to 10^-4 in fluent. The flow is incompressible BTW. Also my yplus is below 1 everywhere on the wing.

 March 24, 2012, 01:13 #4 Senior Member   Lucky Tran Join Date: Apr 2011 Location: Orlando, FL USA Posts: 1,880 Rep Power: 26 Nick, how are you determining that the lift and drag are not stablized and that the flow is not developed? Is there any other issue besides continuity residuals? Continuity residuals, or any residuals for that matter are poor estimators of convergence. I would double check the way residuals are defined for each program (scaled vs normalized, and their actual definitions) if you are concerned about residuals in different programs.

 March 25, 2012, 03:45 #5 Senior Member   Nick Join Date: Nov 2010 Posts: 126 Rep Power: 8 Thanks for your answer. I can't quite follow what you mean by residuals aren't a good criterion for convergence. I monitor lift and drag for judging convergence as well as the residuals. Continuity is the only equation which I am struggling with. Also I was wondering if anyone has had any experience with the k-kl-w transitional model and come across a similar issue. I'd also appreciate it if you could inform me of the solution method with the highest resolution in Fluent.

March 25, 2012, 03:54
#6
Senior Member

Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 1,880
Rep Power: 26
Quote:
 Originally Posted by Nick R Thanks for your answer. I can't quite follow what you mean by residuals aren't a good criterion for convergence. I monitor lift and drag for judging convergence as well as the residuals. Continuity is the only equation which I am struggling with. Also I was wondering if anyone has had any experience with the k-kl-w transitional model and come across a similar issue. I'd also appreciate it if you could inform me of the solution method with the highest resolution in Fluent.
Hi Nick, same question again since you did not quite answer it
how are you determining that the lift and drag are not stabilized and that the flow is not developed?

residuals are not a good criterion to judge convergence. especially continuity because of the way it is calculated.

monitor solution values. it is obvious that in an unsteady simulation, the lift and drag are also time-varying quantities. so how did you determine that they are not converging? did you compare instantaneous lift and drag or averaged? if averaged, how did you averaged and for how long?

also, since it is an unsteady simulation, for which time step is 5e-5? the information you are providing is not very clear and it is not even certain if there is a problem at all with your simulation.

March 25, 2012, 05:32
#7
Member

Join Date: Nov 2011
Posts: 83
Rep Power: 7
Quote:
 Originally Posted by LuckyTran residuals are not a good criterion to judge convergence. especially continuity because of the way it is calculated.
Dear LuckyTran,

I was reading this post and came across this sentence. Could you please argument some more on that? I would be really interested in knowing more. How is mass flow residual calculated? Why it is not reliable? Why generally residuals are not a good criterion to judge convergence?
Thanks a lot.
Cheers,

ROb

March 25, 2012, 14:07
#8
Senior Member

Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 1,880
Rep Power: 26
Quote:
 Originally Posted by robboflea Dear LuckyTran, I was reading this post and came across this sentence. Could you please argument some more on that? I would be really interested in knowing more. How is mass flow residual calculated? Why it is not reliable? Why generally residuals are not a good criterion to judge convergence? Thanks a lot. Cheers, ROb
You can search residuals and get lots of hits.

Residuals are a measure of how much imbalance is left in each cell volume of the quantities, continuity (mass imbalance), momentum (in x,y,z). If additional models are used, there are more residuals for each equation (k,e or k,w for example, and energy).

Usually the raw residuals are not reported. The scaled or normalized residuals are reported so that a relative convergence is reported.

For continuity, the unscaled residual is the sum of mass creation. The reported residual by Fluent is a globally scaled residual, the unscaled residual is normalized by the maximum residual of the first five iterations.

If the max of the continuity residuals during the first five iterations is very small, the scaled residual for continuity will have trouble reaching very small values. In practice, this is achieved by a very good initial guess. Actually, a perfect initialization of a flow with the actual solution will yield a scaled continuity residual of 1. My simulation is perfect, but if I were to use the reduction in residual as my convergence criteria, my simulation would never converge! Even if it did converge, it would converge to the wrong solution, since I know the scaled residual must be 1!

On the other hand, the other scaled (momentum) residuals are calculated differently. The unscaled residuals are calculated the same way but they are normalized by the sum of their convective "speed" of each cell. This actually is not too bad to judge convergence, since if the residual is small then the imbalance is small. But being a large sum, it still does not give too much information about how well converged the solution is at each individual cells. With this method, individual cells can locally have very high imbalances across the cell (even non-physical) without affecting the reported residual. Hence, although the imbalance is small overall, there is potential for individual cells to just be plain wrong. If this individual cell happens to be a critical cell, it can lead to very poor results.

For normalized residuals, the unscaled residuals are normalized by the maximum residual after the first 5 iterations similar to the scaling done by continuity.

These are the defaults for Fluent pressure based solver. There are also normalized residuals and the definition of residual is different for the density based solver but still similar. Also different programs use different methods, but similar discussion applies to them also.

So far I have discussed global scaling, there is also local scaling whereby the unscaled residuals are normalized by the different in max and min residual of the current iteration.

Notice I have not mentioned anything about solution convergence. In other words, there is no guarantee that any of my solution values have converged. Large oscillations in the solution are possible even while the residuals are small and decreasing. Also note that, the continuity equation is normalized differently. This is also the reason for the well-observed result that, typically, the last residual to converge is the continuity equation / the continuity equation typically has the highest residuals.

Last edited by LuckyTran; March 27, 2012 at 11:15.

 March 26, 2012, 04:33 #9 Senior Member   Nick Join Date: Nov 2010 Posts: 126 Rep Power: 8 OK. I'm beginning to understand what you mean...so in my case a look at the unscaled residuals would be a better way to judge convergence in addition to the lift coefficient since the initial guess may mess it up for the scaled residuals ..is this correct? also what i meant before was the flow is supposed to reach a steady state where lift doesnt change (according to experimental results) so perhaps I should only run it in steady mode and monitor the unscaled residuals and lift coefficient

March 26, 2012, 04:49
#10
Senior Member

Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 1,880
Rep Power: 26
Quote:
 Originally Posted by Nick R OK. I'm beginning to understand what you mean...so in my case a look at the unscaled residuals would be a better way to judge convergence in addition to the lift coefficient since the initial guess may mess it up for the scaled residuals ..is this correct? also what i meant before was the flow is supposed to reach a steady state where lift doesnt change (according to experimental results) so perhaps I should only run it in steady mode and monitor the unscaled residuals and lift coefficient
just monitor lift and drag coefficient. monitoring residuals does not give any better data.

flow will probably never reach steady state in an unsteady simulation! unless you magically hit a stable steady solution, probably only possible for a handful of laminar and trivial flow cases. Also, experimentally the flow is time-varying in nature. The experimental results are just the average of time-varying quantities to result in a time-averaged and steady state answer. Flows are inherently time-varying, both numerically and experimentally, and physically!

for unsteady simulations, you need to turn on data sampling for time variables. this will take running averages of all your quantities and you can then see averaged velocity, pressure, etc. I am not sure if averaged lift and drag are included, I doubt it.

that is all the unsteady talk. but I don't see any reason why the problem cannot converge in a steady simulation unless you are specifically trying to extract the unsteady quantities. Running an unsteady simulation on what is supposed to be a steady state simulation will only mask the inherent instabilities, it is a crude workaround and does not actually solve the problem of convergence.

 March 28, 2012, 22:08 #12 Senior Member   Nick Join Date: Nov 2010 Posts: 126 Rep Power: 8 Thank you both for your comprehensive responses. Given your explanation my solution has converged at 4degrees. Last edited by Nick R; March 29, 2012 at 01:31.

 March 29, 2012, 12:04 #13 Senior Member   Join Date: Mar 2011 Location: Germany Posts: 431 Rep Power: 13 No problem at all, you are more than welcome at any time and i also want a little help from you regarding the structured hexa meshing on the wing, i know its not the relevant topic to be discussed here but try to tolerate me. I totally know nothing about hexa meshing and blocking strategy in ICEM and looking forward for some help, kindly help me by any means may be through a self made tutorial. thanks in advance Regards

 March 29, 2012, 20:06 #14 Senior Member   Nick Join Date: Nov 2010 Posts: 126 Rep Power: 8 The best tutorial I've come across for aerodynamic purposes is on youtube. Just run a search on ICEM airfoil and you'll find Simon's three-part tutorial. It'll teach you how to do 2d hex meshing around a foil. Anything more complex such as a 3d geometry would be variations on the 2d version.

 March 31, 2012, 15:32 #15 Senior Member   Join Date: Mar 2011 Location: Germany Posts: 431 Rep Power: 13 I know how to do 2d hexa meshing on the aerofoil, i learned it fron that Simon's tutorials but i know nothing about the 3d hexa meshing, kindly help me in that

 April 1, 2012, 06:00 #16 Senior Member   Nick Join Date: Nov 2010 Posts: 126 Rep Power: 8 if you know how to do a 2d hex mesh around a foil, the 3d scheme will be easy for you. Instead of creating a 2d block begin by creating a 3d block and associate edges to curves accordingly in three dimensions. there is also another tutorial from ICEM where a wing is meshed in 3d

June 3, 2012, 07:14
residual oscillation
#17
Member

sagila
Join Date: May 2012
Posts: 39
Rep Power: 7
Quote:
 Originally Posted by LuckyTran .... Large oscillations in the solution are possible even while the residuals are small and decreasing. ....
Hi Lucky Tran,

Could you please explain more about this? Because, in my case, monitor of solution value has gotten constant and the net flux imbalance of total heat transfer rate is also ok. But, the residuals are oscillating although they are decreasing.

Thanks.
Best,
Sagila

 September 5, 2013, 00:16 #18 New Member   Jalay Shukla Join Date: Aug 2013 Posts: 11 Rep Power: 6 Hi all, In my case i am making simulation on gas flow in room. it has velocity inlet and pressure outlet as BC , but i have a question my k-e doesnt converge but all others converge (gas,velocity,energy) So what is the problem? And please give me some info about solution initialization ..i am quite unaware of it

 February 4, 2015, 12:21 #19 New Member   sajad Join Date: Aug 2010 Posts: 17 Rep Power: 8 Hi all, I'm modelling a a 3d rotor blade using density based transient solver in Ansys Fluent. I am using pressure far field BC and Spalart-Allmaras turbulence model. I have problem with residuals. in fact the continuity residual doesn't go lower than 1e-2 and other residuals doesn't go lower than 1e-3 and remain constant or oscillate around. I have tried first order and second order implicit formulation, single and double precision. i also tried varying the initial condition since i thought it might be responsible for the problem. none of these didn't help. I have noticed that when i increase yplus on the blade surface (by increasing first cell width) i get better residuals. for example if i increase yplus to 600, continuity residual decrease to 1e-3. I thought it might be because of aspect ratio but decreasing aspect ratio didn't solve problem. so what should i do? should i just ignore residual values and monitor lift and drag? what is causing this problem? I'd really appreciate any suggestions to resolve this issue. (I should apologize for my bad english. i hope this doesn't stop you helping me)

 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 xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33 colopolo CFX 13 October 4, 2011 22:03 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03 msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58 liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07