Application of the finite Hankel transform to a diffusion problem without azimuthal symmetry

The problem treated in this paper concerns calculating the evolution of the pressure in a single-phase, slightly compressible fluid in a porous medium consisting of communicating layers. The fluid is produced through a point sink located on the side of an otherwise sealed cylindrical wellbore. This location of the sink causes the flow around the wellbore to be azimuthally asymmetric.The problem is solved through successive application of Laplace, finite Fourier and finite Hankel transforms. Although apparently straightforward, this approach leads to serious numerical difficulties. The published form of the inversion formula for the finite Hankel transform leads to inaccurate computation for the higher azimuthal modes even with 128 bit arithmetic. An alternative form is developed which enables accurate evaluation of the solution with the more practical 64 bit arithmetic. The technique for two-layer solution presented here can be directly extended to a problem with a larger number of communicating layers. This is the first instance of successful application of the finite Hankel transform to an azimuthally asymmetric diffusion problem.     

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.     

Deactivation reactions in the modeled 2,2,6,6‐tetramethyl‐1‐piperidinyloxy‐mediated free‐radical polymerization of styrene: A comparative study with the 2,2,6,6‐tetramethyl‐1‐piperidinyloxy/acrylonitrile system

The competitiveness of the combination and disproportionation reactions between a 1‐phenylpropyl radical, standing for a growing polystyryl macroradical, and a 2,2,6,6‐tetramethyl‐1‐piperidinyloxy (TEMPO) radical in the nitroxide‐mediated free‐radical polymerization of styrene was quantitatively evaluated by the study of the transition geometry and the potential energy profiles for the competing reactions with the use of quantum‐mechanical calculations at the density functional theory (DFT) UB3‐LYP/6‐311+G(3df, 2p)//(unrestricted) Austin Model 1 level of theory. The search for transition geometries resulted in six and two transition structures for the radical combination and disproportionation reactions, respectively. The former transition structures, mainly differing in the out‐of‐plane angle of the N? O bond in the transition structure TEMPO molecule, were correlated with the activation energy, which was determined to be in the range of 8.4–19.4 kcal mol^?1 from a single‐point calculation at the DFT UB3‐LYP/6‐311+G(3df, 2p)//unrestricted Austin Model 1 level. The calculated activation energy for the disproportionation reaction was less favorable by a value of more than 30 kcal mol^?1 in comparison with that for the combination reaction. The approximate barrier difference for the TEMPO addition and disproportionation reaction was slightly smaller for the styrene polymerization system than for the acrylonitrile polymerization system, thus indicating that a β‐proton abstraction through a TEMPO radical from the polymer backbone could diminish control over the radical polymerization of styrene with the nitroxide even more than in the latter system. © 2006 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 45: 232–241, 2007     

虚拟技术在实验教学中的应用

介绍了虚拟技术在实验教学改革中的作用与意义，探讨了电子信息系统实验教学改革的新思路与实现措施，并总结了实验教学改革的初步成果．     

Monte Carlo Simulation of Electron Transport and X-Ray Generation. II. Radiative Processes and Examples in Electron Probe Microanalysis

Theoretical models for Monte Carlo simulation of radiative processes, i.e. bremsstrahlung and characteristic x-ray emission, are presented. Possible strategies for simulating electron transport are briefly described. For mechanisms involving energy loss and angular deflections, difficulties for strict implementation of accurate numerical differential cross sections still remain due to the strong correlations between these variables. Practical solutions for the case of inelastic collisions and bremsstrahlung emission are described. Comparisons of simulation results with experimental data for several problems of interest in electron probe microanalysis are presented.     

Software Cold Test Simulation of Coupled Cavity Slow-Wave Structure in Millimeter Wave TWT

In this paper, the software cold-test simulation method to obtain RF dispersion, interaction impedance by MAFIA software is discussed considering effect of the conducting barrel in the millimeter-wave TWT. The dispersion simulation result is obtained more consistent with experiment data than the method that ignoring effect of the conducting barrel. By changing the structure parameter, high frequency characteristic will be affected. The change relationship of high frequency characteristic with structure parameter is obtained. The different calculation methods of the interaction impedance are discussed. These results are received and are consistent with experiment data.