|
[Sponsors] |
[Airfoil] Bad results for LES in comparison to RANS |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 17, 2020, 05:35 |
|
#21 |
New Member
Join Date: Jul 2020
Posts: 25
Rep Power: 6 |
We already said:
- your BCs are correct - to extend spanwise direction - that 2D LES is meaningless because turbulence is 3D so don't expect good results Then you started talking about 2D-3D and other stuff, if this is not relevant to the first question (Why is 2D LES in Fluent not giving good results?) don't talk about it because it's too confusing. Your english is not the problem. |
|
September 17, 2020, 05:41 |
|
#22 | |
Senior Member
Join Date: Jan 2018
Posts: 121
Rep Power: 8 |
Quote:
Best regards |
||
September 17, 2020, 05:49 |
|
#23 |
New Member
Join Date: Jul 2020
Posts: 25
Rep Power: 6 |
||
September 17, 2020, 08:05 |
|
#24 |
Super Moderator
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,426
Rep Power: 49 |
Let me state the obvious, just in case: when expanding the 2D domain to the third dimension, this third dimension needs to be adequately resolved by the mesh for LES to do its magic. Just stretching out 1-5 cells won't work.
|
|
September 17, 2020, 10:38 |
|
#25 | |
Senior Member
Join Date: Jan 2018
Posts: 121
Rep Power: 8 |
Quote:
But today a guy checked the velocity contours and told me that the high drag coefficient is probably due to the premature transition and separation of the flow (you can see that the flow separates very soon close the leading edge) and this is due to the use of central difference scheme upstream of the transition region which causes numerical trip of the boundary layer. Is that right? Kind regards. |
||
September 17, 2020, 10:43 |
|
#26 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 |
Quote:
I don't see the meaning of the statement of this guy. LES requires accurate discretizations and the second order central discretization was the hystorical scheme used. Again, you miss the key point. The 2D or quasi 2D approximation forces the vorticity to have the direction normal to the plane and no action of the stretching appears. Without stretching the energy has a different way of distributing along the length scale. For example, in a 2D flow you can see small vortices merging in larger. Try to extend the third dimension to a reasonable lenght and use a step size dz+=O(20-30). Only then we can discuss. |
||
September 17, 2020, 11:50 |
|
#27 | |||
Senior Member
Join Date: Jan 2018
Posts: 121
Rep Power: 8 |
Quote:
Thank you dear FMDenaro. My friend sent me a PDF which has more details about the LES of this airfoil. You can download it here: http://www.tfd.chalmers.se/~lada/pos.../simon_lic.pdf He said that his comments on the central difference scheme is from this paper you can check it in the following sections that I quoted some of them: 4.4 Discretisation Schemes of the Convective Terms (Page: 17) Quote:
4.8.2 Transition (Page: 27) Quote:
I kbow it's not a pure 3D simulation and the papers claimed this. It's a quasi-3D simulation Now I'm running the 3D LES with 33 nodes in spanwise direction and the length of 3% of the chord length and waiting for the results. Although you told that it needs a wider spanwise extent, but I have to choose this length because I need to compare the results with a reference simulation and they claimed that for 33 nodes, the 3% extent had better results that 8%. Also my computational power is very limited and I can't have more than 4 million cells. The 3% length is adapted from the paper. You can check the following sections: 4.5.1 The resolution and extension in the spanwise direction (Page: 20) 4.8.3 Resolution in the streamwise and spanwise directions (Page: 29) Best Regards |
||||
September 17, 2020, 12:39 |
|
#28 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 |
Quote:
LES is based on a basic assumption, the SGS model must account for the required energy transfer, that is both the inertial and disispative unresolved parts. The oscillation caused by central discretization is a different issue of numerical character. But in no way that implies that a bounded scheme must be used in LES. This is somehow a conflict of two different actions the limiter acting for numerical reason and the SGS model for the unresolved components. Usually the result is a very excess of dissipation. The reason why the bounded scheme is present in Fluent is only one: to satisfy the users to get some solution, no matter of the accurate physical meaning. But if your LES setup is correct and the grid is suitable, the unbounded scheme must work. |
||
September 17, 2020, 12:50 |
|
#29 | |
Senior Member
Join Date: Jan 2018
Posts: 121
Rep Power: 8 |
Quote:
I'm using unbounded central difference scheme now and I'm agree with you that Fluent just wants to give the user a solution while the LES needs accuracy. The funny thing is that by enabling LES and Smagorinsky subgrid scale model, it says that it must be used with bounded scheme! After the new simulation results are ready, I'll post the results. Best regards. |
||
September 17, 2020, 13:15 |
|
#30 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 |
Quote:
You can force Fluent to use unbounded central scheme for LES. However, I suggest also to use the dynamic Smagorinsky model. |
||
September 18, 2020, 05:23 |
|
#31 | |
Senior Member
Join Date: Jan 2018
Posts: 121
Rep Power: 8 |
Quote:
Thank you dear FMDenaro! I'm using dynamic Samgorinsky model. I just noticed something weird! When I use Convective Outlet B.C for velocity + homogeneous Neumann B.C for pressure at the outlet boundary, the outlet mass flow is lower than the inlet mass flow. (The difference is about 0.5 kg/s and the total mass flow is about 50 kg/s) But when I use the same B.C for the velocity, but Dirichlet B.C for pressure, the inlet and outlet mass flows are equal. It seems that there's a problem with the Neumann B.C for the pressure. is that OK? Best Regards |
||
September 19, 2020, 11:33 |
|
#32 |
New Member
|
Hi Moreza,
I think it is important that you look at time-averaged LES results. Your initial post was a snapshot of the LES, which can't be compared to RANS, even after you run the 3D domain. This will at least partially explain the differences in pressure and mass flows you observe. A 1% fluctuation in mass flow is not surprising, but the time-average fluctuation should be 0%. I didn't see where anyone mentioned, but if you want to predict the transition of laminar to turbulent flow, it will be important to check that the eddy viscosity is very low in the laminar region for your mesh. Dynamic Smagorinsky should accomplish this, but I do not think it is a given. Here is a useful paper on the matter: https://web.stanford.edu/group/ctr/R.../09_sayadi.pdf |
|
September 19, 2020, 12:01 |
|
#33 | |
Senior Member
Join Date: Jan 2018
Posts: 121
Rep Power: 8 |
Quote:
And about the mass flow out, it's strange for me why while using the Drichlet boundary condition, mass flow does not oscillate and it only oscillates while using Nuemann B.C. Best regards. |
||
September 19, 2020, 14:15 |
|
#34 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 |
How could you prescribe a Dirichlet BC at the outlet if the flow is unsteady ...
|
|
September 19, 2020, 14:22 |
|
#35 | ||
Senior Member
Join Date: Jan 2018
Posts: 121
Rep Power: 8 |
Quote:
I know the pressure outlet B.C in Fluent as Dirichlet because in the Docs it is given: Quote:
This is what I understand from the Docs. So if i'm wrong or it's a misunderstanding of the text by me, please correct me. |
|||
September 19, 2020, 14:31 |
|
#36 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 |
Quote:
Be aware, the set of BCs for the compressible subsonic flows is specified by one Dirichlet condition at the outlet but the other required BCs cannot be prescribed as Dirichlet. |
||
September 19, 2020, 14:37 |
|
#37 | |
Senior Member
Join Date: Jan 2018
Posts: 121
Rep Power: 8 |
Quote:
Which boundary condition do you prefer for this problem at the outlet? Is Neumann the best choice? |
||
September 19, 2020, 14:39 |
|
#38 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73 |
Quote:
I am also talking about the BCs for the outlet. For compressible subsonic flows you can fix a condition for the pressure (Dirichlet) but the other variables must be let free to develop from the interior, according to the characteristic theory |
||
September 19, 2020, 14:43 |
|
#39 | |
Senior Member
Join Date: Jan 2018
Posts: 121
Rep Power: 8 |
Quote:
The misunderstanding of your message is due to my poor English. Sorry for that. So, I think I must go with Neumann B.Cs... |
||
September 28, 2020, 13:23 |
|
#40 |
New Member
|
I'm not sure what you mean in your first question. But for the z-direction, if a scale-resolving simulation is done (i.e. LES), the length should be sufficient to capture the largest eddies, and the number of cells should be enough to keep the cell aspect ratio as close to 1 as possible as the filter shape can dampen turbulence in certain directions. Ideally z+ = y+ < 1.
|
|
Tags |
airfoil, les, rans, spalart allmaras model |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
bad heat transfer results with low y+ | me45 | OpenFOAM Running, Solving & CFD | 0 | April 29, 2020 11:56 |
Post Results Comparison | RustyCohle | Main CFD Forum | 0 | April 6, 2020 04:46 |
Comparison of 2D simulation results with 3D | Mohit Gupta | Main CFD Forum | 0 | September 29, 2008 13:04 |
Bad results when compared with Wind Tunnel | Luiz | CFX | 7 | October 27, 2006 09:35 |
RANS results for LES | Li Yang | Main CFD Forum | 5 | June 17, 2002 05:34 |