超音速等离子体炬的磁流体动力学数值研究

针对设计的喉径2mm、工作电流为100A的拉瓦尔喷嘴,在二维轴对称模型的基础上,对超音速等离子体炬中的流动及其外部射流进行了数值模拟。通过在阳极喷嘴内部采用基于磁矢量势的磁流体动力学模型,避免了对磁感应强度的复杂积分计算,得到了喷嘴内部多场耦合的结果及外部射流的流动状态,分析了喷嘴内部电磁场对等离子体的加速作用及射流发展过程。结果显示,等离子体经历了亚音速→跨音速→超音速的发展过程,最终获得2.3 Ma的超音速射流。研究结果为超音速等离子体炬的工业应用提供了理论基础。     

Rotational flow of Oldroyd-B nanofluid subject to Cattaneo-Christov double diffusion theory

A nanofluid is composed of a base fluid component and nanoparticles, in which the nanoparticles are dispersed in the base fluid. The addition of nanoparticles into a base fluid can remarkably improve the thermal conductivity of the nanofluid, and such an increment of thermal conductivity can play an important role in improving the heat transfer rate of the base fluid. Further, the dynamics of non-Newtonian fluids along with nanoparticles is quite interesting with numerous industrial applications. The present predominately predictive modeling studies the flow of the viscoelastic Oldroyd-B fluid over a rotating disk in the presence of nanoparticles. A progressive amendment in the heat and concentration equations is made by exploiting the Cattaneo-Christov heat and mass flux expressions. The characteristic of the Lorentz force due to the magnetic field applied normal to the disk is studied. The Buongiorno model together with the Cattaneo-Christov theory is implemented in the Oldroyd-B nanofluid flow to investigate the heat and mass transport mechanism. This theory predicts the characteristics of the fluid thermal and solutal relaxation time on the boundary layer flow. The von K′arm′an similarity functions are utilized to convert the partial differential equations(PDEs) into ordinary differential equations(ODEs). A homotopic approach for obtaining the analytical solutions to the governing nonlinear problem is carried out. The graphical results are obtained for the velocity field, temperature, and concentration distributions. Comparisons are made for a limiting case between the numerical and analytical solutions, and the results are found in good agreement. The results reveal that the thermal and solutal relaxation time parameters diminish the temperature and concentration distributions, respectively. The axial flow decreases in the downward direction for higher values of the retardation time parameter. The impact of the thermophoresis parameter boosts the temperature distribution.     

Global Regularity of the Logarithmically Supercritical MHD System in Two-dimensional Space

In this paper, we study the global regularity of logarithmically supercritical MHD equations in $2$ dimensional, in which the dissipation terms are $-\mu\Lambda^{2\alpha}u$ and $-\nu\mathcal{L}^{2\beta} b$. We show that global regular solutions in the cases $0<\alpha<\frac{1}{2},\beta>1,3\alpha+2\beta>3$.     

Decoupled schemes for unsteady MHD equations. I. time discretization

In this article, a time discretization decoupled scheme for two‐dimensional magnetohydrodynamics equations is proposed. The almost unconditional stability and convergence of this scheme are provided. The optimal error estimates for velocity and magnet are provided, and the optimal error estimate for pressure are deduced as well. Finite element spatial discretization and numerical implementation are considered in our article (Zhang and He, Comput Math Appl 69 (2015), 1390–1406). © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 956–973, 2017     

Local regularity criteria of a suitable weak solution to MHD equations

We present a regularity condition of a suitable weak solution to the MHD equations in three dimensional space with slip boundary conditions for a velocity and magnetic vector fields. More precisely, we prove a suitable weak solution are H¨older continuous near boundary provided that the scaled mixed L_(x,t)~(p,q) -norm of the velocity vector field with 3/p + 2/q ≤ 2,2 q ∞ is sufficiently small near the boundary. Also, we will investigate that for this solution u ∈ L_(x,t)~(p,q) with 1≤3/p+2/q≤3/2, 3 p ∞, the Hausdorff dimension of its singular set is no greater than max{p, q}(3/p+2/q-1).     

GLOBAL EXISTENCE OF WEAK SOLUTIONS FOR GENERALIZED QUANTUM MHD EQUATION

We prove the existence of a weak solution for a generalized quantum MHD equation in a 2-dimensional periodic box for large initial data. The existence of a global weak solution is established through a three-level approximation,energy estimates, and weak convergence for the adiabatic exponent γ 1.     

Global solutions to the damped MHD system

The global existence of solutions to a damped MHD system in Besov spaces is obtained. This research is motivated by the open question [Wu et al., J. Nonlinear Sci. 25 (2015), 157–192, Remark 5.1] on the possibility of extension to the existence of global small solutions from a potential Banach space to a Besov space. The present result positively confirms the open question. What is more, the initial data are not required to be uniformly small.     

A Low Order Nonconforming Mixed Finite Element Method for Non-Stationary Incompressible Magnetohydrodynamics System

A low order nonconforming mixed finite element method (FEM) is established for the fully coupled non-stationary incompressible magnetohydrodynamics (MHD) problem in a bounded domain in 3D. The lowest order finite elements on tetrahedra or hexahedra are chosen to approximate the pressure, the velocity field and the magnetic field, in which the hydrodynamic unknowns are approximated by inf-sup stable finite element pairs and the magnetic field by $H^1(\Omega)$-conforming finite elements, respectively. The existence and uniqueness of the approximate solutions are shown. Optimal order error estimates of $L^2(H^1)$-norm for the velocity field, $L^2(L^2)$-norm for the pressure and the broken $L^2(H^1)$-norm for the magnetic field are derived.     

Nonlinear MHD Kelvin-Helmholtz instability in a pipe

Nonlinear MHD Kelvin-Helmholtz (K-H) instability in a pipe is treated with the derivative expansion method in the present paper. The linear stability problem was discussed in the past by Chandrasekhar (1961)^[1] and Xu et al. (1981).^[6]Nagano (1979)^[3] discussed the nonlinear MHD K-H instability with infinite depth. He used the singular perturbation method and extrapolated the obtained second order modifier of amplitude vs. frequency to seek the nonlinear effect on the instability growth rate γ. However, in our view, such an extrapolation is inappropriate. Because when the instability sets in, the growth rates of higher order terms on the right hand side of equations will exceed the corresponding secular producing terms, so the expansion will still become meaningless even if the secular producing terms are eliminated. Mathematically speaking, it's impossible to derive formula (39) when γ ₀ ² is negative in Nagano's paper.^[3]Moreover, even as early as γ ₀ ² → O⁺, the expansion becomes invalid because the 2nd order modifier γ₂ (in his formula (56)) tends to infinity. This weakness is removed in this paper, and the result is extended to the case of a pipe with finite depth. Theproject is supported by the National Natural Science Foundation of China.     

Magnetohydrodynamic mixed convection in a vertical channel

The problem of combined free and forced convective magnetohydrodynamic flow in a vertical channel is analysed by taking into account the effect of viscous and ohmic dissipations. The channel walls are maintained at equal or at different constant temperatures. The velocity field and the temperature field are obtained analytically by perturbation series method and numerically by finite difference technique. The results are presented for various values of the Brinkman number and the ratio of Grashof number to the Reynolds number for both equal and different wall temperatures. Nusselt number at the walls is determined. It is found that the viscous dissipation enhances the flow reversal in the case of downward flow while it counters the flow in the case of upward flow. It is also found that the analytical and numerical solutions agree very well for small values of ε.