爆轰驱动惰性气体磁流体发电试验研究

磁流体发电装置作为一种特殊的高功率脉冲电源,具有效率高、容量大、启动快的优点,制约其发展的关键在于如何获得高电导率的发电工质.爆轰驱动具有远超常规方式的驱动能力,在提供高温、高电导率气体方面独具优势.将爆轰驱动激波管技术应用于磁流体发电,有利于突破磁流体发电技术瓶颈,故据此开展了基于爆轰驱动激波管技术的惰性气体磁流体发电试验研究.爆轰驱动根据激波管点火位置不同分为反向和正向两种运行模式,反向爆轰驱动可提供时间较长、状态稳定的试验气流,而正向爆轰优势在于产生高焓试验气流.试验系统由爆轰驱动激波管、拉瓦尔喷管、发电通道、电磁铁和真空罐等组成,试验中分别以反向爆轰和正向爆轰驱动激波管产生发电工质,利用激波将惰性气体压缩至高温从而发生电离,形成的等离子体经喷管加速后,最终在法拉第直线型发电机内切割磁感线输出电能.磁场强度0.9 T的条件下,反向爆轰在负载3.5Ω时获得了较稳定的1.9 kW输出功率,持续时间1.5 ms;外接35 mΩ负载时,正向爆轰在0.3 ms内短时输出功率高达212 k W,功率密度为0.2 GW/m³.试验成功验证了基于爆轰驱动激波管技术的惰性气体...     

Assessment of Five SGS Models for Low-Rem MHD Turbulence

Assessment of three regularization-based and two eddy-viscosity-based subgrid-scale (SGS) turbulence models for large eddy simulations (LES) are carried out in the context of magnetohydrodynamic (MHD) decaying homogeneous turbulence (DHT) with a Taylor scale Reynolds number (Re_λ) of 120 and a MHD transition-to-turbulence Taylor-Green vortex (TGV) problems with a Reynolds number of 3000, through direct comparisons to direct numerical simulations (DNS). Simulations are conducted using the low-magnetic Reynolds number approximation (Re_m<<1). LES predictions using the regularization-based Leray- α,LANS- α, and Clark- α SGS models, along with the eddy viscosity-based non-dynamic Smagorinsky and the dynamic Smagorinsky models are compared to in-house DNS for DHT and previous results for TGV. With regard to the regularization models, this work represents their first application to MHD turbulence. Analyses of turbulent kinetic energy decay rates, energy spectra, and vorticity fields made between the varying magnetic field cases demonstrated that the regularization models performed poorly compared to the eddy-viscosity models for all MHD cases, but the comparisons improved with increase in magnitude of magnetic field, due to a decrease in the population of SGS eddies within the flow field.     

Spontaneous rotation in the exact solution of magnetohydrodynamic equations for flow between two stationary impermeable disks

Magnetohydrodynamic (MHD) flow of a viscous electrically conducting incompressible fluid between two stationary impermeable disks is considered. A homogeneous electric current density vector normal to the surface is specified on the upper disk, and the lower disk is nonconducting. The exact von Karman solution of the complete system of MHD equations is studied in which the axial velocity and the magnetic field depend only on the axial coordinate. The problem contains two dimensionless parameters: the electric current density on the upper plate Y and the Batchelor number (magnetic Prandtl number). It is assumed that there is no external source that produces an axial magnetic field. The problem is solved for a Batchelor number of 0–2. Fluid flow is caused by the electric current. It is shown that for small values of Y, the fluid velocity vector has only axial and radial components. The velocity of motion increases with increasing Y, and at a critical value of Y, there is a bifurcation of the new steady flow regime with fluid rotation, while the flow without rotation becomes unstable. A feature of the obtained new exact solution is the absence of an axial magnetic field necessary for the occurrence of an azimuthal component of the ponderomotive force, as is the case in the MHD dynamo. A new mechanism for the bifurcation of rotation in MHD flow is found.     

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.     

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall considering the interaction between squirmer and wall are numerically studied with an immersed boundary-lattice Boltzmann method. The power-law index, Reynolds number, initial orientation angle of squirmer, and initial distance of squirmer from the wall are all taken into account to investigate the swimming characteristics for pusher (β?<?0), neutral squirmer (β?=?0), and puller (β?>?0) (three kinds of swimmer types) near the no-slip boundary. Four new kinds of swimming modes are found. Results show that, for the pushers and pullers, the wall displays an increasing attraction with increasing power-law index n, which differs from the neutral squirmer who always departs from the wall after the first collision with the wall. Both the initial orientation angle and initial distance from the wall only affect the moving situations rather than the moving modes of the squirmers. However, the squirmers depart from the wall as the Reynolds number increases and chaotic orbits appear for some squirmers at Re?=?5. Several typical flow fields are analyzed and the power consumption and torque for different kinds of flows are also studied. It is found that, as the absolute value of β increases, the power consumption generally increases in shear-thinning (n?=?0.4), Newtonian (n?=?1), and shear-thickening (n?=?1.6) fluids. Moreover, the pushers (β?<?0) and the pullers (β?>?0) expend almost the same power if the absolute value of β remains the same. In addition, the power consumption of the squirmers is highly dependent on the power-law index n.     

Influence of heating load on heat transfer characteristics in micro-pin-fin arrays

Experimental investigations were carried out to explore the convective heat transfer in micro pin-fins with different aspect ratios, and the influence of heating load on Nusselt numbers in micro pin-fins with liquid water as working fluid were investigated. The mechanism of convective heat transfer in micro pin-fins at different heating load were studied by 3-D numerical investigations, and the relationships of thermal physical properties change, the end wall effect and axial thermal conduction with Nu numbers in micro pin-fins were analysed. It was found that the thickness of boundary layer was decreased as much as 33.3 % attributed to the destructive effect of thermal physical properties change, and convective heat transfer in the micro pin-fin channel was more than 20 % enhanced by the flow disturbance caused by the increase of temperature difference. The discrepancy of Nu in micro pin-fin channel with different aspect ratios reached 34.59 %, and this discrepancy was reduced by the increase of heating load. The maximum value of impact factors of dynamic viscosity and thermal conductivity on the Nu in micro-pin-fins reached 25.02 and 7.68 %, respectively.     

Eigenfunctionals of Homogeneous Order-Preserving Maps with Applications to Sexually Reproducing Populations

Homogeneous bounded maps B on cones \(X_+\) of ordered normed vector spaces X allow the definition of a cone spectral radius which is analogous to the spectral radius of a bounded linear operator. If \(X_+\) is complete and B is also order-preserving, conditions are derived for B to have a homogeneous order-preserving eigenfunctional \(\theta : X_+ \rightarrow { \mathbb {R}}_+\) associated with the cone spectral radius in analogy to one part of the Krein–Rutman theorem. Since homogeneous B arise as first order approximations at 0 of maps that describe the year-to-year development of sexually reproducing populations, these eigenfunctionals are an important ingredient in the persistence theory of structured populations with mating.     

Effects of Fuel Lewis Number on Localised Forced Ignition of Turbulent Mixing Layers

The influences of fuel Lewis number Le _F on localised forced ignition of inhomogeneous mixtures are analysed using three-dimensional compressible Direct Numerical Simulations (DNS) of turbulent mixing layers for Le _F  = 0.8, 1.0 and 1.2 and a range of different root-mean-square turbulent velocity fluctuation u′ values. For all Le _F cases a tribrachial flame has been observed in case of successful ignition. However, the lean premixed branch tends to merge with the diffusion flame on the stoichiometric mixture fraction isosurface at later stages of the flame evolution. It has been observed that the maximum values of temperature and reaction rate increase with decreasing Le _F during the period of external energy addition. Moreover, Le _F is found to have a significant effect on the behaviours of mean temperature and fuel reaction rate magnitude conditional on mixture fraction values. It is also found that reaction rate and mixture fraction gradient magnitude \(\vert \nabla \xi \vert \) are negatively correlated at the most reactive region for all values of Le _F explored. The probability of finding high values of \(\vert \nabla \xi \vert \) increases with increasing Le _F . For a given value of u′, the extent of burning decreases with increasing Le _F . A moderate increase in u′ gives rise to an increase in the extent of burning for Le _F  = 0.8 and 1.0, which starts to decrease with further increases in u′. For Le _F  = 1.2, the extent of burning decreases monotonically with increasing u′. The extent of edge flame propagation on the stoichiometric mixture fraction ξ = ξ _st isosurface is characterised by the probability of finding burned gas on this isosurface, which decreases with increasing u′ and Le _F . It has been found that it is easier to obtain self-sustained combustion following localised forced ignition in case of inhomogeneous mixtures than that in the case of homogeneous mixtures with the same energy input, energy deposition duration when the ignition centre is placed at the stoichiometric mixture. The difficultly to sustain combustion unaided by external energy addition in homogeneous mixture is particularly prevalent in the case of Le _F  = 1.2.     

Effect of Inflow Turbulence on an Airfoil Flow with Laminar Separation Bubble: An LES Study

The present paper is concerned with numerical investigations on the effect of inflow turbulence on the flow around a SD7003 airfoil. At a Reynolds number Re_c =?60,000, an angle of attack α =?4^° and a low or zero turbulence intensity of the oncoming flow, the flow past the airfoil is known to be dominated by early separation, subsequent transition and reattachment leading to a laminar separation bubble with a distinctive pressure plateau. The objective of the study is to investigate the effect of inflow turbulence on the flow behavior. For this purpose, a numerical methodology relying on a wall-resolved large-eddy simulation, a synthetic turbulence inflow generator and a specific source term concept for introducing the turbulence fluctuations within the computational domain is used. The numerical technique applied allows the variation of the free-stream turbulence intensity (TI) in a wide range. In order to analyze the influence of TI on the arising instantaneous and time-averaged flow field past the airfoil, the present study evaluates the range 0% ≤ TI ≤?11.2%, which covers typical values found in atmospheric boundary layers. In accordance with experimental studies it is shown that the laminar separation bubble first shrinks and finally completely vanishes for increasing inflow turbulence. Consequently, the aerodynamic performance in terms of the lift-to-drag ratio increases. Furthermore, the effect of the time and length scales of the isotropic inflow turbulence on the development of the flow field around the airfoil is analyzed and a perceptible influence is found. Within the range of inflow scales studied decreasing scales augment the receptivity of the boundary layer promoting an earlier transition.     

Effects of Lewis Number on Head on Quenching of Turbulent Premixed Flames: A Direct Numerical Simulation Analysis

The head on quenching of statistically planar turbulent premixed flames by an isothermal inert wall has been analysed using three-dimensional Direct Numerical Simulation (DNS) data for different values of global Lewis number Le(0.8, 1.0 and 1.2) and turbulent Reynolds number Re_t. The statistics of head on quenching have been analysed in terms of the wall Peclet number Pe (i.e. distance of the flame from the wall normalised by the Zel’dovich flame thickness) and the normalised wall heat flux Φ. It has been found that the maximum (minimum) value of Φ(Pe) for the turbulent Le=0.8 cases are greater (smaller) than the corresponding laminar value, whereas both Pe and Φ in turbulent cases remain comparable to the corresponding laminar values for Le=1.0 and 1.2. Detailed physical explanations are provided for the observed Le dependences of Pe and Φ. The existing closure of mean reaction rate \(\overline {\dot {\omega }}\) using the scalar dissipation rate (SDR) in the near wall region has been assessed based on a-priori analysis of DNS data and modifications to the existing closures of mean reaction rate and SDR have been suggested to account for the wall effects in such a manner that the modified closures perform well both near to and away from the wall.     

Hall效应对三维磁流体发生器的影响

应用三维非理想低磁雷诺数磁流体五方程模型发展了对带有强制项的Navier-Stokes方程组采用熵条件格式, 对椭圆型电势方程采用SOR进行迭代的数值方法,研究了Hall效应对磁流体旁路超燃冲压发动机中磁流体发生器流动及性能的影响.磁流体发生器采用电子束获得有效可靠的电导率. 计算结果表明,Hall效应可引起流场和电场的扭曲, 从而诱导出不稳定二次流的发展与演变,并破坏Joule热的分布. 对这些磁流体现象作出了较详细的分析.最后计算了磁流体发生器的性能参数, 说明Hall效应将导致磁流体发生器的性能下降.      

Pre-Chamber Ignition Mechanism: Simulations of Transient Autoignition in a Mixing Layer Between Reactants and Partially-Burnt Products

The structure of autoignition in a mixing layer between fully-burnt or partially-burnt combustion products from a methane-air flame at ? = 0.85 and a methane-air mixture of a leaner equivalence ratio has been studied with transient diffusion flamelet calculations. This configuration is relevant to scavenged pre-chamber natural-gas engines, where the turbulent jet ejected from the pre-chamber may be quenched or may be composed of fully-burnt products. The degree of reaction in the jet fluid is described by a progress variable c (c = taking values 0.5, 0.8, and 1.0) and the mixing by a mixture fraction ξ (ξ = 1 in the jet fluid and 0 in the CH₄-air mixture to be ignited). At high scalar dissipation rates, N₀, ignition does not occur and a chemically-frozen steady-state condition emerges at long times. At scalar dissipation rates below a critical value, ignition occurs at a time that increases with N₀. The flame reaches the ξ = 0 boundary at a finite time that decreases with N₀. The results help identify overall timescales of the jet-ignition problem and suggest a methodology by which estimates of ignition times in real engines may be made.     

Effect of tube diameter and capillary number on platelet margination and near-wall dynamics

I investigate the effect of tube diameter D and red blood cell capillary number Ca (i.e. the ratio of viscous to elastic forces) on platelet margination in blood flow at ≈37 % tube haematocrit. The system is modelled as three-dimensional suspension of deformable red blood cells and nearly rigid platelets using a combination of the lattice-Boltzmann, immersed boundary and finite element methods. Results of simulations during the dynamics before the steady state has been reached show that a non-diffusive radial platelet transport facilitates margination. This non-diffusive effect is important near the edge of the cell-free layer, but only for Ca > 0.2, when red blood cells are tank-treading. I also show that platelet trapping in the cell-free layer is reversible for Ca ≤ 0.2. Margination is essentially independent of Ca only for the smallest investigated tube diameter (D = 10 μm). Once platelets have reached the cell-free layer, they tend to slide rather than tumble. The tumbling rate is essentially independent of Ca but increases with D. Strong confinement suppresses tumbling due to the relatively small cell-free layer thickness at ≈ 37 % tube haematocrit.     

Determination of the Effective Probabilities of Catalytic Reactions on the Surfaces of Heat Shield Materials in Dissociated Carbon Dioxide Flows

The effect of heterogeneous catalysis on the heat transfer to cold and heated surfaces in subsonic dissociated carbon dioxide jet flows is studied experimentally, using a 100 kW inductive plasma generator, and simulated numerically. The effective probabilities of the heterogeneous reactions CO + O → CO₂ and O + O → O₂ on molybdenum (T_w=300 K) and quartz (T_w=470–620 K) surfaces, the Buran heat shield tile coating (T_w = 1470—1670 K), and two oxidation-resistant carbon-carbon coating materials (T_w=1420—1840 K) are determined by comparing the experimental and calculated data on the heat fluxes at the stagnation point of models at a pressure of 0.1 atm.     

Growths of mechanical elasticity and electrical conductance of graphene nanoplatelet/poly(lactic acid) composites under strong electric field: correlation with time evolution of higher order structure of graphene nanoplatelets

For a composite of poly(lactic acid) containing graphene nanoplatelets (GNPs) at a low concentration (0.34 vol%), this study examined growths of mechanical and electrical properties under an alternating current (AC) electric field, focusing on field-induced GNP structures governing those properties. The composite was subjected to the AC field (60 Hz) of various intensities E for various times t _E at 190 °C. A fraction of randomly oriented GNPs was aligned by the field and then connected into columns, as suggested from optical microscopy. This structural evolution led to qualitatively similar growths of low-frequency storage modulus and static electrical conductivity. The key quantity for understanding this growth was a time t _E* for occurrence of short circuit that detected formation of GNP columns conductively bridging the electrodes. The growths of both modulus and conductivity for various E were summarized as functions of a reduced variable, t _E/t _E*, confirming the growths commonly reflected the evolution of the GNP columns. However, the modulus grew fast and leveled off by t _E/t _E* ~ 1, whereas the conductivity kept growing gradually even at t _E/t _E* > 1. This difference was discussed in relation to the matrix chains and leftover GNPs out the column.     

Power-law creep and residual stresses in a carbopol gel

We report on the interplay between creep and residual stresses in a carbopol microgel. When a constant shear stress σ is applied below the yield stress σ _y, the strain is shown to increase as a power law of time, γ(t) = γ ₀ + (t/τ) ^α , with an exponent α = 0.39 ± 0.04 that is strongly reminiscent of Andrade creep in hard solids. For applied shear stresses lower than some typical value σ _c ? 0.2σ _y, the microgel experiences a more complex, anomalous creep behavior, characterized by an initial decrease of the strain, that we attribute to the existence of residual stresses of the order of σ _c that persist after a rest time under a zero shear rate following preshear. The influence of gel concentration on creep and residual stresses are investigated as well as possible aging effects. We discuss our results in light of previous works on colloidal glasses and other soft glassy systems.     

A class of exact solutions for <Emphasis Type="Italic">N</Emphasis>-dimensional incompressible magnetohydrodynamic equations

In this paper, a sufficient and necessary condition is presented for existence of a class of exact solutions to N-dimensional incompressible magnetohydrodynamic (MHD) equations. Such solutions can be explicitly expressed by appropriate formulae. Once the required matrices are chosen, solutions to the MHD equations are directly constructed.     

Operator method for solving the plane elasticity problem for a strip with periodic cuts

In the present paper, we use the conformal mapping z/c = ζ?2a sin ζ (a, c?const, ζ = u + iv) of the strip {|v| ≤ v ₀, |u| < ∞} onto the domain D, which is a strip with symmetric periodic cuts. For the domain D, in the orthogonal system of isometric coordinates u, v, we solve the plane elasticity problem. We seek the biharmonic function in the form F = C ψ ₀ + S ψ*₀ + x(C ψ ₁ ? S ψ ₂) + y(C ψ ₂ + S ψ ₁), where C(v) and S(v) are the operator functions described in [1] and ψ ₀(u), …, ψ ₂(u) are the desired functions. The boundary conditions for the function F posed for v = ±v ₀ are equivalent to two operator equations for ψ ₁(u) and ψ ₂(u) and to two ordinary differential equations of first order for ψ ₀(u) and ψ*₀(u) [2]. By finding the functions ψ _j (u) in the form of trigonometric series with indeterminate coefficients and by solving the operator equations, we obtain infinite systems of linear equations for the unknown coefficients. We present an efficient method for solving these systems, which is based on studying stable recursive relations. In the present paper, we give an example of analysis of a specific strip (a = 1/4, v ₀ = 1) loaded on the boundary v = v ₀ by a normal load of intensity p. We find the particular solutions corresponding to the extension of the strip by the longitudinal force X ^∞ and to the transverse and pure bending of the strip due to the transverse force Y ^∞ and the constant moment M ^∞, respectively. We also present the graphs of normal and tangential stresses in the transverse cross-section x = 0 and study the stress concentration effect near the cut bottom.     

Pressure-Gradient Turbulent Boundary Layers Developing Around a Wing Section

A direct numerical simulation database of the flow around a NACA4412 wing section at R e _c = 400,000 and 5^° angle of attack (Hosseini et al. Int. J. Heat Fluid Flow 61, 117–128, 2016), obtained with the spectral-element code Nek5000, is analyzed. The Clauser pressure-gradient parameter β ranges from ? 0 and 85 on the suction side, and from 0 to ? 0.25 on the pressure side of the wing. The maximum R e _?? and R e _τ values are around 2,800 and 373 on the suction side, respectively, whereas on the pressure side these values are 818 and 346. Comparisons between the suction side with zero-pressure-gradient turbulent boundary layer data show larger values of the shape factor and a lower skin friction, both connected with the fact that the adverse pressure gradient present on the suction side of the wing increases the wall-normal convection. The adverse-pressure-gradient boundary layer also exhibits a more prominent wake region, the development of an outer peak in the Reynolds-stress tensor components, and increased production and dissipation across the boundary layer. All these effects are connected with the fact that the large-scale motions of the flow become relatively more intense due to the adverse pressure gradient, as apparent from spanwise premultiplied power-spectral density maps. The emergence of an outer spectral peak is observed at β values of around 4 for λ _z ? 0.65δ ₉₉, closer to the wall than the spectral outer peak observed in zero-pressure-gradient turbulent boundary layers at higher R e _?? . The effect of the slight favorable pressure gradient present on the pressure side of the wing is opposite the one of the adverse pressure gradient, leading to less energetic outer-layer structures.     

A simple criterion for finite time stability with application to impacted buckling of elastic columns

The existing theories of finite-time stability depend on a prescribed bound on initial disturbances and a prescribed threshold for allowable responses. It remains a challenge to identify the critical value of loading parameter for finite time instability observed in experiments without the need of specifying any prescribed threshold for allowable responses. Based on an energy balance analysis of a simple dynamic system, this paper proposes a general criterion for finite time stability which indicates that finite time stability of a linear dynamic system with constant coefficients during a given time interval [0, t _f ] is guaranteed provided the product of its maximum growth rate (determined by the maximum eigen-root p₁ >0) and the duration t _f does not exceed 2, i.e., p₁t _f <2. The proposed criterion (p₁t _f =2) is applied to several problems of impacted buckling of elastic columns: (i) an elastic column impacted by a striking mass, (ii) longitudinal impact of an elastic column on a rigid wall, and (iii) an elastic column compressed at a constant speed (“Hoff problem”), in which the time-varying axial force is replaced approximately by its average value over the time duration. Comparison of critical parameters predicted by the proposed criterion with available experimental and simulation data shows that the proposed criterion is in robust reasonable agreement with the known data, which suggests that the proposed simple criterion (p₁t _f =2) can be used to estimate critical parameters for finite time stability of dynamic systems governed by linear equations with constant coefficients.