We systematically study the dissipative anomaly in compressible magnetohydrodynamic (MHD) turbulence using direct numerical simulations, and show that the total dissipation remains finite as viscosity diminishes. The dimensionless dissipation rate $\mathcal {C}_{\varepsilon }$ fits well with the model $\mathcal {C}_{\varepsilon } = \mathcal {C}_{\varepsilon,\infty } + \mathcal {D}/R_L^-$ for all levels of flow compressibility considered here, where $R_L^-$ is the generalized large-scale Reynolds number. The asymptotic value $\mathcal {C}_{\varepsilon,\infty }$ describes the total energy transfer flux, and decreases with increase of the flow compressibility, indicating non-universality of the dimensionless dissipation rate in compressible MHD turbulence. After introducing an empirically modified dissipation rate, the data from compressible cases collapse to a form similar to the incompressible MHD case depending only on the modified Reynolds number.

The shock wave accelerating a heavy fluid layer can induce reverberating waves that continuously interact with the first and second interfaces. In order to manipulate the perturbation growths at fluid-layer interfaces, we present a theoretical framework to eliminate the reverberating waves. A model is established to predict the individual freeze-out (i.e. stagnation of perturbation growth) for the first and second interfaces under specific flow conditions determined based on the shock dynamics theory. The theoretical model quantifies the controllable parameters required for freeze-out, including the initial amplitudes of the first and second interfaces, the interface coupling strength and the maximum initial layer thickness preventing the second interface's phase reversal. The effectiveness of the model in predicting individual freeze-out for the first and second interfaces is validated numerically over a wide range of initial conditions. The upper and lower limits of initial amplitudes for the freeze-out of the whole fluid-layer width growth are further predicted. Within this amplitude range, a slightly higher initial amplitude for the second interface is specified, effectively arresting the growth of the entire fluid-layer width before the phase reversal of the second interface.

This work investigates the compressible turbulence induced by Richtmyer–Meshkov (RM) instability using high-resolution Navier–Stokes simulations. Special attention is paid to the characteristics of RM turbulence including the mixing width growth, the turbulent kinetic energy (TKE) decay, the mixing degree, inhomogeneity and anisotropy. Three distinct initial perturbation spectra are designed at the interface to reveal the initial condition imprint on the RM turbulence. Results show that cases with initial large-scale perturbations present a stronger imprint on statistical characteristics and also a quicker growth of mixing width, whereas cases with small-scale perturbations present a faster TKE decay, greater mixing level, higher isotropy and homogeneity. A thorough analysis on the inter-scale energy transfer in RM turbulence is also presented with the coarse-graining approach that exposes the two subfilter-scale (SFS) energy fluxes (i.e. deformation work and baropycnal work). A strong correlation between the nonlinear model of baropycnal work (Fluids, 4(2), 2019) and the simulation results is confirmed for the first time, demonstrating its potential in modelling the RM turbulence. Two primary mechanisms of baropycnal work (the straining and baroclinic generation processes) are explored with this nonlinear model. The evolutions of two SFS energy fluxes exhibit distinct behaviours at various filter scales, in different flow regions and under various flow motions (strain and rotation). It is found that all three cases share the common inter-scale energy transfer dynamics, which is important for modelling the RM turbulence.

Our earth is immersed in the near-earth space plasma environment, which plays a vital role in protecting our planet against the solar-wind impact and influencing space activities. It is significant to investigate the physical processes dominating the environment, for deepening our scientific understanding of it and improving the ability to forecast the space weather. As a crucial part of the National Major Scientific and Technological Infrastructure–Space Environment Simulation Research Infrastructure (SESRI) in Harbin, the Space Plasma Environment Research Facility (SPERF) builds a system to replicate the near-earth space plasma environment in the laboratory. The system aims to simulate the three-dimensional (3-D) structure and processes of the terrestrial magnetosphere for the first time in the world, providing a unique platform to reveal the physics of the 3-D asymmetric magnetic reconnection relevant to the earth's magnetopause, wave–particle interaction in the earth's radiation belt, particles’ dynamics during the geomagnetic storm, etc. The paper will present the engineering design and construction of the near-earth space plasma simulation system of the SPERF, with a focus on the critical technologies that have been resolved to achieve the scientific goals. Meanwhile, the possible physical issues that can be studied based on the apparatus are sketched briefly. The earth-based system is of great value in understanding the space plasma environment and supporting space exploration.

The third-order law links energy transfer rates in the inertial range of magneto- hydrodynamic (MHD) turbulence with third-order structure functions. Anisotropy, a typical property in the solar wind, challenges the applicability of the third-order law with the isotropic assumption. To shed light on the energy transfer process in the presence of anisotropy, we conducted direct numerical simulations of forced MHD turbulence with normal and hyper-viscosity under various strengths of the external magnetic field ($B_0$), and calculated three forms of third-order structure function with or without averaging of the azimuthal or polar angles with respect to $B_0$ direction. Correspondingly, three estimated energy transfer rates were obtained. The result shows that the peak of normalized third-order structure function occurs at larger scales closer to the $B_0$ direction, and the maximum of longitudinal transfer rates shifts away from the $B_0$ direction at larger $B_0$. Compared with normal viscous cases, hyper-viscous cases can attain better separated inertial range, thus facilitating the estimation of the energy cascade rates. We find that the widespread use of the isotropic form of the third-order law in estimating the energy transfer rates is questionable in some cases, especially when the anisotropy arising from the mean magnetic field is inevitable. In contrast, the direction-averaged third-order structure function properly accounts for the effect of anisotropy and predicts the energy transfer rates and inertial range accurately, even at very high $B_0$. With limited statistics, the third-order structure function shows a stronger dependence on averaging of azimuthal angles than the time, especially for high $B_0$ cases. These findings provide insights into the anisotropic effect on the estimation of energy transfer rates.

The scaling law of the structure function of Richtmyer–Meshkov (RM) turbulence is investigated both numerically and theoretically. High-fidelity simulations with a minimum-dispersion, adaptive-dissipation scheme are first performed. Results show that the mixing width experiences an exponential growth and the turbulent kinetic energy has a visible $-3/2$ spectrum. The scalar field exhibits a greater degree of intermittency than the velocity field, and also the small-scale statistics suffer a larger influence of large scales. Visible differences in the scaling law of the structure function among the RM turbulence and other types of turbulence are observed, which reveal the unique characteristic of RM turbulence. A phenomenological theory, which gives the spatial and temporal scaling laws of the structure functions of velocity and scalar of RM turbulence, is developed for the first time by introducing an external agent. The spatial scaling exponents of structure functions from simulation deviate from the Kolmogorov exponents, but are quite close to the RM-modified anomalous exponents. This demonstrates the validity of the present phenomenological theory. The temporal scaling exponents of structure functions first meet the RM-modified anomalous exponents, and then approach the Kolmogorov–Obukhov–Corrsin non-intermittent ones.

Coastal eutrophication and hypoxia remain a persistent environmental crisis despite the great efforts to reduce nutrient loading and mitigate associated environmental damages. Symptoms of this crisis have appeared to spread rapidly, reaching developing countries in Asia with emergences in Southern America and Africa. The pace of changes and the underlying drivers remain not so clear. To address the gap, we review the up-to-date status and mechanisms of eutrophication and hypoxia in global coastal oceans, upon which we examine the trajectories of changes over the 40 years or longer in six model coastal systems with varying socio-economic development statuses and different levels and histories of eutrophication. Although these coastal systems share common features of eutrophication, site-specific characteristics are also substantial, depending on the regional environmental setting and level of social-economic development along with policy implementation and management. Nevertheless, ecosystem recovery generally needs greater reduction in pressures compared to that initiated degradation and becomes less feasible to achieve past norms with a longer time anthropogenic pressures on the ecosystems. While the qualitative causality between drivers and consequences is well established, quantitative attribution of these drivers to eutrophication and hypoxia remains difficult especially when we consider the social economic drivers because the changes in coastal ecosystems are subject to multiple influences and the cause–effect relationship is often non-linear. Such relationships are further complicated by climate changes that have been accelerating over the past few decades. The knowledge gaps that limit our quantitative and mechanistic understanding of the human-coastal ocean nexus are identified, which is essential for science-based policy making. Recognizing lessons from past management practices, we advocate for a better, more efficient indexing system of coastal eutrophication and an advanced regional earth system modeling framework with optimal modules of human dimensions to facilitate the development and evaluation of effective policy and restoration actions.

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

The new mineral, gysinite-(La), with the ideal formula PbLa(CO3)2(OH)⋅H2O, has been discovered in lujavrite from the Saima alkaline complex, Liaoning Province, China. It commonly occurs as subhedral to anhedral, granular and platy crystals of 5 to 50 μm in size, in interstices or enclosed in microcline, aegirine and nepheline. Associated minerals include nepheline, aegirine, microcline, natrolite, eudialyte, lamprophyllite, bastnäsite-(Ce), parasite-(Ce), ancylite-(La), ancylite-(Ce), bobtraillite, britholite-(Ce), thorite, calcite and galena. The crystallisation of gysinite-(La) may be related to the post-magmatic carbonation event. Gysinite-(La) crystals are generally transparent, colourless, or pale yellow, with a vitreous lustre and white streak. It is brittle with an uneven fracture, and the estimated Mohs hardness is 3½ to 4. The calculated density is 5.007 g/cm3. Optically, gysinite-(La) is biaxial (–), α= 1.832(2), β= 1.849(4), γ = 1.862(5) in white light and 2Vmeas = 81.6°. The empirical formula of gysinite-(La) is (La0.93Pb0.61Nd0.23Pr0.14Sr0.04Gd0.02Sm0.01Eu0.01Ca0.01)Σ2(CO3)2(OH)1.34⋅0.66H2O, which is calculated on the basis of general formula (REExM2+2–x)(CO3)2(OH)x⋅(2–x)H2O. The strongest eight lines of its powder X-ray diffraction pattern [d, Å (I, %) (hkl)] are: 5.596 (21) (011), 4.349 (100) (110), 3.732 (68) (111), 2.984 (61) (121), 2.667 (21) (031), 2.363 (48) (131), 2.090 (29) (221) and 2.028 (21) (212). Gysinite-(La) is orthorhombic, in the space group Pmcn, and unit-cell parameters refined from single-crystal X-ray diffraction data are: a = 5.0655(2) Å, b = 8.5990(3) Å, c = 7.3901(4) Å, V = 321.90(2) Å3 and Z = 2. It is a new member of the ancylite group and isostructural with gysinite-(Nd), but with La and Pb dominant in the metal cation sites in the structure.

During late Carboniferous time, the residual ocean basin gradually closed in West Junggar and only a small amount of seawater remained in the Hala’alat Mountain area, where discussions of provenance and tectonics are limited. In this study, LA-ICP-MS U–Pb dating and heavy mineral identification are conducted on the upper Carboniferous tuffaceous sandstones from the Hala’alat and Aladeyikesai formations in the Hala’alat Mountain area. The results reveal the low maturity of the clastic sediments, indicating proximal deposition. The Hala’alat Formation detrital zircons present a single peak (c. 330 Ma). Speculatively, the primary provenance is the Boshchekul–Chingiz Arc, and the secondary sources are the Darbut Tectono-Magmatic Belt and island arcs in the basin. The main peak and provenance of the Aladeyikesai Formation are similar to those of the Hala’alat Formation. Moreover, several age groups, namely, 370–344 Ma, 427–404 Ma and 478–476 Ma, potentially correspond to provenances of the Darbut Tectono-Magmatic Belt, the Boshchekul–Chingiz Arc and the Kujibai–Hongguleleng Ophiolitic Mélange Belt. In addition, the maximum depositional ages of the Hala’alat and Aladeyikesai formations calculated are 314.6 ± 1.54 Ma and 330.8 ± 0.61 Ma, respectively. Comprehensive analysis shows a relatively singular provenance of the Hala’alat Formation. While the provenance of the Aladeyikesai Formation shows inheritance, the provenance area extends northwards to the Kujibai–Hongguleleng Ophiolitic Mélange Belt. Furthermore, the closure of the Junggar Ocean during Carboniferous time caused the potential source region of the Hala’alat Mountain area to migrate northeastwards from Barleik Mountain to Xiemisitai Mountain. This study provides a basis for the analysis of regional geological evolution.

Fluorsigaiite, ideally Ca2Sr3(PO4)3F, is a new Sr analogue of fluorphosphohedyphane and a new member of the apatite supergroup. It was discovered in lujavrite from the Saima alkaline complex, Liaoning Province, China. Fluorsigaiite commonly occurs as individual prismatic, columnar and platy crystals of 10 to 50 μm in size, associated with microcline, nepheline, aegirine, natrolite, eudialyte, fluorapatite, a fluorstrophite-like mineral, stronadelphite and calcite. Occasionally, crystals of fluorsigaiite form prismatic aggregates in the interstices of lujavrite. Fluorsigaiite is translucent to transparent, colourless to yellowish white with a vitreous lustre and without fluorescence. The estimated Mohs hardness is 5, and the tenacity is brittle with uneven fractures. The calculated density is 3.842 g/cm3. Optically, fluorsigaiite is uniaxial (–) with ω = 1.64(1) and ɛ = 1.63(1) in white light and without dispersion. The mean chemical composition (in wt.%) of fluorsigaiite is Na2O 0.75, CaO 15.17, SrO 44.44, La2O3 3.64, Ce2O3 2.22, Pr2O3 0.19, Nd2O3 0.13, Sm2O3 0.05, Gd2O3 0.23, P2O5 31.87, F 1.91, H2O 0.46, sum 100.26, giving the empirical formula (Sr2.82Ca1.79Na0.16La0.15Ce0.09Pr0.01Nd0.01Gd0.01)Σ5.04P2.97O12[F0.66(OH)0.34]Σ1, which is calculated on the basis of 13 total anions and F+(OH) = 1. The strongest eight lines of its powder X-ray diffraction pattern [d, Å (I, %) (hkl)] are: 3.563 (15) (002), 3.275 (15) (102), 3.144 (19) (120), 2.876 (100) (121), 2.861 (96) (112), 2.772 (27) (300), 1.991 (17) (222) and 1.895 (23) (213). Fluorsigaiite is hexagonal, in the space group P63/m and unit-cell parameters refined from single-crystal X-ray diffraction data are: a = 9.6101(2) Å, c = 7.1311(1) Å, V = 570.35(3) Å3 and Z = 2. It is isostructural with hedyphane-group minerals, and contains different prevailing (species-defining) Ca and Sr cations at the Ca1 and Ca2 sites, respectively. Fluorsigaiite was probably formed from Sr-rich fluids at late-magmatic or hydrothermal stage of the Saima lujavrite.

We perform a theoretical investigation of the instability of a helical vortex filament beneath a free surface in a semi-infinite ideal fluid. The focus is on the leading-order free-surface boundary effect upon the equilibrium form and instability of the vortex. This effect is characterised by the Froude number $F_r = U(gh^*)^{-{1}/{2}}$ where $g$ is gravity, and $U = \varGamma /(2{\rm \pi} b^*)$ with $\varGamma$ being the strength, $2{\rm \pi} b^*$ the pitch and $h^*$ the centre submergence of the helical vortex. In the case of $F_r \rightarrow 0$ corresponding to the presence of a rigid boundary, a new approximate equilibrium form is found if the vortex possesses a non-zero rotational velocity. Compared with the infinite fluid case (Widnall, J. Fluid Mech., vol. 54, no. 4, 1972, pp. 641–663), the vortex is destabilised (or stabilised) to relatively short- (or long-)wavelength sub-harmonic perturbations, but remains stable to super-harmonic perturbations. The wall-boundary effect becomes stronger for smaller helix angle and could dominate over the self-induced flow effect depending on the submergence. In the case of $F_r > 0$, we obtain the surface wave solution induced by the vortex in the context of linearised potential-flow theory. The wave elevation is unbounded when the $m$th wave mode becomes resonant as $F_r$ approaches the critical Froude numbers ${\mathcal {F}} (m) = (C_0^*/U)^{-1} (mh^*/b^*)^{-{1}/{2}}$, $m=1, 2, \ldots,$ where $C_0^*$ is the induced wave speed. We find that the new approximate equilibrium of the vortex exists if and only if $F_r < {\mathcal {F}}(1)$. Compared with the infinite fluid and $F_r \rightarrow 0$ cases, the wave effect causes the vortex to be destabilised to super-harmonic and long-wavelength sub-harmonic perturbations with generally faster growth rate for greater $F_r$ and smaller helix angle.

Dietary modification plays a vital role in the treatment of non-alcoholic liver diseases. We investigated the effects of the consumption of a different amount of dehulled adlay, which has hypolipidaemic and anti-inflammatory properties, on non-alcoholic fatty liver disease (NAFLD). We fed rats a high-fat-high-fructose liquid diet for 16 weeks to induce NAFLD. The rats were divided into three groups fed the NAFLD diet only (NN) or a diet containing 44·9 or 89·8 g/l of dehulled adlay (NA and NB groups, respectively). After 8 weeks, the NA and NB groups had lower C-reactive protein levels and improvement in insulin resistance. In addition, the NB group had lower liver weight and hepatic TAG and cholesterol concentrations than did the NN group. Compared with the NN group, the high-dose NB group had improved steatosis, lower hepatic TNF-α, IL-1β and IL-6 levels and lower adipose leptin levels. Our results suggest that a diet containing dehulled adlay can ameliorate NAFLD progression by decreasing of insulin resistance, steatosis and inflammation.

We perform large-eddy simulations of turbulent flow in a channel constricted by streamwise periodically distributed hill-shaped protrusions. Two Reynolds number cases, i.e. $Re_h=10\,595$ (Fröhlich et al., J. Fluid Mech., vol. 526, 2005, pp. 19–66) and $Re_h=33\,000$ (Kähler et al., J. Fluid Mech., vol. 796, 2016, pp. 257–284), are repeated and utilized to verify and validate our numerical results, including the pressure and skin friction coefficients on bottom and top walls of the channel, mean velocity profiles and Reynolds stresses. All comparisons show reasonable agreement, providing a measure of validity that enables us to further probe simulation results at higher Reynolds number ($Re_h=10^5$) into aspects of flow physics that are not available from experiments. Effects of variation of Reynolds number are studied, with emphasis on the mean skin friction coefficients, separation bubble size and pressure fluctuations that are related to separation and reattachment. In addition, the main large-scale features of the separation behind the hill, including the scaling of the mean velocity profiles, are discussed. Furthermore, the instantaneous near-wall flow field is analysed in terms of skin friction portraits, and we confirm the existence of the local very small separation bubble on the hill crest as observed in experimental and numerical investigations. The flow field at the top wall, which is generally not given sufficient attention, is evaluated with the empirical friction law and universal logarithmic law as in planar channel flows. It is found that these empirical laws compare well with the large-eddy simulation results, although the hill constrictions behave as a perturbation source and the developed shear layer has some effects on the flow field near the top wall.

We describe an ultra-wide-bandwidth, low-frequency receiver recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band ( ${\sim}60\%$ ), the system temperature is approximately 22 K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver, including its astronomical objectives, as well as the feed, receiver, digitiser, and signal processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified, including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration, and timing stability.

Consumption of a high-fat diet increases fat accumulation and may further lead to inflammation and hepatic injuries. The aim of the study was to investigate the effects of Camellia oleifera seed extract (CSE) on non-alcoholic fatty liver disease (NAFLD). After a 16-week NAFLD-inducing period, rats were assigned to experimental groups fed an NAFLD diet with or without CSE. At the end of the study, we found that consuming CSE decreased the abdominal fat weight and hepatic fat accumulation and modulated circulating adipokine levels. We also found that CSE groups had lower hepatic cytochrome P450 2E1 and transforming growth factor (TGF)-β protein expressions. In addition, we found that CSE consumption may have affected the gut microbiota and reduced toll-like receptor (TLR)-4, myeloid differentiation primary response gene 88, toll/IL-1 receptor domain-containing adaptor-inducing interferon-β (TRIF) expression and proinflammatory cytokine concentrations in the liver. Our results suggest that CSE may alleviate the progression of NAFLD in rats with diet-induced steatosis through reducing fat accumulation and improving lipid metabolism and hepatic inflammation.