Home     News     Team     Areas     Services     Research     Links  

Mathematical Engineering Associate Consultants S.A., INGMAT S.A., is a technological management company.

Our purpose is to generate and integrate scientific-technological developments to produce applied science with high added value for national and international industry.



 

José Miguel de la Barra 412 piso 4
Santiago, Chile
56 2 2474905
contactos@ingmat.com

RESEARCH

Explicit Representation for the Infinite-Depth Two-Dimensional Free-Surface Green's Function in Linear Water-Wave Theory

In this paper we derive an explicit representation for the two-dimensional free-surface Green's function in water of infinite depth, based on a finite combination of complex-valued exponential integrals and elementary functions. This representation can easily and accurately be evaluated in a numerical manner, and its main advantage over other representations lies in its simplicity and in the fact that it can be extended towards the complementary half-plane in a straightforward manner. It seems that this extension has not been studied rigorously until now, and it is required when boundary integral equations are extended in the same way. For the computation of the Green's function, the limiting absorption principle and a partial Fourier transform along the free surface are used. Some of its properties are also discussed, and an expression for its far field is developed, which allows us to state appropriately the involved radiation condition. This Green's function is then used to solve the two-dimensional infinite-depth water-wave problem by developing a corresponding boundary integral equation, whose solution is determined by means of the boundary element method. To validate the computations, a benchmark problem based on a half-circle is presented and solved numerically.


FULL TEXT

Theoretical aspects and numerical computation of the time-harmonic Green's function for an isotropic elastic half-plane with an impedance boundary condition

This work presents an effective and accurate method for determining, from a theoretical and computational point of view, the time-harmonic Green's function of an isotropic elastic half-plane where an impedance boundary condition is considered. This method, based on the previous work done by Durán et al. (cf. [Numer. Math.107 (2007) 295–314; IMA J. Appl. Math. 71 (2006) 853–876]) for the Helmholtz equation in a half-plane, combines appropriately analytical and numerical techniques, which has an important advantage because the obtention of explicit expressions for the surface waves. We show, in addition to the usual Rayleigh wave, another surface wave appearing in some special cases. Numerical results are given to illustrate that. This is an extended and detailed version of the previous article by Durán et al. [C. R. Acad. Sci. Paris, Ser. IIB 334 (2006) 725–731].


FULL TEXT

On the Green’s function for the Helmholtz operator in an impedance circular cylindrical waveguide

This paper addresses the problem of finding a series representation for the Green’s function of the Helmholtz operator in an infinite circular cylindrical waveguide with impedance boundary condition. Resorting to the Fourier transform, complex analysis techniques and the limiting absorption principle (when the undamped case is analyzed), a detailed deduction of the Green’s function is performed, generalizing the results available in the literature for the case of a complex impedance parameter. Procedures to obtain numerical values of the Green’s function are also developed in this article.


FULL TEXT

Radiation condition and uniqueness for the outgoing elastic wave in a half-plane with free boundary

In this Note we deduce an explicit Sommerfeld-type radiation condition which is convenient to prove the uniqueness for the time-harmonic outgoing wave problem in an isotropic elastic half-plane with free boundary condition. The expression is obtained from a rigorous asymptotic analysis of the associated Green's function. The main difficulty is that the free boundary condition allows the propagation of a Rayleigh wave which cannot be neglected in the far field expansion. We also give the existence result for this problem. To cite this article: M. Durán et al., C. R. Acad. Sci. Paris, Ser. I 347 (2009).


FULL TEXT

An Efficient Galerkin BEM to Compute High Acoustic Eigenfrequencies

An efficient numerical method, using integral equations, is developed to calculate precisely the acoustic eigenfrequencies and their associated eigenvectors, located in a given high frequency interval. It is currently known that the real symmetric matrices are well adapted to numerical treatment. However, we show that this is not the case when using integral representations to determine with high accuracy the spectrum of elliptic, and other related operators. Functions are evaluated only in the boundary of the domain, so very fine discretizations may be chosen to obtain high eigenfrequencies. We discuss the stability and convergence of the proposed method. Finally we show some examples.


FULL TEXT

The Helmholtz Equation in a Locally Perturbed Half-Space with Non-Absorbing Boundary

We obtain uniqueness and existence results of an outgoing solution for the Helmholtz equation in a half-space, or in a compact local perturbation of it, with an impedance boundary condition. It is worth noting that these kinds of domains have unbounded boundaries which lead to a non-classical exterior problem. The established radiation condition is somewhat different from the usual Sommerfeld’s one, due to the appearance of surface waves (in the case of a non-absorbing boundary). A half-space Green’s function framework is used to carry out our computations. This is an extended and detailed version of the previous article “The Helmholtz equation with impedance in a half-space,” Duran et al. (CR Acad Sci Paris Ser I 341:561–566, 2005).


FULL TEXT

Computing numerically the Green’s function of the half-plane Helmholtz operator with impedance boundary conditions

In this article we compute numerically the Green’s function of the half-plane Helmholtz operator with impedance boundary conditions. A compactly perturbed half-plane Helmholtz problem is used to motivate this calculation, by treating it through integral equation techniques. These require the knowledge of the calculated Green’s function, and lead to a boundary element discretization. The Green’s function is computed using the inverse Fourier operator of its spectral transform, applying an inverse FFT for the regular part, and removing the singularities analytically. Finally, some numerical results for the Green’s function and for a benchmark resonance problem are shown.


FULL TEXT

Computing Green's function of elasticity in a half-plane with impedance boundary condition

This Note presents an effective and accurate method for numerical calculation of the Green's function G associated with the time harmonic elasticity system in a half-plane, where an impedance boundary condition is considered. The need to compute this function arises when studying wave propagation in underground mining and seismological engineering. To theoretically obtain this Green's function, we have drawn our inspiration from the paper by Durán et al. (2005), where the Green's function for the Helmholtz equation has been computed. The method consists in applying a partial Fourier transform, which allows an explicit calculation of the so-called spectral Green's function. In order to compute its inverse Fourier transform, we separate  as a sum of two terms. The first is associated with the whole plane, whereas the second takes into account the half-plane and the boundary conditions. The first term corresponds to the Green's function of the well known time-harmonic elasticity system in (cf. J. Dompierre, Thesis). The second term is separated as a sum of three terms, where two of them contain singularities in the spectral variable (pseudo-poles and poles) and the other is regular and decreasing at infinity. The inverse Fourier transform of the singular terms are analytically computed, whereas the regular one is numerically obtained via an FFT algorithm. We present a numerical result. Moreover, we show that, under some conditions, a fourth additional slowness appears and which could produce a new surface wave. To cite this article: M. Durán et al., C. R. Mecanique 334 (2006).


FULL TEXT

The Helmholtz equation in a locally perturbed half-plane with passive boundary

In this article, we study the existence and uniqueness of outgoing solutions for the Helmholtz equation in locally perturbed half-planes with passive boundary. We establish an explicit outgoing radiation condition which is somewhat different from the usual Sommerfeld's one due to the appearance of surface waves. We work with the help of Fourier analysis and a half-plane Green's function framework. This is an extended and detailed version of the previous article Durán et al. (2005, The Helmholtz equation with impedance in a half-plane. C. R. Acad. Sci. Paris, Ser. I, 340, 483–488).


FULL TEXT

Numerical study of the spectral 3-D Green's function singularities for piezoelectric SAW components

This article deals with the numerical study of the singularities appearing in the spectral 3-D Green's function associated with the piezoelectric surface acoustic wave components (so-called SAW components). These electrical units are currently used today in several devices produced by the telecommunications industry (radio, TV, radar, arid digital telecommunication systems). The need to improve their performance has motivated engineers arid researchers to use mathematical modeling intensively, in particular the integral equations technique here used, which requires the computing of the associated Green's function and the study of its properties.


FULL TEXT

The Helmholtz equation with impedance in a half-space
In this Note we obtain existence and uniqueness results for the Helmholtz equation in the half-space View the MathML source with an impedance or Robin boundary condition. Basically, we follow the procedure we have already used in the bi-dimensional case (the half-plane). Thus, we compute the associated Green's function with the help of a double Fourier transform and we analyze its far field in order to obtain radiation conditions that allow us to prove the uniqueness of an outgoing solution. Again, these radiation conditions are somewhat unusual due to the appearance of a surface wave guided by the boundary. An integral representation of the solution is presented by means of the Green's function and the boundary data.
FULL TEXT

A numerical study of a spectral problem in solid-fluid type structures
This article presents a numerical study of a spectral problem that models the vibrations of a solid–fluid structure. It is a quadratic eigenvalue problem involving incompressible Stokes equations. In its numerical approximation we use Lagrange finite elements. To approximate the velocity, degree 2 polynomials on triangles are used, and for the pressure, degree 1 polynomials. The numerical results obtained confirm the theory, as they show in particular that the known theoretical bound for the maximum number of nonreal eigenvalues admitted by such a system is optimal. The results also take account of the dependence of vibration frequencies with respect to determined physical parameters, which have a bearing on the model.
FULL TEXT

Numerical stability in the calculation of eigenfrequencies using integral equations
We comment on a phenomenon of instability that appears while computing eigenfrequencies using the integral equation framework. More precisely, it is currently known that the real symmetric matrices are well, and sometimes the best, adapted to numerical treatment. However, we show that this is not the case, if we wish to determine with high accuracy the spectrum of elliptic, and other related operators, using integral representations.
FULL TEXT

Rate of convergence estimates for the spectral approximation of a generalized eigenvalue problem
The aim of this work is to derive rate of convergence estimates for the spectral approximation of a mathematical model which describes the vibrations of a solid-fluid type structure. First, we summarize the main theoretical results and the discretization of this variational eigenvalue problem. Then, we state some well known abstract theorems on spectral approximation and apply them to our specific problem, which allow us to obtain the desired spectral convergence. By using classical regularity results, we are able to establish estimates for the rate of convergence of the approximated eigenvalues and for the gap between generalized eigenspaces.
FULL TEXT

A quadratic eigenvalue problem involving Stokes equations
An existence theorem for the eigenvalues of a spectral problem is studied in this paper. The physical situation behind this mathematical problem is the determination of the eigenfrequencies and eigenmotions of a fluid-solid structure. The liquid part in this structure is represented by a viscous incompressible fluid, while the solid part is a set of parallel rigid tubes. The spectral problem governing this system is a quadratic eigenvalue problem which involves Stokes equations with a non-local boundary condition. The strategy for tackling the question of existence of eigenvalues consists of proving that the original problem is equivalent to that of determining the characteristic values of a linear (non-selfadjoint) compact operator. Sharp estimates for the eigenvalues give precise information about the region of ω where the eigenvalues are located. In particular, we prove that this problem admits a countable set of eigenvalues in which only a finite number of them have a non-zero imaginary part.
FULL TEXT