No description available.

<p>Boundary integral equation techniques are useful in thenumerical simulation of scattering problems for wave equations.Their advantage over methods based on partial di.erentialequations comes from the lack of phase errors in the wavepropagation and from the fact that only the boundary of thescattering object needs to be discretized. Boundary integraltechniques are often applied in frequency domain but recentlyseveral time domain integral equation methods are beingdeveloped.</p><p>We study time domain integral equation methods for thescalar wave equation with a Galerkin discretization of twodi.erent integral formulations for a Dirichlet scatterer. The.rst method uses the Kirchho. formula for the solution of thescalar wave equation. The method is prone to get unstable modesand the method is stabilized using an averaging .lter on thesolution. The second method uses the integral formulations forthe Helmholtz equation in frequency domain, and this method isstable. The Galerkin formulation for a Neumann scattererarising from Helmholtz equation is implemented, but isunstable.</p><p>In the discretizations, integrals are evaluated overtriangles, sectors, segments and circles. Integrals areevaluated analytically and in some cases numerically. Singularintegrands are made .nite, using the Du.y transform.</p><p>The Galerkin discretizations uses constant basis functionsin time and nodal linear elements in space. Numericalcomputations verify that the Dirichlet methods are stable, .rstorder accurate in time and second order accurate in space.Tests are performed with a point source illuminating a plateand a plane wave illuminating a sphere.</p><p>We investigate the On Surface Radiation Condition, which canbe used as a medium to high frequency approximation of theKirchho. formula, for both Dirichlet and Neumann scatterers.Numerical computations are done for a Dirichlet scatterer.</p>

No description available.

Boundary integral equation techniques are useful in thenumerical simulation of scattering problems for wave equations.Their advantage over methods based on partial di.erentialequations comes from the lack of phase errors in the wavepropagation and from the fact that only the boundary of thescattering object needs to be discretized. Boundary integraltechniques are often applied in frequency domain but recentlyseveral time domain integral equation methods are beingdeveloped. We study time domain integral equation methods for thescalar wave equation with a Galerkin discretization of twodi.erent integral formulations for a Dirichlet scatterer. The.rst method uses the Kirchho. formula for the solution of thescalar wave equation. The method is prone to get unstable modesand the method is stabilized using an averaging .lter on thesolution. The second method uses the integral formulations forthe Helmholtz equation in frequency domain, and this method isstable. The Galerkin formulation for a Neumann scattererarising from Helmholtz equation is implemented, but isunstable. In the discretizations, integrals are evaluated overtriangles, sectors, segments and circles. Integrals areevaluated analytically and in some cases numerically. Singularintegrands are made .nite, using the Du.y transform. The Galerkin discretizations uses constant basis functionsin time and nodal linear elements in space. Numericalcomputations verify that the Dirichlet methods are stable, .rstorder accurate in time and second order accurate in space.Tests are performed with a point source illuminating a plateand a plane wave illuminating a sphere. We investigate the On Surface Radiation Condition, which canbe used as a medium to high frequency approximation of theKirchho. formula, for both Dirichlet and Neumann scatterers.Numerical computations are done for a Dirichlet scatterer. / NR 20140805

Fluid phenomena dominated by viscous effects can, in many cases, be modeled by the Stokes equations. The boundary integral form of the Stokes equations reduces the number of degrees of freedom in a numerical discretization by reformulating the three-dimensional problem to two-dimensional integral equations to be discretized over the boundaries of the domain. Hence for the study of objects immersed in a fluid, such as drops or elastic/solid particles, integral equations are to be discretized over the surfaces of these objects only. As outer boundaries or confinements are added these must also be included in the formulation. An inherent difficulty in the numerical treatment of boundary integrals for Stokes flow is the integration of the singular fundamental solution of the Stokes equations, e.g. the so called Stokeslet. To alleviate this problem we developed a set of high-order quadrature rules for the numerical integration of the Stokeslet over a flat surface. Such a quadrature rule was first designed for singularities of the type <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?1/%7C%5Cmathbf%7Bx%7D%7C" />. To assess the convergence properties of this quadrature rule a theoretical analysis has been performed. The slightly more complicated singularity of the Stokeslet required certain modifications of the integration rule developed for <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?1/%7C%5Cmathbf%7Bx%7D%7C" />. An extension of this type of quadrature rule to a cylindrical surface is also developed. These quadrature rules are tested also on physical problems that have an analytic solution in the literature. Another difficulty associated with boundary integral problems is introduced by periodic boundary conditions. For a set of particles in a periodic domain periodicity is imposed by requiring that the motion of each particle has an added contribution from all periodic images of all particles all the way up to infinity. This leads to an infinite sum which is not absolutely convergent, and an additional physical constraint which removes the divergence needs to be imposed. The sum is decomposed into two fast converging sums, one that handles the short range interactions in real space and the other that sums up the long range interactions in Fourier space. Such decompositions are already available in the literature for kernels that are commonly used in boundary integral formulations. Here a decomposition in faster decaying sums than the ones present in the literature is derived for the periodic kernel of the stress tensor. However the computational complexity of the sums, regardless of the decomposition they stem from, is <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cmathcal%7BO%7D(N%5E%7B2%7D)" />. This complexity can be lowered using a fast summation method as we introduced here for simulating a sedimenting fiber suspension. The fast summation method was initially designed for point particles, which could be used for fibers discretized numerically almost without any changes. However, when two fibers are very close to each other, analytical integration is used to eliminate numerical inaccuracies due to the nearly singular behavior of the kernel and the real space part in the fast summation method was modified to allow for this analytical treatment. The method we have developed for sedimenting fiber suspensions allows for simulations in large periodic domains and we have performed a set of such simulations at a larger scale (larger domain/more fibers) than previously feasible. / <p>QC 20121122</p>

The importance of the Stokes system stems from the fact that the Stokes system is the stationary linearised form of the Navier Stokes system [Te01, Chapter1]. This linearisation is allowed when neglecting the inertial terms at a low Reinolds numbers Re << 1. The Stokes system essentially models the behaviour of a non - turbulent viscous fluid. The mixed interior boundary value problem related to the compressible Stokes system is reduced to two different BDIES which are equivalent to the original boundary value problem. These boundary-domain integral equation systems (BDIES) can be expressed in terms of surface and volume parametrix-based potential type operators whose properties are also analysed in appropriate Sobolev spaces. The invertibility and Fredholm properties related to the matrix operators that de ne the BDIES are also presented. Furthermore, we also consider the mixed compressible Stokes system with variable viscosity in unbounded domains. An analysis of the similarities and differences with regards to the bounded domain case is presented. Furthermore, we outline the mapping properties of the surface and volume parametrix-based potentials in weighted Sobolev spaces. Equivalence and invertibility results still hold under certain decay conditions on the variable coeffi cient The last part of the thesis refers to the mixed boundary value problem for the stationary heat transfer partial di erential equation with variable coe cient. This BVP is reduced to a system of direct segregated parametrix-based Boundary-Domain Integral Equations (BDIEs). We use a parametrix different from the one employed by Chkadua, Mikhailov and Natroshvili in the paper [CMN09]. Mapping properties of the potential type integral operators appearing in these equations are presented in appropriate Sobolev spaces. We prove the equivalence between the original BVP and the corresponding BDIE system. The invertibility and Fredholm properties of the boundary-domain integral operators are also analysed in both bounded and unbounded domains.

No description available.

The forced water entry and exit of two-dimensional bodies through a free surface is computed for various 2D bodies (symmetric wedges, asymmetric wedges, truncated wedges and boxes). These bodies enter or exit water with constant velocity or constant acceleration. The calculations are based on the fully non-linear timestepping complex-variable method of Vinje and Brevig. The model was formulated as an initial boundary-value problem with boundary conditions specified on the boundaries (dynamic and kinematic free-surface boundary conditions) and initial conditions at time zero (initial velocity and position of the body and free-surface particles). The formulated problem was solved by means of a boundary-element method using collocation points on the boundary of the domain and solutions at each time were calculated using time stepping (Runge-Kutta and Hamming predictor corrector) methods. Numerical results for the deformed free-surface profile, the speed of the point at the intersection of the body and free surface, the pressure along the wetted region of the bodies and force experienced by the bodies, are given for the entry and exit. To verify the results, various tests such as convergence checks, self-similarity for entry (gravity-free solutions) and Froude number effect for constant velocity entry and exit (half-wedge angles 5 up to 55 degrees) are investigated. The numerical results are compared with Mackie's analytical theory for water entry and exit with constant velocities, and the analytical added mass force computed for water entry and exit of symmetric wedges and boxes with constant acceleration and velocity using conformal mapping. Finally, numerical results showing the effect of finite depth are investigated for entry and exit.

This thesis is concerned with the study of transient scattering of acoustic waves by an obstacle in an infinite domain, where the scattered wave is represented in terms of time domain boundary layer potentials. The problem of finding the unknown solution of the scattering problem is thus reduced to the problem of finding the unknown density of the time domain boundary layer operators on the obstacle’s boundary, subject to the boundary data of the known incident wave. Using a Galerkin approach, the unknown density is replaced by a piecewise polynomial approximation, the coefficients of which can be found by solving a linear system. The entries of the system matrix of this linear system involve, for the case of a two dimensional scattering problem, integrals over four dimensional space-time manifolds. An accurate computation of these integrals is crucial for the stability of this method. Using piecewise polynomials of low order, the two temporal integrals can be evaluated analytically, leading to kernel functions for the spatial integrals with complicated domains of piecewise support. These spatial kernel functions are generalised into a class of admissible kernel functions. A quadrature scheme for the approximation of the two dimensional spatial integrals with admissible kernel functions is presented and proven to converge exponentially by using the theory of countably normed spaces. A priori error estimates for the Galerkin approximation scheme are recalled, enhanced and discussed. In particular, the scattered wave’s energy is studied as an alternative error measure. The numerical schemes are presented in such a way that allows the use of non-uniform meshes in space and time, in order to be used with adaptive methods that are based on a posteriori error indicators and which modify the computational domain according to the values of these error indicators. The theoretical analysis of these schemes demands the study of generalised mapping properties of time domain boundary layer potentials and integral operators, analogously to the well known results for elliptic problems. These mapping properties are shown for both two and three space dimensions. Using the generalised mapping properties, three types of a posteriori error estimators are adopted from the literature on elliptic problems and studied within the context of the two dimensional transient problem. Some comments on the three dimensional case are also given. Advantages and disadvantages of each of these a posteriori error estimates are discussed and compared to the a priori error estimates. The thesis concludes with the presentation of two adaptive schemes for the two dimensional scattering problem and some corresponding numerical experiments.

Une technique numérique basée sur le couplage d’une approximation par éléments finis et d’une méthode intégrale a été développée pour le calcul du signal induit dans les antennes I.R.M. Ce signal est issu du mouvement de précession libre de l’aimantation transversale du corps à explorer pendant la relaxation. Dans notre modélisation, l’aimantation transversale représente le champ magnétique source. Celui-ci induit dans l’antenne un courant d’une durée très brève (quelques millisecondes) ; il représente le signal contenant toutes les informations de l’échantillon. Notre modélisation des antennes d’I.R.M de type circuit imprimé a été validée par comparaison avec des mesures expérimentales ainsi qu’avec une méthode analytique. Nous l’avons développée en tenant compte de leurs géométries et de leurs caractéristiques électromagnétiques afin d’avoir un meilleur rapport Signal/Bruit. Nous avons pris en considération des principaux facteurs tels que la distance entre l’antenne et l’échantillon à explorer ainsi que les caractéristiques électromagnétiques de l’antenne. / A numerical technique, based on the combination of a finite element method and a boundary integral method, has been developed to compute the induced signal in MRI antennas. This signal rises from a free movement of precession of the transverse magnetization of the sample to explore. In our modeling, the transverse magnetization represents the magnetic source field. Its flux embraces the antenna to give rise to a sinusoidal current which is very quickly attenuated in time (a few ms); it represents the signal containing all the information of the sample. We here want to find the geometrical and electromagnetic characteristics of the antennas which permit to have a signal to noise ratio as great as possible. In our computation, we have taken into account leading factors such as the distance between the probe and the organ to be explored and also the geometrical and electromagnetic characteristics of the probe. Our modeling of printed circuits MRI antenna has been validated by comparing with experimental measurements and also with an anlytical method. We have developped it by taking into account their geometries and their electromagnetical characteristics in order to have a better signal/noise ratio. We have considered principal factors such as the distance between the antenna and the organ to explore and also the electromagnetic characteristics of the antenna.