数值分析
In this paper we present a new boundary integral equation formulation for the solution of the elastostatic traction boundary value problem in two and three dimensions. The approach relies on the introduction of new layer potentials, called…
We propose several linear, fully decoupled numerical schemes with first- and second-order temporal accuracy for a novel Q-tensor-based two-phase hydrodynamic model describing the coupling of active nematic liquid crystal solutions with…
In this paper, Cimmino's classical reflection algorithm for solving the $n\times n$ nonsingular linear system $A\bx=\bb$ is analysed through the lens of spectral theory. Reformulating the weighted iteration as…
We propose a Weak-form Physics-Informed Neural Operator (WINO), a data-free framework that combines the efficiency of neural operators with the geometric flexibility of the $\varphi$-finite element method ($\varphi$-FEM). $\varphi$-FEM is…
We present a theoretical and numerical analysis of Monte Carlo methods for the estimation of statistical moments of random variables $X:\Omega\rightarrow E$ taking values in a Banach space $E$. For practical computation, we consider…
Return panels, covariances, and large feature matrices evolve one observation or one entry at a time, yet downstream models require an up-to-date low-rank factorization $A_t \approx U_t \Sigma_t V_t^\top$ on every tick -- a regime where…
The nudged elastic band (NEB) and Dimer methods are standard tools for computing minimum-energy paths and index-one saddle points in atomistic transition problems. They are increasingly driven by surrogate or learned force models, whose…
A linear one-dimensional singularly perturbed convection-diffusion problem is solved numerically after its solution is decomposed as $u_0+w$, where $u_0$, the corresponding reduced solution, is treated as a function known exactly or…
For some kernel matrices, low-rank approximations can be quickly obtained via analytic techniques. One important class of analytic methods that has received attention in recent years is based on the use of proxy points. Accuracy analysis…
In this paper, we study Runge--Kutta methods for the computation of ruin probabilities in the classical risk model through the associated Volterra integro-differential equation. The proposed framework combines fourth-order one-step and…
Holland-Simpson thin-wire finite-difference time-domain (FDTD) simulations of obliquely oriented closed-loop antennas exhibit persistent low-frequency parasitic currents because the current-deposition operator fails to conserve charge. This…
This paper presents a novel multilevel projection-based stabilization method for advection-dominated convection--diffusion problems within the framework of Isogeometric Analysis. The proposed approach extracts and penalizes fine-scale…
Berger & Giuliani (2024) developed a provably stable weighted state redistribution (SRD) algorithm for cut-cell meshes. A key limitation of their method is that, although flux redistribution naturally vanishes when updates are small, SRD…
We consider an unregularized optimal control problem subject to the steady-state Navier-Stokes equations. We derive the existence of optimal solutions and prove first- and second-order optimality conditions. To approximate solutions to the…
The Cahn-Hilliard equation is a fundamental model for describing phase separation phenomena in binary mixtures. Traditional numerical methods, such as finite difference and finite element methods, often incur substantial computational cost,…
The quest for an algorithm that solves an $n\times n$ linear system in $O(n^2)$ time complexity, or $O(n^2 \text{poly}(1/\epsilon))$ when solving up to $\epsilon$ relative error, is a long-standing open problem in numerical linear algebra…
The predictive simulation of fluid dynamics in densely packed microfluidic devices, such as Deterministic Lateral Displacement (DLD) arrays, stagnates with standard iterative solvers. We show that this failure is not algorithmic but rooted…
Extracting approximate eigenpairs from a prescribed subspace is of fundamental importance in eigenvalue computation. While projecting the target eigenvector onto the subspace yields satisfactory accuracy, extracting an approximate eigenpair…
This paper introduces inexact versions of several block-splitting preconditioners for solving the three-by-three block linear systems arising from a special class of indefinite least squares problems. We first establish the convergence…
We introduce a new method, dubbed Geometric Structure-Preserving Interpolation ($\Gamma$-SPIN) to preserve physics-constraints inherent in the material parameter limits of the finite-strain Cosserat micropolar model. The method advocates to…