数值分析
Multigrid methods have been a popular approach for solving linear systems arising from the discretization of partial differential equations (PDEs) for several decades. They are particularly effective for accelerating convergence rates with…
Compact Runge-Kutta (cRK) Flux Reconstruction (FR) methods are a variant of RKFR methods for hyperbolic conservation laws with a compact stencil including only immediate neighboring finite elements. We extend cRKFR methods to handle…
The Kirchhoff plate model plays a vital role in modeling, computing and analyzing the mechanical behaviors of thin plate structures. This study propose a novel fourth-order multi-scale (FOMS) computational method for high-accuracy and…
Multigrid solvers are among the most efficient methods for solving the Poisson equation, which is ubiquitous in computational physics. For example, in the context of incompressible flows, it is typically the costliest operation. The present…
We describe a methodology to build vectorial kinetic schemes, targetting the numerical solution of linear symmetric-hyperbolic systems of conservation laws -a minimal application case for those schemes. Precisely, we fully detail the…
We study the $L_p$-discrepancy of random point sets in high dimensions, with emphasis on small values of $p$. Although the classical $L_p$-discrepancy suffers from the curse of dimensionality for all $p \in (1,\infty)$, the gap between…
3D-0D coupled flow models are widely used across many application fields but remain challenging to solve. Implicit coupling introduces non-local terms, whereas explicit coupling results in only conditionally stable schemes. Furthermore,…
We introduce a novel monolithic formulation that employs Lagrange multipliers (LMs) to couple a fluid flow governed by the time-dependent Stokes equations with a poroelastic structure described by the Biot equations. The formulation is…
For stochastic differential equations (SDEs) with Markovian switching, whose drift and diffusion coefficients are allowed to contain superlinear terms, the backward Euler-Maruyama (BEM) method is proposed to approximate the invariant…
In this paper we develop the formalism of rational complex Bezier curves. This framework is a simple extension of the CAD paradigm, since it describes arc of curves in terms of control polygons and weights, which are extended to complex…
We derive and analyze numerical methods for underdamped (kinetic) Langevin dynamics in a domain with elastic reflection at the boundary. First-order approximations are based on an Euler-type scheme incorporating collision-handling at the…
In this paper, we study numerical approximations for stochastic differential equations (SDEs) that use adaptive step sizes. In particular, we consider a general setting where decisions to reduce step sizes are allowed to depend on the…
Due to the COVID-19 pandemic, there is an increasing demand for portable CT machines worldwide in order to diagnose patients in a variety of settings. This has led to a need for CT image reconstruction algorithms that can produce high…
We consider the Boltzmann equation with random uncertainties arising from the initial data and collision kernel in the {\it whole space}, along with their stochastic Galerkin (SG) approximations. By employing Green's function method, we…
This paper extends previous work on finitedifference schemes over staggered grids for infinite-dimensional port-Hamiltonian systems. In the one-dimensional setting, it generalizes the discretization approach originally developed for the…
We present a fully non-intrusive parametric reduced-order modeling (PROM) framework for geometrically nonlinear structures subject to geometric variations. The method builds upon equation-driven Galerkin ROMs constructed from vibration…
The theory of spectral methods for partial differential equations leads to infinite-dimensional matrices which represent the derivative operator with respect to an underlying orthonormal basis. Favourable properties of such differentiation…
In this work, we study the dynamics of similarity classes of tetrahedra generated by the longest-edge bisection (LEB) algorithm. Building on the normalization strategy introduced by Perdomo and Plaza for triangles, we construct a canonical…
Classical convergence analysis for kernel interpolation typically assumes that the target function $f$ lies in the reproducing kernel Hilbert space $\mathcal{H}_k\!\left(\Omega\right)$ induced by a kernel on a domain…
We present a two-level overlapping Schwarz preconditioner for three-dimensional problems discretized with the Virtual Element Method. Our approach handles general polyhedral meshes and irregular subdomains, extending the applicability of…