BNCT优化网格设计及相关算法研究

用MCNP蒙特卡罗程序模拟了硼中子俘获治疗(BNCT)3种国际基准网格模型, 并与 修正的Snyder椭球模型进行了比较. 在此基础上, 给出了一种保质量守恒、内存量少、易于产生输入文件的4种基本材料成分的BNCT网格模型. 计算结果表明, 在4mm网格下, 新模型可以达到基准模型的精度; 根据解析模型剂量随深度的变化规律, 研究构造了多网格组合模型, 在重要区域计算精度不损失的条件下, 计算时间大大缩短. 最后研究给出了一个既保证精度、又在可接受的时间内完成剂量计算的模型、样本数和相应的算法, 它基本上满足临床BNCT的要求.     

BNCT治疗规划系统MCDB算法及测试

给出BNCT治疗计划系统软件MCDB的算法及测试结果,通过使用材料矩阵描述BNCT的网格模型,配合快速径迹算法和计数矩阵,进行粒子输运模拟及计数,模拟结果显示,MCDB的计算精度与MCNP相当,计算速度较MCNP提高3.1~3.4倍.MCDB模拟1 000万粒子的时间不足2个CPU小时,可保证大部分网格剂量误差在5%以内,基本上达到了BNCT临床要求,可用于30 kW的医院中子照射器的临床治疗.     

BNCT蒙特卡洛剂量计算的混合网格算法研究

在硼中子俘获治疗(BNCT)的蒙特卡洛(MC)剂量计算中,通常使用单一的网格模式,如16mm,8mm,4mm.使用细网格计算资源太大,使用粗网格,计算精度不够,为此,根据粒子穿透深度和计数量的变化梯度,采用混合网格模拟计算,达到了细网格的精度,时间仅为细网格的37%.     

5×5棒束CFD模拟计算域和网格标准

采用不可压N-S方程、标准k-ε湍流模型、SIMPLE算法,对某压水堆的两跨距5×5燃料组件的水循环区域流动进行数值分析.对比流道横截面上的压降、速度及其横向分量的相对误差,得到计算域的进出口段长度、网格形式、网格尺寸的选取标准.结果显示：计算域的进口段、出口段长度不宜过短或过长,分别为2.0倍、2.3~2.8倍格架中心距长度;格架区宜采用非结构蜂窝网格,光棒区宜采用结构网格;格架区网格尺寸可取0.3 mm,最大不宜超过0.35 mm,而结构网格横向尺寸可取0.3 mm~0.5 mm,纵向稀疏比应小于1.5为宜.满足上述标准,数值分析能保证较高的可信度,同时保持可接受的计算量.     

α模型的改进和Levy指数的解析计算

为了用Levy稳定性理论研究高能碰撞多重产生中的间歇行为,除了采用蒙特卡罗模拟外,有必要找到一种能对Levy指数作解析计算的方法.为此,改进了随机级联α模型,使其中基元分割几率ω的表达式是线性的.用这一模型,解析地计算了Levy指数,得到了与理论预言相一致的结果;并推出了ω的任意阶矩的解析表达式.为了将这一模型和原有的α模型进行比较,用两种模型分别做MC模拟,计算了间歇指数,证明两种模型都是实用的、可靠的模型.     

基于前沿推进的改进型PEBI网格生成方法

提出一种基于前沿推进的改进型PEBI网格生成方法,将其应用于油藏渗流模拟中,编制相应的油藏数值模拟器.与商业数模软件CMG的计算结果对比,验证了方法的准确性.以某不规则油藏模型为例开展数值模拟研究.结果表明,与常规笛卡尔网格相比,PEBI网格能减弱网格取向的影响.本文提出的PEBI网格生成方法在复杂油藏边界的精确描述、井位及断层内部约束条件的植入等方面具有较强的可应用性.     

可压缩流计算的九分笛卡尔网格技术

简要介绍了一种新的二维网格技术——九分笛卡尔网格技术.与传统方法将一个网格等分为4个不一样,该方法将一个网格单元等分为9个.同时考虑各向同性和各向异性加密,能够对激波和壁面附近进行快速加密,并且有效节省了网格点数目.用此方法计算了单个翼型、多个翼型及直管道下壁面带凸包的无粘带激波的流动.     

四阶紧致差分格式对静力模式中垂直网格计算频散性的影响

为了说明四阶紧致差分格式在大气和海洋数值模式中的潜在价值,提出一种通用方法,推导静力线性斜压适应方程组在微分和差分情况下的频散关系,水平尺度分100 km,10 km和1 km三种情况,从频率、水平群速和垂直群速方面,对采用二阶中央差和四阶紧致差分格式情况下,非跳点网格(N网格)、Lorenz网格(L网格)、Charney-Phillips网格(CP网格)、Lorenz时间跳点网格(LTS网格)和Charney-Phillips时间跳点网格(CPTS网格)的计算特性进行比较,发现采用高精度的四阶紧致差分格式总体上可以明显减少上述三种水平尺度波动在N网格、CP网格、L网格和CPTS网格上的频率、水平群速和垂直群速误差,但对LTS网格,采用四阶紧致差分格式,会使得计算水平群速和垂直群速误差变大.     

嵌套网格技术中的Collar网格和虚拟网格方法

针对嵌套网格技术中的结合部问题,提出了Collar网格和虚拟网格方法.用双曲型微分方程生成的处于结合部的有两个边界分别处于不同固定曲面上的Collar网格,在保证计算网格的生成方便快捷而且网格质量高的前提下,解决了为物体结合部的内外边界点提供有效插值单元的问题.虚拟网格为紧贴物面的面网格,它的作用是将物面边界条件传递给其它网格的边界面,而其本身不作流场计算.计算实践表明,将Collar网格和虚拟网格结合起来应用在嵌套网格技术中能保证几何外形不发生变化,有效地处理各种复杂组合体外形的结合部问题.     

基于Voronoi网格的扩散方程差分格式

在Voronoi网格上利用一种基于回路积分法的有限体积法构造扩散方程的的差分格式.在这种特殊的网格上离散扩散方程比通常在四边形网格上离散的格式要简单,不会引进角点未知量,提高了对网格边上的流的离散精度,及差分格式整体精度.这种Voronoi网格上的扩散计算也可以与单元中心流体力学计算耦合.数值算例表明这种格式比四边形网格上的格式精度高,且能更好的应对网格扭曲情形.     

Numerical Simulation of Electromagnetic Waves Scattering by Discrete Exterior Calculus

We show how to construct discrete Maxwell equations by discrete exterior calculus. The new scheme has many virtues compared to the traditional Yee＇s scheme： it is a multisymplectic scheme and keeps geometric properties. Moreover, it can be applied on triangular mesh and thus is more adaptive to handle domains with irregular shapes. We have implemented this scheme on a Java platform successfully and our experimental results show that this scheme works well.     

Dynamics of ablation plasma induced by a femtosecond laser pulse in electric fields

Direct observations of ablation plasma dynamics in electric field is presented. A time-resolved spatial profile of the ablation plasma induced by femtosecond laser ablation (fsLA) with high fluence is visualized using a planar-laser-induced fluorescence (P-LIF) method. The external electric field is produced by installing a mesh electrode at 6 mm from a Samarium solid target. The Sm ion plasma created by the fsLA showed collective motion regardless of the external electric field, until they reached close to the electrode. When the accelerating and decelerating field was applied, the ions almost disappeared behind the electrode from the field of view. The observations are understood utilizing a SIMION simulation with a conceivable potential gradient caused by Debye shield effect, which is that the ablation plasma keeps the same potential as the target voltage and follows electric potential gradient near the mesh electrode. It is also revealed that this effect degrades time-of-flight resolution at high fluence irradiation. This work gives a new direction for further developments of a fsLA time-of-flight spectrometer.     

Application of three-dimensional magnetic resonance image registration for monitoring hip joint diseases

The purpose of this study was to estimate the accuracy of a method in which three-dimensional (3D) magnetic resonance (MR) volume registration is used for monitoring hip joint disease. Data were analyzed using a normalized cross-correlation (NCC) algorithm involving a user-selected 3D box including the proximal femur. Most of the femoral head was not included in the 3D box because it can become deformed during the course of disease. The accuracy of registration around the femoral head was evaluated using five phantoms and clinical MR data of 17 patients with hip joint disease. In the phantom experiment, registration accuracy was evaluated using four fiducial markers attached to the femoral head. In the experiment using clinical data, registration accuracy was evaluated using a landmark in the femoral head. The registration accuracy in the phantom and clinical experiment was 0.43+/-0.18 mm (S.D.) and 1.12+/-0.46 mm (S.D.), respectively. The former is a value less than half the minimum dimension of a voxel (1.25 x 1.25 x 1.0 mm). Although the latter is slightly larger than the minimum dimension of a voxel, actual errors would be smaller because of the uncertainty in landmark localization. In conclusion, the present method based on an NCC algorithm can be used to accurately register serial MR images of the femoral heads with an error on the order of a voxel. We believe that this method is sufficiently accurate for monitoring hip joint diseases.     

General probabilistic approach to the filtration process

We show experimentally that clogging is basically a matter of the probability of the presence of particles. We describe this process as a function of the main variables of the process, namely, the ratio of particle to mesh hole size, the solid fraction, and the number of grains arriving at each mesh hole during one test, with the help of a simple model, the predictions of which are in very good agreement with our experimental data.     

The role of voxel aspect ratio in determining apparent vascular phase behavior in susceptibility weighted imaging

Susceptibility weighted imaging (SWI) uses apparent phase contrast to enhance the contrast-to-noise ratio (CNR) in the magnitude image. In theory, the apparent phase will depend on the aspect ratio when both venous blood and parenchyma occupy the same voxel. To demonstrate the maximal expected effect of the external field from a vein, we model the vein as an infinitely long cylinder perpendicular to the main magnetic field. The results show that the apparent phase of a voxel in the image is a function of resolution, vessel size and, to a lesser degree, vessel center within the voxel. The simulations explain why a negative-phase mask has worked in SWI processing of high-resolution images collected in the transverse direction, despite the expected positive-phase behavior for vessels perpendicular to the main field. The predicted phase behavior from the simulations is in good agreement with that observed from human brain datasets.     

A non-equilibrium Ising model of turbulence

We introduce a model of interacting lattices at different resolutions driven by the two-dimensional Ising dynamics with a nearest-neighbor interaction. We study this model both with tools borrowed from equilibrium statistical mechanics as well as non-equilibrium thermodynamics. Our findings show that this model keeps the signature of the equilibrium phase transition. The critical temperature of the equilibrium models corresponds to the state maximizing the entropy and delimits two out-of-equilibrium regimes, one satisfying the Onsager relations for systems close to equilibrium and one resembling convective turbulent states. Since the model preserves the entropy and energy fluxes in the scale space, it seems a good candidate for parametric studies of out-of-equilibrium turbulent systems.     

Calculated signal-to-noise ratio of MRI detected with SQUIDs and Faraday detectors in fields from 10 microT to 1.5 T

We examine the calculated signal-to-noise ratio (SNR) achievable with different MRI detection modalities in precession fields ranging from 10 microT to 1.5 T. In particular, we compare traditional Faraday detectors with both tuned and untuned detectors based on superconducting quantum interference devices (SQUIDs). We derive general expressions for the magnetic field noise due to the samples and the detectors, and then calculate the SNR achievable for a specific geometry with each modality with and without prepolarization. We show that each of the three modalities is superior in one of the three field ranges. SQUID-based detection is superior to conventional Faraday detection for MRI in precession fields below 250 mT for a 65 mm diameter surface coil placed a distance of 25 mm from the voxel of interest embedded in a cylinder of tissue 50 mm tall and of radius 50 mm. This crossover field, however, is sensitive to the geometry.     

Role of apparent diffusion coefficient values for the differentiation of viable and necrotic areas of breast cancer and its potential utility to guide voxel positioning for MRS in the absence of dynamic contrast-enhanced MRI data

We carried out retrospective analysis of apparent diffusion coefficient (ADC) values in 48 infiltrating ductal breast cancer patients who had dynamic contrast-enhanced magnetic resonance imaging (DCEMRI; Group I) and in 53 patients (Group II) for whom DCEMRI data were not available. Twenty-three patients of Group I showed no necrosis (Group Ia), while in 25 patients, both viable (nonnecrotic) and necrotic tumor areas (Group Ib) were observed on DCEMRI. T1-weighted, fat-suppressed and short inversion recovery images were used to identify the viable and necrotic tumor areas in Group II patients, and necrosis was not seen in 11 patients (Group IIa), while 42 (Group IIb) showed both viable and necrotic tumor areas. The ADCs of the necrotic area of Group Ib (1.79±0.30 ×10(-3) mm(2)/s) and Group IIb (1.83±0.40 ×10(-3) mm(2)/s) patients were similar and significantly higher (P<.01) compared to the ADCs of the viable tumor area of Group Ia (0.96±0.21 ×10(-3) mm(2)/s) and Group IIa (0.90±0.17 ×10(-3) mm(2)/s) patients. Proton MR spectroscopy (MRS) data were also available in these patients, and the ADC values were retrospectively determined from the voxel from which MR spectrum was obtained. These values were compared with the ADC obtained for the viable and necrotic areas of the tumor. ADC of the MRS voxel was similar to that obtained for the viable tumor area in patients of both groups. This interesting observation reveals the potential utility of using ADC values to identify viable tumor area for positioning of voxel for MRS in the absence of DCEMRI data.