数值分析
In this work, we extend the meshfree generalized multiscale exponential integration framework introduced in Nikiforov et al. (2025) to the simulation of three-dimensional advection--diffusion problems in heterogeneous and high-contrast…
This paper develops a family of fast, structure-preserving numerical algorithms for the nonlinear Maxwell-Ampere Nernst-Planck equations. For the first-order scheme, the Slotboom transformation rewrites the Nernst-Planck equation to enable…
In this work, we develop and analyze a higher-order finite element method for the multidimensional fragmentation equation. To the best of our knowledge, this is the first study to establish a rigorous, conforming finite element framework…
A key problem in approximation theory is the recovery of high-dimensional functions from samples. In many cases, the functions of interest exhibit anisotropic smoothness, and, in many practical settings, the nature of this anisotropy may be…
This paper provides approximation orders for a class of nonlinear interpolation procedures for univariate data sampled over $\sigma$ quasi-uniform grids. The considered interpolation is built using both essentially nonoscillatory (ENO) and…
We develop a hybrid conservative finite-volume / bounded-interval multiwavelet formulation for the deterministic one-dimensional Buckley--Leverett equation. Because Buckley--Leverett transport is a nonlinear hyperbolic conservation law with…
We consider adaptive finite element methods (AFEMs) with inexact algebraic solvers for second-order symmetric linear elliptic diffusion problems. Optimal complexity of AFEM, i.e., optimal convergence rates with respect to the overall…
We propose a novel non-iterative domain decomposition time integrator for acoustic wave equations using a discontinuous Galerkin discretization in space. It is based on a local Crank-Nicolson approximation combined with a suitable local…
We introduce new adaptive schemes for the one- and two-dimensional hyperbolic systems of conservation laws. Our schemes are based on an adaption strategy recently introduced in [{\sc S. Chu, A. Kurganov, and I. Menshov}, Appl. Numer. Math.,…
We propose and analyze a non-iterative domain decomposition integrator for the linear acoustic wave equation. The core idea is to combine an implicit Crank-Nicolson step on spatial subdomains with a local prediction step at the subdomain…
We consider a prototypical parabolic SPDE with finite-dimensional multiplicative noise, which, subject to a nonnegative initial datum, has a unique nonnegative solution. Inspired by well-established techniques in the deterministic case, we…
The hybrid-high order (HHO) scheme has many successful applications including linear elasticity as the first step towards computational solid mechanics. The striking advantage is the simplicity among other higher-order nonconforming schemes…
Wave propagation problems governed by the Helmholtz equation remain among the most challenging in scientific computing, due to their indefinite nature. Domain decomposition methods with spectral coarse spaces have emerged as some of the…
This paper derives a discrete dual problem for a prototypical hybrid high-order method for convex minimization problems. The discrete primal and dual problem satisfy a weak convex duality that leads to a priori error estimates with…
This paper proposes a novel technique for the approximation of strong solutions $u \in C(\overline{\Omega}) \cap W^{2,n}_\mathrm{loc}(\Omega)$ to uniformly elliptic linear PDE of second order in nondivergence form with continuous leading…
Definite integrals with parameters of holonomic functions satisfy holonomic systems of linear partial differential equations. When we restrict parameters to a one dimensional curve, the system becomes a linear ordinary differential equation…
This article presents a computational framework for determining the optimal slip velocity of a microswimmer with arbitrary three-dimensional geometry suspended in a viscous fluid. The objective is to minimize the hydrodynamic power…
We present parameter-robust preconditioners for linear systems that arise after applying static condensation to a hybridizable discontinuous Galerkin (HDG) discretization of the time-dependent Stokes problem. Building upon the theoretical…
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…
We develop a cut finite element method (CutFEM) for convection-diffusion problems posed on mixed-dimensional domains, i.e., unions of manifolds of different dimensions arranged in a hierarchical structure where lower-dimensional components…