铝丝阵列Z箍缩的辐射磁流体动力学过程

在拉格朗日坐标系下,建立了描述Z箍缩等离子体内爆动力学过程的一维三温辐射磁流体力学方程组,并编制了相应的拉氏程序。利用该程序对美国Sandia实验室Saturn装置上的铝丝阵列内爆实验进行了数值模拟,得到了内爆等离子体各参量的时空分布。其中内爆到心时刻、X光峰值功率、X光总能量等计算结果与实验结果基本一致。表明所建立的物理模型和编制的程序是合理和可靠的。     

一种适用于磁流体力学切向间断的HLLC黎曼算子

磁场的存在使得磁流体力学特征波模不同于流体力学, 因此直接由流体力学HLLC黎曼算子导出的HLLC双中间态在交界的间断面会出现不守恒的问题。通常降级采用HLL单磁场中间态代替HLLC双磁场中间态以实现守恒和计算稳定, 代价是切向间断的模拟精度不足。本文对此进行改进, 在模拟切向间断时仍然保留原有的HLLC双磁场中间态, 同时各守恒量仍然能够满足Toro相容条件; 改进型HLLC算子在间断两侧的磁场分量存在差异, 因此能够更精确还原切向间断面。基于数值测试, 包括一维激波管和切向间断的时变模拟, 以及地球磁层三维数值模拟, 将模拟结果进行对比, 结果表明: 相比于已发展的HLLC算子, 改进型HLLC算子对切向间断具有更好的捕捉精度, 能够达到或接近耗时更多的HLLD算子的模拟精度。     

New Approaches and Solutions of Nonlinear Force-Free Magnetic Field

Some new approaches for nonlinear force-free magnetic field are presented and new exact solutions are found analytically. Examples are given and some implications of results to astrophysical solar plasmas as well as tokamak or/and spheromak plasmas are discussed.PACS numbers: 52.30.Cv, 52.55.-s, 95.30.Qd     

Heat and Mass Transfer in MHD Micropolar Flow Over a Vertical Moving Porous Plate in a Porous Medium

An analysis is presented for the problem of free convection with mass transfer flow for a micropolar fluid via a porous medium bounded by a semi-infinite vertical porous plate in the presence of a transverse magnetic field. The plate moves with constant velocity in the longitudinal direction, and the free stream velocity follows an exponentially small perturbation law. A uniform magnetic field acts perpendicularly to the porous surface in which absorbs the micropolar fluid with a suction velocity varying with time. Numerical results of velocity distribution of micropolar fluids are compared with the corresponding flow problems for a Newtonian fluid. Also, the results of the skin-friction coefficient, the couple stress coefficient, the rate of the heat and mass transfers at the wall are prepared with various values of fluid properties and flow conditions.     

Solution to the MHD flow over a non-linear stretching sheet by homotopy perturbation method

In this study,by means of homotopy perturbation method(HPM) an approximate solution of the magnetohydrodynamic(MHD) boundary layer flow is obtained.The main feature of the HPM is that it deforms a difficult problem into a set of problems which are easier to solve.HPM produces analytical expressions for the solution to nonlinear differential equations.The obtained analytic solution is in the form of an infinite power series.In this work,the analytical solution obtained by using only two terms from HPM soluti...     

Analysis of a kinematic dynamo model with FEM–BEM coupling

We consider a kinematic dynamo model in a bounded interior simply connected region Ω and in an insulating exterior region . In the so‐called direct problem, the magnetic field B and the electric field E are unknown and are driven by a given incompressible flow field w . After eliminating E , a vector and a scalar potential ansatz for B in the interior and exterior domains, respectively, are applied, leading to a coupled interface problem. We apply a finite element approach in the bounded interior domain Ω, whereas a symmetric boundary element approach in the unbounded exterior domain Ω^c is used. We present results on the well‐posedness of the continuous coupled variational formulation, prove the well‐posedness and stability of the semi‐discretized and fully discretized schemes, and provide quasi‐optimal error estimates for the fully discretized scheme. Copyright © 2013 John Wiley & Sons, Ltd.     

MHD free convective flow past semi-infinite vertical permeable wall

In this paper,the basic equations governing the flow and heat transfer of an incompressible viscous and electrically conducting fluid past a semi-infinite vertical permeable plate in the form of partia...     

Local and parallel finite element algorithm for stationary incompressible magnetohydrodynamics

This article presents a local and parallel finite element method for the stationary incompressible magnetohydrodynamics problem. The key idea of this algorithm comes from the two‐grid discretization technique. Specifically, we solve the nonlinear system on a global coarse mesh, and then solve a series of linear problems on several subdomains in parallel. Furthermore, local a priori estimates are obtained on a general shape regular grid. The efficiency of the algorithm is also illustrated by some numerical experiments.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1513–1539, 2017