
[Sponsors] 
May 9, 2010, 01:48 
BezierBernstein approximation for airfoil shape

#1 
Senior Member
Join Date: Nov 2009
Posts: 411
Rep Power: 13 
Hello,
I'm interested if it is possible to approximate an airfoil shape with 2 or more Bezier curves of a moderate order (say 1012). What I want to know is if it is possible that given a certain number of points (say 360 points) of an airfoil to fit through these a 12 order Bezier curve. I've read an introductory chapter in Bezier curves but I didn't found any clear explanations of how you can approximate a given data set with a Bezier curve. The classical approach seems to be to use a Bezier curve of order N1 where N is the number of your data points. I would appreciate any insight in the subject. Thanks, Do 

May 9, 2010, 04:20 
what is Shadows?

#2 
New Member
Uday
Join Date: Mar 2010
Posts: 3
Rep Power: 9 
Hi all of you, I want to do meshing at the inetrface of two surfaces using gambit meshing. But at the interface i find error of Shadows when the geometry is exported to Fluent.
So, what is the meaning of Shadows? how to tackel the problems with these Shadows? please reply..... 

May 9, 2010, 07:38 

#3 
Super Moderator

I am interested in the same question at present. I have been reading up on bsplines and nurbs. But the books never tell you how to get a curve to fit given data.
I guess you just have to do some least squares fit. Thats what I am going to try. A new and better way to represent airfoil shapes has been proposed by Brenda Kulfan. Search for her name in aiaa papers and journals. This is especially interesting if you are going to do shape optimization. 

May 9, 2010, 08:51 

#4 
Senior Member
Join Date: Nov 2009
Posts: 411
Rep Power: 13 
Thanks Praveen,
least square approximation was what I have had in mind, basically you will determine by least squares the control points. I've already implemented similar codes for Chebyshev polynomials and other orthogonal polynomials, so (in a first approx.) I will just replace the Chebyshev function with one that calculates a Bernestein polynomial! Sure in order to have a Bezier curve you need to interpolate exactly the first and the last point of your data set, so my code isn't exactly right, but it can be easily modified. Shape optimization is what I'm trying to do. Do 

May 9, 2010, 08:53 

#5 
Senior Member
Join Date: Nov 2009
Posts: 411
Rep Power: 13 
It works! except the for the first and last points of course.
Do 

May 9, 2010, 08:59 

#6 
Super Moderator

I am assuming you are fitting two bezier curves, one for the upper and another for the lower curve. For each curve you can fix the first and last control points to be equal to the first and last point coordinates. Then determine the remaining control points by least squares. You can even put in information about the tangents at the LE/TE points into the bezier curve.
If you are using a high degree bezier curve, you should not directly evaluate the bernstein polynomials as that leads to loss of accuracy. A better way is to use the de Casteljau algorithm. In that case, you must do the least squares fitting by using an optimization method. 

May 9, 2010, 09:11 

#7 
Senior Member
Join Date: Nov 2009
Posts: 411
Rep Power: 13 
Yes, I use 2 Bezier curves of 12 order. I will check the accuracy of my Bernstein polynomials, but I suppose for the order 12 I will have no particular accuracy loss. Anyway I can check the error switching to quad precision versus double precision as I usually use in my code.
In case I will notice some degradation in accuracy I will use 4 Bezier curves of order 6, 2 for the upper side and 2 for the lower side of my airfoil. Do 

May 9, 2010, 13:46 

#8 
Senior Member
Join Date: Nov 2009
Posts: 411
Rep Power: 13 
@Praveen
Just to be sure I'm not in the wrong direction, when you do shape optimization of an airfoil using a Bezier curve, the design variables are the control points of the Bezier curve. Correct ? Thanks, Do 

May 10, 2010, 08:10 

#10 
Senior Member
Join Date: Nov 2009
Posts: 411
Rep Power: 13 
Thanks.
Do 

April 23, 2014, 06:21 
Airfoil representation by bezier curves

#11 
New Member
shikhar
Join Date: Apr 2010
Posts: 22
Rep Power: 9 
Hi
I am trying to represent airfoil by using 2 bezier curves of order 12. what i understand is i have to find control points of bezier curve by using the coordinates of the airfoil and doing least squares approximation. My question is in typical least squares approx. we minimize vertical offset (i.e. we minimize {y  y(t)}^2 ). So if i do this i will find the y coordinates of my control points. Now, how do i find the x coordinates of my control points??? Can anyone please help me. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Calculation of the Governing Equations  Mihail  CFX  7  September 7, 2014 06:27 