The scaling relations mapping the turbulence statistics in compressible turbulent boundary layers (TBLs) onto their incompressible counterparts are of fundamental significance for turbulence modelling, such as the Morkovin scaling for velocity fields, while for pressure fluctuation fields, a corresponding scaling relation is currently absent. In this work, the underlying scaling relations of pressure fluctuations about Mach number ($M$) contained in their generation mechanisms are explored by analysing a series of direct numerical simulation data of compressible TBLs over a wide Mach number range $(0.5\leq M \leq 8.0)$. Based on the governing equation of pressure fluctuations, they are decomposed into components according to the properties of source terms. It is notable that the intensity of the compressible component, predominantly originating from the acoustic mode, obeys a monotonic distribution about the Mach number and wall distance; further, the intensity of the rest of the pressure components, which are mainly generated by the vorticity mode, demonstrates a uniform distribution consistent with its incompressible counterpart. Moreover, the coupling between the two components is negligibly weak. Based on the scaling relations, semiempirical models for the fluctuation intensity of both pressure and its components are constructed. Hence, a mapping relation is obtained that the profiles of pressure fluctuation intensities in compressible TBLs can be mapped onto their incompressible counterparts by removing the contribution from the acoustic mode, which can be provided by the model. The intrinsic scaling relation can provide some basic insight for pressure fluctuation modelling.

The efficacy of steady large-amplitude blowing/suction on instability and transition control for a hypersonic flat plate boundary layer with Mach number 5.86 is investigated systematically. The influence of the blowing/suction flux and amplitude on instability is examined through direct numerical simulation and resolvent analysis. When a relatively small flux is used, the two-dimensional instability critical frequency that distinguishes the promotion/suppression mode effect closely aligns with the synchronisation frequency. For the oblique wave, as the spanwise wavenumber increases, the suppression effects would become weaker and the mode suppression bandwidth diminishes/increases in general in the blowing/suction control. Increasing the blowing/suction flux can effectively broaden the frequency bandwidth of disturbance suppression. The influence of amplitude on disturbance suppression is weak in a scenario of constant flux. To gain a deeper insight into disturbance suppression mechanism, momentum potential theory (MPT) and kinetic energy budget analysis are further employed in analysing disturbance evolution with and without control. When the disturbance is suppressed, the blowing induces the transport of certain acoustic components along the compression wave out of the boundary layer, whereas the suction does not. The velocity fluctuations are derived from the momentum fluctuations of the MPT. Compared with the momentum fluctuations, the evolutions indicated by each component's velocity fluctuations greatly facilitate the investigations of the acoustic nature of the second mode. The rapid variation of disturbance amplitude near the blowing is caused by the oscillations of the acoustic component and phase speed differences between vortical and thermal components. Kinetic energy budget analysis is performed to address the non-parallel effect of the boundary layer introduced by blowing/suction, which tends to suppress disturbances near the blowing. Moreover, viscous effects leading to energy dissipation are identified to be stronger in regions where the boundary layer is rapidly thickening. Finally, it is demonstrated that a flat plate boundary layer transition triggered by a random disturbance can be delayed by a blowing/suction combination control. The resolvent analysis further demonstrates that disturbances with frequencies that dominate the early transition stage are dampened in the controlled base flow.

Achieving optimal nutritional status in patients with penetrating Crohn’s disease is crucial in preparing for surgical resection. However, there is a dearth of literature comparing the efficacy of total parenteral nutrition (TPN) v. exclusive enteral nutrition (EEN) in optimising postoperative outcomes. Hence, we conducted a case-matched study to assess the impact of preoperative EEN v. TPN on the incidence of postoperative adverse outcomes, encompassing overall postoperative morbidity and stoma formation, among penetrating Crohn’s disease patients undergoing bowel surgery. From 1 December 2012 to 1 December 2021, a retrospective study was conducted at a tertiary centre to enrol consecutive patients with penetrating Crohn’s disease who underwent surgical resection. Propensity score matching was utilised to compare the incidence of postoperative adverse outcomes. Furthermore, univariate and multivariate logistic regression analyses were conducted to identify the risk factors associated with adverse outcomes. The study included 510 patients meeting the criteria. Among them, 101 patients in the TPN group showed significant improvements in laboratory indicators at the time of surgery compared with pre-optimisation levels. After matching, TPN increased the occurrence of postoperative adverse outcomes (92·2 % v. 64·1 %, P = 0·001) when compared with the EEN group. In the multivariate analysis, TPN showed a significantly higher OR for adverse outcomes than EEN (OR = 4·241; 95 % CI 1·567–11·478; P = 0·004). The study revealed that penetrating Crohn’s disease patients who were able to fulfil their nutritional requirements through EEN exhibited superior nutritional and surgical outcomes in comparison with those who received TPN.

This study presents a comprehensive analysis on the extreme positive and negative events of wall shear stress and heat flux fluctuations in compressible turbulent boundary layers (TBLs) solved by direct numerical simulations. To examine the compressibility effects, we focus on the extreme events in two representative cases, i.e. a supersonic TBL of Mach number $M=2$ and a hypersonic TBL of $M=8$, by scrutinizing the coherent structures and their correlated dynamics based on conditional analysis. As characterized by the spatial distribution of wall shear stress and heat flux, the extreme events are indicated to be closely related to the structural organization of wall streaks, in addition to the occurrence of the alternating positive and negative structures (APNSs) in the hypersonic TBL. These two types of coherent structures are strikingly different, namely the nature of wall streaks and APNSs are shown to be related to the solenoidal and dilatational fluid motions, respectively. Quantitative analysis using a volumetric conditional average is performed to identify and extract the coherent structures that directly account for the extreme events. It is found that in the supersonic TBL, the essential ingredients of the conditional field are hairpin-like vortices, whose combinations can induce wall streaks, whereas in the hypersonic TBL, the essential ingredients become hairpin-like vortices as well as near-wall APNSs. To quantify the momentum and energy transport mechanisms underlying the extreme events, we proposed a novel decomposition method for extreme skin friction and heat flux, based on the integral identities of conditionally averaged governing equations. Taking advantage of this decomposition method, the dominant transport mechanisms of the hairpin-like vortices and APNSs are revealed. Specifically, the momentum and energy transports undertaken by the hairpin-like vortices are attributed to multiple comparable mechanisms, whereas those by the APNSs are convection dominated. In that, the dominant transport mechanisms in extreme events between the supersonic and hypersonic TBLs are indicated to be totally different.

Polydatin is an active polyphenol displaying multifaceted benefits. Recently, growing studies have noticed its potential therapeutic effects on bone and joint disorders (BJDs). Therefore, this article reviews recent in vivo and in vitro progress on the protective role of polydatin against BJDs. An insight into the underlying mechanisms is also presented. It was found that polydatin could promote osteogenesis in vitro, and symptom improvements have been disclosed with animal models of osteoporosis, osteosarcoma, osteoarthritis and rheumatic arthritis. These beneficial effects obtained in laboratory could be mainly attributed to the bone metabolism-regulating, anti-inflammatory, antioxidative, apoptosis-regulating and autophagy-regulating functions of polydatin. However, studies on human subjects with BJDs that can lead to early identification of the clinical efficacy and adverse effects of polydatin have not been reported yet. Accordingly, this review serves as a starting point for pursuing clinical trials. Additionally, future emphasis should also be devoted to the low bioavailability and prompt metabolism nature of polydatin. In summary, well-designed clinical trials of polydatin in patients with BJD are in demand, and its pharmacokinetic nature must be taken into account.

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.

The purpose of this retrospective study was to optimize the transplantation strategy for women of advanced maternal age to achieve live births within the shortest time. Data were collected from patients older than 40 years who underwent assisted reproductive therapy at our centre from 1 January 2009 to 31 December 2019. In total, 1023 cases of fresh cleavage embryo transfer (CET) cycles, 280 cases of frozen–thawed blastocyst transfer (FBT) cycles, and 26 cases of frozen–thawed CET (FCET) cycles were included. The main outcome was the live birth rate (LBR). The secondary outcomes were the clinical pregnancy rate (CPR) and neonatal outcomes. Multivariable logistic regression was performed to adjust for confounding factors. The blastocyst formation rate of patients older than 40 years was 23.5%, the freezing cycle rate was 19.8%, and the fresh-embryo transfer rate was 83.0%. The implantation rate, CPR, and LBR were significantly different among the CET, FCET, and FBT groups. There were no significant differences in multiple pregnancies and abortion rates among the groups, and neonatal outcomes were similar. Multivariate logistic regression analysis showed that, compared with the CET group, LBR did not increase in the FCET group, whereas LBR increased in the FBT group. For patients older than 40 years when having approximately eight embryos after fertilization, blastocyst transfer can be considered after fully discussing the advantages and disadvantages of blastocyst culture. Alternatively, CET can be performed first, followed by FBT if the cleavage embryo transfer is unsuccessful.

To address the shortcomings of existing methods for rotorcraft searching, positioning, tracking and landing on a ship at sea, a dual-channel LIDAR searching, positioning, tracking and landing system (DCLSPTLS) is proposed in this paper, which utilises the multi-pulse laser echoes accumulation method and the physical phenomenon that the laser reflectivity of the ship deck in the near-infrared band is four orders of magnitude higher than that of the sea surface. The DCLSPTLS searching and positioning model, tracking model and landing model are established, respectively. The searching and positioning model can provide estimates of the azimuth angle, the distance of the ship relative to the rotorcraft and the ship's course. With the above parameters as inputs, the total tracking time and the direction of the rotorcraft tracking speed can be obtained by using the tracking model. The landing model can calculate the pitch and the roll angles of the ship's deck relative to the rotorcraft by using the least squares method and the laser irradiation coordinates. The simulation shows that the DCLSPTLS can realise the functions of rotorcraft searching, positioning, tracking and landing by using the above parameters. To verify the effectiveness of the DCLSPTLS, a functional test is performed using a rotorcraft and a model ship on a lake. The test results are consistent with the results of the simulation.

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.

It is generally accepted that high-oleic crops have at least 70% oleate. As compared to their normal-oleic counterparts, oil and food products made from high-oleic peanut have better keeping quality and are much healthier. Therefore, high-oleic peanut is well recognized by processors and consumers. However, owing to the limited availability of high-oleic donors, most present-day high-oleic peanut varietal releases merely have F435 type FAD2 mutations. Through screening of a mutagenized peanut population of 15L46, a high-yielding peanut line with desirable elliptical oblong large seeds, using near infrared model for predicting oleate content in individual single seeds, high-oleic peanut mutants were identified. Sequencing FAD2A and FAD2B of the mutants along with the wild type revealed that these mutants possessed G448A FAD2A (F435 type FAD2A mutation) and G558A FAD2B (non-F435 type FAD2B mutation). Expression of the wild and mutated type FAD2B in yeast verified that the functional mutation contributed to the high-oleic phenotype in these mutants. The mutants provided additional high-oleic donors to peanut quality improvement.

Based on the fully compressible Navier–Stokes equations, the linear stability of thermal convection in rapidly rotating spherical shells of various radius ratios $\eta$ is studied for a wide range of Taylor number $Ta$, Prandtl number $Pr$ and the number of density scale height $N_\rho$. Besides the classical inertial mode and columnar mode, which are widely studied by the Boussinesq approximation and anelastic approximation, the quasi-geostrophic compressible mode is also identified in a wide range of $N_\rho$ and $Pr$ for all $\eta$ considered, and this mode mainly occurs in the convection with relatively small $Pr$ and large $N_\rho$. The instability processes are classified into five categories. In general, for the specified wavenumber $m$, the parameter space ($Pr, N_\rho$) of the fifth category, in which the base state loses stability via the quasi-geostrophic compressible mode and remains unstable, shrinks as $\eta$ increases. The asymptotic scaling behaviours of the critical Rayleigh numbers $Ra_c$ and corresponding wavenumbers $m_c$ to $Ta$ are found at different $\eta$ for the same instability mode. As $\eta$ increases, the flow stability is strengthened. Furthermore, the linearized perturbation equations and Reynolds–Orr equation are employed to quantitatively analyse the mechanical mechanisms and flow instability mechanisms of different modes. In the quasi-geostrophic compressible mode, the time-derivative term of disturbance density in the continuity equation and the diffusion term of disturbance temperature in the energy equation are found to be critical, while in the columnar and inertial modes, they can generally be ignored. Because the time-derivative term of the disturbance density in the continuity equation cannot be ignored, the anelastic approximation fails to capture the instability mode in the small-$Pr$ and large-$N_\rho$ system, where convection onset is dominated by the quasi-geostrophic compressible mode. However, all the modes are primarily governed by the balance between the Coriolis force and the pressure gradient, based on the momentum equation. Physically, the most important difference between the quasi-geostrophic compressible mode and the columnar mode is the role played by the disturbance pressure. The disturbance pressure performs negative work for the former mode, which appears to stabilize the flow, while it destabilizes the flow for the latter mode. As $\eta$ increases, in the former mode the relative work performed by the disturbance pressure increases and in the latter mode decreases.

The staining procedure is critical for investigating intra- and extra-cellular ultrastructure of microorganisms by transmission electron microscopy (TEM). Here, we propose a new ultra-low lead staining (ULLS) technique for preparing the ultrathin sections for TEM analysis. Sections of Enterobacter sp. (bacteria), Aspergillus niger (filamentous fungi), Rhodotorula mucilaginosa (fungi), and Chlamydomonas reinhardtii (microalgae) were tested. Compared with the sections prepared by the typical double-staining technique, ULLS-based sections showed evident advantages: (i) the staining process only required the addition of Pb(NO3)2; (ii) the Pb level during incubation was set as low as 1 mg/L, which had negligible toxicity to most microbial cells; (iii) the Pb cations were added during microbial culture, which avoided complicated sample preparation as in typical double staining. Taking C. reinhardtii as an example, the ULLS technique allowed fine investigation of microbial ultrastructure, e.g., starch granule, mitochondrion, Golgi apparatus, vacuole, and vesicle. Meanwhile, the physiological processes of the cells such as cell lysis and exocytosis were successfully captured, with relatively high contrast. This study hence shows a bright future on preparation of the high-quality ultrathin sections of microbial cells by the ULLS technique.

GBF1 [Golgi brefeldin A (BFA) resistance factor 1] is a member of the guanine nucleotide exchange factors Arf family. GBF1 localizes at the cis-Golgi and endoplasmic reticulum (ER)-Golgi intermediate compartment where it participates in ER-Golgi traffic by assisting in the recruitment of the coat protein COPI. However, the roles of GBF1 in oocyte meiotic maturation are still unknown. In the present study, we investigated the regulatory functions of GBF1 in mouse oocyte organelle dynamics. In our results, GBF1 was stably expressed during oocyte maturation, and GBF1 localized at the spindle periphery during metaphase I. Inhibiting GBF1 activity led to aberrant accumulation of the Golgi apparatus around the spindle. This may be due to the effects of GBF1 on the localization of GM130, as GBF1 co-localized with GM130 and inhibiting GBF1 induced condensation of GM130. Moreover, the loss of GBF1 activity affected the ER distribution and induced ER stress, as shown by increased GRP78 expression. Mitochondrial localization and functions were affected, as the mitochondrial membrane potential was altered. Taken together, these results suggest that GBF1 has wide-ranging effects on the distribution and functions of Golgi apparatus, ER, and mitochondria as well as normal polar body formation in mouse oocytes.

The upsurge in the number of people affected by the COVID-19 is likely to lead to increased rates of emotional trauma and mental illnesses. This article systematically reviewed the available data on the benefits of interventions to reduce adverse mental health sequelae of infectious disease outbreaks, and to offer guidance for mental health service responses to infectious disease pandemic. PubMed, Web of Science, Embase, PsycINFO, WHO Global Research Database on infectious disease, and the preprint server medRxiv were searched. Of 4278 reports identified, 32 were included in this review. Most articles of psychological interventions were implemented to address the impact of COVID-19 pandemic, followed by Ebola, SARS, and MERS for multiple vulnerable populations. Increasing mental health literacy of the public is vital to prevent the mental health crisis under the COVID-19 pandemic. Group-based cognitive behavioral therapy, psychological first aid, community-based psychosocial arts program, and other culturally adapted interventions were reported as being effective against the mental health impacts of COVID-19, Ebola, and SARS. Culturally-adapted, cost-effective, and accessible strategies integrated into the public health emergency response and established medical systems at the local and national levels are likely to be an effective option to enhance mental health response capacity for the current and for future infectious disease outbreaks. Tele-mental healthcare services were key central components of stepped care for both infectious disease outbreak management and routine support; however, the usefulness and limitations of remote health delivery should also be recognized.

We perform a numerical study of the heat transfer and flow structure of Rayleigh–Bénard (RB) convection in (in most cases regular) porous media, which are comprised of circular, solid obstacles located on a square lattice. This study is focused on the role of porosity $\unicode[STIX]{x1D719}$ in the flow properties during the transition process from the traditional RB convection with $\unicode[STIX]{x1D719}=1$ (so no obstacles included) to Darcy-type porous-media convection with $\unicode[STIX]{x1D719}$ approaching 0. Simulations are carried out in a cell with unity aspect ratio, for Rayleigh number $Ra$ from $10^{5}$ to $10^{10}$ and varying porosities $\unicode[STIX]{x1D719}$, at a fixed Prandtl number $Pr=4.3$, and we restrict ourselves to the two-dimensional case. For fixed $Ra$, the Nusselt number $Nu$ is found to vary non-monotonically as a function of $\unicode[STIX]{x1D719}$; namely, with decreasing $\unicode[STIX]{x1D719}$, it first increases, before it decreases for $\unicode[STIX]{x1D719}$ approaching 0. The non-monotonic behaviour of $Nu(\unicode[STIX]{x1D719})$ originates from two competing effects of the porous structure on the heat transfer. On the one hand, the flow coherence is enhanced in the porous media, which is beneficial for the heat transfer. On the other hand, the convection is slowed down by the enhanced resistance due to the porous structure, leading to heat transfer reduction. For fixed $\unicode[STIX]{x1D719}$, depending on $Ra$, two different heat transfer regimes are identified, with different effective power-law behaviours of $Nu$ versus $Ra$, namely a steep one for low $Ra$ when viscosity dominates, and the standard classical one for large $Ra$. The scaling crossover occurs when the thermal boundary layer thickness and the pore scale are comparable. The influences of the porous structure on the temperature and velocity fluctuations, convective heat flux and energy dissipation rates are analysed, further demonstrating the competing effects of the porous structure to enhance or reduce the heat transfer.

Clonorchis sinensis (C. sinensis) is one of the most serious food-borne parasites, which can lead to liver fibrosis or cholangiocarcinoma. Effective measures for clonorchiasis prevention are still urgently needed. Bacillus subtilis (B. subtilis) is an effective antigen delivery platform for oral vaccines. Chonorchis sinensis serpin (CsSerpin) was proved to be potential vaccine candidates. In this study, CsSerpin3 was displayed on the surface of B. subtilis spore and recombinant spores were orally administrated to BALB/C mice. CsSerpin3-specific IgA levels in faecal, bile and intestinal mucous increased at 4–8 weeks after the first administration compared with those in control groups. The mucus production and the number of goblet cells in intestinal mucosa elevated in B.s-CotC-CsSerpin3 (CotC, coat protein of B. subtilis spore) spores treated group compared to those in blank control. No significant difference in the activities of glutamic-pyruvic transaminase/ alanine aminotransferase and glutamic oxalacetic transaminase/aspartate aminotransferase were observed between groups. There was no side effect inflammation and observable pathological damage in the liver tissue of mice after administration. Moreover, collagen deposition and Ishak score were statistically reduced in B.s-CotC-CsSerpin3 spores treated mice. In conclusion, B. subtilis spores displaying CsSerpin3 could be investigated further as an oral vaccine against clonorchiasis.

We present direct numerical simulations of non-Oberbeck–Boussinesq (NOB) Rayleigh–Bénard (RB) convection due to large temperature differences in two-dimensional (2-D) and three-dimensional (3-D) cells. Perfect air is chosen as the operating fluid and the Prandtl number ($Pr$) is fixed to 0.71 for the reference state $\hat{T}_{0}=300~\text{K}$. In the present system, we consider large temperature differences ranging from 60 K to 240 K, and relatively strong NOB effects are induced at moderate Rayleigh numbers ($Ra$) in the range $3\times 10^{6}\leqslant Ra\leqslant 5\times 10^{9}$. The large temperature difference also induces the turbulence system with large density variation. Due to top-down symmetry breaking under NOB conditions, an increase of the centre temperature $T_{c}$ is found compared to the arithmetic mean temperature $T_{m}$ of the top and bottom plates, and the shift of $T_{c}$ is strongly dependent on Rayleigh number $Ra$ and temperature differential $\unicode[STIX]{x1D716}$. The NOB effects on the Nusselt number ($Nu$) are quite small (${\lesssim}2\,\%$). The power-law scalings of $Nu$ versus $Ra$ are robust against NOB effects, even for the extremely large temperature difference 240 K, which has never been reached in previous experiments and simulations. The Reynolds numbers $Re$, as well as the scalings of $Re$ versus $Ra$, are also insensitive to NOB effects. It is noteworthy that the influence of NOB effects on $Nu$ and $Re$ in 3-D RB flow are weaker than its 2-D counterpart. Furthermore, the extended laminar boundary layer (BL) equations are developed based on the low-Mach-number Navier–Stokes equations, which qualitatively predicts the NOB effects on velocity profiles. Direct numerical simulation results indicate that the top and bottom thermal BLs can compensate each other much better than the velocity BLs under NOB conditions, which contribute to the robustness of $Nu$.