Related papers: FC-based shock-dynamics solver with neural-network…
We address a numerical instability that arises in the directionally split computation of hydrodynamic flows when shock fronts are parallel to a grid plane. Transverse oscillations in pressure, density and temperature are produced that are…
We propose a high order discontinuous Galerkin (DG) method for solving nonlinear Fokker-Planck equations with a gradient flow structure. For some of these models it is known that the transient solutions converge to steady-states when time…
Various forms of numerical shock instabilities are known to plague many contact and shear preserving approximate Riemann solvers, including the popular Harten-Lax-van Leer with Contact (HLLC) scheme, during high speed flow simulations. In…
The Boltzmann equation, a fundamental model in kinetic theory, describes the evolution of particle distribution functions through a nonlinear, high-dimensional collision operator. However, its numerical solution remains computationally…
Preserving scalar boundedness is important for numerical schemes used in turbulent compressible multi-component flow simulations to prevent unphysical results and unstable simulations. However, ensuring scalar boundedness for high-order,…
This paper proposes an approach, Spectral Dynamics Embedding Control (SDEC), to optimal control for nonlinear stochastic systems. This method reveals an infinite-dimensional feature representation induced by the system's nonlinear…
Smoothed Particle Hydrodynamics (SPH) is a popular numerical technique developed for simulating complex fluid flows. Among its key ingredients is the use of nonlocal integral relaxations to local differentiations. Mathematical analysis of…
Accurately characterizing non-linear functional manifolds with singularities is a fundamental challenge in scientific computing. While Multi-Layer Perceptrons (MLPs) dominate, their spectral bias hinders resolving high-curvature features…
This paper presents a novel structure-preserving scheme for Euler equations, focusing on the numerical conservation of entropy and kinetic energy. Explicit flux functions engineered to conserve entropy are introduced within the…
We introduce an accurate and robust technique for accessing causality of network transfer functions given in the form of bandlimited discrete frequency responses. These transfer functions are commonly used to represent the electrical…
In this work, we are concerned with the Fokker-Planck equations associated with the Nonlinear Noisy Leaky Integrate-and-Fire model for neuron networks. Due to the jump mechanism at the microscopic level, such Fokker-Planck equations are…
A fully discrete Crank--Nicolson Leap--Frog (CNLF) scheme is proposed and analyzed for the unsteady bioconvection flow problem with concentration-dependent viscosity. Spatial discretization is handled via the Galerkin finite element method…
In this work it is shown how the immersed boundary method of (Peskin2002) for modeling flexible structures immersed in a fluid can be extended to include thermal fluctuations. A stochastic numerical method is proposed which deals with…
Efficient and accurate spectral solvers for nonlocal models in any spatial dimension are presented. The approach we pursue is based on the Fourier multipliers of nonlocal Laplace operators introduced in a previous work. It is demonstrated…
In smooth-particle hydrodynamics (SPH), artificial viscosity is necessary for the correct treatment of shocks, but often generates unwanted dissipation away from shocks. We present a novel method of controlling the amount of artificial…
Physics-informed neural networks (PINNs) have great potential for flexibility and effectiveness in forward modeling and inversion of seismic waves. However, coordinate-based neural networks (NNs) commonly suffer from the "spectral bias"…
A novel algorithm for the direct numerical simulation of the variable-density, low-Mach Navier-Stokes equations extending the method of Kim, Moin, and Moser (1987) for incompressible flow is presented here. A Fourier representation is…
Hyperbolic conservation laws are conventionally solved by evolving reconstructed floating-point fields, incurring both computational overhead and structural diffusion near discontinuities. Here we introduce the Fast Quantised Numerical…
We solve the convection-diffusion equation using a coupling of cell-centered finite volume (FV) and discontinuous Galerkin (DG) methods. The domain is divided into disjoint regions assigned to FV or DG, and the two methods are coupled…
We propose a new entropy-compatible neural network method for scalar hyperbolic conservation laws and establish, to our knowledge, the first explicit \(L^1\) convergence rates in this setting that apply to piecewise smooth entropy…