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.     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

Jurassic igneous rocks of the central Sanandaj–Sirjan zone (Iran) mark a propagating continental rift,not a magmatic arc

Jurassic igneous bodies of the Sanandaj–Sirjan zone (SaSZ) in SW Iran are generally considered as a magmatic arc but critical evaluation of modern geochronology, geochemistry and radiogenic isotopes challenges this conclusion. There is no evidence for sustained igneous activity along the ~1,200 km long SaSZ, as expected for a convergent plate margin; instead activity was brief at most sites and propagated NW at ~20 mm/a. Jurassic igneous rocks define a bimodal suite of gabbro‐diorite and granite. Chemical and isotopic compositions of mafic rocks indicate subcontinental lithospheric mantle sources that mostly lacked subduction‐related modifications. The arc‐like features of S‐type granites reflect massive involvement of Cadomian crust and younger sediments to generate felsic melts in response to mafic intrusions. We conclude that Jurassic SaSZ igneous activity occurred in a continental rift, not an arc. SaSZ igneous rocks do not indicate that subduction along the SW margin of Eurasia began in Jurassic time.     

Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.     

Active capture and anomalous adsorption: New mechanisms for the incorporation of heavy noble gases

Abstract— Active capture is a new process for the incorporation of large quantities of heavy noble gases into growing surfaces. Adsorption in the conventional sense involves surface bonding by polarization (Van der Waals forces). What is referred to as “anomalous adsorption” of heavy noble gases involves chemical bonds and can occur when other (more chemically active) species are not available to preempt sites with unfilled bonds. Anomalous adsorption has been observed under conditions of fracture, vacuum deposition and ionizing radiation. Active capture depends upon anomalous adsorption to retain noble gases on a surface long enough to be captured in a growing surface film as it is deposited. The fundamental principle may be the impingement onto the growing film with sufficient energy to liberate surface electrons (work function energy of a few electronvolts) so that they are retained by anomalous adsorption long enough to be entrapped in the growing surface. Trapping efficiencies of ?1% have been observed for Kr and Xe in laboratory experiments, implying a fundamentally new mechanism for the incorporation of heavy noble gases onto surfaces. It may play a role in explaining the large concentrations of planetary noble gases contained in phase‐Q.     

A spectroscopic study of the peculiar galaxy UGC 5600

We present our observations of the galaxy UGS 5600 with a long-slit spectrograph (UAGS) and a multipupil field spectrograph (MPFS) attached to the 6-m Special Astrophysical Observatory telescope. Radial-velocity fields of the stellar and gaseous components were constructed for the central region and inner ring of the galaxy. We proved the existence of two nearly orthogonal kinematic subsystems and conclude that UGC 5600 is a galaxy with an inner polar ring. In the circumnuclear region, we detected noncircular stellar motions and suspected the existence of a minibar. The emission lines are shown to originate in H II regions. We estimated the metallicity from the intensity ratio of the [N II]λ6583 and Hα lines to be nearly solar, which rules out the possibility that the polar ring was produced by the accretion of gas from a dwarf companion.