CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

Source-Doublets panel methods MatLAB

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   May 4, 2022, 12:47
Default Source-Doublets panel methods MatLAB
New Member
Sebastiano Vianello
Join Date: Mar 2021
Posts: 1
Rep Power: 0
Wian is on a distinguished road
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 \alpha=0, 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 \alpha = 5, 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
File Type: png 0006_0deg_Cp.png (13.0 KB, 16 views)
File Type: png 0006_0deg_doublet.png (40.1 KB, 20 views)
File Type: png 0006_5deg_Cp.png (24.9 KB, 17 views)
File Type: png 0006_5deg_doublet.png (37.1 KB, 20 views)
Wian is offline   Reply With Quote

Old   May 13, 2022, 12:50
james nathman
Join Date: May 2011
Posts: 62
Rep Power: 14
blackjack is on a distinguished road
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?
blackjack is offline   Reply With Quote

Old   May 17, 2022, 12:43
Senior Member
Will Kernkamp
Join Date: Jun 2014
Posts: 307
Rep Power: 12
wkernkamp is on a distinguished road
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.
wkernkamp is offline   Reply With Quote


doublets, matlab, panel, source

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
[Other] Adding solvers from DensityBasedTurbo to foam-extend 3.0 Seroga OpenFOAM Community Contributions 9 June 12, 2015 17:18
[swak4Foam] swak4foam building problem GGerber OpenFOAM Community Contributions 54 April 24, 2015 16:02
Problem compiling a custom Lagrangian library 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.