Related papers: Efficient analytical implementation of the DOT Rie…
The stability and convergence analysis of high-order numerical approximations for the one- and two-dimensional nonlocal wave equations on unbounded spatial domains are considered. We first use the quadrature-based finite difference schemes…
We consider partitioned time integration for heterogeneous coupled heat equations. First and second order multirate, as well as time-adaptive Dirichlet-Neumann Waveform relaxation (DNWR) methods are derived. In 1D and for implicit Euler…
Consider the elastic scattering of an incident wave by a rigid obstacle in three dimensions, which is formulated as an exterior problem for the Navier equation. By constructing a Dirichlet-to-Neumann (DtN) operator and introducing a…
We introduce a new class of integrators for stiff ODEs as well as SDEs. These integrators are (i) {\it Multiscale}: they are based on flow averaging and so do not fully resolve the fast variables and have a computational cost determined by…
Resolvent analysis identifies the most responsive forcings and most receptive states of a dynamical system, in an input--output sense, based on its governing equations. Interest in the method has continued to grow during the past decade due…
We study the Stokes--Poisson--Boltzmann equations with Dirichlet and Navier boundary conditions. The system consists of the incompressible Stokes equations coupled with a nonlinear Poisson--Boltzmann equation through electrostatic forcing…
This paper develops the numerical inverse scattering transform (NIST) framework for the coupled modified Korteweg-de Vries (mKdV) equation based on its associated Riemann-Hilbert problem. The coupled system gives rise to a $3\times3$…
We evaluate the performance of novel numerical methods for solving one-dimensional nonlinear fractional dispersive and dissipative evolution equations. The methods are based on affine combinations of time-splitting integrators and…
A new implementation of the Roe scheme for solving two-layer shallow-water equations is presented in this paper. The proposed A-Roe scheme is based on the analytical solution to the characteristic quartic of the flux matrix, which is an…
We consider time dependent thermal fluid structure interaction. The respective models are the compressible Navier-Stokes equations and the nonlinear heat equation. A partitioned coupling approach via a Dirichlet-Neumann method and a fixed…
Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in…
A new Active Flux method for the multi-dimensional Euler equations is based on an additive operator splitting into acoustics and advection. The acoustic operator is solved in a locally linearized manner by using the exact evolution…
This work aims to construct an efficient and highly accurate numerical method to address the time singularity at $t=0$ involved in a class of time-fractional parabolic integro-partial differential equations in one and two dimensions. The…
Semi-discrete optimal transport (SOT), which maps a continuous probability measure to a discrete one, is a fundamental problem with wide-ranging applications. Entropic regularization is often employed to solve the SOT problem, leading to a…
Nonlinear parametric inverse problems appear in many applications. Here, we focus on diffuse optical tomography (DOT) in medical imaging to recover unknown images of interest, such as cancerous tissue in a given medium, using a mathematical…
We develop robust and scalable fully implicit nonlinear finite element solvers for the simulations of biological transportation networks driven by the gradient flow minimization of a non-convex energy cost functional. Our approach employs a…
We present the combination of a complex-time tensor-network impurity solver with an analytic continuation scheme based on exponential fitting as an efficient framework for single and multi-orbital dynamical mean-field calculations. By…
The Diffuse Optical Tomography (DOT) has received considerable attention in the recent years in the field of biomedical imaging and disease detection. However, imaging through highly diffusive medium is a challenge and stability is always…
Two reconstruction methods of Electrical Impedance Tomography (EIT) are numerically compared for nonsmooth conductivities in the plane based on the use of complex geometrical optics (CGO) solutions to D-bar equations involving the global…
The Riemann problem for first-order hyperbolic systems of partial differential equations is of fundamental importance for both theoretical and numerical purposes. Many approximate solvers have been developed for such systems; exact solution…