MC程序并行设计及提高加速比措施

MC程序的并行设计涉及算法及模块划分,它直接关系到并行加速效率的高低,中子－γ耦合输运蒙特卡罗程序MCNP经过行改造,实现了PVM和MPI两种系统下的并行化,由于作了模块化设计,并行加速效率极佳,PVM版和MPI版并行程序在多个处理器下的加速比均呈线性增长,相比PVM,MPI的适应性列强,多数情况下其效率高于OPVM,并行MCNP程序的计算结果可靠,MPI并行程序在16、32和64个处理器上的并行效率分别达到99％、97％和89％。     

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

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

三维光滑粒子流体动力学并行计算程序

介绍了基于消息传递并行程序设计平台研制的三维光滑粒子流体动力学并行程序CSPH3D.包括计算格式、并行方案、并行程序逻辑,以及加快邻域粒子搜索的处理方法.对三维微喷射和斜侵彻的计算表明:CSPH3D程序可以较好地计算这类问题.并且程序具有较高的并行效率.对于粒子总数为1527402的微喷射算例和粒子总数为1454225的斜侵彻算例,使用100个处理器时,并行效率可以达到80%.     

三维多物质弹塑性流体动力学Euler方法的并行算法研究及程序测试

 并行计算是解决爆炸与冲击问题大规模数值模拟最有效的手段之一。针对Euler方法并行程序设计的复杂性,阐述了三维多物质弹塑性流体动力学程序MMIC-3D并行设计的总体策略,基于消息传递接口(MPI)设计出相应的PMMIC-3D并行程序,并提出了一套实用的程序测试方案。结合聚能射流形成过程的数值模拟算例,在八节点的集群上测试了加速比、并行效率及可扩放性,分析了影响并行性能的因素。     

二维多群辐射输运程序LARED-R-1的并行化

利用有向图描述数据依赖关系,应用已有的并行流水线通量扫描算法,实现基于非协调网格的二维辐射输运程序LARED-R-1的并行化.同时,采用消息缓冲技术提高并行程序的性能.经测试,对于典型的问题规模(100群、3800个网格单元、40个方向),在某并行机的64个和128个处理器上,并行程序分别获得80%和53%的并行效率.     

并行平台上的CFD第一部分:标准的建立

为解决分布式网络并行技术在计算流体力学应用中的主要问题,提出了"并行平台上的CFD通用界面标准”.应用该CFD通用界面标准简化了传统的CFD程序改为并行程序的工作,使不同的CFD核心程序可以方便而快捷的改为并行程序,应用于不同的算例中;提高了CFD并行程序的移植性,CFD核心层同并行边界层无需作任何改动就可应用到不同的并行平台上.在这个标准基础上开发了基于WINDOWSNT系统的WINSOCK并行平台,搭建了基于LINUX系统的PVM并行平台.     

多核环境下潮汐分析辅助软件并行绘制技术研究

为方便科研人员使用Baytap-G潮汐分析软件,基于C#.NET程序设计语言设计了一套潮汐分析辅助软件,实现了Baytap-G软件的封装,完成了输入数据格式的自动转换、输出数据的自动提取,以及水位、振幅、相位数据的可视化展示;在介绍C#.NET语言中任务(TASK)并行编程模型的基础上,详细讨论了绘制振幅、相位图任务分解、并行化程序设计考虑及其实现;实验结果表明,合理设计并行程序可以充分利用多核计算机的计算资源,提升程序运行效率,但是过多的任务数、不均匀的工作负载通常会影响并行程序的效率。     

结构材料辐照损伤的分子动力学程序GPU并行化及优化

基于NIVIDIA公司的CUDA架构对结构材料辐照损伤的分子动力学程序在单个GPU上进行并行化,并对影响程序运行效率的相关因素进行分析和测试.经过一系列优化,当粒子数为两百万时,对比单CPU的执行时间,优化后的GPU程序其双精度加速比可达112倍,单精度加速比达到了三百倍,为后续扩展多GPU结构材料辐照损伤的分子动力学程序奠定基础.     

2(1/2)维等离子体粒子模拟分布式并行程序设计

设计了两种粒子模拟的并行算法,并对其进行了比较,基于消息传递环境开发2(1/2)维粒子模拟并行程序,测试并分析了并行性能.粒子模拟算法中,给出了一个初始化粒子三维Maxwell速度分布的算法,并对常用的电磁场的Lindman吸收边界作了推广.最后对激光钻孔问题作了粒子模拟计算,验证了该并行程序.     

光互连与电互连的并行运算性能比较

光互连在短距离传输中的应用是未来大规模计算领域的关键技术。针对并行运算系统对高速高效传输的应用需求,研究了不同结构和不同传输方式下并行运算系统的性能,并从并行运算性能这一角度比较了光互连系统与电互连系统的优劣。采用网孔结构和超立方结构两种模型,用加速比和效率两个指标来评价光互连系统和电互连系统在并行运算方面的性能差异。得到光互连系统及电互连系统在进行并行计算时的加速比和效率与计算规模之间的关系,分析了计算规模对计算速比和效率的影响。通过加速比和效率这两个量的比较得到结论,在大规模并行计算中,光互连系统有着电互连系统不可比拟的优势。     

Self-organization of five species in a cyclic competition game

Cyclic competition game models, particularly the “rock–paper–scissors” model, play important roles in exploring the problem of multi-species coexistence in spatially ecological systems. We propose an extended “rock–paper–scissors” game to model cyclic interactions among five species, and find that two of the five can coexistent when biodiversity disappears, which is different from the “rock–paper–scissors” game. As the number of fingers is five, we named the new model the “fingers” game, where the thumb, forefinger, middle finger, ring finger, and little finger cyclically dominate their subsequent species and are dominated by their former species. We investigate the “fingers” model in two ways: direct simulations and nonlinear partial differential equations. An important finding is that the number of species in a cyclic competition game has an influence on the emergence of biodiversity. To be specific, the “rock–paper–scissors” model is in favor of maintaining biodiversity in comparison with the “fingers” model when the variables (population size, reproduction rate, selection rate, and migration rate) are the same. It is also shown that the mobility and reproduction rate can promote or jeopardize biodiversity.     

Contrastive analysis on tone between Putonghua and Taiwan Mandarin

The differences of the pitch and duration of Chinese syllables between Putonghua （PTH） and Taiwan Mandarin （TM） were studied. The speech materials to be used are not only isolated syllables, but also sentences. The results reveal that： For the isolated syllables, T1 and T2 in TM are influenced by Minnan dialect, therefore their pitch are lower than those in PTH. T3 is fall-rise in PTH, while it is fall in TM. Moreover, the syllable duration sequence for different tone is T3〉T2〉T1〉T4 in PTH, while it is T1〉T2〉T3〉T4 in TM. For the syllables in sentences, T2 is mid-rise in PTH, while it is mid-level in TM. And the T3 is longer than T4 but shorter than T1 or T2 in PTH, while it is the shortest in TM. Furthermore the effects of prosodic phrase boundary on duration for different tones are almost the same in PTH, but the lengthening part of T1 or T2 is longer than that of T3 or T4 in TM.     

Experimental study on downstep in Chinese intonation

Downstep in pitch contour of Chinese Putonghua is examined using subtly designed sentences by controlling tone combinations. The results show both automatic and nonautomatic downstep phenomena exist in Chinese. In non-automatic downstep, low tones compress downwards the pitch range of the following syllables. and the main influence of downstep is on topline. Low tone not only lower the topline behind it, but also raise the high tones before it, the effects are compatible with each other. In automatic downstep, the topline of pitch contour in intonational phrase is presented as a linear downtrend, but it differs among speakers due to the effect of personal stress practice. In comparison with downstep phenomenon in other tone or non-tone languages, the downstep ratio in Chinese is not constant, and the domain of downstep is not limited within the adjacent tones.     

Abelian-Higgs phase of SU(2) QCD and glueball energy

It is shown that SU(2) QCD admits an dual Abelian-Higgs phase, with a Higgs vacuum of a type-Ⅱ superconductor. This is done by using a connection decomposition for the gluon field and the random-direction approximation. Using a bag picture with soft wall, we presented a calculational procedure for the glueball energy based on the recent proof for wall-vortices [Nucl. Phys. B 741(2006)1].     

Theoretical analysis and simulation of conjugate heights for dual-conjugate AO system in lidar

A multi-conjugate adaptive optics （MCAO） can offer a possibility of widening field of view （FOV） characterized by the isoplanatic angle, and the choose of conjugate height becomes a basic problem for MCAO, which influences the size of iosplanatic angle. Considering the application of lidar, the isoplanatic angle＇s expressions of two deformable mirrors （DMs） MCAO for uplink and downlink are deduced. The effects of conjugate heights for dual-conjugate AO are thoughtfully discussed, and the isoplanatic angles are further analyzed. The results show that the isopanatic angle varies with the conjugate height and reaches the maximum as the conjugate height is at the optimal altitude. Moreover, the optimal conjugate height changes with the propagation distance.     

Quick single-photon detector with many avalanche photo diodes working on the time division

Due to the limit of response speed of the present single-photon detector, the code rate is still too low to come into practical use for the present quantum key distribution (QKD) system.A new idea is put up to design a quick single-photon detector.This quick single-photon detector is composed of a multi-port optic-fiber splitter and many avalanche photo diodes (APDs).Au of the ports with APDs work on the time division and cooperate with a logic discriminating and deciding unit driven by the clock signal.The operation frequency lies on the number N of ports, and can reach N times of the conventional single-photon detector.The single-photon prompt detection can come true for high repetition-rate pulses.The applying of this detector will largely raise the code rate of the QKD, and boost the commercial use.     

Signal generation in an isotropic medium in scanning electron acoustic microscope

Based on the research in Ref. [5][Materials Science and Engineering, 1989; A122： 57 63], an improved model of heat source is set up, the different modes of Lamb wave in an isotropic sample generated by a chopped electron beam at frequency f are obtained with integral transform and normal function expansion method, and the output signal of PZT coupled at the back surface of the sample is found out. The generation mechanism of SEAM （Scanning Electron Acoustic Microscopy） signal is discussed. It shows that the SEAM is a near field imaging technique with high spatial resolution and its best lateral spatial resolution is about 2√2α （α is the radius of the focused electron beam）. Some of experimental results of SEAM images are presented in the paper and it shows that the spatial resolution of SEAM is better than 0.5 μm and smaller than the thermal diffusion length of the sample. Therefore the character of near field imaging in SEAM is also proved experimentally.     

High efficiency filtering using two-dimensional photonic crystal coupled cavity structures

Transmission spectra of coupled cavity structures (CCSs) in two-dimensional (2D) photonic crystals (PCs) are investigated using a coupled mode theory, and an optical filter based on CCS is proposed. The performance of the filter is investigated using finite-difference time-domain (FDTD) method, and the results show that within a very short coupling distance of about 3λ, where λ is the wavelength of signal in vacuum, the incident signals with different frequencies are separated into different channels with a contrast ratio of 20 dB. The advantages of this kind of filter are small size and easily tunable operation frequencies.