Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Global height datum unification: a new approach in gravity potential space

The problem of “global height datum unification” is solved in the gravity potential space based on: (1) high-resolution local gravity field modeling, (2) geocentric coordinates of the reference benchmark, and (3) a known value of the geoid’s potential. The high-resolution local gravity field model is derived based on a solution of the fixed-free two-boundary-value problem of the Earth’s gravity field using (a) potential difference values (from precise leveling), (b) modulus of the gravity vector (from gravimetry), (c) astronomical longitude and latitude (from geodetic astronomy and/or combination of (GNSS) Global Navigation Satellite System observations with total station measurements), (d) and satellite altimetry. Knowing the height of the reference benchmark in the national height system and its geocentric GNSS coordinates, and using the derived high-resolution local gravity field model, the gravity potential value of the zero point of the height system is computed. The difference between the derived gravity potential value of the zero point of the height system and the geoid’s potential value is computed. This potential difference gives the offset of the zero point of the height system from geoid in the “potential space”, which is transferred into “geometry space” using the transformation formula derived in this paper. The method was applied to the computation of the offset of the zero point of the Iranian height datum from the geoid’s potential value W ₀=62636855.8 m²/s². According to the geometry space computations, the height datum of Iran is 0.09 m below the geoid.     

A new GIS-based data mining technique using an adaptive neuro-fuzzy inference system (ANFIS) and k-fold cross-validation approach for land subsidence susceptibility mapping

In this paper, we evaluate the predictive performance of an adaptive neuro-fuzzy inference system (ANFIS) using six different membership functions (MF). In combination with a geographic information system (GIS), ANFIS was used for land subsidence susceptibility mapping (LSSM) in the Marand plain, northwest Iran. This area is prone to droughts and low groundwater levels and subsequent land subsidence damages. Therefore, a land subsidence inventory database was created from an extensive field survey. Areas of land subsidence or areas showing initial signs of subsidence were used for training, while one-third of inventory database were reserved for testing and validation. The inventory database randomly divided into three different folds of the same size. One of the folds was chosen for testing and validation. Other two folds was used for training. This process repeated for every fold in the inventory dataset. Thereafter, land subsidence related factors, such as hydrological and topographical factors, were prepared as GIS layers. Areas susceptible to land subsidence were then analyzed using the ANFIS approach, and land subsidence susceptibility maps were created, whereby six different MFs were applied. Lastly, the results derived from each MF were validated with those areas of the land subsidence database that were not used for training. Receiver operating characteristics (ROC) curves were drawn for all LSSMs, and the areas under the curves were calculated. The ROC analyses for the six LSSMs yielded very high prediction values for two out of the six methods, namely the difference of DsigMF (0.958) and GaussMF (0.951). The integration of ANFIS and GIS generally led to high LSSM prediction accuracies. This study demonstrated that the choice of training dataset and the MF significantly affects the results.     

The effect of non-uniform magnetic field on the slow mode oscillations

In this paper, the slow MHD mode oscillations of the coronal plasma are studied. The aim is to identify the effect of structuring (such as magnetic field, temperature, density, and pressure) on the frequencies of oscillations. We modelled the coronal medium as a low-β plasma with longitudinally density and pressure stratifications and a weakly inhomogeneous magnetic field varied slowly with height and radial directions. The linearized ideal MHD equations reduced to a single Klein–Gordon differential equation for square of oscillatory frequencies. The eigenfunctions and analytical dispersion relations are derived. The dispersion relations were solved numerically. In the case of uniform magnetic field, the previous studies verified. Our numerical results show that, the frequencies and their ratios are very sensitive functions of pressure scale height, and slightly varying functions of inhomogeneity parameter of magnetic field. By changing the magnetic field strength between the apex and footpoints of the loop about 50%, the frequencies ratio are changed about 5%. We concluded that, the pressure scale height and temperature gradient are first order effects and inhomogeneity of magnetic field is a second order effect on the slow mode oscillations.     

Risk assessment of polycyclic aromatic hydrocarbons in the West Port semi-enclosed basin (Malaysia)

Sediment samples collected from the West Port, the west coastal waters of Malaysia, were analyzed by standard methods to determine the degree of hydrocarbon contamination and identify the sources of polyaromatic hydrocarbons (PAHs). Concentrations of PAHs in the port sediments ranged from 100.3 to 3,446.9 μg/kg dw. The highest concentrations were observed in stations close to the coastline, locations affected by intensive shipping activities and industrial input. These were dominated by high-molecular-weight PAHs (4–6 rings). Source identification showed that PAHs originated mostly pyrogenically, from the combustion of fossil fuels, grass, wood, and coal or from petroleum combustion. Regarding ecological risk estimation, only station 7 was moderately polluted, the rest of the stations suffered rare or slight adverse biological effects with PAH exposure in surface sediment, suggesting that PAHs are not considered as contaminants of concern in the West Port.     

Landslide detection using deep learning and object-based image analysis

Landslides - Recent landslide detection studies have focused on pixel-based deep learning (DL) approaches. In contrast, intuitive annotation of landslides from satellite imagery is based on...     

A shrinking particle—shrinking core model for leaching of a zinc ore containing silica

A new mathematical model was developed for leaching of zinc ores containing silicates such as hemimorphite which produce a gel during leaching with sulfuric acid. This model is based on the shrinking core model in which the particle size and the reacting core shrink simultaneously. It was shown that the actual dissolution time of the ore particles is longer than the time corresponding to the dissolution of chemical zinc oxide itself. It was suggested that because of the existence of silicates in the ore, a gelatinous layer was formed around the reacting core. Since the gel product is soft, it breaks apart when the particles collide and as a result, the particles shrink. However, a thin gelatinous layer always covers the reacting core which increases the mass transfer resistance and increases the leaching time. This model was applied to leaching of a zinc-rich tailing containing hemimorphite and the thickness of the gelatinous layer as well as the diffusion coefficient in this layer was determined.