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的模拟相对误差的关系进行了探讨,指出双对数模型具有放大和缩小指数模型相对误差的特性.对二者的关系进行了理论推导和实例验证,并给出了二者的定量关系式.     

Memory management optimization problems for integrated circuit simulators

In hardware design, it is necessary to simulate the anticipated behavior of the integrated circuit before it is actually cast in silicon. As simulation procedures are long due to the great number of tests to be performed, optimization of the simulation code is of prime importance. This paper describes two mathematical models for the minimization of the memory access times for a cycle-based simulator.An integrated circuit being viewed as a directed acyclic graph, the problem consists in building a graph order on the vertices, compatible with the relation order induced by the graph, in order to minimize a cost function that represents the memory access time. For both proposed cost functions, we show that the corresponding problems are NP-complete. However, we show that the special cases where the graphs are in-trees or out-trees can be solved in polynomial time.     

钱塘江河口弯道潮流场的数值模拟

钱塘江河口潮流错综复杂．依靠零星的实测潮流资料无法全面描述其流场特征．本文采用Dcldt3D模型对钱塘江河口闻家堰至仓前河段的潮流场进行二维数值模拟．着重研究了河口弯道的流场特性．结果表明．河口弯道流场具有跟无潮河流弯道显著不同的特性．如弯道流速向下游递增；涨、落急主流线分歧；涨、落潮流对丁坝工程的响应不同．以及下游弯道水面横比降大于上游弯道．     

各向同性湍流中颗粒引起湍流变动的直接模拟

本文通过直接数值模拟对均匀各向同性湍流中颗粒对湍流的变动作用进行了研究.颗粒相的体积分数很小而质量载荷足够大,以至于颗粒之间的相互作用可以忽略不计,而重点考虑颗粒与湍流间能量的交换。颗粒对湍流的反向作用使得湍动能的耗散率增强,以至于湍动能的衰减速率增大.湍动能的衰减速率随颗粒惯性的增大而增大。三维湍动能谱显示,颗粒对湍动能的影响在不同的尺度上是不均匀的。在低波数段,流体带动颗粒,而高波数段则相反.     

外加磁场微波等离子推力器内流场数值模拟

为了提高微波等离子推力器性能，改善等离子体对电磁波能量的吸收状况，提高核心区温度，提出外加磁场的方案，并对热等离子体进行了数值模拟.假设局域热平衡条件，采用Navier-Stokes，Maxwell和Saha方程，利用压力修正的半隐格式和时域有限差分求解方法，建立了径向磁镜场下推力器内等离子体流场的数值计算模型.数值模拟结果表明：外加磁场后的磁感应强度小于0.5 T时，推力器内热等离子体核心区最高温度随磁感应强度的增加而迅速提高.外加磁场后的磁感应强度大于0.5 T时，核心区最高温度随磁感应强度的增加而缓慢提高.磁感应强度为0.5 T时，热等离子体核心区最高温度与不加磁场相比提高了24%.外加磁场对等离子体流场速度分布影响不大. 关键词： 等离子体模拟 等离子体相互作用 等离子体流动