We study the diffusive instability subject to a solid–liquid interface and a Kolmogorov flow using modal, non-modal analyses and energy analysis. It is found that the phase boundary has different effects on the threshold of diffusive convection for weak and strong salinity stratification. In the context of oceanography where the salinity Rayleigh number RS is very high, the ice–water interface has negligible influence on the onset of diffusive convection. In the presence of shear, the diffusive convection for RS < 106 tends to be inhibited and with the increase of shear intensity, the oscillatory, steady convective and Kelvin–Helmholtz instabilities will be successively dominant. For RS > 106, the shear has a destabilizing effect on the diffusive convection, due to the generation of thermohaline-sheared instability found by Radko (J. Fluid Mech., vol. 805, 2016, pp. 147–170). Non-modal analysis indicates that for realistic parameters of high-latitude oceans, with the transition of the ultimate energy source of instability from the density gradient to background current, the thermohaline-shear instability is expected to transition from oscillatory to steady instability. The initial transient amplification due to double diffusion has also been studied, which is due to the generation of overstable instability at the initial phase. For RS < 106, the optimal initial condition to achieve the maximum transient growth favours longitudinal rolls. For thermohaline-shear instability, however, it favours transverse rolls and specifically, in oscillatory thermohaline-shear instability, the transient amplification can be enhanced by the shear by one order of magnitude, thus having important influence on the stability of the system.

We provide behavioral insights into the economic substitution phenomenon among front-loaded and back-loaded tax-sheltered savings plans. We conduct three behavioral studies with adult Canadian participants to show experimentally that substitution can occur and to explain why substitution can occur. The first study shows that participants transfer savings from a front-loaded plan to a back-loaded plan when the latter becomes available, consistent with a substitution effect. The second study examines how participants trade-off two unique features of back-loaded and front-loaded savings plans. Our results indicate that participants favor the back-loaded tax feature and a variable contribution limit (offered in a front-loaded plan). As participants prefer one feature from each type of plan, this finding can help explain why substitution occurs. The third study provides participants a categorization task with various household budgeting items, including savings items. Results show that 68.1% of participants categorize multiple tax-sheltered savings plans in the same mental account, again consistent with a substitution effect under a budget constraint. As both tax-sheltered savings plans in Canada are used for different purposes, this finding shows that participants tend not to distinguish between the purpose of saving in each account, consistent with a substitution phenomenon.

As optical parametric chirped pulse amplification has been widely adopted for the generation of extreme intensity laser sources, nonlinear crystals of large aperture are demanded for high-energy amplifiers. Yttrium calcium oxyborate (YCa4O(BO3)3, YCOB) is capable of being grown with apertures exceeding 100 mm, which makes it possible for application in systems of petawatt scale. In this paper, we experimentally demonstrated for the first time to our knowledge, an ultra-broadband non-collinear optical parametric amplifier with YCOB for petawatt-scale compressed pulse generation at 800 nm. Based on the SG-II 5 PW facility, amplified signal energy of approximately 40 J was achieved and pump-to-signal conversion efficiency was up to 42.3%. A gain bandwidth of 87 nm was realized and supported a compressed pulse duration of 22.3 fs. The near-field and wavefront aberration represented excellent characteristics, which were comparable with those achieved in lithium triborate-based amplifiers. These results verified the great potential for YCOB utilization in the future.

In the present study, we investigated the influence of different mid-stage N compensation timings on agronomic and physiological traits associated with grain yield and quality in field experiments. Two japonica rice cultivars with a good tasting quality (Nangeng 9108 and Nangeng 5055) were examined under eight N compensation timings (N1–N6: one-time N compensation at 7-2 weeks before heading; N7: split N compensation at 5 and 3 weeks before heading; N8: split N compensation at 4 and 2 weeks before heading) and a control with no N compensation. The highest yield was obtained with N7, followed by N3. The yield advantage is mainly attributable to the improved population structure (higher productive tiller rate with a stable number of effective panicles), higher total number of spikelets per unit area (large panicles with more grains per panicle), larger leaf area index in the late period and higher photosynthetic production capacity (more dry matter accumulation and transportation in the middle and late periods). Delaying N compensation timing improved the processing and nutritional quality of rice, but decreased the quality of appearance and cooking/eating traits. Our results suggest that, from the perspective of achieving relative coordination between high yield and high quality of japonica rice, the optimal N compensation should be divided equally at 5 and 3 weeks before heading. However, if simplifying the number of operations and the pursuit of eating quality were considered, one-time N compensation should be conducted at 5 weeks before heading.

Disequilibrating macro shocks affect different firms’ prospects differently, increasing idiosyncratic variation in forward-looking stock returns before affecting economic growth. Consistent with most such shocks from 1947 to 2020 enhancing productivity, increased idiosyncratic stock return variation forecasts next-quarter real GDP growth, industrial production growth, and consumption growth both in-sample and out-of-sample. These effects persist after controlling for other leading economic indicators.

We present an experimental study of Rayleigh–Bénard convection using liquid metal alloy gallium-indium-tin as the working fluid with a Prandtl number of $Pr=0.029$. The flow state and the heat transport were measured in a Rayleigh number range of $1.2\times 10^{4} \le Ra \le 1.3\times 10^{7}$. The temperature fluctuation at the cell centre is used as a proxy for the flow state. It is found that, as $Ra$ increases from the lower end of the parameter range, the flow evolves from a convection state to an oscillation state, a chaotic state and finally a turbulent state for $Ra>10^5$. The study suggests that the large-scale circulation in the turbulent state is a residual of the cell structure near the onset of convection, which is in contrast with the case of $Pr\sim 1$, where the cell structure is transiently replaced by high order flow modes before the emergence of the large-scale circulation in the turbulent state. The evolution of the flow state is also reflected by the heat transport characterised by the Nusselt number $Nu$ and the probability density function (p.d.f.) of the temperature fluctuation at the cell centre. It is found that the effective local heat transport scaling exponent $\gamma$, i.e. $Nu\sim Ra^{\gamma }$, changes continuously from $\gamma =0.49$ at $Ra\sim 10^4$ to $\gamma =0.25$ for $Ra>10^6$. Meanwhile, the p.d.f. at the cell centre gradually evolves from a Gaussian-like shape before the transition to turbulence to an exponential-like shape in the turbulent state. For $Ra>10^6$, the flow shows self-similar behaviour, which is revealed by the universal shape of the p.d.f. of the temperature fluctuation at the cell centre and a $Nu=0.19Ra^{0.25}$ scaling for the heat transport.

The mitochondrial genome provides important information for phylogenetic analysis and an understanding of evolutionary origin. In this study, the mitochondrial genomes of Ilisha elongata and Setipinna tenuifilis were sequenced, which are typical circular vertebrate mitochondrial genomes composed of 16,770 and 16,805 bp, respectively. The mitogenomes of I. elongata and S. tenuifilis include 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA), two ribosomal RNA (rRNA) genes and one control region (CR). Both two species' genome compositions were highly A + T biased and exhibited positive AT-skews and negative GC-skews. The genetic distance and Ka/Ks ratio analyses indicated that 13 PCGs were affected by purifying selection and the selection pressures were different from certain deep-sea fishes, which were most likely due to the difference in their living environment. Results of phylogenetic analysis support close relationships among Chirocentridae, Denticipitidae, Clupeidae, Engraulidae and Pristigasteridae based on the nucleotide and amino acid sequences of 13 PCGs. Within Clupeoidei, I. elongata and S. tenuifilis were most closely related to the family Pristigasteridae and Engraulidae, respectively. These results will help to better understand the evolutionary position of Clupeiformes and provide a reference for further phylogenetic research on Clupeiformes species.

It has been suggested that added sugar intake is associated with non-alcoholic fatty liver disease (NAFLD). However, previous studies only focused on sugar-sweetened beverages; the evidence for associations with total added sugars and their sources is scarce. This study aimed to examine the associations of total added sugars, their physical forms (liquid v. solid) and food sources with risk of NAFLD among adults in Tianjin, China. We used data from 15 538 participants, free of NAFLD, other liver diseases, CVD, cancer or diabetes at baseline (2013–2018 years). Added sugar intake was estimated from a validated 100-item FFQ. NAFLD was diagnosed by ultrasonography after exclusion of other causes of liver diseases. Multivariable Cox proportional hazards models were fitted to calculate hazard ratios (HR) and corresponding 95 % CI for NAFLD risk with added sugar intake. During a median follow-up of 4·2 years, 3476 incident NAFLD cases were documented. After adjusting for age, sex, BMI and its change from baseline to follow-up, lifestyle factors, personal and family medical history and overall diet quality, the multivariable HR of NAFLD risk were 1·18 (95 % CI 1·06, 1·32) for total added sugars, 1·20 (95 % CI 1·08, 1·33) for liquid added sugars and 0·96 (95 % CI 0·86, 1·07) for solid added sugars when comparing the highest quartiles of intake with the lowest quartiles of intake. In this prospective cohort of Chinese adults, higher intakes of total added sugars and liquid added sugars, but not solid added sugars, were associated with a higher risk of NAFLD.

We study the two-dimensional Rayleigh–Bénard instability subject to the combined effects of a solid–liquid phase boundary and shear using linear stability theory and energy analysis. We consider two thermal states of the solid (isothermal and conducting), and two types of shear that can arise in different contexts. When the melting temperature is equal to that held at the top boundary, three instability modes can arise with the increase of Reynolds number Re that characterizes the shear intensity: a boundary mode, a mixed boundary–bulk mode and a bulk flow mode. When the melting temperature lies between the top and the bottom boundaries, the introduction of Couette flow, independent of its intensity, always leads to the mixed mode, whereas the instability with Poiseuille flow is dominated by the bulk flow mode once Re exceeds a critical value, below which the mixed mode dominates. The energy analysis suggests that there exist two mechanisms by which the shear flow affects the system: one is by inhibiting the upward heat flux and another is by absorbing energy from the perturbed hydrodynamic field. These two mechanisms can play totally different roles in different cases. Results in the high-Re regime indicate that, when Re exceeds its classical threshold, i.e. Re = 5772.2 for Poiseuille flow, the Tollmien–Schlichting instability will be dominant in the present system.

We report an experimental study of the viscous boundary layer (BL) properties of turbulent Rayleigh–Bénard convection in a cylindrical cell. The velocity profile with all three components was measured from the centre of the bottom plate by an integrated home-made particle image velocimetry system. The Rayleigh number $Ra$ varied in the range $1.82 \times 10^8 \le Ra \le 5.26 \times 10^9$ and the Prandtl number $Pr$ was fixed at $Pr = 4.34$. The probability density function of the wall-shear stress indicates that using the velocity component in the mean large-scale circulation (LSC) plane alone may not be sufficient to characterise the viscous BL. Based on a dynamic wall-shear frame, we propose a method to reconstruct the measured full velocity profile which eliminates the effects of complex dynamics of the LSC. Various BL properties including the eddy viscosity are then obtained and analysed. It is found that, in the dynamic wall-shear frame, the eddy viscosity profiles along the centre line of the convection cell at different $Ra$ all collapse on a single master curve described by $\nu _t^d / \nu = 0.81 (z / \delta _u^d) ^{3.10 \pm 0.05}$. The Rayleigh number dependencies of several BL quantities are also determined in the dynamic frame, including the BL thickness $\delta _u^d$ (${\sim } Ra^{-0.21}$), the Reynolds number $Re^d$ (${\sim }Ra^{-0.46}$) and the shear Reynolds number $Re_s^d$ (${\sim } Ra^{0.24}$). Within the experimental uncertainty, these scaling exponents are the same as those obtained in the static laboratory frame. Finally, with the measured full velocity profile, we obtain the energy dissipation rate at the centre of the bottom plate $\varepsilon _{w}$, which is found to follow $\langle \varepsilon _{w} \rangle _t \sim Ra^{1.25}$.

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.

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.

This study is performed to figure out how the presence of diabetes affects the infection, progression and prognosis of 2019 novel coronavirus disease (COVID-19), and the effective therapy that can treat the diabetes-complicated patients with COVID-19. A multicentre study was performed in four hospitals. COVID-19 patients with diabetes mellitus (DM) or hyperglycaemia were compared with those without these conditions and matched by propensity score matching for their clinical progress and outcome. Totally, 2444 confirmed COVID-19 patients were recruited, from whom 336 had DM. Compared to 1344 non-DM patients with age and sex matched, DM-COVID-19 patients had significantly higher rates of intensive care unit entrance (12.43% vs. 6.58%, P = 0.014), kidney failure (9.20% vs. 4.05%, P = 0.027) and mortality (25.00% vs. 18.15%, P < 0.001). Age and sex-stratified comparison revealed increased susceptibility to COVID-19 only from females with DM. For either non-DM or DM group, hyperglycaemia was associated with adverse outcomes, featured by higher rates of severe pneumonia and mortality, in comparison with non-hyperglycaemia. This was accompanied by significantly altered laboratory indicators including lymphocyte and neutrophil percentage, C-reactive protein and urea nitrogen level, all with correlation coefficients >0.35. Both diabetes and hyperglycaemia were independently associated with adverse prognosis of COVID-19, with hazard ratios of 10.41 and 3.58, respectively.

Lower-crust-derived adakitic rocks in the Gangdese belt provide important constraints on the timing of Tibetan crustal thickening and on the relative contributions of magmatic and tectonic processes. Here we present geochronological and geochemical data for the Wangdui porphyritic monzogranites in the western Gangdese belt. Zircon U–Pb dating yields emplacement ages of 46–44 Ma. All samples have high Sr (321–599 ppm), low Yb (0.76–1.33 ppm) and Y (10.6–18.3 ppm) contents, with high La/Yb (51.1–72.3) and Sr/Y (21.0–51.4) ratios, indicating adakitic affinities. The low MgO (0.97–1.76 wt %), Cr (7.49–53.6 ppm) and Ni (4.75–29.1 ppm) contents, as well as high 87Sr/86Sr(i) (0.7143–0.7145), low ϵNd(t) (−10.4 to −9.8) and zircon ϵHf(t) (−17.7 to 0.4) values, suggest that the Wangdui pluton most likely originated from partial melting of the thickened ancient lower crust. In combination with previously published data, despite the east–west-trending heterogeneity of crustal composition in the Gangdese belt, the La/Yb ratios of magmatic rocks reveal that both western and eastern segments experienced remarkable crustal thickening in the Eocene. However, in contrast to the thickened juvenile lower crust in the eastern segment formed by the underplating of mantle-derived magmas, tectonic shortening plays a more crucial role in thickening of the ancient basement in western Gangdese. In fact, such Eocene-thickened ancient lower-crust-derived adakitic rocks are widely distributed in the central Himalayan–Tibetan orogen. This, together with the extensive development of fold–thrust belts, suggests that tectonic shortening might be the main mechanism accounting for the crustal thickening associated with the India–Asia collision.