We present evidence revealing that an object with specific properties can exhibit multiple stable falling postures at low Reynolds numbers. By scrutinizing the force equilibrium relationship of a fixed object at various attack angles and Reynolds numbers, we introduce a methodology that can obtain the stable falling postures of the object. This method saves computational resources and more intuitively presents the results in the full parameter domain. Our findings are substantiated by free-fall tests conducted through both physical experiments and numerical simulations, which validate the existence of multiple stable solutions in accordance with the interpolation results obtained with fixed objects. Additionally, we quantify the abundance and distribution patterns of stable falling postures for a diverse range of representative shapes. This discovery highlights the existence of multiple stable solutions that are universally present across objects of different shapes. The implications of this research extend to the design, stability control and trajectory prediction of all free and controlled flights in both air and water.

Wall temperature has a significant effect on shock wave/turbulent boundary layer interactions (STBLIs) and has become a non-negligible factor in the design process of hypersonic vehicles. In this paper, direct numerical simulations are conducted to investigate the wall temperature effects on STBLIs over a 34° compression ramp at Mach number 6. Three values of the wall-to-recovery-temperature ratio (0.50, 0.75 and 1.0) are considered in the simulations. The results show that the size of the separation bubble declines significantly as the wall temperature decreases. This is because the momentum profile of the boundary layer becomes fuller with wall cooling, which means the near-wall fluid has a greater momentum to suppress flow separation. An equation based on the free-interaction theory is proposed to predict the distributions of the wall pressure upstream of the corner at different wall temperatures. The prediction results are generally consistent with the simulation results (Reynolds number Reτ ranges from 160 to 675). In addition, the low-frequency unsteadiness is studied through the weighted power spectral density of the wall pressure and the correlation between the upstream and downstream. The results indicate that the low-frequency motion of the separation shock is mainly driven by the downstream mechanism and that wall cooling can significantly suppress the low-frequency unsteadiness, including the strength and streamwise range of the low-frequency motions.

The turbulent boundary layer (TBL) is a widely existing flow phenomenon in nature and engineering applications. Its strong mixing effect can achieve more sufficient material mixing, heat transport, etc. The understanding of the entrainment process and mechanism of irrotational fluids entering the turbulent region can be promoted by studying the geometric and dynamic characteristics of turbulent${/}$non-turbulent interfaces (TNTI). In compressible flow, it is unclear whether the properties of TNTI will change and whether the entrainment will show different features due to the influence of compressibility. Based on the direct numerical simulation results of supersonic compressible plate TBLs with Mach number of 2.9, the geometric and dynamic characteristics of TNTI are investigated in this paper. The interface is identified by the enstrophy method, and the height, thickness, fractal dimension, enstrophy transportation and entrainment characteristics of the interface are investigated. It is found that for the enstrophy transportation in a TBL, the contribution of compressibility-related terms accounts for approximately 13.4 % of the total enstrophy transportation, which tends to transfer the enstrophy of turbulence near the interface to both directions vertical to the interface. This promotes the expansion of the turbulent region towards the non-turbulent region, and the mean height, thickness and entrainment velocity are increased by approximately 3.7 %, 7.0 % and 8.5 %, respectively, while the fractal dimension is basically unaffected. Different from the incompressible flow, the contribution of the compressibility-related terms to the entrainment velocity is independent of the local curvature, and the intense entrainment process is more likely to occur on a highly curved concave surface.

In this paper, direct numerical simulations in hypersonic turbulent boundary layers over a $24^{\circ }$ compression ramp at Mach 6.0 are performed. The wall skin friction and its spanwise non-homogeneity in the interaction region are analysed via the spectral analysis and drag decomposition method. On the compression ramp, the premultiplied spanwise energy spectrum of wall shear stress $\tau _{w}$ reveals two energetic spanwise length scales. One occurs in the region of $x/\delta _{ref}=0\unicode{x2013}3$ ($x=0$ lies in the compression corner; $\delta _{ref}$ is the boundary layer thickness upstream of the interaction region) and is consistent with that of the large-scale streamwise vortices, indicating that the fluctuation intensity of $\tau _{w}$ is associated with the Görtler-type structures. The other one is observed downstream of $x/\delta _{ref}=3.0$ and corresponds to the regenerated elongated streaky structures. The fluctuation intensity of $\tau _{w}$ peaks at $x/\delta _{ref}=3.0$, where both the above energetic length scales are observed. The drag decomposition method proposed by Li et al. (J. Fluid Mech., vol. 875, 2019, pp. 101–123) is extended to include the effects of spanwise non-homogeneity so that it can be used in the interaction region where the mean flow field and the mean skin friction $C_f$ exhibit an obvious spanwise heterogeneity. The results reveal that, in the upstream turbulent boundary layer, the drag contribution arising from the spanwise heterogeneity can be neglected, while this value on the compression ramp is up to 20.7 % of $C_f$, resulting from the Görtler-type vortices. With the aid of the drag decomposition method, it is found that the main flow features that contribute positively to the amplification of $C_f$ and its rapid increase on the compression ramp includes: the density increase across the shock, the high mean shear stress and turbulence amplification around the detached shear layer and the Favre-averaged downward velocity towards the ramp wall. Compared with the spanwise-averaged value, $C_f$ and its components at the spanwise station where the downwash and upwash of the Görtler-type vortices occur reveal a spanwise variation exceeding 10 %.

Oocytes with excessively large first polar bodies (PB1) often occur in assisted reproductive procedures. Many times these oocytes are discarded without insemination and, as a result, the application of this portion of oocytes has scarcely been reported to date. Few studies have examined large PB1 oocytes in infertile women and have virtually entirely studied genetic variations for large PB1 oocyte abnormalities. Here, we describe an unusual case of a live birth from a remarkably large PB1 oocyte in a frozen embryo transfer (FET) cycle. This is the first instance of a successful live birth resulting from a PB1 oocyte with an extremely large polar body measuring 80 μM × 40 μM in size. The large PB1 oocyte was performed by an early rescue intracytoplasmic sperm injection (r-ICSI) and was formed into a blastocyst on day 5. Following FET, a healthy boy baby weighing 3100 g was finally delivered by caesarean section at 37 weeks and 5 days after conception. Additionally, there were no complications throughout the antenatal period or the perinatal phase of this following full-term delivery. In this study, it is revealed for the first time that a huge PB1 oocyte can be fertilized, resulting in the growth of a blastocyst, a subsequent pregnancy, and a live birth. This new information prompts us to reconsider the use of large PB1 oocytes. More insightful talks should be given attention to prevent the waste of embryos because not all oocytes with aberrant morphology are unavailable.

Helicity plays a key role in the evolution of vortex structures and turbulent dynamics. The helicity dynamics and vortex structures in streamwise-rotating channel turbulence are discussed in this paper using the helicity budget equation and the differentiated second-order structure function equation of helicity. Generally, rotation and Reynolds numbers exhibit opposing effects on the interscale helicity dynamics and the vortices. Under the buffer layer, the positions of the helicity peaks are proportional to the ratio between the Reynolds and rotation numbers. The mechanism is related to the opposing effects of convection and rotation. Rotation directly affects the helicity balance through the Coriolis term and corresponding pressure term. In the buffer layer, the scale helicity is negative at small scales but positive at large scales, which is mainly induced by the spatial effects (the production and the spatial turbulent convection) but reduced by interscale cascades. Examination of structures reveals the close association between scale helicity and streaks, with streak lift angles exhibiting an increase with rotation and a decrease with Reynolds numbers. In the log-law layer, the Coriolis terms and corresponding pressure terms are proportional to the rotation numbers but remain independent of the Reynolds numbers. The negative scale helicity is forward cascaded towards small scales. Generally, spanwise vortices in the log-law layer are related to sweep events and forward cascades. Our findings indicate that these spanwise vortices are suppressed by rotation but recover with increasing Reynolds numbers, aligning with the effects observed in the scale helicity balance.

In this paper the multiscale dynamics of streamwise-rotating channel turbulence is studied through direct numerical simulations. Using the generalized Kolmogorov equation, we find that as rotation becomes stronger, the turbulence in the buffer layer is obviously reduced by the intense spatial turbulent convection. On the contrary, in other layers, the turbulence is strengthened mainly by the modified production peak, the intense spatial turbulent convection and the suppressed forward energy cascades. It is also discovered that under a system rotation, small- and large-scale inclined structures have different angles with the streamwise direction, and the difference is strengthened with increasing rotation rates. The multiscale inclined structures are further confirmed quantitatively through a newly defined angle based on the velocity vector. Through the budget balance of Reynolds stresses and the hairpin vortex model, it is discovered that the Coriolis force and the pressure–velocity correlation are responsible for sustaining the inclined structures. The Coriolis force directly decreases the inclination angles but indirectly induces inclined structures in a more predominant way. The pressure–velocity correlation term is related to the strain rate tensor. Finally, the anisotropic generalized Kolmogorov equation is used to validate the above findings and reveals that the multiscale behaviours of the inclined structures are mainly induced by the mean spanwise velocity gradients.

The development of high-intensity ultrafast laser facilities provides the possibility to create novel physical phenomena and matter states. The timing fluctuation of the laser pulses is crucial for pump–probe experiments, which is one of the vital means to observe the ultrafast dynamics driven by intense laser pulses. In this paper, we demonstrate the timing fluctuation characterization and control of the front end of a 100-PW laser that is composed of a high-contrast optical parametric amplifier (seed) and a 200-TW optical parametric chirped pulse amplifier (preamplifier). By combining the timing jitter measurement with a feedback system, the laser seed and preamplifier are synchronized to the reference with timing fluctuations of 1.82 and 4.48 fs, respectively. The timing system will be a key prerequisite for the stable operation of 100-PW laser facilities and provide the basis for potential pump–probe experiments performed on the laser.

We demonstrate an ultra-broadband high temporal contrast infrared laser source based on cascaded optical parametric amplification, hollow-core fiber (HCF) and second harmonic generation processes. In this setup, the spectrum of an approximately 1.8 μm laser pulse has near 1 μm full bandwidth by employing an argon gas-filled HCF. Subsequently, after frequency doubling with cascaded crystals and dispersion compensation by a fused silica wedge pair, 9.6 fs (~3 cycles) and 150 μJ pulses centered at 910 nm with full bandwidth of over 300 nm can be generated. The energy stability of the output laser pulse is excellent with 0.8% (root mean square) over 20 min, and the temporal contrast is >1012 at –10 ps before the main pulse. The excellent temporal and spatial characteristics and stability make this laser able to be used as a good seed source for ultra-intense and ultrafast laser systems.

Direct numerical simulations are performed to investigate the wake transitions of the flow normal to a circular rotating disk. The diameter-thickness aspect ratio of the disk is $\chi =50$. The Reynolds number of the free stream is defined as $Re_s=U_\infty D/\nu$, with incoming flow velocity $U_\infty$, disk diameter $D$, and kinematic viscosity of the fluid $\nu$. The rotational motion of the disk is described by the Reynolds number of rotation $Re_r=\varOmega Re_s$, with non-dimensional rotation rate $\varOmega =\frac {1}{2}\omega D/U_\infty$, where $\omega$ is the angular rotation speed of the disk. Extensive numerical simulations are performed in the parameter space $50 \leqslant Re_s \leqslant 250$ and $0 \leqslant Re_r \leqslant 250$, in which six flow regimes are identified as follows: the axisymmetric state, the low-speed steady rotation (LSR) state, the high-speed steady rotation (HSR) state, the low-speed unsteady rotation (LUR) state, the rotational vortex shedding state, and the chaotic state. Although plane symmetry exists in the wake when the disk is stationary, a small rotation will immediately destroy its symmetry. However, the vortex shedding frequencies and wake patterns of the stationary disk are inherited by the unsteady rotating cases at low $Re_r$. A flow rotation rate jump is observed at $Re_s\approx 125$. The LUR state is intermediate between the LSR and HSR states. Due to the rotational motion, the wake of the disk enters the steady rotation state earlier at large $Re_r$, and is delayed into the vortex shedding state in the whole range of $Re_r$. In the steady rotation states (LSR and HSR), the steady flow rotation rate is linearly correlated with the disk rotation rate. It is found that the rotation of the disk can restrain the vortex shedding. The chaotic state can be regularized by the medium rotation speed of the disk.

We report dispersion management based on a mismatched-grating compressor for a 100 PW level laser, which utilizes optical parametric chirped pulse amplification and also features large chirped pulse duration and an ultra-broadband spectrum. The numerical calculation indicates that amplified pulses with 4 ns chirped pulse duration and 210 nm spectral bandwidth can be directly compressed to sub-13 fs, which is close to the Fourier-transform limit (FTL). More importantly, the tolerances of the mismatched-grating compressor to the misalignment of the stretcher, the error of the desired grating groove density and the variation of material dispersion are comprehensively analyzed, which is crucially important for its practical application. The results demonstrate that good tolerances and near-FTL compressed pulses can be achieved simultaneously, just by keeping a balance between the residual second-, third- and fourth-order dispersions in the laser system. This work can offer a meaningful guideline for the design and construction of 100 PW level lasers.

A quasi-dynamic subgrid-scale (SGS) kinetic energy one-equation eddy-viscosity model is introduced in this paper for large-eddy simulation (LES) of compressible flows. With the additional SGS kinetic energy equation, the SGS kinetic energy can be predicted properly. Then, using the dual constraints of SGS kinetic energy and the SGS kinetic energy flux, the eddy-viscosity model can be determined exactly. Taking a similar scheme as the expansion of the SGS stress, other unclosed quantities in the equations to be solved could be well modelled separately. Therefore, with the advance of the equations, all the model coefficients can be determined dynamically. Differing from the classic dynamic procedure, the new methodology needs no test filtering, and thus it could also be called a quasi-dynamic procedure. Using direct numerical simulation of compressible turbulent channel flow, the $a$ $priori$ test shows that the key modelled quantities of the suggested model display high correlations with the real values. In LES of compressible turbulent channel flows of the Mach number being $1.5$ and 3.0, the proposed model can precisely predict some important quantities, including the mean velocity, Reynolds stress and turbulent flux, and it can also supply more abundant turbulent structures. For the compressible flat-plate boundary layer, the new model can correctly predict the transition process, mean velocity and turbulence intensities in the turbulent region. The results show that the proposed model has the advantage of scale adaptivity. Finally, the new model is applied to LES of turbulent mixing in spherical converging Richtmyer–Meshkov instability, and the accurate results show that the new model has a good ability for LES of complex fluids.

Helical rotating turbulence is a chiral and anisotropic flow. The energy and helicity transfers of helical rotating turbulence are studied in this paper. First, we discuss the antisymmetry and conservation of energy and helicity transfers. There are three expressions for helicity transfers due to the commutability of differential operators. The first expression is derived here. The second expression violates antisymmetry, and the third one introduces non-physical effects. The relations of these expressions are discussed in detail, including those about the sum of all triads and partial triads, as well as those about helical wave decomposition. Through direct numerical simulations, we find that helicity can reduce inverse energy cascades. The inhibition is mainly associated with transhelical energy fluxes and the interactions of two-dimensional modes. The inverse cascades of decomposed energy fluxes are related to the two-dimensionalization. For helicity, rotation does not affect the total helicity flux but generally suppresses the decomposed helicity fluxes. Positive homochiral and negative heterochiral helicity fluxes are associated with corresponding positive anisotropic transfers. Notably, the transhelical helicity fluxes increase the amplitudes of both positive and negative helicity, which is related to the chirality polarization.

Direct numerical simulations of a three-dimensional cylindrical Richtmyer–Meshkov instability with and without chemical reactions are carried out to explore the chemical reaction effects on the statistical characteristics of transition and turbulent mixing. We adopt 9-species and 19-reaction models of non-premixed hydrogen and oxygen separated by a multimode perturbed cylindrical interface. A new definition of mixing width suitable for a chemical reaction is introduced, and we investigate the spatio-temporal evolution of typical flow parameters within the mixing regions. After reshock with a fuller mixing of fuels and oxygen, the chemical reaction becomes sufficiently apparent at affecting the evolution of the flow fields. Because of the generation of a combustion wave within the combustion regions and propagation, the growth of the mixing width with a chemical reaction is accelerated, especially around the outer radius with large temperature gradient profiles. However, the viscous dissipation rate in the early stage of the chemical reaction is greater because of heat release, which results in weakened turbulent mixing within the mixing regions. We confirm that small-scale structures begin to develop after reshock and then decay over time. During the developing process, helicity also begins to develop, in addition to kinetic energy, viscous dissipation rate, enstrophy, etc. In the present numerical simulations with cylindrical geometry, the fluctuating flow fields evolve from quasi-two-dimensional perturbations, and the generations of helicity can capture this transition process. The weakened fluctuations during shock compression can be explained as the inverse energy cascade, and the chemical reaction can promote this inverse energy cascade process.

Helicity is a quadratic inviscid conservative quantity in three-dimensional turbulent flows and is crucial for turbulent system evolution. Helicity effects have mainly been highlighted over the past few decades to explore the intrinsic mechanism of turbulent flows, while the statistical characteristics of helicity itself are nearly absent in general anisotropic turbulent flows. In this paper, we investigate the helicity statistics in turbulent channel flows with streamwise rotation at moderate rotation numbers ($Ro_{\tau }=7.5,15$ and 30) and Reynolds numbers ($Re_{\tau }=180$ and 395), including their spatial and scale distributions, anisotropy and cross-scale transfer. The appearance of a mean secondary flow in the spanwise direction corresponds to a mean streamwise vorticity, which indicates the presence of a high-helicity distribution. Numerical results reveal a regular helicity profile along the wall-normal direction, and a new peak is found in the near-wall region around $y^+=6$ of the streamwise or spanwise helicity profiles. The inter-scale helicity transfer is analysed by the filtering method, and the numerical consequences reveal that the second channel of the helicity cascade we proposed previously is dominant in contrast to the first channel. The rotation effects are explored by comparing the numerical results obtained under different rotation numbers. With increasing rotation number, more helical structures in the near-wall regions are present, with peaks of helicity profiles and fluxes coming closer to the wall. With a higher Reynolds number, their amplitudes are larger and scale-space transfer is strengthened. These systematic numerical analyses uncover the helicity distributions and transfer in wall-bounded turbulent flows.

Kinetic energy flux (KEF) is an important physical quantity that characterizes cascades of kinetic energy in turbulent flows. In large-eddy simulation (LES), it is crucial for the subgrid-scale (SGS) model to accurately predict the KEF in turbulence. In this paper, we propose a new eddy-viscosity SGS model constrained by the properly modelled KEF for LES of compressible wall-bounded turbulence. The new methodology has the advantages of both accurate prediction of the KEF and strong numerical stability in LES. We can obtain an approximate KEF by the tensor-diffusivity model, which has a high correlation with the real value. Then, using the artificial neural network method, the local ratios between the real KEF and the approximate KEF are accurately modelled. Consequently, the SGS model can be improved by the product of that ratio and the approximate KEF. In LES of compressible turbulent channel flow, the new model can accurately predict mean velocity profile, turbulence intensities, Reynolds stress, temperature–velocity correlation, etc. Additionally, for the case of a compressible flat-plate boundary layer, the new model can accurately predict some key quantities, including the onset of transitions and transition peaks, the skin-friction coefficient, the mean velocity in the turbulence region, etc., and it can also predict the energy backscatters in turbulence. Furthermore, the proposed model also shows more advantages for coarser grids.

The effect of wall temperature on the transfer of kinetic energy in a hypersonic turbulent boundary layer for different Mach numbers and wall temperature ratios is studied by direct numerical simulation. A cold wall temperature can enhance the compressibility effect in the near-wall region through increasing the temperature gradient and wall heat flux. It is shown that the cold wall temperature enhances the local reverse transfer of kinetic energy from small scales to large scales, and suppresses the local direct transfer of kinetic energy from large scales to small scales. The average filtered spatial convection and average filtered viscous dissipation are dominant in the near-wall region, while the average subgrid-scale flux of kinetic energy achieves its peak value in the buffer layer. It is found that the wall can suppress the inter-scale transfer of kinetic energy, especially for the situation of a cold wall. A strong local reverse transfer of fluctuating kinetic energy is identified in the buffer layer in the inertial range. Helmholtz decomposition is applied to analyse the compressibility effect on the subgrid-scale flux of kinetic energy. A strong transfer of the solenoidal component of fluctuating kinetic energy is identified in the buffer layer, while a significant transfer of the dilatational component of fluctuating kinetic energy is observed in the near-wall region. It is also shown that compression motions have a major contribution to the direct transfer of fluctuating kinetic energy, while expansion motions play a marked role in the reverse transfer of fluctuating kinetic energy.

In this paper, the Richtmyer–Meshkov instabilities in spherical and cylindrical converging geometries with a Mach number of approximately 1.5 are investigated by using the high resolution implicit large eddy simulation method, and the influence of the geometric effect on the turbulent mixing is investigated. The heavy fluid is sulphur hexafluoride (SF6), and the light fluid is nitrogen (N2). The shock wave converges from the heavy fluid into the light fluid. The Atwood number is 0.678. The total structured and uniform Cartesian grid node number in the main computational domain is 20483. In addition, to avoid the influence of boundary reflection, a sufficiently long sponge layer with 50 non-uniform coarse grids is added for each non-periodic boundary. Present numerical simulations have high and nonlinear initial perturbation levels, which rapidly lead to turbulent mixing in the mixing layers. Firstly, some physical-variable mean profiles, including mass fraction, Taylor Reynolds number, turbulent kinetic energy, enstrophy and helicity, are provided. Second, the mixing characteristics in the spherical and cylindrical turbulent mixing layers are investigated, such as molecular mixing fraction, efficiency Atwood number, turbulent mass-flux velocity and density self-correlation. Then, Reynolds stress and anisotropy are also investigated. Finally, the radial velocity, velocity divergence and enstrophy in the spherical and cylindrical turbulent mixing layers are studied using the method of conditional statistical analysis. Present numerical results show that the geometric effect has a great influence on the converging Richtmyer–Meshkov instability mixing layers.

Helicity, as one of only two inviscid invariants in three-dimensional turbulence, plays an important role in the generation and evolution of turbulent flows. Through theoretical analyses, we find that there are two channels in the helicity cascade process, which differs dramatically from the traditional viewpoint. In this paper, we have conducted important research on the newly proposed dual-channel helicity cascade theory, including vortex dynamic processes, intermittent discrepancies, tensor geometries, etc. The first channel mainly originates from the vortex twisting process, and the second channel mainly originates from the vortex stretching process. Antisymmetric tensors are introduced to the derivations of dual-channel helicity cascade theory, and a complex rotation frame leads to a higher helicity transfer efficiency. By analysing data from direct numerical simulations of typical turbulent flows, we find that these two channels behave differently. The ensemble averages of helicity flux in different channels are equal in homogeneous and isotropic turbulence, while they are different in other types of turbulent flows. The intermittency of the second channel is stronger than that of the first channel. In addition, we find a novel mechanism of hindered or even inverse energy cascades, which could be attributed to the second-channel helicity flux.

Water-based polyurethane/alumina hollow microsphere (WPU-hAl2O3) composite films were prepared via a facile spin coating method. The pristine WPU, as the matrix of the composite films, was tailor-made by hAl2O3 with the diameter of 2–5 μm to improve the mechanical and physical properties of the films. The hardness, surface morphology, infrared emissivity, wettability, and light transmittance of the WPU-hAl2O3 films with different hAl2O3 contents were investigated. The results indicate that the Vickers hardness, coefficient of friction, infrared emissivity at the wavelength of 2–22 μm, and wetting angle of the WPU-hAl2O3 films (30 wt%) increased by 53.6%, 51.7%, 21.1%, and 19.0%, respectively, compared with the pristine WPU films. Meanwhile, with the rising of hAl2O3 content, the light transmittance decreased by 75.3% at the wavelength of 400–800 nm. This work not only designs a kind of lightweight multifunctional composite film but also provides an effective route for extending further applications of hAl2O3 in the field of composite films.