Numerical investigation of particle velocity distributions in aeolian sand transport

Particle velocity distribution in a blowing sand cloud is a reflection of saltation movement of many particles. Numerical analysis is performed for particle velocity distribution with a discrete particle model. The probability distributions of resultant particle velocity in the impact-entrainment process, particle horizontal and vertical velocities at different heights and the vertical velocity of ascending particles are analyzed. The probability distributions of resultant impact and lift-off velocities of saltating particles can be expressed by a log-normal function, and that of impact angle comply with an exponential function. The probability distribution of particle horizontal and vertical velocities at different heights shows a typical single-peak pattern. In the lower part of saltation layer, the particle horizontal velocity distribution is positively skewed. Further analysis shows that the probability density function of the vertical velocity of ascending particles is similar to the right-hand part of a normal distribution function, and a general equation is acquired for the probability density function of non-dimensional vertical velocity of ascending particles which is independent of diameter of saltating particles, wind strength and height. These distributions in the present numerical analysis are consistent with reported experimental results. The present investigation is important for understanding the saltation state in wind-blown sand movement. (C) 2009 Elsevier B.V. All rights reserved.

CFD-DEM simulation of three-dimensional aeolian sand movement

A three-dimensional CFD-DEM model is proposed to investigate the aeolian sand movement. The results show that the mean particle horizontal velocity can be expressed by a power function of heights. The probability distribution of the impact and lift-off velocities of particles can be described by a log-normal function, and that of the impact and lift-off angles can be expressed by an exponential function. The probability distribution of particle horizontal velocity at different heights can be described as a lognormal function, while the probability distribution of longitudinal and vertical velocity can be described as a normal function. The comparison with previous two-dimensional calculations shows that the variations of mean particle horizontal velocity along the heights in two-dimensional and three-dimensional models are similar. However, the mean particle density of the two-dimensional model is larger than that in reality, which will result in the overestimation of sand transportation rate in the two-dimensional calculation. The study also shows that the predicted probability distributions of particle velocities are in good agreement with the experimental results.

Slow isocharged sequence ions with helium collisions: Projectile core dependence

The collisions of the isocharged sequence ions of q=6 (C6+, N6+, O6+, F6+, Ne6+, Ar6+, and Ca6+), q=7 (F7+, Ne7+, S7+, Ar7+, and Ca7+), q=8 (F8+, Ne8+, Ar8+, and Ca8+), q=9 (F9+, Ne9+, Si9+, S9+, Ar9+, and Ca9+) and q=11 (Si11+, Ar11+, and Ca11+) with helium at the same velocities were investigated. The cross-section ratios of the double-electron transfer (DET) to the single-electron capture (SEC) sigma(DET)/sigma(SEC) and the true double-electron capture (TDC) to the double-electron transfer sigma(TDC)/sigma(DET) were measured. It shows that for different ions in an isocharged sequence, the experimental cross-section ratio sigma(DET)/sigma(SEC) varies by a factor of 3. The results confirm that the projectile core is another dominant factor besides the charge state and the collision velocity in slow (0.35-0.49v(0); v(0) denotes the Bohr velocity) highly charged ions (HCIs) with helium collisions. The experimental cross-section ratio sigma(DET)/sigma(SEC) is compared with the extended classical over-barrier model (ECBM) [A. Barany , Nucl. Instrum. Methods Phys. Res. B 9, 397 (1985)], the molecular Coulombic barrier model (MCBM) [A. Niehaus, J. Phys. B 19, 2925 (1986)], and the semiempirical scaling laws (SSL) [N. Selberg , Phys. Rev. A 54, 4127 (1996)]. It also shows that the projectile core properties affect the initial capture probabilities as well as the subsequent relaxation of the projectiles. The experimental cross-section ratio sigma(TDC)/sigma(DET) for those lower isocharged sequences is dramatically affected by the projectile core structure, while for those sufficiently highly isocharged sequences, the autoionization always dominates, hence the cross-section ratio sigma(TDC)/sigma(DET) is always small.

Entrance channel dependence of the correlation function of IMFs in Ar-36+Sn-112,Sn-124 at 35MeV/u

The reduced velocity correlation functions of the Intermediate Mass Fragments (IMFs) were measured in the reactions of Ar-36+ Sn-112,Sn-124 at 35MeV/u. The anti-correlation at small reduced velocities is more pronounced in Ar-36+ Sn-124 system than that in Ar-36+ Sn-112 system. The difference of the correlation functions between the two reactions is mainly contributed by the particle pairs with high momenta. A three-body Coulomb repulsive trajectory code (MENEKA) is employed to calculate the emission time scale of IMFs for-the both systems. The time scale is 150fm/c in the Ar-36+ Sn-112 system and 120fm/c in the Ar-36+ Sn-124 system, respectively. A calculation based on an Isospin dependence Quantum Molecular Dynamics code (IQMD) reveals that the emission time spectrum of IMFs is shifted slightly leftwards in Ar-36+ Sn-124 compared with that in the Ar-16+ Sn-112 system, indicating a shorter emission time scale. Correspondingly, the central density of the hot nuclei decreases faster in Ar-36+ Sn-124 than in Ar-36+ Sn-112

Applicability of the reaction window theory

Based on the molecular Coulombic over barrier model for description of slow ion-atom collisions, the reaction window theory related to projectile velocity is presented briefly. According to the theory, the state-selective differential cross sections of single electron capture in O8+ -H, A(8+) -H, Ar8+-He, Ne10+-He and Ar18+-He collisions at different collision velocities are calculated and compared with experimental results. Calculations are also done for single, double, and triple electron capture in N-15(7+)-Ne collisions at fixed velocity of 0.53 a.u., and are compared with experimental data. It is found that the predictions of the final electronic state distribution of captured electron(s) are in agreement with experimental data, and both theory and experiments show that the widths of the reaction window increase with the projectile velocity. The differential cross sections predicted by the theory are larger for smaller Q-values, vice versa, when compared with experimental data.

Largely deformed dinuclear system formed in F-19+Al-27 dissipative collision

Excitation functions are measured for different charge products of the F-19+(27) Al reaction in the laboratory energy range 110.25-118.75MeV in steps of 250keV at theta(lab) = 57 degrees, 31 degrees and -29 degrees. The coherence rotation angular velocities of the intermediate dinuclear systems formed in the reaction are extracted from the cross section energy autocorrelation functions. Compared the angular velocity extracted from the experimental data with the ones deduced from the sticking limit, it is indicated that a larger deformation of the intermediate dinuclear system exists.

Electron impact excitation rate coefficients of N II ion

This paper calculates the electron impact excitation rate coefficients from the ground term 2s(2)2p(2) P-3 to the excited terms of the 2s(2)2p(2), 2s2p(3), 2s(2)2p3s, 2s(2)2p3p, and 2s(2)2p3d configurations of N II. In the calculations, rnulticonfiguration Dirac-Fork wave functions have been applied to describe the target-ion states and relativistic distorted-wave calculation has been performed to generate fine-structure collision strengths. The collision strengths are then averaged over a Maxwellian distribution of electron velocities in order to generate the effective collision strengths. The calculated rate coefficients are compared with available experimental and theoretical data, and some good agreements are found for the outer shell electron excitations. But for the inner shell electron excitations there are still some differences between the present calculations and available experiments.

近玻尔速度高电荷态离子与表面作用诱发的溅射离子能谱

近玻尔速度高电荷态离子在物质表面诱发的溅射离子能谱研究在理论和实验上尚是一片空白。本工作针对溅射离子能谱测量，设计、建立了一台径向位置灵敏127o柱形静电离子能谱仪。该谱仪结合位置灵敏探测器，克服了传统静电谱仪只能通过扫描电场单次、单能量点取谱的缺陷，实现了分段取谱——单次可获取宽为0.85Ec ≤ E ≤ 1.15Ec的能谱，极大提高了实验测量效率。利用该谱仪，首次获得了0.8至1.8倍玻尔速度Ar7，8，9+离子与金属铍靶和高定向性石墨靶碰撞产生的溅射离子能谱（入射角45o ，出射角135o）。结果表明：1）能谱分布很宽，其峰值位于130eV至600eV之间，远大于线性级联碰撞区的几到十几eV；且峰位随入射离子速度增加向高能方向略有移动。2）铍靶的溅射能谱呈一个极不对称的峰分布——高能部分有很长的尾巴，下降趋势服从E-n分布；其能谱与入射离子速度关联很大，小于vBohr时指数n与线性级联碰撞接近，而大于vBohr时非线性贡献明显。3）石墨靶的能谱峰位与入射离子动能、势能相关；其谱形在上升到最大值后均出现一个很宽的平台，表明该系列碰撞系统中，非线性效应占主导

Migration of neutral solutes by double stepwise gradient elution in capillary electrochromatography

Characteristics of electroosmotic flow (EOF) and the migration of neutral solutes under double stepwise gradient elution in capillary electrochromatography were studied systematically. EOF velocity proved to be the function of operation time changing with the introduction of the second mobile phase. Accordingly, the retention of components also changed. The migration of neutral solutes was studied under the following three situations; A, components eluted when the column was filled only with the first kind of mobile phase; B, solutes eluted still in the first kind of mobile phase while at that time two kinds of mobile phase coexisted in the column and C, samples eluted in the second kind of mobile phase. Equations to describe the retention times of components under these three kinds of conditions were deduced and applied to predict the retention times of 12 aromatic compounds. Relative errors between experimental and calculated values were below 5.0%, which proved the reliability of the equations. In addition, parameters that might affect the retention time of solutes, such as the transferring time of mobile phase vials, the capacity factors of components and EOF velocities two steps were studied systematically (C) 2001 Elsevier Science B.V. All rights reserved.

2-D and 3-D processing and interpretation of multi-fold ground penetrating radar data:a case history from an archaeological site

A 2.5-D and 3-D multi-fold GPR survey was carried out in the Archaeological Park of Aquileia (northern Italy). The primary objective of the study was the identification of targets of potential archaeological interest in an area designated by local archaeological authorities. The second geophysical objective was to test 2-D and 3-D multi-fold methods and to study localised targets of unknown shape and dimensions in hostile soil conditions. Several portions of the acquisition grid were processed in common offset (CO), common shot (CSG) and common mid point (CMP) geometry. An 8×8 m area was studied with orthogonal CMPs thus achieving a 3-D subsurface coverage with azimuthal range limited to two normal components. Coherent noise components were identified in the pre-stack domain and removed by means of FK filtering of CMP records. Stack velocities were obtained from conventional velocity analysis and azimuthal velocity analysis of 3-D pre-stack gathers. Two major discontinuities were identified in the area of study. The deeper one most probably coincides with the paleosol at the base of the layer associated with activities of man in the area in the last 2500 years. This interpretation is in agreement with the results obtained from nearby cores and excavations. The shallow discontinuity is observed in a part of the investigated area and it shows local interruptions with a linear distribution on the grid. Such interruptions may correspond to buried targets of archaeological interest. The prominent enhancement of the subsurface images obtained by means of multi-fold techniques, compared with the relatively poor quality of the conventional single-fold georadar sections, indicates that multi-fold methods are well suited for the application to high resolution studies in archaeology.

侵蚀细沟水力学特性及细沟侵蚀与形态特征的试验研究

本文进行了一系列水槽试验 ,研究不同长度的侵蚀细沟在坡度及入流量影响下的浅水水流水力学特性、侵蚀量及细沟的形态特征 .所采用的土壤材料为砂壤土 ,试验使用四种入流流量 ,四种坡度 ,四种细沟长度 .记录水流速度、沟宽及径流中的泥沙含量 ,并用多元回归方法分析了它们与坡度、入流量或它们的交互项之间的关系 .提出了一种数学模型来描述试验中沟宽的周期性形态变化 .提出了通过试验得出的不同细沟长度的产沙量来确定水流输沙能力的方法 ,并得出出测定输沙能力的采样长度为 2～ 4m ,通过对细沟形态变化分析进一步说明了所确定的水流输沙能力采样长度的正确性

Structures and vibrational spectra of C-2 and LaC2+ clusters

C-2 and LaC2+ were studied using Hartree-Fock(HF), B3LYP (Becke 3-paremeter-Lee-Yang-Parr) density functional method, second-order Moller-Plesset perturbation (MP2) and coupled cluster singles and doubles with non-iterative triples(CCSD(T)) methods. The basis set employed was LANL1DZ. Geometries, vibrational frequencies and other quantities were reported. The results showed that for C-2, all the methods performed well for low spin state (singlet), while only HF and B3LYP remained so for high spin state (triplet). For LaC2+, four isomers were presented and fully optimized. The results suggested that linear isomers with C-infinity v and D-infinity h symmetries were predicted to be saddle points on the energy surface for all the methods, while for isomers with C-2 upsilon and C-s symmetries, they were local minima except C-2 upsilon at B3LYP level, and were isoenergetic at HF, MP2 and CCSD(T) levels, near isoenergetic at B3LYP level. From the differences between HOMO and LUMO, it is also known that the isomers with C-2 upsilon and C-s symmetries offer the largest values and therefore correspond to the most stable structure. For La-C bond lengths, B3LYP gives the shortest, the order is B3LYP

Geomorphic and hydrodynamic responses in salt marsh-tidal creek systems, Jiangsu, China

Salt marsh-tidal creek systems as a coastal geomorphological unit represent an important natural resource. The present study on Jiangsu salt marshes, eastern China, shows that variations in tidal current velocities in salt marsh creeks are controlled by the local tidal wave characteristics and the bed slope and elevation of the salt marshes and creeks. Likewise, the tidal currents modify the geomorphology of the salt marsh-tidal creek systems by transporting sediments and causing erosion/deposition. Storm events, which appear to have cyclical changes in their intensity relating to sunspot activities, can affect the geomorphic evolution of such systems. Further, in response to accelerated sea-level rise, accretional rates on salt marshes may increase. The tidal creeks have the function of transporting water and sediment onto the salt marsh surface; further, the energy of tidal currents and waves are dissipated within the salt marsh-tidal creek system. Hence, this coastal system has a potential value for coastal protection.

Higher-order Boussinesq-type equations for interfacial waves in a two-fluid system

Interfacial waves propagating along the interface between a three-dimensional two-fluid system with a rigid upper boundary and an uneven bottom are considered. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. A set of higher-order Boussinesq-type equations in terms of the depth-averaged velocities accounting for stronger nonlinearity are derived. When the small parameter measuring frequency dispersion keeping up to lower-order and full nonlinearity are considered, the equations include the Choi and Camassa's results (1999). The enhanced equations in terms of the depth-averaged velocities are obtained by applying the enhancement technique introduced by Madsen et al. (1991) and Schaffer and Madsen (1995a). It is noted that the equations derived from the present study include, as special cases, those obtained by Madsen and Schaffer (1998). By comparison with the dispersion relation of the linear Stokes waves, we found that the dispersion relation is more improved than Choi and Camassa's (1999) results, and the applicable scope of water depth is deeper.

Strongly nonlinear long interfacial waves in uniform basic flows and their dispersion relationship

In this paper, long interfacial waves of finite amplitude in uniform basic flows are considered with the assumption that the aspect ratio between wavelength and water depth is small. A new model is derived using the velocities at arbitrary distances from the still water level as the velocity variables instead of the commonly used depth-averaged velocities. This significantly improves the dispersion properties and makes them applicable to a wider range of water depths. Since its derivation requires no assumption on wave amplitude, the model thus can be used to describe waves with arbitrary amplitude.