|
[Sponsors] |
Wing in compressible flow convergence problem |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Member
Ferruccio Rossi
Join Date: Jun 2016
Location: Melbourne, FL USA
Posts: 91
Rep Power: 10 ![]() |
Hi,
I am simulating compressible flow over a 737 wing approximately 20 m long (see pictures) at cruise conditions: Pressure Farfield boundary conditions Mach = 0.745 Angle of attack = 2.5 deg Absolute pressure = 24998 Pa Temperature = 220 k Turbulence = Spalart-Allmaras Coupled scheme Second order for all the spatial discretizations Pseudo Transient My solution only converges to 3e-1. It doesn't even get to my target broad convergence value, which I set to 1e-3. Moreover, my results seem to be off. Although the 737 wing is real scale (and thus large), how can I get such high values for the coefficients? Lift coeff is 29! Drag coeff = 2 Lift Coeff = 29 I am using a mesh which I believe is fine (see pictures). Total number of elements = 6,203,332 Total number of nodes = 2,597,576 Wing face elements size = 5e-2 m Inflation layers = 36 layers, 5e-4 m first layer thickness, 1.05 growth rate So why my simulation is not converging? Do I need a finer mesh or finer inflation layers? Or something else could be causing the problem? If you look at my residuals plot, you see that the values are stable although not converged to my target. Are the results still reliable if they are stable but not converged? I am very confused. Anything helps. Thank you all! |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
Ferruccio Rossi
Join Date: Jun 2016
Location: Melbourne, FL USA
Posts: 91
Rep Power: 10 ![]() |
here are my plots
|
|
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,728
Rep Power: 66 ![]() ![]() ![]() |
For coefficient values: Make sure you have set all the reference values.
For convergence: Your solution looks fine. Only the residual for the turbulent viscosity is high. Run the mesh check utility. A fine mesh is not necessarily a good mesh. It looks like there is a sharp transition between your inflation layers and core mesh, so that you end up with a very large volume change. This can limit your convergence. Also check if the skewness is high in these regions. This could easily be fixed by using a larger growth rate. e.g. 1.1 instead of 1.05. Also plot a contour of velocities and see if your prism layers cover the entire boundary layer, they should. But when generating fine meshes it's easy to make the prism layer thickness too thin. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Member
Ferruccio Rossi
Join Date: Jun 2016
Location: Melbourne, FL USA
Posts: 91
Rep Power: 10 ![]() |
Thank you for your answer. I made a more balanced mesh with a better transition between inflation layers and core mesh. I will see the results after the simulation is done.
For coeff values: For what concerns the reference values, I am a little confused on how to use them. I used the default compute from farfield mesh, and this is what I got (see picture). The values for pressure, velocity, temperature, etc all makes sense, but the values of area and length are both set to 1. Is this the reason my my CL is off? To what should I set this values? For convergence: I thought the unconverged residual was continuity (white line). Isn't turbulent viscosity already converged (purple line)? Thank you again for your patience. |
|
![]() |
![]() |
![]() |
![]() |
#5 | |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,728
Rep Power: 66 ![]() ![]() ![]() |
Quote:
The continuity residual is not scaled the same way as other variables (it's dependent on the initial guess). If you have a good initial guess, you'll have high continuity residuals. I would ignore the continuity residual. If you're okay with the turbulent viscosity then you're already done. |
||
![]() |
![]() |
![]() |
![]() |
#6 | |
Member
Ferruccio Rossi
Join Date: Jun 2016
Location: Melbourne, FL USA
Posts: 91
Rep Power: 10 ![]() |
Thank you so much!!! This is the best answer I ever got on this topic and it finally clarified a lot!!!
So, I define lift coeff with the standard CL formula: CL = LIFT / (0.5 * density * velocity^2 * WingArea) where WingArea is the projected area of my wing = 52 m^2 Based on this definition of CL, my CL is calculated using wing area instead of wing length. So for the reference values I would input 52 m^2 for area (since that is my wing area) and 0 for length (since I don't need it for my CL). Is this reasoning correct? For the residuals: When you say Quote:
To test my hypothesis, I ran a simulation using a density based solver, but this time continuity converged, while z-velocity, x-velocity and turbulence didn't converge (see picture). 1. Should I ignore continuity in both solver types when it comes to assess convergence? 2. And what are the most important plots that I must look at to consider my simulation reliably converged? 3. What solver (pressure or density based) would you recommend for solving high speed subsonic flows (Mach = 0.745)? Sorry for the many questions but you are very clear and I am understanding a lot from your explanations. Thanks! |
||
![]() |
![]() |
![]() |
![]() |
#7 | |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,728
Rep Power: 66 ![]() ![]() ![]() |
Quote:
If you initialize your flow with the exact solution, the continuity residuals get stuck at 1. Or if you reset your residuals right now, you'll see something similar. Your plot of solution vs iteration is the best indicator of convergence. You can use the pressure based solver, and I recommend it. It's much easier to use and pretty robust. The density based solver really saves you when you have shockwaves, but the pressure based solver can also handle those problems as long as you are careful. Start with the pressure based solver. If you want to switch to the density based solver after, you can. |
||
![]() |
![]() |
![]() |
![]() |
#8 | ||
Member
Ferruccio Rossi
Join Date: Jun 2016
Location: Melbourne, FL USA
Posts: 91
Rep Power: 10 ![]() |
Ok, got it.
Last thing ![]() Quote:
and when you say Quote:
And what residual value is considered solid convergence (1e-4,1e-5...)? |
|||
![]() |
![]() |
![]() |
![]() |
#9 | |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,728
Rep Power: 66 ![]() ![]() ![]() |
People always ask for black and white answers and I am afraid I must never make them happy. I always answer yes and no.
It's a big mistake in the way we do education because people are always used to a problem with only one method and they know the solution and there's exactly 1 solution. In real life, you have a problem, no method, and no known solution. Quote:
Residuals tell you the imbalances leftover when you say stuff = 0, and the stuff is not zero. Nothing is exact in numerics, so the user is always the one that defines the tolerance within which something is satisfied. You get to choose the residual levels. Residuals don't tell actually anything about convergence (because they don't tell you what the solution is). But you would think (intuitively) that if the solution is satisfying my governing equations more and more, that it's converging to a solution. This is a nice intuition to have, until you one day run into a problem that's bi-stable (a problem with two solutions). The problem is the continuity residual is scaled in a way that's not very meaningful. If you have a horrible initial guess, the continuity residual drops rapidly. If you have a beautiful initial guess, it's stuck. It's a rather dumb residual unless you intentionally make a very bad initial guess. So if I make a an awesome initial guess, I can get to the answer in 200 iterations with a continuity residual of 1E-01; whereas you make can a horrible initial guess, arrive at the same answer in 20,000 iterations, and end up with a continuity residual of 1E-10. We both have the same solution but we have different continuity residuals! And if you look at how much the continuity residual has dropped, you would think that the person who makes the worst initial guess and wastes the most time computing has the "more converged" solution. Who would you rather hire? For what purpose are you doing CFD? Are you doing CFD for the sake of watching residuals drop or are you doing CFD to get some sort of result? Some sort of sought-after-parameter? Whatever that reason is, monitor that variable because that's what you care about most. You want that result to converge. But if that is too painful, then the next best things to check the basic variables, the variables that Fluent is actually solving for: x-velocity, y-velocity, z-velocity, pressure, temperature, turbulence, etc. But then you have to get creative and decide where to monitor these things. Again, you should decide. It needs to make sense. Obviously don't monitor the solution on your neighbor's machine. Monitor something in your solution. If you don't know where to look, then look anywhere (or everywhere). My best advice is to just let Fluent iterate forever. Disable all the checks for convergence and just let it iterate forever. Go away from the computer. Your first few simulations should be done this way. Come back when you're bored (at least a whole day later, over a weekend is a good time to let it do this). Hopefully your simulation did not diverge and blow up. If it did, fix it. If it didn't, now you can sit down and view the results. Look at what happened to your monitors and residuals. And now decide where do you think the solution was converged? |
||
![]() |
![]() |
![]() |
![]() |
#10 |
Member
Ferruccio Rossi
Join Date: Jun 2016
Location: Melbourne, FL USA
Posts: 91
Rep Power: 10 ![]() |
Thank you for your time. Your posts have been really useful.
Much appreciated! |
|
![]() |
![]() |
![]() |
Tags |
convergence, drag, lift, residuals, wing |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Two Mass flow inlet bc convergence problem | nabidinhomessi | Main CFD Forum | 5 | December 14, 2015 07:11 |
Convergence Problem in Axisymmetric Periodic Flow | atheresia | FLUENT | 3 | February 10, 2014 03:00 |
momentum & x-vel convergence, and reverse flow problem | cinwendy | FLUENT | 0 | April 17, 2013 02:19 |
Convergence Problem with two phase flow | challenger85 | CFX | 0 | December 29, 2009 08:01 |
convergence problem of three phase flow | nick | CFX | 3 | April 25, 2008 09:00 |