Frictional contact algorithms in SPH for the simulation of soil–structure interaction

Simulation of frictional contact between soils and rigid or deformable structure in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two algorithms are implemented into the SPH code to describe contact behavior, where the contact forces are calculated using the law of conservation of momentum based on ideal plastic collision or using the criteria of partial penetrating. In both algorithms, the problem of boundary deficiency inherited from SPH is properly handled so that the particles located at contact boundary can have precise acceleration, which is critical for contact detection. And the movement and rotation of the rigid structure are taken into account so that it is easy to simulate the process of pile driving or movement of a retaining wall in geotechnical engineering analysis. Furthermore, the capability of modeling deformability of a structure during frictional contact simulations broadens the fields of SPH application. In contrast to previous work dealing with contact in SPH, which usually use particle‐to‐particle contact or ignoring sliding between particles and solid structure, the method proposed here is more efficient and accurate, and it is suitable to simulate interaction between soft materials and rigid or deformable structures, which are very common in geotechnical engineering. A number of numerical tests are carried out to verify the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions or results from finite element method analysis. Good agreement obtained from these comparisons suggests that the proposed algorithms are robust and can be applied to extend the capability of SPH in solving geotechnical problems. Copyright © 2013 John Wiley & Sons, Ltd.     

Practical second‐order reliability analysis applied to foundation engineering

A practical and efficient approach of implementing second‐order reliability method (SORM) is presented and illustrated for cases related to foundation engineering involving explicit and implicit limit state functions. The proposed SORM procedure is based on an approximating paraboloid fitted to the limit state surface in the neighborhood of the design point and can be easily carried out in a spreadsheet. Complex mathematical operations are relegated to relatively simple user‐created functions. The failure probability is calculated automatically based on the reliability index and principal curvatures of the limit state surface using established closed‐form SORM formulas. Four common foundation engineering examples are analyzed using the proposed method and discussed: immediate settlement of a flexible rectangular foundation, bearing capacity of a shallow footing, axial capacity of a vertical single pile, and deflection of a pile under lateral load. Comparisons with Monte Carlo simulations are made. In the case of the laterally loaded pile, the friction angle of the soil is represented as a one‐dimensional random field, and pile deflections are computed based on finite element analysis on a stand‐alone computer package. The implicit limit state function is approximated via the response surface method using two quadratic models. Copyright © 2011 John Wiley & Sons, Ltd.     

Evaluation of multiple regression models using spatial variables to predict nitrate concentrations in volcanic aquifers

Multiple linear regression of spatial variables including land use, soil type, and topography was applied to predict nitrate concentration and evaluate major factors affecting nitrate occurrence in springs and wells in the southern and northern areas of Jeju volcanic island, Korea. Three types of contributing area surrogates (CAS), namely circle, semicircle, and wedge, were employed to calculate the spatial variables. The regression results showed R² of 0.81–0.84 for springs and 0.74–0.77 for wells; R² values for wedge and semicircular CAS were more than 10% higher than those for circular CAS. The R² of spring models was significantly affected by both the shape and size of CAS, with optimal radii of 150–250 m and 300–400 m in the southern and northern areas, respectively, corresponding to thinner upper basaltic aquifers, and implying shorter flow paths in the southern area. The most influential variables in springs were orchards and soil types related to agriculture including silty loam and silty clay loam, indicating that nitrate levels are strongly affected by N fertilization in cultivated areas. In contrast, wells showed much less sensitivity to both shapes and sizes of CAS, with less contribution of land use and soil type to the regression, which could be attributed to a mix of multiple aquifer zones and widely different factors in the installation and operation of wells. Field parameters of electrical conductivity (EC) and pH increased the R² up to 10%, suggesting that these can be useful when regression with spatial variables yields a lower R². The optimal spatial scales for prediction of nitrate concentration and spatial variables that significantly contribute to nitrate contamination can provide relevant criteria for establishing groundwater management policies, considering the increasing anthropogenic land‐use trends on the island, where groundwater is highly sensitive to changes in spatial variables. Copyright © 2015 John Wiley & Sons, Ltd.     

Seasonal variations of imposex indices and butyltin concentrations in the rock shell Thais clavigera collected from Hong Kong waters

This study investigated butyltin contamination in Hong Kong waters using the rock shell Thais clavigera as a biomonitor. We collected rock shells from nine coastal sites in both winter 2005 and summer 2006 to examine the spatial and seasonal patterns. The relative penis size index (RPSI) and vas deferens sequence index (VDSI) were determined following a standard protocol, while the tissue concentrations of butyltins (i.e., mono-, di- and tri-butyltin) were measured using gas chromatography. The results showed that imposex was particularly severe in T. clavigera collected from Victoria Harbour, Pak Sha Wan and Sai Kung Pier where shipping activities were frequent. Both imposex indices had a positive correlation with tissue concentration of tributyltin in T. clavigera. While VDSI exhibited little seasonal variation, RPSI of winter specimens was significantly higher than that of summer specimens. Thus, sampling season should be standardised when using RPSI as a biomarker.     

Molecular staging of marine medaka: a model organism for marine ecotoxicity study

Oryzias melastigma, also called O. dancena, is becoming a very useful model for estuarine and marine ecotoxicity studies. With O. melastigma being adopted by ILSI Health and Environmental Science Institute (HESI) for embryo toxicity testing, improved knowledge of biomarker based embryonic development becomes especially important for mechanism-based toxicity evaluations. Using whole mount in situ hybridization and immunostaining techniques together with widely used molecular markers, this study describes the molecular development of marine medaka embryos, focusing on the brain, eye, heart, pectoral fin, pancreas, liver, muscle and neuron system. These organs are targets of environmental pollutants that disrupt normal embryonic development in medaka and other fish.     

Numerical simulation of vertical ground heat exchangers: Intermittent operation in unsaturated soil conditions

The effect of varying the thermal properties of inhomogeneous unsaturated soil on the intermittent operation of a vertical ground heat exchanger (GHE) was simulated by a conjugate heat transfer simulation using a transient conductive heat transfer model. A three-phase soil model was used to introduce soil properties that vary with depth. The performance during the first few hours was significantly different from that of an analytical infinite line source model that assumes steady-state borehole conditions, although relatively good agreement was obtained thereafter. Unsaturated soil conditions afforded a 40% lower mean heat exchange rate than saturated conditions. This demonstrates the importance of considering unsaturated conditions in the design and performance evaluation of GHEs.     

The distribution of annual maximum earthquake magnitude around Taiwan and its application in the estimation of catastrophic earthquake recurrence probability

The annual maximum earthquake magnitudes around Taiwan from 1900 to 2009 are presented in this paper. Using the distribution of the AMEM, a probabilistic framework to estimate the recurrence probability of a large-size earthquake is also proposed and an illustration was made in this paper. The mean value of the 110-AMEM is 6.433, and the coefficient of variation is around 10%. The results of two goodness-of-fit tests show that the Gamma and lognormal distributions are relatively suitable to represent the AMEM around Taiwan among five common probability distributions. Using the proposed approach, the recurrence probability is 4% for an earthquake with magnitude greater than 7.5 in a 1-year period around Taiwan. More site-specifically, the probability is around 5% in Central Taiwan for such an earthquake to occur in a 50-year period.     

Monitoring Residual Spatial Patterns using Bayesian Hierarchical Spatial Modelling for Exploring Unknown Risk Factors

This article studies Bayesian hierarchical spatial modelling that monitors the changes of residual spatial pattern (structure) of the outcome variable for exploring unknown risk factors in small‐area analysis. Spatially structured random effects (SRE) and unstructured random effects (URE) terms added to the conventional logistic regression model take into account overdispersion and residual spatial structure, which if unaccounted for could cause incorrect identification of risk factors. Mapping and/or calculating the ratio of random effects that are spatially‐structured monitor the extent of residual spatial structure. The monitoring provides insights into identification of unknown covariates that have similar spatial structures to those of SRE. Adding such covariates to the model has the potential to diminish the residual spatial structure, until possibly all or most of the spatial structure can be explained. Risk factors identified are the added covariates that have statistically significant regression coefficients. We apply the methods to the analysis of domestic burglaries in Cambridgeshire, England. Small‐area analysis of crime where data often display apparent spatial structure would particularly benefit from the methodologies. We discuss the methodologies, their relevancy in our analysis of domestic burglaries, their limitations, and possible paths for future research.     

Thermohaline structures across the Luzon Strait from seismic reflection data

A legacy seismic section across the Luzon Strait was reprocessed using the pre-stack depth migration method to reveal thermohaline structures in the water column. Distinct finescale and mesoscale features can be seen from the seismic image. Vertically, reflective patterns are associated with three water layers: upper layer, intermediate layer, and deep layer. The upper layer can be divided into three areas by two opposing weak-to-transparent zones from west to east with individual thicknesses: ～400, 550, and 700 m. They are interpreted as the South China Sea upper water, the Kuroshio water, and the Pacific upper water, respectively, separated by two Kuroshio frontal zones. Internal waves are ubiquitous in the Kuroshio water fluctuating in different amplitudes and wavelengths laterally. The finestructure shows the western Kuroshio front zone is composed of three subparallel blanking zones with an average eastward dipping angle of 2.3°. They are regarded as the well-mixed frontal interfaces developed in and along a multi-frontal frontal zone. The transparent reflection of the intermediate water above the Hengchun is suggestive of a well-mixed layer by the near-bottom turbulence. In the Pacific intermediate water near the eastern mouth of the Bashi Channel, there is an inverted crescent-shaped structure, which is likely associated with the intrusion of the South China Sea intermediate water through the Bashi Channel. These imaged features are consistent with the observed oceanographic phenomena in the region.