Source author record

Laurence Halpern

Laurence Halpern appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA math.AP Numerical Analysis

Catalog footprint

What is connected

6works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Perfectly Matched Layers on Cubic Domains forPauli's Equations

This article proves the well posedness of the boundary value problemthat arises when PML algorithms are applied to Pauli's equationswith a three dimensional rectangle as computational domain. The absorptionsare positive near the boundary and zero far from the boundary so are always x-dependent. At the flat parts of the boundary of the rectangle, the natural absorbing boundary conditions are imposed.The difficulty addressed is the analysis of the resulting variable coeffi-cient problem on the rectanglar solid with its edges and corners. TheLaplace transform is analysed. It turns on the analysis of a boundaryvalue problem formally obtained by complex stretching. Existence isproved by deriving a boundary value problems for a complex stretchedHelmholtz equation on smoothed domains. This is the first stabilityproof with x-dependent absorptions on a bounded domain whoseboundary is not smooth.

preprint2014arXiv

A new Algorithm Based on Factorization for Heterogeneous Domain Decomposition

Often computational models are too expensive to be solved in the entire domain of simulation, and a cheaper model would suffice away from the main zone of interest. We present for the concrete example of an evolution problem of advection reaction diffusion type a heterogeneous domain decomposition algorithm which allows us to recover a solution that is very close to the solution of the fully viscous problem, but solves only an inviscid problem in parts of the domain. Our new algorithm is based on the factorization of the underlying differential operator, and we therefore call it factorization algorithm. We give a detailed error analysis, and show that we can obtain approximations in the viscous region which are much closer to the viscous solution in the entire domain of simulation than approximations obtained by other heterogeneous domain decomposition algorithms from the literature.

preprint2014arXiv

Optimized Schwarz Waveform Relaxation for Advection Reaction Diffusion Equations in Two Dimensions

Optimized Schwarz Waveform Relaxation methods have been developed over the last decade for the parallel solution of evolution problems. They are based on a decomposition in space and an iteration, where only subproblems in space-time need to be solved. Each subproblem can be simulated using an adapted numerical method, for example with local time stepping, or one can even use a different model in different subdomains, which makes these methods very suitable also from a modeling point of view. For rapid convergence however, it is important to use effective transmission conditions between the space-time subdomains, and for best performance, these transmission conditions need to take the physics of the underlying evolution problem into account. The optimization of these transmission conditions leads to a mathematically hard best approximation problem of homographic type. We study in this paper in detail this problem for the case of linear advection reaction diffusion equations in two spatial dimensions. We prove comprehensively best approximation results for transmission conditions of Robin and Ventcel type. We give for each case closed form asymptotic values for the parameters, which guarantee asymptotically best performance of the iterative methods. We finally show extensive numerical experiments, and we measure performance corresponding to our analysis

preprint2011arXiv

The Analysis of Matched Layers

A systematic analysis of matched layers is undertaken with special attention to better understand the remarkable method of Bérenger. We prove that the Bérenger and closely related layers define well posed transmission problems in great generality. When the Bérenger method or one of its close relatives is well posed, perfect matching is proved. The proofs use the energy method, Fourier-Laplace transform, and real coordinate changes for Laplace transformed equations. It is proved that the loss of derivatives associated with the Bérenger method does not occur for elliptic generators. More generally, an essentially necessary and sufficient condition for loss of derivatives in Bérenger's method is proved. The sufficiency relies on the energy method with pseudodifferential multiplier. Amplifying and nonamplifying layers are identified by a geometric optics computation. Among the various flavors of Bérenger's algorithm for Maxwell's equations our favorite choice leads to a strongly well posed augmented system and is both perfect and nonamplifying in great generality. We construct by an extrapolation argument an alternative matched layer method which preserves the strong hyperbolicity of the original problem and though not perfectly matched has leading reflection coefficient equal to zero at all angles of incidence.

preprint2010arXiv

Optimized Schwarz waveform relaxation and discontinuous Galerkin time stepping for heterogeneous problems

We design and analyze a Schwarz waveform relaxation algorithm for domain decomposition of advection-diffusion-reaction problems with strong heterogeneities. The interfaces are curved, and we use optimized Robin or Ventcell transmission conditions. We analyze the semi-discretization in time with Discontinuous Galerkin as well. We also show two-dimensional numerical results using generalized mortar finite elements in space.

preprint2007arXiv

Nonlinear Nonoverlapping Schwarz Waveform Relaxation for Semilinear Wave Propagation

We introduce a non-overlapping variant of the Schwarz waveform relaxation algorithm for semilinear wave propagation in one dimension. Using the theory of absorbing boundary conditions, we derive a new nonlinear algorithm. We show that the algorithm is well-posed and we prove its convergence by energy estimates and a Galerkin method. We then introduce an explicit scheme. We prove the convergence of the discrete algorithm with suitable assumptions on the nonlinearity. We finally illustrate our analysis with numerical experiments.