Magnetic properties and cation ordering of nanopowders of the synthetic analogue of kuramite, Cu3SnS4

An extensive characterisation of the magnetic properties of synthetic powders of kuramite, with formal composition Cu₃SnS₄, was performed. Powders were investigated through superconducting quantum interference device (SQUID) magnetometry, electron paramagnetic resonance (EPR) spectroscopy, X-ray powder diffraction (XRPD), scanning and transmission electron microscopies (SEM and TEM) and microanalysis. SEM and TEM reveal the presence of nanodimensioned particles. XRPD clearly shows that Cu₃SnS₄ crystallised in a cubic sphalerite-type structural model, in spite of the stannite-type tetragonal structure described for the natural phase. This difference arises from a full random distribution of cations. Synthetic kuramite nanopowders exhibit a marked paramagnetism, originated by the presence of Cu(II), definitely assessed by EPR measurements. Moreover, the overall magnetic behaviour of the sample cannot be simply ascribed to diluted paramagnetism, and this suggests the presence of strong superexchange interactions among Cu(II) ions even at room temperature. The main consequences of these results are the definitive assessment of the chemical formula Cu(I)₂Cu(II)SnS₄ and of a random distribution of Cu(II), Cu(I) and Sn(IV) ions within the available tetrahedral sites.     

Spatial modelling and uncertainty assessment of pyroclastic cover thickness in the Sorrento Peninsula

Assessing spatial variability of soil thickness is a critical issue for understanding and predicting slope processes. The present work was aimed at estimating the spatial scales at which the variation of pyroclastic cover thickness occurs in a sample area in the Sorrento Peninsula (Italy). Stochastic simulation was used to understand the spatial variability of pyroclastic cover thickness on Mount Pendolo and to assess its spatial uncertainty. In the study area, covering about 0.7 km², thickness measurements were collected using electrical resistivity tomography profiles, continuous core drillings and steel rod penetrometric tests. Variographic analysis revealed the occurrence of an anisotropic behaviour along the N50 and N140 directions. In the latter anisotropic direction, a nested variogram was fitted including (1) a long-range component which could be related to large-scale factors, like the curvature of the slope and contributing area and (2) a shorter scale variation which is probably associated with the occurrence of denudation processes or to the articulate cover/bedrock interface. To assess the spatial variability and uncertainty of pyroclastic cover thickness, a stochastic simulation algorithm was used and 500 equally probable images of cover thickness were yielded. The results showed that a better thickness distribution map can be drawn by simulating the data collected on the slope and at the footslope separately. The approach also allowed delineating the areas characterized by greater uncertainty, suggesting supplementary measurements to further improve the cover thickness distribution model, thus reducing the uncertainty.     

Numerical homogenization-based seismic assessment of an English-bond masonry prototype: Structural level application

A multiscale strategy is evaluated at a structural level for the analysis of unreinforced masonry structures. The mechanical characterization of the masonry is deduced from homogenization-based micro-scale finite element (FE) models. The derived data are here employed at a structural level via a discrete FE model. The discrete FE model is composed of quadrilateral rigid plates interconnected through vertical and horizontal interfaces. On the interfaces, between adjoining discrete elements, a model that accounts for the in- and out-of-plane behavior of masonry, with damage and plasticity, is adopted. Such interfaces represent the material pre- and post-peak regimes, its orthotropy, and, depending on the micro-model assumed, account by three-dimensional shear effects that are especially important for multi-leaf walls and complex regular textures. The discrete model has been implemented in an advanced structural analysis software where powerful built-in features as the arc-length method, line-search algorithm, and implicit or explicit solver schemes are available. The multi-scale model is applied for the dynamic study of a small English-bond masonry house prototype subjected to a series of consecutive earthquake records. Detailed comparisons between the experimental and numerical data are presented, including the results obtained through a continuous total strain rotating crack model. Quasi-static and dynamic analyses are conducted. Results demonstrate that when enough experimental information is available on the masonry components under tension, shear, and compression regimes, the approach predicts well the seismic structural response in terms of time-history displacements, seismic capacity, and damage patterns. The required computational cost (CPU time) is very attractive.     

Ensemble modelling of ice cover for a reservoir affected by pumped‐storage operation and climate change

Ensemble modelling was used to assess the robustness of projected impacts of pumped‐storage (PS) operation and climate change on reservoir ice cover. To this end, three one‐dimensional and a two‐dimensional laterally averaged hydrodynamic model were set up. For the latter, the strength of the impacts with increasing distance from the dam was also investigated. Climate change effects were simulated by forcing the models with 150 years of synthetic meteorological time series created with a weather generator based on available air temperature scenarios for Switzerland. Future climate by the end of the 21^st century was projected to shorten the ice‐covered period by ~2 months and decrease ice thicknesses by ~13 cm. Under current climate conditions, the ice cover would already be affected by extended PS operation. For example, the average probability of ice coverage on a specific day was projected to decrease by ~13% for current climate and could further be reduced from ~45% to ~10% for future climate. Overall, the results of all models were consistent. Although the number of winters without ice cover was projected to increase for all one‐dimensional models, studying individual segments of the two‐dimensional model showed that the impact was pronounced for segments close to the PS intake/outlet. In summary, the reservoir's ice cover is expected to partially vanish with higher probability of open water conditions closer to the PS intake/outlet.     

Capture into first-order resonances and long-term stability of pairs of equal-mass planets

Massive planets form within the lifetime of protoplanetary disks, and therefore, they are subject to orbital migration due to planet–disk interactions. When the first planet reaches the inner edge of the disk, its migration stops and consequently the second planet ends up locked in resonance with the first one. We detail how the resonant trapping works comparing semi-analytical formulae and numerical simulations. We restrict to the case of two equal-mass coplanar planets trapped in first-order resonances, but the method can be easily generalized. We first describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using series expansions up to different orders in eccentricity as well as a non-expanded Hamiltonian. Then we show that during convergent migration the planets evolve along these families of equilibrium points. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict precisely analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that in most cases the resonant configurations achieved by migration become unstable after the removal of the protoplanetary disk. Here we probe the stability of the resonances as a function of planetary mass. For this purpose, we fictitiously increase the masses of resonant planets, adiabatically maintaining the low-amplitude libration regime until instability occurs. We discuss two hypotheses for the instability, that of a low-order secondary resonance of the libration frequency with a fast synodic frequency of the system, and that of minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter nevertheless play the decisive role in the destabilization of resonant pairs. We show evidence that as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.     

Upscaling permeability for three-dimensional fractured porous rocks with the multiple boundary method

Upscaling permeability of grid blocks is crucial for groundwater models. A novel upscaling method for three-dimensional fractured porous rocks is presented. The objective of the study was to compare this method with the commonly used Oda upscaling method and the volume averaging method. First, the multiple boundary method and its computational framework were defined for three-dimensional stochastic fracture networks. Then, the different upscaling methods were compared for a set of rotated fractures, for tortuous fractures, and for two discrete fracture networks. The results computed by the multiple boundary method are comparable with those of the other two methods and fit best the analytical solution for a set of rotated fractures. The errors in flow rate of the equivalent fracture model decrease when using the multiple boundary method. Furthermore, the errors of the equivalent fracture models increase from well-connected fracture networks to poorly connected ones. Finally, the diagonal components of the equivalent permeability tensors tend to follow a normal or log-normal distribution for the well-connected fracture network model with infinite fracture size. By contrast, they exhibit a power-law distribution for the poorly connected fracture network with multiple scale fractures. The study demonstrates the accuracy and the flexibility of the multiple boundary upscaling concept. This makes it attractive for being incorporated into any existing flow-based upscaling procedures, which helps in reducing the uncertainty of groundwater models.