This article presents a data-driven model based on modal decomposition, applied to approximate the low-order statistics of the spatially averaged wall-shear stress in a turbulent channel flow over a porous wall with two anisotropic permeabilities, producing drag increase or reduction when compared with the case of an isotropic porous wall. The model is comparable to a neural network architecture using a linear map to a classification. To create this model, we use high-order dynamic mode decomposition (DMD) to identify the structures describing the main flow dynamics, and then test different linear combinations of these modes to estimate the time evolution of the stress at the porous interface. The coefficients of the model are obtained by training the model against the results of direct numerical simulations over different time intervals. Depending on the number and the way of combining the DMD modes, the reduced-order models presented can reconstruct the wall-shear stress with relative error smaller than 0.01 % and reproduce its statistical variations for at least 1500 time units with relative error in the standard deviation or the mean smaller than 5 %. The model has also been tested to approximate the statistics of the wall-shear stress over the whole wall, showing that the regeneration of the flow structures can be reproduced by the nonlinear interaction of modes. Finally, considering the DMD modes as communities in a neural network, we examine the influence of the mode-to-mode interaction on the nonlinear flow dynamics, which explains the performance of the different models.

We analysed both structural and functional aspects of sentences containing the four adverbials “after”, “before”, “because”, and “if” in two dense corpora of parent-child interactions from two British English-acquiring children (2;00–4;07). In comparing mothers’ and children's usage we separate out the effects of frequency, cognitive complexity and pragmatics in explaining the course of acquisition of adverbial sentences. We also compare these usage patterns to stimuli used in a range of experimental studies and show how differences may account for some of the difficulties that children have shown in experiments. In addition, we report descriptive data on various aspects of adverbial sentences that have not yet been studied as a resource for future investigations.

In this numerical and theoretical work, we study the turbulent channel flow of Newtonian and elastoviscoplastic fluids. The coherent structures in these flows are identified by means of higher order dynamic mode decomposition (HODMD), applied to a set of data non-equidistant in time, to reveal the role of the near-wall streaks and their breakdown, and the interplay between turbulent dynamics and non-Newtonian effects. HODMD identifies six different high-amplitude modes, which either describe the yielded flow or the yielded–unyielded flow interaction. The structure of the low- and high-frequency modes suggests that the interaction between high- and low-speed streamwise velocity structures is one of the mechanisms triggering the streak breakdown, dominant in Newtonian turbulence where we observe shorter near-wall streaks and a more chaotic dynamics. As the influence of elasticity and plasticity increases, the flow becomes more correlated in the streamwise direction, with long streaks disrupted for short times by localised perturbations, reflected in reduced drag. Finally, we present streamwise-periodic dynamic mode decomposition modes as a viable tool to describe the highly complex turbulent flows, and identify simple well-organised groups of travelling waves.

The Empire of Aksum was one of Africa's most influential ancient civilisations. Traditionally, most archaeological fieldwork has focused on the capital city of Aksum, but recent research at the site of Beta Samati has investigated a contemporaneous trade and religious centre located between Aksum and the Red Sea. The authors outline the discovery of the site and present important finds from the initial excavations, including an early basilica, inscriptions and a gold intaglio ring. From daily life and ritual praxis to international trade, this work illuminates the role of Beta Samati as an administrative centre and its significance within the wider Aksumite world.

The COllaborative project of Development of Anthropometrical measures in Twins (CODATwins) project is a large international collaborative effort to analyze individual-level phenotype data from twins in multiple cohorts from different environments. The main objective is to study factors that modify genetic and environmental variation of height, body mass index (BMI, kg/m2) and size at birth, and additionally to address other research questions such as long-term consequences of birth size. The project started in 2013 and is open to all twin projects in the world having height and weight measures on twins with information on zygosity. Thus far, 54 twin projects from 24 countries have provided individual-level data. The CODATwins database includes 489,981 twin individuals (228,635 complete twin pairs). Since many twin cohorts have collected longitudinal data, there is a total of 1,049,785 height and weight observations. For many cohorts, we also have information on birth weight and length, own smoking behavior and own or parental education. We found that the heritability estimates of height and BMI systematically changed from infancy to old age. Remarkably, only minor differences in the heritability estimates were found across cultural–geographic regions, measurement time and birth cohort for height and BMI. In addition to genetic epidemiological studies, we looked at associations of height and BMI with education, birth weight and smoking status. Within-family analyses examined differences within same-sex and opposite-sex dizygotic twins in birth size and later development. The CODATwins project demonstrates the feasibility and value of international collaboration to address gene-by-exposure interactions that require large sample sizes and address the effects of different exposures across time, geographical regions and socioeconomic status.

We study single-phase and particle-laden turbulent channel flows bounded by two incompressible hyper-elastic walls with different deformability at bulk Reynolds number $5600$. The solid volume fraction of finite-size neutrally buoyant rigid spherical particles considered is $10\,\%$. The elastic walls are assumed to be of a neo-Hookean material. A fully Eulerian formulation is employed to model the elastic walls together with a direct-forcing immersed boundary method for the coupling between the fluid and the particles. The data show a significant drag increase and the enhancement of the turbulence activity with growing wall elasticity for both the single-phase and particle-laden flows when compared with the single-phase flow over rigid walls. Drag reduction and turbulence attenuation is obtained, on the other hand, with highly elastic walls when comparing the particle-laden flow with the single-phase flow for the same wall properties; the opposite effect, drag increase, is observed upon adding particles to the flow over less elastic walls. This is explained by investigating the near-wall turbulence, where the strong asymmetry in the magnitude of the wall-normal velocity fluctuations (favouring positive $v^{\prime }$), is found to push the particles towards the channel centre. The particle layer close to the wall contributes to turbulence production by increasing the wall-normal velocity fluctuations, so that in the absence of this layer, smaller wall deformations and in turn turbulence attenuation is observed. For a moderate wall elasticity, we increase the particle volume fraction up to $20\,\%$ and find that particle migration away from the wall is the cause of turbulence attenuation with respect to the flow over rigid walls. However, for this higher volume fractions, the particle induced stress compensates for the decreasing Reynolds shear stress, resulting in a higher overall drag for the case with elastic walls. The effect of the wall elasticity on the overall drag reduces significantly with increasing particle volume fraction.

We present a study of Lagrangian intermittency and its characteristic time scales. Using the concepts of flying and diving residence times above and below a given threshold in the magnitude of turbulence quantities, we infer the time spectra of the Lagrangian temporal fluctuations of dissipation, acceleration and enstrophy by means of a direct numerical simulation in homogeneous and isotropic turbulence. We then relate these time scales, first, to the presence of extreme events in turbulence and, second, to the local flow characteristics. Analyses confirm the existence in turbulent quantities of holes mirroring bursts, both of which are at the core of what constitutes Lagrangian intermittency. It is shown that holes are associated with quiescent laminar regions of the flow. Moreover, Lagrangian holes occur over few Kolmogorov time scales while Lagrangian bursts happen over longer periods scaling with the global decorrelation time scale, hence showing that loss of the history of the turbulence quantities along particle trajectories in turbulence is not continuous. Such a characteristic partially explains why current Lagrangian stochastic models fail at reproducing our results. More generally, the Lagrangian dataset of residence times shown here represents another manner for qualifying the accuracy of models. We also deliver a theoretical approximation of mean residence times, which highlights the importance of the correlation between turbulence quantities and their time derivatives in setting temporal statistics. Finally, whether in a hole or a burst, the straining structure along particle trajectories always evolves self-similarly (in a statistical sense) from shearless two-dimensional to shear bi-axial configurations. We speculate that this latter configuration represents the optimum manner to dissipate locally the available energy.

Finite-size particles modulate wall-bounded turbulence, leading, for the case of spherical particles, to increased drag also owing to the formation of a particle wall layer. Here, we study the effect of particle shape on the turbulence in suspensions of spheroidal particles at volume fraction $\unicode[STIX]{x1D719}=10\,\%$ and show how the near-wall particle dynamics deeply changes with the particle aspect ratio and how this affects the global suspension behaviour. Direct numerical simulations are performed using a direct-forcing immersed boundary method to account for the dispersed phase, combined with a soft-sphere collision model and lubrication corrections for short-range particle–particle and particle–wall interactions. The turbulence reduces with the aspect ratio of oblate particles, leading to drag reduction with respect to the single-phase flow for particles with aspect ratio ${\mathcal{A}}{\mathcal{R}}\leqslant 1/3$, when the significant reduction in Reynolds shear stress is more than the compensation by the additional stresses, induced by the solid phase. Oblate particles are found to avoid the region close to the wall, travelling parallel to it with small angular velocities, while preferentially sampling high-speed fluid in the wall region. Prolate particles also tend to orient parallel to the wall and avoid its vicinity. Their reluctance to rotate around the spanwise axis reduces the wall-normal velocity fluctuation of the flow and therefore the turbulence Reynolds stress, similar to oblates; however, they undergo rotations in wall-parallel planes which increase the additional solid stresses due to their relatively larger angular velocities compared to the oblates. These larger additional stresses compensate for the reduction in turbulence activity and lead to a wall drag similar to that of single-phase flows. Spheres on the other hand, form a layer close to the wall with large angular velocities in the spanwise direction, which increases the turbulence activity in addition to exerting the largest solid stresses on the suspension, in comparison to the other studied shapes. Spherical particles therefore increase the wall drag with respect to the single-phase flow.

We study heat transfer in plane Couette flow laden with rigid spherical particles by means of direct numerical simulations. In the simulations we use a direct-forcing immersed boundary method to account for the dispersed phase together with a volume-of-fluid approach to solve the temperature field inside and outside the particles. We focus on the variation of the heat transfer with the particle Reynolds number, total volume fraction (number of particles) and the ratio between the particle and fluid thermal diffusivity, quantified in terms of an effective suspension diffusivity. We show that, when inertia at the particle scale is negligible, the heat transfer increases with respect to the unladen case following an empirical correlation recently proposed in the literature. In addition, an average composite diffusivity can be used to approximate the effective diffusivity of the suspension in the inertialess regime when varying the molecular diffusion in the two phases. At finite particle inertia, however, the heat transfer increase is significantly larger, smoothly saturating at higher volume fractions. By phase-ensemble-averaging we identify the different mechanisms contributing to the total heat transfer and show that the increase of the effective conductivity observed at finite inertia is due to the increase of the transport associated with fluid and particle velocity. We also show that the contribution of the heat conduction in the solid phase to the total wall-normal heat flux reduces when increasing the particle Reynolds number, so that particles of low thermal diffusivity weakly alter the total heat flux in the suspension at finite particle Reynolds numbers. On the other hand, a higher particle thermal diffusivity significantly increases the total heat transfer.

Phytoplankton are the foundation of aquatic food webs. Through photosynthesis, phytoplankton draw down $\text{CO}_{2}$ at magnitudes equivalent to forests and other terrestrial plants and convert it to organic material that is then consumed by other planktonic organisms in higher trophic levels. Mechanisms that affect local concentrations and velocities are of primary significance to many encounter-based processes in the plankton, including prey–predator interactions, fertilization and aggregate formation. We report results from simulations of sinking phytoplankton, considered as elongated spheroids, in homogenous isotropic turbulence to answer the question of whether trajectories and velocities of sinking phytoplankton are altered by turbulence. We show in particular that settling spheroids with physical characteristics similar to those of diatoms weakly cluster and preferentially sample regions of downwelling flow, corresponding to an increase of the mean settling speed with respect to the mean settling speed in quiescent fluid. We explain how different parameters can affect the settling speed and what underlying mechanisms might be involved. Interestingly, we observe that the increase in the aspect ratio of the prolate spheroids can affect the clustering and the average settling speed of particles by two mechanisms: first is the effect of aspect ratio on the rotation rate of the particles, which saturates faster than the second mechanism of increasing drag anisotropy.