Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.     

An analytical solution for deforming twin‐parallel tunnels in an elastic half plane

The interaction between twin‐parallel tunnels affects the tunnelling‐induced ground deformation, which may endanger the nearby structures. In this paper, an analytical solution is presented for problems in determining displacements and stresses around deforming twin‐parallel tunnels in an elastic half plane, on the basis of complex variable theory. As an example, a uniform radial displacement was assumed as the boundary condition for each of the two tunnels. Special attention was paid to the effects of tunnel depth and spacing between the two tunnels on the surface movement to gain deep insight into the effect of the interaction between twin‐parallel tunnels using the proposed analytical approach. It is revealed that the influence of twin tunnel interaction on surface movements diminishes with both the increase of the tunnel depth and the spacing between the two tunnels. The presented analytical solution manifests that, similar to most of the existing numerical results, the principle of superposition can be applied to determine ground deformation of twin‐parallel tunnels with a certain large depth and spacing; otherwise, the interaction effect between the two tunnels should be taken into account for predicting reliable ground movement. Copyright © 2014 John Wiley & Sons, Ltd.     

Empirical assessment of rockfall and debris flow risk along the Karakoram Highway,Pakistan

Natural Hazards - The Karakoram Highway links north Pakistan with southwest China. It passes through unique geomorphological, geological and tectonic setting. This study focused 200-km-long section...     

Investigation of Failed Slopes Based on a New Method of Back Analysis in Rock Masses: A Case Study

Geotechnical and Geological Engineering - Small variation in shear strength parameters results in remarkable changes in the safety factor (SF) of a rock slope. In this regard, rock mass strength of...     

Geochemical Prospectivity Mapping Through a Feature Extraction–Selection Classification Scheme

Natural Resources Research - Machine learning (ML) schemes can enhance success in geochemical prospectivity mapping. This study has examined the effectiveness of several feature extraction or...     

An integrated object-based image analysis and CA-Markov model approach for modeling land use/land cover trends in the Sarab plain

The present paper is an attempt to integrate a semi-automated object-based image analysis (OBIA) classification framework and a cellular automata-Markov model to study land use/land cover (LULC) changes. Land use maps for the Sarab plain in Iran for the years 2000, 2006, and 2014 were created from Landsat satellite data, by applying an OBIA classification using the normalized difference vegetation index, salinity index, moisture stress index, soil-adjusted vegetation index, and elevation and slope indicators. The classifications yielded overall accuracies of 91, 93, and 94% for 2000, 2006, and 2014, respectively. Finally, using the transition matrix, the spatial distribution of land use was simulated for 2020. The results of the study revealed that the number of orchards with irrigated agriculture and dry-farm agriculture in the Sarab plain is increasing, while the amount of bare land is decreasing. The results of this research are of great importance for regional authorities and decision makers in strategic land use planning.     

The removal of color and COD from wastewater containing water base color by coagulation process

In this study acrylic water base color was removed from synthetic wastewater using coagulation process. Experiments were conducted on the sample containing 100 mg/L and 400 mg/L of acrylic water base color. Destruction of color by means of coagulation/flocculation techniques using ferrous sulfate, alum, lime and polyelectrolyte (cationic, anionic and non ionic). The study was performed in a systematic approach searching optimum values of alum and FeSO₄ concentration, pH and temperature. All the experiments were run in a laboratory scale.The obtained results show that treatment with alum and ferrous sulfate alone proved to be very effective in removing the color (> 99 %) and part of COD (60–70 %) from aqueous solution. Lime alone did not significant change on COD and color removal.     

Higher-order compositional modeling with Fickian diffusion in unstructured and anisotropic media

We present advances in compositional modeling of two-phase multi-component flow through highly complex porous media. Higher-order methods are used to approximate both mass transport and the velocity and pressure fields. We employ the Mixed Hybrid Finite Element (MHFE) method to simultaneously solve, to the same order, the pressure equation and Darcy's law for the velocity. The species balance equation is approximated by the discontinuous Galerkin (DG) approach, combined with a slope limiter. In this work we present an improved DG scheme where phase splitting is analyzed at all element vertices in the two-phase regions, rather than only as element averages. This approximation is higher-order than the commonly employed finite volume method and earlier DG approximations. The method reduces numerical dispersion, allowing for an accurate capture of shock fronts and lower dependence on mesh quality and orientation. Further new features are the extension to unstructured grids and support for arbitrary permeability tensors (allowing for both scalar heterogeneity, and shear anisotropy). The most important advancement in this work is the self-consistent modeling of two-phase multi-component Fickian diffusion. We present several numerical examples to illustrate the powerful features of our combined MHFE–dg method with respect to lower-order calculations, ranging from simple two component fluids to more challenging real problems regarding CO₂ injection into a vertical domain saturated with a multi-component petroleum fluid.     

Numerical modeling of two-phase flow in heterogeneous permeable media with different capillarity pressures

Contrast in capillary pressure of heterogeneous permeable media can have a significant effect on the flow path in two-phase immiscible flow. Very little work has appeared on the subject of capillary heterogeneity despite the fact that in certain cases it may be as important as permeability heterogeneity. The discontinuity in saturation as a result of capillary continuity, and in some cases capillary discontinuity may arise from contrast in capillary pressure functions in heterogeneous permeable media leading to complications in numerical modeling. There are also other challenges for accurate numerical modeling due to distorted unstructured grids because of the grid orientation and numerical dispersion effects. Limited attempts have been made in the literature to assess the accuracy of fluid flow modeling in heterogeneous permeable media with capillarity heterogeneity. The basic mixed finite element (MFE) framework is a superior method for accurate flux calculation in heterogeneous media in comparison to the conventional finite difference and finite volume approaches. However, a deficiency in the MFE from the direct use of fractional flow formulation has been recognized lately in application to flow in permeable media with capillary heterogeneity. In this work, we propose a new consistent formulation in 3D in which the total velocity is expressed in terms of the wetting-phase potential gradient and the capillary potential gradient. In our formulation, the coefficient of the wetting potential gradient is in terms of the total mobility which is smoother than the wetting mobility. We combine the MFE and discontinuous Galerkin (DG) methods to solve the pressure equation and the saturation equation, respectively. Our numerical model is verified with 1D analytical solutions in homogeneous and heterogeneous media. We also present 2D examples to demonstrate the significance of capillary heterogeneity in flow, and a 3D example to demonstrate the negligible effect of distorted meshes on the numerical solution in our proposed algorithm.