声流现象的研究及其应用

本文从声流产生的机理出发，综合了影响声流的诸因素，介绍了声流现象的检测方法及声流效应的应用。     

作物覆盖条件下土壤中水热分布的稳态数值比较

本文应用非饱和多孔介质中传热传质的数学模型,考虑到作物根系和冠层的影响,对稳态情况下圆柱形士壤床中的水热分布进行了数值计算模拟。在不同环境条件下对土壤中的水热分布也进行了模拟计算。同时,在相同环境条件下对有作物覆盖的土壤床和无作物覆盖的土壤床中的水热分布也进行了模拟比较。文中数值模拟的结果较好地反映了实际情况。     

Velocity spectrum for non-Markovian Brownian motion in a periodic potential

Non-Markovian Brownian motion in a periodic potential is studied by means of an electronic analogue simulator. Velocity spectra, the Fourier transforms of velocity autocorrelation functions, are obtained for three types of random force, that is, a white noise, an Ornstein—Uhlenbeck process, and a quasimonochromatic noise. The analogue results are in good agreement both with theoretical ones calculated with the use of a matrix-continued-fraction method, and with the results of digital simulations. An unexpected extra peak in the velocity spectrum is observed for Ornstein-Uhlenbeck noise with large correlation time. The peak is attributed to a slow oscillatory motion of the Brownian particle as it moves back and forth over several lattice spaces. Its relationship to an approximate Langevin equation is discussed.     

系统动力学模型结构分析中表函数的改进方法

分析了非线性系统动力学(SD)模型中表函数环节的传统表示方法,指出该方法对于进行SD模型的结构分析所造成的障碍,并提出了用牛顿插值的承袭性算法来改进表函数环节的表示方法的思路.接着分别介绍了用差商和差分牛顿插值法改进表函数的过程.最后,通过一个实例说明了表函数改进的方法,并展示了以此为基础进行结构分析的思路.     

EAST装置等离子体放电位形快速识别研究

为了适应超导托卡马克装置EAST位形控制运行模式的需要，研究了等离子体放电位形快速识别算法，给出了等离子体放电位形重建的模拟计算结果，并与美国GA平衡反演程序EFIT计算的平衡结果进行了详细的误差对比分析.结果表明，通过外部磁测量和合适的电流剖面模型，并结合实时平衡重建算法可以快速地对EAST等离子体放电位形进行识别. 关键词： 托卡马克 等离子体 平衡重建 数值模拟     

Influence of Relative Permeability on the Stability Characteristics of Immiscible Flow in Porous Media

Relative permeability functions for immiscible displacements in porous media show a wide range of profiles. Although, this behavior is well known, its impact on the stability of the displacement process is unexplored. Our analysis clearly demonstrates for the first time that the viscous instability characteristics of two-phase flows are governed not only by their end point values, but are strongly dependent on the actual profile of relative permeability functions. Linear stability analysis predicts the capacity of the flow to develop large scale fingers which can result in substantial bypassing of the resident fluid. It is observed that relative permeability functions attributed to drainage processes yield a more unstable displacement as compared to functions related to imbibition processes. Moreover, instability is observed to increase for those relative permeability functions which result from increased wettability of the wetting fluid. High accuracy numerical simulations show agreement with these predictions and demonstrate how large amplitude viscous fingers result in significant bypassing for certain relative permeability functions. In the nonlinear regime, the finger amplitude grows at a rate ∝ t^1/2 initially, drops to t^1/4 at a later time and finally grows ∝ t. The basic mechanisms of finger interaction, however, are not substantially influenced by relative permeability functions.     

Measurements and simulations of particle dispersion in a turbulent flow

A particle imaging technique has been used to collect droplet displacement statistics in a round turbulent jet of air. Droplets are injected on the jet axis, and a laser sheet and position-sensitive photomultiplier tube are used to track their radial displacement and time-of-flight. Dispersion statistics can be computed which are Lagrangian or Eulerian in nature. The experiments have been simulated numerically using a second-order closure scheme for the jet and a stochastic simulation for the particle trajectories. Results are presented for non-vaporizing droplets of sizes from 35 to 160 μm. The simulations have underscored the importance of initial conditions and early droplet displacement history on the droplet trajectory for droplets with large inertia relative to the turbulence. Estimates of initial conditions have been made and their effect on dispersion is quantified.     

Equivariant mappings: a new approach in stochastic simulations

Stochastic simulations on manifolds usually are traced back to ⁿ via charts. If a group G is acting on a manifold M and if the respective distribution v is invariant under this group action then in many cases of practical interest there exists a more convenient approach which uses equivariant mappings. The concept of equivariant mappings will be discussed intensively at the instance of the Grassman manifold in which case G equals the orthogonal group. Further advantages of this concept will be demonstrated by applying it to a probabilistic problem from the field of combinatorial geometry.     

硼中子俘获治疗的蒙特卡罗方法模拟

用通用蒙特卡罗程序MCNP模拟了粒子在人脑中的输运过程. 吸收剂量率主要来自以下四个反应:¹⁰B(n,α)⁷Li,¹⁴N(n,p)¹⁴C,¹H(n,γ)²D,快中子弹性散射反应.对肿瘤区的贡献主要来自硼中子吸收反应.结果表明,超热中子比热中子适合于深肿瘤的治疗,而热中子对浅肿瘤的治疗有优越性,比如皮肤癌.同确定论方法的结果相比,蒙特卡罗方法不失为一种模拟中子俘获治疗的好工具.     

A recursive-faulting model of distributed damage in confined brittle materials

We develop a model of distributed damage in brittle materials deforming in triaxial compression based on the explicit construction of special microstructures obtained by recursive faulting. The model aims to predict the effective or macroscopic behavior of the material from its elastic and fracture properties; and to predict the microstructures underlying the microscopic behavior. The model accounts for the elasticity of the matrix, fault nucleation and the cohesive and frictional behavior of the faults. We analyze the resulting quasistatic boundary value problem and determine the relaxation of the potential energy, which describes the macroscopic material behavior averaged over all possible fine-scale structures. Finally, we present numerical calculations of the dynamic multi-axial compression experiments on sintered aluminum nitride of Chen and Ravichandran [1994. Dynamic compressive behavior of ceramics under lateral confinement. J. Phys. IV 4, 177-182; 1996a. Static and dynamic compressive behavior of aluminum nitride under moderate confinement. J. Am. Soc. Ceramics 79(3), 579-584; 1996b. An experimental technique for imposing dynamic multiaxial compression with mechanical confinement. Exp. Mech. 36(2), 155-158; 2000. Failure mode transition in ceramics under dynamic multiaxial compression. Int. J. Fracture 101, 141-159]. The model correctly predicts the general trends regarding the observed damage patterns; and the brittle-to-ductile transition resulting under increasing confinement.