# Source-Doublets panel methods MatLAB

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

May 4, 2022, 12:47
Source-Doublets panel methods MatLAB
#1
New Member

Sebastiano Vianello
Join Date: Mar 2021
Posts: 1
Rep Power: 0
Good evening everyone.
I'm developing a 3D panel method in MatLAB for aerodynamic analysis of a generic wing. This will be later on coupled with a structural solver to conduct some aeroelasic analysis.
The problem setup follows the method proposed by Katz and Plotkin in the well known 'Low Speed Aerodynamics'.
The geometry is discretized using a constant source-doublets distribution on each panel: the former provides the thickness contribution, the latter the vorticity contribution.
The analyzed airfoil is a NACA0006 straight wing, with chord = 1m and span = 10m. There are 10 spanwise sections and 100 panels per section.
For , the pressure distribution (Cp), is close to the one computed using XFOIL (albeit the Cp variation is not as abrupt as though close to the leading edge).
When imposing an angle to the profile( in this case , however, the situation is much different, with the pressure distribution being much different.
You can find attached the doublets and sources intensities distributed on the airfoil, which seems reasonable to me.
I've tried to look at other implementations of the method but from that I couldn't find what was wrong with my code.
Did anyone had a similar experience? Or is willing to help me with the solution of the problem?
Thank you everyone in advance for your help.
Attached Images
 0006_0deg_Cp.png (13.0 KB, 16 views) 0006_0deg_doublet.png (40.1 KB, 20 views) 0006_5deg_Cp.png (24.9 KB, 17 views) 0006_5deg_doublet.png (37.1 KB, 20 views)

 May 13, 2022, 12:50 #2 Member   james nathman Join Date: May 2011 Posts: 62 Rep Power: 14 The trailing-edge pressure at 5 degrees indicates your method has no wake. What are the lift coefficients at the centerline of the wing for 0 & 5 degrees?

 May 17, 2022, 12:43 #3 Senior Member   Will Kernkamp Join Date: Jun 2014 Posts: 307 Rep Power: 12 Blackjack is correct. You need to create wake panels with a constant doublet strength equal to the sum of the doublet strength of the trailing edge panels of upper and lower surface. In that way, you ensure that there is no vortex at the trailing edge. This is the so called Kutta condition at the trailing edge. It causes the flow to smoothly leave the trailing edge. This is necessary for lift to occur. The wake panels should be in the direction of the flow, but it doesn't have to be very precise for a decent result. A good start would be to extend the edges from the trailing edge straight in the flow direction. The length should be a several times the span of the wing to ensure that the vortex at the end of your wake panels does not influence the flow at the wing.

 Tags doublets, matlab, panel, source

 Thread Tools Search this Thread Search this Thread: Advanced Search Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Seroga OpenFOAM Community Contributions 9 June 12, 2015 17:18 [swak4Foam] swak4foam building problem GGerber OpenFOAM Community Contributions 54 April 24, 2015 16:02 brbbhatti OpenFOAM Programming & Development 2 July 7, 2014 11:32 [swak4Foam] Swak4FOAM 0.2.3 / OF2.2.x installation error FerdiFuchs OpenFOAM Community Contributions 27 April 16, 2014 15:14 [swak4Foam] Error bulding swak4Foam sfigato OpenFOAM Community Contributions 18 August 22, 2013 12:41

All times are GMT -4. The time now is 05:54.

 Contact Us - CFD Online - Privacy Statement - Top