Abstract:This study presents two-dimensional direct numerical simulations for sediment-laden current with higher density propagating forward through a lighter ambient water. The incompressible NavierStokes equations including the buoyancy force for the density difference between the light and heavy fluids are solved by a finite difference scheme based on a structured mesh. The concentration transport equations are used to explore such rich transport phenomena as gravity and turbidity currents. Within the framework of an Upwinding Combined Compact finite Difference (UCCD) scheme, rigorous determination of weighting coefficients underlies the modified equation analysis and the minimization of the numerical modified wavenumber. This sixth-order UCCD scheme is implemented in a four-point grid stencil to approximate advection and diffusion terms in the concentration transport equations and the first-order derivative terms in the Navier-Stokes equations, which can greatly enhance convective stability and increase dispersive accuracy at the same time. The initial discontinuous concentration field is smoothed by solving a newly proposed Heaviside function to prevent numerical instabilities and unreasonable concentration values. A two-step projection method is then applied to obtain the velocity field. The numerical algorithm shows a satisfying ability to capture the generation, development, and dissipation of the Kelvin-Helmholz instabilities and turbulent billows at the interface between the current and the ambient fluid. The simulation results also are compared with the data in published literatures and good agreements are found to prove that the present numerical model can well reproduce the propagation, particle deposition, and mixing processes of lock-exchange gravity and turbidity currents.

Simulation of Kelvin-Helmholtz instability with high-order finite difference method

A numerical method for two-phase flow with its application to the Kelvin Helmholtz instability problem

Simulation of Kelvin-Helmholtz Instability with Flux-Corrected Transport Method

High-fidelity Simulation of Regular Waves Based on Multi-Moment Finite Volume Formulation and THINC Method

Development of a Volume of Fluid Method for Computing Interfacial Incompressible Fluid Flows

Numerical Modeling of Lock-Exchange Gravity/turbidity Currents by a High-Order Upwinding Combined Compact Difference Scheme

Velocity Gradient in Kelvin-Helmholtz Instability for Supersonic Fluid

The influence of the numerical viscosity of high-order weighted essentially non-oscillatory (WENO) schemes on simulation of Rayleigh-Taylor instability

Study on the Kelvin-Helmholtz Instability in Two-Dimensional Incompressible Fluid

A Weakly Nonlinear Model for Kelvin–Helmholtz Instability in Incompressible Fluids

Lattice Boltzmann Study On Kelvin-Helmholtz Instability: Roles Of Velocity And Density Gradients

Weakly Nonlinear Analysis on the Kelvin-Helmholtz Instability

Numerical simulation and mutual conversion of the temporal-spatial modes of hydrodynamic stability in boundary layer flow

Numerical Investigations of the Kelvin-Helmholtz Instability Based on a Multi-Resolution Moving Particle Semi-Implicit Method

A Well-Balanced Numerical Model for the Simulation of Long Waves over Complex Domains

Nonlinear Simulations of Coalescence Instability Using a Flux Difference Splitting Method

Weakly nonlinear analysis on the Kelvin-Helmholtz instability (vol 86, pg 15002, 2009)

Numerical Investigation on the Ablative Kelvin-Helmholtz Instability

A Newly Improved WENO Scheme and Its Application to the Simulation of Richtmyer-Meshkov Instability

Most-unstable Wave Method for Inlet Conditions of Large Eddy Simulations of Spatially Developing Channel Flows

NUMERICAL SIMULATIONS OF RICHTMYER-MESHKOV INSTABILITY