数值分析

Modern high performance computers are massively parallel; for many PDE applications spatial parallelism saturates long before the computer's capability is reached. Parallel-in-time methods enable further speedup beyond spatial saturation by…

FIAT (the FInite element Automatic Tabulator) provides a powerful Python library for the generation and evaluation of finite element basis functions on a reference element. This release paper describes recent improvements to FIAT aimed at…

This paper proposes a simple technique of curve and surface construction with B-splines. Given a control polygon or a control mesh together with node ordinates corresponding to all control points, a rational curve or surface is obtained by…

Modelling the evolution process for the growth of microalgae in an artificial pond is a huge challenge, given the complex interaction between hydrodynamics and biological processes occurring across various timescales. In this paper, we…

For a positive integer $d$, the $d$-dimensional Chebyshev-Frolov lattice is the $\mathbb{Z}$-lattice in $\mathbb{R}^d$ generated by the Vandermonde matrix associated to the roots of the $d$-dimensional Chebyshev polynomial. It is important…

We investigate multivariate integration for a space of infinitely times differentiable functions $\mathcal{F}_{s, \boldsymbol{u}} := \{f \in C^\infty [0,1]^s \mid \| f \|_{\mathcal{F}_{s, \boldsymbol{u}}} < \infty \}$, where $\| f…

We establish formulas for the $b$-adic Walsh coefficients of functions in $C^\alpha[0,1]$ for an integer $\alpha \geq 1$ and give upper bounds on the Walsh coefficients of these functions. We also study the Walsh coefficients of periodic…

In this paper, we study quasi-Monte Carlo (QMC) rules for numerical integration. J. Dick proved a Koksma-Hlawka type inequality for $\alpha$-smooth integrands and gave an explicit construction of QMC rules achieving the optimal rate of…

Reachable Minimally supported (RM) B-splines have been recently introduced as a novel B-spline--like basis. They feature local linear independence and admit a fast de Boor--like evaluation algorithm. These properties make them particularly…

Physics-informed neural networks (PINNs) often struggle with multi-scale PDEs featuring sharp gradients and nontrivial boundary conditions, as the physics residual and boundary enforcement compete during optimization. We present a…

Partial differential equations (PDEs) describing thermodynamically isolated systems typically possess conserved quantities (like mass, momentum, and energy) and dissipated quantities (like entropy). Preserving these conservation and…

Originally introduced to describe a transition region in stars, the magnetic rotating shallow water (MRSW) model is now used in many solar physics and geophysical applications. Derived from the 3-D incompressible magnetohydrodynamic system,…

When extended to two-phase flows, weakly compressible models lead to a non-conservative system, which precludes its treatment using standard finite volume techniques. In this paper, a novel HLLC-type path-conservative scheme is formulated…

Graph autoencoders have gained attention in nonlinear reduced-order modeling of parameterized partial differential equations defined on unstructured grids. Despite they provide a geometrically consistent way of treating complex domains,…

In this paper, we extend the previous work on absolutely convergent fixed-point fast sweeping WENO methods by Li et al. (J. Comput. Phys. 443: 110516, 2021) and design a fifth-order hybrid fast sweeping scheme for solving steady state…

In this paper, we propose a high-order extension of the multiscale method introduced by the authors in [SIAM J. Numer. Anal., 63(4) (2025), pp. 1617--1641] for heterogeneous Stokes problems, while also providing several other improvements,…

The paper considers function-valued tensors, viewed as multidimensional arrays with entries in an abstract Hilbert space. Despite the absence of the algebraic structure of a field, the geometric inner-product structure suffices to introduce…

Numerical investigations of partial differential equations with hysteresis have largely focused on simulations, leaving numerical error analysis unexplored and relying mainly on derivative-free nonlinear solvers. This work establishes…

In this work we construct a high-order Asymptotic-Preserving (AP) Implicit-Explicit (IMEX) scheme for the ES-BGK model for gas mixtures introduced in [Brull, Commun. Math. Sci., 2015]. The time discretization is based on the IMEX strategy…

This article introduces a novel approach for broken-FEEC (Finite Element Exterior Calculus), extending its application to locally refined spline spaces with non-matching interfaces. Traditional broken-FEEC allows for discontinuous…