基于GPU的快速有限元法求解密度场

doi:10.13700/j.bh.1001-5965.2020.0346

北京航空航天大学学报 ›› 2021, Vol. 47 ›› Issue (10): 2088-2096.doi: 10.13700/j.bh.1001-5965.2020.0346

基于GPU的快速有限元法求解密度场

李果阳¹, 严华¹, 张征宇², 陈沁梅¹, 祝福顺¹

1. 四川大学电子信息学院, 成都 610065;
2. 中国空气动力研究与发展中心高速空气动力研究所, 绵阳 621000

收稿日期:2020-07-21 发布日期:2021-11-08
通讯作者: 严华 E-mail:yanhua@scu.edu.cn
基金资助:
国家自然科学基金（11872069）

A fast finite element method based on GPU to solve density field

LI Guoyang¹, YAN Hua¹, ZHANG Zhengyu², CHEN Qinmei¹, ZHU Fushun¹

1. College of Electronics Information Engineering, Sichuan University, Chengdu 610065, China;
2. High Speed Aerdynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China

Received:2020-07-21 Published:2021-11-08

摘要/Abstract

摘要： 为了快速计算分析利用视频测量方法测得的高速风洞试验密度场在扰动流场作用下的实验数据，针对密度场的数值求解问题，经过光线偏折理论分析密度场得到的二阶偏微分方程，对其研究实现了CPU串行有限元法求解。在此基础上提出了基于GPU的快速有限元求解密度场的方法，该方法经过对串行有限元法求解过程效率分析后，将耗时的神经网络拟合、总刚度矩阵和总载荷向量的求解进行了基于GPU的并行加速。实验结果表明：在精度满足实际工程要求的前提下，相对于CPU串行求解方法，所提方法可大大提高求解效率，且随着网格剖分成倍加密，其加速比成倍增加。

关键词: 密度投影场, 有限元, 并行加速, 二阶偏微分方程, 刚度矩阵

Abstract: In order to quickly calculate and analyze the experimental data under the action of turbulence field, which are measured by video measurement method in high-speed wind tunnel test density projection field, and aimed at the numerical solution of density projection field, the second-order partial differential equation of density projection field is analyzed by ray deflection theory. And the CPU serial finite element method is realized to solve the problem. On this basis, a GPU-based fast finite element method for solving density field is proposed. After analyzing the efficiency of the serial finite element method solving process, the time-consuming neural network fitting, total stiffness matrix and total load vector are solved and then GPU-based parallel acceleration are carried out. The experimental results show that, under the premise that the accuracy meets the actual engineering requirements, the proposed method can greatly improve the solving efficiency compared with the CPU serial solution method, and the acceleration ratio is multiplied with the grid subdivision and multiplied encryption.

Key words: density projection field, finite element, parallel acceleration, second-order partial differential equation, stiffness matrix

中图分类号:

TP301.6

李果阳, 严华, 张征宇, 陈沁梅, 祝福顺. 基于GPU的快速有限元法求解密度场[J]. 北京航空航天大学学报, 2021, 47(10): 2088-2096.

LI Guoyang, YAN Hua, ZHANG Zhengyu, CHEN Qinmei, ZHU Fushun. A fast finite element method based on GPU to solve density field[J]. Journal of Beijing University of Aeronautics and Astronautics, 2021, 47(10): 2088-2096.

参考文献

[1] 郭隆德, 杨建军, 吴运刚, 等.利用光学层析技术重构超声速绕流流场密度分布[J].实验流体力学, 2009, 23(2):68-72.GUO L D, YANG J J, WU Y G, et al.Reconstruction of density distribution in the supersonic flow field with optical computerized tomography[J].Journal of Experiments in Fluid Mechanics, 2009, 23(2):68-72(in Chinese).
[2] MEIER G E A.Cornputerized background oriented schlieren[J].Experiments in Fluids, 2002, 33:181-187.
[3] 张征宇, 王显圣, 黄叙辉, 等.高速复杂流动结构的视频测量[J].航空学报, 2017, 38(8):120989.ZHANG Z Y, WANG X S, HUANG X H, et al.Videogrammetry measurement for high-speed complex flow structures[J].Acta Aeronautica et Astronautica Sinica, 2017, 38(8):120989(in Chinese).
[4] 张正贺, 黄贞, 陈汝婷, 等.基于投影式背景纹影技术的火焰温度场测量仪设计[J].激光与光电子学进展, 2019, 56(5):052803.ZHANG Z H, HUANG Z, CHEN R T, et al.Design of flame temperature measurement instrument based on projective background-oriented schlieren technique[J].Laser & Optoelectronics Progress, 2019, 56(5):052803(in Chinese).
[5] OTA M, KURIHARA K, ARIMOTO H, et al.Background oriented schlieren (BOS) measurement in supersonic flow with 4K high-speed camera[C]//Proc SPIE 10328, Selected Papers from the 31st International Congress on High-Speed Imaging and Photonics, 2017, 1032:103281.
[6] VERSO L, LIBERZON A.Background oriented schlieren in a density stratified fluid[J].The Review of Scientific Instruments, 2015, 86(10):103705.
[7] 张卫庆, 王成亮, 徐洪, 等.基于有限元法的油浸式变压器绕组温度场数值模拟分析[J].自动化与仪器仪表, 2020(4):56-59.ZHANG W Q, WANG C L, XU H, et al.Numerical simulation analysis of winding temperature field of oil-immersed transformer based on finite element method[J].Automation & Instrumentation, 2020(4):56-59(in Chinese).
[8] 夏一帆.面向航空发动机燃烧室点火问题的数值计算方法研究[D].杭州:浙江大学, 2019.XIA Y F.A study of numerical methods for the ignition process in aeroengine combustors[D].Hangzhou:Zhejiang University, 2019(in Chinese).
[9] FEVGAS A, DALOUKAS K, TSOMPANOPOULOU P, et al.A study of sparse matrix methods on new hardware[J].International Journal of Monitoring and Surveillance Technologies Research, 2015, 3(3):1-19.
[10] BARAN B, KÖHLER M, PRASAD N, et al.Numerical solution of large scale sparse matrix equations in python[J].Proceedings in Applied Mathematics and Mechanics, 2014, 14(1):959-960.
[11] 贾光辉, 于云瑞, 王丹.卷积神经网络求解有限元单元刚度矩阵[J].北京航空航天大学学报, 2020, 46(3):481-487.JIA G H, YU Y R, WANG D.Solving finite element stiffness matrix based on convolutional neural network[J].Journal of Beijing University of Aeronautics and Astronautics, 2020, 46(3):481-487(in Chinese).
[12] 刘尧喜, 唐进元, 周炜, 等.扩展有限元刚度矩阵的CSR存储实现[J].机械强度, 2019, 41(6):1384-1390.LIU Y X, TANG J Y, ZHOU W, et al.CSR storage method of extended finite element stiffness matrix[J].Journal of Mechanical Strength, 2019, 41(6):1384-1390(in Chinese).
[13] 曹建伟, 徐翔, 王友年.基于GPU求解椭圆型偏微分方程的并行算法[J].计算物理, 2015, 32(4):475-481.CAO J W, XU X, WANG Y N.Parallel algorithms for separable elliptic equation based on GPU[J].Chinese Journal of Computational Physics, 2015, 32(4):475-481(in Chinese).
[14] LOH W K, KIM Y K.A GPU-accelerated density-based clustering algorithm[C]//2014 IEEE Fourth International Conference on Big Data and Cloud Computing.Piscataway:IEEE Press, 2014:775-776.
[15] 周树荃.有限元结构分析并行计算[M].北京:科学出版社, 1994.ZHOU S Q.Parallel calculation of finite element structural analysis[M].Beijing:Science Press, 1994(in Chinese).
[16] 张俊, 胥頔, 张龙.基于BOS技术的密度场测量研究[J].实验流体力学, 2015, 29(1):77-82.ZHANG J, XU D, ZHANG L.Research on density measurement based on background oriented schlieren method[J].Journal of Experiments in Fluid Mechanics, 2015, 29(1):77-82(in Chinese).
[17] 张俊, 胥頔, 朱博.背景纹影技术算法研究[C]//中国空气动力学会测控技术专委会第六届四次学术交流会, 2013:68-73.ZHANG J, XU D, ZHU B.Background schlieren technique algorithm research[C]//The Chinese Society of Aerodynamic Control Technology of the 6th Four Academic Conference Proceedings, 2013:68-73(in Chinese).
[18] 陈国荣.有限单元法原理及应用[M].2版.北京:科学出版社, 2016:91-116.CHEN G R.Principle and application of finite element method[M].2nd ed.Beijing:Science Press, 2016:91-116(in Chinese).
[19] 朱道建.利用有限元方法求解流体动力学问题[D].武汉:武汉大学, 2017.ZHU D J.To solve the fluid dynamics problems by finite element method[D].Wuhan:Wuhan University, 2017(in Chinese).
[20] 张健飞, 沈德飞.有限元GPU加速计算的实现方法[J].计算机辅助工程, 2014, 23(2):41-45.ZHANG J F, SHEN D F.Implementation method of GPU-accelerated finite element calculation[J].Computer Aided Engineering, 2014, 23(2):41-45(in Chinese).
[21] KIRK D B, HWU W W.Programming massively parallel processors:A hands-on approach[M].2nd ed.New York:Morgan Kaufmann, 2010:19-130.

基于GPU的快速有限元法求解密度场

A fast finite element method based on GPU to solve density field

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	何存富, 李永坤, 吕炎, 宋国荣. 厚胶层复合材料黏接结构中超声反射/透射特性的有限元仿真[J]. 北京航空航天大学学报, 2021, 47(9): 1721-1728.
[2]	喻健, 何宇廷, 冯宇, 马斌麟, 张腾, 杨飞. 环氧树脂基复合材料加筋板结构吸湿行为研究[J]. 北京航空航天大学学报, 2021, 47(9): 1908-1917.
[3]	盖玉新, 李敏. 深海采矿扬矿管的浮力提升方案分析[J]. 北京航空航天大学学报, 2021, 47(8): 1679-1686.
[4]	姜周, 李琳, 范雨, 王文君. 压电周期板中耦合禁带影响规律分析[J]. 北京航空航天大学学报, 2021, 47(7): 1422-1437.
[5]	冯兴, 姚仰平, 李汝宁, 张献民. 山区机场高填方地基变形分析[J]. 北京航空航天大学学报, 2021, 47(10): 2013-2023.
[6]	张婕, 程小全, 程羽佳, 冯振宇. 湿热环境对CCF800/环氧挖补板拉压性能的影响[J]. 北京航空航天大学学报, 2020, 46(6): 1116-1124.
[7]	王保兴, 孟刚, 林苗, 李巍, 曹毅. 3-PPP型柔性并联微定位平台的设计与分析[J]. 北京航空航天大学学报, 2020, 46(4): 798-807.
[8]	贾光辉, 于云瑞, 王丹. 卷积神经网络求解有限元单元刚度矩阵[J]. 北京航空航天大学学报, 2020, 46(3): 481-487.
[9]	李亮, 吴颂平. 矩阵带宽缩减技术在隐式间断有限元中的应用[J]. 北京航空航天大学学报, 2020, 46(3): 532-540.
[10]	张华, 白俊强, 乔磊, 刘艳. 一种适用于浸入有限元方法的网格自适应方法[J]. 北京航空航天大学学报, 2020, 46(3): 588-597.
[11]	岳喜山, 闫群, 赵伟, 谢宗蕻. 单侧面板裂纹损伤对钛合金蜂窝夹层结构弯曲性能影响[J]. 北京航空航天大学学报, 2020, 46(11): 2018-2025.
[12]	赵瑞杰, 王春洁, 阎肃. 潜望式激光通信粗指向装置的多工况拓扑优化[J]. 北京航空航天大学学报, 2020, 46(11): 2114-2120.
[13]	周伟, 张超, 甘淑媛. 基于XFEM-CZM耦合法的胶接接头裂纹扩展分析及强度预测[J]. 北京航空航天大学学报, 2020, 46(11): 2121-2130.
[14]	董虎, 林苗, 顾苏程, 曹毅, 李巍. 多向气动驱动器软体仿生舌弯曲状态的研究[J]. 北京航空航天大学学报, 2019, 45(9): 1882-1893.
[15]	靳涛, 胡小平, 于保华. 大端接圆柱杆的复合圆锥形变幅杆设计及应用[J]. 北京航空航天大学学报, 2019, 45(8): 1630-1638.