Shock Mach Number Research Articles

Self-organization of collisionless shocks: from a laminar profile to a rippled time-dependent structure

A collisionless shock structure results from the nonlinear interaction between charged particles and electromagnetic fields. Yet, a collisionless shock is globally governed by the mass, momentum and energy conservation requirements. A stable shock structure must ensure that the fluxes of the conserved quantities are constant on average, and, therefore, is determined by this necessity. Here, we study an observed high upstream temperature high Mach number shock and show that the conservation laws cannot be fulfilled unless the shock is spatially inhomogeneous along the shock front and time-dependent.

Electron Acceleration at Quasi-parallel Nonrelativistic Shocks: A 1D Kinetic Survey

We present a survey of 1D kinetic particle-in-cell simulations of quasi-parallel nonrelativistic shocks to identify the environments favorable for electron acceleration. We explore an unprecedented range of shock speeds v sh ≈ 0.067–0.267c, Alfvén Mach numbers MA=5–40 , sonic Mach numbers Ms=5–160 , as well as the proton-to-electron mass ratios m i/m e = 16–1836. We find that high Alfvén Mach number shocks can channel a large fraction of their kinetic energy into nonthermal particles, self-sustaining magnetic turbulence and acceleration to larger and larger energies. The fraction of injected particles is ≲0.5% for electrons and ≈1% for protons, and the corresponding energy efficiencies are ≲2% and ≈10%, respectively. The extent of the nonthermal tail is sensitive to the Alfvén Mach number; when MA≲10 , the nonthermal electron distribution exhibits minimal growth beyond the average momentum of the downstream thermal protons, independently of the proton-to-electron mass ratio. Acceleration is slow for shocks with low sonic Mach numbers, yet nonthermal electrons still achieve momenta exceeding the downstream thermal proton momentum when the shock Alfvén Mach number is large enough. We provide simulation-based parameterizations of the transition from thermal to nonthermal distribution in the downstream (found at a momentum around pi,e/mivsh≈3mi,e/mi ), as well as the ratio of nonthermal electron to proton number density. The results are applicable to many different environments and are important for modeling shock-powered nonthermal radiation.

Numerical study of bifurcation structures in reflected shock-wave/laminar-boundary-layer interaction within an end-wall tube

Abstract The primary aim of this study is to analyze the unsteady characteristics of the interaction between a reflected shock wave and a laminar boundary layer in an end-wall shock tube. Our direct numerical simulations at shock Mach numbers of Ms=1.9, 2.5, and 3.5 using a fifth-order WENO scheme and three-step Runge-Kutta time integration method revealed inhomogeneity and anisotropy in the shock bifurcation. Surprisingly, the upper and lower bifurcated structures maintain a notably asymmetric flow during the forward propagation of the reflected shock bifurcation. The inverse flow in the bifurcation resembles a crooked earthworm structure, exhibiting high-frequency oscillations indicative of instability. However, at higher shock intensities, the earthworm transforms into a stable strip-like configuration, facilitating the entrapment of inverse flow and leading to rapid bifurcation height growth and early convergence. Additionally, isolated islands with high density, temperature, and pressure emerge in the transitional region behind the bifurcated shocks, due to variations in wave propagation speed.

Computational Study of Shocked V-Shaped N2/SF6 Interface across Varying Mach Numbers

The Mach number effect on the Richtmyer–Meshkov instability (RMI) evolution of the shocked V-shaped N2/SF6 interface is numerically studied in this research. Four distinct Mach numbers are taken into consideration for this purpose: Ms=1.12,1.22,1.42, and 1.62. A two-dimensional space of compressible two-component Euler equations is simulated using a high-order modal discontinuous Galerkin approach to computational simulations. The numerical results show good consistency when compared to the available experimental data. The computational results show that the RMI evolution in the shocked V-shaped N2/SF6 interface is critically dependent on the Mach number. The flow field, interface deformation, intricate wave patterns, inward jet development, and vorticity generation are all strongly impacted by the shock Mach number. As the Mach number increases, the V-shaped interface deforms differently, and the distance between the Mach stem and the triple points varies depending on the Mach number. Compared to lower Mach numbers, higher ones produce larger rolled-up vortex chains. A thorough analysis of the Mach number effect identifies the factors that propel the creation of vorticity during the interaction phase. Moreover, kinetic energy and enstrophy both dramatically rise with increasing Mach number. Lastly, a detailed analysis is carried out to determine how the Mach number affects the temporal variations in the V-shaped interface’s features.

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. Steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

Comprehensive analysis of normal shock wave propagation in turbulent non-ideal gas flows with analytical and neural network methods

Shock wave propagation in gases through turbulent flow has wide-reaching implications for both theoretical research and practical applications, including aerospace engineering, propulsion systems, and industrial gas processes. The study of normal shock propagation in turbulent flow over non-ideal gas investigates the changes in pressure, density, and flow velocity across the shock wave. The Mach number is derived for the system and explored across various gas molecule quantities and turbulence intensities. This study analytically investigated the normal shock wave propagation in turbulent flow of adiabatic gases with modified Rankine–Hugoniot conditions. Artificial neural network (ANN) techniques are used to estimate the solutions for shock strength and Mach number training validation phases of back-propagated neural networks with the Levenberg–Marquardt algorithm. The results reveal that pressure ratio with density ratio increase for higher values of increase in the turbulence level as well as intermolecular forces. A reverse trend is observed in velocity coefficient after shock in the presence of adiabatic gas. The regression coefficient values obtained using the network model ranged from 0.999 99 to 1, indicating an almost perfect correlation. These findings demonstrate that the ANN can predict the Mach number with high accuracy.

Numerical investigation of the reflection of unsteady decelerating incident shock waves

The incident shock attenuation phenomenon in shock tube has received widespread interest because of its inevitable influence on experimental gas properties. However, few studies have investigated the cascading effects of the resulting nonuniformity on shock reflection in shock tunnels before the nozzle. This paper describes a numerical study on the unsteady reflection of a decelerating incident shock driven by a decelerating piston, as a simplification of the complex nonideal factors. The initial decay and uniform parameter distributions are generated based on the non-inertial frame. The results indicate that the existence of the nonuniform area forces the reflected shock wave region to undergo a transition process of attenuation and then stability. The final values of the gas parameters in zone 5 will, therefore, deviate from those given by the traditional relation for an ideal shock tube, which are only dependent on the terminal shock Mach number. This affects the determination of the total temperature T5, which is difficult to measure directly. We discuss the reconstruction of the nonuniform region with the attenuation trajectory of the incident shock wave and find that there is no one-to-one correspondence between this trajectory and the resulting nonuniformity, which introduces additional uncertainties to the predictions. Thus, an analogy method is developed through the multilevel block building algorithm, allowing the total temperature T5 to be determined using the total pressure p5 considering the effect of nonuniformity. Further assessments verify the applicability of this method in cases with multidimensional and viscous interactions.

Analyzing Richtmyer–Meshkov Phenomena Triggered by Forward-Triangular Light Gas Bubbles: A Numerical Perspective

In this paper, we present a numerical investigation into elucidating the complex dynamics of Richtmyer–Meshkov (RM) phenomena initiated by the interaction of shock waves with forward-triangular light gas bubbles. The triangular bubble is filled with neon, helium, or hydrogen gas, and is surrounded by nitrogen gas. Three different shock Mach numbers are considered: Ms=1.12,1.21, and 1.41. For the numerical simulations, a two-dimensional system of compressible Euler equations for two-component gas flows is solved by utilizing the high-fidelity explicit modal discontinuous Galerkin technique. For validation, the numerical results are compared with the existing experimental results and are found to be in good agreement. The numerical model explores the impact of the Atwood number on the underlying mechanisms of the shock-induced forward-triangle bubble, encompassing aspects such as flow evolution, wave characteristics, jet formation, generation of vorticity, interface features, and integral diagnostics. Furthermore, the impacts of shock strengths and positive Atwood numbers on the flow evolution are also analyzed. Insights gained from this numerical perspective enhance our understanding of RM phenomena triggered by forward-triangular light gas bubbles, with implications for diverse applications in engineering, astrophysics, and fusion research.

Hypersonic inlet flow field reconstruction dominated by shock wave and boundary layer based on small sample physics-informed neural networks

High Mach number inlets face complex challenges such as shock wave/boundary layer interference, adversely impacting the aerodynamic characteristics and stable operating boundaries. Traditional computational fluid dynamics (CFD) simulations for performance and flow field analysis are slow, and data-driven deep learning technologies struggle with predicting flow fields in complex aerodynamic phenomena like shock wave/boundary layer interactions, boundary layer separation, and Mach disks. This study introduces a rapid, robust method for reconstructing hypersonic viscous inlet flow fields using inlet geometry design parameters. This method is called shock wave and boundary layer dominated two-dimensional multiphysical field reconstruction model based on multi-scale sensory field fusion residual physical information neural network (PIMSRF_ResNet), ensuring high-precision training data. Using the optimal Pareto solution set from multi-objective optimization with a small sample intelligence method, 150 sets of multi-physics fields under Ma10 conditions with varying geometric design parameters are calculated. Compared to traditional pure data-driven models, PIMSRF_ResNet's structural similarity in the test set reaches 0.9773, with a peak signal-to-noise ratio close to 30 dB and a correlation coefficient exceeding 98%. These experimental results demonstrated that the proposed method is a promising tool for real-time prediction of complex internal flow fields dominated by high Mach number shock waves and boundary layers.

Quasi-perpendicular shocks of galaxy clusters in hybrid kinetic simulations

Context. Cosmic ray acceleration in galaxy clusters is still an ongoing puzzle, with relativistic electrons forming radio relics at merger shocks and emitting synchrotron radiation. These shocks are also potential sources of ultra-high-energy cosmic rays, gamma rays, and neutrinos. Our recent work focuses on electron acceleration at low Mach number merger shocks in the hot intracluster medium which is characterized by high plasma beta. Using particle-in-cell (PIC) simulations, we previously showed that electrons are energized through the stochastic shock-drift acceleration process, which is facilitated by multi-scale turbulence, including ion-scale shock surface rippling. For the present work, we performed hybrid-kinetic simulations in a range of various quasi-perpendicular foreshock conditions, including plasma beta, magnetic obliquity, and the shock Mach number. Aims. We study the ion kinetic physics, which is responsible for the shock structure and wave turbulence, that in turn affects the particle acceleration processes. We cover the spatial and temporal scales, which allow the development of large-scale ion turbulence modes in the system. Methods. We applied a recently developed generalized fluid-particle hybrid numerical code that can combine fluid modeling for both electrons and ions with an arbitrary number of kinetic species. We limited this model to a standard hybrid simulation configuration with kinetic ions and fluid electrons. The model utilizes the exact form of the generalized Ohm’s law, allowing for an arbitrary choice of mass and energy densities, as well as the charge-to-mass ratio of the kinetic species. Results. We show that the properties of ion-driven multi-scale magnetic turbulence in merger shocks are in agreement with the ion structures observed in PIC simulations. In typical shocks with the sonic Mach number Ms = 3, the magnetic structures and shock front density ripples grow and saturate at wavelengths reaching approximately four ion Larmor radii. Only shocks with Ms ≳ 2.3 develop ripples. At very weak shocks with Ms ≲ 2.3, weak turbulence is formed downstream of the shock. We observed a moderate dependence of the strength of magnetic field fluctuations on the quasi-perpendicular magnetic field obliquity. However, as the field obliquity decreases, the shock front ripples exhibit longer wavelengths. Finally, we note that the steady-state structure of Ms = 3 shocks in high-beta plasmas shows evidence that there is little difference between 2D and 3D simulations. The turbulence near the shock front seems to be a 2D-like structure in 3D simulations.

Numerical simulation of the interaction between a planar shock wave and a backward-facing triangular bubble containing gases with different Atwood numbers

The interaction between a shock wave and an interface delineating two gases engenders intricate flow physics, with particular attention drawn to the hydrodynamic instability due to its practical significance. Previous studies have primarily focused on elucidating different wave patterns and instabilities evolution at the interface during the initial phase of shock interaction with cylindrical or spherical bubbles. However, scant literature has shifted its focus toward exploring the long-term morphology of bubbles, especially those characterized by polygonal interfaces. Notably, the detailed examination of shock interaction with a polygonal interface, such as a triangular one with a constant incident angle, remains largely unexplored in existing literature. Recently, the longtime evolution of detailed flow structures across the interface of shock-forward-facing triangular bubbles was captured by Kundu et al. [“A study on dynamics of shock-accelerated forward-facing triangular bubbles at different Atwood numbers,” Phys. Fluids 36, 016110 (2024)] through numerical simulation. In this study, the dynamics of a shock-accelerated backward-facing triangular interface containing various gases, namely, Sulfur Hexafluoride, Refrigerant-22, Argon, Neon, and Helium, is studied for a shock Mach number of 1.21. Simulations were performed by solving the two-dimensional Euler equation using low-dissipative advection upwind splitting methods (AUSMD), in conjunction with a derived ninth-order upwind scheme and a four-stage third-order Runge–Kutta scheme for temporal integration. The development of Richtmyer–Meshkov (RM) and Kelvin–Helmholtz (K–H) instabilities at the interface, mixing, and normalized movements of backward-facing triangular bubbles is captured at different Atwood (At) numbers.

Formation and propagation characteristics of a weak shock wave in maglev tube

The propagation of the weak shock wave (WSW) to the tunnel exits and their radiation as micro-pressure waves (MPWs) may cause sonic booms or structural resonance of buildings, posing potential hazards to humans, animals, and buildings in the exit's environment. The characteristics of the WSW and sonic booms of a maglev train/tube coupling model were studied based on the two-dimensional axisymmetric unsteady Reynolds average Navier–Stokes turbulence model. In the later stage of a MPW, the formation mechanism, geometry, and kinematic characteristics of compressible vortex rings (CVRs) were systematically analyzed. The inertial effect causes the initial wavefront to gradually transition from a Gaussian-shape waveform to a triangular waveform during its propagation, eventually coalescing into a WSW. The overpressure, density jump, and shock Mach number at the WSW location all increase with the increasing train speed, while the WSW thickness decreases accordingly. The formation distance of the WSW is inversely proportional to the amplitude of the initial wavefront gradient, and the WSW directly causes the occurrence of the exit sonic boom. The MPW amplitude has significant directionality with a largest value in the axial direction. Within the speed range of 450–700 km/h, the sound pressure level of the MPW exceeds the hearing threshold and even reaches the feeling threshold. The evolution of CVRs includes primary CVR, secondary CVR, and Kelvin–Helmholtz vortices. Primary CVR has the greatest impact on the axial MPW among them. The occurrence of CVRs will cause a second small noise level other than the sonic boom.

A Study of Shock-Induced Transient Jets: The Impact of Exit Conditions on Mach Disk and Vortex Ring

Abstract We investigate the effect of exit pressure history on the flow characteristics of underexpanded transient jets. Using both experiments and numerical simulations, we study the dynamics of shock-cell and vortex structures within these jets. A shock tube with an open-ended configuration allows us to generate transient jets by adjusting the diaphragm pressure ratio and the length of the driver section. Our results indicate that when the shock Mach number exceeds 1.6, a Mach disk forms, indicating a highly underexpanded transient jet at the exit of the shock tube. A distinguishing feature of this jet is the emergence of counter-rotating vortex rings (CRVRs) alongside the initial primary vortex ring. Our findings reveal a substantial influence of both the amplitude and duration of the peak exit pressure on the characteristics of the Mach disk and vortex ring. Notably, the characteristics of the primary vortex ring exhibit significant sensitivity to the formation and evolution of CRVRs. In cases of continuously decreasing exit pressure, the Mach disk follows a consistent self-similar decay pattern, regardless of the peak exit pressure magnitude. Finally, we present an empirical relationship between exit pressure and the characteristics of the Mach disk. In summary, this research provides insight into the complex interaction between the exit pressure history and the flow characteristics in underexpanded transient jets.

Computational study of shock diffraction over convex edges

Numerical analysis of normal shock diffraction over various step corners is presented using an implicit second-order upwind least squares cell-based method. Slipstream angle calculation for various corner angles and Mach numbers, variation of pressure load and specific heat flux with incident shock Mach numbers and different step corners along the step wall, separation point movement, and reattachment of the streamlines are some of the key features of the present study. Detailed studies have been conducted about the perturbed region enclosed within the diffracted shock wave and the last running expansion wave with the increase in the incident shock Mach numbers. The finite Volume Method is used to solve the governing equations numerically for the simulation of the moving shock. Various flow characteristics such as secondary shock, contact surface, expansion wave, slipstream, vortex, etc are well captured. Apart from isopycnics; Mach contours, isobars, and isotherms are plotted as well. The reattachment region formed after the flow separation near the step edge is mentioned and corresponding lengths over various steps are presented. Separation point movement along the step wall is highlighted and a relevant increase in reattachment length is reported. The slipstream angles get increased with incident shock Mach numbers and with corner angles. However, the rate of increment of slipstream angle decreases gradually. The core of the vortices generated and the reattachment of the separated streamlines are indicated by the two downfalls in the heat flux values.

Parametric effects on Richtmyer–Meshkov instability of a V-shaped gaseous interface within linear stage

The Richtmyer–Meshkov instability of a V-shaped air/SF6 gaseous interface is numerically studied via a high-order finite difference scheme and a localized artificial diffusivity method. The oblique angle of the interface ranges from 20° to 75°, and the incident shock Mach number varies from 1.05 to 1.75. The wave patterns and the vortex structures are visualized during the interface evolution. A cavity is observed at the spike fingertip when the oblique angle decreases, which proves to be formed due to Mach reflection of the transmitted shock through velocity decomposition. By analyzing the linear growth rates of the interface, a modified empirical model for the reduction factor is suggested with model coefficients acquired by linear fitting for different Mach numbers. With shock polar analysis (SPA) method and visualization of the wave configuration, a criterion is proposed to explain the non-monotonic dependence of the linear growth rate on the oblique angle. In addition, Mach number effects on the linear growth rate are discussed by the SPA method, especially the anomalous behavior of the Mach 1.05 case.

A study on dynamics of shock-accelerated forward-facing triangular bubbles at different Atwood numbers

The complexity of flow physics and the associated hydrodynamic instability arising out of interactions of a shock wave with forward and backward-facing triangular interfaces drew the attention of researchers around the globe. In earlier studies, many researchers focused on the formation of different wave patterns, the development of instabilities at the interface, and the flow morphology during the initial phase of shock wave interacting with light and heavier bubbles. However, limited studies are available in the literature on the interaction of shock with a polygonal interface. Furthermore, it is difficult to capture the complex flow physics of a polygonal interface accelerated by shock waves at later time instants. In the present study, the dynamics of shock-accelerated forward-facing triangular interface containing various gases, namely, sulfur hexafluoride, refrigerant-22, argon, neon, and helium, are examined numerically for a longer time duration for a shock Mach number of 1.21. The simulations were performed by solving two-dimensional Euler equations using a low-dissipative advection upstream splitting method algorithm coupled with a derived ninth-order upwind scheme and a four-stage third-order Runge–Kutta scheme. The numerical results demonstrated the influence of the Atwood number on vorticity generation, bubble deformation, mixing, and the development of Kelvin Helmholtz instabilities on the bubble interface up to long instants, which are not available in the literature. The Fourier spectra of the streamwise kinetic energy showed the distribution of energy in the larger and smaller scale vortical structures.

A New Enigmatic Radio Relic in the Low-mass Cluster Abell 2108

ABSTRACT We report the discovery of a radio relic in the north-eastern periphery of the galaxy cluster Abell 2108 (A2108). A2108 is part of the uGMRT LOw-MAss Galaxy Cluster Survey (GLOMACS), where our main aim is to search for diffuse radio emission signatures in very sparsely explored low-mass galaxy clusters using uGMRT band-3 (central frequency 400 MHz). We used our uGMRT band-3 observations along with the existing archival band-3 uGMRT data to improve image sensitivity. Along with the previously reported south-western relic, the discovery of the new relic makes A2108 one of the few low-mass clusters hosting a double relic. The new relic spans over a region of 610 kpc × 310 kpc and, interestingly, differs considerably in size and morphology from the other relic. Using XMM–Newton science archive data, we also report the tentative detection of a mildly supersonic shock of Mach number $\mathcal {M}_\mathrm{SB}=1.42$ and $\mathcal {M}_\mathrm{T} = 1.43$ from the surface brightness and temperature discontinuities, respectively, near this newly found relic. Both the relics in A2108 are found to be significantly under-luminous compared with other double-relic systems in the mass–luminosity plane. Moreover, the north-eastern relic is remarkably under-luminous in the size–luminosity plane. Although mild supersonic shocks resulting from an off-axis merger could have influenced the origin of both relics, we hypothesize that local environments have played a crucial role in shaping their morphologies.

The dispersion and ignition behavior of zirconium particles under the effect of a shock wave

Experiments of the interaction process between shock waves and zirconium particles have been carried out using a self-designed shock tube. The interaction process was recorded using a high-speed shadowgraph imaging system, pressure testing system, and infrared thermal imager. The wave-systems structures resulting from the interaction between shock waves and accumulated zirconium powders have been characterized. Generally, the axial distance of the powder cloud reduces with the increase in accumulation thickness, which in turn increases with the addition of shock Mach number. The incident shock wave with an intensity of 1.53 Ma was enough to ignite zirconium powders under these experimental conditions. The ignition delay time for zirconium powders decreases with the addition of shock Mach number, while it decreases for thinner accumulations. Besides, the flame propagation characteristics and the ignition mechanism of zirconium powders are discussed in detail. The results can serve as a supplementary guide for the utilization and safe operation of zirconium powders.

Experimental investigation on initiation mechanism, overpressure, and flame propagation characteristics of methane-air mixtures explosion induced by hexogen in a closed pipeline

Gas explosion accidents under shock wave stimulation cause more severe hazards than electric spark ignition. To explore the initiation mechanism and explosion characteristics of methane induced by shock waves. A series of experimental investigations of methane-air mixtures initiated by hexogen (RDX) was performed based on a self-developed closed pipeline. The effects of incident shock wave intensity and methane concentration on the characteristics of methane explosion were studied. The results indicate that methane participating in the reaction can significantly increase the positive pressure duration compared with RDX explosion. With the increase of incident shock Mach number (Ma), methane deflagration pressure increases and counteracts the attenuation of RDX detonation pressure. When the ignition distance is 0 cm, the overpressure inside the pipeline is directly proportional to the methane concentration from 5 vol% to 13 vol%. When the concentration approaches the upper explosion limit, the flame velocity exceeds the theoretical detonation velocity, resulting in the unstable detonation of methane in the pipeline.

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths, vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn., vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids, vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. (Phys. Fluids, vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.