面向结构网格自适应并行计算的矩形区域求差集快速算法

结构网格自适应程序需要使用矩形区域求差集算法计算网格层间数据依赖关系和网格层嵌套关系.原有的矩形区域求差集算法时间复杂度较高,成为该类应用大规模并行计算可扩展性能瓶颈.本文利用分而治之的方法,构造近似线性时间复杂度的矩形区域求差集快速算法,并利用区域分解实现该算法的并行计算.分别针对规则矩形区域和多层自适应网格的非规则矩形区域求差集问题,验证该算法的效率.结果表明,该算法具有近似线性计算复杂度,对于大规模计算问题,加速效果显著.     

光码分多址并行图像传输系统的最佳光正交签名图形构造

图像传输和接入要求系统的吞吐量超过Gb/s,若采用常规光码分多址系统来实现图像传输,对光电器件的要求高、实现难度大,解决方案是采用光码分多址并行图像传输系统.针对并行图像传输的签名序列是一种二维光正交签名图形,是一种全新形式的签名码,构造最大容量的光正交签名图形(最佳光正交签名图形)的方法很少这一问题,提出一种区组设计构造算法.算法的思想是根据差集合和光正交码的关系,将差集合和区组设计的概念推广到二维空间,根据二维差集合和光正交签名图形的关系从而构造出最佳光正交签名图形.详细阐述了构造最佳光正交签名图形的步骤,对算法进行了仿真.仿真结果说明该算法可以有效地构造出最佳光正交签名图形.     

基于快速W变换-Ⅱ计算偶离散sine变换-Ⅳ的新算法

偶离散sine变换(EDST)在声信号处理中有广泛的应用.而EDST—IV的快速算法是计算各类EDST的关键.本文导出了用离散W变换-Ⅱ(DWT-Ⅱ)表示EDST—Ⅳ的关系式,由此构造了基于快速W变换-Ⅱ(FWT-Ⅱ)计算EDST-Ⅳ的新算法,并与其它快速算法作了运算量的比较.新算法结构规则所需运算量最少.     

量子计算算法介绍

量子计算机利用量子力学原理进行计算,具有量子并行计算的优势,能够超越经典计算1990年中期,量子算法取得突破,舒尔(Shor)构造了大数质因子的量子算法,葛洛沃(Grover)构造了无序数据库的量子搜索算法,引起了人们对量子计算的重视,极大地推动了量子计算的研究.文章简单介绍了几个典型的量子算法以及量子算法研究的一些新进展.     

基于位相抽取的三维信息加密算法研究

本文提出了以位相抽取为基础的三维信息加密算法.构造由纯振幅和纯位相物体组成的简单三维信息作为加密对象.先用标量衍射理论计算其复振幅并抽取相位分布,再结合双随机位相编码完成加密.解密算法为上述算法的逆过程.计算机模拟结果证实了该算法的有效性、鲁棒性和安全性,并揭示了位相抽取算法用于加密更大信息量三维信息的潜力. 关键词： 傅里叶光学 光学信息安全 三维信息加密 位相抽取     

二维Lagrangian坐标系下可压气动方程组的间断Petrov-Galerkin方法

构造矩形网格下求解Lagrangian坐标系下气动方程组的单元中心型格式. 空间离散采用控制体积间断Petrov-Galerkin方法,时间离散采用二阶TVD Runge-Kutta方法. 利用限制器来抑制非物理震荡并保证RKCV算法的稳定性. 构造的算法可以保证物理量的局部守恒. 与Runge-Kutta间断Galerkin（RKDG）方法相比较,RKCV方法的计算公式少一项积分项使得计算较简单. 给出一些数值算例验证了算法的可靠性及效率.     

基于自适应提升小波变换的图像压缩

介绍一种用基于提升算法的中值滤波构造自适应小波变换方法.其思想是设计自适应的“预测算子”及自适应的“更新”运算以更好的逼近信号,从而达到更“紧凑”的信号表示.实验表明,在图像压缩应用中,在峰值信噪比(PSNR)上,虽然本文算法比不上D9/7小波,但在主观评价上更好.     

多介质流模拟的Runge-Kutta控制体积间断有限元方法（英文）

构造可用于多介质流数值模拟的Runge-Kutta控制体积(RKCV)间断有限元方法.对于多介质流模拟,使用线性和非线性的Riemann问题解法器计算界面处的数值流通量.该方法是一种高精度的数值方法且可以保证流体的局部守恒.数值结果表明,即使是利用线性Riemann问题解法器的计算格式也可获得较好的数值结果.与Runge-kutta间断Galerkin方法的比较展示了本文构造算法的优势.     

基于融合与高斯加权暗通道的单幅图像去雾算法

针对图像去雾算法在景深突变处出现光晕现象和远景区域去雾不足的问题,提出了一种基于融合与高斯加权暗通道的单幅图像去雾算法.利用图像形态学梯度的特点,将形态学梯度图像与暗通道图像线性融合获取融合暗通道,构造自适应的高斯权重参数对融合的暗通道图像逐像素处理获取粗透射率,在使用L1正则化优化透射率,通过大气散射模型与修复的大气光值恢复无雾图像.仿真实验表明,本文算法可以较好地恢复出图像的细节并抑制光晕现象,与几种典型的图像去雾算法客观对比,证实了本文算法的可行性.     

一种基于分级分块处理的背景估计算法

视频序列中连续两帧纯背景图像的帧间差服从高斯分布.利用高斯随机变量的峰态和偏态系数的有关性质构造了一种背景子块与目标子块的分块检测算法,并将多分辨率分析的思想引入其中.对帧差图像进行分级分块处理,构造并完善了分级分块背景重建算法,完成了相应的计算机仿真实验.实验结果表明,该算法能快速、有效的获得背景图像.     

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

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

Lattice Boltzmann Solver for Multiphase Flows: Application to High Weber and Reynolds Numbers

The lattice Boltzmann method, now widely used for a variety of applications, has also been extended to model multiphase flows through different formulations. While already applied to many different configurations in low Weber and Reynolds number regimes, applications to higher Weber/Reynolds numbers or larger density/viscosity ratios are still the topic of active research. In this study, through a combination of a decoupled phase-field formulation—the conservative Allen–Cahn equation—and a cumulant-based collision operator for a low-Mach pressure-based flow solver, we present an algorithm that can be used for higher Reynolds/Weber numbers. The algorithm was validated through a variety of test cases, starting with the Rayleigh–Taylor instability in both 2D and 3D, followed by the impact of a droplet on a liquid sheet. In all simulations, the solver correctly captured the flow dynamics andmatched reference results very well. As the final test case, the solver was used to model droplet splashing on a thin liquid sheet in 3D with a density ratio of 1000 and kinematic viscosity ratio of 15, matching the water/air system at We = 8000 and Re = 1000. Results showed that the solver correctly captured the fingering instabilities at the crown rim and their subsequent breakup, in agreement with experimental and numerical observations reported in the literature.     

Schur-decomposition for 3D matrix equations and its application in solving radiative discrete ordinates equations discretized by Chebyshev collocation spectral method

The Schur-decomposition for three-dimensional matrix equations is developed and used to directly solve the radiative discrete ordinates equations which are discretized by Chebyshev collocation spectral method. Three methods, say, the spectral methods based on 2D and 3D matrix equation solvers individually, and the standard discrete ordinates method, are presented. The numerical results show the good accuracy of spectral method based on direct solvers. The CPU time cost comparisons against the resolutions between these three methods are made using MATLAB and FORTRAN 95 computer languages separately. The results show that the CPU time cost of Chebyshev collocation spectral method with 3D Schur-decomposition solver is the least, and almost only one thirtieth to one fiftieth CPU time is needed when using the spectral method with 3D Schur-decomposition solver compared with the standard discrete ordinates method.     

多级轴流压气机三维气动设计的一种快速方法

本文提出了一种适用于多级压气机的快速三维数值气动设计方法,该方法的核心是应用一种快速网格生成技术和NS方程求解器之间不断的相互迭代,最终得到各叶片排的三维叶片造型。首先给定流量,压比和压气机子午通道的几何形状。在初始设计阶段通过计算得到各叶片排沿叶高周向平均的进出口气流角分布,并把其作为计算的目标参数。然后通过网格生成和NS求解之间的迭代,不断调整目标参数,直到计算收敛。应用本方法设计了一台具有三个重复级的实验用轴流低速压气机。     

一种基于HLLC算法的Mie-Grüneisen多介质混合模型

建立一种修正的HLLC （Harten-Lax-Van Leer-Contact）格式下稳定的Mie-Grüneisen多介质混合计算模型.Mie-Grüneisen混合模型中的通量包括守恒和非守恒两个部分，原始的HLLC格式对守恒部分适用，但是原始的HLLC格式直接用于非守恒部分，很难控制住数值振荡的产生.在原始格式中，间断面的移动速度为间断网格的左侧或右侧速度，修正后替换为网格中的平均速度，经过修正后，对HLLC格式重新进行推导，并随之扩展到二维问题.数值实验表明，利用修改后的HLLC格式，Mie-Grüneisen混合模型可以取得较好的稳定性和准确性.     

3D quantum transport solver based on the perfectly matched layer and spectral element methods for the simulation of semiconductor nanodevices

A 3D quantum transport solver based on the spectral element method (SEM) and perfectly matched layer (PML) is introduced to solve the 3D Schrödinger equation with a tensor effective mass. In this solver, the influence of the environment is replaced with the artificial PML open boundary extended beyond the contact regions of the device. These contact regions are treated as waveguides with known incident waves from waveguide mode solutions. As the transmitted wave function is treated as a total wave, there is no need to decompose it into waveguide modes, thus significantly simplifying the problem in comparison with conventional open boundary conditions. The spectral element method leads to an exponentially improving accuracy with the increase in the polynomial order and sampling points. The PML region can be designed such that less than −100 dB outgoing waves are reflected by this artificial material. The computational efficiency of the SEM solver is demonstrated by comparing the numerical and analytical results from waveguide and plane-wave examples and its utility is illustrated by multiple-terminal devices and semiconductor nanotube devices.     

Experimental and Computational Study of the Shock-Wave Loading of Optically Transparent Objects

Technical Physics - Shock-wave loading of transparent objects is studied with the aid of a theoretical method based on numerical simulation on a 3D regular mesh using an explicit solver in the...     

An efficient mixed-precision,hybrid CPU–GPU implementation of a nonlinearly implicit one-dimensional particle-in-cell algorithm

Recently, an implicit, nonlinearly consistent, energy- and charge-conserving one-dimensional (1D) particle-in-cell method has been proposed for multi-scale, full-f kinetic simulations [G. Chen et al., J. Comput. Phys. 230 (18) (2011)]. The method employs a Jacobian-free Newton–Krylov (JFNK) solver, capable of using very large timesteps without loss of numerical stability or accuracy. A fundamental feature of the method is the segregation of particle-orbit computations from the field solver, while remaining fully self-consistent. This paper describes a very efficient, mixed-precision hybrid CPU–GPU implementation of the 1D implicit PIC algorithm exploiting this feature. The JFNK solver is kept on the CPU in double precision (DP), while the implicit, charge-conserving, and adaptive particle mover is implemented on a GPU (graphics processing unit) using CUDA in single-precision (SP). Performance-oriented optimizations are introduced with the aid of the roofline model. The implicit particle mover algorithm is shown to achieve up to 400 GOp/s on a Nvidia GeForce GTX580. This corresponds to 25% absolute GPU efficiency against the peak theoretical performance, and is about 100 times faster than an equivalent single-core CPU (Intel Xeon X5460) compiler-optimized execution. For the test case chosen, the mixed-precision hybrid CPU–GPU solver is shown to over-perform the DP CPU-only serial version by a factor of ～100, without apparent loss of robustness or accuracy in a challenging long-timescale ion acoustic wave simulation.     

基于黎曼解的移动最小二乘粒子动力学数值方法

叙述移动最小二乘粒子动力学(MLSPH)的基本原理,讨论一维MLSPH计算方法,提出两种不同精度的基于无网格的构造黎曼问题初值的方法.通过黎曼解近似粒子间的相互作用,减小接触间断附近的震荡,并给出多介质一维激波管问题的算例,验证此方法的有效性.     

微射流放大器内部三维流动数值模拟

采用三维流动求解器CFD-ACE ,以水作为工作介质,对不同流动工况的偏向型微射流放大器内部三维流场进行了数值模拟,得到了流场的压力分布图和速度矢量图,并通过计算得到平均流量增益曲线,将偏向型微射流放大器内部流场分析与外部运行特性分析有机地结合起来,为实验方案的优化提供了依据。