2025-03-28 | | Total: 15
This paper presents a novel inverse Lax-Wendroff (ILW) boundary treatment for finite difference Hermite weighted essentially non-oscillatory (HWENO) schemes to solve hyperbolic conservation laws on arbitrary geometries. The complex geometric domain is divided by a uniform Cartesian grid, resulting in challenge in boundary treatment. The proposed ILW boundary treatment could provide high order approximations of both solution values and spatial derivatives at ghost points outside the computational domain. Distinct from existing ILW approaches, our boundary treatment constructs the extrapolation via optimized through a least squares formulation, coupled with the spatial derivatives at the boundary obtained via the ILW procedure. Theoretical analysis indicates that compared with other ILW methods, our proposed one would require fewer terms by using the relatively complicated ILW procedure and thus improve computational efficiency while preserving accuracy and stability. The effectiveness and robustness of the method are validated through numerical experiments.
We present a rational filter for computing all eigenvalues of a symmetric definite eigenvalue problem lying in an interval on the real axis. The linear systems arising from the filter embedded in the subspace iteration framework, are solved via a preconditioned Krylov method. The choice of the poles of the filter is based on two criteria. On the one hand, the filter should enhance the eigenvalues in the interval of interest, which suggests that the poles should be chosen close to or in the interval. On the other hand, the choice of poles has an important impact on the convergence speed of the iterative method. For the solution of problems arising from vibrations, the two criteria contradict each other, since fast convergence of the eigensolver requires poles to be in or close to the interval, whereas the iterative linear system solver becomes cheaper when the poles lie further away from the eigenvalues. In the paper, we propose a selection of poles inspired by the shifted Laplace preconditioner for the Helmholtz equation. We show numerical experiments from finite element models of vibrations. We compare the shifted Laplace rational filter with rational filters based on quadrature rules for contour integration.
Volterra's integral equations with local and nonlocal loads represent the novel class of integral equations that have attracted considerable attention in recent years. These equations are a generalisation of the classic Volterra integral equations, which were first introduced by Vito Volterra in the late 19th century. The loaded Volterra integral equations are characterised by the presence of a load which complicates the process of their theoretical and numerical study. Sometimes these equation are called the equations with ``frozen'' argument. The present work is devoted to the study of Volterra equations with locally loaded integral operators. The existence and uniquness theorems are proved. Among the main contributions is the collocation method for approximate solution of such equations based on the piecewise linear approximation. To confirm the convergence of the method, a number of numerical results for solving model problems are provided.
This paper considers the problem of detecting adjoint mismatch for two linear maps. To clarify, this means that we aim to calculate the operator norm for the difference of two linear maps, where for one we only have a black-box implementation for the evaluation of the map, and for the other we only have a black-box for the evaluation of the adjoint map. We give two stochastic algorithms for which we prove the almost sure convergence to the operator norm. The algorithm is a random search method for a generalization of the Rayleigh quotient and uses optimal step sizes. Additionally, a convergence analysis is done for the corresponding singular vector and the respective eigenvalue equation.
The Koopman-Hill projection method is used to approximate the fundamental solution matrix of linear time-periodic ordinary differential equations, possibly stemming from linearization around a periodic solution of a nonlinear dynamical system. By expressing both the true fundamental solution and its approximation as series, we derive an upper bound for the approximation error that decays exponentially with the size of the Hill matrix. Exponential decay of the Fourier coefficients of the system dynamics is key to guarantee convergence. The paper also analyzes a subharmonic formulation that improves the convergence rate. Two numerical examples, including a Duffing oscillator, illustrate the theoretical findings.
This paper focuses on continuous data assimilation (CDA) for the Navier-Stokes equations with nonlinear slip boundary conditions. CDA methods are typically employed to recover the original system when initial data or viscosity coefficients are unknown, by incorporating a feedback control term generated by observational data over a time period. In this study, based on a regularized form derived from the variational inequalities of the Navier-Stokes equations with nonlinear slip boundary conditions, we first investigate the classical CDA problem when initial data is absent. After establishing the existence, uniqueness and regularity of the solution, we prove its exponential convergence with respect to the time. Additionally, we extend the CDA to address the problem of missing viscosity coefficients and analyze its convergence order, too. Furthermore, utilizing the predictive capabilities of partial evolutionary tensor neural networks (pETNNs) for time-dependent problems, we propose a novel CDA by replacing observational data with predictions got by pETNNs. Compared with the classical CDA, the new one can achieve similar approximation accuracy but need much less computational cost. Some numerical experiments are presented, which not only validate the theoretical results, but also demonstrate the efficiency of the CDA.
The inverse reflector problem aims to design a freeform reflecting surface that can direct the light from a specified source to produce the desired illumination in the target area, which is significant in the field of geometrical non-imaging optics. Mathematically, it can be formulated as an optimization problem, which is exactly the optimal transportation problem (OT) when the target is in the far field. The gradient of OT is governed by the generalized Monge-Amp`ere equation that models the far-field reflector system. Based on the gradient, this work presents a Sobolev gradient descent method implemented within a finite element framework to solve the corresponding OT. Convergence of the method is established and numerical examples are provided to demonstrate the effectiveness of the method.
A Taylor method for solving an ordinary differential equation initial-value problem $\dot x = f(t,x)$, $x(t_0) = x_0$, computes the Taylor series (TS) of the solution at the current point, truncated to some order, and then advances to the next point by summing the TS with a suitable step size. A standard ODE method (e.g. Runge-Kutta) treats function $f$ as a black box, but a Taylor solver requires $f$ to be preprocessed into a code-list of elementary operations that it interprets as operations on (truncated) TS. The trade-off for this extra work includes arbitrary order, typically enabling much larger step sizes. For a standard function, such as $\exp$, this means evaluating $v(t)=\exp(u(t))$, where $u(t),v(t)$ are TS. The sub-ODE method applies the ODE $d v/d u=v$, obeyed by $v=\exp(u)$, to in-line this operation as $\dot v=v\dot u$. This gives economy of implementation: each function that satisfies a simple ODE goes into the "Taylor library" with a few lines of code--not needing a separate recurrence relation, which is the typical approach. Mathematically, however, the use of sub-ODEs generally transforms the original ODE into a differential-algebraic system, making it nontrivial to ensure a sound system of recurrences for Taylor coefficients. We prove that, regardless of how many sub-ODEs are incorporated into $f$, this approach guarantees a sound system. We introduce our sub-ODE-based Matlab ODE solver and show that its performance compares favorably with solvers from the Matlab ODE suite.
ConicCurv is a new derivative-free algorithm to estimate the curvature of a plane curve from a sample of data points. It is based on a known tangent estimator method grounded on classic results of Projective Geometry and Bézier rational conic curves. The curvature values estimated by ConicCurv are invariant to Euclidean changes of coordinates and reproduce the exact curvature values if the data are sampled from a conic. We show that ConicCurv< has convergence order $3$ and, if the sample points are uniformly arc-length distributed, the convergence order is $4$. The performance of ConicCurv is compared with some of the most frequently used algorithms to estimate curvatures and its performance is illustrated in the calculation of the elastic energy of subdivision curves and the location of L-curves corners.
We extend the Geometric Refinement Transform (GRT) by introducing centroidal Voronoi tessellations (CVTs) into the refinement process, enhancing symmetry, reconstruction accuracy, and numerical stability. By applying Lloyds algorithm at each refinement level, we minimize centroidal energy and generate Voronoi regions that better align with the functions underlying structure. This approach reduces geometric distortion, suppresses reconstruction error, and provides a natural framework for adaptive refinement. We analyze convergence properties, quantify the reduction in reconstruction error using Taylor-based estimates and Lipschitz continuous functions, and propose perturbation strategies to escape symmetry-preserving local minima. The resulting transform offers improved accuracy for applications in medical imaging, signal processing, and physics simulations, while preserving the theoretical completeness and stability guarantees of the original GRT framework.
We present and discuss a generalization of the popular MINI mixed finite element for the 2D Stokes equation by means of conforming virtual elements on polygonal meshes. We prove optimal error estimates for both velocity and pressure. Theoretical results are confirmed by several numerical tests performed with different choices of polynomial accuracy and meshes.
The non-Hermitian Bethe-Salpeter eigenvalue problem is a structured eigenproblem, with real eigenvalues coming in pairs $\{\lambda,-\lambda\}$ where the corresponding pair of eigenvectors are closely related, and furthermore the left eigenvectors can be trivially obtained from the right ones. We exploit these properties to devise three variants of structure-preserving Lanczos eigensolvers to compute a subset of eigenvalues (those of either smallest or largest magnitude) together with their corresponding right and left eigenvectors. For this to be effective in real applications, we need to incorporate a thick-restart technique in a way that the overall computation preserves the problem structure. The new methods are validated in an implementation within the SLEPc library using several test matrices, some of them coming from the Yambo materials science code.
In this paper, the stochastic theta (ST) method is investigated for a class of stochastic differential equations driven by a time-changed Brownian motion, whose coefficients are time-space-dependent and satisfy the local Lipschitz condition. It is proved that under the local Lipschitz and some additional assumptions, the ST method with $\theta\in[1/2,1]$ is strongly convergent. It is also obtained that, for all positive stepsizes, the ST method with $\theta\in[1/2,1]$ is asymptotically mean square stable under a coercivity condition. With some restrictions on the stepsize, the ST method with $\theta\in[0,1/2)$ is asymptotically mean square stable under a stronger assumption. Some numerical simulations are presented to illustrate the theoretical results.
This article builds on the recently proposed RB-ML-ROM approach for parameterized parabolic PDEs and proposes a novel hierarchical Trust Region algorithm for solving parabolic PDE constrained optimization problems. Instead of using a traditional offline/online splitting approach for model order reduction, we adopt an active learning or enrichment strategy to construct a multi-fidelity hierarchy of reduced order models on-the-fly during the outer optimization loop. The multi-fidelity surrogate model consists of a full order model, a reduced order model and a machine learning model. The proposed hierarchical framework adaptively updates its hierarchy when querying parameters, utilizing a rigorous a posteriori error estimator in an error aware trust region framework. Numerical experiments are given to demonstrate the efficiency of the proposed approach.
Message-Passing Monte Carlo (MPMC) was recently introduced as a novel low-discrepancy sampling approach leveraging tools from geometric deep learning. While originally designed for generating uniform point sets, we extend this framework to sample from general multivariate probability distributions with known probability density function. Our proposed method, Stein-Message-Passing Monte Carlo (Stein-MPMC), minimizes a kernelized Stein discrepancy, ensuring improved sample quality. Finally, we show that Stein-MPMC outperforms competing methods, such as Stein Variational Gradient Descent and (greedy) Stein Points, by achieving a lower Stein discrepancy.