On the analysis of sloshing of water in rectangular containers: Numerical study and experimental validation

In this work the analysis of sloshing of water in rectangular open tanks has been extensively carried out. Two mathematical models are employed, respectively the Reynolds Averaged Navier Stokes Equations (RANSE) and the Shallow Water Equations (SWE). The RANSE are solved using a modified form of the well established MAC method (SIMAC) able to treat both the free surface motion and the viscous stresses over the rigid walls accurately. The Shallow Water Equations are solved by means of a simple and powerful algorithm (CE-SE) able to deal with large impacting waves over the tank walls.Successively, in order to validate the mentioned algorithms and for a better understanding of the sloshing phenomenon, experimental tests have been carried out using a 0.5 m breadth rectangular tank in periodic roll motion.It has been shown that RANSE provide more accurate solutions than SWE for small or moderate amplitudes of excitation. In particular in this paper it is proved that the shallow water approximation can be efficiently adopted within liquid depth to tank breadth RATIO = 0.15, when examining the sloshing problem. By increasing the water level inside the tank, results by SWE show large qualitative and quantitative disagreement with experiments. Nevertheless, in the case of large amplitude excitation, when sprays and large breaking waves are expected, SWE provide a fairly good estimate of the sloshing induced waves.Finally a simple baffle configuration inside the tank has been considered. By the analysis of numerical results, it has been observed that the presence of a vertical baffle at the middle of the tank dramatically changes the sloshing response compared to the unbaffled configuration. It produces a jump-like effect, resulting in a weak magnification of the dynamic loads on the vertical walls out of resonance, and a strong reduction of the dynamic loads in the resonance condition.     

Influence of topsoil of pyroclastic origin on microbial contamination of groundwater in fractured carbonate aquifers

The aim of the research was to analyse the influence of a topsoil of pyroclastic origin on microbial contamination of groundwater in a carbonate aquifer and verify the reliability of thermotolerant coliforms and fecal enterococci as bacterial indicators. The research was carried out through hydrogeological and microbiological monitoring at an experimental field site in Italy during two hydrologic years and through column tests in a laboratory. The taxonomic classification of fecal indicators detected in spring water samples was performed using API20 galleries. Fecal enterococci were also identified by means of 16S rRNA gene sequencing. The topsoil of pyroclastic origin significantly retains both thermotolerant coliforms and fecal enterococci. Results of column tests carried out in soil blocks collected randomly within the test site suggest that Escherichia coli was more retained than Enterococcus faecalis, even though this difference is statistically significant in only two out of six soil samples. Thus, a non-uniform difference in retention is expected at field scale. This suggestion is in agreement with the results of the microbiological monitoring. In fact, fecal enterococci were a more reliable indicator than thermotolerant coliforms for detecting contamination at both seasonal springs of the aquifer system, while no significant differences were observed at the perennial spring.     

Petrology of Karoo volcanic rocks in the southern Lebombo monocline, Mozambique

The Karoo volcanic sequence in the southern Lebombo monocline in Mozambique contains different silicic units in the form of pyroclastic rocks, and two different basalt types. The silicic units in the lower part of the Lebombo sequence are formed by a lower unit of dacites and rhyolites (67–80 wt.% SiO₂) with high Ba (990–2500 ppm), Zr (800–1100 ppm) and Y (130–240 ppm), which are part of the Jozini–Mbuluzi Formation, followed by a second unit, interlayered with the Movene basalts, of high-SiO₂ rhyolites (76–78 wt.%; the Sica Beds Formation), with low Sr (19–54 ppm), Zr (340–480 ppm) and Ba (330–850 ppm) plus rare quartz-trachytes (64–66 wt.% SiO₂), with high Nb and Rb contents (240–250 and 370–381 ppm, respectively), and relatively low Zr (450–460 ppm). The mafic rocks found at the top of the sequence are basalts and ferrobasalts belonging to the Movene Formation. The basalts have roughly flat mantle-normalized incompatible element patterns, with abundances of the most incompatible elements not higher than 25 times primitive mantle. The ferrobasalt has TiO₂  4.7 wt.%, Fe₂O_3t = 16 wt.%, and high Y (100 ppm), Zr (420 ppm) and Ba (1000 ppm). The Movene basalts have initial (at 180 Ma) ⁸⁷Sr/⁸⁶Sr = 0.7052–0.7054 and ¹⁴³Nd/¹⁴⁴Nd = 0.51232, and the Movene ferrobasalt has even lower ⁸⁷Sr/⁸⁶Sr (0.70377) and higher ¹⁴³Nd/¹⁴⁴Nd (0.51259). The silicic rocks show a modest range of initial Sr-(⁸⁷Sr/⁸⁶Sr = 0.70470–0.70648) and Nd-(¹⁴³Nd/¹⁴⁴Nd = 0.51223–0.51243) isotope ratios. The less evolved dacites could have been formed after crystal fractionation of oxide-rich gabbroic cumulates from mafic parental magmas, whereas the most silica-rich rhyolites could have been formed after fractional crystallization of feldspars, pyroxenes, oxides, zircon and apatite from a parental dacite magma. The composition of the Movene basalts imply different feeding systems from those of the underlying Sabie River basalts.     

Unidirectional cyclic resistance of Ticino and Toyoura sands from centrifuge cone penetration tests

The evaluation of the undrained cyclic resistance of sandy deposits is required to forecast the soil behaviour during an earthquake (liquefaction, cyclic mobility); due to the difficulties in obtaining undisturbed samples of most liquefiable soils, it is usually deduced from field test results such as cone penetration tests. This paper proposes a methodology to evaluate the undrained cyclic resistance from normalised cone resistance of two well-studied silica sands (Ticino and Toyoura), with different mineralogy, one mainly composed of feldspar, the other of quartz. The determination of the cyclic resistance of Ticino and Toyoura sands was achieved through undrained cyclic triaxial tests on reconstituted specimens. The tip resistance was deduced from CPTs performed in centrifuge with a miniaturised piezocone on homogeneous reconstituted models. Both the undrained cyclic and tip resistances were correlated with the state parameter ψ. Results of centrifuge and triaxial tests were combined through ψ to deduce the cyclic resistance ratio CRR directly from the normalised cone resistance. The shape of the curve relating CRR to the normalised cone resistance resulted unusual respect to all the recognised curves widespread in the geotechnical literature. The aim of the proposed correlations is to provide a useful instrument to improve the actual knowledge on liquefaction and to give a contribution based on the critical state soil mechanics framework to the development of refined correlations between the cyclic resistance of a sand and the results of cone penetration tests.     

A new version of the USLE‐MM for predicting bare plot soil loss at the Sparacia (South Italy) experimental site

Improving empirical prediction of plot soil erosion at the event temporal scale has both scientific and practical importance. In this investigation, 492 runoff and soil loss data from plots of different lengths, λ (11 ≤ λ ≤ 44 m), and steepness, s (14.9 ≤ s ≤ 26.0%), established at the Sparacia experimental station, in Sicily, South Italy, were used to derive a new version of Universal Soil Loss Equation (USLE)‐MM model, by only assuming a value of one for the topographic length, L, and steepness, S, factors for λ = 22 m and s = 9%, respectively. An erosivity index equal to (Q_REI₃₀)^b1, Q_R and EI₃₀ being the runoff coefficient and the event rainfall erosivity index, respectively, with b₁ > 1 was found to be an appropriate choice for the Sparacia area. The specifically developed functions for L and S did not differ appreciably from other, more widely accepted relationships (maximum differences by a factor of 1.22 for L and 1.09 for S). The new version of the USLE‐MM performed particularly well for highly erosive events, because predicted soil loss differed by not more than a factor of 1.19 from the measured soil loss for measured values of more than 100 Mg ha^?1. The choice of the relationships to predict topographic effects on plot soil loss should not represent a point of particular concern in the application of the USLE‐MM in other environments. However, tests of the empirical approach should be carried out in other experimental areas in an attempt to develop analytical tools, usable at the event temporal scale, reasonably simple and of wide validity. Copyright © 2015 John Wiley & Sons, Ltd.     

Measuring rill erosion at plot scale by a drone‐based technology

The traditional direct method (i.e. metric ruler and rillmeter) of monitoring rill erosion at plot scale is time consuming and invasive because it modifies the surface of the rilled area. Measuring rill features using a drone‐based technology is considered a non‐invasive method allowing a fast field relief. In the experimental Sparacia area a survey by a quadricopter Microdrones md4‐200 was carried out, and this relief allowed the generation of a Digital Elevation Model (DEM), with a mesh size of 1 cm and a resolution elevation equal to 2 mm, for three plots (L, G and C) affected by rill erosion. At first for the experimental L plot, which is 44 m long, the rill features were surveyed by a ‘manual’ method which was carried out by drawing on the PC screen the rill paths obtained by a visual orthophoto interpretation. This manual method was not applicable for the plots in which rills of limited depth occurred and were not detectable. Then, for both L plot and the other experimental plots having a length ranging from 22 to 44 m, an ‘automatic’ extraction method of rills from DEM was applied. Using an appropriate calculation routine, a vector coverage of transects orthogonal to the main flow direction (i.e. the maximum slope steepness path) was generated. The intersection of each plot DEM with the transect coverage allowed to obtain both the cross sections and the main rill morphological features. For the L plot the comparison between the rill morphological features obtained by the two different methods (manual, automatic) was carried out. Finally, the length–volume relationship and a dimensionless relationship proposed in literature were tested for all studied experimental plots. Copyright © 2015 John Wiley & Sons, Ltd.     

A shell model for anisotropic magnetohydrodynamic turbulence

Abstract

In this paper, starting from the spectral DIA equations obtained by Veltri et al. (1982), describing the spectral dynamical evolution of magnetohydrodynamic (MHD) turbulence in the presence of a background magnetic field B ₀, we have derived an approximate form of these equations (shell model) more appropriate for numerical integration at high Reynolds numbers.

We have studied the decay of an initially isotropic state, with an initial imbalance between the energies for the two signs of the cross-helicity. Reynolds numbers up to 10⁵ have been considered.

Numerical results show that the nonlinear energy cascade behaves anisotropically in the k-space, i.e. in the spectra there is a prevalence of the wavevectors perpendicular to B ₀ with respect to the parallel wavevectors. This anisotropic effect, which is due to the presence of the background magnetic field, can be understood in terms of the so-called ‘‘Alfvén effect''.

A different source of anisotropy, due to the difference of the energy transfer for the two polarizations perpendicular to k, is recovered, but its effect is found to be mainly concentrated in the injection range.

Only little differences have been found, in the inertial range, in the spectral indices from the Kraichnan 3/2 value, which is valid for an isotropic spectrum. A form for the anisotropic spectrum can be recovered phenomenologically from our results. Values of the spectral indices quite different from the Kraichnan 3 2 value are obtained only when we consider stationary states with different forcing terms for the two modes of Alfvén wave propagation.

The comparison of our results with the observations of the v and B fluctuations in the interplanatery space shows that the anisotropy found in interplanetary fluctuations might be attributed only partially to the result of a nonlinear energy cascade.