The presence of dispersed-phase droplets can result in a notable increase in a system's drag. However, our understanding of the mechanism underlying this phenomenon remains limited. In this study, we use three-dimensional direct numerical simulations with a modified multi-marker volume-of-fluid method to investigate liquid–liquid two-phase turbulence in a Taylor–Couette geometry. The dispersed phase has the same density and viscosity as the continuous phase. The Reynolds number $Re\equiv r_i\omega _i d/\nu$ is fixed at 5200, the volume fraction of the dispersed phase is up to $40\,\%$, and the Weber number $We\equiv \rho u^2_\tau d/\sigma$ is approximately 8. It is found that the increase in the system's drag originates from the contribution of interfacial tension. Specifically, droplets experience significant deformation and stretching in the streamwise direction due to shear near the inner cylinder. Consequently, the rear end of the droplets lags behind the fore head. This causes opposing interfacial tension effects on the fore head and rear end of the droplets. For the fore head of the droplets, the effect of interfacial tension appears to act against the flow direction. For the rear end, the effect appears to act in the flow direction. The increase in the system's drag is attributed primarily to the effect of interfacial tension on the fore head of the droplets which leads to the hindering effect of the droplets on the surrounding continuous phase. This hindering effect disrupts the formation of high-speed streaks, favouring the formation of low-speed ones, which are generally associated with higher viscous stress and drag of the system. This study provides new insights into the mechanism of drag enhancement reported in our previous experiments.

Pore-resolved direct numerical simulations have been performed to investigate the turbulent open-channel flow over a rough and permeable sediment bed, represented by a mono-disperse random sphere pack. After a careful validation, eleven cases were simulated to systemically sample a parameter space spanned by a friction Reynolds number $Re_\tau \in [150, 500]$ and a permeability Reynolds number $Re_K \in [0, 2.8]$. By varying the ratio of flow depth to sphere diameter within a range of $h/D \in \{ 3,5,10,\infty \}$, the influence of both Reynolds numbers on the flow field and the turbulence structure could be investigated independently. The simulation results are analysed within a time–space double-averaging framework, whereas flow visualizations provide insight into instantaneous fields. Based on the drag distribution, we propose a consistent interface description, which can be used to define both near-interface and outer-flow coordinates. In these near-interface coordinates, the profiles of the mean velocity and the total shear stress collapse. Furthermore, the proposed interface definition yields outer-layer coordinates, in which the flow and turbulence statistics over a rough and permeable bed reveal similarity to a smooth-wall flow at a similar $Re_\tau$. Within the parameter space, $Re_\tau$ has a strong influence on the wake region of the velocity profile. In contrast, $Re_K$ changes the wall-blocking effect and the shear intensity, which is reflected by the turbulence structure and vortex orientation in the near-interface region. As streamwise velocity streaks disappear and the vortex inclination increases with higher $Re_K$, differences between near-interface and outer-layer turbulence structure are reduced.

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative $\boldsymbol {\varPhi }_{I}(\boldsymbol {x},t)$ and residual $\boldsymbol {\varPhi }_{R}(\boldsymbol {x},t)$ components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in $\boldsymbol {\varPhi }_{I}(\boldsymbol {x},t)$, whereas $\boldsymbol {\varPhi }_{R}(\boldsymbol {x},t)$ offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

We explore the application of the reference map technique, originally developed for Eulerian simulation of solid mechanics, in Lagrangian kinematics of turbulent flows. Unlike traditional methods based on explicit particle tracking, the reference map facilitates the calculation of flow maps and gradients without the need for particles. This is achieved through an Eulerian update of the reference map, which records the take-off positions of fluid particles. This approach is found to be mathematically equivalent to the work of Leung (J. Comput. Phys., vol. 230, issue 9, 2011, pp. 3500–3524), who computed the flow map of simple two-dimensional flows using an Eulerian approach. We discuss important modifications necessary for its first application to complex three-dimensional turbulent flows, including the conservative, low-dissipation update of the flow map and the treatment of periodic boundary conditions. We first demonstrate the accuracy of finite-time Lyapunov exponent (FTLE) calculations based on the reference map against the standard particle-based approach in a two-dimensional Taylor–Green vortex. Then we apply it to turbulent channel flow at $Re_\tau =180$, where Lagrangian coherent structures identified as ridges of the backward-time FTLE are found to bound vortical regions of flow, consistent with Eulerian coherent structures from the $Q$-criterion. The reference map also proves suitable for material surface tracking despite not explicitly tracking particles. This capability can provide valuable insights into the Lagrangian landscape of turbulent momentum transport, complementing Eulerian velocity field analysis. The evolution of initially wall-normal material surfaces in the viscous sublayer, buffer layer and log layer sheds light on the Reynolds stress-generating events from a Lagrangian perspective.

We introduce a novel recursive procedure to a neural-network-based subgrid-scale (NN-based SGS) model for large eddy simulation (LES) of high-Reynolds-number turbulent flow. This process is designed to allow an SGS model to be applicable to a hierarchy of different grid sizes without requiring expensive filtered direct numerical simulation (DNS) data: (1) train an NN-based SGS model with filtered DNS data at a low Reynolds number; (2) apply the trained SGS model to LES at a higher Reynolds number; (3) update this SGS model with training data augmented with filtered LES (fLES) data, accommodating coarser filter size; (4) apply the updated NN to LES at a further higher Reynolds number; (5) go back to Step (3) until a target (very coarse) filter size divided by the Kolmogorov length scale is reached. We also construct an NN-based SGS model using a dual NN architecture whose outputs are the SGS normal stresses for one NN and the SGS shear stresses for the other NN. The input is composed of the velocity gradient tensor and grid size. Furthermore, for the application of an NN-based SGS model trained with one flow to another flow, we modify the NN by eliminating bias and introducing a leaky rectified linear unit function as an activation function. The present recursive SGS model is applied to forced homogeneous isotropic turbulence (FHIT) and successfully predicts FHIT at high Reynolds numbers. The present model trained from FHIT is also applied to decaying homogeneous isotropic turbulence and shows an excellent prediction performance.

We provide scaling relations for the Nusselt number $Nu$ and the friction coefficient $C_{S}$ in sheared Rayleigh–Bénard convection, i.e. in Rayleigh–Bénard flow with Couette- or Poiseuille-type shear forcing, by extending the Grossmann & Lohse (J. Fluid Mech., vol. 407, 2000, pp. 27–56, Phys. Rev. Lett., vol. 86, 2001, pp. 3316–3319, Phys. Rev. E, vol. 66, 2002, 016305, Phys. Fluids, vol. 16, 2004, pp. 4462–4472) theory to sheared thermal convection. The control parameters for these systems are the Rayleigh number $Ra$, the Prandtl number $Pr$ and the Reynolds number $Re_S$ that characterises the strength of the imposed shear. By direct numerical simulations and theoretical considerations, we show that, in turbulent Rayleigh–Bénard convection, the friction coefficients associated with the applied shear and the shear generated by the large-scale convection rolls are both well described by Prandtl's (Ergeb. Aerodyn. Vers. Gött., vol. 4, 1932, pp. 18–29) logarithmic friction law, suggesting some kind of universality between purely shear-driven flows and thermal convection. These scaling relations hold well for $10^6 \leq Ra \leq 10^8$, $0.5 \leq Pr \leq 5.0$, and $0 \leq Re_S \leq 10^4$.

We employ direct numerical simulations to investigate the heat transfer and flow structures in turbulent Rayleigh–Bénard convection in both cylindrical cells and laterally periodic domains, spanning an unprecedentedly wide range of aspect ratios $0.075 \leqslant \varGamma \leqslant 32$. We focus on Prandtl number ${Pr}=1$ and Rayleigh numbers ${{Ra}}=2\times 10^7$ and ${{Ra}}=10^8$. In both cases, with increasing aspect ratio, the heat transfer first increases, then reaches a maximum (which is more pronounced for the cylindrical case due to confinement effects), and then slightly goes down again before it finally saturates at the large aspect ratio limit, which is achieved already at $\varGamma \approx 4$. Already for $\varGamma \gtrsim 0.75$, the heat transfers in both cylindrical and laterally periodic domains become identical. The large-$\varGamma$ limit for the volume-integrated Reynolds number and the boundary layer thicknesses are also reached at $\varGamma \approx 4$. However, while the integral flow properties converge at $\varGamma \approx 4$, the confinement of a cylindrical domain impacts the temperature and velocity variance distributions up to $\varGamma \approx 16$, as thermal superstructures cannot form close to the sidewall.

This DNS study considers transition to turbulence in plane Couette flow (pCf) with a rough stationary wall and a smooth moving wall. The roughness elements are square ribs of height $k=0.2h$ (where $h$ is the half-channel height). Two different pitch separations, $\lambda =2k$ and $10k$, are considered, i.e. d-type and k-type roughness, respectively. The transition in both rough pCf cases takes place through a stage of alternate laminar–turbulent bands aligned in an oblique fashion. However, roughness causes a shift in the transitional Reynolds number ($Re$) range. In the k-type roughness, stable bands are observed in the range $Re \in [300, 325]$, which is a downwards shift from the transitional $Re$ range for the smooth pCf ($Re \in [325,400]$). The d-type roughness, on the other hand, surprisingly shifts the transitional $Re$ range upwards to $Re \in [350,425]$. This peculiar behaviour is associated with the ability of the ribs to shed and regenerate vorticity. Large-scale components extracted using a filtering process relate to the transition bands and flow parallel to the oblique laminar–turbulent boundaries. Counter-rotating vortices are present in the turbulent regions of the flow field, which exist in tandem with the high- and low-velocity streaks. Another interesting observation is the secondary Reynolds shear stresses, $-\overline {v^{\prime }u^{\prime }}$ and $-\overline {w^{\prime }v^{\prime }}$, which are non-zero in the transitional regime, in contrast to the turbulent regime where they are negligible.

We investigate the coupling effects of the two-phase interface, viscosity ratio and density ratio of the dispersed phase to the continuous phase on the flow statistics in two-phase Taylor–Couette turbulence at a system Reynolds number of $6\times 10^3$ and a system Weber number of 10 using interface-resolved three-dimensional direct numerical simulations with the volume-of-fluid method. Our study focuses on four different scenarios: neutral droplets, low-viscosity droplets, light droplets and low-viscosity light droplets. We find that neutral droplets and low-viscosity droplets primarily contribute to drag enhancement through the two-phase interface, whereas light droplets reduce the system's drag by explicitly reducing Reynolds stress due to the density dependence of Reynolds stress. In addition, low-viscosity light droplets contribute to greater drag reduction by further reducing momentum transport near the inner cylinder and implicitly reducing Reynolds stress. While interfacial tension enhances turbulent kinetic energy (TKE) transport, drag enhancement is not strongly correlated with TKE transport for both neutral droplets and low-viscosity droplets. Light droplets primarily reduce the production term by diminishing Reynolds stress, whereas the density contrast between the phases boosts TKE transport near the inner wall. Therefore, the reduction in the dissipation rate is predominantly attributed to decreased turbulence production, causing drag reduction. For low-viscosity light droplets, the production term diminishes further, primarily due to their greater reduction in Reynolds stress, while reduced viscosity weakens the density difference's contribution to TKE transport near the inner cylinder, resulting in a more pronounced reduction in the dissipation rate and consequently stronger drag reduction. Our findings provide new insights into the physics of turbulence modulation by the dispersed phase in two-phase turbulence systems.

We study the effect of surface texture on an overlying turbulent flow for the case of textures made of an alternating slip/no-slip pattern, a common model for superhydrophobic surfaces, but also a particularly simple form of texture. For texture sizes $L^+ \gtrsim 25$, we have previously reported that, even though the texture effectively imposes homogeneous slip boundary conditions on the overlying, background turbulence, this is not its sole effect. The effective conditions only produce an origin offset on the background turbulence, which remains otherwise smooth-wall-like. For actual textures, however, as their size increases from $L^+ \gtrsim 25$ the flow progressively departs from this smooth-wall-like regime, resulting in additional shear Reynolds stress and increased drag, in a non-homogeneous fashion that could not be reproduced by the effective boundary conditions. This paper focuses on the underlying physical mechanism of this phenomenon. We argue that it is caused by the nonlinear interaction of the texture-coherent flow, directly induced by the surface topology, and the background turbulence, as it acts directly on the latter and alters it. This does not occur at the boundary where effective conditions are imposed, but within the overlying flow itself, where the interaction acts as a forcing on the governing equations of the background turbulence, and takes the form of cross-advective terms between the latter and the texture-coherent flow. We show this by conducting simulations where we remove the texture and introduce additional, forcing terms in the Navier–Stokes equations, in addition to the equivalent homogeneous slip boundary conditions. The forcing terms capture the effect of the nonlinear interaction on the background turbulence without the need to resolve the surface texture. We show that, when the forcing terms are derived accounting for the amplitude modulation of the texture-coherent flow by the background turbulence, they quantitatively capture the changes in the flow up to texture sizes $L^+ \approx 70{-}100$. This includes not just the roughness function but also the changes in the flow statistics and structure.

Turbulent circular pipe flows subjected to axial system rotation are studied using direct numerical simulations (DNS) for a wide range of rotation numbers of $Ro_b = 0\unicode{x2013}20$ at a fixed Reynolds number. To ensure that energetic turbulent eddy motions are captured at high rotation numbers, long pipes up to $L_z = 180{\rm \pi} R$ are used in DNS. Two types of energy-containing flow structures have been observed. The first type is hairpin structures that are characteristic of the turbulent boundary layer developing over the pipe wall for both non-rotating and axially rotating flows. The second type is Taylor columns forming at moderate and high rotation numbers. Based on the study of two-point autocorrelation coefficients, it is observed that Taylor columns exhibit quasi-periods in both axial and azimuthal directions. According to the premultiplied spectra, Taylor columns feature one single characteristic axial length scale at the moderate rotation numbers but two at high rotation numbers. It is discovered that the axial system rotation suppresses the sweep events systematically and impedes the formation of hairpin structures. As the rotation number is increased, the turbulence kinetic energy held by Taylor columns enhances rapidly associated with significant increases in their axial length scales.

The friction drag of the axial flow along the outer surface of a cylinder varies with the cylinder radius and flow conditions. This study included direct numerical simulations of the axial turbulent flow along a circular cylinder under different conditions for obtaining the turbulence statistics and wall friction coefficient. Then the characteristics of velocity streaks were observed from a geometrical perspective of turbulence structures around the circular cylinder, and compared with the characteristics of the turbulence structures in a boundary layer on a flat plate. The results showed that the velocity streak spacing and the distance between the velocity streak and the cylinder surface in the viscous length scale do not vary substantively with the radius of the cylinder, and are the same as those of the turbulent flow along a flat plate. Therefore, they can be considered geometrical characteristics of the turbulence structure independent of the cylinder radius. Moreover, the friction coefficient per pair of high- and low-speed velocity streaks is the same as that of flat-plate turbulent flow, independent of the cylinder radius, and can be regarded as a dynamical characteristic for a pair of velocity streaks. Two equations were derived based on the characteristics of wall turbulence. The characteristics of the turbulence predicted by the two formulae were consistent with the simulation results. Consequently, we showed that the wall friction coefficient and number of the velocity streak pairs, which are statistical and structural characteristics of wall turbulence, can be predicted appropriately by specifying the radius Reynolds number.

Simulations of elastic turbulence, the chaotic flow of highly elastic and inertialess polymer solutions, are plagued by numerical difficulties: the chaotically advected polymer conformation tensor develops extremely large gradients and can lose its positive-definiteness, which triggers numerical instabilities. While efforts to tackle these issues have produced a plethora of specialized techniques – tensor decompositions, artificial diffusion, and shock-capturing advection schemes – we still lack an unambiguous route to accurate and efficient simulations. In this work, we show that even when a simulation is numerically stable, maintaining positive-definiteness and displaying the expected chaotic fluctuations, it can still suffer from errors significant enough to distort the large-scale dynamics and flow structures. We focus on two-dimensional simulations of the Oldroyd-B and FENE-P equations, driven by a large-scale cellular body forcing. We first compare two positivity-preserving decompositions of the conformation tensor: symmetric square root (SSR) and Cholesky with a logarithmic transformation (Cholesky-log). While both simulations yield chaotic flows, only the latter preserves the pattern of the forcing, i.e. its fluctuating vortical cells remain ordered in a lattice. In contrast, the SSR simulation exhibits distorted vortical cells that shrink, expand and reorient constantly. To identify the accurate simulation, we appeal to a hitherto overlooked mathematical bound on the determinant of the conformation tensor, which unequivocally rejects the SSR simulation. Importantly, the accuracy of the Cholesky-log simulation is shown to arise from the logarithmic transformation. We also consider local artificial diffusion, a potential low-cost alternative to high-order advection schemes. Unfortunately, the artificially enhanced diffusive smearing of polymer stress in regions of intense stretching substantially modifies the global dynamics. We then show how the spurious large-scale motions, identified here, contaminate predictions of scalar mixing. Finally, we discuss the effects of spatial resolution, which controls the steepness of gradients in a non-diffusive simulation.

Non-helical turbulence within a linear shear flow has demonstrated efficient amplification of large-scale magnetic fields in numerical simulations, but its precise mechanism remains elusive. The incoherent $\alpha$ mechanism proposes that a zero-mean fluctuating transport coefficient $\alpha$ (linked to kinetic helicity) in the shear flow is a candidate driver. Previous renovating-flow models have proposed that the correlation time of helicity fluctuations must be sufficiently extended to overcome turbulent magnetic diffusivity, yet only empirical validation of this concept has been obtained. In this study, we conduct direct numerical simulations of weakly compressible non-helical hydrodynamic turbulence. We scrutinize the correlation times of velocity and kinetic helicity fluctuations in distinct flow configurations, including rotation, shearing and Keplerian flows, as well as the shearing burgulence counterpart. Our findings indicate that rotation contributes to a prolonged correlation time of helicity compared with velocity, particularly notable in auto-correlations of both volume-averaged quantities and individual Fourier modes due to the formation of large-scale vortices. In contrast, moderate shear strength does not exhibit significant scale separation, with shear flows elongating vortices in the shear direction. Shearing burgulence, characterized by shorter helicity correlation times, appears less conducive to hosting the incoherent $\alpha$ effect. Notably, at modest shear rates, only Keplerian flows exhibit sufficiently coherent helicity fluctuations, in contrast to shearing flows. However, the relative strength of helicity fluctuations compared with turbulent diffusivity is significantly lower, raising doubts about the viability of the incoherent $\alpha$ effect as a potential dynamo driver in the subsonic flows examined in this study.

The dimensional transition in turbulent jets of a shear-thinning fluid is studied via direct numerical simulations. Our findings reveal that under vertical confinement, the flow exhibits a unique mixed-dimensional (or 2.5-dimensional) state, where large-scale two-dimensional and small-scale three-dimensional structures coexist. This transition from three-dimensional turbulence near the inlet to two-dimensional dynamics downstream is dictated by the level of confinement: weak confinement guarantees turbulence to remain three-dimensional, whereas strong confinement forces the transition to two dimensions; the mixed-dimensional state is observed for moderate confinement and it emerges as soon as flow scales are larger than the vertical length. In this scenario, we observed that the mixed-dimensional state is an overall more energetic state, and it shows a multi-cascade process, where the direct cascade of energy at small scales and the direct cascade of enstrophy at large scales coexist. The results provide insights into the complex dynamics of confined turbulent flows, relevant in both natural and industrial settings.

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan (J. Fluid Mech., vol. 974, 2023, A35) and Jacobs & Dalziel (J. Fluid Mech., vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

The research on elasto-inertial turbulence (EIT), a new type of turbulent flow, has reached the stage of identifying the minimal flow unit (MFU). On this issue, direct numerical simulations of FENE-P fluid flow in two-dimensional channels with variable sizes are conducted in this study. We demonstrate with the increase of channel length that the simulated flow experiences several different flow patterns, and there exists an MFU for EIT to be self-sustained. At Weissenberg number ($Wi$) higher than the one required to excite EIT, when the channel length is relatively small, a steady arrowhead regime (SAR) flow structure and a laminar-like friction coefficient is achieved. However, as the channel length increases, the flow can fully develop into EIT characterized with high flow drag. Close to the size of the MFU, the simulated flow behaves intermittently between the SAR state with low drag and EIT state with high drag. The flow falling back to ‘laminar flow’ is caused by the insufficient channel size below the MFU. Furthermore, we give the relationship between the value of the MFU and the effective $Wi$, and explain its physical reasons. Moreover, the intermittent flow regime obtained based on the MFU gives us an opportunity to look into the origin and exciting process of EIT. Through capturing the onset process of EIT, we observed that EIT originates from the sheet-like extension structure located near the wall, which is maybe related to the wall mode rather than the centre mode. The fracture and regeneration of this sheet-like structure is the key mechanism for the self-sustaining of EIT.

We report direct numerical simulations results of the rough-wall channel, focusing on roughness with high $k_{rms}/k_a$ statistics but small to negative $Sk$ statistics, and we study the implications of this new dataset on rough-wall modelling. Here, $k_{rms}$ is the root mean square, $k_a$ is the first-order moment of roughness height, and $Sk$ is the skewness. The effects of packing density, skewness and arrangement of roughness elements on mean streamwise velocity, equivalent roughness height ($z_0$) and Reynolds and dispersive stresses have been studied. We demonstrate that two-point correlation lengths of roughness height statistics play an important role in characterizing rough surfaces with identical moments of roughness height but different arrangements of roughness elements. Analysis of the present as well as historical data suggests that the task of rough-wall modelling is to identify geometric parameters that distinguish the rough surfaces within the calibration dataset. We demonstrate a novel feature selection procedure to determine these parameters. Further, since there is no finite set of roughness statistics that distinguish between all rough surfaces, we argue that obtaining a universal rough-wall model for making equivalent sand-grain roughness ($k_s$) predictions would be challenging, and that each rough-wall model would have its applicable range. This motivates the development of group-based rough-wall models. The applicability of multi-variate polynomial regression and feedforward neural networks for building such group-based rough-wall models using the selected features has been shown.

The present study offers a twofold contribution on counter-gradient transport (CGT) of turbulent scalar flux. First, by examining turbulent scalar mixing through synchronized particle image velocimetry and planar laser-induced fluorescence on an inclined jet in cross-flow, we clarify the previously unexplained phenomenon of CGT, revealing key flow structures, their spatial distribution and modelling implications. Statistical analysis identifies two distinct CGT regions: local cross-gradient transport in the windward shear layer and non-local effects near the wall after injection. These behaviours are driven by specific flow structures, namely Kelvin–Helmholtz vortices (local) and wake vortices (non-local), suggesting that scalar flux can be decomposed into a gradient-type term for gradient diffusion and a term for large-eddy stirring. Second, we propose a new approach for reconstruction of turbulent mean flow and scalar fields using continuous adjoint data assimilation (DA). By rectifying model-form errors through anisotropic correction under observational constraints, our DA model minimizes discrepancies between experimental measurements and numerical predictions. As expected, the introduced forcing term effectively identifies regions where traditional models fall short, particularly in the jet centreline and near-wall regions, thereby enhancing the accuracy of the mean scalar field. These enhancements occur not only within the observation region but also in unseen regions, underscoring present DA approach's reliability and practicality for reproducing mean flow behaviours from limited data. These findings lay a solid foundation for adjoint-based model-consistent data-driven methods, offering promising potential for accurately predicting complex flow scenarios like film cooling.