数值分析
We present an anisotropic mesh adaptation procedure based on Riemannian metrics for the simulation of two-phase incompressible flows with non-matching densities. The system dynamics are governed by the Cahn-Hilliard Navier-Stokes (CHNS)…
This paper investigates the simultaneous identification of a spatially dependent potential and the initial condition in a subdiffusion model based on two terminal observations. The existence, uniqueness, and conditional stability of the…
Rooted trees are essential for describing numerical schemes via the so-called B-series. They have also been used extensively in rough analysis for expanding solutions of singular Stochastic Partial Differential Equations (SPDEs). When one…
We present a data-assisted iterative regularization method for solving ill-posed inverse problems. The proposed approach, termed \texttt{IRMGL+\(\Psi\)}, integrates classical iterative techniques with a data-driven regularization term…
This report showcases the role of, and future directions for, the field of Randomized Numerical Linear Algebra (RNLA) in a selection of scientific applications. These applications span the domains of imaging, genomics and dynamical systems,…
We analyze an advection-diffusion-reaction problem with non-homogeneous boundary conditions that models the chromatography process, a vital stage in bioseparation. We prove stability and error estimates for both constant and affine…
We discuss the order, efficiency, stability and positivity of several meshless schemes for linear scalar hyperbolic equations. Meshless schemes are Generalised Finite Difference Methods (GFDMs) for arbitrary irregular grids in which there…
The canonical polyadic (CP) decomposition is one of the most widely used tensor decomposition techniques. The conventional CP decomposition algorithm combines alternating least squares (ALS) with the normal equation. However, the normal…
We use linear programming bounds to analyze point sets in the torus with respect to their optimality for problems in discrepancy theory and quasi-Monte Carlo methods. These concepts will be unified by introducing tensor product energies. We…
Optimised lightweight structures, such as shallow domes and slender towers, are prone to sudden buckling failure because geometric uncertainties/imperfections can lead to a drastic reduction in their buckling loads. We introduce a framework…
We develop the novel method of artificial barriers for scalar stochastic differential equations (SDEs) and use it to construct boundary-preserving numerical schemes for strong approximation of scalar SDEs, possibly with non-globally…
Block classical Gram-Schmidt (BCGS) is commonly used for orthogonalizing a set of vectors $X$ in distributed computing environments due to its favorable communication properties relative to other orthogonalization approaches, such as…
In this paper we design, analyse and test domain decomposition methods for linear systems of equations arising from conforming finite element discretisations of positive Maxwell-type equations, namely for $\mathbf{H}(\mathbf{curl})$…
This paper addresses the analysis and numerical assessment of a computational method for solving the Cahn--Hilliard equation defined on a surface. The proposed approach combines the stabilized trace finite element method for spatial…
This work describes and analyzes the domain derivative for a time-dependent acoustic scattering problem. We study the nonlinear operator that maps a sound-soft scattering object to the solution of the time-dependent wave equation evaluated…
In this manuscript, we propose matrix- and tensor-oriented methods for the numerical solution of the multidimensional evolutionary space-fractional complex Ginzburg--Landau equation. After a suitable spatial semidiscretization, the…
Exponential Runge-Kutta methods are a well-established tool for the numerical integration of parabolic evolution equations. However, these schemes are typically developed under the assumption of homogeneous boundary conditions. In this…
This work deals with two groups of spectral analysis results for matrices arising in fully implicit Runge-Kutta methods used for linear time-dependent partial differential equations. These were applied for different formulations of the same…
We present a framework for solving a broad class of ill-posed inverse problems governed by partial differential equations (PDEs), where the target coefficients of the forward operator are recovered through an iterative regularization scheme…
This paper studies the numerical approximation of the ground state of rotating Bose--Einstein condensates, formulated as the minimization of the Gross--Pitaevskii energy functional under a mass conservation constraint. To solve this…