Related papers: Implicit Asymptotic Preserving Method for Linear T…
An operator-splitting finite element scheme for the time-dependent, high-dimensional radiative transfer equation is presented in this paper. The streamline upwind Petrov-Galerkin finite element method and discontinuous Galerkin finite…
We derive a grey linear diffusion equation for photons with respect to inertial (or lab-frame) space and time, using asymptotic analysis in 1D planar geometry. The solution of the equation is the comoving radiation energy density. Our…
We develop a semi-implicit algorithm for time-accurate simulation of the compressible Navier-Stokes equations, with special reference to wall-bounded flows. The method is based on linearization of the partial convective fluxes associated…
Quasi-linear hyperbolic systems with source terms introduce significant computational challenges due to the presence of a stiff source term. To address this, a finite volume Nessyahu-Tadmor (NT) central numerical scheme is explored and…
A transfer matrix method is presented for solving the scattering problem for the quasi one-dimensional massless Dirac equation applied to graphene in the presence of an arbitrary inhomogeneous electric and perpendicular magnetic field. It…
We propose an inexact low-rank source iteration with diffusion synthetic acceleration (SI-DSA) for solving the multidimensional steady-state radiative transfer equation (RTE) in the second-order formulation. The angular flux is represented…
We present compact semi-implicit finite difference schemes on structured grids for numerical solutions of the advection by an external velocity and by a speed in normal direction that are applicable in level set methods. The most involved…
Since the early 1970s, inversion techniques have become the most useful tool for inferring the magnetic, dynamic, and thermodynamic properties of the solar atmosphere. The intrinsic model dependence makes it necessary to formulate specific…
Dispersion of low-density rigid particles with complex geometries is ubiquitous in both natural and industrial environments. We show that while explicit methods for coupling the incompressible Navier-Stokes equations and Newton's equations…
In this paper, we study the remote estimation problem of a Markov process over a channel with a cost. We formulate this problem as an infinite horizon optimization problem with two players, i.e., a sensor and a monitor, that have distinct…
Low-rank methods for kinetic equations have attracted increasing attention due to their effectiveness in reducing the high dimensionality of phase space. In our previous work [G. Wang & J. Hu, J. Comput. Phys. 558 (2026) 114884], we…
This work aims to extend the residual distribution (RD) framework to stiff relaxation problems. The RD is a class of schemes which is used to solve hyperbolic system of partial differential equations. Up to our knowledge, it was used only…
We present a preconditioner based on spectral projection that is combined with a deflated Krylov subspace method for solving ill conditioned linear systems of equations. Our results show that the proposed algorithm requires many fewer…
In this paper, we revisit asymptotic stability for the two-dimensional incompressible porous media equation and the Stokes transport system in a periodic channel. It is well-known that a stratified density, which strictly decreases in the…
We develop a semi-discrete optimal transport scheme for the compressible semi-geostrophic equations, a system that plays an important role in modelling large-scale atmospheric dynamics and frontogenesis. Unlike the incompressible case, the…
We consider the inverse problem of reconstructing the optical parameters of the radiative transfer equation (RTE) from boundary measurements in the diffusion limit. In the diffusive regime (the Knudsen number $\mathsf{Kn}\ll 1$), the…
We develop an unconditionally energy-stable tensor-product space-time discretization framework for the solution of a linear kinetic transport equation in one space dimension. The kinetic equation is a simplified model of radiative transfer…
This article introduces an iterative method for solving nonsingular non-Hermitian positive semidefinite systems of linear equations. To construct the iteration process, the coefficient matrix is split into two non-Hermitian positive…
In the first part of the paper we consider accelerated first order optimization method for convex functions with $L$-Lipschitz-continuous gradient, that is able to automatically adapts to problems which satisfies Polyak-{\L}ojasiewicz…
The stable principal component pursuit (SPCP) problem is a non-smooth convex optimization problem, the solution of which has been shown both in theory and in practice to enable one to recover the low rank and sparse components of a matrix…