Related papers: A fast direct solver for two dimensional quasi-per…
For a model convection-diffusion problem, we address the presence of oscillatory discrete solutions, and study difficulties in recovering standard approximation results for its solution. We justify the presence of non-physical oscillations…
This paper aims to solve numerically the two-dimensional inverse medium scattering problem with far-field data. This is a challenging task due to the severe ill-posedness and strong nonlinearity of the inverse problem. As already known, it…
In this paper, we study the inverse medium scattering problem to reconstruct unknown inhomogeneous medium from far field patterns of scattered waves. In the first part of our work, the linear inverse scattering problem was discussed, while…
In a previous work the authors described a fast high-fidelity computer model for acoustic scattering from multi-layered elastic spheres. This work is now extended with a scaling strategy significantly mitigating the problem of overflow and…
The paper describes a sparse direct solver for the linear systems that arise from the discretization of an elliptic PDE on a two dimensional domain. The scheme decomposes the domain into thin subdomains, or ``slabs'' and uses a two-level…
Ultrafast scattering using X-rays or electrons is an emerging method to obtain structure dynamics at the atomic length and time scales. However, directly resolving in real-space atomic motions is inherently limited by the finite detector…
Fast and high-order accurate algorithms for three dimensional elastic scattering are of great importance when modeling physical phenomena in mechanics, seismic imaging, and many other fields of applied science. In this paper, we develop a…
We propose a second order, fully semi-Lagrangian method for the numerical solution of systems of advection-diffusion-reaction equations, which employs a semi-Lagrangian approach to approximate in time both the advective and the diffusive…
We focus here on a class of fourth-order parabolic equations that can be written as a system of second-order equations by introducing an auxiliary variable. We design a novel second-order fully discrete mixed finite element method to…
We present a fast algorithm for evaluating the (non-smooth) solution of the free-space two-dimensional (2D) scalar wave equation with many point sources, each with a high-frequency band-limited time signature. Such an algorithm is key to an…
This article is a follow up of our submitted paper [11] in which a decomposition of the Richards equation along two soil layers was discussed. A decomposed problem was formulated and a decoupling and linearisation technique was presented to…
This work develops a nonlinear multigrid method for diffusion problems discretized by cell-centered finite volume methods on general unstructured grids. The multigrid hierarchy is constructed algebraically using aggregation of degrees of…
The computational cost of fluid simulations increases rapidly with grid resolution. This has given a hard limit on the ability of simulations to accurately resolve small scale features of complex flows. Here we use a machine learning…
We present a fast direct solution method for the advection-diffusion equation in one and two dimensions with non-periodic boundaries. Computational cost is reduced to $\mathcal O(N)$ by making a low-rank approximation of the Green's…
We present a class of new explicit and stable numerical algorithms to solve the spatially discretized linear heat or diffusion equation. After discretizing the space and the time variables like conventional finite difference methods, we do…
Although high-order Maxwell integral equation solvers provide significant advantages in terms of speed and accuracy over corresponding low-order integral methods, their performance significantly degrades in presence of non-smooth…
We consider the direct electromagnetic scattering problem of time-harmonic obliquely incident waves by a infinitely long, homogeneous and doubly-connected cylinder in three dimensions. We apply a hybrid integral equation method (combination…
Predicting the optical response of macroscopic arrangements of individual scatterers is a computational challenge, as the problem involves length scales across multiple orders of magnitude. We present a full-wave optical method to highly…
Multiple light scattering hampers imaging objects in complex scattering media. Approaches used in real practices mainly aim to filter out multiple scattering obscuring the ballistic waves that travel straight through the scattering medium.…
We study two-dimensional wave propagation in materials whose properties vary periodically in one direction only. High order homogenization is carried out to derive a dispersive effective medium approximation. One-dimensional materials with…