Controllability method for acoustic scattering with spectral elements

We formulate the Helmholtz equation as an exact controllability problem for the time-dependent wave equation. The problem is then discretized in time domain with central finite difference scheme and in space domain with spectral elements. This approach leads to high accuracy in spatial discretization. Moreover, the spectral element method results in diagonal mass matrices, which makes the time integration of the wave equation highly efficient. After discretization, the exact controllability problem is reformulated as a least-squares problem, which is solved by the conjugate gradient method. We illustrate the method with some numerical experiments, which demonstrate the significant improvements in efficiency due to the higher order spectral elements. For a given accuracy, the controllability technique with spectral element method requires fewer computational operations than with conventional finite element method. In addition, by using higher order polynomial basis the influence of the pollution effect is reduced.     

Preconditioning techniques for the iterative solution of scattering problems

We consider a time-harmonic electromagnetic scattering problem for an inhomogeneous medium. Some symmetry hypotheses on the refractive index of the medium and on the electromagnetic fields allow to reduce this problem to a two-dimensional scattering problem. This boundary value problem is defined on an unbounded domain, so its numerical solution cannot be obtained by a straightforward application of usual methods, such as for example finite difference methods, and finite element methods. A possible way to overcome this difficulty is given by an equivalent integral formulation of this problem, where the scattered field can be computed from the solution of a Fredholm integral equation of second kind. The numerical approximation of this problem usually produces large dense linear systems. We consider usual iterative methods for the solution of such linear systems, and we study some preconditioning techniques to improve the efficiency of these methods. We show some numerical results obtained with two well known Krylov subspace methods, i.e., Bi-CGSTAB and GMRES.     

Domain decomposition iterative procedures for solving scalar waves in the frequency domain

The propagation of dispersive waves can be modeled relevantly in the frequency domain. A wave problem in the frequency domain is difficult to solve numerically. In addition to having a complex–valued solution, the problem is neither Hermitian symmetric nor coercive in a wide range of applications in Geophysics or Quantum–Mechanics. In this paper, we consider a parallel domain decomposition iterative procedure for solving the problem by finite differences or conforming finite element methods. The analysis includes the decomposition of the domain into either the individual elements or larger subdomains ( of finite elements). To accelerate the speed of convergence, we introduce relaxation parameters on the subdomain interfaces and an artificial damping iteration. The convergence rate of the resulting algorithm turns out to be independent on the mesh size and the wave number. Numerical results carried out on an nCUBE2 parallel computer are presented to show the effectiveness of the method. Received October 30, 1995 / Revised version received January 10, 1997     

A finite element method for approximating the time-harmonic Maxwell equations

Summary In this paper we study the use of Nédélec's curl conforming finite elements to approximate the time-harmonic Maxwell equations on a bounded domain. The analysis is complicated by the fact that the bilinear form is not coercive, and the principle part has a very large null-space. This difficulty is circumvented by using a discrete Helmholtz decomposition of the error vector. Numerical results are presented that compare two different linear elements.Research supported in part by grants from AFOSR and NSF     

Locally linearized fractional step methods for nonlinear parabolic problems

This work deals with the efficient numerical solution of a class of nonlinear time-dependent reaction-diffusion equations. Via the method of lines approach, we first perform the spatial discretization of the original problem by applying a mimetic finite difference scheme. The system of ordinary differential equations arising from that process is then integrated in time with a linearly implicit fractional step method. For that purpose, we locally decompose the discrete nonlinear diffusion operator using suitable Taylor expansions and a domain decomposition splitting technique. The totally discrete scheme considers implicit time integrations for the linear terms while explicitly handling the nonlinear ones. As a result, the original problem is reduced to the solution of several linear systems per time step which can be trivially decomposed into a set of uncoupled parallelizable linear subsystems. The convergence of the proposed methods is illustrated by numerical experiments.     

A GLOBAL-LOCAL FINITE ELEMENT METHOD IN SPACE-TIME FOR A HYPERBOLIC PROBLEM

ABSTRACT

A finite element method for solving the wave equation with couples boundary conditions is presented. In this approach finite elements are applied globally with respect to space and simultaneously but locally with respect to time. This gives rise to a single-step method in time. The method is a practical and economic one and the numerical results obtained compare favourably with the available analytic solution.     

Modified block preconditioners for the discretized time-harmonic Maxwell equations in mixed form

In this paper, building on the previous work by Greif and Schötzau [Preconditioners for the discretized time-harmonic Maxwell equations in mixed form, Numer. Linear Algebra Appl. 14 (2007) 281–297] and Benzi and Olshanskii [An augmented lagrangian-based approach to the Oseen problem, SIAM J. Sci. Comput. 28 (2006) 2095–2113], we present the improved preconditioning techniques for the iterative solution of the saddle point linear systems, which arise from the finite element discretization of the mixed formulation of the time-harmonic Maxwell equations. The modified block diagonal and triangular preconditioners considered are based on augmentation with using the symmetric nonsingular weighted matrix. We discuss the spectral properties of the preconditioned matrix in detail and generalize the results of the above-mentioned paper by Greif and Schötzau. Numerical experiments are given to demonstrate the efficiency of the presented preconditioners.     

Finite element analysis of varitional crimes for a quasilinear elliptic problem in 3D

Summary. We examine a finite element approximation of a quasilinear boundary value elliptic problem in a three-dimensional bounded convex domain with a smooth boundary. The domain is approximated by a polyhedron and a numerical integration is taken into account. We apply linear tetrahedral finite elements and prove the convergence of approximate solutions on polyhedral domains in the -norm to the true solution without any additional regularity assumptions. Received May 23, 1997 / Published online December 6, 1999     

A collocation method for high-frequency scattering by convex polygons

We consider the problem of scattering of a time-harmonic acoustic incident plane wave by a sound soft convex polygon. For standard boundary or finite element methods, with a piecewise polynomial approximation space, the computational cost required to achieve a prescribed level of accuracy grows linearly with respect to the frequency of the incident wave. Recently Chandler–Wilde and Langdon proposed a novel Galerkin boundary element method for this problem for which, by incorporating the products of plane wave basis functions with piecewise polynomials supported on a graded mesh into the approximation space, they were able to demonstrate that the number of degrees of freedom required to achieve a prescribed level of accuracy grows only logarithmically with respect to the frequency. Here we propose a related collocation method, using the same approximation space, for which we demonstrate via numerical experiments a convergence rate identical to that achieved with the Galerkin scheme, but with a substantially reduced computational cost.     

Fast parallel solvers for symmetric boundary element domain decomposition equations

Summary. The boundary element method (BEM) is of advantage in many applications including far-field computations in magnetostatics and solid mechanics as well as accurate computations of singularities. Since the numerical approximation is essentially reduced to the boundary of the domain under consideration, the mesh generation and handling is simpler than, for example, in a finite element discretization of the domain. In this paper, we discuss fast solution techniques for the linear systems of equations obtained by the BEM (BE-equations) utilizing the non-overlapping domain decomposition (DD). We study parallel algorithms for solving large scale Galerkin BE–equations approximating linear potential problems in plane, bounded domains with piecewise homogeneous material properties. We give an elementary spectral equivalence analysis of the BEM Schur complement that provides the tool for constructing and analysing appropriate preconditioners. Finally, we present numerical results obtained on a massively parallel machine using up to 128 processors, and we sketch further applications to elasticity problems and to the coupling of the finite element method (FEM) with the boundary element method. As shown theoretically and confirmed by the numerical experiments, the methods are of algebraic complexity and of high parallel efficiency, where denotes the usual discretization parameter. Received August 28, 1996 / Revised version received March 10, 1997     

Solution of Maxwell equation in axisymmetric geometry by Fourier series decompostion and by use of H (rot) conforming finite element

Summary. This study deals with the mathematical and numerical solution of time-harmonic Maxwell equation in axisymmetric geometry. Using Fourier decomposition, we define weighted Sobolev spaces of solution and we prove expected regularity results. A practical contribution of this paper is the construction of a class of finite element conforming with the H (rot) space equipped with the weighted measure rdrdz. It appears as an extension of the well-known cartesian mixed finite element of Raviart-Thomas-Nédélec [11]–[15]. These elements are built from classical lagrangian and mixed finite element, therefore no special approximations functions are needed. Finally, following works of Mercier and Raugel [10], we perform an interpolation error estimate for the simplest proposed element. Received March 15, 1996 / Revised version received November 30, 1998 / Published online December 6, 1999     

Schwarz methods for inequalities with contraction operators

We prove the convergence of some multiplicative and additive Schwarz methods for inequalities which contain contraction operators. The problem is stated in a reflexive Banach space and it generalizes the well-known fixed-point problem in the Hilbert spaces. Error estimation theorems are given for three multiplicative algorithms and two additive algorithms. We show that these algorithms are in fact Schwarz methods if the subspaces are associated with a decomposition of the domain. Also, for the one- and two-level methods in the finite element spaces, we write the convergence rates as functions of the overlapping and mesh parameters. They are similar with the convergence rates of these methods for linear problems. Besides the direct use of the five algorithms for the inequalities with contraction operators, we can use the above results to obtain the convergence rate of the Schwarz method for other types of inequalities or nonlinear equations. In this way, we prove the convergence and estimate the error of the one- and two-level Schwarz methods for some inequalities in Hilbert spaces which are not of the variational type, and also, for the Navier–Stokes problem. Finally, we give conditions of existence and uniqueness of the solution for all problems we consider. We point out that these conditions and the convergence conditions of the proposed algorithms are of the same type.     

A fictitious domain approach to the numerical solution of PDEs in stochastic domains

We present an efficient method for the numerical realization of elliptic PDEs in domains depending on random variables. Domains are bounded, and have finite fluctuations. The key feature is the combination of a fictitious domain approach and a polynomial chaos expansion. The PDE is solved in a larger, fixed domain (the fictitious domain), with the original boundary condition enforced via a Lagrange multiplier acting on a random manifold inside the new domain. A (generalized) Wiener expansion is invoked to convert such a stochastic problem into a deterministic one, depending on an extra set of real variables (the stochastic variables). Discretization is accomplished by standard mixed finite elements in the physical variables and a Galerkin projection method with numerical integration (which coincides with a collocation scheme) in the stochastic variables. A stability and convergence analysis of the method, as well as numerical results, are provided. The convergence is “spectral” in the polynomial chaos order, in any subdomain which does not contain the random boundaries.     

Two-scale composite finite element method for Dirichlet problems on complicated domains

In this paper, we define a new class of finite elements for the discretization of problems with Dirichlet boundary conditions. In contrast to standard finite elements, the minimal dimension of the approximation space is independent of the domain geometry and this is especially advantageous for problems on domains with complicated micro-structures. For the proposed finite element method we prove the optimal-order approximation (up to logarithmic terms) and convergence estimates valid also in the cases when the exact solution has a reduced regularity due to re-entering corners of the domain boundary. Numerical experiments confirm the theoretical results and show the potential of our proposed method.     

Estimates from above and from below for the homogenized coefficients

Summary The computation of the homogenized coefficients involves the solution of a minimum problem for a quadratic functional. The use of conforming finite elements provides estimates from above. In this paper we prove that, under realistic assumptions, the use of linear non-conforming finite elements provides estimates from below. Hence, a-posteriori error estimates for the homogenized coefficients can be obtained.     

Numerical simulation of tridimensional electromagnetic shaping of liquid metals

Summary We describe a numerical method to compute free surfaces in electromagnetic shaping and levitation of liquid metals. We use an energetic variational formulation and optimization techniques to compute, a critical point. The surfaces are represented by piecewise linear finite elements. Each step of the algorithm requires solving an elliptic boundary value problem in the exterior of the intermediate surfaces. This is done by using an integral representation on these surfaces.     

A preconditioned domain decomposition algorithm for the solution of the elliptic Neumann problem

A new preconditioned conjugate gradient (PCG)-based domain decomposition method is given for the solution of linear equations arising in the finite element method applied to the elliptic Neumann problem. The novelty of the proposed method is in the recommended preconditioner which is constructed by using cyclic matrix. The resulting preconditioned algorithms are well suited to parallel computation.     

An elliptic singularly perturbed problem with two parameters II: Robust finite element solution

In this paper we consider a singularly perturbed elliptic model problem with two small parameters posed on the unit square. The problem is solved numerically by the finite element method using piecewise linear or bilinear elements on a layer-adapted Shishkin mesh. We prove that method with bilinear elements is uniformly convergent in an energy norm. Numerical results confirm our theoretical analysis.     

An ultra-weak method for acoustic fluid–solid interaction

We introduce the ultra-weak variational formulation (UWVF) for fluid–solid vibration problems. In particular, we consider the scattering of time-harmonic acoustic pressure waves from solid, elastic objects. The problem is modeled using a coupled system of the Helmholtz and Navier equations. The transmission conditions on the fluid–solid interface are represented in an impedance-type form after which we can employ the well known ultra-weak formulations for the Helmholtz and Navier equations. The UWVF approximation for both equations is computed using a superposition of propagating plane waves. A condition number based criterion is used to define the plane wave basis dimension for each element. As a model problem we investigate the scattering of sound from an infinite elastic cylinder immersed in a fluid. A comparison of the UWVF approximation with the analytical solution shows that the method provides a means for solving wave problems on relatively coarse meshes. However, particular care is needed when the method is used for problems at frequencies near the resonance frequencies of the fluid–solid system.     

A cascadic multigrid algorithm for semilinear elliptic problems

Summary. We propose a cascadic multigrid algorithm for a semilinear elliptic problem. The nonlinear equations arising from linear finite element discretizations are solved by Newton's method. Given an approximate solution on the coarsest grid on each finer grid we perform exactly one Newton step taking the approximate solution from the previous grid as initial guess. The Newton systems are solved iteratively by an appropriate smoothing method. We prove that the algorithm yields an approximate solution within the discretization error on the finest grid provided that the start approximation is sufficiently accurate and that the initial grid size is sufficiently small. Moreover, we show that the method has multigrid complexity. Received February 12, 1998 / Revised version received July 22, 1999 / Published online June 8, 2000