No Arabic abstract
We evaluate the effect of boundary layer losses on two-dimensional H2/O2/Ar cellular detonations obtained in narrow channels. The experiments provide the details of the cellular structure and the detonation speed deficits from the ideal CJ speed. We model the effect of the boundary layer losses by incorporating the flow divergence in the third dimension due to the negative boundary layer displacement thickness, modeled using Mirels theory. The cellular structures obtained numerically with the resulting quasi-2D formulation of the reactive Euler equations with two-step chain-branching chemistry are found in excellent agreement with experiment, both in terms of cell dynamics and velocity deficits, provided the boundary layer constant of Mirels is modified by a factor of 2. A significant increase in the cell size is found with increasing velocity deficit. This is found to be very well captured by the induction zone increase in slower detonations due to the lower temperatures in the induction zone.
Numerical work on shockwave/boundary-layer interactions (SBLIs) to date has largely focused on span-periodic quasi-2D configurations that neglect the influence lateral confinement has on the core flow. The present study is concerned with the effect of flow confinement on Mach 2 laminar SBLI in rectangular ducts. An oblique shock generated by a 2 degree wedge forms a conical swept SBLI with sidewall boundary layers before reflecting from the bottom wall of the domain. Multiple large regions of flow-reversal are observed on the sidewalls, bottom wall and at the corner intersection. The main interaction is found to be strongly three-dimensional and highly dependent on the geometry of the duct. Comparison to quasi-2D span-periodic simulations showed sidewalls strengthen the interaction by 31% for the baseline configuration with an aspect ratio of one. The length of the shock generator and subsequent trailing edge expansion fan position was shown to be a critical parameter in determining the central separation length. By shortening the length of the shock generator, control of the interaction and suppression of the central interaction is demonstrated. Parametric studies of shock strength and duct aspect ratio were performed to find limiting behaviours. For the largest aspect ratio of four, three-dimensionality was visible across 30% of the span width away from the wall. Topological features of the three-dimensional separation are identified and shown to be consistent with `owl-like separations of the first kind. The reflection of the initial conical swept SBLIs is found to be the most significant factor determining the flow structures downstream of the main interaction.
We perform direct numerical simulations of rotating Rayleigh--Benard convection of fluids with low ($Pr=0.1$) and high ($Pr=5$) Prandtl numbers in a horizontally periodic layer with no-slip top and bottom boundaries. At both Prandtl numbers, we demonstrate the presence of an upscale transfer of kinetic energy that leads to the development of domain-filling vortical structures. Sufficiently strong buoyant forcing and rotation foster the quasi-two-dimensional turbulent state of the flow, despite the formation of plume-like vertical disturbances promoted by so-called Ekman pumping from the viscous boundary layer.
This fluid dynamics video submitted to the Gallery of Fluid motion shows a turbulent boundary layer developing under a 5 metre-long flat plate towed through water. A stationary imaging system provides a unique view of the developing boundary layer as it would form over the hull of a ship or fuselage of an aircraft. The towed plate permits visualisation of the zero-pressure-gradient turbulent boundary layer as it develops from the trip to a high Reynolds number state ($Re_tau approx 3000$). An evolving large-scale coherent structure will appear almost stationary in this frame of reference. The visualisations provide an unique view of the evolution of fundamental processes in the boundary layer (such as interfacial bulging, entrainment, vortical motions, etc.). In the more traditional laboratory frame of reference, in which fluid passes over a stationary body, it is difficult to observe the full evolution and lifetime of turbulent coherent structures. An equivalent experiment in a wind/water-tunnel would require a camera and laser that moves with the flow, effectively `chasing eddies as they advect downstream.
The present study addresses the reaction zone structure and burning mechanism of unstable detonations. Experiments investigated mainly two-dimensional methane-oxygen cellular detonations in a thin channel geometry. The sufficiently high temporal resolution permitted to determine the PDF of the shock distribution, a power-law with an exponent of -3, and the burning rate of unreacted pockets from their edges - through surface turbulent flames with a speed approximately 3-7 times larger than the laminar one at the local conditions. Numerical simulations were performed using a novel Large Eddy Simulation method where the reactions due to both auto-ignition and turbulent transport and treated exactly at the sub-grid scale in a reaction-diffusion formulation. The model is an extension of Kerstein & Menons Linear Eddy Model for Large Eddy Simulation to treat flows with shock waves and rapid gasdynamic transients. The two-dimensional simulations recovered well the amplification of the laminar flame speed owing to the turbulence generated mainly by the shear layers originating from the triple points and subsequent Richtmyer-Meshkov instability associated with the internal pressure waves. The simulations clarified how the level of turbulence generated controlled the burning rate of the pockets, the hydrodynamic thickness of the wave, the cellular structure and its distribution. Three-dimensional simulations were found in general good agreement with the two-dimensional ones, in that the sub-grid scale model captured the ensuing turbulent burning once the scales associated with the cellular dynamics, where turbulent kinetic energy is injected, are well resolved.
We investigate the structures of the near-plate velocity and temperature profiles at different horizontal positions along the conducting bottom (and top) plate of a Rayleigh-B{e}nard convection cell, using two-dimensional (2D) numerical data obtained at the Rayleigh number Ra=10^8 and the Prandtl number Pr=4.4 of an Oberbeck-Boussinesq flow with constant material parameters. The results show that most of the time, and for both velocity and temperature, the instantaneous profiles scaled by the dynamical frame method [Q. Zhou and K.-Q. Xia, Phys. Rev. Lett. 104, 104301 (2010) agree well with the classical Prandtl-Blasius laminar boundary layer (BL) profiles. Therefore, when averaging in the dynamical reference frames, which fluctuate with the respective instantaneous kinematic and thermal BL thicknesses, the obtained mean velocity and temperature profiles are also of Prandtl-Blasius type for nearly all horizontal positions. We further show that in certain situations the traditional definitions based on the time-averaged profiles can lead to unphysical BL thicknesses, while the dynamical method also in such cases can provide a well-defined BL thickness for both the kinematic and the thermal BLs.