Numerical analysis of nitrogen-mixed argon plasma characteristicsand injected particle behavior in an ICP torch for ultrafine powdersynthesis

A numerical model is presented for the analysis of plasma characteristics of an ICP torch and gas mixing effects on the plasma when a nitrogen gas is added into the argon plasma as a carrier or sheath gas at the torch inlet, The fluid equations describing the plasma flow and temperature fields and the diffusions between two different gases are solved along with a magnetic vector potential equation for electromagnetic fields. The trajectory and the temperature change with time for a particle injected into the plasma are also investigated by a plasma-particle interaction model to find out optimum injection conditions for the synthesis of ultrafine nitride ceramic powders, It is found from the calculations that the nitrogen-mixed argon plasma with a nitrogen sheath gas is more favorable than the plasma with a nitrogen carrier gas for the reaction kinetics of nitride synthesis. It is also found that the radial injection through the holes of the tube wall Is preferable to the axial injection at the torch inlet for the complete evaporation of injected particle and the effective chemical reaction of reactant vapor with nitrogen. For the radial injection in an ICP torch of 20 cm in axial length, the optimum injection locations and initial velocities of 50 μm aluminum particles are found for synthesizing aluminum nitride are in the range of 6~12 cm apart from the torch inlet and over 15 m/s, respectively     

����FLUENT�Ķ�ά��Ӧ��ϵ����������ֵģ��

将等离子体作为磁流体,考虑其流体属性和电磁属性,介绍了利用FLUENT软件包并将其进行二次开发,解算电磁场方程、质量连续性方程、动量守恒方程、以及能量守恒方程的数值模拟方法,得到了以磁矢势为表达形式的电磁场分布、温度分布和速度分布.数值模拟了粉末球化所用的感应耦合等离子体炬电磁场分布、温度分布、速度分布.分析了温度分布、速度分布产生的物理原因,为感应耦合等离子体炬球化粉末颗粒提供理论性指导.     

Fluid dynamics analysis of a gas attenuator for X‐ray FELs under high‐repetition‐rate operation

Newtonian fluid dynamics simulations were performed using the Navier–Stokes–Fourier formulations to elucidate the short time‐scale (µs and longer) evolution of the density and temperature distributions in an argon‐gas‐filled attenuator for an X‐ray free‐electron laser under high‐repetition‐rate operation. Both hydrodynamic motions of the gas molecules and thermal conductions were included in a finite‐volume calculation. It was found that the hydrodynamic wave motions play the primary role in creating a density depression (also known as a filament) by advectively transporting gas particles away from the X‐ray laser–gas interaction region, where large pressure and temperature gradients have been built upon the initial energy deposition via X‐ray photoelectric absorption and subsequent thermalization. Concurrent outward heat conduction tends to reduce the pressure in the filament core region, generating a counter gas flow to backfill the filament, but on an initially slower time scale. If the inter‐pulse separation is sufficiently short so the filament cannot recover, the depth of the filament progressively increases as the trailing pulses remove additional gas particles. Since the rate of hydrodynamic removal decreases while the rate of heat conduction back flow increases as time elapses, the two competing mechanisms ultimately reach a dynamic balance, establishing a repeating pattern for each pulse cycle. By performing simulations at higher repetition rates but lower per pulse energies while maintaining a constant time‐averaged power, the amplitude of the hydrodynamic motion per pulse becomes smaller, and the evolution of the temperature and density distributions approach asymptotically towards, as expected, those calculated for a continuous‐wave input of the equivalent power.     

Near‐field optical enhancement by lead‐sulfide quantum dots and metallic nanoparticles for SERS

There is a growing interest in using quantum dots (QDs) and metallic nanoparticles (NPs), both for luminescence enhancement and surface‐enhanced Raman scattering (SERS). Here, we study the electromagnetic‐field enhancement that can be generated by lead‐sulfide (PbS) QDs using three‐dimensional finite‐element simulations. We investigate the field enhancement associated with combinations of PbS QDs with metallic NPs and substrates. The results show that high field enhancement can be achieved by combining PbS QDs with metallic NPs of larger sizes. The ideal size for Ag NPs is 25 nm, providing a SERS enhancement factor of ~5*10⁸ for light polarization parallel to the NP dimer axis and a gap of 0.6 nm. For Au NPs, the bigger the size, the higher is the field for the studied diameters, up to 50 nm. The near‐field values for PbS QDs above metallic substrates were found to be lower compared to the case of PbS QD‐metal NP dimers. This study provides the understanding for the design and application of QDs for the enhancement of near‐field phenomena. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical modelling of a direct current non-transferred thermal plasma torch for optimal performance

A laminar steady-state 2D axisymmetric model of a direct current (DC) thermal plasma torch using a magneto-hydrodynamic approach has been developed. The model takes into account the entire torch system comprising the plasma gas injection, the inner region of the torch, and the jet exiting into the ambient environment. Numerical results are obtained for two different power inputs chosen from published experimental data. The temperature predictions at the torch exit are found in good agreement with experimental results. Velocity analysis of the plasma jet has been presented and the impact of electromagnetic force on jet velocity is analysed. The Lorentz force arising due to the coupling of fluid motion and electromagnetic forces shoots up the jet velocity to significantly high values near the cathode tip. Temperature and velocity profiles are in good agreement with the characteristics of a long laminar plasma jet. An operating value of heat transfer coefficient (h) has been suggested for optimal torch operation, thus ensuring a low anode erosion rate and acceptable thermal efficiency. The argon torch has the maximum temperature and longest jet length among the plasma gases considered.     

A high‐accuracy complex‐phase method of simulating X‐ray propagation through a multi‐lens system

The propagation of X‐ray waves through an optical system consisting of many X‐ray refractive lenses is considered. For solving the problem for an electromagnetic wave, a finite‐difference method is applied. The error of simulation is analytically estimated and investigated. It was found that a very detailed difference grid is required for reliable and accurate calculations of the propagation of X‐ray waves through a multi‐lens system. The reasons for using a very detailed difference grid are investigated. It was shown that the wave phase becomes a function, very quickly increasing with increasing distance from the optical axis, after the wave has passed through the multi‐lens system. If the phase is a quickly increasing function of the coordinates perpendicular to the optical axis, then the electric field of the wave is a quickly oscillating function of these coordinates, and thus a very detailed difference grid becomes necessary to describe such a wavefield. To avoid this difficulty, an equation for the phase function is proposed as an alternative to the equation of the electric field. This allows reliable and accurate simulations to be carried out when using the multi‐lens system. An equation for the phase function is derived and used for accurate simulations. The numerical error of the suggested method is estimated. It is shown that the equation for the phase function allows efficient simulations to be fulfilled for the multi‐lens system.     

Current-voltage characteristics of a non-transferred plasma spray torch

I–V characteristics of a non-transferred DC plasma spray torch operating on argon and argon + nitrogen mixtures are reported. Arc voltage is decreased with increase in arc current and increased with increase in electrode gap. Arc power is higher at higher percentage of nitrogen in argon. Nottingham co-efficients were calculated using numerical method.     

High resolution finite volume scheme for the quantum hydrodynamic equations

The theory of quantum fluid dynamics (QFD) helps nanotechnology engineers to understand the physical effect of quantum forces. Although the governing equations of quantum fluid dynamics and classical fluid mechanics have the same form, there are two numerical simulation problems must be solved in QFD. The first is that the quantum potential term becomes singular and causes a divergence in the numerical simulation when the probability density is very small and close to zero. The second is that the unitarity in the time evolution of the quantum wave packet is significant. Accurate numerical evaluations are critical to the simulations of the flow fields that are generated by various quantum fluid systems.A finite volume scheme is developed herein to solve the quantum hydrodynamic equations of motion, which significantly improve the accuracy and stability of this method. The QFD equation is numerically implemented within the Eulerian method. A third-order modified Osher–Chakravarthy (MOC) upwind-centered finite volume scheme was constructed for conservation law to evaluate the convective terms, and a second-order central finite volume scheme was used to map the quantum potential field. An explicit Runge–Kutta method is used to perform the time integration to achieve fast convergence of the proposed scheme.In order to meet the numerical result can conform to the physical phenomenon and avoid numerical divergence happening due to extremely low probability density, the minimum value setting of probability density must exceed zero and smaller than certain value. The optimal value was found in the proposed numerical approach to maintain a converging numerical simulation when the minimum probability density is 10^?5 to 10^?12. The normalization of the wave packet remains close to unity through a long numerical simulation and the deviations from 1.0 is about 10^?4.To check the QFD finite difference numerical computations, one- and two-dimensional particle motions were solved for an Eckart barrier and a downhill ramp barrier, respectively. The results were compared to the solution of the Schrödinger equation, using the same potentials, which was obtained using by a finite difference method. Finally, the new approach was applied to simulate a quantum nanojet system and offer more intact theory in quantum computational fluid dynamics.     

Effect of the Electron‐Beam Self‐Fields on Gain in an Optical Wiggler Pumped Free‐Electron Laser

The effects of self fields on gain for a free‐electron lasers (FELs) with electromagnetic‐wave wiggler and an axial guide magnetic field is presented. The relativistic equation of motion for a single electron for all relevant fields, including wiggler, self‐fields and axial guide magnetic field has been solved. Two classes of possible single‐particle trajectories in this configuration are found. Result of the numerical calculation shown that the relativistic part of group I (group II) orbits decreases (increases) monotonically with the axial field. The gain equations for the FEL configuration by adding the effect of self‐fields have been derived. The numerical calculation has been employed to analysis the gain induced by the effects of the self‐fields. It is shown that, for group I orbits the gain decreases in the presence of self‐fields and the gain decrement increases with increasing axial guide magnetic field, while for group II orbits the self‐fields enhances the gain. The gain decrement and enhancement are due to diamagnetic and paramagnetic effects of the self‐magnetic field, respectively. The comparison of the gain for electromagnetic‐wave wiggler with the gain in helical wiggler has been done (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modelling of an argon plasma flow

One-fluid MHD equations are numerically solved for an axisymmetric flow of argon inside and outside a discharge chamber of a cascaded plasma torch. Arc currents of 100 and 150 A, flow rates of 10 and 30 slpm, and two different optical thicknesses are assumed. The flow is shown to be only weakly compressible, of the Mach number below 0.3. The calculated torch powers range from 3 to 7 kW, the energy lost by radiation is (1–3) kW. Axial velocities of (200–600) m/s and temperatures of (12 500–14 400) K are found at the torch nozzle exit and compared with available experimental data. Inside the chamber, several recirculation zones are predicted. The physical model and the chamber shape used can readily be extended for a hybrid plasma torch with the simultaneous stabilization of the electric arc inside the chamber by argon and, e.g., water swirl.Dedicated to Professor Jií Horáek on the occasion of his 60th birthday.     

直流非转移弧热等离子体喷枪加热效率分析

利用计算机模拟技术,采用有限容积方法,针对自制的低能耗、高效率内送粉等离子喷涂设备,对喷枪内部影响喷涂热效率的因素进行了分析。研究结果表明通过适当提高进气流量、提高混合气体中的氮气或氢气含量和减少喷枪通道长度,均能提高等离子喷枪热效率。而电流的大小以及喷枪的通道直径对热效率无显著影响。     

基于格子Boltzmann方法的超声速非平衡流模拟

基于D1Q4可压缩格子Boltzmann模型,按照流通矢量分裂方法的思路,采用坐标旋转技术构造求解三维带化学反应Navier-Stokes方程对流通量求解器.结合有限体积法求解三维化学非平衡流Navier-Stokes方程,采用时间算子分裂算法解决化学反应刚性问题,数值模拟超声速化学非平衡流的三个经典算例.数值结果表明:在高马赫数下,采用D1Q4可压缩格子Boltzmann模型构造的三维对流通量求解器数值模拟中没有出现非物理解,同时在超声速化学非平衡流场中正确分辨激波、燃烧波等物理现象,精度和分辨率均较高,验证了本文构造的三维对流通量求解器的可靠性,拓宽了D1Q4可压缩格子Boltzmann模型的应用范围,为计算超声速化学非平衡流提供一种新方法.     

Free Mesh Method: fundamental conception, algorithms and accuracy study

The finite element method (FEM) has been commonly employed in a variety of fields as a computer simulation method to solve such problems as solid, fluid, electro-magnetic phenomena and so on. However, creation of a quality mesh for the problem domain is a prerequisite when using FEM, which becomes a major part of the cost of a simulation. It is natural that the concept of meshless method has evolved. The free mesh method (FMM) is among the typical meshless methods intended for particle-like finite element analysis of problems that are difficult to handle using global mesh generation, especially on parallel processors. FMM is an efficient node-based finite element method that employs a local mesh generation technique and a node-by-node algorithm for the finite element calculations. In this paper, FMM and its variation are reviewed focusing on their fundamental conception, algorithms and accuracy.     

Ar介质感应耦合等离子加热器能量转化与流动特性

对感应耦合等离子(inductively coupled plasma,ICP)加热器内能量转化过程与分布规律、 流动特性的研究和认识能够为高频等离子加热器的设计提供理论指导,同时能够为加热器向大功率、多介质、广适用方向的发展提供支撑.基于二维轴对称、层流流动和局部热力学平衡等假设条件,利用COMSOL对百千瓦级Ar介...     

Viscous quark‐gluon plasma in the early universe

In the present work a study is given for the evolution of a flat, isotropic and homogeneous Universe, which is filled with a causal bulk viscous cosmological fluid. We describe the viscous properties by an ultra‐relativistic equation of state, and bulk viscosity coefficient obtained from recent lattice QCD calculations. The basic equation for the Hubble parameter is derived by using the energy equation obtained from the assumption of the covariant conservation of the energy‐momentum tensor of the matter in the Universe. By assuming a power law dependence of the bulk viscosity coefficient, temperature and relaxation time on the energy density, we derive the evolution equation for the Hubble function. By using the equations of state from recent lattice QCD simulations and heavy‐ion collisions we obtain an approximate solution of the field equations. In this treatment for the viscous cosmology, no evidence for singularity is observed. For example, both the Hubble parameter and the scale factor are finite at t = 0, where t is the comoving time. Furthermore, their time evolution essentially differs from the one associated with non‐viscous and ideal gas. Also it is noticed that the thermodynamic quantities, like temperature, energy density and bulk pressure remain finite. Particular solutions are also considered in order to prove that the free parameter in this model does qualitatively influence the final results.     

Numerical and experimental investigation of seam welding with a pulsed laser

The present work develops two numerical models to compute thermal phenomena during pulsed laser welding. The first one which is based on finite difference model calculates the welding width and its penetration by utilizing heat transfer equations. Parametric design capabilities of the finite element code ANSYS were also employed for the simulation of the second model. The transient temperature profiles, the fusion dimensions and the heat affected zones (HAZ) have been calculated here. The thermo-physical properties are dependent on temperature and so a nonlinear solution is employed. It is found that the temperature profile and penetration depth are strongly dependent on the pulse parameters of laser beam. Finally, the results of the two models and the experimental outcomes are compared.     

Interaction phenomena of ion‐acoustic quasi‐solitons and classical‐solitons in three component plasma

The space–time evolution of the cnoidal‐soliton solution, characteristics of the quasi‐soliton solution of Korteweg‐de‐Vries (KdV) equation, and the interaction phenomena of ion‐acoustic waves (IAWs) are investigated in a plasma system consisting of positive and negative ions with superthermal electrons. To do this, and (Ar⁺, F^?) plasmas are considered and two‐sided KdV equations (KdVEs) are derived applying the extended Poincaré‐Lighthill‐Kuo (ePLK) method. The effects on wave structures, potential profiles, and propagation characteristics with plasma parameters of the cnoidal‐wave, quasi‐soliton solution, and head‐on collision phenomena of IAWs are presented graphically. It was found that the superthermality parameter and the mass ratio of ions play a significant role in the head‐on collision between soliton and standing cnoidal wave and reveal that the collision is elastic and both waves change their phase shifts due to collision. Moreover, the superthermality parameters are also responsible for the production of compressive and rarefactive phase shifts in overtaking collision processes between right travelling classical soliton (CS) and cnoidal wave (CW) and reduced the amplitudes of IAWs. It was also found that a new wave is created with a high amplitude in the interacting region during collision depending on the plasma parameters.     

Implementation of the moving particle semi-implicit method on GPU

The Moving Particle Semi-implicit (MPS) method performs well in simulating violent free surface flow and hence becomes popular in the area of fluid flow simulation. However, the implementations of searching neighbouring particles and solving the large sparse matrix equations (Poisson-type equation) are very time-consuming. In order to utilize the tremendous power of parallel computation of Graphics Processing Units (GPU), this study has developed a GPU-based MPS model employing the Compute Unified Device Ar...