Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Infiltration at the Tai rain forest (Ivory Coast): Measurements and modelling

Infiltration experiments have been performed at three sites along a well-known catena under virgin tropical rain forest using a portable sprinkling infiltrometer. Experimentally determined infiltration curves are presented. Infiltration curves are also simulated on the basis of the Mein-Larson equation. The parameters for this model have been obtained from the infiltration curves (saturated conductivity) and simple soil moisture determinations (fillable porosity). The agreement between experimentally determined and modelled infiltration is reasonable, provided (a) saturated conductivity as derived from the experimental data is corrected, (b) a storage parameter, also derived from the experimental data, is added to the Mein-Larson model, and (c) the decline in soil porosity with depth is either small or occurs abruptly at shallow depth. Comparison of observed infiltration rates with rainfall intensity shows that Horton Overland Flow has to occur naturally at least on the middle and lower section of the catena. Despite the fact that most parameters can be estimated in principle from basic soil data, it remains advisable to obtain sprinkling infiltrometer field measurements, because of soil variability due to dynamic surface conditions, macroporosity, air entrapment, and irregularity of the wetting front.     

Doline probability map using logistic regression and GIS technology in the central Ebro Basin (Spain)

In the surroundings of Zaragoza, karstification processes are especially intense in covered karst areas where fluvial terraces lie directly on Tertiary evaporites. Since the beginning of Quaternary, these processes have lead to the development of collapse and subsidence dolines with a wide range of sizes, which have significant economic impacts. To reduce economic impact and increase safety, a regional analysis of this phenomenon is needed for spatial management. Therefore, a probability map of dolines was developed using logistic regression and geographic information system (GIS) techniques. This paper covers the selection of input data, manipulation of data using the GIS technology, and the use of logistic regression to generate a doline probability map. The primary variable in the doline development in this area is geomorphology, represented by the location of endorheic areas and different terrace levels. Secondary variables are the presence of irrigation and the water table gradient.     

Modelling of periodic wave transformation in the inner surf zone

In this paper, we analyse the ability of the nonlinear shallow-water (NSW) equations to predict wave distortion and energy dissipation of periodic broken waves in the inner surf zone. This analysis is based on the weak-solution theory for conservative equations. We derive a new one-way model, which applies to the transformation of non-reflective periodic broken waves on gently sloping beaches. This model can be useful to develop breaking-wave parameterizations (in particular broken-wave celerity expression) in both time-averaged wave models and time-dependent Boussinesq-type models. We also derive a new wave set-up equation which provides a simple and explicit relation between wave set-up and energy dissipation. Finally, we compare numerical simulations of both, the NSW model and the simplified one-way model, with spilling wave breaking experiments and we find a good agreement.     

关于具二次约束最小二乘问题的敏度分析

本文系统地研究具二次约束最小二乘问题的敏度分析。首先给出长期方程唯一正根的上界和下界；然后证明割线法用于计算拉格朗日乘子时全局收敛；最后给出解的扰动界。     

青岛港风暴潮经验统计预报

本文利用青岛港多年实测资料,分析了该港风暴潮概况。而后通过多元回归技术,求取了该港极值增减水的预报公式。经非独立和独立检验,结果令人满意。     

Numerical modelling of deformation of multi-directional random waves over a varying topography

Two numerical formulations of the breaking phenomenon were implemented in a numerical model for random wave propagation based on the elliptic formulation of the mild-slope equation. The randomness of the wave field was simulated based on a spectral component method, in which the 3-D spectrum is discretised in components of equal energy. One of the breaking process formulations is based on the concept of breaking each independent spectral component. The other is based on the distribution of the local amount of energy dissipated through the independent spectral components. The model based on the concept of breaking each independent spectral component produces the best estimates of the wave field, when the numerical results are compared with laboratory data.     

Comparative tolerance of two estuarine annelids to fluoranthene under normoxic and moderately hypoxic conditions

The tolerance of the oligochaete Monopylephorus rubroniveus and the polychaete Streblospio benedicti to sediment-associated fluoranthene was characterized under normoxic (>80% dissolved oxygen saturation) and moderately hypoxic ( approximately 50% dissolved oxygen saturation) conditions. Under both conditions, M. rubroniveus was highly tolerant of fluoranthene. Streblospio benedicti was considerably less tolerant of fluoranthene compared with M. rubroniveus. In addition, S. benedicti was less tolerant to fluoranthene under moderately hypoxic conditions, although no differences in sensitivity between the two oxygen conditions were observed based upon median lethal tissue residues. Bioaccumulation factors were higher for S. benedicti exposed to moderate hypoxia, suggesting that behavioral adaptations to compensate for the lower dissolved oxygen increased its bioaccumulation of fluoranthene. The results of the present laboratory study demonstrate that (1). changes in annelid tolerance to fluoranthene under varying oxygen conditions is a species-dependent phenomenon and (2). the differential tolerance of these two annelids to these combined stressors is consistent with their relative abundances in the field.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.