Related papers: A Sampling-Based Adaptive Rank Approach to the Wig…
In this paper, we propose a numerical method for solving weakly compressible fluid flow based on a dynamical low-rank projector splitting. The low-rank splitting scheme is applied to the Boltzmann equation with BGK collision term, which…
In this paper, we propose and analyze first-order time-stepping pressure-correction projection scheme for the Navier-Stokes-Planck-Nernst-Poisson equations. By introducing a governing equation for the auxiliary variable through the ionic…
For the problems of low-rank matrix completion, the efficiency of the widely-used nuclear norm technique may be challenged under many circumstances, especially when certain basis coefficients are fixed, for example, the low-rank correlation…
In this paper, we construct a robust adaptive central-upwind scheme on unstructured triangular grids for two-dimensional shallow water equations with variable density. The method is well-balanced, positivity-preserving, and oscillation-free…
We generalize the interpolative separable density fitting (ISDF) method, used for compressing the four-index electron repulsion integral (ERI) tensor, to incorporate adaptive real space grids for potentially highly localized single-particle…
In this paper, we study the polynomial optimization problem of multi-forms over the intersection of the multi-spheres and the nonnegative orthants. This class of problems is NP-hard in general, and includes the problem of finding the best…
We present a simple model of alignment of a large number of rigid bodies (modeled by rotation matrices) subject to internal rotational noise. The numerical simulations exhibit a phenomenon of first order phase transition with respect the…
We present a comprehensive discussion of the so-called asymmetric Wigner bilayer system, where mobile point charges, all of the same sign, are immersed into the space left between two parallel, homogeneously charged plates (with possibly…
Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase space due to their better scalability than continuum approaches with respect to dimension. Complex processes…
We consider a coupled system of Navier-Stokes and Nernst-Planck equations, describing the evolution of the velocity and the concentration fields of dissolved constituents in an electrolyte solution. Motivated by recent applications in the…
We present a fully iterative adaptive algorithm for the numerical minimization of strongly convex energy functionals in Hilbert spaces. The proposed approach, which we first present in abstract form, generates a hierarchical sequence of…
We develop new adaptive alternative weighted essentially non-oscillatory (A-WENO) schemes for hyperbolic systems of conservation laws. The new schemes employ the recently proposed local characteristic decomposition based central-upwind…
In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker-Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source.…
Randomized numerical linear algebra is proved to bridge theoretical advancements to offer scalable solutions for approximating tensor decomposition. This paper introduces fast randomized algorithms for solving the fixed Tucker-rank problem…
The existence of large-data weak solutions to a steady compressible Navier-Stokes-Fourier system for chemically reacting fluid mixtures is proved. General free energies are considered satisfying some structural assumptions, with a pressure…
In this work, we derive a reliable and efficient residual-typed error estimator for the finite element approximation of a 2d cathodic protection problem governed by a steady-state diffusion equation with a nonlinear boundary condition. We…
We present a dynamical low-rank (DLR) method for the Vlasov-Poisson-Fokker-Planck (VPFP) system. Our main contributions are two-fold: (i) a conservative spatial discretization of the Fokker-Planck operator that factors into velocity-only…
The low-rank approximation is a complexity reduction technique to approximate a tensor or a matrix with a reduced rank, which has been applied to the simulation of high dimensional problems to reduce the memory required and computational…
Efficient sampling in biomolecular simulations is critical for accurately capturing the complex dynamical behaviors of biological systems. Adaptive sampling techniques aim to improve efficiency by focusing computational resources on the…
In the field of model compression, choosing an appropriate rank for tensor decomposition is pivotal for balancing model compression rate and efficiency. However, this selection, whether done manually or through optimization-based automatic…