Search

The highly nonlinear evolution of the single-mode stratified compressible Rayleigh–Taylor instability (RTI) is investigated via direct numerical simulation over a range of Atwood numbers ( $A_T=0.1$ – $0.9$ ) and Mach numbers ( $Ma=0.1$ – $0.7$ ) for characterising the isothermal background stratification. After the potential stage, it is found that the bubble is accelerated to a velocity which is well above the saturation value predicted in the potential flow model. Unlike the bubble re-acceleration behaviour in quasi-incompressible RTI with uniform background density, the characteristics in the stratified compressible RTI are driven by not only vorticity accumulation inside the bubble but also flow compressibility resulting from the stratification. Specifically, in the case of strong stratification and high $A_T$ , the flow compressibility dominates the bubble re-acceleration characters. To model the effect of flow compressibility, we propose a novel model to reliably describe the bubble re-acceleration behaviours in the stratified compressible RTI, via introducing the dilatation into the classical model that takes into account only vorticity accumulation.

Three-dimensional elastic turbulence in Taylor–Couette flows of dilute polymer solutions has been realized and thoroughly investigated via direct numerical simulations. A novel flow transition pathway from elastically dominated turbulence to solitary vortex pairs (or diwhirls) and eventually to elastic turbulence is observed by decreasing the fluid inertia ( $Re$ ) over seven orders of magnitude, i.e. from $Re=1000$ to $0.0001$ . The dominant spatio-temporal flow features in the elastic turbulence regime are those of large-scale unsteady diwhirls and small-scale axial and azimuthal travelling waves in the outer and inner halves of the gap, respectively. Moreover, it is conclusively shown that production of turbulent kinetic energy in purely elastic turbulence solely arises due to the stochastic nature of polymer stretch/relaxation. Overall, based on this comprehensive numerical investigation, the differences in the underlying fluid physics that give rise to turbulent fluctuations in elastically dominated and purely elastic turbulence have been delineated.

The width effect on contact angle hysteresis in a microchannel with patterned heterogeneous surfaces is systematically investigated. In the model, identical defects periodically appear on the background surface. To this end, a droplet's evaporation and condensation processes inside the microchannel are studied by theoretical analysis and numerical simulation based on a diffuse-interface lattice Boltzmann method. The microchannel width effect on the system's equilibrium properties is studied. The results demonstrate that the number of equilibrium configurations increases linearly with the microchannel width ( $b$ ), and has a quadratic relationship with the cosine of the reference contact angle and the heterogeneity strength ( $\varepsilon$ ). The average most stable contact angle is independent of $b$ and is always equal to the contact angle predicted by the Cassie–Baxter equation. For contact angle hysteresis ( $H$ ), when the microchannels are narrow and wide, there are individual-effect-dominated hysteresis (IDH) and collective-effect-dominated hysteresis (CDH), respectively. The IDH and CDH are hysteresis modes corresponding to the jumping behaviour of contact lines affected by individual defects and two neighbouring defects, respectively. Based on the graphical force balance approach, we establish a scaling law to quantify the connection between $H$ , $b$ and $\varepsilon$ . Specifically, in the IDH mode, $H\sim b \varepsilon ^2$ , while in the CDH mode, $H$ increases linearly with $\varepsilon$ but nonlinearly with $b$ .

Pressure fluctuations play an essential role in the transport of turbulent kinetic energy and vibrational loading. This study focuses on examining the effect of wall cooling on pressure fluctuations in compressible turbulent boundary layers by high-fidelity direct numerical simulations. Pressure fluctuations result from the vorticity mode and the acoustic mode that are both closely dependent on compressibility. To demonstrate the effects of wall cooling at various compressibility intensities, three free-stream Mach numbers are investigated, i.e. $M_\infty =0.5$ , 2.0 and 8.0, with real gas effects being absent for $M_\infty =8.0$ due to a low enthalpy inflow. Overall, opposite effects of wall cooling on pressure fluctuations are found between the subsonic/supersonic cases and the hypersonic case. Specifically, the pressure fluctuations normalized by wall shear stress $p^\prime _{rms}/\tau _w$ are suppressed in the subsonic and supersonic cases, while enhanced in the hypersonic case near the wall. Importantly, travelling-wave-like alternating positive and negative structures (APNS), which greatly contribute to pressure fluctuations, are identified within the viscous sublayer and buffer layer in the hypersonic cases. Furthermore, generating mechanisms of pressure fluctuations are explored by extending the decomposition based on the fluctuating pressure equation to compressible turbulent boundary layers. Pressure fluctuations are decomposed into five components, in which rapid pressure, slow pressure and compressible pressure are dominant. The suppression of pressure fluctuations in the subsonic and supersonic cases is due to both rapid pressure and slow pressure being suppressed by wall cooling. In contrast, wall cooling strengthens compressible pressure for all Mach numbers, especially in the hypersonic case, resulting in increased wall pressure fluctuations. Compressible pressure plays a leading role in the hypersonic case, mainly due to the APNS. Essentially, the main effects of wall cooling can be interpreted by the suppression of the vorticity mode and the enhancement of the acoustic mode.

Trailing-edge serrations inspired by owls are capable of reducing broadband noise. In this study, the wall-resolved large-eddy simulations (LES) are carried out on the flow over NACA-0012 airfoil with additional serrated trailing edges. The computations are conducted with the high-order flux reconstruction method on unstructured meshes. Three kinds of serrations with different lengths are studied and compared with the straight trailing-edge case, and all three types of serration achieved a certain degree of noise reduction. Presently, the medium-length serration achieves the best noise reduction effect. The maximum decrease of overall sound pressure level is approximately 2.4 dB, implying that the length of serration has a substantial impact on the noise reduction effect. The serration has no significant effect on the upstream turbulence statistics, but it changes the flow structure near the serration, such as inducing side vortex pairs attached to the serration edges. Moreover, dynamic mode decomposition shows that the pressure structures vary with the serration length. For the most unstable hydrodynamic wave, the spanwise coherence of the mode structure of pressure in the upstream boundary layer is weakened. In addition, serrations can redistribute the dipole sources on the surfaces of airfoil and serrations. The destructive interference is enhanced to some extent, which is favourable for noise reduction. In contrast with LES simulations, the pure dipole analysis shows that the longest serration case seems to be the best. Furthermore, a recently developed noise theory is used to evaluate the influence of serrations on the flow noise sources qualitatively and quantitatively. It is found that the serrations can mitigate noise source intensity near the serration edges but increase the source intensity in the near wake. The combined effect of serration on the dipole source and flow noise source determines the overall noise reduction effect. To conclude, destructive interference plays a primary role in suppressing noise radiation by serrated trailing edges, and the dual effect of flow noise sources should be considered in future serration designs. As the influence of turbulence structure will make it more difficult to find the optimal serration parameters, the position of high-fidelity simulation will become increasingly important.

We investigate by direct numerical simulation the active control of laminar-turbulent transition in a hypersonic flat-plate boundary layer at a freestream Mach number of 5.86. The control mechanism is a synthetic jet. Based upon the linear stability theory of Mack, in hypersonic flow the important path to transition involves a high-frequency, second-mode fundamental resonance. Through systematic investigation, we reveal that the forcing the boundary layer with a synthetic jet at appropriate combinations of amplitude and frequency suppresses the second mode and delays transition. To gain physical insights into the major control mechanism, we employ the momentum potential theory (MPT) to analyze the flows with and without control. Essentially, the underlying control mechanism relies on an intriguing effect of the synthetic jet via generating the outward radiated wave structures, which are identified to split the upstream acoustic and vortical components. The splitting treatment presents the second-mode energy to drop sharply after the flow passes through the synthetic jet slot. The MPT source-term analysis reveals that the significantly suppressed near-wall source terms are responsible for suppressing the second mode downstream. Compared with the vortical and thermal source terms, the acoustic source term is found to be suppressed most. The kinetic budget analysis further reveals that the splitting treatment is related to the non-parallel effect and the nonlinear interaction.

A direct numerical simulation database of a weakly compressible turbulent channel flow with bulk Mach number 1.56 is studied in detail, including the geometrical relationships between the pressure-Hessian tensor and the vorticity/strain-rate tensor, as well as the mechanism of the pressure-Hessian tensor contributing to the evolution of invariants of the velocity gradient tensor. The results show that the geometrical relationships between the pressure-Hessian tensor and the vorticity/strain-rate tensor in the central region of the channel are consistent with that of isotropic turbulence. However, in the buffer layer with relatively stronger inhomogeneity and anisotropy, the vorticity tends to be aligned with the first or second eigenvector of the pressure-Hessian tensor in the unstable focus/compressing topological region, and tends to be aligned with the first eigenvector of the pressure-Hessian tensor in the stable focus/stretching topological region. In the unstable node/saddle/saddle and stable node/saddle/saddle topological regions, the vorticity prefers to lie in the plane of the first and second eigenvectors of the pressure-Hessian tensor. The strain-rate and the pressure-Hessian tensors tend to share their second principal direction. Moreover, for the coupling between the pressure-Hessian tensor and the principal strain rates, we clarify the influence on dissipation, the nonlinear generation of dissipation and the enstrophy generation. The decomposition of the pressure-Hessian tensor further shows that the slow pressure-related term dominates the pressure-Hessian tensor's contribution, and the influence of inhomogeneity and anisotropy mainly originates from the inhomogeneity and anisotropy of the fluctuating velocity. These statistical properties would be instructive in formulating dynamical models of the velocity gradient tensor for wall turbulence.

The evaporation mechanism of miscible binary nanodroplets from heated homogeneous surfaces was studied by molecular dynamics simulations, which has never been studied before. The binary droplets contain a hydrophilic component (type-2 particles) and a hydrophobic component (type-3 particles). It is shown that liquid-liquid interaction strength (e23) and hydrophilic particle number fraction (φ) have great influence on the surface tension, wetting characteristics, evaporation patterns, evaporation rate, and local mass flux. It is observed that when e23 ≥ 1, or φ ≈ 0.5, the evaporation mode is the constant-contact-angle mode. Otherwise, it is the mixed mode. We found that the evaporation rate becomes faster when φ and e23 increase. The droplets become more hydrophilic when φ increases, which promotes heat transfer efficiency between the liquid-solid interface. Besides, a larger e23 promotes the heat transfer inside the droplet. The mass transfer to the vapor phase occurs preferentially in the vicinity of TPCL (three phase contact line) in the hydrophilic systems (θ θc), the mass flux close to the TPCL is suppressed. We found that θc ∈ (102°-106°), which is different from the theoretical one, θc = 90°. The discrepancy is attributed to the existence of the adsorption layer near the TPCL.

Direct numerical simulation of spanwise-rotation-driven flow transitions in viscoelastic plane Couette flow from a drag-reduced inertial to a drag-enhanced elasto-inertial turbulent flow state followed by full relaminarization is reported for the first time. Specifically, this novel flow transition begins with a drag-reduced inertial turbulent flow state at a low rotation number $0\leqslant Ro \leqslant 0.1$ , and then transitions to a rotation/polymer-additive-driven drag-enhanced inertial turbulent regime, $0.1\leqslant Ro \leqslant 0.3$ . In turn, the flow transitions to a drag-enhanced elasto-inertial turbulent state, $0.3\leqslant Ro \leqslant 0.9$ , and eventually relaminarizes at $Ro=1$ . In addition, two novel rotation-dependent drag enhancement mechanisms are proposed and substantiated. (1) The formation of large-scale roll cells results in enhanced convective momentum transport along with significant polymer elongation and stress generated in the extensionally dominated flow between adjacent roll cells at $Ro\leqslant 0.2$ . (2) Coriolis-force-generated turbulent vortices cause strong incoherent transport and homogenization of significant polymer stress in the bulk via their vortical circulations at $Ro=0.5 - 0.9$ .

The nonlinear evolution of mixing layer in cylindrical Rayleigh-Taylor (RT) turbulence is studied theoretically and numerically. The scaling laws including the hyperbolic cosine growth for outward mixing layer and the cosine growth for inward mixing layer of the cylindrical RT turbulence are proposed for the first time and verified reliably by direct numerical simulation of the Navier-Stokes equations. It is identified that the scaling laws for the cylindrical RT turbulence transcend the classical power law for the planar RT turbulence and can be recovered to the quadratic growth as cylindrical geometry effect vanishes. Further, characteristic time- and length scales are reasonably obtained based on the scaling laws to reveal the self-similar evolution features for the cylindrical RT turbulence.

The stability of the two-layer film flow driven by an oscillatory plate under long-wave disturbances is studied. The influence of key factors, such as thickness ratio ( $n$ ), viscosity ratio ( $m$ ), density ratio ( $r$ ), oscillatory frequency ( $\beta$ ) and insoluble surfactants on the stability behaviours is studied systematically. Four special Floquet patterns are identified, and the corresponding growth rates are obtained by solving the eigenvalue problem of the fourth-order matrix. A small viscosity ratio ( $m\le 1$ ) may stabilize the flow but it depends on the thickness ratio. If the viscosity ratio is not very small ( $m>0.1$ ), in the $(\beta ,n)$ -plane, stable and unstable curved stripes appear alternately. In other words, under the circumstances, if the two-layer film flow is unstable, slightly adjusting the thickness of the upper film may make it stable. In particular, if the upper film is thin enough, even under high-frequency oscillation, the flow is always stable. The influence of density ratio is similar, i.e. there are curved stable and unstable stripes in the $(\beta ,r)$ -planes. Surface surfactants generally stabilize the flow of the two-layer oscillatory membrane, while interfacial surfactants may stabilize or destabilize the flow but the effect is mild. It is also found that gravity can generally stabilize the flow because it narrows the bandwidth of unstable frequencies.

The stability of the inverted flexible plate with non-uniform stiffness distribution in a free stream is studied by numerical simulation and mathematical theory. In our study, the bending stiffness distribution is expressed as the function of the leading edge's bending stiffness [Formula: see text] and the polynomial of the plate's coordinate. Based on the former theoretical work on the stability of inverted plates with uniform stiffness distribution, we derive the upper limit value of [Formula: see text] at which the zero-deflection equilibrium loses its stability for the plate with non-uniform stiffness distribution. The critical [Formula: see text] derived from the mathematical theory agrees well with that obtained from the numerical simulation. An effective bending stiffness is defined, which can be used to unify the regimes of the motion modes between uniform plates and non-uniform plates. Moreover, three orders of mass ratio [[Formula: see text], and [Formula: see text]] are investigated, and the underlying mechanism for large amplitude flapping is clarified for the inverted plate with different mass ratios. An appropriate bending stiffness distribution can greatly improve the deformation of the plate. The findings shed some light on the energy harvesting of the inverted plate.

When controlling the trailing-edge (TE) interference noise of airfoil, the design of the TE serration shape is still an open issue. To this end, the flow and noise generation for different TE serration shapes are explored by the wall-resolved implicit large-eddy simulation and acoustic analogy. The feather-like serrations are found to achieve the most prominent noise reduction among the four types of curved serrations, especially in the low-frequency range. With the aid of acoustic analogy, the coherence analysis of far-field noise produced by the dipole sources on the airfoil surface is performed. The results show that destructive interference is still the critical mechanism responsible for noise reduction. Considering only the dipole sources, we find that the feather-like serrated TE shape can obtain the best noise reduction performance among all the serrated cases. Furthermore, since flow structures are reorganized near the TE serrations, we investigated the flow noise sources quantitatively in the near field. In these cases, the noise source due to flow structures is suppressed to a greater extent in the feather-like serrated case near the TE serration roots. Consequently, the above findings indicate that the feather-like serration favors suppressing dipole and flow noise sources in the near field, which makes it an efficient configuration for reducing airfoil noise.

A high-order transition route from inertial to elasticity-dominated turbulence (EDT) in Taylor–Couette flows of polymeric solutions has been discovered via direct numerical simulations. This novel two-step transition route is realized by enhancing the extensional viscosity and hoop stresses of the polymeric solution via increasing the maximum chain extension at a fixed polymer concentration. Specifically, in the first step inertial turbulence is stabilized to a laminar flow much like the modulated wavy vortex flow. The second step destabilizes this laminar flow state to EDT, i.e. a spatially smooth and temporally random flow with a $-3.5$ scaling law of the energy spectrum reminiscent of elastic turbulence. The flow states involved are distinctly different to those observed in the reverse transition route from inertial turbulence via a relaminarization of the flow to elasto-inertial turbulence in parallel shear flows, underscoring the importance of polymer-induced hoop stresses in realizing EDT that are absent in parallel shear flows.

The flow physics of inertio-elastic turbulent Taylor–Couette flow for a radius ratio of $0.5$ in the Reynolds number ( $Re$ ) range of $500$ to $8000$ is investigated via direct numerical simulation. It is shown that as $Re$ is increased the turbulence dynamics can be subdivided into two distinct regimes: (i) a low $Re \leqslant 1000$ regime where the flow physics is essentially dominated by nonlinear elastic forces and the main contribution to transport and mixing of momentum, stress and energy comes from large-scale flow structures in the bulk region and (ii) a high $Re \geqslant 5000$ regime where inertial forces govern the flow physics and the flow dynamics is mainly governed by small-scale flow structures in the near-wall region. Flow–microstructure coupling analysis reveals that the elastic Görtler instability in the near-wall region is triggered via significant polymer extension and commensurately high hoop stresses. This instability gives rise to small-scale elastic vortical structures identified as elastic Görtler vortices which are present at all $Re$ considered. In fact, these vortices develop herringbone streaks near the inner wall that have a longer average life span than their Newtonian counterparts due to their elastic origin. Examination of the budgets of mean streamwise enstrophy, mean kinetic energy, turbulent kinetic energy and Reynolds shear stress demonstrates that increasing fluid inertia hinders the generation of elastic stresses, leading to a monotonic reduction of the elastic-related effects on the flow physics.

It is fascinating that fish groups spontaneously form different formations. The collective locomotions of two and multiple undulatory self-propelled foils swimming in a fluid are numerically studied and the deep reinforcement learning (DRL) is applied to control the locomotion. We explored whether typical patterns emerge spontaneously under the driven two DRL strategies. One strategy is that only the following fish gets hydrodynamic advantages. The other is that all individuals in the group take advantage of the interaction. In the DRL strategy, we use swimming efficiency as the reward function, and the visual information is included. We also investigated the effect of involving hydrodynamic force information, which is an analogy to that detected by the lateral line of fish. Each fish can adjust its undulatory phase to achieve the goal. Under the two strategies, collective patterns with different characteristics, i.e., the staggered-following, tandem-following phalanx and compact modes emerge. They are consistent with the results in the literature. The hydrodynamic mechanism of the above high-efficiency collective traveling modes is analyzed by the vortex-body interaction and thrust. We also found that the time sequence feature and hydrodynamic information in the DRL are essential to improve the performance of collective swimming. Our research can reasonably explain the controversial issue observed in the relevant experiments. The paper may be helpful for the design of bionic fish.

In human beings and experimental animals, maladaptive impulsivity is manifested by the acute injection of psychostimulants, such as amphetamine. Cannabinoid CB1 receptors have been implicated in the regulation of stimulant-induced impulsive action, but the role of CB1 receptors in timing-related impulsive action by amphetamine remains unknown. Male rats were used in evaluating the effects of CB1 receptor antagonist and agonist (SR141716A and WIN55,212–2, respectively) systemically administered individually and combined with d-amphetamine on a differential reinforcement of low-rate response (DRL) task, an operant behavioral test of timing and behavioral inhibition characterized as a type of timing impulsive action. A distinct pattern of DRL behavioral changes was produced by acute d-amphetamine (0, 0.5, 1.0, and 1.5 mg/kg) treatment in a dose-dependent fashion, whereas no significant dose effect was detected for acute SR141716A (0, 0.3, 1, and 3 mg/kg) or WIN55,212–2 (0, 0.5, 1, and 2 mg/kg) treatment. Furthermore, DRL behavior altered by 1.5 mg/kg d-amphetamine was reversed by a noneffective dose of SR141716A (3 mg/kg) pretreatment. The minimally influenced DRL behavior by 0.5 mg/kg d-amphetamine was affected by pretreatment with a noneffective dose of WIN55,212–2 (1 mg/kg). These findings reveal that the activation and blockade of CB1 receptors can differentially modulate the timing impulsive action of DRL behavior induced by acute amphetamine treatment. Characterizing how CB1 receptors modulate impulsive behavior will deepen our understanding of the cannabinoid psychopharmacology of impulsivity and may be helpful in developing an optimal pharmacotherapy for reducing maladaptive impulsivity in patients with some psychiatric disorders.

Structures of wall turbulence due to the mean shear created by the wall are generated. Numerical tests are performed with the rough-wall-like mean shear imposed in the near-wall region without resolving the surface roughness in the constrained large-eddy simulation. The results indicate that the major effects of roughness on wall turbulence can be well reproduced.

Pool boiling heat transfer enhancement on structured surfaces with columns is investigated by a thermal two-phase lattice Boltzmann model. The effects of geometrical parameters, including column height (H), width (W) and gap spacing (D), are discussed in detail. The size of columns is comparable to the diameter of a typical detached bubble ( d b ). It is found that the heat transfer performance mainly depends on two factors: the heated surface area A and local convective flow field. When W and D are properly chosen as about W = D ≈ 4 d b , increasing H enhances the heat transfer solely due to the increase of the surface area. When W and D are small, it is observed that the top of the columns and the channels are mostly covered by a layer of vapor, respectively, which weaken the convection. Under the circumstances, although surface area A increases, heat transfer enhancement is not so significant. For the surface tension effect, the enhancement of the structured surface decreases with the increase of capillary number compared to that of the plain surface. The mechanism is that at larger capillary number, the convection close to the channel may be partially blocked, which hinders the heat transfer. An optimal enhancement could be achieved when A is as large as possible and meanwhile the convection is not hindered by the bubbles behaviors. The finding here may shed some light on mechanism of heat transfer enhancement on structured surfaces with columns.

The topological evolution near the turbulent/non-turbulent interface (TNTI) in turbulent mixing layer is studied by means of statistical analysis of the invariants of velocity gradient tens...

Droplet evaporation on heterogeneous or patterned surfaces has numerous potential applications, for example, inkjet printing. The effect of surface heterogeneities on the evaporation of a nanometer-sized cylindrical droplet on a solid surface is studied using molecular dynamics simulations of Lennard-Jones particles. Different heterogeneities of the surface were achieved through alternating stripes of equal width but two chemical types, which lead to different contact angles. The evaporation induced by the heated substrate instead of the isothermal evaporation is investigated. It is found that the whole evaporation process is generally dominated by the nonuniform evaporation effect. However, at the initial moment, the volume expansion and local evaporation effects play important roles. From the nanoscale point of view, the slow movement of the contact line during the pinning process is observed, which is different from the macroscopic stationary pinning. Particularly, we found that the speed of the contact line may be not only affected by the intrinsic energy barrier between the two adjacent stripes ( ũ) but also relevant to the evaporation rate. Generally speaking, the larger the intrinsic energy barrier, the slower the movement of the contact line. At the specified temperature, when ũ is less than a critical energy barrier ( ũ*), the speed of the contact line would increase with the evaporate rate. When ũũ*, the speed of the contact line is determined only by ũ and no longer affected by the evaporation rate at different stages (the first stick and the second stick).

The effective viscosities of dilute and semi-dilute suspensions in a two-dimensional shear flow are studied using the lattice Boltzmann method. The suspensions contain non-Brownian hard circular buoyant porous particles. Here a more accurate formula for intrinsic viscosity as a function of Darcy number (Da) for the whole Da regime is proposed through our numerical result. The effects of fluid inertia, permeability of the particle, and confinement of the bounding walls are investigated. It is found that for the cases with a small Da, the effective viscosity significantly increases with confinement and fluid inertia. However, for the cases with a large Da, the confinement ratio and fluid inertia have very minor effect. Moreover, for semi-dilute suspensions, the permeability of the particle weakens the effect of the hydrodynamic interactions between particles on the relative viscosity ηr and makes ηr decrease. The above phenomena can be well understood through quantifying the disturbance of the porous particle to the flow.

Liquid films can be entrained when the dewetting velocity attains a threshold, and this dynamical wetting transition has been well studied in the situation of plane substrates. We investigate the forced dewetting in a capillary tube using diffuse-interface simulations and lubrication analysis, focusing on the onset of wetting transition and subsequent interface evolution. Results show that the meniscus remains stable when the displacing rate is below a threshold, beyond which film entrainment occurs and eventually leads to the formation of Taylor bubbles separated by liquid slugs, as has also been observed in the recent experiments of Zhao et al. (Phys. Rev. Lett., vol. 120, 2018, 084501). We derive an analytical solution of the critical capillary number, and demonstrate that the wetting transition is accompanied by a vanishing apparent contact angle and an abrupt drop of the contact-line velocity. Both the bubble and slug lengths are found to depend on the capillary number and the wettability of the wall. A theoretical formula for the bubble length is also proposed and compares favourably with numerical and experimental results.

In this work we extend the method of the constrained large-eddy simulation (CLES) to simulate the turbulent flow over inhomogeneous rough walls. In the original concept of CLES, the subgrid-scale (SGS) stress is constrained so that the mean part and the fluctuation part of the SGS stress can be modelled separately to improve the accuracy of the simulation result. Here in the simulation of the rough-wall flows, we propose to interpret the extra stress terms in the CLES formulation as the roughness-induced stress so that the roughness inhomogeneity can be incorporated by modifying the formulation of the constrained SGS stress. This is examined with the simulations of the channel flow with the spanwise alternating high/low roughness strips. Then the CLES method is employed to investigate the temporal response of the turbulence to the change of the wall condition from rough to smooth. We demonstrate that the temporal development of the internal boundary layer is just similar to that in a spatial rough-to-smooth transition process, and the spanwise roughness inhomogeneity has little impact on the transition process.

Direct numerical simulation of polymer-induced flow relaminarization of turbulent spanwise-rotating plane Couette flow (RPCF) is reported for the first time. Specifically, the reverse transition pathway from a Newtonian turbulent RPCF to a fully relaminarized drag enhanced viscoelastic flow has been elucidated. Evidently, this transition occurs gradually by weakening and eventual elimination of small-scale vortices as the Weissenberg number ( ) is enhanced, paving the way for a two-dimensional laminar flow consisting of large-scale and highly organized roll cells. The influence of polymer additives on convective momentum exchange by large-scale roll cells and small-scale turbulent vortices, namely, the drag reduction (DR) realized by elimination of turbulent vortices and the significant drag enhancement (DE) that results from polymer roll cell interactions has been identified as the mechanism of DE. The observed vortical changes point to a universal mechanism for the coupling of polymer chains and turbulent vortices in wall-bounded viscoelastic DE and DR flows.

We present an analytical model of nonlinear evolution of two-dimensional single-mode Rayleigh–Taylor instability (RTI) in cylindrical geometry at arbitrary Atwood number for the first time. Our model covers a full scenario of bubble evolution from the earlier exponential growth to the nonlinear regime with the bubbles growing in time as is formulated as a simplified function of the external acceleration, Atwood number and number of perturbation waves. This model's predictions are in good agreement with data from direct numerical simulations.

Intermittent locomotion is a widely used behavioral strategy for fish and birds to reduce the cost of movement. The intermittent locomotion performance of a self-propelled flapping plate is investigated numerically. Two intermittent swimming modes, namely, the multiple-tail-beat mode (MT mode) and the half-tail-beat mode (HT mode), as well as the continuous swimming mode (CT mode), are considered. Performance is evaluated from propulsive speed, efficiency, and cost of transport. The hydrodynamic performances of the intermittent modes are found to be better than the hydrodynamic performance of the CT mode when the bending stiffness $K$ is moderate [i.e., $K\ensuremath{\approx}O(1)]$ and the duty cycle is not too small. For the two intermittent modes, the performance of the HT mode is better than that of the MT mode when $K$ is small or moderate, while the situation is opposite when $K$ is large. It is found that compared to the asymmetric wake of the MT mode, the symmetric wake of the HT mode is favorable to generate more thrust force and therefore achieve better performance. Besides, at moderate $K$, the largest bending deformation of the plate in the HT mode, as well as the large normal force, produces the largest thrust during the flapping. The present results can help us to better understand the intermittent locomotion of animals and may be helpful for bionic design.

The subgrid effects on the dynamics of the filtered velocity gradient tensor (VGT) in compressible turbulence are studied by means of statistical analysis of the invariants of the filtered VGT in compressible mixing layers. The evolution of the filtered VGT is determined by the interaction among the invariants, the pressure effects, the viscous effects and the subgrid effects. Based on the probability fluxes in the plane of the second ( plane is identified in the locally expanded regions, which is determined by the non-normal effect. Further, an SGS model with the non-local effect is proposed to give a better prediction of the SGS energy dissipation.

Droplets interacting with deformable moving boundaries is ubiquitous. The flexible boundaries may dramatically affect the hydrodynamic behavior of droplets. A numerical method for simulating droplet impact on flexible substrates is developed. The effect of flexibility is investigated. To reduce the contact time and increase the remaining upward momentum in the flexible cases, the Weber number should be larger than a critical value. Moreover, the ratio of the natural frequency of the plate to that of the droplet ${F}_{r}$ should approximately equal to the reciprocal of the contact time of droplets impact on the rigid surfaces $({t}_{\text{ctr}})$ at the same We, e.g., ${F}_{r}\ensuremath{\approx}1/{t}_{\text{ctr}}$. Only under this circumstance would the kinetic energy convert into the surface energy of the droplet and the elastic energy of the plate simultaneously, and vice versa. Moreover, based on a double spring model, we proposed scaling laws for the maximal deflection of the plate and spreading diameter of the drop. Finally, the droplet impact under different wettability is qualitatively studied. We found that the flexibility may contribute to the droplet bouncing at a smaller contact angle.

The effect of trailing-edge shape on the self-propulsive performance of three-dimensional flexible plates is studied numerically. In our study, the trailing edges of the plates are symmetric chevron shapes, and the trailing-edge angle varies from (concave plate) to (convex plate). Under different bending stiffnesses , three regimes of the propulsive performance in terms of propulsive velocity and efficiency as a function of are identified. When is small, moderate and large, the square, convex and concave plate achieves the best performance, respectively. Analyses of vortical structures and velocity fields show that usually the jet behind the plate with the best performance is longest. Besides, the inclination angle of the jet may be small. The different propulsive performances at small and moderate are mainly attributed to the phase lag of the trailing edge. The force acting on the plate is analysed and it is found that the thrust force is mainly contributed by the normal force. If , and are rescaled by the normal force and the area moment of the plate, the curves for different almost collapse into a single curve when the bending stiffness coefficient is small or moderate. The scaling confirms that the normal force should be the characteristic fluid force at small or moderate and the effect is governed by the area moment. The findings may shed some light on the propulsive performance of aquatic animals.

Large-eddy simulation of sonic coaxial jets issuing into a quiescent environment was carried for two total pressure ratios of the inner to outer nozzle, i.e. 5/3 and 3/5. The effects of the total pressure ratio on the flow characteristics were mainly investigated. Various fundamental mechanisms dictating the complex flow characteristics, including jet shear layer evolution, shock system formation, shock/shear-layer interaction, turbulence behavior, and mixing property, have been studied. It is found that the total pressure ratio has an important influence on the flow evolution of coaxial jets as well as turbulence behavior and jet mixing property. The fluid-dynamic shearing and compressing processes are analyzed based on the Lamb vector curl and divergence. The multi-layer structures of the shear layers and shock waves are reasonably captured and the relevant fluid-dynamic processes are clearly clarified. It is also identified that the turbulence behavior and mixing property of the coaxial jets are mainly associated with the shearing effect in the outer shear layer region and the shearing and compressing coupled effect in the jet core region.

Energetic benefit and enhanced performance are considered among the most fascinating achievements of collective behaviours, e.g. fish schools and flying formations. The collective locomotion of two self-propelled flapping plates initially in a side-by-side arrangement is investigated numerically. Both in-phase and antiphase oscillations for the two plates are considered. It is found that the plates will spontaneously form some stable configurations as a result of the flow-mediated interaction, specifically, the staggered-following (SF) mode and the alternate-leading (AL) mode for the in-phase scenario and the moving abreast (MA) mode and the AL mode for the antiphase scenario. In the SF mode, the rear plate follows the front one with a staggered configuration. In the AL mode, the plates chase each other side-by-side alternately. In terms of propulsive speed and efficiency, the performance of the plates in the SF mode with small lateral spacing $H$ is found to be better than those in the tandem following case ($H=0$) and the side-by-side case (i.e. the AL mode). To achieve higher propulsive efficiency, no matter in-phase or antiphase oscillations, the two plates with moderate bending stiffness, e.g. $K\approx O(1)$, are preferred and they should be close enough in the lateral direction. For the side-by-side configuration, the performance of each plate in the antiphase and in-phase scenarios is enhanced and weakened in comparison with that of the isolated plate, respectively. Besides the pressure and vorticity contours, the normal force and thrust acting on the plates are also analysed. It is revealed that the thrust is mainly contributed by the normal force at moderate bending stiffness. The normal force and thrust are critical to the propulsive speed and efficiency. For two self-propelled plates, in view of hydrodynamics, to achieve higher performance the in-phase SF mode and antiphase flappings in the side-by-side configuration are preferred.

Drag reduction (DR) in plane Couette flow (PCF) induced by the addition of flexible polymers has been studied via direct numerical simulation (DNS). The similarities and differences in the drag reduction features of PCF and plane Poiseuille flow (PPF) have been examined in detail, particularly in regard to the polymer-induced modification of large-scale structures (LSSs) in the near-wall turbulence. Specifically, it has been demonstrated that in the near-wall region, drag-reduced PCF has features similar to those of drag-reduced PPF; however, in the core region, intriguing differences are found between these two drag-reduced shear flows. Chief among these differences is the significant polymer stretch that arises from the enhanced exchanges between elastic potential energy and turbulent kinetic energy and the commensurate observation of peak values of the conformation tensor components $\unicode[STIX]{x1D60A}_{yy}$ and $\unicode[STIX]{x1D60A}_{zz}$ in this region. This finding is in stark contrast to that of drag-reduced PPF where the polymer stretch and the exchanges between elastic potential energy and turbulent kinetic energy in the core region are insignificant; to this end, in drag-reduced PPF, peak values of the conformation tensor components appear in the near-wall region. Therefore, this study paves the way for understanding the underlying flow physics in drag-reduced PCF, particularly in the context of elastic theory. Moreover, the longitudinal large-scale streaks at the channel centre of drag-reduced PCF are greatly strengthened due to the increased production/dissipation ratio; the LSS imprint effects on the near-wall flow of drag-reduced PCF monotonically increase as the Weissenberg number is enhanced.

We numerically investigate the mechanism leading to the entrapment of spheres at the gas–liquid interface after impact. Upon impact onto a liquid pool, a hydrophobic sphere is seen to follow one of the three regimes identified in the experiment (Lee & Kim, Langmuir, vol. 24, 2008, pp. 142–145): sinking, bouncing or being entrapped at the interface. It is important to understand the role of wettability in this process of flow–structure interaction with dynamic wetting, and in particular, to what extent the wettability can determine whether the sphere is entrapped at the interface. For this purpose, a diffuse-interface immersed boundary method is adopted in the numerical simulations. We expand the parameter space considered previously, provide the phase diagrams and identify the key phenomena in the impact dynamics. Then, we propose the scaling models to interpret the critical conditions for the occurrence of sphere entrapment, accounting for the wettability of the sphere. The models are shown to provide a good correlation among the impact inertia of the drop, the surface tension, the wettability and the density ratio of the sphere to the liquid.

The interaction of tandem inverted flexible flags in a uniform flow is investigated. For the inverted flags, their ends are fixed with their heads freely flapping. A direct numerical simulation is performed for which the Reynolds number is of order 200. Large flapping amplitude as well as large drag force is preferred because more energy may be harvested if more bending energy is generated. For the simple case of two tandem inverted flags, the drag force and flapping amplitude of the rear flag are found to be smaller than those of an isolated inverted flag due to the destructive merging mode of vortices. However, it is still unknown whether more bending energy can be generated when coupled inverted flags are arranged properly. To explore the possibility, inverted flags are proposed to be arranged as two rows, which indicate two lines of inverted flags perpendicular to the direction of the incoming flow, and flags in the front and rear rows are in-line or staggered. First the results for infinite flags with periodic boundary condition are presented. In both the in-line and the staggered arrangements, due to the interactions between the front–rear flags, the flapping amplitude or the maximum bending deformation and bending energy of a flag in the rear row can be enhanced, which may be significantly higher than those of an isolated case. Meanwhile, the bending energy of a flag in the front row is close to that of an isolated case. Second, results for finite inverted flag groups show that antiphase synchronization is preferred. When the group number is large enough, the bending energies of the front and rear flags in the inner groups are close to those in the infinite case. This finding may be helpful for the designing of an efficient energy harvesting device using inverted flags.

The complete growth process of a single bubble in quiescent liquid is simulated using a three-dimensional hybrid thermal lattice Boltzmann model. The non-equilibrium extrapolation pressure boundary condition is extended to handle the thermal multiphase flow. Unfavorable spurious currents are usually generated in the vicinity of curved interfaces when two-phase lattice Boltzmann methods are applied. Here a level-set scheme is incorporated into the simulations to accurately represent interfacial dynamics. The phase change is controlled by an equation of state automatically instead of any artificial phase change model. Hence the present simulation is more accurate and thermodynamically consistent. The temperature, velocity fields during the bubble growth are consistent with relevant theories. The bubble growth rate obtained from the lattice Boltzmann simulations agree well with the analytical solutions. The result shows that the present scheme is able to simulate the relevant thermal bubble dynamics quantitatively.