|
[Sponsors] |
March 27, 2014, 18:16 |
Is there a problem in the Euler solver?
|
#1 |
Member
Join Date: Sep 2013
Posts: 43
Rep Power: 13 |
I am doing some validation cases in Euler and Navier-Stokes with the naca0012 airfoil.
The Navier-Stokes solver seems to give good results compared to the ones provided by NASA (http://turbmodels.larc.nasa.gov/naca0012_val_sa.html) In my opinion, at low Mach number, the Euler solver should give zero drag (in 2D). Unfortunately, it is not what I get with the naca0012 airfoil. For example, at zero angle of attack, I get Cd = 0.0004 (even if it is not huge, it is not negligible). When I zoom on the flow around the airfoil, there is a loss of speed near the airfoil, and so a loss of total pressure, which could explain the drag I get. My question is: is it a problem of numerical scheme, or a problem of mesh (it is the same mesh as the one used for Navier-Stokes, so with very thin cells near the boundary), or a problem in the Euler solver? I put the some files in attachment to illustrate what I said. Thank you for your answers Laurent |
|
March 27, 2014, 21:17 |
|
#2 |
Senior Member
Heather Kline
Join Date: Jun 2013
Posts: 309
Rep Power: 14 |
Yes, for a 2-D body in inviscid flow there should be zero drag (d'Alembert's paradox).
The small finite drag is likely due to numerical dissipation, which is the same mechanism that enforces the Kutta condition here. Most if not all CFD codes have this effect to some extent. Mesh refinement can reduce the magnitude, as would reduced numerical error within the scheme. |
|
March 28, 2014, 05:31 |
|
#3 |
Member
Join Date: Sep 2013
Posts: 43
Rep Power: 13 |
Thank you for your answer, but what seems strange to me is that with the same mesh and the same scheme but in Navier-Stokes, results are much better (I get exactly the drag given by NASA on http://turbmodels.larc.nasa.gov/naca0012_val_sa.html so a numerical error on Cd below 0.0001)
And the loss of velocity near the wall should not happen: it looks like a boundary layer with a very low friction coefficient... That is why I am wondering if the Euler solver does not have a problem... Laurent |
|
March 28, 2014, 14:49 |
|
#4 | |
Super Moderator
Francisco Palacios
Join Date: Jan 2013
Location: Long Beach, CA
Posts: 404
Rep Power: 15 |
Quote:
My answer is related with Heather's reply. The key is the artificial dissipation introduced by the numerical scheme. Navier-Stokes grids have been created to simulate problems with boundary layers, in other words they have an a-priori grid adaptation that has beneficial effect when you run a Navier-Stokes problems. But these grids are a very bad choice when you try to run an Euler problem (completely different behavior on the wall: hyperbolic vs. elliptic). So… with high-stretching grids, the convective scheme is not adding the optimal dissipation (in the normal direction the wall) and you detect that strange effect. Obviously, if you solve a NS problem, the viscous effects are dominant close to the wall and the solution is correct. Conclusion, the a-priori adaptation should be done based on the equation that you are solving. Otherwise it will deteriorate the convergence and the solution quality because you are adapting some features that are not present in your continuous model. Thanks for using SU2, Best, Francisco |
||
March 28, 2014, 16:48 |
|
#5 |
Member
Join Date: Sep 2013
Posts: 43
Rep Power: 13 |
Thank you very much for this answer that solves my problem. So I am going to try better meshes for Euler cases.
Thank you for your great work! Best regards, Laurent |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Abaqus/cfd steady-state solver problem | skuznet | Main CFD Forum | 2 | June 11, 2016 16:46 |
Problem with compiling new solver | palazi88 | OpenFOAM Programming & Development | 2 | December 24, 2013 20:52 |
Problem with implicit unsteady solver | CCMuser | STAR-CCM+ | 2 | March 3, 2010 12:20 |
Problem Interface Solid Fluid with wall velocity Solver v12 | hills1 | CFX | 2 | October 12, 2009 06:36 |
compressible two phase flow in CFX4.4 | youngan | CFX | 0 | July 2, 2003 00:32 |