数值分析
We design and investigate efficient multigrid solvers for multiphase Stokes problems discretised via mixed-degree local discontinuous Galerkin methods. Using the template of a standard multigrid V-cycle, we develop a smoother analogous to…
This work is devoted to the numerical approximation of high-dimensional advection-diffusion equations. It is well-known that classical methods, such as the finite volume method, suffer from the curse of dimensionality, and that their time…
As a combination of the logarithmic transformation with the truncated Euler-Maruyama (TEM) scheme, the positivity-preserving logarithmic truncated Euler-Maruyama (LTEM) scheme has been generally developed for scalar stochastic differential…
The Riemann problem is fundamental in the computational modeling of hyperbolic partial differential equations, enabling the development of stable and accurate upwind schemes. While exact solvers provide robust upwinding fluxes, their high…
Our focus is on simulating the dynamics of non-interacting particles including the effects of an external potential, which, under certain assumptions, can be formally described by the Dean-Kawasaki equation. The Dean-Kawasaki equation can…
We present a general framework to construct symmetric, well-conditioned, cross-element compatible nodal distributions that can be used for high-order and high-dimensional finite elements. Starting from the inherent symmetries of an element…
We present an extension of an adaptive, partially matrix-free, Hierarchically Semi-Separable (HSS) matrix construction algorithm by Gorman et al. [SIAM J. Sci. Comput. 41(5), 2019] which uses Gaussian sketching operators to a broader class…
We develop finite element spaces of symmetric tensor products of two-forms with polynomial coefficients. In three dimensions, these give higher order finite element spaces of matrix fields with normal-normal continuity, which have…
Purpose of writing this paper is to solve a transcendental function containing a product of a variable and its double exponential by a unique method of approximation. If the value of the said product is given, then its inverse function is…
Scientific calculations involving multiplication, division, exponents, inverse exponents of real numbers, geometric mean, reciprocal, Euler number, logarithm, and antilogarithm are generally carried out using battery operated electronic…
As the discretization error for the solution of a partial differential equation (PDE) decreases, the precision required to store the corresponding coefficients naturally increases. Storing the solution's finite element coefficients…
We introduce a novel quadrature strategy for Isogeometric Analysis (IgA) boundary element discretizations, specifically tailored to collocation methods. Thanks to the dimensionality reduction and the natural handling of unbounded domains,…
We present a MATLAB package for reconstructing sound-speed images from transmission ultrasound data. The package is based on two-point ray tracing and implements two complementary inversion strategies for image reconstruction. The first is…
Due to added numerical stabilization (diffusion), the stationary states of numerical methods for hyperbolic problems need not be consistent discretizations of those of the PDEs. A closely related phenomenon is the lack of consistency of…
This paper introduces a new numerical method for approximating the Lambert W function in the real domain. The method transforms the function into a simpler form that allows iterative refinement of an initial guess. Two iterative strategies…
The trigonometric interpolation has been recently applied to solve a second-order Fredholm integro-differentiable equation (FIDE). It achieves high accuracy with a moderate size of grid points and effectively addresses singularities of…
A dynamic linear thermo-poroelasticity model, containing inertial and relaxation terms with second-order time derivatives, is investigated in this paper. The mathematical and numerical analysis of this model is performed in the frequency…
This paper presents the symmetric wave interpolation method for stable global interpolation using readily available equidistant points. Its key achievement is the integration of the practical utility of such points with the numerical…
We generalize the closest point method (CPM) to solve surface partial differential equations with general boundary conditions. The proposed extrapolation method provides a unified framework for treating a broad class of inhomogeneous…
Multidimensional continuous-domain inverse problems are often solved by the minimization of a loss functional, formed as the sum of a data fidelity and a regularization. In this work, we present a new construction where the regularization…