数值分析
Time-harmonic wave propagation problems, especially those governed by Maxwell's equations, pose significant computational challenges due to the non-self-adjoint nature of the operators and the large, non-Hermitian linear systems resulting…
Diffusion kernels over graphs have been widely utilized as effective tools in various applications due to their ability to accurately model the flow of information through nodes and edges. However, there is a notable gap in the literature…
We are interested in the Euler-Maruyama dicretization of the SDE dXt =b(t,Xt)dt+ dZt, X0 =x$\in$Rd, where Zt is a symmetric isotropic d-dimensional $\alpha$-stable process, $\alpha$ $\in$ (1, 2] and the drift b $\in$ L$\infty$…
We examine nonlinear Kolmogorov partial differential equations (PDEs). Here the nonlinear part of the PDE comes from its Hamiltonian where one maximizes over all possible drift and diffusion coefficients which fall within a…
We construct a family of finite element sub-complexes of the conformal complex on tetrahedral meshes and show their exactness on contractible domains. This complex includes vector fields and symmetric and traceless tensor fields, connected…
We study time integration schemes for $\dot H^1$-solutions to the energy-(sub)critical semilinear wave equation on $\mathbb{R}^3$. We show first-order convergence in $L^2$ for the Lie splitting and convergence order $3/2$ for a corrected…
Time-dependent convection-dominated convection-diffusion problems are considered. We develop a moving mesh streamline upwind Petrov-Galerkin (MM-SUPG) method by combining residual-based SUPG stabilization with a metric-based moving mesh PDE…
We propose a new random sketching approach for embedding high-dimensional Hilbert-Schmidt operators, using random input-output pairs. Such operator can then be approximated in a low-dimensional subspace of operators by solving a small…
We consider the compressible Euler system with anelastic scaling, modeling isentropic flows under the influence of gravity. In the zero-Mach-number limit, the solution of the compressible Euler system converges to a variable density…
In this paper, we study a tensor-based method for the numerical solution of a class of diffusion--reaction equations defined on spatial domains that admit common curvilinear coordinate representations. Typical examples in 2D include disks…
Time-Resolved Spectroscopy (TRS) is a powerful modality for non-invasive characterization of turbid media. However, extracting optical properties, absorption $\mu_a$ and reduced scattering $\mu_s'$, from 3D stochastic measurements remains…
The Stokes problem with non-homogeneous Dirichlet boundary condition is solved numerically using conforming discretizations and an approximation of the boundary datum in the corresponding trace space. Optimal discretization error estimates…
In this contribution we present recent developments in the formulation and solution of Initial Boundary Value Problems (IBVPs). Building upon a modern variational action formulation of classical dynamics, we treat Initial Boundary Value…
We study the numerical approximation of a class of degenerate parabolic stochastic partial differential equations on non-compact metric graphs, which naturally arise in the asymptotic analysis of Hamiltonian flows under small noise…
The high-precision solution of the radial Dirac equation is fundamental to relativistic quantum chemistry, essential for reliable pseudopotential generation and all-electron electronic structure methods. However, standard basis-set…
This paper presents a convergence analysis for the Hessian Discretisation Method (HDM) applied to fourth-order semilinear elliptic equations involving a trilinear nonlinearity and general source, based on two complementary approaches. The…
Sparse grids are popular tools for high-dimensional function approximation. In this work, we study the use of sparse grids for interpolation using separable Mat\'ern kernels…
We consider an elastic model for a circular arch that incorporates membrane, transverse shear, and bending effects. The central line of the arch is partitioned into elements, and an ultra-weak variational formulation is developed alongside…
Particle breakage due to collisional interactions plays a vital role in the development of several phenomena in science and engineering. The nonlinear collisional breakage equations (NCBEs) are a significant set of equations in this…
We construct and analyse a $hp$-FE/BE coupling on non-matching meshes, based on Nitsche's method. Both the mesh size and the polynomial degree are changed to improve accuracy. Nitsche's method leads to a positive definite formulation, thus,…