二维不可压磁流体力学方程组与高精度数值模拟.pdf

捅斐 的流动进行了数值模拟． 首先，基于二维不可压流的涡量．流函数方法，结合感应电流密度定义，给出二维 不可压磁流体的涡量．流函数．电流密度形式方程组．涡量方程的对流项采用五阶迎 风紧致差分格式离散，控制方程组中的其它一阶导数项和扩散项采用四阶对称紧致差 分格式离散．流函数和势函数的Poisson方程采用四阶对称紧致差分格式离散，并结合 超松弛迭代方法求解． 学精确解问题，验证了高精度紧致差分格式对不可压磁流问题的适用性和有效性．并 考察外加磁场作用下的平面驱动方腔流问题，验证磁雷诺数和哈特曼数的的物理性质． 在中高雷诺数下，观察外加磁场受到的感应磁场的扰动影响、流场的结构随哈特曼数 增大发生的变化及时间延续磁场的稳定性；不同的磁雷诺数下，讨论了磁粘滞作用及 磁对流、磁扩散效应对磁场结构产生的影响． 在上述研究工作的基础上，总结磁流相互作用机理并提出了下一步研究计划和展 望． 关键词：不可压磁流体力学方程组；高精度紧致差分格式：涡量．流函数一电流密度 方法；导电流体；数值模拟：外加磁场 Abstract the the conductive of fluidsin Magnetohydrodynamics(MHD)isstudy electrically mag— netic is ofthe and fields，the governingequmionscomposed theMaxwell this schemeiSusedtosire． equations．Indissertation．High．Precisioncompact of ulatetheflow conductivefluidsinthe of fields． electrically presencemagnetic accountofthe of dimensional formationtwo Firstly,on stream—vorticity N．S and withthedefinitionofinduced stream function． equationscombining electricity．the formoftwodimensional MHD is vorticity-electricdensity incompressibleequationsgiven． Theconvectivetermsof arediscretizedwiththefjfth．ordercom— vorticity equation upwind FDschemes．thediffusivetermsandothertermsarediscretizedwiththefourth．order pact FD schemesofsecondderivativeandfirstderivative symmetricallycompact