A new operational vector approach for time-fractional subdiffusion equations of distributed order based on hybrid functions

This research study deals with the numerical solutions of linear and nonlinear time-fractional subdiffusion equations of distributed order. The main aim of our approach is based on the hybrid of block-pulse functions and shifted Legendre polynomials. We produce a novel and exact operational vector for the fractional Riemann–Liouville integral and use it via the Gauss–Legendre quadrature formula and collocation method. Consequently, we reduce the proposed equations to systems of equations. The convergence and error bounds for the new method are investigated. Six problems are tested to confirm the accuracy of the proposed approach. Comparisons between the obtained numerical results and other existing methods are provided. Numerical experiments illustrate the reliability, applicability, and efficiency of the proposed method.     

Migration of seismic data

Prospecting for oil and gas resources poses the problem of determining the geological structure of the earth's crust from indirect measurements. Seismic migration is an acoustic image reconstruction technique based on the inversion of the scalar wave equation. Extensive computation is necessary before reliable information can be extracted from large sets of recorded data. In this paper a collection of industrial migration techniques, each giving rise to a data parallel algorithm, is outlined. Computer simulations on synthetic seismic data illustrate the problem and the approach.     

基于Sagnac/双Mach-Zehnder原理的分布式光纤传感系统的偏振态分析

运用琼斯矩阵分析法,建立了该分布式传感器系统在传感光纤双折射和外部事件共同作用下,传输光的偏振态对系统功率传输系数影响的数学模型.仿真分析了在传感光纤线性双折射和圆双折射共同作用下,传输光的偏振态对系统功率传输系数的影响,并提出了利用偏振控制器对输入偏振态进行控制,从而改善干涉信号输出的方法.     

论浮力的计算

通过讨论物体下表面不与液体接触时浮力的计算,总结浮力计算的一般方法.     

基于联合变换相关器的Dempster-Shafer理论计算

提出了用联合变换相关器(JTC)来实现证据理论概率分配函数正交和的光学计算,在理论上对于如何用联合变像相关器实现证据理论分配函数正交和作了详细的讨论,并作了相应的仿真测试,结果显示该方法是可行的。与John Caulfield提出的用声光器件矢量外积实现的正交和计算相比,在该结构中由于采用二进制编码的数值计算,因此其计算精度得到了提高,同时对联合变换相关器输入端二进制编码的数值空间位置的适当调整可以直接得到所需要的证据理论正交和矢量,在处理步骤上得到了简化。     

Accelerating image registration of MRI by GPU-based parallel computation

Automatic image registration for MRI applications generally requires many iteration loops and is, therefore, a time-consuming task. This drawback prolongs data analysis and delays the workflow of clinical routines. Recent advances in the massively parallel computation of graphic processing units (GPUs) may be a solution to this problem. This study proposes a method to accelerate registration calculations, especially for the popular statistical parametric mapping (SPM) system. This study reimplemented the image registration of SPM system to achieve an approximately 14-fold increase in speed in registering single-modality intrasubject data sets. The proposed program is fully compatible with SPM, allowing the user to simply replace the original image registration library of SPM to gain the benefit of the computation power provided by commodity graphic processors. In conclusion, the GPU computation method is a practical way to accelerate automatic image registration. This technology promises a broader scope of application in the field of image registration.     

Using self-similar properties of turbulent premixed flames to downsize chemical tables in high-performance numerical simulations

Detailed chemical mechanisms have to be incorporated in turbulent combustion modelling to predict flame propagation, ignition, extinction or pollutant formation. Unfortunately, hundreds of species and thousands of elementary reactions are involved in hydrocarbon chemical schemes and cannot be handled in practical simulations, because of the related computational costs and the need to model the complexity of their interaction with turbulent motions. Detailed chemistry may be handled using look-up tables, where chemical parameters such as reaction rates and/or species mass fractions are determined from a reduced set of coordinates, progress variables or mixture fractions, as proposed in ILDM, FPI or FGM methods. Nevertheless, these tables may require large computer memory spaces and non-negligible access times. This issue becomes of crucial importance when running on massively parallel computers: to implement these databases in shared memories would induce a large number of data exchanges, reducing the overall code performance; on the other hand duplicating databases in every local processor memory may become impossible either for large databases or small local memories. This work proposes to take advantage of the self-similar behaviour of turbulent premixed flames to reduce the size of these chemical databases, specifically when running on massively parallel machines, under the FPI (Flame Prolongation of ILDM) framework. Several approaches to reduce the database are investigated and discussed both in terms of memory requirements and access times. A very good compromise is obtained for methane–air turbulent premixed flames, where the size of the database is decreased by a factor of 1000, while the access time is reduced by about 60%.     

LBM伪势MRT三维模型GPU并行计算的性能优化

格子Boltzmann方法伪势模型算法中的格点间计算未完全局部化,因此在并行计算时需要更多次的全局内存读写、使用更多数量的寄存器和线程同步操作,从而导致GPU并行计算效率下降.本文针对伪势模型并行计算的局限性,基于三维十五速格子结构的多松弛时间伪势模型,以气液相分离为算例,通过合并访问的方式提高全局内存的读写效率;并提出一种"定向转移"算法,提高格子边界格点获取邻居格点数据的效率;最后探索不同资源分配中各种因素对计算效率的影响,总结最优资源分配的方法.     

一种复杂多边形最大内圆的快速查找算法

任意多边形的最大内圆的快速查找在多个领域有着重要的应用。首先给出了任意多边形的Voronoi图的生成方法,以及基于Voronoi图实现该多边形的最大内圆的寻找算法,在此基础上分别采用中线简化及并行化处理对算法进行了效率方面的改进,并重点研究了矢量数据并行化过程中的数据分配策略,给出了其最优数据划分策略——算法复杂度均衡法。最后对北美阿拉斯加地区的多个湖泊进行了最大内圆的查找实验,相应的实验结果表明,改进的方法能够很大程度地提高最大内圆的查找效率。     

基于JASMIN的地下水流大规模并行数值模拟

针对具有精细网格剖分、长时间跨度特征的地下水流模拟中计算时间长、存储开销大等瓶颈问题,基于MODFLOW三维非稳定流计算方法,提出基于网格片的核心算法以及基于影像区的通信机制,并在JASMIN框架上研制了大规模地下水流并行数值模拟程序JOGFLOW.通过河南郑州市中牟县雁鸣湖水源地地下水流的模拟,对程序正确性和性能进行了验证;通过建立一个具有精细网格剖分的假想地下水概念模型对可扩展性进行测试.相对于32核的并行程序,在512以及1 024个处理机上的并行效率分别可达77.2%和67.5%.数值模拟结果表明,JOGFLOW具有较好的计算性能与可扩展性,能够有效使用数百上千计算核心,支持千万量级以上网格剖分的地下水流模型的大规模并行计算.