A new numerical scheme of 3rd order Weighted Essentially Non-Oscillatory (WENO) type for 2.5D mixed GLM-MHD in Cartesian coordinates is proposed. The MHD equations are modified by combining the arguments as by Dellar ...A new numerical scheme of 3rd order Weighted Essentially Non-Oscillatory (WENO) type for 2.5D mixed GLM-MHD in Cartesian coordinates is proposed. The MHD equations are modified by combining the arguments as by Dellar and Dedner et al to couple the divergence constraint with the evolution equations using a Generalized Lagrange Multiplier (GLM). Moreover, the magnetohydrodynamic part of the GLM-MHD system is still in conservation form. Meanwhile, this method is very easy to add to an existing code since the underlying MHD solver does not have to be modified. To show the validation and capacity of its application to MHD problem modelling, interaction between a magnetosonic shock and a denser cloud and magnetic reconnection problems are used to verify this new MHD code. The numerical tests for 2D Orszag and Tang's MHD vortex, interaction between a magnetosonic shock and a denser cloud and magnetic reconnection problems show that the third order WENO MHD solvers are robust and yield reliable results by the new mixed GLM or the mixed EGLM correction here even if it can not be shown that how the divergence errors are transported as well as damped as done for one dimensional ideal MHD by Dedner et al.展开更多
帕克太阳探针(Parker Solar Probe,PSP)在太阳附近发现大量磁力线回弯结构,通常还伴随有太阳风速度增加.这些磁力线回弯的产生机制到目前为止有多种解释,其中有代表性的一种是由慢太阳风中的喷流引起的.我们首先对PSP的就地观测数据进...帕克太阳探针(Parker Solar Probe,PSP)在太阳附近发现大量磁力线回弯结构,通常还伴随有太阳风速度增加.这些磁力线回弯的产生机制到目前为止有多种解释,其中有代表性的一种是由慢太阳风中的喷流引起的.我们首先对PSP的就地观测数据进行了统计分析并给出了发生率和空间尺度随径向距离的演化情况,然后使用简化的1.5维磁流体动力学(magnetohydrodynamics,MHD)模型对喷流在太阳风中的演化进行了模拟,其中太阳风被简化为位于黄道面的球对称流.模拟结果表明喷流的确可以导致太阳附近磁力线发生偏转,验证了喷流可以对磁场方向改变有贡献的图景.不过喷流形成的原因还需要进一步研究.展开更多
An analysis was made to study the steady momentum and heat transfer characteristics of a viscous electrically conducting fluid near a stagnation point due to a stretching/shrinking sheet in the presence of a transvers...An analysis was made to study the steady momentum and heat transfer characteristics of a viscous electrically conducting fluid near a stagnation point due to a stretching/shrinking sheet in the presence of a transverse magnetic field and generalized slip condition. Two flow problems corresponding to the planar and axisymmetric stretching/shrinking sheet were considered. By means of similarity transformations, the obtained resultant nonlinear ordinary differential equations were solved numerically using a shooting method for dual solutions of velocity and temperature profiles. Some important physical features of the flow and heat transfer in terms of the fluid velocity, the temperature distribution, the skin friction coefficient and the local Nusselt number for various values of the controlling governing parameters like velocity slip parameter, critical shear rate, magnetic field, ratio of stretching/shrinking rate to external flow rate and Prandtl number were analyzed and discussed. An increase of the critical shear rate decreases the fluid velocity whereas the local Nusselt number increases. The comparison of the present numerical results with the existing literature in a limiting case is given and found to be in an excellent agreement.展开更多
基金Supported by the National Natural Science Foundation of China (40374056, 40536029, 40574068)the International Collaboration Research Team Program of the Chinese Academy of Sciences
文摘A new numerical scheme of 3rd order Weighted Essentially Non-Oscillatory (WENO) type for 2.5D mixed GLM-MHD in Cartesian coordinates is proposed. The MHD equations are modified by combining the arguments as by Dellar and Dedner et al to couple the divergence constraint with the evolution equations using a Generalized Lagrange Multiplier (GLM). Moreover, the magnetohydrodynamic part of the GLM-MHD system is still in conservation form. Meanwhile, this method is very easy to add to an existing code since the underlying MHD solver does not have to be modified. To show the validation and capacity of its application to MHD problem modelling, interaction between a magnetosonic shock and a denser cloud and magnetic reconnection problems are used to verify this new MHD code. The numerical tests for 2D Orszag and Tang's MHD vortex, interaction between a magnetosonic shock and a denser cloud and magnetic reconnection problems show that the third order WENO MHD solvers are robust and yield reliable results by the new mixed GLM or the mixed EGLM correction here even if it can not be shown that how the divergence errors are transported as well as damped as done for one dimensional ideal MHD by Dedner et al.
文摘帕克太阳探针(Parker Solar Probe,PSP)在太阳附近发现大量磁力线回弯结构,通常还伴随有太阳风速度增加.这些磁力线回弯的产生机制到目前为止有多种解释,其中有代表性的一种是由慢太阳风中的喷流引起的.我们首先对PSP的就地观测数据进行了统计分析并给出了发生率和空间尺度随径向距离的演化情况,然后使用简化的1.5维磁流体动力学(magnetohydrodynamics,MHD)模型对喷流在太阳风中的演化进行了模拟,其中太阳风被简化为位于黄道面的球对称流.模拟结果表明喷流的确可以导致太阳附近磁力线发生偏转,验证了喷流可以对磁场方向改变有贡献的图景.不过喷流形成的原因还需要进一步研究.
文摘An analysis was made to study the steady momentum and heat transfer characteristics of a viscous electrically conducting fluid near a stagnation point due to a stretching/shrinking sheet in the presence of a transverse magnetic field and generalized slip condition. Two flow problems corresponding to the planar and axisymmetric stretching/shrinking sheet were considered. By means of similarity transformations, the obtained resultant nonlinear ordinary differential equations were solved numerically using a shooting method for dual solutions of velocity and temperature profiles. Some important physical features of the flow and heat transfer in terms of the fluid velocity, the temperature distribution, the skin friction coefficient and the local Nusselt number for various values of the controlling governing parameters like velocity slip parameter, critical shear rate, magnetic field, ratio of stretching/shrinking rate to external flow rate and Prandtl number were analyzed and discussed. An increase of the critical shear rate decreases the fluid velocity whereas the local Nusselt number increases. The comparison of the present numerical results with the existing literature in a limiting case is given and found to be in an excellent agreement.
