An explicit high-order single-stage single-step positivity-preserving finite difference WENO method for the compressible Euler equations

In this work we construct a high-order, single-stage, single-step positivity-preserving method for the compressible Euler equations. Space is discretized with the finite difference weighted essentially non-oscillatory (WENO) method. Time is discretized through a Lax-Wendroff procedure that is constructed from the Picard integral formulation (PIF) of the partial differential equation. The method can be viewed as a modified flux approach, where a linear combination of a low- and high-order flux defines the numerical flux used for a single-step update. The coefficients of the linear combination are constructed by solving a simple optimization problem at each time step. The high-order flux itself is constructed through the use of Taylor series and the Cauchy-Kowalewski procedure that incorporates higher-order terms. Numerical results in one- and two-dimensions are presented

A positivity-preserving scheme for fluctuating hydrodynamics

A finite-difference hybrid numerical method for the solution of the isothermal fluctuating hydrodynamic equations is proposed. The primary focus is to ensure the positivity-preserving property of the numerical scheme, which is critical for its functionality and reliability especially when simulating fluctuating vapour systems. Both cases of single- and two-phase flows are considered by exploiting the van der Waals' square-gradient approximation to model the fluid (often referred to as “diffuse-interface” model). The accuracy and robustness of the proposed scheme is verified against several benchmark theoretical predictions for the statistical properties of density, velocity fluctuations and liquid-vapour interface, including the static structure factor of the density field and the spectrum of the capillary waves excited by thermal fluctuations at interface. Finally, the hybrid scheme is applied to the challenging bubble nucleation process, and is shown to capture the salient features of the phenomenon, namely nucleation rate and subsequent bubble-growth dynamics