2025-04-18 | | Total: 10
Physics-Informed Neural Networks (PINNs) have emerged as a powerful framework for solving partial differential equations (PDEs). However, their performance heavily relies on the strategy used to select training points. Conventional adaptive sampling methods, such as residual-based refinement, often require multi-round sampling and repeated retraining of PINNs, leading to computational inefficiency due to redundant points and costly gradient computations-particularly in high-dimensional or high-order derivative scenarios. To address these limitations, we propose RL-PINNs, a reinforcement learning(RL)-driven adaptive sampling framework that enables efficient training with only a single round of sampling. Our approach formulates adaptive sampling as a Markov decision process, where an RL agent dynamically selects optimal training points by maximizing a long-term utility metric. Critically, we replace gradient-dependent residual metrics with a computationally efficient function variation as the reward signal, eliminating the overhead of derivative calculations. Furthermore, we employ a delayed reward mechanism to prioritize long-term training stability over short-term gains. Extensive experiments across diverse PDE benchmarks, including low-regular, nonlinear, high-dimensional, and high-order problems, demonstrate that RL-PINNs significantly outperforms existing residual-driven adaptive methods in accuracy. Notably, RL-PINNs achieve this with negligible sampling overhead, making them scalable to high-dimensional and high-order problems.
We propose a characterization of a $p$-Laplace higher eigenvalue based on the inverse iteration method with balancing the Rayleigh quotients of the positive and negative parts of solutions to consecutive $p$-Poisson equations. The approach relies on the second eigenvalue's minimax properties, but the actual limiting eigenvalue depends on the choice of initial function. The well-posedness and convergence of the iterative scheme are proved. Moreover, we provide the corresponding numerical computations. As auxiliary results, which also have an independent interest, we provide several properties of certain $p$-Poisson problems.
Magnetic Particle Imaging (MPI) is a tomographic imaging modality capable of real-time, high-sensitivity mapping of superparamagnetic iron oxide nanoparticles. Model-based image reconstruction provides an alternative to conventional methods that rely on a measured system matrix, eliminating the need for laborious calibration measurements. Nevertheless, model-based approaches must account for the complexities of the imaging chain to maintain high image quality. A recently proposed direct reconstruction method leverages weighted Chebyshev polynomials in the frequency domain, removing the need for a simulated system matrix. However, the underlying model neglects key physical effects, such as nanoparticle anisotropy, leading to distortions in reconstructed images. To mitigate these artifacts, an adapted direct Chebyshev reconstruction (DCR) method incorporates a spatially variant deconvolution step, significantly improving reconstruction accuracy at the cost of increased computational demands. In this work, we evaluate the adapted DCR on six experimental phantoms, demonstrating enhanced reconstruction quality in real measurements and achieving image fidelity comparable to or exceeding that of simulated system matrix reconstruction. Furthermore, we introduce an efficient approximation for the spatially variable deconvolution, reducing both runtime and memory consumption while maintaining accuracy. This method achieves computational complexity of O(N log N ), making it particularly beneficial for high-resolution and three-dimensional imaging. Our results highlight the potential of the adapted DCR approach for improving model-based MPI reconstruction in practical applications.
The existence of explicit symplectic integrators for general nonseparable Hamiltonian systems is an open and important problem in both numerical analysis and computing in science and engineering, as explicit integrators are usually more efficient than the implicit integrators of the same order of accuracy. Up to now, all responses to this problem are negative. That is, there exist explicit symplectic integrators only for some special nonseparable Hamiltonian systems, whereas the universal design involving explicit symplectic integrators for general nonseparable Hamiltonian systems has not yet been studied sufficiently. In this paper, we present a constructive proof for the existence of explicit symplectic integrators for general nonseparable Hamiltonian systems via finding explicit symplectic mappings under which the special submanifold of the extended phase space is invariant. It turns out that the proposed explicit integrators are symplectic in both the extended phase space and the original phase space. Moreover, on the basis of the global modified Hamiltonians of the proposed integrators, the backward error analysis is made via a parameter relaxation and restriction technique to show the linear growth of global errors and the near-preservation of first integrals. In particular, the effective estimated time interval is nearly the same as classical implicit symplectic integrators when applied to (near-) integrable Hamiltonian systems. Numerical experiments with a completely integrable nonseparable Hamiltonian and a nonintegrable nonseparable Hamiltonian illustrate the good long-term behavior and high efficiency of the explicit symplectic integrators proposed and analyzed in this paper.
Stochastic differential equations (SDEs) on Riemannian manifolds have numerous applications in system identification and control. However, geometry-preserving numerical methods for simulating Riemannian SDEs remain relatively underdeveloped. In this paper, we propose the Exponential Euler-Maruyama (Exp-EM) scheme for approximating solutions of SDEs on Riemannian manifolds. The Exp-EM scheme is both geometry-preserving and computationally tractable. We establish a strong convergence rate of $\mathcal{O}(\delta^{\frac{1 - \epsilon}{2}})$ for the Exp-EM scheme, which extends previous results obtained for specific manifolds to a more general setting. Numerical simulations are provided to illustrate our theoretical findings.
In this paper, we study the numerical simulation of stochastic differential equations (SDEs) on the special orthogonal Lie group $\text{SO}(n)$. We propose a geometry-preserving numerical scheme based on the stochastic tangent space parametrization (S-TaSP) method for state-dependent multiplicative SDEs on $\text{SO}(n)$. The convergence analysis of the S-TaSP scheme establishes a strong convergence order of $\mathcal{O}(\delta^{\frac{1-\epsilon}{2}})$, which matches the convergence order of the previous stochastic Lie Euler-Maruyama scheme while avoiding the computational cost of the exponential map. Numerical simulation illustrates the theoretical results.
We construct an efficient primal-dual forward-backward (PDFB) splitting method for computing a class of minimizing movement schemes with nonlinear mobility transport distances, and apply it to computing Wasserstein-like gradient flows. This approach introduces a novel saddle point formulation for the minimizing movement schemes, leveraging a support function form from the Benamou-Brenier dynamical formulation of optimal transport. The resulting framework allows for flexible computation of Wasserstein-like gradient flows by solving the corresponding saddle point problem at the fully discrete level, and can be easily extended to handle general nonlinear mobilities. We also provide a detailed convergence analysis of the PDFB splitting method, along with practical remarks on its implementation and application. The effectiveness of the method is demonstrated through several challenging numerical examples.
Approximating a manifold-valued function from samples of input-output pairs consists of modeling the relationship between an input from a vector space and an output on a Riemannian manifold. We propose a function approximation method that leverages and unifies two prior techniques: (i) approximating a pullback to the tangent space, and (ii) the Riemannian moving least squares method. The core idea of the new scheme is to combine pullbacks to multiple tangent spaces with a weighted Fréchet mean. The effectiveness of this approach is illustrated with numerical experiments on model problems from parametric model order reduction.
This paper is concerned with non-uniform fully-mixed FEMs for dynamic coupled Stokes-Darcy model with the well-known Beavers-Joseph-Saffman (BJS) interface condition. In particular, a decoupled algorithm with the lowest-order mixed non-uniform FE approximations (MINI for the Stokes equation and RT0-DG0 for the Darcy equation) and the classical Nitsche-type penalty is studied. The method with the combined approximation of different orders is commonly used in practical simulations. However, the optimal error analysis of methods with non-uniform approximations for the coupled Stokes-Darcy flow model has remained challenging, although the analysis for uniform approximations has been well done. The key question is how the lower-order approximation to the Darcy flow influences the accuracy of the Stokes solution through the interface condition. In this paper, we prove that the decoupled algorithm provides a truly optimal convergence rate in L^2-norm in spatial direction: O(h^2) for Stokes velocity and O(h) for Darcy flow in the coupled Stokes-Darcy model. This implies that the lower-order approximation to the Darcy flow does not pollute the accuracy of numerical velocity for Stokes flow. The analysis presented in this paper is based on a well-designed Stokes-Darcy Ritz projection and given for a dynamic coupled model. The optimal error estimate holds for more general combined approximations and more general coupled models, including the corresponding model of steady-state Stokes-Darcy flows and the model of coupled dynamic Stokes and steady-state Darcy flows. Numerical results confirm our theoretical analysis and show that the decoupled algorithm is efficient.
This paper presents a generalized energy-based modeling framework extending recent formulations tailored for differential-algebraic equations. The proposed structure, inspired by the port-Hamiltonian formalism, ensures passivity, preserves the power balance, and facilitates the consistent interconnection of subsystems. A particular focus is put on low-frequency power applications in electrical engineering. Stranded, solid, and foil conductor models are investigated in the context of the eddy current problem. Each conductor model is shown to fit into the generalized energy-based structure, which allows their structure-preserving coupling with electrical circuits described by modified nodal analysis. Theoretical developments are validated through a numerical simulation of an oscillator circuit, demonstrating energy conservation in lossless scenarios and controlled dissipation when eddy currents are present.