BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Christiaan C. Stolk

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

ResearcherUnclaimed source record
math.NAmath.APmath-phmath.MPphysics.comp-ph

Catalog footprint

What is connected

8works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

8 published item(s)

preprint2016arXiv

A dispersion minimizing scheme for the 3-D Helmholtz equation based on ray theory

We develop a new dispersion minimizing compact finite difference scheme for the Helmholtz equation in 2 and 3 dimensions. The scheme is based on a newly developed ray theory for difference equations. A discrete Helmholtz operator and a discrete operator to be applied to the source and the wavefields are constructed. Their coefficients are piecewise polynomial functions of $hk$, chosen such that phase and amplitude errors are minimal. The phase errors of the scheme are very small, approximately as small as those of the 2-D quasi-stabilized FEM method and substantially smaller than those of alternatives in 3-D, assuming the same number of gridpoints per wavelength is used. In numerical experiments, accurate solutions are obtained in constant and smoothly varying media using meshes with only five to six points per wavelength and wave propagation over hundreds of wavelengths. When used as a coarse level discretization in a multigrid method the scheme can even be used with downto three points per wavelength. Tests on 3-D examples with up to $10^8$ degrees of freedom show that with a recently developed hybrid solver, the use of coarser meshes can lead to corresponding savings in computation time, resulting in good simulation times compared to the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An improved sweeping domain decomposition preconditioner for the Helmholtz equation

In this paper we generalize and improve a recently developed domain decomposition preconditioner for the iterative solution of discretized Helmholtz equations. We introduce an improved method for transmission at the internal boundaries using perfectly matched layers. Simultaneous forward and backward sweeps are introduced, thereby improving the possibilities for parallellization. Finally, the method is combined with an outer two-grid iteration. The method is studied theoretically and with numerical examples. It is shown that the modifications lead to substantial decreases in computation time and memory use, so that computation times become comparable to that of the fastests methods currently in the literature for problems with up to 10^8 degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A multigrid method for the Helmholtz equation with optimized coarse grid corrections

We study the convergence of multigrid schemes for the Helmholtz equation, focusing in particular on the choice of the coarse scale operators. Let G_c denote the number of points per wavelength at the coarse level. If the coarse scale solutions are to approximate the true solutions, then the oscillatory nature of the solutions implies the requirement G_c > 2. However, in examples the requirement is more like G_c >= 10, in a trade-off involving also the amount of damping present and the number of multigrid iterations. We conjecture that this is caused by the difference in phase speeds between the coarse and fine scale operators. Standard 5-point finite differences in 2-D are our first example. A new coarse scale 9-point operator is constructed to match the fine scale phase speeds. We then compare phase speeds and multigrid performance of standard schemes with a scheme using the new operator. The required G_c is reduced from about 10 to about 3.5, with less damping present so that waves propagate over > 100 wavelengths in the new scheme. Next we consider extensions of the method to more general cases. In 3-D comparable results are obtained with standard 7-point differences and optimized 27-point coarse grid operators, leading to an order of magnitude reduction in the number of unknowns for the coarsest scale linear system. Finally we show how to include PML boundary layers, using a regular grid finite element method. Matching coarse scale operators can easily be constructed for other discretizations. The method is therefore potentially useful for a large class of discretized high-frequency Helmholtz equations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A rapidly converging domain decomposition method for the Helmholtz equation

A new domain decomposition method is introduced for the heterogeneous 2-D and 3-D Helmholtz equations. Transmission conditions based on the perfectly matched layer (PML) are derived that avoid artificial reflections and match incoming and outgoing waves at the subdomain interfaces. We focus on a subdivision of the rectangular domain into many thin subdomains along one of the axes, in combination with a certain ordering for solving the subdomain problems and a GMRES outer iteration. When combined with multifrontal methods, the solver has near-linear cost in examples, due to very small iteration numbers that are essentially independent of problem size and number of subdomains. It is to our knowledge only the second method with this property next to the moving PML sweeping method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A mathematical framework for inverse wave problems in heterogeneous media

This paper provides a theoretical foundation for some common formulations of inverse problems in wave propagation, based on hyperbolic systems of linear integro-differential equations with bounded and measurable coefficients. The coefficients of these time-dependent partial differential equations respresent parametrically the spatially varying mechanical properties of materials. Rocks, manufactured materials, and other wave propagation environments often exhibit spatial heterogeneity in mechanical properties at a wide variety of scales, and coefficient functions representing these properties must mimic this heterogeneity. We show how to choose domains (classes of nonsmooth coefficient functions) and data definitions (traces of weak solutions) so that optimization formulations of inverse wave problems satisfy some of the prerequisites for application of Newton's method and its relatives. These results follow from the properties of a class of abstract first-order evolution systems, of which various physical wave systems appear as concrete instances. Finite speed of propagation for linear waves with bounded, measurable mechanical parameter fields is one of the by-products of this theory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Linearized inverse scattering based on seismic Reverse Time Migration

In this paper we study the linearized inverse problem associated with imaging of reflection seismic data. We introduce an inverse scattering transform derived from reverse-time migration (RTM). In the process, the explicit evaluation of the so-called normal operator is avoided, while other differential and pseudodifferential operator factors are introduced. We prove that, under certain conditions, the transform yields a partial inverse, and support this with numerical simulations. In addition, we explain the recently discussed 'low-frequency artifacts' in RTM, which are naturally removed by the new method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the generalization of wavelet diagonal preconditioning to the Helmholtz equation

We present a preconditioning method for the multi-dimensional Helmholtz equation with smoothly varying coefficient. The method is based on a frame of functions, that approximately separates components associated with different singular values of the operator. For the small singular values, corresponding to propagating waves, the frame functions are constructed using ray theory. A series of 2-D numerical experiments demonstrates that the number of iterations required for convergence is small and independent of the frequency. In this sense the method is optimal.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Preconditioners based on Windowed Fourier Frames Applied to Elliptic Partial Differential Equations

We investigate the application of windowed Fourier frames (WFFs) to the numerical solution of partial differential equations, focussing on elliptic equations. The action of a partial differential operator (PDO) on a windowed plane wave is close to a multiplication, where the multiplication factor is given by the symbol of the PDO evaluated at the wave number and central position of the windowed plane wave. This can be exploited in a preconditioning method for use in iterative inversion. For domains with periodic boundary conditions we find that the condition number with the preconditioning becomes bounded and the iteration converges well. For problems with a Dirichlet boundary condition, some large and small singular values remain. However the iterative inversion still appears to converge well.

0 citations0 reviewsNo code linkedOpen access