Three‐dimensional finite element study of arching behavior in slope/drilled shafts system

Two‐dimensional slope stability analysis for a slope with a row of drilled shafts needs a mechanism to take into account the three‐dimensional effect of the soil arching due to the spaced drilled shafts on slope. To gain a better understanding of the arching mechanisms in a slope with evenly spaced drilled shafts socketed into a stable stratum (or a rock layer), the three‐dimensional finite element modelling technique was used for a comprehensive parametric study, where the nonlinear and plastic nature of the soil and the elastic behavior of the drilled shafts as well as the interface frictions were modelled. Various factors were varied in the parametric study to include (1) the rigidity of the drilled shafts as influenced by its diameter, modulus of elasticity, and total length; (2) shafts spacing and location on the slope; (3) the material properties of rock and the socket length of shaft; and (4) the soil movement and strength parameters. Evidences of soil arching and reduction in the stresses and displacements through the load transfer mechanisms due to the presence of the drilled shafts were elucidated through the finite element method (FEM) simulation results. Design charts based on regression analysis of FEM simulation results were created to obtain a numerical value of the load transfer factor for the arching mechanism provided by the drilled shafts on the slope. Observations of the arching behavior learned from the FEM simulations provide an insight into the behavior of drilled shafts stabilized slope. Copyright © 2009 John Wiley & Sons, Ltd.     

A Unified Computational Method for Simulating Dynamic Behavior of Planing Vessels

High speed planing hulls have complex hydrodynamic behaviors. The trim angle and drafts are very sensitive to speed and location of the center of gravity. Therefore, motion simulation for such vessels needs a strong coupling between rigid body motions and hydrodynamic analysis. In addition, free surface should be predicted with good accuracy for each time step. In this paper, velocity and pressure fields are coupled by use of the fractional step method. On the basis of integration of the two-phase viscous f...     

Towards simulation of 3D nonlinear high-speed vessels motion

A numerical simulation algorithm based on the finite volume discretisation is presented for analyzing ship motions. The algorithm employs a fractional step method to deal with the coupling between the pressure and velocity fields. The free surface capturing is fulfilled by using a volume of fluid method in which the interface between the liquid (water) and gas (air) phases are computed by solving a scalar transport equation for the volume fraction of the liquid phase. The computed velocity field is employed to evaluate the acting forces and moments on the vessel. Using the strategy of boundary-fitted body-attached mesh and calculating all six degrees-of-freedom of motion in each time step, time history of ship motions including displacements, velocities and accelerations are evaluated.To verify the proposed algorithm, a series of verification tests are conducted. First, a two-dimensional asymmetrical wedge slamming is simulated as a simple type of a common case for high-speed vessels. Then, the steady-state forward motion of a high-speed planing catamaran is investigated. Results of both test cases show good agreement with experimental data. It is concluded that the proposed algorithm can be a promising strategy for both performance prediction and design of high-speed vessels.     

Modeling the near-field effects of the worst-case tsunami in the Makran subduction zone

As a first step towards the development of inundation maps for the northwestern Indian Ocean, we simulated the near-field inundation of two large tsunami in the Makran subduction zone (MSZ). The tsunami scenarios were based on large historical earthquakes in the region. The first scenario included the rupture of about 500 km of the plate boundary in the eastern MSZ, featuring a moment magnitude of M_w 8.6. The second scenario involved the full rupture of the plate boundary resulting from a M_w 9 earthquake. For each scenario, the distribution of tsunami wave height along the coastlines of the region is presented. Also, detailed runup modeling was performed at four main coastal cities in the region for the second scenario. To investigate the possible effect of splay fault branching on tsunami wave height, a hypothetical splay fault was modeled which showed that it can locally increase the maximum wave height by a factor of 2. Our results showed that the two tsunami scenarios produce a runup height of 12-18 m and 24-30 m, respectively. For the second scenario, the modeled inundation distance was between 1 and 5 km.     

Population ecology of fiddler crab, Uca lactea annulipes (Decapoda: Ocypodidae) in Sirik mangrove estuary,Iran

Fiddler crabs are known as the most conspicuous bioturbating animals in mangrove forest, which by their burrowing activities promote nutrient cycling in ecosystems. Mangroves, crabs and bacteria are connected in the form of positive feedback loop, whereby an increase in the activity of any participant will tend to increase the activity of all others as well. Population structure, growth, mortality, recruitment, biomass and production of fiddler crabs as dominant species of deposit-feeder crabs have been investigated from May 2005 to March 2006 in the Sirik mangrove estuary in the eastern part of Strait of Hormuz. Ten 1 m² squares of the substrate were randomly sampled during low tide. There was a significant change in organic content in the sampling area during study period. The sex ratio was significantly deviated from 1:1 ratio in large sizes, where males were dominated year around. Fiddler crabs showed seasonal growth pattern during a year, which the maximum growth was observed during autumn and early summer and ceased during winter by low temperature. Growth parameters were estimated by Appeldoorn's method and fitted with seasonalized von Bertalanffy growth models. The Kimura likelihood ratio test revealed significant differences in the overall growth curves between males and females, originated from the differences in k$k$ and L_∞$L_{\mathrm{\infty }}$ between sexes. Males were found to exhibit higher growth rate and larger size than females, this could be attributed to the higher investment of energy for reproduction in females. Reproduction was restricted to spring and early summer, while recruitment occurs at the end of summer. Total annual production of males and females was 37.90 and 10.05 g dry weight/10 m² respectively and major production took place during early summer. Because of impact of temperature on the activity of fiddler crabs in subtropical regions, density of crabs has significant positive correlation with temperature (P < 0.01), while growth rate shows significant negative correlation with organic matter of sediment (P < 0.05). In conclusion these crabs showed lower mortality and production rate with the seasonal growth and recruitment in comparison with tropical species.     

Assessing the effect of temporal dynamics on urban growth simulation: Towards an asynchronous cellular automata

Time is a fundamental dimension in urban dynamics, but the effect of various definitions of time on urban growth models has rarely been evaluated. In urban growth models such as cellular automata (CA), time has typically been defined as a sequence of discrete time steps. However, most urban growth processes such as land‐use changes are asynchronous. The aim of this study is to examine the effect of various temporal dynamics scenarios on urban growth simulation, in terms of urban land‐use planning, and to introduce an asynchronous parcel‐based cellular automata (AParCA) model. In this study, eight different scenarios were generated to investigate the impact of temporal dynamics on CA‐based urban growth models, and their outputs were evaluated using various urban planning indicators. The obtained results show that different degrees of temporal dynamics lead to various patterns appearing in urban growth CA models, and the application of asynchronous (event‐driven) CA models achieves better simulation results than synchronous models.     

Bearing capacity performance of soft cohesive soil treated by kenaf limited life geotextile

Abstract

This article presents the performance of the short-term bearing capacity on soft clay soil treated by Kenaf geotextile under vertical loading via a small-scale modelling test at unit gravity. The ground model was formulated by consolidating kaolin in a rigid testing compartment. In the loading test, the strip footing was represented by a rigid footing. For the treated case, a series of tests were performed to examine the effects of the burial depth of the Kenaf geotextile on the bearing capacity of the soft soil. The Kenaf geotextile was laid beneath the rigid footing (at the ground surface) and at 50, 75 and 100?mm depth from the soil surface. All the measured results of the Kenaf geotextile treated ground were compared with the untreated ground. The incorporation of Kenaf fibre geotextile was observed to enhance the bearing capacity of soft cohesive clay up to 281% depending upon the depth of the installed geotextile. The geotextile at the surface provided the highest bearing capacity and sustained the highest displacement at failure. The outcome of this research will promote the use of natural fibre geotextiles as sustainable earth reinforcement in temporary earthwork applications.     

Combined Gravimetric–Seismic Crustal Model for Antarctica

The latest seismic data and improved information about the subglacial bedrock relief are used in this study to estimate the sediment and crustal thickness under the Antarctic continent. Since large parts of Antarctica are not yet covered by seismic surveys, the gravity and crustal structure models are used to interpolate the Moho information where seismic data are missing. The gravity information is also extended offshore to detect the Moho under continental margins and neighboring oceanic crust. The processing strategy involves the solution to the Vening Meinesz-Moritz’s inverse problem of isostasy constrained on seismic data. A comparison of our new results with existing studies indicates a substantial improvement in the sediment and crustal models. The seismic data analysis shows significant sediment accumulations in Antarctica, with broad sedimentary basins. According to our result, the maximum sediment thickness in Antarctica is about 15 km under Filchner-Ronne Ice Shelf. The Moho relief closely resembles major geological and tectonic features. A rather thick continental crust of East Antarctic Craton is separated from a complex geological/tectonic structure of West Antarctica by the Transantarctic Mountains. The average Moho depth of 34.1 km under the Antarctic continent slightly differs from previous estimates. A maximum Moho deepening of 58.2 km under the Gamburtsev Subglacial Mountains in East Antarctica confirmed the presence of deep and compact orogenic roots. Another large Moho depth in East Antarctica is detected under Dronning Maud Land with two orogenic roots under Wohlthat Massif (48–50 km) and the Kottas Mountains (48–50 km) that are separated by a relatively thin crust along Jutulstraumen Rift. The Moho depth under central parts of the Transantarctic Mountains reaches 46 km. The maximum Moho deepening (34–38 km) in West Antarctica is under the Antarctic Peninsula. The Moho depth minima in East Antarctica are found under the Lambert Trench (24–28 km), while in West Antarctica the Moho depth minima are along the West Antarctic Rift System under the Bentley depression (20–22 km) and Ross Sea Ice Shelf (16–24 km). The gravimetric result confirmed a maximum extension of the Antarctic continental margins under the Ross Sea Embayment and the Weddell Sea Embayment with an extremely thin continental crust (10–20 km).     

Application of geostatistical simulation to compile seismotectonic provinces based on earthquake databases (case study: Iran)

This article is devoted to application of a simulation algorithm based on geostatistical methods to compile and update seismotectonic provinces in which Iran has been chosen as a case study. Traditionally, tectonic maps together with seismological data and information (e.g., earthquake catalogues, earthquake mechanism, and microseismic data) have been used to update seismotectonic provinces. In many cases, incomplete earthquake catalogues are one of the important challenges in this procedure. To overcome this problem, a geostatistical simulation algorithm, turning band simulation, TBSIM, was applied to make a synthetic data to improve incomplete earthquake catalogues. Then, the synthetic data was added to the traditional information to study the seismicity homogeneity and classify the areas according to tectonic and seismic properties to update seismotectonic provinces. In this paper, (i) different magnitude types in the studied catalogues have been homogenized to moment magnitude (M_w), and earthquake declustering was then carried out to remove aftershocks and foreshocks; (ii) time normalization method was introduced to decrease the uncertainty in a temporal domain prior to start the simulation procedure; (iii) variography has been carried out in each subregion to study spatial regressions (e.g., west-southwestern area showed a spatial regression from 0.4 to 1.4 decimal degrees; the maximum range identified in the azimuth of 135?±?10); (iv) TBSIM algorithm was then applied to make simulated events which gave rise to make 68,800 synthetic events according to the spatial regression found in several directions; (v) simulated events (i.e., magnitudes) were classified based on their intensity in ArcGIS packages and homogenous seismic zones have been determined. Finally, according to the synthetic data, tectonic features, and actual earthquake catalogues, 17 seismotectonic provinces were introduced in four major classes introduced as very high, high, moderate, and low seismic potential provinces. Seismotectonic properties of very high seismic potential provinces have been also presented.     

Multicomponent competitive monovalent cation exchange in hierarchical porous media with multimodal reactive mineral facies

This paper presents a stochastic model for multicomponent competitive monovalent cation exchange in hierarchical porous media. Reactive transport in porous media is highly sensitive to heterogeneities in physical and chemical properties, such as hydraulic conductivity (K), and cation exchange capacity (CEC). We use a conceptual model for multimodal reactive mineral facies and develop a Eulerian-based stochastic theory to analyze the transport of multiple cations in heterogeneous media with a hierarchical organization of reactive minerals. Numerical examples investigate the retardation factors and dispersivities in a chemical system made of three monovalent cations (Na⁺, K⁺, and Cs⁺). The results demonstrate how heterogeneity influences the transport of competitive monovalent cations, and highlight the importance of correlations between K and CEC. Further sensitivity analyses are presented investigating how the dispersion and retardation of each cation are affected by the means, variances, and integral scales of K and CEC. The volume fraction of organic matter is shown to be another important parameter. The Eulerian stochastic framework presented in this work clarifies the importance of each system parameters on the migration of cation plumes in formations with hierarchical organization of facies types. Our stochastic approach could be used as an alternative to numerical simulations for 3D reactive transport in hierarchical porous media, which become prohibitively expensive for the multicomponent applications considered in this work.