Scattering Calculations with Wavelets

We show that the use of wavelet bases for solving the momentum-space scattering integral equation leads to sparse matrices which can simplify the solution. Wavelet bases are applied to calculate the K-matrix for nucleon-nucleon scattering with the s-wave Malfliet-Tjon V potential. We introduce a new method, which uses special properties of the wavelets, for evaluating the singular part of the integral. Analysis of this test problem indicates that a significant reduction in computational size can be achieved for realistic few-body scattering problems.Comment: 26 pages, Latex, 6 eps figure

Wavelet Methods in the Relativistic Three-Body Problem

In this paper we discuss the use of wavelet bases to solve the relativistic three-body problem. Wavelet bases can be used to transform momentum-space scattering integral equations into an approximate system of linear equations with a sparse matrix. This has the potential to reduce the size of realistic three-body calculations with minimal loss of accuracy. The wavelet method leads to a clean, interaction independent treatment of the scattering singularities which does not require any subtractions.Comment: 14 pages, 3 figures, corrected referenc

Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D

We present an effective harmonic density interpolation method for the numerical evaluation of singular and nearly singular Laplace boundary integral operators and layer potentials in two and three spatial dimensions. The method relies on the use of Green's third identity and local Taylor-like interpolations of density functions in terms of harmonic polynomials. The proposed technique effectively regularizes the singularities present in boundary integral operators and layer potentials, and recasts the latter in terms of integrands that are bounded or even more regular, depending on the order of the density interpolation. The resulting boundary integrals can then be easily, accurately, and inexpensively evaluated by means of standard quadrature rules. A variety of numerical examples demonstrate the effectiveness of the technique when used in conjunction with the classical trapezoidal rule (to integrate over smooth curves) in two-dimensions, and with a Chebyshev-type quadrature rule (to integrate over surfaces given as unions of non-overlapping quadrilateral patches) in three-dimensions

Planewave density interpolation methods for 3D Helmholtz boundary integral equations

This paper introduces planewave density interpolation methods for the regularization of weakly singular, strongly singular, hypersingular and nearly singular integral kernels present in 3D Helmholtz surface layer potentials and associated integral operators. Relying on Green's third identity and pointwise interpolation of density functions in the form of planewaves, these methods allow layer potentials and integral operators to be expressed in terms of integrand functions that remain smooth (at least bounded) regardless the location of the target point relative to the surface sources. Common challenging integrals that arise in both Nystr\"om and boundary element discretization of boundary integral equation, can then be numerically evaluated by standard quadrature rules that are irrespective of the kernel singularity. Closed-form and purely numerical planewave density interpolation procedures are presented in this paper, which are used in conjunction with Chebyshev-based Nystr\"om and Galerkin boundary element methods. A variety of numerical examples---including problems of acoustic scattering involving multiple touching and even intersecting obstacles, demonstrate the capabilities of the proposed technique

Development of an integrated BEM approach for hot fluid structure interaction

The progress made toward the development of a boundary element formulation for the study of hot fluid-structure interaction in Earth-to-Orbit engine hot section components is reported. The convective viscous integral formulation was derived and implemented in the general purpose computer program GP-BEST. The new convective kernel functions, in turn, necessitated the development of refined integration techniques. As a result, however, since the physics of the problem is embedded in these kernels, boundary element solutions can now be obtained at very high Reynolds number. Flow around obstacles can be solved approximately with an efficient linearized boundary-only analysis or, more exactly, by including all of the nonlinearities present in the neighborhood of the obstacle. The other major accomplishment was the development of a comprehensive fluid-structure interaction capability within GP-BEST. This new facility is implemented in a completely general manner, so that quite arbitrary geometry, material properties and boundary conditions may be specified. Thus, a single analysis code (GP-BEST) can be used to run structures-only problems, fluids-only problems, or the combined fluid-structure problem. In all three cases, steady or transient conditions can be selected, with or without thermal effects. Nonlinear analyses can be solved via direct iteration or by employing a modified Newton-Raphson approach

Aerodynamic influence coefficient method using singularity splines

A numerical lifting surface formulation, including computed results for planar wing cases is presented. This formulation, referred to as the vortex spline scheme, combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of loading function methods. The formulation employes a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfied all the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise. Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral. The current formulation uses the elementary horseshoe vortex as the basic singularity and is therefore restricted to linearized potential flow. As part of the study, a non planar development was considered, but the numerical evaluation of the lifting surface concept was restricted to planar configurations. Also, a second order sideslip analysis based on an asymptotic expansion was investigated using the singularity spline formulation

Fredholm factorization of Wiener-Hopf scalar and matrix kernels

A general theory to factorize the Wiener-Hopf (W-H) kernel using Fredholm Integral Equations (FIE) of the second kind is presented. This technique, hereafter called Fredholm factorization, factorizes the W-H kernel using simple numerical quadrature. W-H kernels can be either of scalar form or of matrix form with arbitrary dimensions. The kernel spectrum can be continuous (with branch points), discrete (with poles), or mixed (with branch points and poles). In order to validate the proposed method, rational matrix kernels in particular are studied since they admit exact closed form factorization. In the appendix a new analytical method to factorize rational matrix kernels is also described. The Fredholm factorization is discussed in detail, supplying several numerical tests. Physical aspects are also illustrated in the framework of scattering problems: in particular, diffraction problems. Mathematical proofs are reported in the pape