Numerical method for hypersingular integrals of highly oscillatory functions on the positive semiaxis

This paper deals with a quadrature rule for the numerical evaluation of hypersingular integrals of highly oscillatory functions on the positive semiaxis. The rule is of product type and consists in approximating the density function f by a truncated interpolation process based on the zeros of generalized Laguerre polynomials and an additional point. We prove the stability and the convergence of the rule, giving error estimates for functions belonging to weighted Sobolev spaces equipped with uniform norm. We also show how the proposed rule can be used for the numerical solution of hypersingular integral equations. Numerical tests which confirm the theoretical estimates and comparisons with other existing quadrature rules are presented

Approximation of Hilbert and Hadamard transforms on (0, +∞)

The authors propose a numerical method for computing Hilbert and Hadamard transforms on (0, +∞) by a simultaneous approximation process involving a suitable Lagrange polynomial of degree s and “truncated” Gaussian rule of order m, with s&lt

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

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.Comment: 34 page

Modelling Fluid Structure Interaction problems using Boundary Element Method

This dissertation investigates the application of Boundary Element Methods (BEM) to Fluid Structure Interaction (FSI) problems under three main different perspectives. This work is divided in three main parts: i) the derivation of BEM for the Laplace equation and its application to analyze ship-wave interaction problems, ii) the imple- mentation of efficient and parallel BEM solvers addressing the newest challenges of High Performance Computing, iii) the developing of a BEM for the Stokes system and its application to study micro-swimmers.First we develop a BEM for the Laplace equation and we apply it to predict ship-wave interactions making use of an innovative coupling with Finite Element Method stabilization techniques. As well known, the wave pattern around a body depends on the Froude number associated to the flow. Thus, we throughly investigate the robustness and accuracy of the developed methodology assessing the solution dependence on such parameter. To improve the performance and tackle problems with higher number of unknowns, the BEM developed for the Laplace equation is parallelized using OpenSOURCE tech- nique in a hybrid distributed-shared memory environment. We perform several tests to demonstrate both the accuracy and the performance of the parallel BEM developed. In addition, we explore two different possibilities to reduce the overall computational cost from O(N2) to O(N). Firstly we couple the library with a Fast Multiple Method that allows us to reach for higher order of complexity and efficiency. Then we perform a preliminary study on the implementation of a parallel Non Uniform Fast Fourier Transform to be coupled with the newly developed algorithm Sparse Cardinal Sine De- composition (SCSD).Finally we consider the application of the BEM framework to a different kind of FSI problem represented by the Stokes flow of a liquid medium surrounding swimming micro-organisms. We maintain the parallel structure derived for the Laplace equation even in the Stokes setting. Our implementation is able to simulate both prokaryotic and eukaryotic organisms, matching literature and experimental benchmarks. We finally present a deep analysis of the importance of hydrodynamic interactions between the different parts of micro-swimmers in the prevision of optimal swimming conditions, focusing our attention on the study of flagellated \u201crobotic\u201d composite swimmers

Numerical computation of hypersingular integrals on the real semiaxis

In this paper we propose some different strategies to approximate hypersingular integrals. Hadamard Finite Part integrals (shortly FP integrals), regarded as p th derivative of Cauchy principal value integrals, are of interest in the solution of hypersingular BIE, which model many different kind of Physical and Engineering problems (see [1] and the references therein, [2], [3, 4]). The procedure we employ here is based on a simple tool like the “truncated”Gaussian rule (see [5]), conveniently modified to remove numerical cancellation. We will consider density functions having different decays at infinity. The method is shown to be numerically stable and convergent and some error estimates in suitable Zygmund-type spaces are proved. Finally, some numerical tests which confirm the efficiency of the proposed procedures are presented