IMP 重离子治癌中的剂量计算方法

Basic algorithms of biological effective dose optimization and dose distribution on CT image for the heavy ion therapy project at the Institute of Modern Physics(IMP),Chinese Academy of Sciences(CAS) are reported in this paper.Firstly,biological effective dose optimization is conducted in water.According to the relationship between CT number and water equivalent path length,an integral algorithm is used to calculate the average dose within a pixel and then the dose distribution in tissue is derived.Secondly...中文文摘：针对深部肿瘤重离子治疗临床试验的需求,首先在水介质中进行生物有效剂量的优化计算,然后根据CT图像中像素CT值与水等效长度转换系数之间的关系,结合水中的深度剂量分布曲线对每个像素进行积分得到CT图像上的生物有效剂量分布。同时介绍了基于被动式束流配送系统适形照射时的剂量确定方式,并提出二维适形放疗也应使用分层照射方式以适应治疗时的不同要求。这些方法适合目前及今后在IMP进行的重离子治癌临床试验研究中治疗计划系统的需要。

X射线在介质中的剂量分布计算

研究了放射治疗中X射线在介质中的输运过程,编程实现了基于蒙特卡罗方法的剂量计算．并在便于图形处理的软件Matlab中对光子输运结果进行了可视化处理．对X射线在均匀介质和菲均匀介质中的蒙特卡罗模拟结果与实测结果、其他蒙特卡罗软件模拟结果进行了比较,结果符合较好．实验结果表明该方法既可以获得很快的仿真速度,又能得到精确直观的剂量计算结果,为提高放射治疗水平具有重要的指导意义和应用价值．

Rotating wheel filter design and optimization for a uniform depth dose distribution in heavy ion therapy

Treatment planning of heavy-ion radiotherapy involves predictive calculation of not only the physical dose but also the biological dose in a patient body. The goal in designing beam-modulating devices for heavy ion therapy is to achieve uniform biological effects across the spread-out Bragg peak (SOBP). To achieve this, a mathematical model of Bragg peak movement is presented. The parameters of this model have been resolved with Monte Carlo method. And a rotating wheel filter is designed basing on the velocity of the Bragg peak movement.

Numerical approximation improves well survey calculation

It is now possible to improve the precision of well survey calculations by order of magnitude with numerical approximation.

Although the most precise method of simulating and calculating a wellbore trajectory generally requires more calculation than other, less-accurate methods, the wider use of computers in oil fields now eliminates this as an obstacle.

The results of various calculations show that there is a deviation of more than 10 m among the different methods of calculation for a directional well of 3,000 m.¹ Consequently, it is important to improve the precision and reliability of survey calculation-the fundamental, necessary work of quantitatively monitoring and controlling wellbore trajectories.

Hybrid Finite Compact-WENO Schemes for Shock Calculation

Hybrid finite compact (FC)-WENO schemes are proposed for shock calculations. The two sub-schemes (finite compact difference scheme and WENO scheme) are hybridized by means of the similar treatment as in ENO schemes. The hybrid schemes have the advantages of FC and WENO schemes. One is that they possess the merit of the finite compact difference scheme, which requires only bi-diagonal matrix inversion and can apply the known high-resolution flux to obtain high-performance numerical flux function; another is that they have the high-resolution property of WENO scheme for shock capturing. The numerical results show that FC-WENO schemes have better resolution properties than both FC-ENO schemes and WENO schemes. In addition, some comparisons of FC-ENO and artificial compression method (ACM) filter scheme of Yee et al. are also given.

Latin hypercube sampling used in the calculation of the fracture probability

The paper proposes Latin hypercube sampling combined with the stratified sampling of variance reduction technique to calculate accurate fracture probability. In the compound sampling, the number of simulations is relatively small and the calculation error is satisfactory.

An Approximate Method For Calculation Of Fluid Force And Response Of A Circular Cylinder At Lock-In

In this paper, equations calculating lift force of a rigid circular cyclinder at lock-in uniform flow are deduced in detail. Besides, equations calculating the lift force on a long flexible circular cyclinder at lock-in are deduced based on mode analysis of a multi-degree freedom system. The simplified forms of these equations are also given. Furthermore, an approximate method to predict the forces and response of rigid circular cyclinders and long flexible circular cyclinders at lock-in is introduced in the case of low mass-damping ratio. A method to eliminate one deficiency of these equations is introduced. Comparison with experimental results show the effectiveness of this approximate method.

Variation of added mass and its application to the calculation of amplitude response for a circular cylinder

In the present study, analyzed are the variation of added mass for a circular cylinder in the lock-in ( synchronization) range of vortex-induced vibration (VIV) and the relationship between added mass and natural frequency. A theoretical minimum value of the added mass coefficient for a circular cylinder at lock-in is given. Developed are semi-empirical formulas for the added mass of a circular cylinder at lock-in as a function of flow speed and mass ratio. A comparison between experiments and numerical simulations shows that the semi-empirical formulas describing the variation of the added mass for a circular cylinder at lock-in are better than the ideal added mass. In addition, computation models such as the wake oscillator model using the present formulas can predict the amplitude response of a circular cylinder at lock-in more accurately than those using the ideal added mass.

Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.

Numeric Calculation for Transient Thermal Stress Field in Wall-Shape Metal Part During Laser Direct Forming

In order to investigate the transient thermal stress field in wall-shape metal part during laser direct forming, a FEM model basing on ANSYS is established, and its algorithm is also dealt with. Calculation results show that while the wall-shape metal part is being deposited, in X direction, the thermal stress in the top layer of the wall-shape metal part is tensile stress and in the inner of the wall-shape metal part is compressive stress. The reason causing above-mentioned thermal stress status in the wall-shape metal part is illustrated, and the influence of the time and the processing parameters on the thermal stress field in wall-shape metal part is also studied. The calculation results are consistent with experimental results in tendency.

Measurement and Calculation of Vibrational Temperature behind Strong Shock Waves

In the present paper the measured values of vibrational temperature behind strong shock waves are compared with theoretical ones. The histories of vibrational temperature behind strong shock waves in a shock tube were measured using two monochromators. The test gas was pure nitrogen at 100-300Pa, and the speeds of shock waves were 5.0-6.0km/s. The electronic temperature of N-2(+) was also approximately determined from experiment and compared with the experimental vibrational temperature. The results show that the presented calculational method is effective, and the electronic energy of N2+ is excited much faster than its vibrational energy. One Langmuir probe was used to determine the effective time of region 2. The influence of viscosity in the shock tube is also analyzed.