数值分析
Electroencephalography (EEG) source imaging aims to infer brain activity from electrical potentials measured on the scalp. This is a difficult problem because many different source patterns can explain the same measurements. The result…
This paper extends and analyzes the high-order kernel regularization framework of Beale & Tlupova (arXiv:2510.13639) to all four on-surface boundary integral operators of the Helmholtz Calderon calculus in three dimensions: the…
A mixed accuracy framework for Runge--Kutta methods presented in [Grant, JSC 2022] has been shown to speed up the computation in diagonally implicit Runge--Kutta (DIRK) methods by using less expensive low accuracy approaches for the…
The double-exponential Sinc-collocation method is known as a super-accurate method for solving initial value problems of ordinary differential equations, for which the error decreases almost exponentially as a function of the number of…
This work presents the discontinuous Galerkin discretization of the consistent splitting scheme proposed by Liu [J. Liu, J. Comp. Phys., 228(19), 2009]. The method enforces the divergence-free constraint implicitly, removing…
Yee's finite-difference method preserves two crucial properties of Maxwell's equations -- locality and symplecticity -- and thereby enjoys two computational advantages: scalability on high-performance architectures and long-time numerical…
In this contribution we introduce a mixed interpolation-regression operator for functions defined in some domains of the plane. We focus the attention on the ellipse, an annulus and a polygon. An upper bound for such an operator is…
We study a model describing the slow flow of a fluid through a deformable, porous, elastic solid undergoing small deformations. The stress-strain relationship of the solid incorporates nonlinear effects, formulated as a perturbation of the…
We propose, analyze, and numerically validate a correction adaptive two-grid finite element method (CAT-GFEM) for nonselfadjoint or indefinite elliptic problems. In contrast to the adaptive two-grid finite element method (ATGFEM) of Li and…
We study the temporal-spatial regularity properties of tamed Euler approximations for L\'evy-driven SDEs with superlinearly growing drift and diffusion coefficients. We first introduce a novel tamed Euler-type scheme and establish its…
This paper proposes a new randomized design of digital nets in which the generating matrices are chosen to be random Hankel matrices. Compared with previous randomized designs of digital nets, this approach simplifies the construction…
The long-standing hypothesis formulated by Goreinov, Tyrtyshnikov and Zamarashkin \cite{GTZ1997} has recently been solved affirmatively in the case of real two-column matrices by Sengupta and Pautov \cite{SP2026}. In this paper, we consider…
Randomized neural networks (RaNNs) are attractive for partial differential equations (PDEs) because they replace expensive end-to-end training with a linear least-squares solve over randomized hidden features. Their practical performance,…
Numerically solving magnetohydrodynamic (MHD) equations faces many challenges: avoiding divergence error, maintaining positivity, and satisfying entropy conditions. Among discontinuous Galerkin (DG) schemes, there has been a modal version…
In our previous work [29], we proposed a class of high-order asymptotic preserving (AP) finite difference weighted essentially non-oscillatory (WENO) schemes for solving the shallow water equations (SWEs) with bottom topography and Manning…
We formulate a well posed interface formulation for canonical one-dimensional evaporation two-phase model problems (the Stefan and Sucking problems) commonly used to validate production codes. We focus on the interface between the vapor and…
A residual-based a posteriori error estimator is proposed for the incompressible Oseen problem in the convection-dominated regime. The SUPG/PSPG/grad-div stabilized finite element method is used as discretization. The error estimator…
Deep generative models and neural operators have demonstrated significant potential for 3D aerodynamic inference. However, they often face inherent challenges in maintaining physical consistency and preserving high-frequency features,…
Gauss-Lobatto quadrature nodes and weights are optimal for closed summation-by-parts (SBP) formulations based on polynomial approximation spaces in the sense that for a prescribed function space they yield an SBP operator of minimal…
Conventional numerical solvers for the radiative transfer equation (RTE) exhibit severe sensitivity to medium parameters. To address this, we propose an operator learning framework that approximates the RTE solution map as a function of…