二维三温流体力学数值模拟程序的并行化

基于消息传递,通过设计并行算法和组织网格划分,实现了二维三温流体力学Lagrange数值模拟串行程序(Lared-Ⅰ)的并行化。结合物理问题特性和程序的执行特征,设计了动态负载平衡方法,进一步提高并行计算性能。并给出两个并行计算环境上的数值实验结果。     

偏振相移与相位共轭结合的散射光聚焦技术研究

散射介质对光的随机散射作用是制约其光学聚焦和成像的重要因素，光学相位共轭技术能够通过对散射光场共轭还原实现透过散射介质的光学聚焦和成像，其中散射光场相位的获取是共轭还原的核心。阐述了偏振相移的基本原理，将偏振相移与相位共轭技术相结合，设计了新的基于偏振相移的数字光学相位共轭系统。采用633 nm的HeNe激光器作为系统光源，毛玻璃作为散射介质开展散射光聚焦实验研究。实验结果表明:该装置能够成功实现透过散射介质的光学聚焦，其中聚焦点与背景光强的比值可达约400倍。     

基于JASMIN框架的FFT并行解法器及其应用

为解决并行应用程序使用FFTW(Fastest Fourier Transform in the West)并行软件包所面临的计算规模难以扩展、数据结构变动大、实现不同数据结构间通信难度大、接口不确定等问题,在JASMIN框架内设计实现FFT并行解法器.该解法器封装了数据分布存储、数据通信等并行计算细节,通过重新分布存储数据,调用一维FFT变换实现高维FFT的并行计算,并提供规范接口,支撑用户简便地实现FFT的并行计算.数值测试表明,该解法器具有很好的并行性能.该解法器已应用于激光等离子体成丝不稳定性的数值模拟并行程序,它在2048个处理器上的并行效率可达80%以上.     

3维全电磁粒子软件NEPTUNE中的并行计算方法

 介绍了NEPTUNE软件采用的一些并行计算方法：采用“块-网格片”二层并行区域分解方法,使计算规模能够扩展到上千个处理器核。基于复杂几何特征采用自适应技术并行生成结构网格,在原有规则区域的基础上剔除无效网格,大幅降低了存储量和并行执行时间。在经典的Boris和SOR迭代方法基础上,采用红黑排序和几何约束,提出了非规则区域上的Poisson方程并行求解方法。采用这些方法后,当使用NEPTUNE软件模拟MILO器件时,可在1 024个处理器核上获得51.8%的并行效率。     

闪光X射线照相蒙特卡罗程序的并行化

研究了闪光X射线辐射照像蒙特卡罗程序(FXRMC)在MPI平台下的并行计算实现,给出了实现过程中并行随机数的产生方法。对并行程序的测试结果表明:并行程序与串行程序结果一致,加速比比较理想,呈线性增长,并行效率在16个处理器上可达80%以上。算例的结果说明了并行化可有效地解决程序计算散射技术性能时的耗时问题,从而有效化解FXRMC耗时和大规模计算的困难,提高了FXRMC程序的计算规模和计算速度,达到了研究要求。 (Institute of Fluid Physics, CAEP, P. O. Box 919-105, Mianyang 621900, China)     

基于时域多分辨率的非球形气溶胶散射并行计算模型设计及验证EI北大核心CSCD

为解决串行时域多分辨率(MRTD)散射模型运行时间长和内存消耗大的问题,基于消息传递接口(MPI)技术设计了一种非球形气溶胶散射并行计算模型。介绍了MRTD散射模型的基本框架和2种并行数据通信方案,并基于MPI重复非阻塞通信技术实现了MRTD散射模型的并行化设计;搭建了网络并行计算平台,实现了模型的并行化计算。将MRTD散射模型与Mie散射模型、T矩阵法进行了对比,验证了并行MRTD散射模型的计算准确性。结果表明,MRTD模型可较准确地模拟非球形粒子散射特性,并行计算技术可显著提高计算效率;电磁场分量同时交换的并行设计方案的计算效率略高于仅交换磁场分量的方案;通过增加中央处理器核数,程序的并行加速比随之增大,但单核运行效率却略有降低。随着粒子尺度参数的增大,单核计算效率随之增加,复折射率的改变并不会显著影响并行计算效率。     

基于时域多分辨率的非球形气溶胶散射并行计算模型设计及验证

为解决串行时域多分辨率(MRTD)散射模型运行时间长和内存消耗大的问题,基于消息传递接口(MPI)技术设计了一种非球形气溶胶散射并行计算模型。介绍了MRTD散射模型的基本框架和2种并行数据通信方案,并基于MPI重复非阻塞通信技术实现了MRTD散射模型的并行化设计;搭建了网络并行计算平台,实现了模型的并行化计算。将MRTD散射模型与Mie散射模型、T矩阵法进行了对比,验证了并行MRTD散射模型的计算准确性。结果表明,MRTD模型可较准确地模拟非球形粒子散射特性,并行计算技术可显著提高计算效率;电磁场分量同时交换的并行设计方案的计算效率略高于仅交换磁场分量的方案;通过增加中央处理器核数,程序的并行加速比随之增大,但单核运行效率却略有降低。随着粒子尺度参数的增大,单核计算效率随之增加,复折射率的改变并不会显著影响并行计算效率。     

半导体多物理效应并行计算程序JEMS-CDS-Device设计与实现

针对复杂电磁环境下器件多物理效应机理研究需求,研发了半导体多物理效应并行计算程序JEMS-CDS-Device。介绍了JEMS-CDS-Device的架构设计与实现技术。程序基于非结构网格并行框架JAUMIN实现,采用有限体积法(FVM)离散,使用牛顿法全耦合求解"电-载流子输运-热"问题。程序采用"内核+算法库"形式架构,支持2维和3维非结构网格、千万自由度问题并行求解,支持物理方程、离散算法、材料物理模型等的扩展开发。     

用基于方位分辨的漫后向散射光确定粒子尺寸分布及相对折射率

基于偏振门技术,提出了利用漫后向散射光随线偏振光入射方位的变化来测量双层模型表层的粒子尺寸分布及相对折射率。双层物理模型用来模拟上皮组织,上层分别由聚苯乙烯小球或HeLa细胞悬浮液组成,下层由脂肪乳溶液组成,每层的光学特性与上皮组织光学特性匹配。根据米氏(Mie)散射理论,建立了与散射体形态参量相关的偏振后向漫散射强度随入射偏振方位变化的理论模型,并采用浮点遗传算法反演程序获得表层粒子的尺寸分布(PSD)及相对折射率。将理论值和实验测量值进行比较。研究结果表明通过测量随方位变化的后向漫散射光的偏振差分信号,能够获得类上皮组织模型表层的形态结构信息。     

基于Tahoe框架的某夹具并行计算

在开源软件Tahoe框架基础上,结合有限元前后处理程序MSC.Patran及Tecplot,对某复杂夹具进行建模.通过区域分解、编制接口和采用PHG中提供的PCG(preconditioned conjugate gradient,预处理共轭梯度法)迭代解法成功实现262×10⁴自由度模型的串、并行计算.结果表明,并行计算收敛速度更快,4进程并行计算时间不到串行计算时间的1/4.通过与商用程序MSC.Nastran比较,验证计算结果的正确性.利用大型并行计算机对该模型并行计算性能进行研究,获得最高32进程的并行计算加速比.研究表明,改进后的Tahoe计算框架对于开展大规模自由度下的结构并行计算分析研究是可行的,并且随计算节点增加,并行计算过程基本呈线性加速.     

Simple parametric representations of the polarization optical properties of birefringent biological tissues in reflection polarization spectroscopy

The use of polarization reflection spectroscopy for investigating complex scattering media, including orientationally ordered optically anisotropic elements, requires a special choice of output parameters (represented by spectra) that characterize polarized light scattering by objects under study because the standard characteristics, such as the Mueller matrices, depend on the azimuthal orientation of the medium. We propose compact and convenient sets of experimentally determined output parameters that involve separate detection of co- and cross-polarized components of light scattered by a sample irradiated by a normally incident linearly polarized beam and include characteristics invariant with respect to the azimuthal rotation of the sample.     

Measurements of diffuse backscattering Mueller matrices of turbid media

We report on the development of a method that records spatially dependent intensity patterns of polarized light that arise from illuminating a turbid media with a polarized laser beam and being diffusely backscattered. Our technique employs polarized light from a He-Ne laser () which is focused onto the surface of the scattering medium. A surface area of approximately 2×2 cm centered on the light input point is imaged through polarization analysis optics onto a CCD camera. The Mueller matrix is reconstructed by 49 intensity measurements with various orientations of polarizer and analyzer. The measured Mueller matrix of polystyrene spheres was compared with the theory result of incoherent scattering of light by spheres. The experimental and theory results are in excellent agreement. It appears that the azimuthal patterns of the Mueller matrix are determined by the symmetry of the turbid media.     

浑浊散射介质中被测成分线性吸光度提取与模型转换研究

近红外漫反射光谱法可用于物质成分定性、定量分析,具有快速、无损检测等优点,在食品、制药、环境监测、生命科学等领域有广泛应用。采用近红外漫反射光谱对浑浊散射介质进行成分定量分析时,由于光散射作用,介质吸光度随测量距离不再呈线性变化,不同测量距离下测量模型也很难转换。从扩散方程出发,提出了从漫射光谱中提取线性吸光度的方法,获得了与测量距离无关的介质光学参数--有效衰减系数的光谱,该光谱可反映介质吸收变化,并用于物质成分测量。具体采用了双位置差分法,可在线性吸光度的范围内任意选取两个光源-探测器距离测量,通过差分获得与测量位置无关的有效衰减系数光谱,其中测量距离可依据波段灵活地进行选择。另外,经理论推导可知,差分运算的同时还大大削减了光子扩散对测量的影响,有利于不同散射特性的介质间光谱模型的相互借鉴。经理论计算、蒙特卡洛模拟和实验验证,在1 000~1 360 nm波段,该方法适用于强散射的浑浊水溶液。对散射系数处于较大范围（28.53～87.47 cm-1）的intralipid水溶液进行了测试,获得了线性吸光度的光源-探测器距离范围。进一步以葡萄糖为被测成分为例,获得了葡萄糖在intralipid水溶液中的有效衰减系数光谱。采用测量距离连续可调的漫反射光谱测量系统进行了实验,测量范围位于0～0.3 cm,测试了3种典型的介质（3%,5%,10%Intralipid水溶液）。结果可知,3种介质中葡萄糖光谱呈现一致性,它们之间的测量模型可以方便地进行线性转换。研究对象intralipid水溶液常作为生物组织的仿体,其光学参数也覆盖了常见的牛奶、果汁等被测浑浊溶液,因此测量结果具有广泛的适应性。综上,该研究实现了对浑浊水溶液中物质成分线性吸光度的提取,该方法不仅使得测量位置的选择更加灵活方便,还有利于实现不同散射介质间模型的线性转换,特别适用于介质成分复杂多变的场合下的测量应用,如对人体组织、牛奶或食品中的成分测量等。     

New technique to simulate angular patterns of polarized signals from multiply scattering media in MFOV and CCD lidars

We present a semi-analytical approach to determine angular patterns of a polarized intensity that form lidar signals from multiply scattering media registered by a monostatic multiple-field-of-view receiver or CCD camera for any polarization state of incident light and any analyzer state at a receiver. Multiply scattering media with a highly forward elongated phase functions are considered. The model of lidar signal formation includes the single near-backscattering and small-angle multiple scattering of light. The developed approach allows computation and analysis of polarized images, including the Mueller matrix images, formed by the laser light backscattered by multiply scattering media, for example, from atmospheric clouds, ocean waters, tissue, etc.     

Random walk of polarized light in turbid media

We study the propagation of polarized light in turbid media as a random walk of vector photons. Both propagation and polarization directions of light are found to isotropize, following a power law of the number of scattering events. The characteristic length scale governing light isotropization and linear depolarization, the isotropization length , is derived using the exact Mie scattering for spherical particles. A simple relation is obtained for Rayleigh-Gans scatterers where is the transport mean free path and is the mean cosine of scattering angles.     

A Monte Carlo study of penetration depth and sampling volume of polarized light in turbid media

Detection depth and sampling volume of polarized light in highly turbid, cylindrically-shaped samples are estimated using pathlength distributions calculated from a polarization-sensitive Monte Carlo model. Due to defined ranges of the polarized light pathlength distribution, the estimated penetration depth and the interrogated volume of the polarization-maintaining photon subpopulation are smaller than those of the whole collected photon population, the latter exhibiting a wider pathlength distribution resulting from multiple scattering. It is also demonstrated that the spatial interrogation extent of polarized light in turbid media is greatly affected by the experimental detection geometry.     

Ultimate degree of residual polarization of incoherently backscattered light for multiple scattering of linearly polarized light

By invoking ideas about the distribution of the optical paths of partial components of the scattered field, we obtain an expression for estimating the degree of residual polarization of light that is incoherently backscattered from a disordered multiply scattering semi-infinite medium illuminated by linearly polarized light. In the backscattering regime, the depolarization length of the linearly polarized light in the disordered medium becomes smaller with the passage from the isotropic to anisotropic scattering. Experiments with model media featuring substantially anisotropic scattering (the anisotropy parameter of 0.90 ≤ g ≤ 0.95) demonstrated that for backscattering of linearly polarized light, the depolarization length is close to the transport length of the scattering medium.     

On the derivation of vector radiative transfer equation for polarized radiative transport in graded index media

Light transport in graded index media follows a curved trajectory determined by Fermat's principle. Besides the effect of variation of the refractive index on the transport of radiative intensity, the curved ray trajectory will induce geometrical effects on the transport of polarization ellipse. This paper presents a complete derivation of vector radiative transfer equation for polarized radiation transport in absorption, emission and scattering graded index media. The derivation is based on the analysis of the conserved quantities for polarized light transport along curved trajectory and a novel approach. The obtained transfer equation can be considered as a generalization of the classic vector radiative transfer equation that is only valid for uniform refractive index media. Several variant forms of the transport equation are also presented, which include the form for Stokes parameters defined with a fixed reference and the Eulerian forms in the ray coordinate and in several common orthogonal coordinate systems.     

利用表面散射光偏振差异的目标识别技术

 分析了不同物体散射场偏振特性的差异,依据消光定理的矢量微扰解方法将物体散射场分为零阶和高阶解,零阶反射光可完全保持入射光的偏振度,高阶散射则会导致偏振度的降低,因而总散射光的偏振度依赖于散射表面的粗糙程度。提出了利用斯托克斯-穆勒体系测量物体消偏特性的方法并通过实验对一些物体作了测量,实验结果表明：光滑表面可以较好地保留入射光偏振度,而粗糙表面则有很强的消偏作用,因此偏振成像方法可有效地提高目标探测和识别效率。