|
[Sponsors] |
December 13, 2021, 15:13 |
Grid refinement or Higher order scheme
|
#1 |
Member
Yusuf Elbadry
Join Date: Sep 2018
Posts: 65
Rep Power: 8 |
Hello,
I hope you are all fine, I am solving the heaving airfoil problem while using the Eulerian grid (Level Set Method). The 2D Incompressible N-S equations are solved using GLS FEM. I am simulating the flow at Re=800, Angle of Attack =20 Deg. I am using bilinear shape function for geometry and solution and I am using an explicit scheme. I have a problem with the accuracy of the result, I am comparing my result (Cl, CD) with the literature and the error is not small, and I was thinking to use a second-order shape function for both the geometry and solution to get a better result and better shape for the airfoil instead of stairs. Would this be better or should I use a more refined grid with a very small time step? |
|
December 13, 2021, 16:49 |
|
#2 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73 |
Quote:
Before introducing a higher order scheme, you must verify that your code works well for via via refined grids. If your solution does not improve under grid refinement, your problem is not the low accuracy order but some different issues. |
||
December 13, 2021, 17:04 |
|
#3 | |
Member
Yusuf Elbadry
Join Date: Sep 2018
Posts: 65
Rep Power: 8 |
Quote:
Thank you for your reply. |
||
December 13, 2021, 17:19 |
|
#4 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73 |
Quote:
For an arbitrary fixed grid of size h, you cannot be sure that the solution obtained with a second order scheme is better than a first order scheme. You need to have a sufficient grid resolution to get a real improving. That is a classical topic in CFD. |
||
December 13, 2021, 17:58 |
|
#5 | |
Member
Yusuf Elbadry
Join Date: Sep 2018
Posts: 65
Rep Power: 8 |
Quote:
For example, if I have an airfoil with a chord of 1 (nondimensional length), I will need elements with a side length (in y-dir) of 3e-3 or less to have a smooth airfoil shape, If the airfoil will move in the y-direction with an amplitude of 0.5c (from -0.5c to 0.5c), then I need at least 333 divisions in this direction, and let us say I will have one-third of this number in the x-direction, this may result in around 350,000 elements near the airfoil where it moves, not mentioning the rest of the domain. This may result in a grid of more than 600,000 elements. That is why I thought of using high order shape function, to reduce the number of elements and nodes. I really appreciate your help and guidance and thank you for your patience. |
||
December 14, 2021, 05:24 |
|
#6 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73 |
Quote:
I don't understand, the shape functions defines a certain functional relation for the unknown variables, what do you mean now for the highlighted statement? The better description of the airfoil geometry is a different issue. However, grid size and accuracy order defines different topics and different limits.The former fix the range of resolved wavenumber in the solution, the latter is the theoretical and asymptotical velocity rate of the error decreasing. You could have also a spectral accuracy but if you have fixed the grid size you cannot resolve nothing greater than the Nyquist frequency. |
||
December 14, 2021, 14:59 |
|
#7 | |
Member
Yusuf Elbadry
Join Date: Sep 2018
Posts: 65
Rep Power: 8 |
Quote:
Please correct me if I am wrong or if there is something I am missing. |
||
December 14, 2021, 16:35 |
|
#8 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73 |
Quote:
But this is a further topic. Think about a BEM formulation (panel method), you can describe the airfoilf by using straight segment or curved step. Then you prescribe a functional relation for the distribution of singularity. In you case you want to use elements that are not linear triangles but curved triangles and then use high order shape functions. Again, on the same and arbitrary mesh size h you cannot say a-priori that the error is for sure lesser for the quadratic shape function. You need to assess that by a grid refinement study. Only at a certain value h you get a monotonically expected rate for the error. |
||
December 15, 2021, 15:02 |
|
#9 | |
Member
Yusuf Elbadry
Join Date: Sep 2018
Posts: 65
Rep Power: 8 |
Quote:
I really appreciate your patience and your reply. |
||
December 15, 2021, 15:59 |
|
#10 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73 |
Quote:
|
||
December 16, 2021, 14:26 |
|
#11 | |
Member
Yusuf Elbadry
Join Date: Sep 2018
Posts: 65
Rep Power: 8 |
Quote:
I will try using a higher-order shape function for both and solution and geometry. Thank you |
||
December 16, 2021, 14:33 |
|
#12 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73 |
Quote:
You could search for suitable isoparametric elements. On the other hand, the airfoil is better described by a triangular tesselation. |
||
December 16, 2021, 15:15 |
|
#13 |
Member
Yusuf Elbadry
Join Date: Sep 2018
Posts: 65
Rep Power: 8 |
I know that there is a Lagrangian Element(9 nodes in case of bi-quadratic shape function) and Serendipity Element(8 nodes in case of bi-quadratic shape function), can I find a different element than these two types of the element?
|
|
December 16, 2021, 15:25 |
|
#14 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73 |
Quote:
Lagrangian simplex, that is triangular elements (or tetrahedron in 3D), having 3, 6, 10, ... nodes |
||
Tags |
fem, high order scheme., lsm |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Third order MUSCL Scheme | 94bb494nd41f | FLUENT | 4 | January 25, 2024 13:23 |
First Order to Higher Order Blending Factor | NormalVector | FLUENT | 4 | November 13, 2023 08:06 |
Implementation of higher order FD scheme near boundries | pawank | Main CFD Forum | 16 | September 23, 2016 09:52 |
Higher order downwind scheme | jelmer | OpenFOAM Running, Solving & CFD | 4 | August 9, 2006 07:43 |
Combustion Convergence problems | Art Stretton | Phoenics | 5 | April 2, 2002 06:59 |