Effects of grain morphology on critical state: a computational analysis

We introduce a new DEM scheme (LS-DEM) that takes advantage of level sets to enable the inclusion of real grain shapes into a classical discrete element method. Then, LS-DEM is validated and calibrated with respect to real experimental results. Finally, we exploit part of LS-DEM potentiality by using it to study the dependency of critical state (CS) parameters such as critical state line (CSL) slope \(\lambda \), CSL intercept \(\varGamma \), and CS friction angle \(\varPhi _{\mathrm{CS}}\) on the grain’s morphology, i.e., sphericity, roundness, and regularity. This study is carried out in three steps. First, LS-DEM is used to capture and simulate the shape of five different two-dimensional cross sections of real grains, which have been previously classified according to the aforementioned morphological features. Second, the same LS-DEM simulations are carried out for idealized/simplified grains, which are morphologically equivalent to their real counterparts. Third, the results of real and idealized grains are compared, so the effect of “imperfections” on real particles is isolated. Finally, trends for the CS parameters (CSP) dependency on sphericity, roundness, and regularity are obtained as well as analyzed. The main observations and remarks connecting particle’s morphology, particle’s idealization, and CSP are summarized in a table that is attempted to help in keeping a general picture of the analysis, results, and corresponding implications.     

Strength criterion for cross-anisotropic sand under general stress conditions

By incorporating the fabric effect and Lode’s angle dependence into the Mohr–Coulomb failure criterion, a strength criterion for cross-anisotropic sand under general stress conditions was proposed. The obtained criterion has only three material parameters which can be specified by conventional triaxial tests. The formula to calculate the friction angle under any loading direction and intermediate principal stress ratio condition was deduced, and the influence of the degree of the cross-anisotropy was quantified. The friction angles of sand in triaxial, true triaxial, and hollow cylinder torsional shear tests were obtained, and a parametric analysis was used to detect the varying characteristics. The friction angle becomes smaller when the major principal stress changes from perpendicular to parallel to the bedding plane. The loading direction and intermediate principal stress ratio are unrelated in true triaxial tests, and their influences on the friction angle can be well captured by the proposed criterion. In hollow cylinder torsional shear tests with the same internal and external pressures, the loading direction and intermediate principal stress ratio are related. This property results in a lower friction angle in the hollow cylinder torsional shear test than that in the true triaxial test under the same intermediate principal stress ratio condition. By comparing the calculated friction angle with the experimental results under various loading conditions (e.g., triaxial, true triaxial, and hollow cylinder torsional shear test), the proposed criterion was verified to be able to characterize the shear strength of cross-anisotropic sand under general stress conditions.     

Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

Gravimetric and magnetic fabric study of the Sintra Igneous complex: laccolith-plug emplacement in the Western Iberian passive margin

International Journal of Earth Sciences - The geometry and emplacement of the ~ 96 km2, Late Cretaceous Sintra Igneous complex (SIC, ca. 80 Ma) into the West Iberian passive...     

The Juchatengo complex: an upper-level ophiolite assemblage of late Paleozoic age in Oaxaca,southern Mexico

International Journal of Earth Sciences - The Juchatengo complex (JC) suite is located between the Proterozoic Oaxacan complex to the north and the Xolapa complex to the south, and was amalgamated...     

Rainfall thresholds for the triggering of landslides considering previous soil moisture conditions (Asturias,NW Spain)

Given its geological and climatic conditions and its rugged orography, Asturias is one of the most landslide prone areas in the North of Spain. Most of the landslides occur during intense rainfall episodes. Thus, precipitation is considered the main triggering factor in the study area, reaching average annual values of 960 mm. Two main precipitation patterns are frequent: (i) long-lasting periods of moderate rainfall during autumn and winter and (ii) heavy short rainfall episodes during spring and early summer. In the present work, soil moisture conditions in the locations of 84 landslides are analysed during two rainfall episodes, which represent the most common precipitation patterns: October–November 2008 and June 2010. Empirical data allowed the definition of available water capacity percentages of 99–100% as critical soil moisture conditions for the landslide triggering. Intensity-duration rainfall thresholds were calculated for each episode, considering the periods with sustained high soil moisture levels before the occurrence of each analysed landslide event. For this purpose, data from daily water balance models and weather stations were used. An inverse relationship between the duration of the precipitation and its intensity, consistent with published intensity-duration thresholds, was observed, showing relevant seasonal differences.     

Analysis of physical factors of climate and bioclimate and their effects on almonds production to increase the economy in Hebron area of Palestine

For the 1993–2009 period, we analyzed the relationship between almond yield and three climatic variables (mean annual temperature, soil water reserve, and precipitation), and four bioclimatic variables (annual ombrothermic index, water deficit, simple continentality index, and compensated thermicity index), for one major Hebron crop (soft and hard almonds). Moreover, we obtained data almond production from the Palestinian Central Bureau of Statistics, while the climate data from the Palestinian meteorological station during the study period from 1993 to 2009, and analysis is it by using bioclimatic classification of the Earth of Salvador Rivas-Martinez to study the relationship between the almond yield and climate and bioclimate factors (variables). The climatic and bioclimate variables of greatest importance to almond were used to develop regressions analysis relating yield to climatic conditions. Hebron was positively affected by annual ombrothermic index, simple continentality index, precipitation, water soil reserve, and mean annual temperature, but negatively affected by water deficit, with a large proportion of the variance explained by axis F1 (72.48%), F2 (22.38%), and axes F1and F2 (94. 86%). However, in order to produce a high amount of almonds and quality, it can be grown in the regions of the mesomediterranean region, with the value of annual ombrothemic index more than 3, compensated thermicity index between 220/220 to 350/350, simple continentality index between 14 and 20, and in areas where the average annual temperature is between 15 and 20 °C.