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.     

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.     

任意曲线曲面上重磁位场转换的B样条函数法

提出用B样条函数求解曲线、曲面上重磁位场的向上延拓,水平、垂向导数计算,磁异常分量互换的方法。该方法的特点是:原理简明,程序通用性强,计算精度高。     

In situ and satellite observations of the eastward migration of the Western Alboran Sea Gyre

During October 2003 an intensive oceanographic survey (BIOMEGA) was carried out in the Alboran Sea, coinciding with a migration event of the Western Alboran Sea Gyre (WAG). The observations gathered during that cruise constitute the first field evidence of a migrated stage of the WAG. In this work we present the main differences between the 3D hydrodynamic fields observed during BIOMEGA and those corresponding to a WAG located at its usual position. The migration of the gyre was followed by satellite (altimetry and sea surface temperature) imagery. The causes of the gyre migration are explored in terms of the quasi-geostrophic tendency equation, in particular of the dynamics governing scales larger than the Rossby radius of deformation. It is shown that the steady state gyre must be almost equivalent barotropic and that the key process to break down the stationarity would be a density advection at gyre scale. The mechanisms to explain the migration of the WAG proposed by previous authors are discussed in light of the explanation proposed in this work.