CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   FLUENT (http://www.cfd-online.com/Forums/fluent/)
-   -   Continuity convergence issue on a 3d wing (http://www.cfd-online.com/Forums/fluent/98918-continuity-convergence-issue-3d-wing.html)

 Nick R March 21, 2012 23:25

Continuity convergence issue on a 3d wing

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

 cfd seeker March 22, 2012 04:50

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

 Nick R March 22, 2012 08:46

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.

 LuckyTran March 24, 2012 01:13

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.

 Nick R March 25, 2012 03:45

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.

 LuckyTran March 25, 2012 03:54

Quote:
 Originally Posted by Nick R (Post 351307) 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.

 robboflea March 25, 2012 05:32

Quote:
 Originally Posted by LuckyTran (Post 351308) 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

 LuckyTran March 25, 2012 14:07

Quote:
 Originally Posted by robboflea (Post 351320) 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.

 Nick R March 26, 2012 04:33

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

 LuckyTran March 26, 2012 04:49

Quote:
 Originally Posted by Nick R (Post 351443) 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.

 cfd seeker March 27, 2012 11:07

kindly tell these and if still you are not satisfied and has any confusion feel free to ask. Regards

 Nick R March 28, 2012 22:08

Thank you both for your comprehensive responses. Given your explanation my solution has converged at 4degrees.

 cfd seeker March 29, 2012 12:04

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

 Nick R March 29, 2012 20:06

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.

 cfd seeker March 31, 2012 15:32

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

 Nick R April 1, 2012 06:00

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

 sagila June 3, 2012 07:14

residual oscillation

Quote:
 Originally Posted by LuckyTran (Post 351372) .... 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

 jalay September 5, 2013 00:16

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