Related papers: A fast and robust computational method for the ion…
We introduce a simple and stable computational method for ill-posed partial differential equation (PDE) problems. The method is based on Schr\"odingerization, introduced in [S. Jin, N. Liu and Y. Yu, arXiv:2212.13969][S. Jin, N. Liu and Y.…
We propose an iterative solution method for the 3D high-frequency Helmholtz equation that exploits a contour integral formulation of spectral projectors. In this framework, the solution in certain invariant subspaces is approximated by…
We implement the Numerical Unified Transform Method to solve the Nonlinear Schr\"odinger equation on the half-line. For so-called linearizable boundary conditions, the method solves the half-line problems with comparable complexity as the…
This paper presents a linear, decoupled, mass- and energy-conserving numerical scheme for the multi-dimensional coupled nonlinear Schr\"odinger (CNLS) system. The scheme combines the fourth-order compact difference approximation in space…
Atomic and molecular breakup reactions, such as multiple-ionisation, are described by a driven Schr\"odinger equation. This equation is equivalent to a high-dimensional Helmholtz equation and it has solutions that are outgoing waves,…
This paper is to give an overview of AMG methods for solving large scale systems of equations such as those from the discretization of partial differential equations. AMG is often understood as the acronym of "Algebraic Multi-Grid", but it…
We present a numerical formulation for the solution of non-isothermal, compressible, Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical…
The nonlinear Schr\"{o}dinger (NLS) equation possesses an infinite hierarchy of conserved densities and the numerical preservation of some of these quantities is critical for accurate long-time simulations, particularly for multi-soliton…
We present an approach to numerically solving the time-dependent Schroedinger equation and other parabolic equations by the split-step technique with fast Fourier transform, which suppresses the backreflection of waves from the grid…
We study the suppression of reflections in the numerical simulation of the time-dependent Schr\"odinger equation for strong-field problems on a grid using exterior complex scaling (ECS) as absorbing boundary condition. It is shown that the…
This paper discusses the multiscale approach and the convergence of the time-dependent Maxwell-Schr\"{o}dinger system with rapidly oscillating discontinuous coefficients arising from the modeling of a heterogeneous nanostructure with a…
Spectral deferred correction (SDC) methods are an attractive approach to iteratively computing collocation solutions to an ODE by performing so-called sweeps with a low-order time stepping method. SDC allows to easily construct high order…
A new iterative method for non-LTE multilevel polarized radiative transfer in hydrogen lines is presented. Iterative methods (such as the Jacobi method) tend to damp out high-frequency components of the error fast, but converges poorly due…
We present a hybrid method for time-dependent particle transport that combines Monte Carlo (MC) estimation with a deterministic discrete ordinates (\(S_N\)) solve, augmented by quasi-Monte Carlo (QMC) sampling. For spatial discretizations,…
One might argue that solving a trajectory optimization problem over a million grid points is preposterous. How about solving such a problem at an incredibly fast computational time? On a small form-factor processor? Algorithmic details that…
A review of a recent method is presented to construct certain exact solutions to the focusing nonlinear Schr\"odinger equation on the line with a cubic nonlinearity. With motivation by the inverse scattering transform and help from the…
In this work, we propose an efficient and robust multigrid method for solving the time-fractional heat equation. Due to the nonlocal property of fractional differential operators, numerical methods usually generate systems of equations for…
This paper introduces a second-order method for solving general elliptic partial differential equations (PDEs) on irregular domains using GPU acceleration, based on Ying's kernel-free boundary integral (KFBI) method. The method addresses…
We introduce a numerical solver for the steady-state Boltzmann equation based on the symmetric Gauss-Seidel (SGS) method. To solve the nonlinear system on each grid cell derived from the SGS method, a fixed-point iteration preconditioned…
In this paper, we present a data-driven approach to iteratively solve the discrete heterogeneous Helmholtz equation at high wavenumbers. In our approach, we combine classical iterative solvers with convolutional neural networks (CNNs) to…