Direct numerical simulations are performed of gravity-current fronts in the lock-exchange configuration. The case of small density differences is considered, where the Boussinesq approximations can be adopted. The key objective of the investigation is a detailed analysis of the flow structure at the foremost part of the front, where no previous high-resolution data were available. For the simulations, high-order numerical methods are used, based on spectral and spectral-element discretizations and compact finite differences. A three-dimensional simulation is conducted of a front spreading along a no-slip boundary at a Reynolds number of about 750. The simulation exhibits all features typically observed in experimental flows near the gravity-current head, including the lobe-and-cleft structure at the leading edge. The results reveal that the flow topology at the head differs from what has been assumed previously, in that the foremost point is not a stagnation point in a translating system. Rather, the stagnation point is located below and slightly behind the foremost point in the vicinity of the wall. The relevance of this finding for the mechanism behind the lobe-and-cleft instability is discussed. In order to explore the high-Reynolds-number regime, and to assess potential Reynolds-number effects, two-dimensional simulations are conducted for Reynolds numbers up to about 30 000, for both no-slip and slip (i.e. shear-stress free) boundaries. It is shown that although quantitative Reynolds-number effects persist over the whole range examined, no qualitative changes in the flow structure at the head can be observed. A comparison of the two-dimensional results with laboratory data and the three-dimensional simulation provides evidence that a two-dimensional model is able to capture essential features of the flow at the head. The simulations also show that for the free-slip case the shape of the head agrees closely with the classical inviscid theory of Benjamin.
High-resolution simulations are presented of particle-driven gravity currents in the lock-exchange configuration. The study concentrates on dilute flows with small density differences between particle-laden and clear fluid. Moreover, particles are considered which have negligible inertia, and which are much smaller than the smallest length scales of the buoyancy-induced fluid motion. For the mathematical description of the particulate phase a Eulerian approach is employed with a transport equation for the local particlenumber density. The governing equations are integrated numerically with a high-order mixed spectral/ spectral-element technique. In the analysis of the results, special emphasis is placed on the sedimentation of particles and the influence of particle settling on the flow dynamics. Time-dependent sedimentation profiles at the channel floor are presented which agree closely with available experimental data. A detailed study is conducted of the balance between the various components of the energy budget of the flow, i.e. the potential and kinetic energy, and the dissipative losses. Furthermore, the simulation results, along with a modified Shields criterion, are used to show that resuspension of sediment back into the particle-driven current is unlikely to occur in the cases considered. Two-dimensional (2D) and three-dimensional (3D) computations are compared which reveal that, for the present configuration, a 2D model can predict reliably the flow development at early times. However, concerning the long-time evolution of the flow, more substantial differences exist between 2D and 3D simulations. Ó
Results are presented from a high-resolution computational study of particle-driven gravity currents in a plane channel. The investigation was conducted in order to obtain better insight into the energy budget and the mixing behaviour of such flows. Twoand three-dimensional simulations are discussed, and the effects of different factors influencing the flow are examined in detail. Among these are the aspect ratio of the initial suspension reservoir, the settling speed of the particles, and the initial level of turbulence in the suspension. While most of the study is concerned with the lockexchange configuration, where the initial height of the suspension layer is equal to the height of the channel, part of the analysis is also done for a deeply submerged case. Here, the suspension layer is only one-tenth of the full channel height. Concerning the energy budget, a careful analysis is undertaken of dissipative losses in the flow. Dissipative losses arising from the macroscopic fluid motion are distinguished from those due to the microscopic flow around each sedimenting particle. It is found that over a large range of settling velocities and suspension reservoir aspect ratios, sedimentation accounts for roughly half of all dissipative losses. The analysis of the mixing behaviour of the flow concentrates on the mixing between interstitial and ambient fluid, which are taken to be of identical density. The former is assumed to carry a passive contaminant, whose dispersion with time is analysed qualitatively and quantitatively by means of Lagrangian markers. The simulations show the mixing between interstitial and ambient fluid to be more intense for larger values of the particle settling velocity. Finally, the question is addressed of whether or not initial turbulence in the suspension may exert a significant effect on the flow evolution. To this end, results from three simulations with widely different levels of initial kinetic energy are compared. While the initial turbulence level strongly affects the mixing within the current, it has only a small influence on the front velocity and the overall sedimentation rate.
Direct numerical simulation databases of turbulent channel and pipe flow have been used in order to assess the energy transfer between resolved and unresolved motions in large-eddy simulations. To this end, the velocity fields are split into three parts: a statistically stationary mean flow, the resolved, and the unresolved turbulent fluctuations. The distinction between the resolved and unresolved motions is based on the application of a cutoff filter in spectral space. Within the buffer layer a backward transfer of averaged kinetic energy from subgrid to grid-scale turbulent motions has been found to exist, which is primarily caused by subgrid-scale stresses aligned with the mean rates of strain. Such reverse transfer generally cannot be described by the simple eddy-viscosity-type subgrid models usually applied in large-eddy simulations. The use of a conditional averaging technique revealed that the reverse transfer of energy within the near-wall flow is strongly enhanced by coherent motions, such as the well-known bursting events.
Results are presented from a linear-stability analysis of the flow at the head of two-dimensional gravity-current fronts. The analysis was undertaken in order to clarify the instability mechanism that leads to the formation of the complex lobe-and-cleft pattern which is commonly observed at the leading edge of gravity currents propagating along solid boundaries. The stability analysis concentrates on the foremost part of the front, and is based on direct numerical simulation data of two-dimensional lock-exchange flows which are described in the companion paper, Härtel et al. (2000). High-order compact finite differences are employed to discretize the stability equations which results in an algebraic eigenvalue problem for the amplification rate, that is solved in an iterative fashion. The analysis reveals the existence of a vigorous linear instability that acts in a localized way at the leading edge of the front and originates in an unstable stratification in the flow region between the nose and stagnation point. It is shown that the amplification rate of this instability as well as its spanwise length scale depend strongly on Reynolds number. For validation, three-dimensional direct numerical simulations of the early stages of the frontal instability are performed, and close agreement with the results from the linear-stability analysis is demonstrated.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.