请教LBM中用C++编扰流程序的一些问题
1.我参照何雅玲老师的《格子Boltzmann方法的理论及应用》一书其后的顶盖驱动流的程序,编制的3d单球颗粒扰流程序,网格150*50*50,Re=10, 边界也用的是非平衡外推法,但程序运行到n=10步时,输出的数据中就会出现一些“-1.#IND”,请问这是什么原因呢?是这种方法不适合进行3d计算么?2.在2d圆柱扰流C++中,如果我的曲边界处理还用非平衡外推法的话,曲边界上的受力如何处理啊?写成C++代码的话,应该是什么格式呢? -nan(ind)?
页:
[1]