Related papers: Two-dimensional Gauss--Jacobi Quadrature for Multi…
We propose new iterative methods for computing nontrivial extremal generalized singular values and vectors. The first method is a generalized Davidson-type algorithm and the second method employs a multidirectional subspace expansion…
The Lattice Boltzmann Method (LBM) has emerged as a powerful tool in computational fluid dynamics and material science. However, standard LBM formulation imposes some limitations on the applications of the method, particularly compressible…
Gaussian process is a theoretically appealing model for nonparametric analysis, but its computational cumbersomeness hinders its use in large scale and the existing reduced-rank solutions are usually heuristic. In this work, we propose a…
Simulations of the discrete Boltzmann Bhatnagar-Gross-Krook (BGK) equation are an important tool for understanding fluid dynamics in non-continuum regimes. Here, we introduce a discontinuous Galerkin finite element method (DG-FEM) for…
We present an immersed interface method for the vorticity-velocity form of the 2D Navier Stokes equations that directly addresses challenges posed by multiply connected domains, nonconvex obstacles, and the calculation of force…
We present a simple, parallel and distributed algorithm for setting up and partitioning a sparse representation of a regular discretized simulation domain. This method is scalable for a large number of processes even for complex geometries…
In this paper we propose a novel numerical approach for the Boltzmann equation with uncertainties. The method combines the efficiency of classical direct simulation Monte Carlo (DSMC) schemes in the phase space together with the accuracy of…
A family of weak Galerkin finite element discretization is developed for solving the coupled Darcy-Stokes equation. The equation in consideration admits the Beaver-Joseph-Saffman condition on the interface. By using the weak Galerkin…
Efficient simulation of nonlinear and dispersive free-surface flows governed by the incompressible Navier-Stokes equations remains a central challenge in ocean and coastal engineering. The computational bottleneck arises from solving a…
Algorithms for computing the classical Gaussian quadrature rules (Gauss--Jacobi, Gauss--Laguerre, and Gauss--Hermite) are presented, based on globally convergent fourth-order iterative methods combined with asymptotic approximations, which…
We present a systematic computational framework for generating positive quadrature rules in multiple dimensions on general geometries. A direct moment-matching formulation that enforces exact integration on polynomial subspaces yields…
We develop an efficient algorithm for a spatially inhomogeneous matrix-valued quantum Boltzmann equation derived from the Hubbard model. The distribution functions are $2 \times 2$ matrix-valued to accommodate the spin degree of freedom,…
We develop a new spatial semidiscrete multiscale method based upon the edge multiscale methods to solve semilinear parabolic problems with heterogeneous coefficients and smooth initial data. This method allows for a cheap spatial…
We present and analyze in a unified setting two schemes for the numerical discretization of a Darcy-Forchheimer fluid flow model coupled with an advection-diffusion equation modeling the temperature distribution in the fluid. The first…
This work focuses on the development of efficient solvers for the pseudo-stress formulation of the unsteady Stokes problem, discretised by means of a discontinuous Galerkin method on polytopal grids (PolyDG). The introduction of the…
We develop stochastic mixed finite element methods for spatially adaptive simulations of fluid-structure interactions when subject to thermal fluctuations. To account for thermal fluctuations, we introduce a discrete fluctuation-dissipation…
We consider a mono-dimensional two-velocities scheme used to approximate the solutions of a scalar hyperbolic conservative partial differential equation. We prove the convergence of the discrete solution toward the unique entropy solution…
We present and analyze a new space-time parallel multigrid method for parabolic equations. The method is based on arbitrarily high order discontinuous Galerkin discretizations in time, and a finite element discretization in space. The key…
The numerical approximation of some Boussinesq systems in two spatial dimensions is here considered. The differential systems under study are proposed as asymptotic models for the propagation of waves along the interface of two layers of…
This paper presents an analysis of properties of two hybrid discretization methods for Gaussian derivatives, based on convolutions with either the normalized sampled Gaussian kernel or the integrated Gaussian kernel followed by central…