Melting and solidification in periodically time-modulated thermal convection are relevant for numerous natural and engineering systems, for example, glacial melting under periodic sun radiation and latent thermal energy storage under periodically pulsating heating. It is highly relevant for the estimation of melt rate and melt efficiency management. However, even the dynamics of a solid–liquid interface shape subjected to a simple sinusoidal heating has not yet been investigated in detail. In this paper, we offer a better understanding of the modulation frequency dependence of the melting and solidification front. We numerically investigate periodic melting and solidification in turbulent convective flow with the solid above and the melted liquid below, and sinusoidal heating at the bottom plate with the mean temperature equal to the melting temperature. We investigate how the periodic heating can prevent the full solidification, and the resulting flow structures and the quasi-equilibrium interface height. We further study the dependence on the heating modulation frequency. As the frequency decreases, we found two distinct regimes, which are ‘partially solid’ and ‘fully solid’. In the fully solid regime, the liquid freezes completely, and the effect of the modulation is limited. In the partially solid regime, the solid partially melts, and a steady or unsteady solid–liquid interface forms depending on the frequency. The interface height can be derived based on the energy balance through the interface. In the partially solid regime, the interface height oscillates periodically, following the frequency of modulation. Here, we propose a perturbation approach that can predict the dependency of the oscillation amplitude on the modulation frequency.

For dissolving active oil droplets in an ambient liquid, it is generally assumed that the Marangoni effect results in repulsive interactions, while the buoyancy effects caused by the density difference between the droplets, diffusing product and the ambient fluid are usually neglected. However, it has been observed in recent experiments that active droplets can form clusters due to buoyancy-driven convection (Krüger et al., Eur. Phys. J. E, vol. 39, 2016, pp. 1–9). In this study we numerically analyse the buoyancy effect, in addition to the propulsion caused by Marangoni flow (with its strength characterized by the Péclet number $Pe$). The buoyancy effects have their origin in (i) the density difference between the droplet and the ambient liquid, which is characterized by the Galileo number $Ga$; and (ii) the density difference between the diffusing product (i.e. filled micelles) and the ambient liquid, which can be quantified by a solutal Rayleigh number $Ra$. We analyse how the attracting and repulsing behaviour of neighbouring droplets depends on the control parameters $Pe$, $Ga$ and $Ra$. We find that while the Marangoni effect leads to the well-known repulsion between the interacting droplets, the buoyancy effect of the reaction product leads to buoyancy-driven attraction. At sufficiently large $Ra$, even collisions between the droplets can take place. Our study on the effect of $Ga$ further shows that with increasing $Ga$, the collision becomes delayed. Moreover, we derive that the attracting velocity of the droplets, which is characterized by a Reynolds number $Re_d$, is proportional to $Ra^{1/4}/( \ell /R)$, where $\ell /R$ is the distance between the neighbouring droplets normalized by the droplet radius. Finally, we numerically obtain the repulsive velocity of the droplets, characterized by a Reynolds number $Re_{rep}$, which is proportional to $PeRa^{-0.38}$. The balance of attractive and repulsive effect leads to $Pe\sim Ra^{0.63}$, which agrees well with the transition curve between the regimes with and without collision.

The global transition towards diets high in calories has contributed to 2.1 billion people becoming overweight, or obese, which damages male reproduction and harms offspring. Recently, more and more studies have shown that paternal exposure to stress closely affects the health of offspring in an intergenerational and transgenerational way. SET Domain Containing 2 (SETD2), a key epigenetic gene, is highly conserved among species, is a crucial methyltransferase for converting histone 3 lysine 36 dimethylation (H3K36me2) into histone 3 lysine 36 trimethylation (H3K36me3), and plays an important regulator in the response to stress. In this study, we compared patterns of SETD2 expression and the H3K36me3 pattern in pre-implantation embryos derived from normal or obese mice induced by high diet. The results showed that SETD2 mRNA was significantly higher in the high-fat diet (HFD) group than the control diet (CD) group at the 2-cell, 4-cell, 8-cell, and 16-cell stages, and at the morula and blastocyst stages. The relative levels of H3K36me3 in the HFD group at the 2-cell, 4-cell, 8-cell, 16-cell, morula stage, and blastocyst stage were significantly higher than in the CD group. These results indicated that dietary changes in parental generation (F0) male mice fed a HFD were traceable in SETD2/H3K36me3 in embryos, and that a paternal high-fat diet brings about adverse effects for offspring that might be related to SETD2/H3K36me3, which throws new light on the effect of paternal obesity on offspring from an epigenetic perspective.

We numerically investigate turbulent Rayleigh–Bénard convection with gas bubbles attached to the hot plate, mimicking a core feature in electrolysis, catalysis or boiling. The existence of bubbles on the plate reduces the global heat transfer due to the much lower thermal conductivity of gases as compared with liquids and changes the structure of the boundary layers. The numerical simulations are performed in three dimensions at Prandtl number $\mbox{Pr}=4.38$ (water) and Rayleigh number $10^7\leqslant \mbox{Ra}\leqslant 10^8$. For simplicity, we assume the bubbles to be equally sized and having pinned contact lines. We vary the total gas-covered area fraction $0.18 \leqslant S_0 \leqslant 0.62$, the relative bubble height $0.02\leqslant h/H \leqslant 0.05$ (where $H$ is the height of the Rayleigh–Bénard cell), the bubble number $40 \leqslant n \leqslant 144$ and their spatial distribution. In all cases, asymmetric temperature profiles are observed, which we quantitatively explain based on the heat flux conservation at each horizontal section. We further propose the idea of using an equivalent single-phase set-up to mimic the system with attached bubbles. Based on this equivalence, we can calculate the heat transfer. Without introducing any free parameter, the predictions for the Nusselt number, the upper and lower thermal boundary layer thicknesses and the mean centre temperature agree well with the numerical results. Finally, our predictions also work for the cases with much larger $\mbox{Pr}$ (e.g. $400$), which indicates that our results can also be applied to predict the mass transfer in water electrolysis with bubbles attached to the electrode surface or in catalysis.

We numerically investigate turbulent Rayleigh–Bénard convection through two immiscible fluid layers, aiming to understand how the layer thickness and fluid properties affect the heat transfer (characterized by the Nusselt number $\mbox {Nu}$) in two-layer systems. Both two- and three-dimensional simulations are performed at fixed global Rayleigh number $\mbox {Ra}=10^8$, Prandtl number $\mbox {Pr}=4.38$ and Weber number $\mbox {We}=5$. We vary the relative thickness of the upper layer between $0.01 \le \alpha \le 0.99$ and the thermal conductivity coefficient ratio of the two liquids between $0.1 \le \lambda _k \le 10$. Two flow regimes are observed. In the first regime at $0.04\le \alpha \le 0.96$, convective flows appear in both layers and $\mbox {Nu}$ is not sensitive to $\alpha$. In the second regime at $\alpha \le 0.02$ or $\alpha \ge 0.98$, convective flow only exists in the thicker layer, while the thinner one is dominated by pure conduction. In this regime, $\mbox {Nu}$ is sensitive to $\alpha$. To predict $\mbox {Nu}$ in the system in which the two layers are separated by a unique interface, we apply the Grossmann–Lohse theory for both individual layers and impose heat flux conservation at the interface. Without introducing any free parameter, the predictions for $\mbox {Nu}$ and for the temperature at the interface agree well with our numerical results and previous experimental data.

This numerical study presents a simple but extremely effective way to considerably enhance heat transport in turbulent wall-bounded multiphase flows, namely by using oleophilic walls. As a model system, we pick the Rayleigh–Bénard set-up, filled with an oil–water mixture. For oleophilic walls, using only $10\,\%$ volume fraction of oil in water, we observe a remarkable heat transport enhancement of more than $100\,\%$ as compared to the pure water case. In contrast, for oleophobic walls, the enhancement is only of about $20\,\%$ as compared to pure water. The physical explanation of the heat transport increment for oleophilic walls is that thermal plumes detach from the oil-rich boundary layer and carry the heat with them. In the bulk, the oil–water interface prevents the plumes from mixing with the turbulent water bulk and to diffuse their heat. To confirm this physical picture, we show that the minimum amount of oil necessary to achieve the maximum heat transport is set by the volume fraction of the thermal plumes. Our findings provide guidelines of how to optimize heat transport in wall-bounded thermal turbulence. Moreover, the physical insight of how coherent structures are coupled with one of the phases of a two-phase system has very general applicability for controlling transport properties in other turbulent wall-bounded multiphase flows.

Celestial navigation is an important means of maritime navigation; it can automatically achieve inertially referenced positioning and orientation after a long period of development. However, the impact of different accuracy of observations and the influence of nonstationary states, such as ship speed change and steering, are not taken into account in existing algorithms. To solve this problem, this paper proposes an adaptively robust maritime celestial navigation algorithm, in which each observation value is given an equivalent weight according to the robust estimation theory, and the dynamic balance between astronomical observation and prediction values of vessel motion is adjusted by applying the adaptive factor. With this system, compared with the frequently used least square method and extended Kalman filter algorithm, not only are the real-time and high-precision navigation parameters, such as position, course, and speed for the vessel, calculated simultaneously, but also the influence of abnormal observation and vessel motion status change could be well suppressed.

The present study aimed to explore the association between dietary patterns in abdominal obesity obtained by reduced-rank regression (RRR) with visceral fat index (VFI) as a dependent variable and dyslipidemia in rural adults in Henan, China. A total of 29538 people aged 18–79 were selected from the Henan Rural Cohort Study. RRR analysis was used to identify dietary patterns; logistic regression analysis and restricted cubic spline regression models were applied to analyze the association between dietary patterns in abdominal obesity and dyslipidemia. VFI was used as a mediator to estimate the mediation effect. The dietary pattern in abdominal obesity was characterized by high carbohydrate and red meat intake and low consumption of fresh fruits, vegetables, milk, etc. After full adjustment, the highest quartile of dietary pattern scores was significantly associated with an increased risk of dyslipidemia (OR: 1·33, 95 % CI 1·23–1·44, P trend < 0·001), there was a non-linear dose–response relationship between them (P overall-association < 0·001, P non-lin-association = 0·022). The result was similar in dose-response between the dietary pattern scores and VFI. The indirect effect partially mediated by VFI was significant (OR: 1·07, 95 % CI 1·06–1·08). VIF explained approximately 53·3 % of odds of dyslipidemia related to the dietary pattern. Abdominal obesity dietary pattern scores positively affected VFI and dyslipidemia; there was a dose-response in both relationships. Dyslipidemia progression increased with higher abdominal obesity dietary pattern scores. In addition, VFI played a partial mediating role in relationship between abdominal obesity dietary pattern and dyslipidemia.

The present study evaluated whether fat mass assessment using the triceps skinfold (TSF) thickness provides additional prognostic value to the Global Leadership Initiative on Malnutrition (GLIM) framework in patients with lung cancer (LC). We performed an observational cohort study including 2672 LC patients in China. Comprehensive demographic, disease and nutritional characteristics were collected. Malnutrition was retrospectively defined using the GLIM criteria, and optimal stratification was used to determine the best thresholds for the TSF. The associations of malnutrition and TSF categories with survival were estimated independently and jointly by calculating multivariable-adjusted hazard ratios (HR). Malnutrition was identified in 808 (30·2 %) patients, and the best TSF thresholds were 9·5 mm in men and 12 mm in women. Accordingly, 496 (18·6 %) patients were identified as having a low TSF. Patients with concurrent malnutrition and a low TSF had a 54 % (HR = 1·54, 95 % CI = 1·25, 1·88) greater death hazard compared with well-nourished individuals, which was also greater compared with malnourished patients with a normal TSF (HR = 1·23, 95 % CI = 1·06, 1·43) or malnourished patients without TSF assessment (HR = 1·31, 95 % CI = 1·14, 1·50). These associations were concentrated among those patients with adequate muscle mass (as indicated by the calf circumference). Additional fat mass assessment using the TSF enhances the prognostic value of the GLIM criteria. Using the population-derived thresholds for the TSF may provide significant prognostic value when used in combination with the GLIM criteria to guide strategies to optimise the long-term outcomes in patients with LC.

We theoretically and numerically investigate the instabilities driven by diffusiophoretic flow, caused by a solutal concentration gradient along a reacting surface. The important control parameters are the Péclet number $Pe$, which quantifies the ratio of the solutal advection rate to the diffusion rate, and the Schmidt number $Sc$, which is the ratio of viscosity and diffusivity. First, we study the diffusiophoretic flow on a catalytic plane in two dimensions. From a linear stability analysis, we obtain that for $Pe$ larger than $8{\rm \pi}$ mass transport by convection overtakes that by diffusion, and a symmetry-breaking mode arises, which is consistent with numerical results. For even larger $Pe$, nonlinear terms become important. For $Pe > 16{\rm \pi}$, multiple concentration plumes are emitted from the catalytic plane, which eventually merge into a single larger one. When $Pe$ is even larger ($Pe \gtrsim 603$ for Schmidt number $Sc=1$), there are continuous emissions and merging events of the concentration plumes. The newly found flow states have different flow structures for different $Sc$: for $Sc\geqslant 1$, we observe the chaotic emission of plumes, but the fluctuations of concentration are only present in the region near the catalytic plane. In contrast, for $Sc<1$, chaotic flow motion occurs also in the bulk. In the second part of the paper, we conduct three-dimensional simulations for spherical catalytic particles, and beyond a critical Péclet number again find continuous plume emission and plume merging, now leading to a chaotic motion of the phoretic particle. Our results thus help us to understand the experimentally observed chaotic motion of catalytic particles in the high $Pe$ regime.

It is commonly accepted that the breakup criteria of drops or bubbles in turbulence is governed by surface tension and inertia. However, also buoyancy can play an important role at breakup. In order to better understand this role, here we numerically study two-dimensional Rayleigh–Bénard convection for two immiscible fluid layers, in order to identify the effects of buoyancy on interface breakup. We explore the parameter space spanned by the Weber number $5\leqslant We \leqslant 5000$ (the ratio of inertia to surface tension) and the density ratio between the two fluids $0.001 \leqslant \varLambda \leqslant 1$, at fixed Rayleigh number $Ra=10^8$ and Prandtl number $Pr=1$. At low $We$, the interface undulates due to plumes. When $We$ is larger than a critical value, the interface eventually breaks up. Depending on $\varLambda$, two breakup types are observed. The first type occurs at small $\varLambda \ll 1$ (e.g. air–water systems) when local filament thicknesses exceed the Hinze length scale. The second, strikingly different, type occurs at large $\varLambda$ with roughly $0.5 < \varLambda \leqslant 1$ (e.g. oil–water systems): the layers undergo a periodic overturning caused by buoyancy overwhelming surface tension. For both types, the breakup criteria can be derived from force balance arguments and show good agreement with the numerical results.

Hypertension represents one of the most common pre-existing conditions and comorbidities in Coronavirus disease 2019 (COVID-19) patients. To explore whether hypertension serves as a risk factor for disease severity, a multi-centre, retrospective study was conducted in COVID-19 patients. A total of 498 consecutively hospitalised patients with lab-confirmed COVID-19 in China were enrolled in this cohort. Using logistic regression, we assessed the association between hypertension and the likelihood of severe illness with adjustment for confounders. We observed that more than 16% of the enrolled patients exhibited pre-existing hypertension on admission. More severe COVID-19 cases occurred in individuals with hypertension than those without hypertension (21% vs. 10%, P = 0.007). Hypertension associated with the increased risk of severe illness, which was not modified by other demographic factors, such as age, sex, hospital geological location and blood pressure levels on admission. More attention and treatment should be offered to patients with underlying hypertension, who usually are older, have more comorbidities and more susceptible to cardiac complications.

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.

Rabbits play an important role in people’s lives due to their high nutritional value and high-quality hair that can be used as raw material for textiles. Furthermore, rabbits are an important animal model for human disease, as genome-edited animals are particularly valuable for studying gene functions and pathogenesis. Somatic cell nuclear transfer (SCNT) is an important technique for producing genome-edited animals and it has great value in saving endangered species and in clone stem cell therapy. However, the low efficiency of SCNT limits its application, with the selection of suitable rabbit oocytes being crucial to its success. In the present study, we collected oocytes from ovarian follicles and stained them with 26 μM brilliant cresyl blue (BCB). We then matured the oocytes in vitro and used them for SCNT. Comparison of the BCB-positive oocytes with BCB-negative oocytes and the control group showed that the BCB-positive group had a significantly higher maturation rate (81.4% vs. 48.9% and 65.3% for the negative and control groups, respectively), cleavage rate (86.6% vs. 67.9% and 77.9%), blastocyst rate (30.5% vs. 12.8% and 19.6%), total number of blastocysts (90±7.5 vs. 65.3±6.3 and 67.5±5.7), and inner cell mass (ICM)/ trophectoderm (TE) index (42.3±4.2 vs. 30.2±2.1 and 33.9±5.1) (P<0.05). The BCB-positive group had a significantly lower apoptosis index (2.1±0.6 vs. 8.2±0.9 and 6.7±1.1 for the negative and control groups, respectively) (P<0.05). These findings demonstrate that BCB-positive oocytes have a higher maturation ability and developmental competence in vitro, indicating that BCB staining is a reliable method for selecting oocytes to enhance the efficiency of SCNT.

In this paper, we review the status of the multifunctional experimental platform at the National Laboratory of High Power Laser and Physics (NLHPLP). The platform, including the SG-II laser facility, SG-II 9th beam, SG-II upgrade (SG-II UP) facility, and SG-II 5 PW facility, is operational and available for interested scientists studying inertial confinement fusion (ICF) and a broad range of high-energy-density physics. These facilities can provide important experimental capabilities by combining different pulse widths of nanosecond, picosecond, and femtosecond scales. In addition, the SG-II UP facility, consisting of a single petawatt system and an eight-beam nanosecond system, is introduced including several laser technologies that have been developed to ensure the performance of the facility. Recent developments of the SG-II 5 PW facility are also presented.

The aim of the research reported in this Research Communication was to identify differentially expressed proteins in dairy cows with normal and lutein diet and to elucidate the mechanisms of lutein-induced effects on bovine mammary gland metabolism using a comparative proteomic approach. Thirty-three differentially expressed proteins were identified from mammary gland of control diet-fed and lutein diet-fed dairy cows. Among these proteins, 15 were upregulated and 18 were downregulated in the lutein group. Functional analysis of the differentially expressed proteins showed that increased blood flow, depressed glycolysis, enhanced lactose anabolism, decreased fatty acid oxidation and up-regulated beta lactoglobulin expression were connected with lutein addition. These results suggested that the increased blood flow, reduced glucose catabolism, enhanced capacity for milk lactose synthesis, depressed fatty acid catabolism and increased expression of antioxidantion related protein may be the prime factors contributing to the increased milk production and enhanced immune status in lutein-fed dairy cows. This study provides molecular mechanism of dietary lutein in regulating lactation of dairy cows.

The hot deformation behavior of Nb–V–Ti microalloyed ultra-high strength steel was investigated by isothermal compression at 900–1200 °C with strain rates from 0.01 to 10 s−1. The microstructure evolution and precipitation behavior were studied using an optical microscope and a transmission electron microscope Results indicate that the peak stress of experimental steel increases with increasing the strain rate and decreasing the deformation temperature. The constitutive equation of hot deformation was developed with the activation energy Q being about 407.29 kJ/mol. The processing maps were also obtained to identify the instable regions of the flow behavior and to evaluate the efficiency of hot deformation. The size of dynamically recrystallized grains increases gradually with a decrease in the strain rate. Three types of carbides were identified, namely M3C, rich-Ti MC, and rich-Nb MC. With the increase of the deformation rate, the amounts of carbides increase, and the average sizes of the carbides decrease gradually.