Quadrature methods for integro-differential equations of Prandtl's type in weighted spaces of continuous functions

Maria Carmela De Bonis,Donatella Occorsio

DOI: https://doi.org/10.48550/arXiv.2008.00294

2020-08-02

Abstract:The paper deals with the approximate solution of integro-differential equations of Prandtl's type. Quadrature methods involving ``optimal'' Lagrange interpolation processes are proposed and conditions under which they are stable and convergent in suitable weighted spaces of continuous functions are proved. The efficiency of the method has been tested by some numerical experiments, some of them including comparisons with other numerical procedures. In particular, as an application, we have implemented the method for solving Prandtl's equation governing the circulation air flow along the contour of a plane wing profile, in the case of elliptic or rectangular wing-shape.

Numerical Analysis

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is the numerical solution of Prandtl - type integro - differential equations (IDE). Specifically, the authors are concerned with how to propose stable and convergent quadrature methods in weighted continuous function spaces to approximately solve this type of equations.

### 1. Research Background
Prandtl - type integro - differential equations are widely used in many problems in physics and engineering, especially in aerodynamics to describe the circulation distribution around the wing profile. The general form of this type of equations is:

\[
\sigma(y)\zeta(y)+a\zeta'(y)+\frac{b}{\pi}\int_{ - 1}^{1}\frac{\zeta'(x)}{x - y}dx+\frac{1}{\pi}\int_{ - 1}^{1}\bar{k}(x,y)\zeta(x)dx = g(y),\quad y\in(-1,1)
\]

where $\sigma(y)$, $\bar{k}(x,y)$ and $g(y)$ are given functions, $a$ and $b$ are constants satisfying $a^{2}+b^{2}=1$, $\zeta(y)$ is the unknown solution, and satisfies the boundary conditions $\zeta(-1)=\zeta(1) = 0$.

### 2. Solutions
To solve the above equations, the authors proposed quadrature methods based on the "optimal" Lagrange interpolation process. These methods are stable and convergent in the appropriate weighted continuous function spaces. The specific steps are as follows:

1. **Equation Rewriting**: Use the zero - boundary conditions to rewrite the original equation as:
\[
(M_{\sigma\rho}+DA_{\rho}+K + H)f(y)=g(y)
\]
where,
\[
(M_{\sigma\rho}f)(y)=(\sigma\rho f)(y),\quad(Dq)(y)=\frac{d}{dy}q(y),
\]
\[
(A_{\rho}f)(y)=a(f\rho)(y)+\frac{b}{\pi}\int_{ - 1}^{1}\frac{(f\rho)(x)}{x - y}dx,
\]
\[
(Kf)(y)=\frac{1}{\pi}\int_{ - 1}^{1}k(x,y)(f\rho)(x)dx,\quad(Hf)(y)=\frac{1}{\pi}\int_{ - 1}^{1}h(x,y)(f\rho)(x)dx.
\]

2. **Numerical Methods**: Two quadrature methods were introduced, which use the optimal Lagrange interpolation process, and the conditions to ensure their stability and convergence were determined. Error estimates and the condition number of the final linear system were also studied.

3. **Numerical Experiments**: The effectiveness of the proposed methods was verified through some numerical experiments and compared with other numerical methods.

### 3. Application Examples
As an application example, the authors implemented this method to solve the Prandtl equations describing the shapes of elliptical wings and rectangular wings in order to calculate the circulation distribution along the wing profile.

In conclusion, this paper aims to develop an effective and stable numerical method for solving Prandtl - type integro - differential equations, and verifies the effectiveness of this method through theoretical analysis and numerical experiments.