The mortar spectral element method in domains of operators. Part I: The divergence operator and Darcy's equations

The mortar spectral element method is a domain decompositiontechnique that allows for discretizing second- or fourth-orderelliptic equations when set in standard Sobolev spaces. Theaim of this paper is to extend this method to some problemsformulated in spaces of square-integrable functions with square-integrabledivergence. A discretization of the equations due to Darcy whichmodel the flow in porous media is proposed. The numerical analysisof the discrete problem is performed and 2D numerical experimentsare presented; they turn out to be in good agreement with thetheoretical results. Résumé: La méthode d'élémentsspectraux avec joints est une technique de décompositionde domaine permettant de discrétiser des équationselliptiques d'ordre 2 ou 4 posés dans des espaces deSobolev usuels. Le but de cet article est d'étendre cetteméthode à certains problèmes variationnelsformulés dans des espaces de fonctions de carréintégrable à divergence de carré intégrable.On propose une discrétisation des équations deDarcy qui modélisent l'écoulement dans des milieuxporeux, on en effectue l'analyse numérique et on présentedes expériences numériques en dimension 2 quis'avèrent cohérentes avec les résultatsde l'analyse. Received on 6 April 2004. revised on 24 January 2005.     

Solution of the spectral problem for the curl and Stokes operators with periodic boundary conditions

In this paper, we establish relations between eigenvalues and eigenfunctions of the curl operator and Stokes operator (with periodic boundary conditions). These relations show that the curl operator is the square root of the Stokes operator with ν = 1. The multiplicity of the zero eigenvalue of the curl operator is infinite. The space L ₂(Q, 2π) is decomposed into a direct sum of eigenspaces of the operator curl. For any complex number λ, the equation rot u + λu = f and the Stokes equation −ν(Δv + λ ²v) + ∇p = f, div v = 0, are solved. Bibliography: 15 titles. Dedicated to the memory of Olga Aleksandrovna Ladyzhenskaya __________ Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 318, 2004, pp. 246–276.     

Mortar spectral element discretization of the stokes problem in axisymmetric domains

The Stokes problem in a tri‐dimensional axisymmetric domain results into a countable family of two‐dimensional problems when using the Fourier coefficients with respect to the angular variable. Relying on this dimension reduction, we propose and study a mortar spectral element discretization of the problem. Numerical experiments confirm the efficiency of this method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 44–73, 2014     

The double-layer potential operator over polyhedral domains II: Spline Galerkin methods

We examine the numerical approximation of the integral equation (λ ? K)u =f, where K is the double layer (harmonic) potential operator on a closed polyhedral surface in ?³ and λ, ∣λ∣≥1, is a complex constant. The solution is approximated by Galerkin's method, which is based on piecewise polynomials of arbitrary degree on graded triangulations. By utilizing spline spaces which are modified in that the trial functions vanish on some of the triangles closest to the vertices and edges, we investigate the stability of this method in L². Furthermore, the use of suitably graded meshes leads to the same quasioptimal error estimates as in the case of a smooth surface.     

The Fourier Singular Complement Method for the Poisson problem. Part II: axisymmetric domains

This paper is the second part of a threefold article, aimed at solving numerically the Poisson problem in three-dimensional prismatic or axisymmetric domains. In the first part of this series, the Fourier Singular Complement Method was introduced and analysed, in prismatic domains. In this second part, the FSCM is studied in axisymmetric domains with conical vertices, whereas, in the third part, implementation issues, numerical tests and comparisons with other methods are carried out. The method is based on a Fourier expansion in the direction parallel to the reentrant edges of the domain, and on an improved variant of the Singular Complement Method in the 2D section perpendicular to those edges. Neither refinements near the reentrant edges or vertices of the domain, nor cut-off functions are required in the computations to achieve an optimal convergence order in terms of the mesh size and the number of Fourier modes used.     

Finite element approximation of the elasticity spectral problem on curved domains

We analyze the finite element approximation of the spectral problem for the linear elasticity equation with mixed boundary conditions on a curved non-convex domain. In the framework of the abstract spectral approximation theory, we obtain optimal order error estimates for the approximation of eigenvalues and eigenvectors. Two kinds of problems are considered: the discrete domain does not coincide with the real one and mixed boundary conditions are imposed. Some numerical results are presented.     

Focusing problem and the asymptotics of the spectral function of the Laplace-Beltrami operator. II

In the paper a formal high-frequency solution of the problem of a point source of oscillations near a reflecting boundary is constructed. The boundary is geodesically concave which makes it possible to use methods developed in diffraction theory by V. A. Fock and J. B. Keller. On the basis of the formal solution of the problem of a point source of oscillations, it is possible to construct the asymptotics of the spectral functions of the Laplace-Beltrami operator.Translated fromZapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 89, pp. 14–53, 1979.     

Regularity of a vector potential problem and its spectral curve

In this note we study a minimization problem for a vector of measures subject to a prescribed interaction matrix in the presence of external potentials. The conductors are allowed to have zero distance from each other but the external potentials satisfy a growth condition near the common points.We then specialize the setting to a specific problem on the real line which arises in the study of certain biorthogonal polynomials (studied elsewhere) and we prove that the equilibrium measures solve a pseudo-algebraic curve under the assumption that the potentials are real analytic. In particular, the supports of the equilibrium measures are shown to consist of a finite union of compact intervals.     

A spectral method for the Stokes problem in three-dimensional unbounded domains

We present a method for solving the Stokes problem in unbounded domains. It relies on the coupling of the transparent boundary operator and a spectral method in spherical coordinates. It is done explicitly by the use of vector-valued spherical harmonics. A uniform inf-sup condition is proved, which provides an optimal error estimate.

     

Dynamics in dumbbell domains II. The limiting problem

In this work we continue the analysis of the asymptotic dynamics of reaction-diffusion problems in a dumbbell domain started in [J.M. Arrieta, A.N. Carvalho, G. Lozada-Cruz, Dynamics in dumbbell domains I. Continuity of the set of equilibria, J. Differential Equations 231 (2) (2006) 551-597]. Here we study the limiting problem, that is, an evolution problem in a “domain” which consists of an open, bounded and smooth set Ω⊂RN with a curve R₀ attached to it. The evolution in both parts of the domain is governed by a parabolic equation. In Ω the evolution is independent of the evolution in R₀ whereas in R₀ the evolution depends on the evolution in Ω through the continuity condition of the solution at the junction points. We analyze in detail the linear elliptic and parabolic problem, the generation of linear and nonlinear semigroups, the existence and structure of attractors.     

The pricing problem. Part II: Computational complexity

In this paper, we show that the problem under study is in the class Log-APX, cannot be approximate with absolute error bounded by a constant, and the associated evaluation problem is nontrivial in the class Δ ₂ ^p . The two cases of the problem solvable in polynomial time are provided.     

The problem concerning the epimorphicity of a convolution operator in convex domains

We deduce an epimorphicity criterion for the convolution operator $$(a * x) (z) = \frac{1}{{2\pi i}}\oint {x (t) \tilde a (t - z) dt} ,$$ acting from a space of functions analytic in a convex domain into another such space;a(z) is the Borel transformation of the exponential functiona(z).     

A potential function method for the Cauchy problem of elliptic operators

In this paper, we deal with the Cauchy problem of elliptic operators. Through the use of a single-layer potential function, we devise a numerical method for approximating the solution of the Cauchy problem of elliptic operators, which are well known to be highly ill-posed in nature. The method is based on the denseness of single-layer potential functions. Convergence and stability estimates are then given with some examples for numerical verification on the efficiency of the proposed method. It has been observed that the use of more Cauchy data will greatly improve the accuracy of the approximate solutions.     

A boundary element method for the Dirichlet eigenvalue problem of the Laplace operator

The solution of eigenvalue problems for partial differential operators by using boundary integral equation methods usually involves some Newton potentials which may be resolved by using a multiple reciprocity approach. Here we propose an alternative approach which is in some sense equivalent to the above. Instead of a linear eigenvalue problem for the partial differential operator we consider a nonlinear eigenvalue problem for an associated boundary integral operator. This nonlinear eigenvalue problem can be solved by using some appropriate iterative scheme, here we will consider a Newton scheme. We will discuss the convergence and the boundary element discretization of this algorithm, and give some numerical results.     

A symmetric boundary element method for the Stokes problem in multiple connected domains

We consider a symmetric Galerkin boundary element method for the Stokes problem with general boundary conditions including slip conditions. The boundary value problem is reformulated as Steklov–Poincaré boundary integral equation which is then solved by a standard approximation scheme. An essential tool in our approach is the invertibility of the single layer potential which requires the definition of appropriate factor spaces due to the topology of the domain. Here we describe a modified boundary element approach to solve Dirichlet boundary value problems in multiple connected domains. A suitable extension of the standard single layer potential leads to an operator which is elliptic on the original function space. Copyright © 2003 John Wiley & Sons, Ltd.     

The Fourier Singular Complement Method for the Poisson problem. Part I: prismatic domains

This is the first part of a threefold article, aimed at solving numerically the Poisson problem in three-dimensional prismatic or axisymmetric domains. In this first part, the Fourier Singular Complement Method is introduced and analysed, in prismatic domains. In the second part, the FSCM is studied in axisymmetric domains with conical vertices, whereas, in the third part, implementation issues, numerical tests and comparisons with other methods are carried out. The method is based on a Fourier expansion in the direction parallel to the reentrant edges of the domain, and on an improved variant of the Singular Complement Method in the 2D section perpendicular to those edges. Neither refinements near the reentrant edges of the domain nor cut-off functions are required in the computations to achieve an optimal convergence order in terms of the mesh size and the number of Fourier modes used. This author was supported in part by the France/Hong Kong Joint Research Scheme. This author was supported by DGA/DSP-ENSTA 00.60.075.00.470.75.01 Research Programme This author was fully supported by Hong Kong RGC grants (Project CUHK4048/02P and project 403403).     

Scheme of the finite element method with multiplicative separation of the singularity for a spectral boundary value problem for a degenerate differential operator

The paper deals with the numerical solution of a generalized spectral boundary value problem for an elliptic operator with degenerating coefficients. We suggest an approximate method based on the multiplicative separation of the singularity, whereby the eigenfunctions are approximated by piecewise linear functions multiplied by a weight specially chosen depending on the order of degeneration of the coefficients. For this method, we obtain error estimates justifying its optimality.