磁流体动力学在航空工程中的应用与展望

介绍了磁流体动力学在航空工程中的主要应用方式,主要包括:磁流体冲压组合发动机、磁流体涡轮组合发动机、燃烧室后磁流体发电、表面磁流体发电、磁流体加速风洞、磁流体推力矢量、进气道大尺寸磁流体流动控制、边界层分离流动控制、边界层转捩控制、飞行器头部热流控制等;探讨了磁流体技术在应用中存在的关键科学与技术问题,对导电流体的产生、磁流体实验设备与实验技术、多场耦合机理及数值模拟方法等进行了分析;最后对磁流体技术在航空工程上的应用与发展进行了总结与展望.     

A high‐order flux reconstruction adaptive mesh refinement method for magnetohydrodynamics on unstructured grids

We report our recent development of the high‐order flux reconstruction adaptive mesh refinement (AMR) method for magnetohydrodynamics (MHD). The resulted framework features a shock‐capturing duo of AMR and artificial resistivity (AR), which can robustly capture shocks and rotational and contact discontinuities with a fraction of the cell counts that are usually required. In our previous paper, 36 we have presented a shock‐capturing framework on hydrodynamic problems with artificial diffusivity and AMR. Our AMR approach features a tree‐free, direct‐addressing approach in retrieving data across multiple levels of refinement. In this article, we report an extension to MHD systems that retains the flexibility of using unstructured grids. The challenges due to complex shock structures and divergence‐free constraint of magnetic field are more difficult to deal with than those of hydrodynamic systems. The accuracy of our solver hinges on 2 properties to achieve high‐order accuracy on MHD systems: removing the divergence error thoroughly and resolving discontinuities accurately. A hyperbolic divergence cleaning method with multiple subiterations is used for the first task. This method drives away the divergence error and preserves conservative forms of the governing equations. The subiteration can be accelerated by absorbing a pseudo time step into the wave speed coefficient, therefore enjoys a relaxed CFL condition. The AMR method rallies multiple levels of refined cells around various shock discontinuities, and it coordinates with the AR method to obtain sharp shock profiles. The physically consistent AR method localizes discontinuities and damps the spurious oscillation arising in the curl of the magnetic field. The effectiveness of the AMR and AR combination is demonstrated to be much more powerful than simply adding AR on finer and finer mesh, since the AMR steeply reduces the required amount of AR and confines the added artificial diffusivity and resistivity to a narrower and narrower region. We are able to verify the designed high‐order accuracy in space by using smooth flow test problems on unstructured grids. The efficiency and robustness of this framework are fully demonstrated through a number of two‐dimensional nonsmooth ideal MHD tests. We also successfully demonstrate that the AMR method can help significantly save computational cost for the Orszag‐Tang vortex problem.     

Unconditionally stable numerical simulations of a new generalized reduced resistive magnetohydrodynamics model

Reduced‐resistive magnetohydrodynamics (MHD) models are used in understanding different phenomenon in various domains, for example, astrophysics to model magnetotail or for solar arcades [Finn Bozkaya JM, Guzdar PN. Loss of equilibrium and reconnection in tearing of two dimensional equilibrias. Physics of Fluids B 1993; 5 :2870–2876], modeling plasma confinements in reverse field pinch [Strauss HR. The dynamo effect in fusion plasmas. Physics of Fluids 1985; 28 :2786–2792] and tokamaks [Strauss HR. Reduced MHD in nearly potential magnetic fields. Journal of Plasma Physics 1997; 57 (1):83–87; Freidberg J. Plasma Physics and Fusion Energy. Cambridge University Press: Cambridge, 2008]. In this context, recently, a new generalized reduced‐resistive MHD model, which can make use of an arbitrary density profile was proposed [Després B, Sart R. Reduced resistive MHD in Tokamaks with general density. ESAIM: Mathematical Modelling and Numerical Analysis 2012; 46 (5):1081–1106. EDP Sciences, SMAI, 2012 online 2011]. We in this work show that this proposed theoretical model can be realized numerically as well, and that it is very robust if the equation set is written in a very particular form using the properties of FEM. To illustrate these points, we pick the current hole configuration [Fujita T, Oikawa T, Suzuki T, Ide S, Sakamoto Y, Koide Y, Hatae T, Naito O, Isayama A, Hayashi N, Shirai H. Plasma equilibrium and confinement in a tokamak with nearly zero central current density in JT‐60U. Physical Review Letters 2001; 87 (11):245001; Hawkes NC, Stratton BC, Tala T, Challis CD, Conway G, DeAngelis R, Giroud C, Hobirk J, Joffrin E, Lomas P, Lotte P, Mailloux J, Mazon D, Rachlew E, Reyes‐Cortes S, Solano E, Zastrow K‐D. Observation of zero current density in the core of JET discharges with lower hybrid heating and current drive. Physical Review Letters 2001; 87 (11):115001], which was modeled using reduced‐resistive MHD and remodel it using different combinations of current sources and density profiles. Our model can be implemented with reasonable computational resources at the price of solving a well‐posed global linear system and it is unconditionally stable. These features are also demonstrated as a part of our numerical experiments. Copyright © 2013 John Wiley & Sons, Ltd.     

Parallel finite element computation of incompressible magnetohydrodynamics based on three iterations

Based on local algorithms, some parallel finite element(FE) iterative methods for stationary incompressible magnetohydrodynamics(MHD) are presented. These approaches are on account of two-grid skill include two major phases: find the FE solution by solving the nonlinear system on a globally coarse mesh to seize the low frequency component of the solution, and then locally solve linearized residual subproblems by one of three iterations(Stokes-type, Newton, and Oseen-type) on subdomains with fine...     

Effect of magnetic Reynolds number on the two‐dimensional hydromagnetic flow around a cylinder

Numerical experiments have been conducted to study the effect of magnetic Reynolds number on the steady, two‐dimensional, viscous, incompressible and electrically conducting flow around a circular cylinder. Besides usual Reynolds number Re, the flow is governed by the magnetic Reynolds number R_m and Alfvén number β. The flow and magnetic field are uniform and parallel at large distances from the cylinder. The pressure Poisson equation is solved to find the pressure fields in the entire flow region. The effects of the magnetic field and electrical conductivity on the recirculation bubble, drag coefficient, standing vortex and pressure are presented and discussed. For low interaction parameter (N<1), the suppression of the flow‐separation is nearly independent of the conductivity of the fluid, whereas for large interaction parameters, the conductivity of the fluid strongly influences the control of flow‐separation. Copyright © 2008 John Wiley & Sons, Ltd.     

Similarity analysis in magnetohydrodynamics: Hall effects on free convection flow and mass transfer past a semi-infinite vertical flat plate

Heat and mass transfer along a semi-infinite vertical flat plate under the combined buoyancy force effects of thermal and species diffusion is investigated in the presence of a strong non-uniform magnetic field and the Hall currents are taken into account. The induced magnetic field due to the motion of the electrically conducting fluid is negligible. This assumption is valid for a small magnetic Reynolds number. The similarity solutions are obtained using the scale group of transformations. These are the only symmetry transformations admitted by the field equations. The non-linear boundary layer equations with the boundary conditions are transferred to a system of non-linear ordinary differential equations with the appropriate boundary conditions. Furthermore, the similarity equations are solved numerically by using a fourth order Runge-Kutta scheme with the shooting method. Numerical results for the velocity profiles, the temperature profiles and the concentration profiles are presented graphically for various values of the magnetic parameter M in the range of 0-1 with the Hall parameter m taking the values 0.5, 1, 2, and 3.     

A divergence-free semi-implicit finite volume scheme for ideal,viscous, and resistive magnetohydrodynamics

In this paper, we present a novel pressure-based semi-implicit finite volume solver for the equations of compressible ideal, viscous, and resistive magnetohydrodynamics (MHD). The new method is conservative for mass, momentum, and total energy, and in multiple space dimensions, it is constructed in such a way as to respect the divergence-free condition of the magnetic field exactly, also in the presence of resistive effects. This is possible via the use of multidimensional Riemann solvers on an appropriately staggered grid for the time evolution of the magnetic field and a double curl formulation of the resistive terms. The new semi-implicit method for the MHD equations proposed here discretizes the nonlinear convective terms as well as the time evolution of the magnetic field explicitly, whereas all terms related to the pressure in the momentum equation and the total energy equation are discretized implicitly, making again use of a properly staggered grid for pressure and velocity. Inserting the discrete momentum equation into the discrete energy equation then yields a mildly nonlinear symmetric and positive definite algebraic system for the pressure as the only unknown, which can be efficiently solved with the (nested) Newton method of Casulli et al. The pressure system becomes linear when the specific internal energy is a linear function of the pressure. The time step of the scheme is restricted by a CFL condition based only on the fluid velocity and the Alfvén wave speed and is not based on the speed of the magnetosonic waves. Being a semi-implicit pressure-based scheme, our new method is therefore particularly well suited for low Mach number flows and for the incompressible limit of the MHD equations, for which it is well known that explicit density-based Godunov-type finite volume solvers become increasingly inefficient and inaccurate because of the more and more stringent CFL condition and the wrong scaling of the numerical viscosity in the incompressible limit. We show a relevant MHD test problem in the low Mach number regime where the new semi-implicit algorithm is a factor of 50 faster than a traditional explicit finite volume method, which is a very significant gain in terms of computational efficiency. However, our numerical results confirm that our new method performs well also for classical MHD test cases with strong shocks. In this sense, our new scheme is a true all Mach number flow solver.     

Radiation and magnetic field effects on the unsteady mixed convection flow of a second grade fluid over a vertical stretching sheet

This present analysis discusses the mixed convection boundary layer flow of a magnetohydrodynamic second grade fluid over an unsteady permeable stretching sheet. The time‐dependent stretching velocity and the surface temperature are chosen. Series solutions of the governing boundary value problems are obtained by employing homotopy analysis method. Convergence of the obtained solution is explicitly discussed. The dependence of velocity and temperature profiles on the various quantities is shown and discussed by plotting graphs. Skin friction coefficient and the local Nusselt number tabulated and analyzed. Copyright © 2010 John Wiley & Sons, Ltd.     

High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

Effects of induced magnetic field on peristaltic flow of Johnson-Segalman fluid in a vertical symmetric channel

In this paper, the influence of heat transfer and induced magnetic field on peristaltic flow of a Johnson-Segalman fluid is studied. The purpose of the present investigation is to study the effects of induced magnetic field on the peristaltic flow of non-Newtonian fluid. The two-dimensional equations of a Johnson-Segalman fluid are simplified by assuming a long wavelength and a low Reynolds number. The obtained equations are solved for the stream function, magnetic force function, and axial pressure gradient by using a regular perturbation method. The expressions for the pressure rise, temperature, induced magnetic field, pressure gradient, and stream function are sketched and interpreted for various embedded parameters.     

Simultaneous impacts of MHD and variable wall temperature on transient mixed Casson nanofluid flow in the stagnation point of rotating sphere

A numerical analysis is provided to scrutinize time-dependent magnetohydrodynamics(MHD) free and forced convection of an electrically conducting non-Newtonian Casson nanofluid flow in the forward stagnation point region of an impulsively rotating sphere with variable wall temperature. A single-phase flow of nanofluid model is reflected with a number of experimental formulae for both effective viscosity and thermal conductivity of nanofluid. Exceedingly nonlinear governing partial differential equations(PDEs)subject to their compatible boundary conditions are mutated into a system of nonlinear ordinary differential equations(ODEs). The derived nonlinear system is solved numerically with implementation of an implicit finite difference procedure merging with a technique of quasi-linearization. The controlled parameter impacts are clarified by a parametric study of the entire flow regime. It is depicted that from all the exhibited nanoparticles,Cu possesses the best convection. The surface heat transfer and surface shear stresses in the x-and z-directions are boosted with maximizing the values of nanoparticle solid volume fraction ? and rotation λ. Besides, as both the surface temperature exponent n and the Casson parameter γ upgrade, an enhancement of the Nusselt number is given.     

Resonance behaviour of viscoelastic fluid in Poiseuille flow in the presence of a transversal magnetic field

The oscillatory flow of a viscoelastic fluid in a circular pipe under the influence of a transversal magnetic field is studied. Exact solutions for the axial velocity and flow rate are presented. The velocity enhancement and the resonance behaviour are analysed both numerically and asymptotically in the case of small pipe radii. Approximations for the resonance frequencies and the achievable velocity enhancements are derived. Copyright © 2005 John Wiley & Sons, Ltd.     

The Calculation of Semi-circular Corrugated Tube——An Application of the General Solution of Ring Shell

In this paper, the deformation and stress distribution in semi-circular corrugated tube under axial force are calculated by means of the general solutions of circular ring shell given in previous paper[1].     

复杂激励下平置板式浮筏功率流传递特性研究

针对隔振理论研究领域中舰船的振动和噪声控制问题，建立六维^[1、2]平置板式弹性浮筏的多机组隔振理论分析计算模型，导出了弹性浮筏传递功率流的表达式，从振动能量传输的角度来评价浮筏系统隔振效果，绘制了功率谱曲线，揭示机组的布置，筏架刚度等结构参数对功率流传递的影响，给出了浮筏设计中结构参数选择的一般准则。     

参数灵敏度分析的神经网络方法及其工程应用

在系统分析中,参数灵敏度分析不仅为判断各系统参数的重要性大小提供了依据,量化的灵敏度指标也是后续参数估计的前提。然而,在多效实际系统中,系统参数与系统状态间的显式函数关系不易得到,导致一阶灵敏度指标无法直接求取。简化的单因素分析方法亦存在模型粗糙、精度不高的缺点。本文研究采用人工神经网络的高精度泛化映射,通过少量样本的训练,建立复杂系统中多个系统参数与系统状态间的近似映射关系,继而推导得到统一的灵敏度计算列式。简单结构的神经网络方法和解析方法的对比计算显示了方法的有效性和可靠性。最后,应用该法对某斜拉桥结构的荷载参数和刚度参数进行了考查,得到一般性结论。     

等离子体诱导圆锥前体流场的PIV测量

在封闭光学玻璃箱体内,应用介质阻挡放电等离子体对20°顶角圆锥附近的静止大气进行了定常和脉冲循环控制,对等离子体诱导的圆锥截面绕流速度场进行了二维PIV测量,在定常控制和脉冲循环控制下,比较分析了最大绕流速度及最大轴向涡量,研究了脉冲频率和电压对上述参数的影响。实验结果表明:存在一个交流电压阈值(Vp-p≈14.0 kV),当电压高于这个阈值后,相对于定常控制模式,脉冲循环控制下沿90°方位角径向线上分布的时间平均切向速度和轴向涡量迅速增加;在脉冲循环控制控制下,动量传递主要表现为分离涡的形成而不是气流的加速。     

A plasma-fluid model for EHD flow in DBD actuators and experimental validation

A numerical method to simulate plasma induced electrohydrodynamic flow is proposed in this study. The numerical model consists of three components. Firstly, a potential module to simulate temporal potential and electric field generated in the ionized fluid. Secondly, a plasma module to simulate plasma development and charge particle densities. Finally, a fluid module to simulate the flow affected by the body forces induced by the movement of the charged particles. Fluid flow is modeled using modified predictor-corrector strategy as proposed in the marker and cell method. The velocity field was corrected to achieve incompressible flow by calculating pressure correction factors, considered in all cells. Numerical convergence and time sensitivity analysis were carried to confirm grid independence and determine an efficient time step for simulations. Numerical computations are validated by comparing with experimental results of discharge currents and current densities. They were found to be in very good agreement thus providing an extensive validation. Furthermore, quiescent flow over a dielectric barrier discharge actuator is simulated in the this study, using the proposed plasma-fluid model, to model flow evolution and resolve temporal flow features for detailed analysis. The streamline and vorticity plots were analyzed and compared with experimental results, and flow results were found to be in-line with the experiments.     

Transport and deposition of neutral particles in magnetohydrodynamic turbulent channel flows at low magnetic Reynolds numbers

The effect of Lorentz force on particle transport and deposition is studied by using direct numerical simulation of turbulent channel flow of electrically conducting fluids combined with discrete particle simulation of the trajectories of uncharged, spherical particles. The magnetohydrodynamic equations for fluid flows at low magnetic Reynolds numbers are adopted. The particle motion is determined by the drag, added mass, and pressure gradient forces. Results are obtained for flows with particle ensembles of various densities and diameters in the presence of streamwise, wall-normal or spanwise magnetic fields. It is found that the particle dispersion in the wall-normal and spanwise directions is decreased due to the changes of the underlying fluid turbulence by the Lorentz force, while it is increased in the streamwise direction. The particle accumulation in the near-wall region is diminished in the magnetohydrodynamic flows. In addition, the tendency of small inertia particles to concentrate preferentially in the low-speed streaks near the walls is strengthened with increasing Hartmann number. The particle transport by turbophoretic drift and turbulent diffusion is damped by the magnetic field and, consequently, particle deposition is reduced.     

Role of inter-particle collision on solids acceleration in riser

Collision among particles plays a significant role in governing the structure of gas-solids flow in a riser,especially in the dense and acceleration region.The inter-particle collision is the major cause not only for the kinetic energy dissipation(in terms of additional pressure drop beyond the solids hold-up) but also for the control of solids acceleration(in terms of a balancing force to prevent a free acceleration of solids).A neglect of the balancing force of inter-particle collision against the hydrody...     

Stability of conducting liquid flowing down an inclined plane at moderate Reynolds number in the presence of constant electromagnetic field

The stability of a conducting viscous film flowing down an inclined plane at moderate Reynolds number in the presence of electromagnetic field is investigated under induction-free approximation. Using momentum integral method a non-linear evolution equation for the development of the free surface is derived. The linear stability analysis of the evolution equation shows that the magnetic field stabilizes the flow whereas the electric field stabilizes or destabilizes the flow depending on its orientation with the flow. The weakly non-linear study reveals that both the supercritical stability and subcritical instability are possible for this type of thin film flow. The influence of magnetic field on the different zones is very significant, while the impact of electric field is very feeble in comparison.