Role of Ponded Turbidity Currents in Reservoir Trap Efficiency

The capacity to store water in a reservoir declines as it traps sediment. A river entering a reservoir forms a prograding delta. Coarse sediment (e.g., sand) deposits in the fluvial topset and avalanching foreset of the delta, and is typically trapped with an efficiency near 100%. The trap efficiency of fine sediment (e.g., mud), on the other hand, may be below 100%, because some of this sediment may pass out of the reservoir without settling out. Here, a model of trap efficiency of mud is developed in terms of the mechanics of a turbidity current that plunges on the foreset. The dam causes a sustained turbidity current to reflect and form a muddy pond bounded upstream by a hydraulic jump. If the interface of this muddy pond rises above any vent or overflow point at the dam, the trap efficiency of mud drops below 100%. A model of the coevolution of topset, foreset, and bottomset in a reservoir that captures the dynamics of the internal muddy pond is presented. Numerical implementation, comparison against an experiment, and application to a field-scale case provide the basis for a physical understanding of the processes that determine reservoir trap efficiency.     

Cellular-Automata Model for Dense-Snow Avalanches

This paper introduces a three-dimensional model for simulating dense-snow avalanches, based on the numerical method of cellular automata. This method allows one to study the complex behavior of the avalanche by dividing it into small elements, whose interaction is described by simple laws, obtaining a reduction of the computational power needed to perform a three-dimensional simulation. Similar models by several authors have been used to model rock avalanches, mud and lava flows, and debris avalanches. A peculiar aspect of avalanche dynamics, i.e., the mechanisms of erosion of the snowpack and deposition of material from the avalanche is taken into account in the model. The capability of the proposed approach has been illustrated by modeling three documented avalanches that occurred in Susa Valley (Western Italian Alps). Despite the qualitative observations used for calibration, the proposed method is able to reproduce the correct three-dimensional avalanche path, using a digital terrain model, and the order of magnitude of the avalanche deposit volume.     

Practical Equivalent Continuum Model for Simulation of Jointed Rock Mass Using FLAC3D

A simple practical equivalent continuum numerical model for simulating the behavior of jointed rock mass has been extended to three-dimensional using FLAC3D. This model estimates the properties of jointed rock mass from the properties of intact rock and a joint factor (Jf), which is the integration of the properties of joints to take care of the effects of frequency, orientation, and strength of joint. A new FISH function has been written in FLAC3D specifically for modeling jointed rocks using the Duncan and Chang hyperbolic model. This model has been validated first with simple element tests at different confining pressures for different rocks with different joint configurations. Explicit modeling of the joints has also been done in element tests and results obtained compare well with the results of equivalent continuum model and also with experimental results. Further, this has been applied for a case study of a large underground power house cavern in the Himalayas. The analysis has been done under various stages of excavation, assigning a null model available in FLAC3D for simulating the excavation.     

Two-Dimensional Physical and Numerical Modeling of a Pile-Supported Earth Platform over Soft Soil

This paper focuses on the mechanisms occurring in a granular earth platform over soft ground improved by rigid piles. Two-dimensional physical model experiments were performed using the Schneebeli’s analogical soil to investigate the load transfer mechanisms by arching and the settlement reduction and homogenization. Experimental outputs are compared to results obtained on a numerical model using a plane strain continuum approach. The impact of the constitutive model complexity to simulate the platform material behavior was first assessed, but no large difference was recorded. As far as the proposed model, which takes the main features of the observed behavior satisfactorily into account, the numerical procedure could be validated and the parametric studies extended numerically. Both approaches of this study underlined the main geometrical and geotechnical parameters which should inevitably be taken into account in a simplified design method, namely the capping ratio, the platform height, and the platform material shear strength.     

Applications of the discrete element method in mechanical engineering

Compared to other fields of engineering, in mechanical engineering, the Discrete Element Method (DEM) is not yet a well known method. Nevertheless, there is a variety of simulation problems where the method has obvious advantages due to its meshless nature. For problems where several free bodies can collide and break after having been largely deformed, the DEM is the method of choice. Neighborhood search and collision detection between bodies as well as the separation of large solids into smaller particles are naturally incorporated in the method. The main DEM algorithm consists of a relatively simple loop that basically contains the three substeps contact detection, force computation and integration. However, there exists a large variety of different algorithms to choose the substeps to compose the optimal method for a given problem. In this contribution, we describe the dynamics of particle systems together with appropriate numerical integration schemes and give an overview over different types of particle interactions that can be composed to adapt the method to fit to a given simulation problem. Surface triangulations are used to model complicated, non-convex bodies in contact with particle systems. The capabilities of the method are finally demonstrated by means of application examples. Commemorative Contribution.     

熔窑玻璃液流数值模拟的新发展

熔窑玻璃液流数值模拟的新发展赵国昌（清华大学核能技术设计研究院，１０００８４）ＮｅｗＤｅｖｅｌｏｐｍｅｎｔｏｆｔｈｅＮｕｍｅｒｉｃａｌＳｉｍｕｌａｔｉｏｎｆｏｒＦｌｏｗａｎｄＨｅａｔＴｒａｎｓｆｅｒｏｆＧｌａｓｓＭｅｌｔｉｎＴａｎｋＦｕｒｎａｃｅ￥Ｚ...     

RADIATION FIELD FROM A RECTANGULAR PLAQUE SOURCE FOR RADIATION PROCESSING

Evaluation of the Hubbell integral for an omnidirectional detector response to a rectangular plaque gamma emitting radiation source using a different numerical computation method on a personal computer was made. For a plaque consisting of several strip gamma emitters, exposure rates at different positions above the source were calculated using 'exact' evaluation and an approximate method. Optimization method is used to find source activity of each strip, so as to acquire desirable dose uniformity of the irradiated products. In uses of the Hubbell integral, by introducing appropriate programming, computation time can be saved.     

Comparison of Displacement and Mixing Diffusers

In order to cumpare the peformance of different supply diffuers of ventilation air, the airflow passern, temperature stratifiation and contaminant dispersion in a furnitured office ventilated by three kinds of air diffuer were numerically investigated. The air diffuers studied in this paper are a quarter-cylinder displacement diffuer on the floor and mixing diffuers (linear and vortex diffuers) on the ceiling. The heat sources in the of-fice are considered to be 50% convective and 50% radiative. The k-? two-equatwn model of turbulence is employed to predict the turbulent diffusion. The results show that the displacement diffuser provides a rather uniform flow field with low velocify in most areas, and the vertical temperature difference from floor to ceiling is as high as 6 K. With the linear diffuser, the air velociry is high, and the temperature is uniform both horizontally and vertically. The air velocity generated by the vortex diffuser is moderate. The distributions of the temperature and the contaminant are rather uniform.     

平面分层介质结构全波分析综述

本文综述了广义谱域导抗法，直线法及直线法的快速算法，全波离散镜象法等主要方法，并对各种方法的特点进行了比较，这些方法可有效地处理平面分层介质结构问题，文中包括作者近期的研究成果。     

A numerical model for the dissolution of spherical particles in binary alloys under mixed mode control

A general numerical model is described for the dissolution kinetics of spherical particles in binary systems for any combination of first order reactions at the particle-matrix interface and long distance diffusion within the matrix. The model is applicable to both finite and infinite media and handles both complete and partial particle dissolution. It is shown that interfacial reactions can have a strong effect on the dissolution kinetics, the solute concentration at the particle-matrix interface and the solute concentration profile in the matrix.