三维射频加速腔程序包MAFIA的开发和应用

介绍了三维射频加速器计算机辅助设计(CAD)程序包MAFIA的特点、功能及其应用。用CAD的方法实现了射频加速器特别是非轴对称带耦合腔的加速器的三维腔形设计。同时介绍了节省内存、节约机时的模型方法,使得在中小型机器上运行大程序(八万条)成为可能。文章主要对L波段的加速腔进行了数值模拟,并和实验进行了比较,结果符合较好;通过三维数值模拟及图形模拟的方法比较了各种形状谐振腔的性能;计算分析了输入耦合腔及高阶模耦合器对加速腔电磁场分布等参量的影响;提出了加耦合器后谐振频率改变的补偿方法。     

斜孔气膜冷却数值模拟分析

本文通过数值模拟分析圆柱孔和扩散孔的单斜孔气膜冷却特性,考察复合角、孔型和吹风比对流场和气膜冷却效果的影响。结果表明,复合角的引入使气膜侧向分布更宽,但冷却效果沿主流方向衰减更快。适中吹风比得到的气膜能更有效的保护壁面。在相同吹风比和复合角条件下,扩散孔的气膜冷却效率比圆柱孔更好,且冷却更为均匀持久。     

Molecular dynamics simulation of entropy and surface tension for grain boundary of α-Fe

The grain boundary is an interface and the surface tension is one of its important thermodynamic properties. In this paper, the surface tension of the ∑9 grain boundary for α-Fe at various temperatures and pressures is calculated by means of Computer Molecular Dynamics (CMD). The results agree satisfactorily with the experimental data. It is shown that the contribution of entropy to surface tension of grain boundary can be ignored. The project supported by the National Natural Science Foundation of China and the Science Foundation of Chinese Academy of Sciences.     

CdZnTe核探测器的蒙特卡罗模拟的初步研究

以CdZnTe核探测器的工作原理为依据,探测器内反应的随机性和反应产生的电子空穴对数目的统计规律为物理模型,应用Visual C + + 自行编制了蒙特卡罗模拟软件.模拟了γ射线在CdZnTe探测器中的响应能谱,并将模拟结果与实际器件的测试结果进行了比较讨论.模拟能谱与实际测得的能谱的主峰符合较好.此外,通过分析⁵⁷Co源辐照下探测效率与器件厚度的关系,可以推测探测效率达到最大时所对应CdZnTe探测器的理想厚度     

影响产生均匀等离子体状态的几个因素

 用数值模拟方法研究利用辐射加热来产生均匀等离子体状态,它可被用来测量元素的辐射不透明度,校验辐射不透明度理论。研究了在辐射加热铁的“三明治”型靶时影响生成均匀等离子体状态的几个重要因素（样品厚度、CH膜厚度和辐射源）所起的作用，研究发现，当铁等离子体通过热传导达到均匀状态时，其尺度必须与此时的传热距离相当，从而定出铁样品的厚度；低Z介质CH膜对铁等离子体有明显的箍束作用，调整CH膜的厚度可以调节所产生的等离子体状态；调整不同方向上的CH膜厚度，可以控制铁等离子体的膨胀方向，使它尽可能地达到一维膨胀，使得反推出的等离子体密度可以更加准确；样品的种类、厚度以及外面的低Z介质厚度决定了在某时刻能获得的等离子体状态，以及为产生此等离子体状态所需的最低辐射能。     

Simplified analytic solutions and a novel fast algorithm for Yb-doped double-clad fiber lasers

In this paper, continuous wave Yb³⁺-doped double-clad fiber lasers (DCFLs) with linear-cavity are investigated theoretically and numerically using the rate equations. Under the steady state conditions, the simplified analytic solutions of Yb³⁺-doped DCFLs under considering the scattering loss are deduced in the strongly pump condition. Compared with the known analytic solutions in published literatures, our analytic solutions are more accurate, especially, at higher reflectivity of output mirror. In addition, a fast and stable algorithm based on the Newton-Raphson method is proposed to simulate numerically Yb³⁺-doped DCFLs. The results by simplified analytic solutions are in good agreement with those by the numerical simulation. Moreover, we have performed the optimization of an Yb³⁺-doped DCFL using the simplified analytic solutions and the numerical simulations, respectively.     

Numerical simulation of soft palate movement and airflow in human upper airway by fluid-structure interaction method

In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome (OSAHS), three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed three-dimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosing diseases related to anatomical structure and function of the upper airway. The project supported by the National Natural Science Foundation of China (10672036, 10472025 and 10421002), the Natural Science Foundation of Liaoning Province (20032109). English text was polished by Yunming Chen.     

Estimating tail probabilities of heavy tailed distributions with asymptotically zero relative error

Efficient estimation of tail probabilities involving heavy tailed random variables is amongst the most challenging problems in Monte-Carlo simulation. In the last few years, applied probabilists have achieved considerable success in developing efficient algorithms for some such simple but fundamental tail probabilities. Usually, unbiased importance sampling estimators of such tail probabilities are developed and it is proved that these estimators are asymptotically efficient or even possess the desirable bounded relative error property. In this paper, as an illustration, we consider a simple tail probability involving geometric sums of heavy tailed random variables. This is useful in estimating the probability of large delays in M/G/1 queues. In this setting we develop an unbiased estimator whose relative error decreases to zero asymptotically. The key idea is to decompose the probability of interest into a known dominant component and an unknown small component. Simulation then focuses on estimating the latter ‘residual’ probability. Here we show that the existing conditioning methods or importance sampling methods are not effective in estimating the residual probability while an appropriate combination of the two estimates it with bounded relative error. As a further illustration of the proposed ideas, we apply them to develop an estimator for the probability of large delays in stochastic activity networks that has an asymptotically zero relative error.      

双对数模型对模型模拟误差的放缩问题探讨

对双对数模型lg Y=a0+a1lg X1+a2lg X2+…+anlg Xn与其对应的指数模型y=c0xa11xa22…xann的模拟相对误差的关系进行了探讨,指出双对数模型具有放大和缩小指数模型相对误差的特性.对二者的关系进行了理论推导和实例验证,并给出了二者的定量关系式.     

LES-DSMC方法研究超细颗粒气固两相流动过程

采用考虑颗粒脉动流动对气相湍流流动影响的大涡模拟(LES)研究气相湍流,采用直接模拟蒙特卡罗方法(DSMC)模拟颗粒间的碰撞。单颗粒运动满足牛顿第二定律,颗粒相和气相相间作用的双向耦合由牛顿第三定律确定,考虑超细颗粒间的van der Waals作用力。数值模拟垂直管内超细颗粒气固两相流动,对颗粒相速度、浓度以及团聚物流动过程进行分析。