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

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

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

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

分析端接频变负载的多导体传输线FDTD新方法

提出了一种新的有效求解端接频变负载的时域有限差分法（FDTD）.该方法将频域导纳经矢量匹配后以有理式表示,通过逆傅里叶变换将其转换到时域,基于分段线性递归卷积法,推导出该导纳两端的电压和电流关系,将其代入到传输线FDTD模型中,得到了求解端接频变负载多导体（MTL）传输线FDTD的一般公式.对一典型电路进行了仿真,和以状态变量法求解任意负载的结果进行了对比,两者基本重合,说明该方法用来求解频变负载是有效的. 关键词： 多导体传输线 分段线性递归卷积 时域有限差分法 任意负载     

基于MPI的分子对接并行算法

基于消息传递接口(Message Passing Interface,MPI),用两种不同的并行程序设计方法对Autodock程序进行修改.将修改后的程序应用于HIV-1蛋白酶(Protease)和小分子抑制剂XK263的对接体系,测试了并行程序的加速比和并行效率.结果表明,两种改进的并行Autodock程序都可以很好地完成计算,尤其是方案Ⅱ并行程序的加速比和并行效率更高.     

强光电离团簇并行数值模拟

 为研究强激光电离氢原子团簇,在理论上采用1维氢原子团簇的经典粒子动力学模型,结合粒子对(PP)算法及粒子模拟（PIC）方法,采用自行搭建的9节点并行集群系统,利用消息传递接口（MPI)与OpenMP混合编程模型进行了并行数值模拟计算,获得了较为理想的计算加速比。并且引入了弛豫时间参数,有效地处理了粒子间的碰撞过程,在极大简化计算量的同时,保留了物理本质。所得模拟结果与已有的实验结果符合较好,表明该并行计算模型是稳定、可行的。     

非定常输运问题适应于消息传递并行编程环境的香农熵计算方法

在多计算步的非定常输运问题的蒙特卡罗模拟中,为自动调整每一步的样本数以获得较高的计算效率,可以有多种准则.一种可选的方法是在每一步每隔若干样本监测一次系统中未死亡粒子属性分布对应的香农熵的收敛情况以决定何时停止追加样本,此种方法需要在每一步频繁计算香农熵值.由于在MPI消息传递并行编程环境下香农熵的经典计算方法必须广播大量的数据,导致每一步的计算时间随香农熵计算频率的提高而快速增大,这显然是不能满足实际需求的.本文提出了一种适应于消息传递并行编程环境的香农熵计算新方法,该方法计算得到的香农熵值并不等价于经典方法,但二者之间的差别会随着样本数的增加而趋于零.新方法的最大优势是高频计算香农熵值的时间代价大为降低,为最终实现基于香农熵收敛判断的每步样本数的自动调整奠定了必要的基础.     

利用节点效率评估复杂网络功能鲁棒性

为了克服现有复杂网络鲁棒性研究模型只考虑节点失效的局部影响性和网络拓扑鲁棒性的缺陷, 提出了一种利用节点效率来评估复杂网络功能鲁棒性的方法. 该方法综合考虑节点失效的全局影响性, 利用网络中节点的效率来定义各节点的负载、极限负载和失效模型, 通过打击后网络中最终失效节点的比例来衡量网络的功能鲁棒性, 并给出了其评估优化算法. 实验分析表明该方法对考虑节点负载的复杂网络功能鲁棒性的评定可行有效, 对于大型复杂网络可以获得理想的计算能力.     

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

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

二维分子动力学程序（MDP）的并行与优化

基于消息传递并行程序设计平台MPI,结合当前微处理器的高性能特征,探讨了二维分子动力学数值模拟程序的优化和并行,且具体应用到高速碰撞物理模型问题。其中,串行优化的性能提高了两倍,并行效率在由8台高性能微机构成的微机机群上。均大于90％。     

二维分子动力学程序(MDP)的并行与优化

基于消息传递并行程序设计平台ＭＰＩ，结合当前微处理器的高性能特征，探讨了二维分子动力学数值模拟程序的优化和并行，且具体应用到高速碰撞物理模型问题。其中，串行优化的性能提高了两倍，并行效率在由８台高性能微机构成的微机机群上，均大于９０％。     

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.