数值分析
Jin-Xin relaxation is a method for approximating non-linear hyperbolic conservation laws by a linear system of hyperbolic equations with an $\varepsilon$ dependent stiff source term. The system formally relaxes to the original conservation…
In this paper, we first consider linear 2D and 3D convection-diffusion-reaction equations $-\nabla\cdot (\kappa \nabla u) + {\bm v} \cdot \nabla u + \lambda u = \phi$ and $u_t - \nabla\cdot (\kappa \nabla u) + {\bm v} \cdot \nabla u +…
The immersed peridynamics (IPD) method is a fluid-structure interaction (FSI) model to simulate fluid-driven material damage and failure of an immersed structure, in which a peridynamic (PD) constitutive correspondence model is employed…
We propose first-order pressure-correction scheme for the incompressible Navier-Stokes equations, incorporating the recently developed the Dynamically Regularized Lagrange Multiplier (DRLM) methods. The resulting algorithms are fully…
We present a unified framework for the construction of localized exponential integrators that bypasses the traditional trade-off between the accuracy of global spectral methods and the efficiency of sparse finite differences. By evaluating…
Reynolds-Averaged Navier-Stokes (RANS) models are widely used for turbulent flow simulations due to their computational efficiency, but their accuracy strongly depends on the selected turbulence closure and may vary across the flow domain.…
In this paper, we introduce a Stancu-type generalization of multivariate neural network operators by incorporating two parameters that perturb the sampling nodes. The proposed operators extend the existing neural network operator by…
Two notable examples of dual functionals in approximation theory and computer-aided geometric design are the blossom and the divided difference operator. Both of these dual functionals satisfy a similar set of formulas and identities.…
This paper presents a new Monte Carlo (MC) algorithm for time-dependent particle transport problems with global variance reduction based on automatic weight windows (WWs). The centers of WWs at a time step are defined by the solution of an…
We consider an algebraic multigrid (AMG) scheme for the direct solution of complex- valued square linear systems based on a recursive 2 x 2 block partitioning of the coefficient matrix and study the optimal choices of its components. In…
Stochastic variance reduced gradient (SVRG) is an accelerated version of stochastic gradient descent based on variance reduction, and is promising for solving large-scale inverse problems. In this work, we analyze SVRG and a regularized…
In this paper, we study an evolution equation that involves a half-Laplacian operator derived from the Riesz fractional Laplacian, combined with a differential operator \(\mathcal{L}\). Using the identity $(-\Delta)^{1/2}=\mathcal…
We present new finite elements for solving the Riesz maps of the de Rham complex on triangular and tetrahedral meshes at high order. The finite elements discretize the same spaces as usual, but with different basis functions, so that the…
We are interested in the simulation of open quantum systems governed by the Lindblad master equation in an infinite-dimensional Hilbert space. To simulate the solution of this equation, the standard approach involves two sequential…
We show that, under certain circumstances, it is possible to automatically compute Jacobian-inverse-vector and Jacobian-inverse-transpose-vector products about as efficiently as Jacobian-vector and Jacobian-transpose-vector products. The…
We derive an algorithm for compression of the currents and varifolds representations of shapes, using ridge leverage score (RLS) sampling, and the theory of Nystrom approximation in Reproducing Kernel Hilbert Spaces. Our method is faster…
In this paper analysis is performed on a computational method for thermal radiative transfer (TRT) problems based on the multilevel quasidiffusion (variable Eddington factor) method with the method of long characteristics (ray tracing) for…
In this paper, new implicit methods with reduced memory are developed for solving the time-dependent Boltzmann transport equation (BTE). One-group transport problems in 1D slab geometry are considered. The reduced-memory methods are…
In this paper we present a multilevel projection-based iterative scheme for solving thermal radiative transfer problems that performs iteration cycles on the high-order Boltzmann transport equation (BTE) and low-order moment equations.…
Two combined methods for computing solutions of time-varying semilinear differential-algebraic equations (descriptor systems) are obtained. When constructing the methods, time-varying spectral projectors which can be found numerically are…