Modelling raindrop impact and splash erosion processes within a spatial cell: a stochastic approach

A new approach is proposed to simulate splash erosion on local soil surfaces. Without the effect of wind and other raindrops, the impact of free-falling raindrops was considered as an independent event from the stochastic viewpoint. The erosivity of a single raindrop depending on its kinetic energy was computed by an empirical relationship in which the kinetic energy was expressed as a power function of the equivalent diameter of the raindrop. An empirical linear function combining the kinetic energy and soil shear strength was used to estimate the impacted amount of soil particles by a single raindrop. Considering an ideal local soil surface with size of I m x I m, the expected number of received free-failing raindrops with different diameters per unit time was described by the combination of the raindrop size distribution function and the terminal velocity of raindrops. The total splash amount was seen as the sum of the impact amount by all raindrops in the rainfall event. The total splash amount per unit time was subdivided into three different components, including net splash amount, single impact amount and re-detachment amount. The re-detachment amount was obtained by a spatial geometric probability derived using the Poisson function in which overlapped impacted areas were considered. The net splash amount was defined as the mass of soil particles collected outside the splash dish. It was estimated by another spatial geometric probability in which the average splashed distance related to the median grain size of soil and effects of other impacted soil particles and other free-falling raindrops were considered. Splash experiments in artificial rainfall were carried out to validate the availability and accuracy of the model. Our simulated results suggested that the net splash amount and re-detachment amount were small parts of the total splash amount. Their proportions were 0.15% and 2.6%, respectively. The comparison of simulated data with measured data showed that this model could be applied to simulate the soil-splash process successfully and needed information of the rainfall intensity and original soil properties including initial bulk intensity, water content, median grain size and some empirical constants related to the soil surface shear strength, the raindrop size distribution function and the average splashed distance. Copyright (c) 2007 John Wiley & Sons, Ltd.

The electrochemical behavior of a system with a limited number of molecules

In this paper, based on Einstein relationship between diffusion and random walk, the electrochemical behavior of a system with a limited number of molecules was simulated and explored theoretically. The transition of the current vs time responses from discrete to continuous was clearly obtained as the number of redox molecules increased from 10 to 10(6).

What has brought on the effects of number-average molecular weight on the spinodals in polymer mixtures?

On the basis of the thermodynamics of Gibbs, the spinodal for the quasibinary system was derived in the framework of the Sanchez-Lacombe lattice fluid theory. All of the spinodals were calculated based on a model polydisperse polymer mixture, where each polymer contains three different molecular weight subcomponents. According to our calculations, the spinodal depends on both weight-average ((M) over bar (w)) and number-average ((M) over bar (n)) molecular weights, whereas that of the z-average molecular weight is invisible. Moreover, the extreme of the spinodal decreases when the polydispersity index (eta = (M) over bar (w)/(M) over bar (n)) of the polymer increases. The effect of polydispersity on the spinodal decreases when the molecular weight gets larger and can be negligible at a certain large molecular weight. It is well-known that the influence of polydispersity on the phase equilibrium (coexisting curve, cloud point curves) is much more pronounced than on the spinodal. The effect of M, on the spinodal is discussed as it results from the infuluence of composition temperatures, molecular weight, and the latter's distribution on free volume. An approximate expression, which is in the assumptions of v* v(1)* = v(2)* and 1/r --> 0 for both of the polymers, was also derived for simplification. It can be used in high molecular weight, although it failed to make visible the effect of number-average molecular weight on the spinodal.

A novel electrochemical method for simultaneously measuring diffusion coefficient and ion transference number

This paper introduces a new method to estimate the diffusion coefficient and transference number of a salt or an electroactive ion in a solution with little or no supporting electrolyte. The above two parameters can be obtained from a single potential step experiment without previous knowledge of either one. It would appear that the method could also be used in the study of ion transport in a high viscosity solvent or a solid electrolyte. (C) 1998 Elsevier Science S.A.

Wave generation by the falling rock in the two-dimensional wave tank

Wave generation by the falling rock in the two-dimensional wave tank is experimentally and numerically studied, where the numerical model utilizes the boundary element method to solve the fully nonlinear potential flow theory. The wave profiles at different times are measured in the laboratory, which are also used to test the numerical model. Comparisons show that the experimental and numerical results are in good agreement, and the numerical model can be used to simulate the wave generation due to the submarine rock falling. Further numerical tests on the influences of the rock size, density, initial position and the falling angle on the wave elevation of the generated waves are performed, respectively. The results show that the size and density of the rock have strong effects on the maximum elevation of the generated wave, while the effects of the initial position and the falling angle of the rock are also significant. When the size or the density of the rock increases, the maximum elevation of the generated wave increases. The same effect on the generated wave would be produced if the initial position of the rock becomes closer to the surface, or the falling angle between the falling route and the vertical direction turns larger. In addition, the present numerical tests reveal that the submarine rock falling provides a new generation method for the breaking wave in the wave tank.

Outer wave number spectrum and statistical distribution of apparent wave number vector

Based on the ray theory and Longuet-Higgins's linear,model of sea waves, the joint distribution of wave envelope and apparent wave number vector is established. From the joint distribution, we define a new concept, namely the outer wave number spectrum, to describe the outer characteristics of ocean waves. The analytical form of the outer wave number spectrum, the probability distributions of the apparent wave number vector and its components are then derived. The outer wave number spectrum is compared with the inner wave number spectrum for the average status of wind-wave development corresponding to a peakness factor P = 3. Discussions on the similarity and difference between the outer wave number spectrum and inner one are also presented in the paper. (C) 2002 Elsevier Science Ltd. All rights reserved.

边坡动力响应分析及应用研究

As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum VHlhro) point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fN

Effect of pore ring number of zeolites on catalytic performance of Mo/Zeolite in methane aromatization under non-oxygen condition

The catalytic behavior of Mo-based zeolite catalysts with different pore structure and size, particularly with 8 membered ring ( M R), 10 M R, coexisted 10 and 12 M R, and 12 M R, was studied in methane aromatization under the conditions of SV=1500 ml/(g.h), p=0.1 MPa and T = 973 K. It was found that the catalytic performance is correlated with the pore structure of the zeolite supports. The zeolites that possess 10 MR or 10 and 12 MR pore structure with a pore diameter equal to or slightly larger than the dynamic diameter of benzene molecule, such as ZSM-5, ZSM-11, ZRP-1 and MCM-22, are fine supports. Among the tested zeolite supports, MCM-22 exhibits the highest activity and selectivity for benzene. A methane conversion of 10.5% with benzene selectivity of 80% was achieved over Mo/MCM-22 catalyst. The Mo/ERS-7 catalyst with 8 MR (0.45 nm) does not show any activity in methane dehydro-aromatization, while Mo/JQX-1 and Mo/SBA-15 catalysts with 12 MR pore exhibit little activity in the reaction. It can be concluded that the zeolites with 10 MR pore or coexisted 10 and 12 MR, having pore size equal to or slightly larger than the dynamic diameter of benzene molecule, are fine supports for methane activation and aromatization.