Influence of the tangential velocity on the compressible Kelvin−Helmholtz instability with nonequilibrium effects

Kelvin−Helmholtz (KH) instability is a fundamental fluid instability that widely exists in nature and engineering. To better understand the dynamic process of the KH instability, the influence of the tangential velocity on the compressible KH instability is investigated by using the discrete Boltzmann method based on the nonequilibrium statistical physics. Both hydrodynamic and thermodynamic nonequilibrium (TNE) effects are probed and analyzed. It is found that, on the whole, the global density gradients, the TNE strength and area firstly increase and decrease afterwards. Both the global density gradient and heat flux intensity in the vertical direction are almost constant in the initial stage before a vortex forms. Moreover, with the increase of the tangential velocity, the KH instability evolves faster, hence the global density gradients, the TNE strength and area increase in the initial stage and achieve their peak earlier, and their maxima are higher for a larger tangential velocity. Physically, there are several competitive mechanisms in the evolution of the KH instability. (i) The physical gradients increase and the TNE effects are strengthened as the interface is elongated. The local physical gradients decrease and the local TNE intensity is weakened on account of the dissipation and/or diffusion. (ii) The global heat flux intensity is promoted when the physical gradients increase. As the contact area expands, the heat exchange is enhanced and the global heat flux intensity increases. (iii) The global TNE intensity reduces with the decreasing of physical gradients and increase with the increasing of TNE area. (iv) The nonequilibrium area increases as the fluid interface is elongated and is widened because of the dissipation and/or diffusion.     

Rayleigh-Taylor instability under multi-mode perturbation: Discrete Boltzmann modeling with tracers

The two-dimensional Rayleigh-Taylor Instability (RTI) under multi-mode perturbation in compressible flow is probed via the Discrete Boltzmann Modeling (DBM) with tracers. The distribution of tracers provides clear boundaries between light and heavy fluids in the position space. Besides, the position-velocity phase space offers a new perspective for understanding the flow behavior of RTI with intuitive geometrical correspondence. The effects of viscosity, acceleration, compressibility, and Atwood number on the mixing of material and momentum and the mean non-equilibrium strength at the interfaces are investigated separately based on both the mixedness defined by the tracers and the non-equilibrium strength defined by the DBM. The mixedness increases with viscosity during early stage but decreases with viscosity at the later stage. Acceleration, compressibility, and Atwood number show enhancement effects on mixing based on different mechanisms. After the system relaxes from the initial state, the mean non-equilibrium strength at the interfaces presents an initially increasing and then declining trend, which is jointly determined by the interface length and the macroscopic physical quantity gradient. We conclude that the four factors investigated all significantly affect early evolution behavior of an RTI system, such as the competition between interface length and macroscopic physical quantity gradient. The results contribute to the understanding of the multi-mode RTI evolutionary mechanism and the accompanied kinetic effects.     

Nonequilibrium and morphological characterizations of Kelvin–Helmholtz instability in compressible flows

We investigate the effects of viscosity and heat conduction on the onset and growth of Kelvin–Helmholtz instability (KHI) via an efficient discrete Boltzmann model. Technically, two effective approaches are presented to quantitatively analyze and understand the configurations and kinetic processes. One is to determine the thickness of mixing layers through tracking the distributions and evolutions of the thermodynamic nonequilibrium (TNE) measures; the other is to evaluate the growth rate of KHI from the slopes of morphological functionals. Physically, it is found that the time histories of width of mixing layer, TNE intensity, and boundary length show high correlation and attain their maxima simultaneously. The viscosity effects are twofold, stabilize the KHI, and enhance both the local and global TNE intensities. Contrary to the monotonically inhibiting effects of viscosity, the heat conduction effects firstly refrain then enhance the evolution afterwards. The physical reasons are analyzed and presented.     

燃烧系统的离散Boltzmann建模与模拟研究进展

燃烧系统的诸多模拟依托于流体建模, 离散Boltzmann方法(discrete Boltzmann method, DBM) 是近年来发展起来的一种新的流体介观建模方法. 本文简要评述DBM发展的两个方向——Navier-Stokes等偏微分方程的数值逼近解法和复杂系统的微介观动理学建模. 主要介绍在燃烧系统模拟方面DBM已有的工作、新近的思路、与传统流体建模的异同以及近期的研究成果. 本文重点传递的信息为: 作为复杂系统微介观动理学建模出现的DBM在模拟过程中同时给出“流动”及其相伴随的、关系最密切的那部分“热动”非平衡效应; 它为燃烧等复杂系统中各类非平衡行为的描述、非平衡信息的提取、非平衡程度的度量提供了一种简洁、有效的方法; 它所提供的热动非平衡测量量有两类: 一类是直接比较分布函数和平衡态分布函数的动理学矩关系得到的, 一类是来自于Chapman-Enskog多尺度分析给出的热传导和黏性项. 基于第二类DBM, 可以实现(燃烧等)一大类复杂流体系统的多尺度物理建模.     

轴向剪切流对Z箍缩等离子体瑞利-泰勒不稳定性的抑制

利用理想磁流体力学模型对有轴向剪切流的Z箍缩等离子体不稳定性进行了分析。给出了可压缩模型的色散关系,分别对可压缩及不可压缩模型中轴向剪切流对Z箍缩等离子体瑞利-泰勒不稳定性的抑制作用进行了比较,讨论了可压缩性对含有轴向剪切流系统不稳定性的影响。结果表明,可压缩性能够减缓瑞利-泰勒P凯尔文-亥姆霍兹(RTPKH)模扰动的增长,因而使得轴向剪切流对系统不稳定性的抑制作用表现得更为突出。计算结果还说明,在RT不稳定性线性增长阶段,等离子体温度较低,使用可压缩模型能够更真实地描述系统的状态。     

Pseudo-potential investigations of structural, elastic and thermal properties of tungsten disilicide

The plane-wave pseudopotential method using the generalized gradient approximation within the density functional theory is used to investigate the structure and bulk modulus of WSi₂. The quasi-harmonic Debye model, using a set of total energy versus cell volume obtained with the plane-wave pseudopotential method, is applied to the study of the elastic properties and vibrational effects. We have analysed the bulk modulus of WSi₂ up to 1600~K. The major trend shows that the WSi₂ crystal becomes more compressible when the temperature rises and the increase of compressibility leads to the decrease of Debye temperature. The predicted temperature and pressure effects on the thermal expansion, heat capacity and Debye temperatures are determined from the non-equilibrium Gibbs functions and compared with the data available.     

Nonequilibrium fluctuations in the Rayleigh-Bénard problem for binary fluid mixtures

We have employed a simple Galerkin-approximation scheme to calculate nonequilibrium temperature and concentration fluctuations in a binary fluid subjected to a temperature gradient with realistic boundary conditions. When a fluid mixture is driven outside thermal equilibrium, there are two instability mechanisms, namely a Rayleigh (stationary) and a Hopf (oscillatory) instability, causing long-ranged fluctuations. The competition of these two mechanisms causes the structure factor associated with the temperature fluctuations to exhibit two maxima as a function of the wave number q of the fluctuations, in particular, close to the convective instability. In the presence of thermally conducting but impermeable walls the intensity of the temperature fluctuations vanishes as q goes to zero, while the intensity of the concentration fluctuations remains finite in the limit of vanishing q. Finally, we propose a simpler small-Lewis-number approximation scheme, which is useful to represent nonequilibrium concentration fluctuations for mixtures with positive separation ratio, even close to (but below) the convective instability.     

Critical behaviour of random compressible magnets

The critical properties of a compressible random magnet are studied using renormalization group methods. Then-component orderparameter is coupled to quenched disorder and to the elastic fluctuations of the anisotropic solid. It is shown, that the critical behaviour of a compressible random magnet is in general the same as that of a random magnet on a rigid lattice. However, if the specific heat exponent of the ideal magnet is positive and the disorder is sufficiently weak, a macroscopic instability may prevent the system in reaching the critical point. The resulting first-order transition may be preceded by pseudocritical behaviour characteristic to pure compressible magnets. The effect of random magnetic fields on the critical properties of compressible magnets is also discussed.     

粘性耗散及变物性对多孔介质中对流换热的影响研究

本文通过理论分析、数值模拟与实验研究的综合,研究了粘性耗散及变物性对多孔介质中对流换热的影响。发现：用理想的等热流边界条件可以较好地数值模拟等热流边界条件下多孔介质中的对流换热;在实际可实现的参数条件下,粘性耗散对多孔介质中对流换热的影响不大;油的变物性对换热的影响很大—热流密度越高,对流换热越强。     

Competing mechanisms of plasma transport in inhomogeneous configurations with velocity shear: the solar-wind interaction with earth's magnetosphere

Two-dimensional simulations of the Kelvin-Helmholtz instability in an inhomogeneous compressible plasma with a density gradient show that, in a transverse magnetic field configuration, the vortex pairing process and the Rayleigh-Taylor secondary instability compete during the nonlinear evolution of the vortices. Two different regimes exist depending on the value of the density jump across the velocity shear layer. These regimes have different physical signatures that can be crucial for the interpretation of satellite data of the interaction of the solar wind with the magnetospheric plasma.     

Effects of Compressibility on the Rayleigh-Taylor Instability in Z-pinch Implosions

In this paper we shall consider the effect of compressibility on the RT instability in Z-pinch implosions, importance is the comparing growth rates of the RT instability for two systems of the compressible and incompressible MHD plasma. For which reason, we shall use as simple model as possible. Obviously, slab geometry is the most simple. For example, in the case of annular plasma implosion, during the linear growth phase of the RT instability there are vacuums at both sides of the annular plasma shell and its thickness is sufficiently smaller than the pinch radius, allowing us to use slab geometry instead of the annular one. For simplicity, we do not consider the effects of the finite Larmor radius and the sheared axial flow which are the important physical mechanisms to compress the RT instabilities.     

微通道脉冲管中He气交替振荡的分子动力学模拟

使用非平衡分子动力学模拟方法分析了微通道脉冲管(MPT)中由正弦速度活塞提供驱动力时He气交替振荡的微观动力学过程,并对MPT的冷却机制进行了分析.结果表明,MPT的压缩和膨胀过程之间存在一个交替的振荡过程,两个过程具有不对称的属性分布,膨胀过程具有比压缩过程更大的轴向压力梯度.当充气压力较低时,循环时间对冷端温度的影响很小,但是当充气压力高于20 bar时,冷端温度对时间较为敏感,随着时间的减少,冷端温度进一步降低,而冷端瞬时平均温度随着充气压力的增加而增加.另外,压比随着时间的减少而增加,并且明显不受充气压力的影响,但它会在MPT的轴向上产生较大的温度梯度.综上所述,在热端使用不同形式的换热器和调相元件会释放或回收额外的声功率.固定工作模式和尺寸参数的MPT具有最佳频率,可以在冷端获得最低的空载温度.仿真结果增进了对脉冲管制冷机的认识,并为微通道脉冲管制冷机的优化设计提供理论支持.     

三维激光烧蚀瑞利-泰勒不稳定性并行计算

给出了用三维激光烧蚀不稳定性程序求解的物理方程组,简要介绍了所用的数值方法.根据所模拟问题的物理性质和数值模拟方法探索了在多处理机上实行并行计算的技术和方法,并行效率达75%,在此基础上进行了三维多模激光烧蚀不稳定性的数值模拟,深化了对有关物理规律的认识.     

Resonant interaction of an electromagnetic wave with an inhomogeneous plasma slab

Resonant interaction at oblique incidence of an electromagnetic wave on an inhomogeneous plasma slab is studied. The time evolution of this interaction is solved numerically from two-fluid equations, adiabatic equation for electron pressure and from Maxwell equations. It is shown that the electromagnetic energy of an incident wave is transformed both into the heat energy and into the energy of plasma oscillations in the direction of density gradient. The distribution of the transformed energy between the heat energy and the energy of plasma oscillations is strongly dependent on the plasma temperature. The ratio of heat energy to the energy of plasma oscillations is growing with growing temperature. The plasma oscillations are generated by magnetic induction of the penetrating wave. In a cold plasma they are generated especially in the overdense region and their frequency is equal to local plasma frequency. The electric field in the direction of plasma gradient has a form of a wave packet whose envelope reaches a maximum at resonance. The characteristic wavelength in the wave packet decreases and the amplitude of the packet increases with the time.     

Time-dependent plasma transport simulation for the study of edge impurity radiation dynamics with magnetic island in large helical device

The dynamics of impurity radiation distribution during detachment transition with edge magnetic island induced by the application of resonant magnetic perturbation (RMP) is studied numerically by solving time-dependent plasma fluid equations together with impurity and neutral transports in a 2D grid system. Computations provide the appearance of a macroscopic structure in the heat transfer along the island separatrix and outside the island. The resultant parallel temperature gradient generates the plasma flow and the density gradient according to the parallel momentum balance. The resulting plasma flow effectively transports impurity towards the X-point region. As a result, the impurity radiation is more intense near the X-point than in the vicinity of the O-point. This leads to the predominant cooling of the region around the X-point. These results are in agreement with experimental observation with RMP application in the large helical device (LHD). The time scale of the thermal condensation instability is found to be of the order of 10 ms.     

Influence of pump beam size on laser diode end-pumped solid state lasers

Based on space-dependent rate equation, the lowest threshold input power for diode end-pumped solid state lasers is obtained for pump spot size w_p → 0. However, when the pump beam waist is decreased, the far-field half-angle increases due to the same beam quality, which will influence the match between the pump volume and mode volume. Most importantly, when the pump beam waist is decreased, the temperature and temperature gradient in the laser rod would be very high due to the resulting heat, which would cause material distortion, variation of refractive index, instability of laser resonator etc. We analyze the temperature distributions with different pump beam sizes and conclude that the condition w_p ≈ w₀ (the radius of the laser beam at the laser rod) is often taken as a rough guide to the optimum case.     

Numerical simulation of anode heat transfer of nitrogen arc utilizing two-temperature chemical non-equilibrium model

A detailed understanding of anode heat transfer is important for the optimization of arc processing technology. In this paper, a two-temperature chemical non-equilibrium model considering the collisionless space charge sheath is developed to investigate the anode heat transfer of nitrogen free-burning arc. The temperature, total heat flux and different heat flux components are analyzed in detail under different arc currents and anode materials. It is found that the arc current can affect the parameter distributions of anode region by changing plasma characteristics in arc column. As the arc current increases from 100 A to 200 A, the total anode heat flux increases, however, the maximum electron condensation heat flux decreases due to the arc expansion. The anode materials have a significant effect on the temperature and heat flux distributions in the anode region. The total heat flux on thoriated tungsten anode is lower than that on copper anode, while the maximum temperature is higher. The power transferred to thoriated tungsten anode, ranked in descending order, is heat flux from heavy-species, electron condensation heat, heat flux from electrons and ion recombination heat. However, the electron condensation heat makes the largest contribution for power transferred to copper anode.     

基于拟沸腾理论的超临界CO2管内传热恶化量纲分析

超临界流体广泛应用于工程技术领域,其流动传热特性对工程设计具有重要意义,但是,由于超临界流体的物理微观和宏观行为的机理尚不清晰,所以其异常的流动传热特性并未得到很好的解决.普遍认为超临界流体在分子尺度上可分为类气和类液两种不同的特性,直到最近通过实验在宏观上监测到超临界水类液和类气之间的转变,且这一过程与拟沸腾理论一致,使得问题逐渐变得清晰.本文基于拟沸腾理论对超临界CO2异常流动传热行为进行了研究,在假设类液和类气转换过程不均匀的情况下,从经典的量纲分析和亚临界过冷沸腾理论模型出发,提出了一个适用于超临界流体拟沸腾换热过程的分析方法.通过引入表征类气膜生长速度与流体主流平均流速之比π=(qw·ρ1)/(G·Δi·ρg)和表征近壁区类气膜温度梯度π13=(qw·βpc·di)/λg两个无量纲数,来表征拟沸腾如何导致传热恶化,解释了超临界CO2竖直向上加热流动过程中的异常换热特性,即较大的类气膜生长速度使近壁区快速聚集了较多的高温流体,而较大的类气膜温度梯度使类气膜覆盖在壁面.当核心的冷类液不能充分润湿热壁面时,传热恶化.新无量纲数较好的诠释了超临界流体拟沸腾诱导传热恶化机制,为超临界拟沸腾传热研究提供了理论依据.     

Hydrodynamical field fluctuations in a non-equilibrium quasi-stationary state due to a temperature gradient : II. Application to light scattering

Within the general theory of hydrodynamical field fluctuations in a system at non-equilibrium due to a temperature gradient, the problem of light scattering on fluctuations in number density is discussed. The effect of the gradient as well as fluctuations in sound velocity and heat conductivity on the Rayleigh-Brillouin spectrum is analyzed. The presence of the temperature gradient and fluctuations modifies the scattered-light spectrum leading to asymmetry in the heights of the Brillouin lines.