数值分析
In this paper, we propose an adaptive multilevel preconditioned Helmholtz-Jacobi-Davidson (PHJD) method for the Maxwell eigenvalue problem with singularities. The key idea in this work is to employ the local multilevel method for…
In this paper we introduce a new mathematical model describing the erosion process caused in carbonate stones by the dissolution of the porous matrix due to the penetration of carbonic acid present in the environment. Such model is…
In this article, we solve the Stokes and Navier-Stokes equations with the deal$.$II finite-element library. In particular, we use its multigrid, adaptive-mesh, and matrix-free infrastructures to design efficient linear and nonlinear…
Tumor growth beyond a critical size relies on the development of a functional vascular network, which ensures adequate oxygen and nutrient supply. In this work, we present a modeling framework based on an optimization-based 3D-1D coupling…
The paper deals with the 3D Stokes problem with Navier-Tresca stick-slip boundary conditions. A weak formulation of this problem leads to a variational inequality of the second kind, coupled with an equality constraint. This problem is then…
We present a comparative study of numerical methods for computingelectronic densities of states (DOS) in periodic systems. We provide a detailed analysis of the domain of validity of the Brillouincomplex deformation (BCD), a…
Large-scale linear systems of the form $Ax=b$ are often doubly-noisy, in the sense that both its measurement matrix $A$ and measurement vector $b$ are noisy. In this paper, we extend the relaxed greedy randomized Kaczmarz (RGRK) method to…
This paper introduces a unified model for thermo-poroelasticity and multiple-network poroelasticity, reformulated into a total-pressure-based system. We first establish the well-posedness of the problem via a Galerkin-based argument and…
Singularly perturbed partial differential equations arise in many applications, including magnetohydrodynamic duct flows, chemical reaction transport systems, and Poisson Boltzmann electrostatics. These problems are characterized by sharp…
We study the solution of large symmetric positive-definite linear systems in a matrix-free setting with a limited iteration budget. We focus on the preconditioned conjugate gradient (PCG) method with spectral preconditioning. Spectral…
The paper that follows describes a numerical algorithm to solve the parabolic-parabolic Keller--Segel system characterized by singular sensitivity and signal absorption in such a manner that the numerical approximations converge towards a…
In mathematical physics the Galois top, introduced by S. Adlaj, possesses a fixed point on one of two Galois axes through its center of mass. This heavy top has two algebraic motion invariants and an additional transcendental…
In this paper we investigate the influence of the discretization of PDE constraints on shape and topological derivatives. To this end, we study a tracking-type functional and a two-material Poisson problem in one spatial dimension. We…
We consider the numerical approximation of a nonlinear system of partial differential equations modeling magnetostriction in the small-strain regime consisting of the Landau--Lifshitz--Gilbert equation for the magnetization and the…
Lie groups and their actions are ubiquitous in the description of physical systems, and we explore implications in the setting of model order reduction (MOR). We present a novel framework of MOR via Lie groups, called MORLie, in which…
We propose a numerical method for the valuation of European-style options under two-asset infinite-activity exponential L\'evy models. Our method extends the effective approach developed by Wang, Wan & Forsyth (2007) for the 1-dimensional…
We present a consistent high-order staggered Lagrangian hydrodynamics framework designed to reconcile an underlying disparity in existing curvilinear formulations: the mismatch between quadrature-based "strong" mass conservation and the…
This work describes a new version of the Fast Multipole Method for summing pairwise particle interactions that arise from discretizing integral transforms and convolutions on the sphere. The kernel approximations use barycentric Lagrange…
We perform a backward stability analysis of preconditioned sketched GMRES [Nakatsukasa and Tropp, SIAM J. Matrix Anal. Appl, 2024] for solving linear systems $Ax=b$, and show that the backward stability at iteration $i$ depends on the…
This paper presents stability and accuracy analysis of a high-order explicit time stepping scheme introduced by \cite[Section 2.2]{Buvoli2019}, which exhibits superior stability compared to classical Adams-Bashforth. A conjecture that is…