Showing papers in "Journal of Fluid Mechanics in 2022"

Direct numerical simulation of hypersonic turbulent boundary layers: effect of spatial evolution and Reynolds number

Abstract: Abstract Direct numerical simulations (DNS) are performed to investigate the spatial evolution of flat-plate zero-pressure-gradient turbulent boundary layers over long streamwise domains (${>}300\delta _i$, with $\delta _i$ the inflow boundary-layer thickness) at three different Mach numbers, $2.5$, $4.9$ and $10.9$, with the surface temperatures ranging from quasiadiabatic to highly cooled conditions. The settlement of turbulence statistics into a fully developed equilibrium state of the turbulent boundary layer has been carefully monitored, either based on the satisfaction of the von Kármán integral equation or by comparing runs with different inflow turbulence generation techniques. The generated DNS database is used to characterize the streamwise evolution of multiple important variables in the high-Mach-number, cold-wall regime, including the skin friction, the Reynolds analogy factor, the shape factor, the Reynolds stresses, and the fluctuating wall quantities. The data confirm the validity of many classic and newer compressibility transformations at moderately high Reynolds numbers (up to friction Reynolds number $Re_\tau \approx 1200$) and show that, with proper scaling, the sizes of the near-wall streaks and superstructures are insensitive to the Mach number and wall cooling conditions. The strong wall cooling in the hypersonic cold-wall case is found to cause a significant increase in the size of the near-wall turbulence eddies (relative to the boundary-layer thickness), which leads to a reduced-scale separation between the large and small turbulence scales, and in turn to a lack of an outer peak in the spanwise spectra of the streamwise velocity in the logarithmic region.

Interaction of cavitation bubbles with the interface of two immiscible fluids on multiple time scales

Abstract: We experimentally, numerically and theoretically investigate the nonlinear interaction between a cavitation bubble and the interface of two immiscible fluids (oil and water) on multiple time scales. The underwater electric discharge method is utilized to generate a cavitation bubble near or at the interface. Both the bubble dynamics on a short time scale and the interface evolution on a much longer time scale are recorded via high-speed photography. Two mechanisms are found to contribute to the fluid mixing in our system. First, when a bubble is initiated in the oil phase or at the interface, an inertia-dominated high-speed liquid jet generated from the collapsing bubble penetrates the water–oil interface, and consequently transports fine oil droplets into the water. The critical standoff parameter for jet penetration is found to be highly dependent on the density ratio of the two fluids. Furthermore, the pinch-off of an interface jet produced long after the bubble dynamics stage is reckoned as the second mechanism, carrying water droplets into the oil bulk. The dependence of the bubble jetting behaviours and interface jet dynamics on the governing parameters is systematically studied via experiments and boundary integral simulations. Particularly, we quantitatively demonstrate the respective roles of surface tension and viscosity in interface jet dynamics. As for a bubble initiated at the interface, an extended Rayleigh–Plesset model is proposed that well predicts the asymmetric dynamics of the bubble, which accounts for a faster contraction of the bubble top and a downward liquid jet.

Closure mechanism of the A1 and A2 modes in jet screech

Abstract: This paper explores the screech closure mechanism for different axisymmetric modes in shock-containing jets. While many of the discontinuities in tonal frequency exhibited by screeching jets can be associated with a change in the azimuthal mode, there has to date been no explanation for the existence of multiple axisymmetric modes at different frequencies. This paper provides just such an explanation. As shown in previous works, specific wavenumbers arise from the interaction of waves in the flow with the shocks. This provides new paths for driving upstream-travelling waves that can potentially close the resonance loop. Predictions using locally parallel and spatially periodic linear stability analyses and the wavenumber spectrum of the shock-cell structure suggest that the A1 mode resonance is closed by a wave generated when the Kelvin-Helmholtz mode interacts with the leading wavenumber of the shock-cell structure. The A2 mode is closed by a wave that arises due to interaction between the Kelvin-Helmholtz wave and a secondary wavenumber peak, which arises from the spatial variation of the shock-cell wavelength. The predictions are shown to closely match experimental data, and possible justifications for the dominance of each mode are provided based on the growth rates of the absolute instability.

High-resolution direct simulation of deep water breaking waves: transition to turbulence, bubbles and droplets production

Abstract: We present high-resolution three-dimensional (3-D) direct numerical simulations of breaking waves solving for the two-phase Navier–Stokes equations. We investigate the role of the Reynolds number ( Re , wave inertia relative to viscous effects) and Bond number ( Bo , wave scale over the capillary length) on the energy, bubble and droplet statistics of strong plunging breakers. We explore the asymptotic regimes at high Re and Bo , and compare with laboratory breaking waves. Energetically, the breaking wave transitions from laminar to 3-D turbulent flow on a time scale that depends on the turbulent Re up to a limiting value $Re_\lambda \sim 100$ , consistent with the mixing transition in other canonical turbulent flows. We characterize the role of capillary effects on the impacting jet and ingested main cavity shape and subsequent fragmentation process, and extend the buoyant-energetic scaling from Deike et al. ( J. Fluid Mech. , vol. 801, 2016, pp. 91–129) to account for the cavity shape and its scale separation from the Hinze scale, $r_H$ . We confirm two regimes in the bubble size distribution, $N(r/r_H)\propto (r/r_H)^{-10/3}$ for $r>r_H$ , and $\propto (r/r_H)^{-3/2}$ for $r<r_H$ . Bubbles are resolved up to one order of magnitude below $r_H$ , and we observe a good collapse of the numerical data compared to laboratory breaking waves (Deane & Stokes, Nature , vol. 418 (6900), 2002, pp. 839–844). We resolve droplet statistics at high Bo in good agreement with recent experiments (Erinin et al. , Geophys. Res. Lett. , vol. 46 (14), 2019, pp. 8244–8251), with a distribution shape close to $N_d(r_d)\propto r_d^{-2}$ . The evolution of the droplet statistics appears controlled by the details of the impact process and subsequent splash-up. We discuss velocity distributions for the droplets, finding ejection velocities up to four times the phase speed of the wave, which are produced during the most intense splashing events of the breaking process.

Comprehensive analysis of spherical bubble oscillations and shock wave emission in laser-induced cavitation

Abstract: The dynamics of spherical laser-induced cavitation bubbles in water is investigated by plasma photography, time-resolved shadowgraphs, and single-shot probe beam scattering enabling to portray the transition from initial nonlinear to late linear oscillations. The frequency of late oscillations yields the bubble's gas content. Simulations with the Gilmore model using plasma size as input and oscillation times as fit parameter provide insights into experimentally not accessible bubble parameters and shock wave emission. The model is extended by a term covering the initial shock-driven acceleration of the bubble wall, an automated method determining shock front position and pressure decay, and an energy balance for the partitioning of absorbed laser energy into vaporization, bubble and shock wave energy, and dissipation through viscosity and condensation. These tools are used for analysing a scattering signal covering 102 oscillation cycles. The bubble was produced by a plasma with 1550 K average temperature and had 36 $\mu$m maximum radius. Predicted bubble wall velocities during expansion agree well with experimental data. Upon first collapse, most energy was stored in the compressed liquid around the bubble and then radiated away acoustically. The collapsed bubble contained more vapour than gas, and its pressure was 13.5 GPa. The pressure of the rebound shock wave initially decayed $\propto r^{-1.8}$, and energy dissipation at the shock front heated liquid near the bubble wall above the superheat limit. The shock-induced temperature rise reduces damping during late bubble oscillations. Bubble dynamics changes significantly for small bubbles with less than 10 $\mu$m radius.

Comparative analysis of machine learning methods for active flow control

Abstract: Abstract Machine learning frameworks such as genetic programming and reinforcement learning (RL) are gaining popularity in flow control. This work presents a comparative analysis of the two, benchmarking some of their most representative algorithms against global optimization techniques such as Bayesian optimization and Lipschitz global optimization. First, we review the general framework of the model-free control problem, bringing together all methods as black-box optimization problems. Then, we test the control algorithms on three test cases. These are (1) the stabilization of a nonlinear dynamical system featuring frequency cross-talk, (2) the wave cancellation from a Burgers’ flow and (3) the drag reduction in a cylinder wake flow. We present a comprehensive comparison to illustrate their differences in exploration versus exploitation and their balance between ‘model capacity’ in the control law definition versus ‘required complexity’. Indeed, we discovered that previous RL control attempts of controlling the cylinder wake were performing linear control and that the wide observation space was limiting their performances. We believe that such a comparison paves the way towards the hybridization of the various methods, and we offer some perspective on their future development in the literature of flow control problems.

Scaling laws of the maximum spreading factor for impact of nanodroplets on solid surfaces

Abstract: Abstract This study investigates the dynamics of low-viscosity nanodroplets impacting surfaces with static contact angles from θ = 73° to 180° via molecular dynamics (MD) simulations. Two typical morphologies of impacting nanodroplets are observed at the maximum spreading state, a Hertz-ball-like in a low-Weber-number range and a thin-film-like in a high-Weber-number range. Only inertial and capillary forces dominate the impact for the former, whereas viscous force also becomes dominant for the latter. Regardless of morphologies at the maximum spreading state, the ratio of spreading time to contact time always remains constant on an ideal superhydrophobic surface with θ = 180°. With the help of different kinematic approximations of the spreading time and scaling laws of the contact time, scaling laws of the maximum spreading factor ${\beta _{max}}\sim W{e^{1/5}}$ in the low-Weber-number range (capillary regime) and ${\beta _{max}}\sim W{e^{2/3}}R{e^{ - 1/3}}$ (or ${\beta _{max}}\sim W{e^{1/2}}O{h^{1/3}}$) in the high-Weber-number range (cross-over regime) are obtained. Here, We, Re, and Oh are the Weber number, Reynolds number, and Ohnesorge number, respectively. Although the scaling laws are proposed only for the ideal superhydrophobic surface, they are tested valid for θ over 73° owing to the ignorable zero-velocity spreading effect. Furthermore, combining the two scaling laws leads to an impact number, $W{e^{3/10}}O{h^{1/3}} = 2.1$. This impact number can be used to determine whether viscous force is ignorable for impacting nanodroplets, thereby distinguishing the capillary regime from the cross-over regime.

Wall shear stress from jetting cavitation bubbles: influence of the stand-off distance and liquid viscosity

Abstract: We study systematically the cavitation-induced wall shear stress on rigid boundaries as a function of liquid viscosity is the distance between lower bubble interface and wall prior to impact.

A unifying theory of jet screech

Abstract: Abstract This paper describes the mechanism underpinning modal staging behaviour in screeching jets. An upstream-propagating subsonic guided-jet mode is shown to be active in all stages of screech. Axial variation of shock-cell spacing manifests in the spectral domain as a series of suboptimal peaks. It is demonstrated that the guided-jet mode may be energized by interactions of the Kelvin–Helmholtz wavepacket with not only the primary shock wavenumber peak, but also suboptimal peaks; interaction with suboptimals is shown to be responsible for closing the resonance loop in multiple stages of jet screech. A consideration of the full spectral representation of the shocks reconciles several of the classical models and results for jet screech that had heretofore been paradoxical. It is demonstrated that there are multiple standing waves present in the near field of screeching jets, corresponding to the superposition of the various waves active in these jets. Multimodal behaviour is explored for jets in a range of conditions, demonstrating that multiple peaks in the frequency spectra can be due to either changes in which peak of the shock spectra the Kelvin–Helmholtz wavepacket is interacting with, or a change in azimuthal mode, or both. The absence of modal staging in high-aspect-ratio non-axisymmetric jets is also explained in the context of the aforementioned mechanism. The paper closes with a new proposed theory for frequency selection in screeching jets, based on the observation that these triadic interactions appear to underpin selection of the guided-jet mode wavelength in all measured cases.

Absolute instability in shock-containing jets

Abstract: We present an analysis of the linear stability characteristics of shock-containing jets. The flow is linearised around a spatially periodic mean, which acts as a surrogate for a mean flow with a shock-cell structure, leading to a set of partial differential equations with periodic coefficients in space. Disturbances are written using the Floquet ansatz and Fourier modes in the streamwise direction, leading to an eigenvalue problem for the Floquet exponent. The characteristics of the solution are directly compared with the locally parallel case, and some of the features are similar. The inclusion of periodicity induces minor changes in the growth rate and phase velocity of the relevant modes for small shock amplitudes. On the other hand, the eigenfunctions are now subject to modulation related to the periodicity of the flow. Analysis of the spatio-temporal growth rates led to the identification of a saddle point between the Kelvin–Helmholtz mode and the guided jet mode, characterising an absolute instability mechanism. Frequencies and mode shapes related to the saddle points for two conditions (associated with axisymmetric and helical modes) are compared with screech frequencies and the most energetic coherent structures of screeching jets, resulting in a good agreement for both. The analysis shows that a periodic shock-cell structure has an impulse response that grows upstream, leading to oscillator behaviour. The results suggest that screech can occur in the absence of a nozzle, and that the upstream reflection condition is not essential for screech frequency selection. Connections to previous models are also discussed.

About the influences of compressibility, heat transfer and pressure gradients in compressible turbulent boundary layers

Abstract: This paper presents a comprehensive analysis of the momentum and energy transfer in compressible turbulent boundary layers based on integral identities. By considering data obtained from direct numerical simulations for a wide parameter range, the superordinate influences of compressibility, wall heat transfer and pressure gradient on the terms of the governing equations are identified and visualized. This allows us both to determine to what degree cases corresponding to different Mach number, heat transfer and pressure-gradient conditions have physically comparable behaviour and to design turbulent boundary-layer cases with specific sought-after behaviour. To this end, newly formulated identities for the skin-friction coefficient is proposed as a similarity parameter, mapping cases with different Mach numbers and wall heat transfer conditions.

On the role of nonlinear correlations in reduced-order modelling

Abstract: A major goal for reduced-order models of unsteady fluid flows is to uncover and exploit latent low-dimensional structure. Proper orthogonal decomposition (POD) provides an energy-optimal linear basis to represent the flow kinematics, but converges slowly for advection-dominated flows and tends to overestimate the number of dynamically relevant variables. We show that nonlinear correlations in the temporal POD coefficients can be exploited to identify the underlying attractor, characterized by a minimal set of driving modes and a manifold equation for the remaining modes. By viewing these nonlinear correlations as an invariant manifold reduction, this least-order representation can be used to stabilize POD-Galerkin models or as a state space for data-driven model identification. In the latter case, we use sparse polynomial regression to learn a compact, interpretable dynamical system model from the time series of the active modal coefficients. We demonstrate this perspective on a quasiperiodic shear-driven cavity flow and show that the dynamics evolve on a torus generated by two independent Stuart-Landau oscillators. These results emphasize the importance of nonlinear dimensionality reduction to reveal underlying structure in complex flows.

Residual dynamic mode decomposition: robust and verified Koopmanism

Abstract: Abstract Dynamic mode decomposition (DMD) describes complex dynamic processes through a hierarchy of simpler coherent features. DMD is regularly used to understand the fundamental characteristics of turbulence and is closely related to Koopman operators. However, verifying the decomposition, equivalently the computed spectral features of Koopman operators, remains a significant challenge due to the infinite-dimensional nature of Koopman operators. Challenges include spurious (unphysical) modes and dealing with continuous spectra, which both occur regularly in turbulent flows. Residual dynamic mode decomposition (ResDMD), introduced by Colbrook & Townsend (Rigorous data-driven computation of spectral properties of Koopman operators for dynamical systems. 2021. arXiv:2111.14889), overcomes such challenges through the data-driven computation of residuals associated with the full infinite-dimensional Koopman operator. ResDMD computes spectra and pseudospectra of general Koopman operators with error control and computes smoothed approximations of spectral measures (including continuous spectra) with explicit high-order convergence theorems. ResDMD thus provides robust and verified Koopmanism. We implement ResDMD and demonstrate its application in various fluid dynamic situations at varying Reynolds numbers from both numerical and experimental data. Examples include vortex shedding behind a cylinder, hot-wire data acquired in a turbulent boundary layer, particle image velocimetry data focusing on a wall-jet flow and laser-induced plasma acoustic pressure signals. We present some advantages of ResDMD: the ability to resolve nonlinear and transient modes verifiably; the verification of learnt dictionaries; the verification of Koopman mode decompositions; and spectral calculations with reduced broadening effects. We also discuss how a new ordering of modes via residuals enables greater accuracy than the traditional modulus ordering (e.g. when forecasting) with a smaller dictionary. This result paves the way for more significant dynamic compression of large datasets without sacrificing accuracy.

Vortex dynamics and hydrodynamic performance enhancement mechanism in batoid fish oscillatory swimming

Abstract: The effects of chordwise deformation and the half-amplitude asymmetry on the hydrodynamic performance and vortex dynamics of batoid fish have been numerically investigated, in which the two parameters were represented by the wavenumber ( was also due to the increase of the vortical contributions. This work provides a new insight into the role of vortices and the force enhancement mechanism in aquatic swimming.

On the polydisperse particle migration and formation of chains in a square channel flow of non-Newtonian fluids

Abstract: Abstract The migration of polydisperse particles and the formation of self-organized particle chains in a square channel flow of non-Newtonian fluids is studied. The effects of rheological behaviour of the fluid, solution concentration and flow rate are explored experimentally. The direct forcing/fictitious domain method is adopted to qualitatively verify the experiments and further analyse the mechanisms of particle migration and particle chain self-organization. The results show that only particles in viscoelastic fluids with negligible shear-thinning effect will remain at the channel centreline as the flow rate increases. The monodisperse particles reach the same velocity when migrating to the equilibrium position. However, in polydisperse suspensions, the smaller the particle diameter, the greater the velocity when the particle migrates to the equilibrium position. In a viscoelastic fluid, the polydisperse particles are more likely to self-organize into long particle chains along the channel centreline than the monodisperse particles, where the large and small particles are at the front and end of the chain. The dimensionless alignment factor (Af) is adopted to quantify the formation of particle chains, which is the largest in viscoelastic fluids and rapidly increases before decreasing to a stable value as the flow rate increases. For larger particle diameter ratios and stronger shear-thinning effect, the long particle chain self-organization is less obvious. The self-organizing particle chains at the channel centreline are strongly influenced by the fluid elastic properties and weakly by the inertial effect; however, the shear-thinning effect disperses the particles and prevents the formation of long straight particle chains.

Stability and tip streaming of a surfactant-loaded drop in an extensional flow. Influence of surface viscosity

Abstract: Abstract We study numerically the nonlinear stationary states of a droplet covered with an insoluble surfactant in a uniaxial extensional flow. We calculate both the eigenfunctions to reveal the instability mechanism and the time-dependent states resulting from it, which provides a coherent picture of the phenomenon. The transition is of the saddle-node type, both with and without surfactant. The flow becomes unstable under stationary linear perturbations. Surfactant considerably reduces the interval of stable capillary numbers. Inertia increases the droplet deformation and decreases the critical capillary number. In the presence of the surfactant monolayer, neither the droplet deformation nor the stability is significantly affected by the droplet viscosity. The transient state resulting from instability is fundamentally different for drops with and without surfactant. Tip streaming occurs only in the presence of surfactants. The critical eigenmode leading to tip streaming is qualitatively the same as that yielding the central pinching mode for a clean interface, which indicates that the small local scale characterizing tip streaming is set during the nonlinear droplet deformation. The viscous surface stress does not significantly affect the droplet deformation and the critical capillary number. However, the damping rate of the dominant mode considerably decreases for viscous surfactants. Interestingly, shear viscous surface stress considerably alters the tip streaming arising in the supercritical regime, even for very small surface viscosities. The viscous surface stresses alter the balance of normal interfacial stresses and affect the surfactant transport over the stretched interface.

Analytical solutions for reactive shear dispersion with boundary adsorption and desorption

Abstract: Abstract Surface reactions such as the adsorption and desorption at boundaries are very common for solute dispersion in many applications of chemistry, biology, hydraulics, etc. To study how reversible adsorption affects the transient dispersion, Zhang, Hesse & Wang (J. Fluid Mech., vol. 828, 2017, pp. 733–752) have investigated the temporal evolution of moments using the Laplace transform method. Owing to difficulties introduced by the adsorption–desorption boundary condition, great challenges arise from the inverse Laplace transform: dealing with the singularities by the residue theorem can tremendously increase complexities. This work provides a much simpler analytical method to derive solutions in a more compact form that is valid for the entire range of the reactive transport process. Such a progress demonstrates that the classic framework of separation of variables can be extended and applied to this more general adsorption–desorption condition, based on which higher-order statistics including skewness and kurtosis can be explicitly explored in practice. Also extended is Gill's generalised dispersion model for solute concentration distributions, which can now address the entire transient dispersion characteristics, instead of just applied for the long-time asymptotic reactive process as done previously. Regarding the most classic Taylor dispersion problem, we investigate the influence of the reversible adsorption–desorption on the solute cloud in a tube flow. Not only the transient dispersion characteristics of transverse-average concentration distribution but also those of the bulk, surface and total-average distributions are discussed. We further investigate the influence of initial conditions on the non-uniformity of the transient dispersion over the cross-section.

Swimming dynamics of a self-propelled droplet

Abstract: Abstract Chemically active droplets often show intriguing self-propulsion behaviour in a surfactant solution. The drop motion is controlled by the nonlinear coupling among chemical transport in the bulk fluid, consumption of surfactant at the drop surface, and the fluid flow driven by the self-generated Marangoni stress. To quantify the underlying hydrodynamics, this work investigates the swimming motion of a two-dimensional drop that is determined by two dimensionless parameters, the Péclet number ($Pe$) and Damköhler number ($Da$). The weakly nonlinear analysis shows that near the instability threshold, the drop undergoes a supercritical bifurcation with velocity $U\sim \sqrt {Pe-Pe_c}$, where $Pe_c$ is the critical Péclet number for the onset of dipole mode. In the highly nonlinear regime, the drop transits from steady translation of pusher swimming to unsteady motion of mixed pusher–puller swimming along zigzaging trajectories of quadrangle and/or triangle waves. Mode decomposition shows that the zigzag motion is directly related to the interaction between the secondary wake of low surfactant concentration and the primary wake.

On the role of turbulent large-scale streaks in generating sediment ridges

Abstract: The role of turbulent large-scale streaks or large-scale motions in forming subaqueous sediment ridges on an initially flat sediment bed is investigated with the aid of particle resolved direct numerical simulations of open channel flow at bulk Reynolds numbers up to 9500. The regular arrangement of quasi-streamwise ridges and troughs at a characteristic spanwise spacing between 1 and 1.5 times the mean fluid height is found to be a consequence of the spanwise organisation of turbulence in large-scale streamwise velocity streaks. Ridges predominantly appear in regions of weaker erosion below large-scale low-speed streaks and vice versa for troughs. The interaction between the dynamics of the large-scale streaks in the bulk flow and the evolution of sediment ridges on the sediment bed is best described as ‘top-down’ process, as the arrangement of the sediment bedforms is seen to adapt to changes in the outer flow with a time delay of several bulk time units. The observed ‘top-down’ interaction between the outer flow and the bed agrees fairly well with the conceptual model on causality in canonical channel flows proposed by Jimenez (J. Fluid Mech., vol. 842, 2018, P1, § 5.6). Mean secondary currents of Prandtl's second kind of comparable intensity and lateral spacing are found over developed sediment ridges and in single-phase smooth-wall channels alike in averages over bulk time units. This indicates that the secondary flow commonly observed together with sediment ridges is the statistical footprint of the regularly organised large-scale streaks.

Modulation of homogeneous and isotropic turbulence in emulsions

Abstract: We present a numerical study of emulsions in homogeneous and isotropic turbulence at $Re_\lambda=137$. The problem is addressed via Direct Numerical Simulations (DNS), where the Volume of Fluid (VOF) is used to represent the complex features of the liquid-liquid interface. We consider a mixture of two iso-density fluids, where fluid properties are varied with the goal of understanding their role in turbulence modulation, in particular the volume fraction ($0.03<\alpha<0.5$), viscosity ratio ($0.01<\mu_d/\mu_c<100$) and large scale Weber number ($10.6

Electro-fluid-mechanics of the heart

Abstract: Abstract This article presents an overview of the dynamics of the human heart and the main goal is the discussion of its fluid mechanic features. We will see, however, that the flow in the heart can not be fully described without considering its electrophysiology and elastomechanics as well as the interaction with the systemic and pulmonary circulations with which it is strongly connected. Biologically, the human heart is similar to that of all warm-blooded mammals and it satisfies the same allometric laws. Since the Paleolithic Age, however, humans have improved their living conditions, have modified the environment to satisfy their needs and, more recently, have developed advanced medical knowledge which has allowed triple the number of heartbeats with respect to other mammals. In the last century, effective diagnostic tools, reliable surgical procedures and prosthetic devices have been developed and refined leading to substantial progress in cardiology and heart surgery with routine clinical practice which nowadays cures many disorders, once lethal. Pulse duplicators have been built to reproduce the pulsatile flow and ‘blood analogues’, have been realized. Heart phantoms, can attain deformations similar to the real heart although the active contraction and the tissue anisotropy still can not be replicated. Numerical models have also become a viable alternative for cardiovascular research: they do not suffer from limitations of material properties and device technologies, thus making possible the realization of truly digital twins. Unfortunately, a high-fidelity model for the whole heart consists of a system of coupled, nonlinear partial differential equations with a number of degrees of freedom of the order of a billion and computational costs become the bottleneck. An additional challenge comes from the inherent human variability and the uncertainty of the heart parameters whose statistical assessment requires a campaign of simulations rather than a single deterministic calculation; reduced and surrogate models can be employed to alleviate the huge computational burden and all possibilities are currently being pursued. In the era of big data and artificial intelligence, cardiovascular research is also advancing by exploiting the latest technologies: equation-based augmented reality, virtual surgery and computational prediction of disease progression are just a few examples among many that will become standard practice in the forthcoming years.

Weakly nonlinear analysis of the viscoelastic instability in channel flow for finite and vanishing Reynolds numbers

Abstract: The recently discovered centre-mode instability of rectilinear viscoelastic shear flow (Garg et al. , Phys. Rev. Lett. , vol. 121, 2018, 024502) has offered an explanation for the origin of elasto-inertial turbulence that occurs at lower Weissenberg numbers ( $Wi$ ). In support of this, we show using weakly nonlinear analysis that the subcriticality found in Page et al. ( Phys. Rev. Lett. , vol. 125, 2020, 154501) is generic across the neutral curve with the instability becoming supercritical only at low Reynolds numbers ( $Re$ ) and high $Wi$ . We demonstrate that the instability can be viewed as purely elastic in origin, even for $Re=O(10^3)$ , rather than ‘elasto-inertial’, as the underlying shear does not feed the kinetic energy of the instability. It is also found that the introduction of a realistic maximum polymer extension length, $L_{max}$ , in the FENE-P model moves the neutral curve closer to the inertialess $Re=0$ limit at a fixed ratio of solvent-to-solution viscosities, $\beta$ . At $Re=0$ and in the dilute limit ( $\beta \rightarrow 1$ ) with $L_{max} =O(100)$ , the linear instability can be brought down to more physically relevant $Wi\gtrsim 110$ at $\beta =0.98$ , compared with the threshold $Wi=O(10^3)$ at $\beta =0.994$ reported recently by Khalid et al. ( Phys. Rev. Lett. , vol. 127, 2021, 134502) for an Oldroyd-B fluid. Again, the instability is subcritical, implying that inertialess rectilinear viscoelastic shear flow is nonlinearly unstable – i.e. unstable to finite-amplitude disturbances – for even lower $Wi$ .

The principal stage in wind-wave generation

Abstract: Abstract The dynamics of wind-generated water waves in the principal stage of the Phillips theory (Phillips, J. Fluid Mech., vol. 2, 1957, pp. 417–445) is investigated by a combined numerical and analytical approach. We perform a number of high-resolution direct numerical simulation (DNS) of turbulent wind over initially calm water to capture the multistage generation of water waves. Detailed analyses are conducted to evaluate the Phillips theory in both physical space and wavenumber space. Numerical evidence is obtained for the existence of a principal stage when the surface elevation variance grows linearly with time. We further propose a random sweeping turbulence pressure–wave interaction model by introducing the random sweeping hypothesis of air pressure fluctuations to the Phillips theory, and obtain an asymptotic solution of the mean square of surface wave elevations over time. This asymptotic analysis captures the temporal oscillations of surface elevation variance in the principal stage, which is also confirmed by our DNS results. The wavenumber spectrum of surface wave elevations is analysed using a time-dependent norm to elucidate the role of the resonance mechanism on wave generation. In physical space, we use the random sweeping turbulence pressure–wave interaction model to obtain a quantitative prediction of the growth rate of surface elevation variance in the principal stage, which is found to agree with the DNS results better than the original Phillips model.

Boundary layer transition and linear modal instabilities of hypersonic flow over a lifting body

Abstract: Abstract Boundary layer transition over a lifting body of 1.6 m length at $2^\circ$ angle of attack has been simulated at Mach 6 and a unit Reynolds number $1.0 \times 10^7$ m$^{-1}$. The model geometry is the same as the Hypersonic Transition Research Vehicle designed by the China Aerodynamics Research and Development Center. Four distinct transitional regions are identified, i.e. windward vortex region, shoulder vortex region, windward cross-flow region and shoulder cross-flow region. Multi-dimensional linear stability analyses by solving the two-dimensional eigenvalue problem (spatial BiGlobal approach) and the plane-marching parabolized stability equations (PSE3D approach) are further carried out to uncover the dominant instabilities in the last three regions as well as the shoulder attachment-line region. The shoulder vortex is conducive to both inner and outer modes of shear-layer instability, of which the latter most likely trigger the vortex breakdown. A novel method is presented to substantially reduce the computational cost of BiGlobal and PSE3D in resolving the cross-flow instabilities in cross-flow regions. The peak frequencies of cross-flow modes lie between 15 and 45 kHz. Whereas oblique second Mack modes are marginally unstable in the windward cross-flow region, they could be strong enough to compete with the cross-flow modes in the shoulder cross-flow region. In the shoulder attachment-line region, there exists only one unstable mode of Mack instability, differing from previous studies that show a hierarchy of modes in the context of symmetrical attachment-line flows. Results of the numerical simulation and multi-dimensional stability analyses are compared when possible, showing a fair agreement between the two approaches and highlighting the necessity of considering non-parallel effects.

Pressure fluctuations due to ‘trapped waves’ in the initial region of compressible jets

Abstract: An experimental study is conducted on unsteady pressure fluctuations occurring near the nozzle exit and just outside the shear layer of compressible jets. These fluctuations are related to ‘trapped waves’ within the jet's potential core, as investigated and reported recently by other researchers. Round nozzles of three different diameters and rectangular nozzles of various aspect ratios are studied. The fluctuations manifest as a series of peaks in the spectra of the fluctuating pressure. Usually the first peak at the lowest frequency (fundamental) has the highest amplitude and the amplitude decreases progressively for successive peaks at higher frequencies. These ‘trapped wave spectral peaks’ are found to occur with all jets at high subsonic conditions and persist into the supersonic regime. Their characteristics and variations with axial and radial distances, jet Mach number and aspect ratio of the nozzle are documented. For round nozzles, the frequency of the fundamental is found to be independent of the jet's exit boundary layer characteristics and scales with the nozzle diameter. On a Strouhal number (based on diameter) versus jet Mach number plot it is represented by a unique curve. Relative to the fundamental the frequencies of the successive peaks are found to bear the ratios of 5/3, 7/3, 9/3 and so on, at a given Mach number. For rectangular nozzles, the number of peaks observed on the major axis is found to be greater than that observed on the minor axis by a factor approximately equal to the nozzle's aspect ratio; the fundamental is the same on either edge. For all nozzles the onset of screech tones appears as a continuation of the evolution of these peaks; it is as if one of these peaks abruptly increases in amplitude and turns into a screech tone as the jet Mach number is increased.

Transition, intermittency and phase interference effects in airfoil secondary tones and acoustic feedback loop

Abstract: A large eddy simulation is performed to study secondary tones generated by a NACA0012 airfoil at angle of attack of $\alpha = 3^{\circ}$ with freestream Mach number of $M_{\infty} = 0.3$ and Reynolds number of $Re = 5 \times 10^4$. Laminar separation bubbles are observed over the suction side and near the trailing edge, on the pressure side. Flow visualization and spectral analysis are employed to investigate vortex shedding aft of the suction side separation bubble. Vortex interaction results in merging or bursting such that coherent structures or turbulent packets are advected towards the trailing edge leading to different levels of noise emission. Despite the intermittent occurrence of laminar-turbulent transition, the noise spectrum depicts a main tone with multiple equidistant secondary tones. To understand the role of flow instabilities on the tones, the linearized Navier-Stokes equations are examined in its operator form through bi-global stability and resolvent analyses, and by time evolution of disturbances using a matrix-free method. These linear global analyses reveal amplification of disturbances over the suction side separation bubble. Non-normality of the linear operator leads to further transient amplification due to modal interaction among eigenvectors. Two-point, one time autocovariance calculations of pressure along the spanwise direction elucidate aspects of the acoustic feedback loop mechanism in the non-linear solutions. This feedback process is self-sustained by acoustic waves radiated from the trailing edge, which reach the most sensitive flow location between the leading edge and the separation bubble, as identified by the resolvent analysis. Leading edge disturbances arising from secondary diffraction and phase interference among the most unstable frequencies computed in the eigenspectrum are also shown to have an important role in the feedback loop.

A one-dimensional model for the pulsating flow of cerebrospinal fluid in the spinal canal

Abstract: Abstract The monitoring of intracranial pressure (ICP) fluctuations, which is needed in the context of a number of neurological diseases, requires the insertion of pressure sensors, an invasive procedure with considerable risk factors. Intracranial pressure fluctuations drive the wave-like pulsatile motion of cerebrospinal fluid (CSF) along the compliant spinal canal. Systematically derived simplified models relating the ICP fluctuations with the resulting CSF flow rate can be useful in enabling indirect evaluations of the former from non-invasive magnetic resonance imaging (MRI) measurements of the latter. As a preliminary step in enabling these predictive efforts, a model is developed here for the pulsating viscous motion of CSF in the spinal canal, assumed to be a linearly elastic compliant tube of slowly varying section, with a Darcy pressure-loss term included to model the fluid resistance introduced by the trabeculae, which are thin collagen-reinforced columns that form a web-like structure stretching across the spinal canal. Use of Fourier-series expansions enables predictions of CSF flow rate for realistic anharmonic ICP fluctuations. The flow rate predicted using a representative ICP waveform together with a realistic canal anatomy is seen to compare favourably with in vivo phase-contrast MRI measurements at multiple sections along the spinal canal. The results indicate that the proposed model, involving a limited number of parameters, can serve as a basis for future quantitative analyses targeting predictions of ICP temporal fluctuations based on MRI measurements of spinal-canal anatomy and CSF flow rate.

Characterization of instability mechanisms on sharp and blunt slender cones at Mach 6

Abstract: Abstract Experiments are performed to investigate the effect of nose-tip bluntness on the instability mechanisms leading to boundary-layer transition on a $7^{\circ }$ half-angle cone in a Mach-6 free stream. The development of disturbances is characterized using a combination of high-speed calibrated schlieren images and pressure measurements, and the data are compared with results computed using the parabolized stability equations. The approximately 414 mm long cone model is equipped with an interchangeable nose tip ranging from sharp to 5.08 mm in radius. For nose tips with a radius $R_{N}<2.54\ {\rm mm}$, second-mode instability waves are the dominant mechanism leading to transition. Time-averaged frequency spectra computed from the calibrated schlieren visualizations and pressure measurements are used to compute the second-mode most-amplified frequencies and integrated amplification rates ($N$ factors). Good agreement is observed between the measurements and computations in the linear-growth regime for the sharp-nose configuration at each free-stream condition. Additionally, a bispectral analysis identifies quadratic phase locking of frequency content responsible for the growth of higher harmonics. For nose tips of $R_{N}\geqslant 2.54\ {\rm mm}$, the schlieren visualization region is upstream of the entropy-layer swallowing length, and second-mode waves are no longer visible within the boundary layer; instead, elongated, steeply inclined features believed to be associated with non-modal instability mechanisms develop between the entropy-layer and boundary-layer edges. Simultaneously acquired surface pressure measurements reveal high-frequency pressure oscillations similar to second-mode instability waves associated with the trailing edge of these non-modal features.

A transformer-based synthetic-inflow generator for spatially developing turbulent boundary layers

Abstract: Abstract This study proposes a newly developed deep-learning-based method to generate turbulent inflow conditions for spatially developing turbulent boundary layer (TBL) simulations. A combination of a transformer and a multiscale-enhanced super-resolution generative adversarial network is utilised to predict velocity fields of a spatially developing TBL at various planes normal to the streamwise direction. Datasets of direct numerical simulation (DNS) of flat plate flow spanning a momentum thickness-based Reynolds number, $Re_\theta = 661.5\unicode{x2013}1502.0$, are used to train and test the model. The model shows a remarkable ability to predict the instantaneous velocity fields with detailed fluctuations and reproduce the turbulence statistics as well as spatial and temporal spectra with commendable accuracy as compared with the DNS results. The proposed model also exhibits a reasonable accuracy for predicting velocity fields at Reynolds numbers that are not used in the training process. With the aid of transfer learning, the computational cost of the proposed model is considered to be effectively low. Furthermore, applying the generated turbulent inflow conditions to an inflow–outflow simulation reveals a negligible development distance for the TBL to reach the target statistics. The results demonstrate for the first time that transformer-based models can be efficient in predicting the dynamics of turbulent flows. They also show that combining these models with generative adversarial networks-based models can be useful in tackling various turbulence-related problems, including the development of efficient synthetic-turbulent inflow generators.