Mixing process in estuaries

A theoretical expression for vertical profile-of horizontal velocity in terms of its depth-average is derived based on oscillatory boundary layer theory and estuarine flow characteristics. The derived theoretical profile is then incorporated into a vertical quasi-two-dimensional model, which is proved advantageous in more physical implications and less CPU time demand. To validate the proposed model, the calculated results are compared to the field data in the Yangtze River Estuary, exhibiting good agreement with observations. The proposed quasi-two-dimensional vertical model is used to study mixing process, especially dependence of salinity distribution and salt front strength on runoff and tides in estuaries.

The cellular structure of a two-dimensional H2/O2/Ar detonation wave

In this paper, the cellular structure of a two-dimensional detonation wave in a low pressure H2/O2/Ar mixture calculated with a detailed chemical reaction model, high order scheme and high resolution grids is investigated. The regular cellular structure is produced about 1 ms after introducing perturbations in the reaction zone of a steady one-dimensional detonation wave. It is found from the present resolution study that the discrepancies concerning the structure type arising from the coarser grid employed can be resolved using a sufficiently fine grid size of 0.05 mm and below and shows a double-Mach-like strong-type configuration. During the structure evolution process, the structure configuration does not change much in the periods before and after the triple point collision. Through the triple point collision, three regular collision processes are observed and are followed by a quick change to the double-Mach-like configuration. The simulated structure tracks show that there are three different tracks associated with different triple points or the kink on the transverse wave. Comparisons with previous work and experiments indicate the presence of a strong structure for an ordinary detonation.

Direct numerical test of the B-G-K model equation by the DSMC method

In the present paper the rarefied gas how caused by the sudden change of the wall temperature and the Rayleigh problem are simulated by the DSMC method which has been validated by experiments both in global flour field and velocity distribution function level. The comparison of the simulated results with the accurate numerical solutions of the B-G-K model equation shows that near equilibrium the BG-K equation with corrected collision frequency can give accurate result but as farther away from equilibrium the B-G-K equation is not accurate. This is for the first time that the error caused by the B-G-K model equation has been revealed.

球面聚心气相爆轰波传播过程的数值研究

应用频散可控耗散格式对球面聚心气相爆轰波的传播过程进行了数值模拟研究．通过跟踪波阵面上压力和温度变化，分析这些参数在爆轰波传播过程中的演变规律，及其与几何尺度和初始条件之间的依赖关系．研究结果表明，仅在远场波面压力的变化近似只依赖于r／R，对称中心附近则需要考虑初始半径R的影响；波面压力与初始压力的变化呈线性关系；汇聚过程中温度升高比压力慢得多．

湖北省清江隔河岩水库茅坪滑坡及白岩危岩体简介

茅坪滑坡 茅坪滑坡位于湖北清江隔河岩水库库区近尾部左岸，下距隔河岩水库大坝66km，上距正在施工的水布垭坝址25km。它是目前库区规模最大、变形最严重的一处堆积层滑坡。因该滑坡有发生堵江的极大可能性而倍受有关方面关注。目前，该滑坡已被选作中科院知识创新工程重要方向支持项目“山体滑坡灾害防治中的关键力学问题研究(2001～2005)”的样板地段(封面照片)。

Numerical Study on Rainfall Infiltration in Rock-Soil Blocks

A mathematical model for the rain infiltration in the rock-soil slop has been established and solved by using the finite element method. The unsteady water infiltrating process has been simulated to get water content both in the homogeneous and heterogeneous media. The simulated results show that the rock blocks in the rock-soil slop can cause the wetting front moving fast. If the rain intensity is increased, the saturated region will be formed quickly while other conditions are the same. If the rain intensity keeps a constant, it is possible to accelerate the generation of the saturated region by properly increasing the vertical filtration rate of the rock-soil slop. However, if the vertical filtration rate is so far greater than the rain intensity, it will be difficult to form the saturated region in the rock-soil slop. The numerical method was verified by comparing the calculation results with the field test data.

潘洛铁矿滑坡特征分析及防治

分析了潘洛铁矿露天矿滑坡灾害发生发展的特征,得出该滑坡是一个动态的、多层次的复合型滑坡;找出了影响该边坡稳定的主要因素是大气降雨和采矿切割坡脚;根据这些影响因素和滑坡发展变化规律,提出了科学合理的滑坡治理方案,即动态的综合治理方案.经现场实施后,通过滑坡监测数据表明,该滑坡得到了彻底根治.

The Interaction Between Shock Waves and Foam in a Shock Tube

An experimental study and a numerical simulation were conducted to investigate the mechanical and thermodynamic processes involved in the interaction between shock waves and low density foam. The experiment was done in a stainless shock tube (80mm in inner diameter, 10mm in wall thickness and 5360mm in length). The velocities of the incident and reflected compression waves in the foam were measured by using piezo-ceramic pressure sensors. The end-wall peak pressure behind the reflected wave in the foam was measured by using a crystal piezoelectric sensor. It is suggested that the high end-wall pressure may be caused by a rapid contact between the foam and the end-wall surface. Both open-cell and closed-cell foams with different length and density were tested. Through comparing the numerical and experimental end-wall pressure, the permeability coefficients a and 0 are quantitatively determined.

K418与42CrMo异种金属的激光穿透焊接

实验研究了连续波Nd：YAG激光焊接速度、侧吹保护气流量和离焦量等参量对激光穿透焊接K418和42CrMo焊缝成形的影响。结果表明，K418与42CrMo激光穿透焊接有X形和T形两种典型的焊缝形貌，且焊缝形貌是不对称的。随着焊接速度的提高，焊接线能量降低，焊缝尺寸变小，且焊缝上部尺寸变化比下部尺寸变化慢，焊缝形貌由X形过渡到T形。当离焦量在瑞利长度范围内时，焊缝正面宽度变化很小；当离焦量超出瑞利长度范围时，在足够高的激光功率密度下，焊缝正面宽度快速增加。在激光功率为3kW，侧吹保护气角度为35°条件下，通过优化焊接速度、侧吹保护气流量和离焦量等参量可以得到最佳焊缝质量。

由残余变形反推微构件材料的应力应变关系

通过对微结构破坏模式的测量，反推出微构件材料的应力应变关系。在完成这一分析过程中，首先通过建模与分析，获得了刻划微结构力学行为的几何及物理参数关系，然后结合测量出的数据数值求解得到表征应力应变曲线的特征参量，进而获得应力应变曲线。结合薄膜／基体微结构进行分析。

利用广义应变比能进行结构损伤识别的数值研究

目前结构损伤识别大多是基于位移模态的试验结果进行的，但由于位移模态结果的局限性，很多情况下识别效果并不明显，本文基于应变模态的工作，利用应变模记识别得到的振型和频率结果，提出了广义应变比能（GSED）的概念。将其应用于结构损伤的识别中，具有识别效果明显，物理意义清晰的优点，GSED方法不仅可以判断损伤的存在与否，而且可以对损伤进行较精确的定位，在测试数据充分的情况下还可以大致判断出损伤的程度，最后通过对简支梁的数值仿真计算验证了方法的有效性。

土壤传热特性的实验研究(Experimental Study of Soil Heat Diffusivity)

发展了测定实验室土样热扩散率的方法,介绍了研制的实验装置和建议的操作程序。给出的实验结果表明土壤热扩散率随土壤空隙率、含水量和温度等许多参数而变化。

轴对称钝体前缘分离剪切层的非线性行为

实验研究了锐前缘轴对称钝体有缘分离剪切层发展演化的动力学过程，以及分离剪切层与外加声激励之间的非线性作用。结果表明，自然状态下，存在着亚谐频现象，分离剪切层的演化发展是通过涡合并方式实现的。从分离到再附，层流分离剪切层一般经过2次涡合并过程，湍流分离剪切层一般经过5次涡合并过程。混沌动力学分析表明，层流分离剪切层是经过Hopf分岔和周期倍分岔途径进入混沌状态的。

三峡永久船闸高边坡开挖三维离散元数值模拟:3-D Simulation of the Excavation of High Steep Slope of Three-Gorges Permanent Lock by Distinct Element Method

采用面一面接触的三维离散元刚性块体模型,从实测节理面中取出其中的三组,按照其倾向、倾角和节理间距将三峡永久船闸未开挖的区域划分为10~5个离散单元,通过施加力边界条件,给出了与实测初始地应力场接近的数值模拟结果；然后,分4步模拟了永久船闸的开挖过程。计算结果表明：开挖过程会引起节理面出现张开趋势,个别岩体还会沿着节理面滑移。岩体位移的不对称现象较为自然地说明了由节理引起的岩体各向异性特征。

A Numerical Study of Indentation Using Indenters of Different Geometry

Finite element simulation of the Berkovich, Vickers, Knoop, and cone indenters was carried out for the indentation of elastic-plastic material. To fix the semiapex angle of the cone, several rules of equivalence were used and examined. Despite the asymmetry and differences in the stress and strain fields, it was established that for the Berkovich and Vickers indenters, the load-displacement relation can closely be simulated by a single cone indenter having a semiapex angle equal to 70.3degrees in accordance with the rule of the volume equivalence. On the other hand, none of the rules is applicable to the Knoop indenter owing to its great asymmetry. The finite element method developed here is also applicable to layered or gradient materials with slight modifications.