2025-07-04 | | Total: 15
The current work investigates the effectiveness of block triangular preconditioners in accelerating and stabilizing the numerical solution of inverse source problems governed by time-space fractional diffusion equations (TSFDEs). We focus on the recovery of an unknown spatial source function in a multi-dimensional TSFDE, incorporating Caputo time-fractional derivatives and the fractional Laplacian. The inherent ill-posedness is addressed via a quasi-boundary value regularization, followed by a finite difference discretization that leads to large, structured linear systems. We develop and analyze a block triangular preconditioning strategy that mimics the coefficient matrix, while simplifying its structure for computational efficiency. Numerical experiments using the GMRES solver demonstrate that the proposed preconditioner significantly improve convergence rates, robustness, and accuracy, making it well-suited for large-scale, real-world inverse problems involving fractional modeling.
Efficient simulation of the semiclassical Schrödinger equation has garnered significant attention in the numerical analysis community. While controlling the error in the unitary evolution or the wavefunction typically requires the time step size to shrink as the semiclassical parameter h decreases, it has been observed -- and proved for first- and second-order Trotterization schemes -- that the error in certain classes of observables admits a time step size independent of h. In this work, we explicitly characterize this class of observables and present a new, simple algebraic proof of uniform-in-h error bounds for arbitrarily high-order Trotterization schemes. Our proof relies solely on the algebraic structure of the underlying operators in both the continuous and discrete settings. Unlike previous analyses, it avoids Egorov-type theorems and bypasses heavy semiclassical machinery. To our knowledge, this is the first proof of uniform-in-h observable error bounds for Trotterization in the semiclassical regime that relies only on algebraic structure, without invoking the semiclassical limit.
We propose a new nonconforming finite element method for solving Stokes interface problems. The method is constructed on local anisotropic mixed meshes, which are generated by fitting the interface through simple connection of intersection points on an interface-unfitted background mesh, as introduced in \cite{Hu2021optimal}. For triangular elements, we employ the standard CR element; for quadrilateral elements, a new rotated Q1-type element is used. We prove that this rotated Q1 element remains unisolvent and stable even on degenerate quadrilateral elements. Based on these properties, we further show that the space pair of CR-rotated Q1 elements (for velocity) and piecewise P0 spaces (for pressure) satisfies the inf-sup condition without requiring any stabilization terms. As established in our previous work \cite{Wang2025nonconforming}, the consistency error achieves the optimal convergence order without the need for penalty terms to control it. Finally, several numerical examples are provided to verify our theoretical results.
We address the initial source identification problem for the heat equation, a notably ill-posed inverse problem characterized by exponential instability. Departing from classical Tikhonov regularization, we propose a novel approach based on moment analysis of the heat flow, transforming the problem into a more stable inverse moment formulation. By evolving the measured terminal time moments backward through their governing ODE system, we recover the moments of the initial distribution. We then reconstruct the source by solving a convex optimization problem that minimizes the total variation of a measure subject to these moment constraints. This formulation naturally promotes sparsity, yielding atomic solutions that are sums of Dirac measures. Compared to existing methods, our moment-based approach reduces exponential error growth to polynomial growth with respect to the terminal time. We provide explicit error estimates on the recovered initial distributions in terms of moment order, terminal time, and measurement errors. In addition, we develop efficient numerical discretization schemes and demonstrate significant stability improvements of our approach through comprehensive numerical experiments.
We introduce two multiscale numerical schemes for the time integration of weakly nonlinear Schrödinger equations, built upon the discretization of Picard iterates of the solution. These high-order schemes are designed to achieve high precision with respect to the small nonlinearity parameter under particular CFL condition. By exploiting the scattering properties of these schemes thanks to a low-frequency projected linear flow, we also establish its uniform accuracy over long time horizons. Numerical simulations are provided to illustrate the theoretical results, and these schemes are further applied to investigate dynamics in the framework of wave turbulence.
Many elliptic boundary value problems exhibit an interior regularity property, which can be exploited to construct local approximation spaces that converge exponentially within function spaces satisfying this property. These spaces can be used to define local ansatz spaces within the framework of generalised finite element methods, leading to a better relation between dimensionality and convergence order. In this paper, we present a new technique for the construction of such spaces for Lipschitz domains. Instead of the commonly used approach based on eigenvalue problems it relies on extensions of approximations performed on the boundary. Hence, it improves the influence of the spatial dimension on the exponential convergence and allows to construct the local spaces by solving the original kind of variational problems on easily structured domains.
This paper presents a novel framework for goal-oriented optimal static sensor placement and dynamic sensor steering in PDE-constrained inverse problems, utilizing a Bayesian approach accelerated by low-rank approximations. The framework is applied to airborne contaminant tracking, extending recent dynamic sensor steering methods to complex geometries for computational efficiency. A C-optimal design criterion is employed to strategically place sensors, minimizing uncertainty in predictions. Numerical experiments validate the approach's effectiveness for source identification and monitoring, highlighting its potential for real-time decision-making in crisis management scenarios.
We propose and study a Particle-In-Cell (PIC) method based on the Crank-Nicolson time discretization for the Vlasov-Poisson system with a strong and inhomogeneous external magnetic field with fixed direction, where we focus on the motion of particles in the plane orthogonal to the magnetic field (so-called poloidal directions). In this regime, the time step can be subject to stability constraints related to the smallness of Larmor radius and plasma frequency [21]. To avoid this limitation, our approach is based on numerical schemes [9, 10, 12], providing a consistent PIC discretization of the guiding-center system taking into account variations of the magnetic field. We carry out some theoretical proofs and perform several numerical experiments to validate the method and its underlying concepts.
In this work, we propose, analyze, and test a novel computational low-complexity, linear, second-order, and unconditional energy-stable semi-discrete time-stepping scheme for the Cahn-Hilliard-Navier-Stokes (CHNS) system by employing the time filter technique. Firstly, the first-order semi-implicit backward Euler (BE) method is utilized to discretize the CHNS model; Secondly, the time filter, as a post-processing strategy, is incorporated into the BE scheme, requiring only minimal modifications to the existing BE framework to improve its temporal accuracy from first- to second-order. The unconditional energy stability and second-order temporal error estimations are obtained, and several numerical experiments are conducted to verify the theoretical results.
In photoacoustic tomography (PAT), the computation of the initial pressure distribution within an object from its time-dependent boundary measurements over time is considered. This problem can be approached from two well-established points of view: deterministically using regularisation methods, or stochastically using the Bayesian framework. Both approaches frequently require the solution of a variational problem. In the paper we elaborate the connection between these approaches by establishing the equivalence between a smoothing Matérn class of covariance operators and Sobolev embedding operator Es:Hs↪L2. We further discuss the use of a Wavelet-based implementation of the adjoint operator E∗s which also allows for efficient evaluations for certain Matérn covariance operators, leading to efficient implementations both in terms of computational effort as well as memory requirements. The proposed methods are validated with reconstructions for the photoacoustic problem.
In this paper, we develop an efficient asymptotic-preserving (AP) Monte Carlo (MC) method for frequency-dependent radiative transfer equations (RTEs), which is based on the AP-MC method proposed for the gray RTEs in \cite{shi2023efficient}. We follow the characteristics-based approach by Zhang et al. \cite{zhang2023asymptotic} to get a reformulated model, which couples a low dimension convection-diffusion-type equation for macroscopic quantities with a high dimension transport equation for the radiative intensity. To recover the correct free streaming limit due to frequency-dependency, we propose a correction to the reformulated macroscopic equation. The macroscopic system is solved using a hybrid method: convective fluxes are handled by a particle-based MC method, while diffusive fluxes are treated implicitly with central difference. To address the nonlinear coupling between radiative intensity and the Planck function across multiple frequency groups, we adopt a Picard iteration with a predictor-corrector procedure, which decouples a global nonlinear system into a linear system restricted to spatial dimension (independent of frequency) with scalar algebraic nonlinear equations. Once the macroscopic update is done, the transport equation, with a known emission source provided by the macroscopic variables, is efficiently solved using an implicit MC method. This approach enables larger time steps independent of the speed of light and also the frequency across a wide range, significantly enhancing computational efficiency, especially for frequency-dependent RTEs. Formal AP analysis in the diffusive scaling is established. Numerical experiments are performed to demonstrate the high efficiency and AP property of the proposed method.
This paper focuses on the numerical approximation of the linearized shallow water equations using hybridizable discontinuous Galerkin (HDG) methods, leveraging the Hamiltonian structure of the evolution system. First, we propose an equivalent formulation of the equations by introducing an auxiliary variable. Then, we discretize the space variables using HDG methods, resulting in a semi-discrete scheme that preserves a discrete version of the Hamiltonian structure. The use of an alternative formulation with the auxiliary variable is crucial for developing the HDG scheme that preserves this Hamiltonian structure. The resulting system is subsequently discretized in time using symplectic integrators, ensuring the energy conservation of the fully discrete scheme. We present numerical experiments that demonstrate optimal convergence rates for all variables and showcase the conservation of total energy, as well as the evolution of other physical quantities.
Analyzing electromagnetic fields in complex, multi-material environments presents substantial computational challenges. To address these, we propose a hybrid numerical method that couples discrete exterior calculus (DEC) with surface integral equations (SIE) in the potential-based formulation of Maxwell's equations. The method employs the magnetic vector and electric scalar potentials (A-Φ) under the Lorenz gauge, offering natural compatibility with multi-physics couplings and inherent immunity to low-frequency breakdown. To effectively handle both bounded and unbounded regions, we divide the computational domain: the inhomogeneous interior is discretized using DEC, a coordinate-free framework that preserves topological invariants and enables structure-preserving discretization on unstructured meshes, while the homogeneous exterior is treated using SIEs, which inherently satisfy the radiation condition and eliminate the need for artificial domain truncation. A key contribution of this work is a scalar reformulation of the SIEs, which reduces the number of surface integral operators from fourteen to two by expressing the problem in terms of the Cartesian components of the vector potential and their normal derivatives. This simplification motivates a corresponding adaptation in the DEC domain: each vector potential component is represented as a discrete 0-form, in contrast to the conventional 1-form representation. This novel treatment improves compatibility at the interface and significantly enhances numerical performance. The proposed hybrid method thus offers a unified, efficient, and physically consistent framework for solving electromagnetic scattering and radiation problems in complex geometries and heterogeneous materials
Accurately learning solution operators for time-dependent partial differential equations (PDEs) from sparse and irregular data remains a challenging task. Recurrent DeepONet extensions inherit the discrete-time limitations of sequence-to-sequence (seq2seq) RNN architectures, while neural-ODE surrogates cannot incorporate new inputs after initialization. We introduce NCDE-DeepONet, a continuous-time operator network that embeds a Neural Controlled Differential Equation (NCDE) in the branch and augments the trunk with explicit space-time coordinates. The NCDE encodes an entire load history as the solution of a controlled ODE driven by a spline-interpolated input path, making the representation input-resolution-independent: it encodes different input signal discretizations of the observed samples. The trunk then probes this latent path at arbitrary spatial locations and times, rendering the overall map output-resolution independent: predictions can be queried on meshes and time steps unseen during training without retraining or interpolation. Benchmarks on transient Poisson, elastodynamic, and thermoelastic problems confirm the robustness and accuracy of the framework, achieving almost instant solution prediction. These findings suggest that controlled dynamics provide a principled and efficient foundation for high-fidelity operator learning in transient mechanics.
Motivated by the need for efficient estimation of conditional expectations, we consider a least-squares function approximation problem with heavily polluted data. Existing methods that are powerful in the small noise regime are suboptimal when large noise is present. We propose a hybrid approach that combines Christoffel sampling with certain types of optimal experimental design to address this issue. We show that the proposed algorithm enjoys appropriate optimality properties for both sample point generation and noise mollification, leading to improved computational efficiency and sample complexity compared to existing methods. We also extend the algorithm to convex-constrained settings with similar theoretical guarantees. When the target function is defined as the expectation of a random field, we extend our approach to leverage adaptive random subspaces and establish results on the approximation capacity of the adaptive procedure. Our theoretical findings are supported by numerical studies on both synthetic data and on a more challenging stochastic simulation problem in computational finance.