格子玻尔兹曼方法_格子玻尔兹曼方法(LBM)python程序提速

本文介绍了作者在学习格子玻尔兹曼方法(LBM)编程时,如何通过将Fortran代码转换为Python并利用pypy解释器提升Python代码执行速度的过程。尽管pypy对某些库支持有限,但它显著提高了计算效率,甚至比Matlab更快。作者还讨论了遇到的问题及解决方案,并计划进一步优化代码和学习数值传热学。
摘要由CSDN通过智能技术生成

研究生开学已经两周了,一直在学习跟LBM相关的编程知识。由于自己数值传热学的基础不是太好,为了能够快速地融入到现有的工作当中我将工作重心侧重在了编程方面,而不是相关模型和边界条件等的学习。我的主要参考书是:

中文名:格子玻尔兹曼方法--基础与工程应用(附计算机代码)

英文名:Lattice Boltzmann Method -- Fundamentals and Engineering Applications with Computer Codes

这本书的一个好处在于书后有相关案例对应的Fortran代码可以参考,在一些编程思路方面可以节省很多时间不需要自己思考。但与此同时不得不说这本书的代码里还是存在着一些问题。由于笔者能力尚且有限无法完全辨别相关代码的正确与否,所以在此暂不进行有关错误点的整理。


下面是前一个阶段的进度整理以及后一个阶段打算:

  • 熟悉基础理论与简单模型
  • 读懂书后的Fortran代码并修改为Python(扩散、对流-扩散、等温不可压缩、非等温不可压缩)基本上每章都做了1-2个完整的案例,但与书上的结果多多少少存在一些差距,大体上成功,估计细节还要好好推敲比对。
  • 数据可视化处理(仰仗Tecplot,真的很心累,学习资料不是很多)
  • Python代码提速(最开始写完python代码很慢,近乎崩溃,可能是由于python是解释型语言,像C是编译型语言,在计算上速度差距过大。还好现在已经解决,也是这篇文章要详述的地方)
  • 计算过程生成可视化动画(依然是依靠Tecp
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值