Constant potential method (2D) - OCW UPM恒电位法（2D）- OCW芬欧汇川集团.pptVIP

UPM IDR Dirichlet formulation with constant potential Ana Laverón Simavilla Ma Victoria Lapuerta González Dirichlet Formulation The integral equation must be solved making the point P tend toward the body surface Doublet potential: Fd Source distribution Doublet distribution Doublet potential: Fd Dirichlet Formulation Making the inner potential zero: Doublet distribution =G Constant potential method (2D) The profile’s contour is modeled with N panels, that will be straight segments. Panels are determined by N+1 nodes Collocation points are the midpoints of each panel: Node of trailing edge node panel j node j node j+1 Path followed around the profile Constant potential method (2D) The following conditions are imposed: Constant potential on each panel The previous equation is rewritten for the collocation points k G=FN-F1 Constant potential method (2D) To solve the integrals on each panel it’s suitable to choose axes bound to the panel (u,v) v u C.P. k j j+1 panel j k j j+1 Constant potential method (2D) For the panel, k=j For the layer surface, a semi-infinite panel is used Finally: Constant potential method (2D) Matrix form: with and Constant potential method (2D) Flow around a Karman-Trefftz profile a x y R x0 y0 q0 a a ka We will compare the numerical and the analytic results for a Kármán-Trefftz profile. For n panels calculate the pressure coefficient on the profile’s lower and upper, . For this, calculate the speed on the nodes using the expression: where di is the distance between the collocation points i and i+1. Calculate the global circulation around the airfoil: Compare the results with those obtained exactly using the Kármán-Trefftz transformation. Required calculations Comments for the resolution Function that gives the nodes’ coordinates: function[