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.     

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.     

基于神经网络的平底结构砰击压力预报

对利用神经网络预报平底结构入水砰击压力的方法进行了探讨。首先利用仿真软件计算各种情况下平底结构入水所产生的砰击压力，以此形成训练神经网络的数据集。其次利用数据集对三层反馈式网络进行了训练，讨论了不同隐含层节点数对该非线性系统的拟合能力，并且对梯度下降法、动量修正法和基于优化的LM算法的有效性和精度进行了比较，最后得出了适合平底结构入水砰击系统的网络结构。     

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.