数值分析
The inverse scattering problem for biharmonic waves, governing flexural vibrations of elastic plates, presents fundamental analytical challenges distinct from acoustic inverse problems due to the fourth-order differential operator and…
We extend our two-scale neural-network method for scalar singularly perturbed problems with one small parameter to dynamical systems with multiple small parameters. To accommodate multiple small parameters, we use a single effective scale…
The tear film (TF) plays a critical role in maintaining ocular surface health, and its disruption through tear breakup (TBU) is closely associated with dry eye disease. Evaporation-driven thinning is a primary mechanism underlying TBU, yet…
We consider Gaussian Random Fields on metric graphs defined implicitly as the stationary solution to a fractional SPDE driven by Gaussian white noise. Sampling from the finite element approximation requires the Cholesky factorization of the…
High-order discontinuous Galerkin (DG) methods equipped with subcell finite-volume (FV) limiters provide an efficient framework for the simulation of nonlinear hyperbolic balance laws featuring shocks and complex flow structures. However,…
We develop local discontinuous Galerkin (LDG) methods for conservation laws with heterogeneous stochastic fluxes, where the Stratonovich-driven transport terms may be linear or nonlinear. Such equations arise, for example, in simplified…
This work investigates finite element approximations for a general class of elliptic hemivariational inequalities arising in semipermeable media. The proposed model incorporates non-isotropic and heterogeneous diffusion coefficients,…
Nonlocal MEMS equations exhibit finite-time quenching, or touchdown, which is difficult to capture numerically. We study a stagewise rescaling algorithm for a two-dimensional nonlocal MEMS equation in an asymptotically constant-feedback…
Local linear instability refers to the linearized discrete operator exhibiting perturbation growth exceeding that of the corresponding continuous linearized problem. In the context of nonlinear entropy-stable discretizations, we argue that…
A numerical framework for constructing full-dimensional quasi-periodic solutions in nearly integrable systems was recently developed by Fu and Shi[2026]. Based on an alternating scheme, this approach effectively overcomes the secular drift…
The incompressible Hall-magnetohydrodynamics (Hall--MHD) system presents substantial analytical and computational challenges due to its stiff, highly nonlinear Hall term and the strict requirement that the magnetic field remains solenoidal.…
We study a nonlocal Poisson problem with discontinuous source term and analyze how the regularity of the integral kernel determines the discontinuity structure of the corresponding solution. Under general assumptions on compactly supported…
We propose a family of low-rank, completely positive and trace preserving schemes for the Lindblad equation, a common model for open quantum systems. Low-rank representation is employed at two levels: the density matrix is factorized into…
The computation of generalized inverses of quaternion matrices is a fundamental problem in quaternion linear algebra, with wide-ranging applications in signal processing, image restoration, and multidimensional data analysis. This paper…
Solving inverse problems in cardiac electrophysiology consists in the recovery of physiological parameters from surface electrocardiogram (ECG) measurements, a task which is often computationally unfeasible due to the severe ill-posedness…
Semi-implicit semi-Lagrangian (SISL) methods are commonly used for the shallow water equations (SWE) because they allow for larger time steps than those permitted by the Courant-Friedrichs-Lewy (CFL) stability condition in Eulerian schemes.…
To address the initial singularity inherent in solutions to fractional partial differential equations (fPDEs), we propose an accelerated Alikhanov discretization formulation implemented on nonuniform time grids. Based on the…
We propose a new class of short matrix recurrences for the solution of nonsymmetric linear equations of the type $\mathbf{A}_1\mathbf{X}\mathbf{B}_1+\ldots+\mathbf{A}_p\mathbf{X}\mathbf{B}_p=CD^T$. These iterative methods combine local…
We investigate inverse scattering problems for Dirac equations that arise as continuum models of waveguide arrays. We first establish the well-posedness of the forward models. For the associated inverse problems, we develop the inverse Born…
We present an exact sampling algorithm for Pfaffian point processes based on a skew-symmetric analogue of the Cholesky factorization. This algorithm enables efficient sampling of a wide range of statistics arising in random matrix theory…