BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Wim Vanroose

Wim Vanroose contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.comp-phmath-phmath.DSmath.MPComputational Engineering, Finance, and Sciencemath.APquant-phBiological Physicscond-mat.supr-conmath.OCnucl-thphysics.atom-phphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
15topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

21 published item(s)

preprint2022arXiv

Solving for the low-rank tensor components of a scattering wave function

Atomic and molecular breakup reactions, such as multiple-ionisation, are described by a driven Schrödinger equation. This equation is equivalent to a high-dimensional Helmholtz equation and it has solutions that are outgoing waves, emerging from the target. We show that these waves can be described by a low-rank approximation. For 2D problems this it a matrix product of two low-rank matrices, for 3D problems it is a low-rank tensor decomposition. We propose an iterative method that solves, in an alternating way, for these low-rank components of the scattered wave. We illustrate the method with examples in 2D and 3D.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Krylov-Simplex method that minimizes the residual in $\ell_1$-norm or $\ell_\infty$-norm

The paper presents two variants of a Krylov-Simplex iterative method that combines Krylov and simplex iterations to minimize the residual $r = b-Ax$. The first method minimizes $\|r\|_\infty$, i.e. maximum of the absolute residuals. The second minimizes $\|r\|_1$, and finds the solution with the least absolute residuals. Both methods search for an optimal solution $x_k$ in a Krylov subspace which results in a small linear programming problem. A specialized simplex algorithm solves this projected problem and finds the optimal linear combination of Krylov basis vectors to approximate the solution. The resulting simplex algorithm requires the solution of a series of small dense linear systems that only differ by rank-one updates. The $QR$ factorization of these matrices is updated each iteration. We demonstrate the effectiveness of the methods with numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Projected Newton Method for noise constrained Tikhonov regularization

Tikhonov regularization is a popular approach to obtain a meaningful solution for ill-conditioned linear least squares problems. A relatively simple way of choosing a good regularization parameter is given by Morozov's discrepancy principle. However, most approaches require the solution of the Tikhonov problem for many different values of the regularization parameter, which is computationally demanding for large scale problems. We propose a new and efficient algorithm which simultaneously solves the Tikhonov problem and finds the corresponding regularization parameter such that the discrepancy principle is satisfied. We achieve this by formulating the problem as a nonlinear system of equations and solving this system using a line search method. We obtain a good search direction by projecting the problem onto a low dimensional Krylov subspace and computing the Newton direction for the projected problem. This projected Newton direction, which is significantly less computationally expensive to calculate than the true Newton direction, is then combined with a backtracking line search to obtain a globally convergent algorithm, which we refer to as the Projected Newton method. We prove convergence of the algorithm and illustrate the improved performance over current state-of-the-art solvers with some numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A multi-level preconditioned Krylov method for the efficient solution of algebraic tomographic reconstruction problems

Classical iterative methods for tomographic reconstruction include the class of Algebraic Reconstruction Techniques (ART). Convergence of these stationary linear iterative methods is however notably slow. In this paper we propose the use of Krylov solvers for tomographic linear inversion problems. These advanced iterative methods feature fast convergence at the expense of a higher computational cost per iteration, causing them to be generally uncompetitive without the inclusion of a suitable preconditioner. Combining elements from standard multigrid (MG) solvers and the theory of wavelets, a novel wavelet-based multi-level (WMG) preconditioner is introduced, which is shown to significantly speed-up Krylov convergence. The performance of the WMG-preconditioned Krylov method is analyzed through a spectral analysis, and the approach is compared to existing methods like the classical Simultaneous Iterative Reconstruction Technique (SIRT) and unpreconditioned Krylov methods on a 2D tomographic benchmark problem. Numerical experiments are promising, showing the method to be competitive with the classical Algebraic Reconstruction Techniques in terms of convergence speed and overall performance (CPU time) as well as precision of the reconstruction.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An efficient multigrid calculation of the far field map for Helmholtz and Schrödinger equations

In this paper we present a new highly efficient calculation method for the far field amplitude pattern that arises from scattering problems governed by the d-dimensional Helmholtz equation and, by extension, Schrödinger's equation. The new technique is based upon a reformulation of the classical real-valued Green's function integral for the far field amplitude to an equivalent integral over a complex domain. It is shown that the scattered wave, which is essential for the calculation of the far field integral, can be computed very efficiently along this complex contour (or manifold, in multiple dimensions). Using the iterative multigrid method as a solver for the discretized damped scattered wave system, the proposed approach results in a fast and scalable calculation method for the far field map. The complex contour method is successfully validated on Helmholtz and Schrödinger model problems in two and three spatial dimensions, and multigrid convergence results are provided to substantiate the wavenumber scalability and overall performance of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constrained Runs algorithm as a lifting operator for the Boltzmann equation

Lifting operators play an important role in starting a kinetic Boltzmann model from given macroscopic information. The macroscopic variables need to be mapped to the distribution functions, mesoscopic variables of the Boltzmann model. A well-known numerical method for the initialization of Boltzmann models is the Constrained Runs algorithm. This algorithm is used in literature for the initialization of lattice Boltzmann models, special discretizations of the Boltzmann equation. It is based on the attraction of the dynamics toward the slow manifold and uses lattice Boltzmann steps to converge to the desired dynamics on the slow manifold. We focus on applying the Constrained Runs algorithm to map density, average flow velocity, and temperature, the macroscopic variables, to distribution functions. Furthermore, we do not consider only lattice Boltzmann models. We want to perform the algorithm for different discretizations of the Boltzmann equation and consider a standard finite volume discretization.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Initialization of lattice Boltzmann models with the help of the numerical Chapman-Enskog expansion

We extend the applicability of the numerical Chapman-Enskog expansion as a lifting operator for lattice Boltzmann models to map density and momentum to distribution functions. In earlier work [Vanderhoydonc et al. Multiscale Model. Simul. 10(3): 766-791, 2012] such an expansion was constructed in the context of lifting only the zeroth order velocity moment, namely the density. A lifting operator is necessary to convert information from the macroscopic to the mesoscopic scale. This operator is used for the initialization of lattice Boltzmann models. Given only density and momentum, the goal is to initialize the distribution functions of lattice Boltzmann models. For this initialization, the numerical Chapman-Enskog expansion is used in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A new level-dependent coarsegrid correction scheme for indefinite Helmholtz problems

In this paper we construct and analyse a level-dependent coarsegrid correction scheme for indefinite Helmholtz problems. This adapted multigrid method is capable of solving the Helmholtz equation on the finest grid using a series of multigrid cycles with a grid-dependent complex shift, leading to a stable correction scheme on all levels. It is rigourously shown that the adaptation of the complex shift throughout the multigrid cycle maintains the functionality of the two-grid correction scheme, as no smooth modes are amplified in or added to the error. In addition, a sufficiently smoothing relaxation scheme should be applied to ensure damping of the oscillatory error components. Numerical experiments on various benchmark problems show the method to be competitive with or even outperform the current state-of-the-art multigrid-preconditioned Krylov methods, like e.g. CSL-preconditioned GMRES or BiCGStab.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

GMRES-based multigrid for the complex scaled preconditoner for the indefinite Helmholtz equation

Multigrid preconditioners and solvers for the indefinite Helmholtz equation suffer from non-stability of the stationary smoothers due to the indefinite spectrum of the operator. In this paper we explore GMRES as a replacement for the stationary smoothers of the standard multigrid method. This results in a robust and efficient solver for a complex shifted or stretched Helmholtz problem that can be used as a preconditioner. Very few GMRES iterations are required on each level to build a good multigrid method. The convergence behavior is compared to a theoretically derived stable polynomial smoother. We test this method on some benchmark problems and report on the observed convergence behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Local Fourier Analysis of the Complex Shifted Laplacian preconditioner for Helmholtz problems

In this paper we solve the Helmholtz equation with multigrid preconditioned Krylov subspace methods. The class of Shifted Laplacian preconditioners are known to significantly speed-up Krylov convergence. However, these preconditioners have a parameter beta, a measure of the complex shift. Due to contradictory requirements for the multigrid and Krylov convergence, the choice of this shift parameter can be a bottleneck in applying the method. In this paper, we propose a wavenumber-dependent minimal complex shift parameter which is predicted by a rigorous k-grid Local Fourier Analysis (LFA) of the multigrid scheme. We claim that, given any (regionally constant) wavenumber, this minimal complex shift parameter provides the reader with a parameter choice that leads to efficient Krylov convergence. Numerical experiments in one and two spatial dimensions validate the theoretical results. It appears that the proposed complex shift is both the minimal requirement for a multigrid V-cycle to converge, as well as being near-optimal in terms of Krylov iteration count.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An optimal linear solver for the Jacobian system of the extreme type-II Ginzburg--Landau problem

This paper considers the extreme type-II Ginzburg--Landau equations, a nonlinear PDE model for describing the states of a wide range of superconductors. Based on properties of the Jacobian operator and an AMG strategy, a preconditioned Newton--Krylov method is constructed. After a finite-volume-type discretization, numerical experiments are done for representative two- and three-dimensional domains. Strong numerical evidence is provided that the number of Krylov iterations is independent of the dimension $n$ of the solution space, yielding an overall solver complexity of O(n).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Efficient determination of the energy landscape of nonlinear Schrödinger-type equations

We describe a systematic approach for the efficient numerical solution of nonlinear Schrödinger-type partial differential equations of the form $(K +V + g|ψ|^2)ψ=0$, with an energy operator $K$, a scalar potential $V$, and a scalar parameter $g$. Instrumental to the approach are developments in numerical linear and nonlinear algebra, specifically numerical parameter continuation. We demonstrate how a continuous sequence of solutions can be obtained regardless of their stability, so that finally the spectrum of stable and unstable solutions in the specified parameter range is fully revealed. The method is demonstrated for the GL equation in a three-dimensional superconducting domain with an inhomogeneous magnetic field, a numerically demanding problem known to have an involved solution landscape.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Improved convergence of scattering calculations in the oscillator representation

The Schrödinger equation for two and tree-body problems is solved for scattering states in a hybrid representation where solutions are expanded in the eigenstates of the harmonic oscillator in the interaction region and on a finite difference grid in the near-- and far--field. The two representations are coupled through a high--order asymptotic formula that takes into account the function values and the third derivative in the classical turning points. For various examples the convergence is analyzed for various physics problems that use an expansion in a large number of oscillator states. The results show significant improvement over the JM-ECS method [Bidasyuk et al, Phys. Rev. C 82, 064603 (2010)].

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical bifurcation analysis of pattern formation in a cell based auxin transport model

Transport models of growth hormones can be used to reproduce the hormone accumulations that occur in plant organs. Mostly, these accumulation patterns are calculated using time step methods, even though only the resulting steady state patterns of the model are of interest. We examine the steady state solutions of the hormone transport model of Smith et al (2006) for a one-dimensional row of plant cells. We search for the steady state solutions as a function of three of the model parameters by using numerical continuation methods and bifurcation analysis. These methods are more adequate for solving steady state problems than time step methods. We discuss a trivial solution where the concentrations of hormones are equal in all cells and examine its stability region. We identify two generic bifurcation scenarios through which the trivial solution loses its stability. The trivial solution becomes either a steady state pattern with regular spaced peaks or a pattern where the concentration is periodic in time.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical extraction of a macroscopic pde and a lifting operator from a lattice Boltzmann model

Lifting operators play an important role in starting a lattice Boltzmann model from a given initial density. The density, a macroscopic variable, needs to be mapped to the distribution functions, mesoscopic variables, of the lattice Boltzmann model. Several methods proposed as lifting operators have been tested and discussed in the literature. The most famous methods are an analytically found lifting operator, like the Chapman-Enskog expansion, and a numerical method, like the Constrained Runs algorithm, to arrive at an implicit expression for the unknown distribution functions with the help of the density. This paper proposes a lifting operator that alleviates several drawbacks of these existing methods. In particular, we focus on the computational expense and the analytical work that needs to be done. The proposed lifting operator, a numerical Chapman-Enskog expansion, obtains the coefficients of the Chapman-Enskog expansion numerically. Another important feature of the use of lifting operators is found in hybrid models. There the lattice Boltzmann model is spatially coupled with a model based on a more macroscopic description, for example an advection-diffusion-reaction equation. In one part of the domain, the lattice Boltzmann model is used, while in another part, the more macroscopic model. Such a hybrid coupling results in missing data at the interfaces between the different models. A lifting operator is then an important tool since the lattice Boltzmann model is typically described by more variables than a model based on a macroscopic partial differential equation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A preconditioned iterative solver for the scattering solutions of the Schrödinger equation

The Schrödinger equation defines the dynamics of quantum particles which has been an area of unabated interest in physics. We demonstrate how simple transformations of the Schrödinger equation leads to a coupled linear system, whereby each diagonal block is a high frequency Helmholtz problem. Based on this model, we derive indefinite Helmholtz model problems with strongly varying wavenumbers. We employ the iterative approach for their solution. In particular, we develop a preconditioner that has its spectrum restricted to a quadrant (of the complex plane) thereby making it easily invertible by multigrid methods with standard components. This multigrid preconditioner is used in conjuction with suitable Krylov-subspace methods for solving the indefinite Helmholtz model problems. The aim of this study is to report the feasbility of this preconditioner for the model problems. We compare this idea with the other prevalent preconditioning ideas, and discuss its merits. Results of numerical experiments are presented, which complement the proposed ideas, and show that this preconditioner may be used in an automatic setting.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analyzing the wave number dependency of the convergence rate of a multigrid preconditioned Krylov method for the Helmholtz equation with an absorbing layer

This paper analyzes the Krylov convergence rate of a Helmholtz problem preconditioned with Multigrid. The multigrid method is applied to the Helmholtz problem formulated on a complex contour and uses GMRES as a smoother substitute at each level. A one-dimensional model is analyzed both in a continuous and discrete way. It is shown that the Krylov convergence rate of the continuous problem is independent of the wave number. The discrete problem, however, can deviate significantly from this bound due to a pitchfork in the spectrum. It is further shown in numerical experiments that the convergence rate of the Krylov method approaches the continuous bound as the grid distance $h$ gets small.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical bifurcation study of superconducting patterns on a square

This paper considers the extreme type-II Ginzburg-Landau equations that model vortex patterns in superconductors. The nonlinear PDEs are solved using Newton's method, and properties of the Jacobian operator are highlighted. Specifically, it is illustrated how the operator can be regularized using an appropriate phase condition. For a two-dimensional square sample, the numerical results are based on a finite-difference discretization with link variables that preserves the gauge invariance. For two exemplary sample sizes, a thorough bifurcation analysis is performed using the strength of the applied magnetic field as a bifurcation parameter and focusing on the symmetries of this system. The analysis gives new insight in the transitions between stable and unstable states, as well as the connections between stable solution branches.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Applying numerical continuation to the parameter dependence of solutions of the Schrödinger equation

In molecular reactions at the microscopic level the appearance of resonances has an important influence on the reactivity. It is important to predict when a bound state transitions into a resonance and how these transitions depend on various system parameters such as internuclear distances. The dynamics of such systems are described by the time-independent Schrödinger equation and the resonances are modeled by poles of the S-matrix. Using numerical continuation methods and bifurcation theory, techniques which find their roots in the study of dynamical systems, we are able to develop efficient and robust methods to study the transitions of bound states into resonances. By applying Keller's Pseudo-Arclength continuation, we can minimize the numerical complexity of our algorithm. As continuation methods generally assume smooth and well-behaving functions and the S-matrix is neither, special care has been taken to ensure accurate results. We have successfully applied our approach in a number of model problems involving the radial Schrödinger equation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical Continuation of resonances and bound states in coupled channel Schrödinger equations

In this contribution, we introduce numerical continuation methods and bifurcation theory, techniques which find their roots in the study of dynamical systems, to the problem of tracing the parameter dependence of bound and resonant states of the quantum mechanical Schrödinger equation. We extend previous work on the subject to systems of coupled equations. Bound and resonant states of the Schrödinger equation can be determined through the poles of the S-matrix, a quantity that can be derived from the asymptotic form of the wave function. We introduce a regularization procedure that essentially transforms the S-matrix into its inverse and improves its smoothness properties, thus making it amenable to numerical continuation. This allows us to automate the process of tracking bound and resonant states when parameters in the Schrödinger equation are varied. We have applied this approach to a number of model problems with satisfying results.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the indefinite Helmholtz equation: complex stretched absorbing boundary layers, iterative analysis, and preconditioning

This paper studies and analyzes a preconditioned Krylov solver for Helmholtz problems that are formulated with absorbing boundary layers based on complex coordinate stretching. The preconditioner problem is a Helmholtz problem where not only the coordinates in the absorbing layer have an imaginary part, but also the coordinates in the interior region. This results into a preconditioner problem that is invertible with a multigrid cycle. We give a numerical analysis based on the eigenvalues and evaluate the performance with several numerical experiments. The method is an alternative to the complex shifted Laplacian and it gives a comparable performance for the studied model problems.

0 citations0 reviewsNo code linkedOpen access