Tamer A. Zaki
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Tamer A. Zaki has a total of 119 co-authors. Here are the top co-authors by number of shared publications:
- Charles Meneveau
Johns Hopkins University
9 shared publications Mengze Wang
8 shared publications- Gregory L. Eyink
Johns Hopkins University
6 shared publications - Qi Wang
San Diego State University
4 shared publications - Amir Esteghamatian
Worldwide Research and Development, Pfizer Inc, Groton, CT 06340, USA
3 shared publications Charles Meneveau
Department of Mechanical Engineering, Johns Hopkins University, Baltimore, USA
3 shared publications
Latest Publications
- 14 October 2026
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 11 (10)
Latent-space variational data assimilation in two-dimensional turbulence
- 14 October 2026
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 11 (10)
Latent-space variational data assimilation in two-dimensional turbulence
- Open Access
- 23 July 2026
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1039, A48
Mitigating adjoint chaos in wall turbulence
Estimating past events in wall turbulence based solely on surface measurements and first principles is an ill-posed problem that is complicated by chaos. The sensitivity of a measurement to the earlier flow state is described by the adjoint Navier–Stokes equations, which are solved in reverse time starting from the measurement kernel at the sensing position and time. The resulting adjoint field is the spatio-temporal domain of dependence (DOD) of the sensor, which is a dual to the concept of the domain of influence (DOI) of an actuator in the linearised forward equations. In channel turbulence, the energy of each adjoint realisation grows exponentially in backward time according to the Lyapunov exponent, even though the energy of the ensemble average should decay. We introduce a linear eddy-viscosity closure model in the ensemble-averaged adjoint equations, and directly compute the mean DOD and compare our prediction with the ensemble average. Furthermore, we demonstrate that the DOD of a wall-stress measurement and the DOI resulting from a wall-stress perturbation exhibit respective universal behaviours across Reynolds numbers. However, their spatio-temporal structures differ qualitatively, due to the time asymmetry of the governing equations. The DOD field has a two-part structure: one component is associated with the Orr mechanism, characterised by rapid reorientation under mean shear, and the other is related to self-similar expanding streaky structures. These two components jointly define the sensitivity of the wall-stress measurement to past flow events.
- 23 July 2026
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1039, A48
Open AccessMitigating adjoint chaos in wall turbulence
Estimating past events in wall turbulence based solely on surface measurements and first principles is an ill-posed problem that is complicated by chaos. The sensitivity of a measurement to the earlier flow state is described by the adjoint Navier–Stokes equations, which are solved in reverse time starting from the measurement kernel at the sensing position and time. The resulting adjoint field is the spatio-temporal domain of dependence (DOD) of the sensor, which is a dual to the concept of the domain of influence (DOI) of an actuator in the linearised forward equations. In channel turbulence, the energy of each adjoint realisation grows exponentially in backward time according to the Lyapunov exponent, even though the energy of the ensemble average should decay. We introduce a linear eddy-viscosity closure model in the ensemble-averaged adjoint equations, and directly compute the mean DOD and compare our prediction with the ensemble average. Furthermore, we demonstrate that the DOD of a wall-stress measurement and the DOI resulting from a wall-stress perturbation exhibit respective universal behaviours across Reynolds numbers. However, their spatio-temporal structures differ qualitatively, due to the time asymmetry of the governing equations. The DOD field has a two-part structure: one component is associated with the Orr mechanism, characterised by rapid reorientation under mean shear, and the other is related to self-similar expanding streaky structures. These two components jointly define the sensitivity of the wall-stress measurement to past flow events.
- Open Access
- 7 July 2026
- Published by Elsevier in Proceedings of the Combustion Institute
- Vol. 42
Parasitic combustion model-based data assimilation for predicting wave speeds in rotating detonation engines
- 7 July 2026
- Published by Elsevier in Proceedings of the Combustion Institute
- Vol. 42
Open AccessParasitic combustion model-based data assimilation for predicting wave speeds in rotating detonation engines
- Open Access
- 19 June 2026
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1037, A28
Assimilation of wall-pressure measurements in direct numerical simulations of high-speed flow over a cone–flare geometry
Ensemble-variational (EnVar) assimilation of wall-pressure measurements in direct numerical simulations of Mach 6 flow over a cone–flare is performed. The experimental data include pressure spectra and intensities from seven wall-mounted PCB sensors positioned upstream, within and downstream of the separation region induced by the compression corner. Assimilation of the first two sensors only, all upstream of separation, is insufficient to accurately predict the downstream flow. Assimilating all the sensor data is shown to be essential to correctly predict separation onset and the downstream wall-pressure data. Similar to the experiments, the assimilated flow features intense rope-like structures in the attached region. The simulations additionally predict a localised amplification of disturbances beneath the separation shock, where experimental data are not available. This amplification results from the interaction of the boundary-layer instability modes with the compression shock. The simulations also capture the sharp decrease in wall-pressure intensity across separation, and the amplification of low-frequency three-dimensional disturbances within the recirculation bubble. Additionally, the computations highlight the uncertainty in the post-separation predictions due to the low-frequency unsteadiness of the separation shock. Oscillations of the streamwise velocity modulate the boundary-layer thickness, which in turn introduces variability in disturbance amplification.
- 19 June 2026
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1037, A28
Open AccessAssimilation of wall-pressure measurements in direct numerical simulations of high-speed flow over a cone–flare geometry
Ensemble-variational (EnVar) assimilation of wall-pressure measurements in direct numerical simulations of Mach 6 flow over a cone–flare is performed. The experimental data include pressure spectra and intensities from seven wall-mounted PCB sensors positioned upstream, within and downstream of the separation region induced by the compression corner. Assimilation of the first two sensors only, all upstream of separation, is insufficient to accurately predict the downstream flow. Assimilating all the sensor data is shown to be essential to correctly predict separation onset and the downstream wall-pressure data. Similar to the experiments, the assimilated flow features intense rope-like structures in the attached region. The simulations additionally predict a localised amplification of disturbances beneath the separation shock, where experimental data are not available. This amplification results from the interaction of the boundary-layer instability modes with the compression shock. The simulations also capture the sharp decrease in wall-pressure intensity across separation, and the amplification of low-frequency three-dimensional disturbances within the recirculation bubble. Additionally, the computations highlight the uncertainty in the post-separation predictions due to the low-frequency unsteadiness of the separation shock. Oscillations of the streamwise velocity modulate the boundary-layer thickness, which in turn introduces variability in disturbance amplification.
- Open Access
- 9 December 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1024
Back-in-time analysis of vorticity in viscous separated flows over immersed bodies
We present a back-in-time analysis for the origin of vorticity in viscous separated flows over immersed bodies, using the adjoint-vorticity framework recently introduced by Xiang et al. (2025 J. Fluid Mech. vol. 1011, A33. The solution of the adjoint-vorticity equations yields the volume density of mean deformation, which captures the stretching and tilting of the earlier vorticity that leads to the terminal value. The analysis also takes into account the boundary contributions of vorticity and its flux. Three examples are considered. Steady, axisymmetric separation in the flow over a sphere at Reynolds number $Re=200$ is shown to be established due to wall flux from both upstream and downstream of separation, the latter contribution being absent from the classical description by Lighthill. For unsteady separation at higher $Re=300$ , the streamwise vorticity within the wake hairpin vortex is traced back, quantitatively, to the azimuthal vorticity on the sphere. The third configuration is the flow over a prolate spheroid at $Re=3000$ . The null vorticity at three-dimensional separation originates from the cancellation of opposite interior contributions adjacent to the separation surface. The contribution from the downstream side migrates across the separation surface into the upstream region due to a tilting effect – a fundamental distinction between two- and three-dimensional separation. We also examine the detached vortical structures. The streamwise vorticity in the primary vortex originates from tilting of near-wall azimuthal vorticity, differing from Lighthill’s conjecture that the origin is streamwise near-wall vorticity that arises due to the reduced Coriolis force. Finally, a necklace vortex in the turbulent wake is traced back in time, and is shown to have contributions from the spheroid trailing-edge shed shear layer and the large-scale counter-rotating primary vortices.- 9 December 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1024
Open AccessBack-in-time analysis of vorticity in viscous separated flows over immersed bodies
We present a back-in-time analysis for the origin of vorticity in viscous separated flows over immersed bodies, using the adjoint-vorticity framework recently introduced by Xiang et al. (2025 J. Fluid Mech. vol. 1011, A33. The solution of the adjoint-vorticity equations yields the volume density of mean deformation, which captures the stretching and tilting of the earlier vorticity that leads to the terminal value. The analysis also takes into account the boundary contributions of vorticity and its flux. Three examples are considered. Steady, axisymmetric separation in the flow over a sphere at Reynolds number $Re=200$ is shown to be established due to wall flux from both upstream and downstream of separation, the latter contribution being absent from the classical description by Lighthill. For unsteady separation at higher $Re=300$ , the streamwise vorticity within the wake hairpin vortex is traced back, quantitatively, to the azimuthal vorticity on the sphere. The third configuration is the flow over a prolate spheroid at $Re=3000$ . The null vorticity at three-dimensional separation originates from the cancellation of opposite interior contributions adjacent to the separation surface. The contribution from the downstream side migrates across the separation surface into the upstream region due to a tilting effect – a fundamental distinction between two- and three-dimensional separation. We also examine the detached vortical structures. The streamwise vorticity in the primary vortex originates from tilting of near-wall azimuthal vorticity, differing from Lighthill’s conjecture that the origin is streamwise near-wall vorticity that arises due to the reduced Coriolis force. Finally, a necklace vortex in the turbulent wake is traced back in time, and is shown to have contributions from the spheroid trailing-edge shed shear layer and the large-scale counter-rotating primary vortices. - Open Access
- 3 November 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1022, A30
Lighthill’s mechanism and vorticity cascade in the logarithmic layer of wall turbulence
We investigate Lighthill’s proposed turbulent mechanism for near-wall concentration of spanwise vorticity by calculating mean flows conditioned on motion away from or toward the wall in an (friction Reynolds number) ${\textit{Re}}_\tau =1000$ database of plane-parallel channel flow. Our results corroborate Lighthill’s proposal throughout the entire logarithmic layer, but extended by counter-flows that help explain anti-correlation of vorticity transport by advection and by stretching/tilting. We present evidence also for Lighthill’s hypothesis that the vorticity transport in the log layer is a ‘cascade process’ through a scale hierarchy of eddies, with intense competition between transport outward from and inward to the wall. Townsend’s model of attached eddies of hairpin-vortex type accounts for half of the vorticity cascade, whereas we identify necklace type or ’shawl vortices’ that envelop turbulent sweeps as supplying the other half.- 3 November 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1022, A30
Open AccessLighthill’s mechanism and vorticity cascade in the logarithmic layer of wall turbulence
We investigate Lighthill’s proposed turbulent mechanism for near-wall concentration of spanwise vorticity by calculating mean flows conditioned on motion away from or toward the wall in an (friction Reynolds number) ${\textit{Re}}_\tau =1000$ database of plane-parallel channel flow. Our results corroborate Lighthill’s proposal throughout the entire logarithmic layer, but extended by counter-flows that help explain anti-correlation of vorticity transport by advection and by stretching/tilting. We present evidence also for Lighthill’s hypothesis that the vorticity transport in the log layer is a ‘cascade process’ through a scale hierarchy of eddies, with intense competition between transport outward from and inward to the wall. Townsend’s model of attached eddies of hairpin-vortex type accounts for half of the vorticity cascade, whereas we identify necklace type or ’shawl vortices’ that envelop turbulent sweeps as supplying the other half. - Open Access
- 23 September 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1019, A54
Analysis of flow-wall deformation coupling in high Reynolds number compliant wall boundary layers
Interactions of turbulent boundary layers with a compliant surface are investigated experimentally at Re τ = 3300–8900. Integrating tomographic particle tracking with Mach–Zehnder interferometry enables simultaneous mapping of the compliant wall deformation and the three-dimensional velocity and pressure fields. Our initial study (J. Fluid. Mech. vol. 980, R2) shows that the flow–deformation correlations decrease with increasing Re τ , despite an order of magnitude increase in deformation amplitude. To elucidate the mechanisms involved, the same velocity, pressure and kinetic energy fields are decomposed to ‘wave-coherent’ and ‘stochastic’ parts using a Hilbert projection method. The phase dependent coherent variables, especially the pressure, are highly correlated with the wave, but decrease with increasing Re τ . While the coherent energy is 6 %–10 % of the stochastic level, the pressure root mean square is comparable near the wall. The energy flux between the coherent and stochastic parts and the pressure diffusion reverse sign at the critical layer. To explain the Re τ dependence, the characteristic deformation wavelength (three times the thickness) is compared with the scales of the energy-containing eddies in the boundary layer represented by the k−1 range in the energy spectrum. When the deformation wavelength is matched with the k x E uu peak at the present lowest Re τ , the flow–deformation correlations and coherent pressure become strong, even for submicron deformations. In this case, the flow and wall motion become phase locked, suggesting resonant behaviours. As Re τ increases, the wall wavelengths and spectral range of attached eddies are no longer matched, resulting in reduced correlations and lower coherent energy and pressure, despite larger deformation.
- 23 September 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1019, A54
Open AccessAnalysis of flow-wall deformation coupling in high Reynolds number compliant wall boundary layers
Interactions of turbulent boundary layers with a compliant surface are investigated experimentally at Re τ = 3300–8900. Integrating tomographic particle tracking with Mach–Zehnder interferometry enables simultaneous mapping of the compliant wall deformation and the three-dimensional velocity and pressure fields. Our initial study (J. Fluid. Mech. vol. 980, R2) shows that the flow–deformation correlations decrease with increasing Re τ , despite an order of magnitude increase in deformation amplitude. To elucidate the mechanisms involved, the same velocity, pressure and kinetic energy fields are decomposed to ‘wave-coherent’ and ‘stochastic’ parts using a Hilbert projection method. The phase dependent coherent variables, especially the pressure, are highly correlated with the wave, but decrease with increasing Re τ . While the coherent energy is 6 %–10 % of the stochastic level, the pressure root mean square is comparable near the wall. The energy flux between the coherent and stochastic parts and the pressure diffusion reverse sign at the critical layer. To explain the Re τ dependence, the characteristic deformation wavelength (three times the thickness) is compared with the scales of the energy-containing eddies in the boundary layer represented by the k−1 range in the energy spectrum. When the deformation wavelength is matched with the k x E uu peak at the present lowest Re τ , the flow–deformation correlations and coherent pressure become strong, even for submicron deformations. In this case, the flow and wall motion become phase locked, suggesting resonant behaviours. As Re τ increases, the wall wavelengths and spectral range of attached eddies are no longer matched, resulting in reduced correlations and lower coherent energy and pressure, despite larger deformation.
- 22 September 2025
- Published by Elsevier in Computer Methods in Applied Mechanics and Engineering
- Vol. 447
Accelerating Hamiltonian Monte Carlo for Bayesian inference in neural networks and neural operators
- 22 September 2025
- Published by Elsevier in Computer Methods in Applied Mechanics and Engineering
- Vol. 447
Accelerating Hamiltonian Monte Carlo for Bayesian inference in neural networks and neural operators
- 17 September 2025
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 10 (9)
Assimilation of wall-pressure measurements in high-speed boundary layers using a Bayesian optimization with DeepONet
- 17 September 2025
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 10 (9)
Assimilation of wall-pressure measurements in high-speed boundary layers using a Bayesian optimization with DeepONet