FLUENT软件在射流模拟中的应用

介绍了利用CFD软件对二维射流作用下的吸气流动的模拟,通过与实验研究得出的结果进行比较,结果吻合得很好。同时,与理论分析方法进行了比较,证明利用商业软件对射流作用的吸气运动进行研究是可行的。

静电雾化中圆锥-射流模式的数值模拟

在液体高压静电雾化过程中,如何匹配众多参数得到稳定的均匀液滴是该技术的关键和难点。为此,针对以上问题,建立了均匀液滴喷射过程的流场计算模型,应用FLUENT软件,通过参数的合理匹配,模拟得到圆锥-射流模式下的均匀液滴流和稳定雾场,并对静电雾化过程的流场和雾场的相关特性进行分析。结果表明,在场强很大的区域内,雾滴轴向速度也较大;在接收板附近,外部流动的空气对雾滴产生加速作用,使雾场区域发散。研究结果对静电雾化效果的优化有一定的参考价值。

一种新的模拟固体火箭发动机射流铝颗粒燃烧的方法

在总结以往火箭射流数值模拟的基础上,采用多组分有限速率化学反应模型模拟气体在喷管后的复燃,并引入拉格朗日颗粒轨道模型来模拟由Al和Al_2O_3组成的混合颗粒。通过快速欧拉模型得到Al_2O_3烟尘的速度场,然后通过质量、动量、能量源项将气体相、颗粒相和烟尘相三者联立耦合计算来模拟铝颗粒的燃烧和燃烧产生的烟尘颗粒对整个射流结构的影响,为以后更准确地计算射流红外辐射特性提供了理论工具。

大气压非平衡等离子体射流活性粒子产生机制

大气压非平衡等离子体中的高能电子和大量活性粒子如O、OH等能够同污染物发生反应,达到促进污染物降解的目的。因此通过数值模拟和实验的方法对大气压脉冲直流等离子体射流推进以及OH产生机理进行了研究。模拟结果表明:介质管内表面附近形成的鞘层将等离子体通道和介质隔离开来,介质边缘的强电场增强了气体混合区域的电离,因此射流头部呈现出环状结构。Penning电离在电离过程中起主导作用并提高了等离子体通道的电导率。电子分解水分子,电子和水离子的复合以及O(1D)分解水分子是产生OH最主要的3个反应,其中电子分解水分子在3者中作用最大。工作气体中掺入额外的N2、O2和水蒸气能够小幅度提高OH的密度是因为Penning电离增加了电子密度。由于空气的扩散,尽管掺入水蒸气能够使管内OH密度有显著提高,该方法对提高管外OH密度效果并不明显。选择合适的脉冲周期数是控制射流活性的有效手段。

一种新型串联战斗部前级装药结构的数值模拟研究

针对串联战斗部前级装药结构兼顾大的侵深与扩孔,设计了一种阶梯状药型罩。采用数值模拟的方法研究了阶梯状药型罩射流的形成机理以及对45号钢的侵彻能力,并与锥角药型罩形成的射流参数以及对45号钢的侵彻深度、扩孔大小进行对比。研究结果表明:相同装药条件下,阶梯状药型罩比锥角药型罩的侵彻深度提升约25%,对45号钢板的扩孔直径扩大约1.5倍;并且形成杵体较小,更加符合串联战斗部前级装药结构的特性。

等离子喷涂-物理气相沉积高温防护涂层研究进展

等离子喷涂-物理气相沉积(PS-PVD)是基于低压等离子喷涂发展起来的一种新型多功能薄膜及涂层制备技术。由于其独特的等离子射流特征,可实现气液固多相涂层沉积,获得非视线沉积。文中首先介绍了国内外PS-PVD技术等离子体数值模拟和在线检测技术的研究现状,其次讨论了PS-PVD羽-柱状结构热障涂层的形成机制及与传统热障涂层在热导率、抗冲蚀等性能方面的差异,阐述了PS-PVD技术制备环境障涂层的研究进展,最后对PS-PVD技术沉积高温防护涂层的优势和存在的问题进行了总结。

Numerical simulation of fluid dynamics in rare earth chloride solution in jet-flow pyrolysis reactor

Rare earth oxide was prepared via direct pyrolysis of rare earth chloride solution. Based on this technique, a new-type jet-flow pyrolysis reactor was designed, and then the fluid dynamics (pressure and velocity) inside the reactor was numerically simulated using a computational fluid dynamics method. The self-produced pressure (p) and the fuel inlet velocity (v) satisfied a quadratic function,p=0.06v2+0.23v?4.49. To fully utilize the combustion-generated heat in pyrolysis of rare earth chloride, an appropriate external pressure p=v2+3v?4.27 should be imposed at the feed inlet. The 1.25- and 1.5-fold increase of feed inlet diameter resulted in decline of adsorption dynamic pressure, but the intake of rare earth chloride increased by more than 30% and 60%, respectively. The fluid flow in the reactor was affected by the feeding rate; the fluid flow peaked near the throat of venturi and gradually smoothed down at the jet-flow reactor’s terminal along with the sharp decline of feeding rate.

Mathematical modelling and numerical optimization of particle heating process in copper flash furnace

A mathematical model of the particle heating process in the reaction shaft of flash smelting furnace was established and the calculation was performed.The results indicate that radiation plays a significant role in the heat transfer process within the first 0.6 m in the upper part of the reaction shaft,whilst the convection is dominant in the area below 0.6 m for the particle heating.In order to accelerate the particle ignition,it is necessary to enhance the convection,thus to speed up the particle heating.A high-speed preheated oxygen jet technology was then suggested to replace the nature gas combustion in the flash furnace,aiming to create a lateral disturbance in the gaseous phase around the particles,so as to achieve a slip velocity between the two phases and a high convective heat transfer coefficient.Numerical simulation was carried out for the cases with the high-speed oxygen jet and the normal nature gas burners.The results show that with the high-speed jet technology,particles are heated up more rapidly and ignited much earlier,especially within the area of the radial range of R=0.3−0.6 m.As a result,a more efficient smelting process can be achieved under the same operational condition.

Using a Time-domain Higher-order Boundary Element Method to Simulate Wave and Current Diffraction from a 3-D Body

To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.

Flow Behaviors of Gas-Solid Injector by 3D Simulation with Kinetic Theory of Granular Flow

A computational study on the flow behavior of a gas-solid injector by Eulerian approach was carried out. The gas phase was modeled with k-ε turbulent model and the particle phase was modeled with kinetic theory of granular flow. The simulations by Eulerian two-fluid model （TFM） were compared with the corresponding results by discrete element method （DEM） and experiments. It was showed that TFM simulated results were in reasonable agreement with the experimental and DEM simulated results. Based on TFM simulations, gas-solid flow pattern, gas velocity, particle velocity and the static pressure under different driving jet velocity, backpressure and convergent section angle were obtained. The results showed that the time average axial gas velocity sharply decreased and then slightly increased to a constant value in the horizontal conveying pipe. The time average axial particle velocity increased initially and then decreased, but in the outlet region of the convergent section the particle velocity remarkably increased once more to the maximal value. As a whole, the static pressure distribution change trends were found to be independent on driving gas velocity, backpressure and convergent section angle. However, the static pressure increased with increase of convergent section angle and gas jet velocities. The difference of static pressure to backpressure increased with increasing backpressure.

Simulation of turbulent combustion in DLR Scramjet

Turbulent combustion in a DLR (German Aerospace Center) Scramjet engine was simulated using the newly-proposed Partially Resolved Numerical Simulation (PRNS) procedure. The PRNS procedure uses temporal filtering to define large-scale turbulence, and the model developed to account for unresolved scales is grid independent. No problem about inner commutation error and inconsistencies will arise from the PRNS, while such issues are of concern in traditional Large Eddy Simulation (LES) methods. The mean results have good agreement with the experiment data and the flow structures with small scales are well resolved.

Numerical Simulation of Direct-contact Condensation from a Supersonic Steam Jet in Subcooled Water

The phenomenon of direct-contact condensation,used in steam driven jet injectors,nuclear reactor emergency core cooling systems and direct-contact heat exchangers,was investigated computationally by introducing a thermal equilibrium model for direct-contact condensation of steam in subcooled water.The condensation model presented was a two resistance model which takes care of the heat transfer process on both sides of the interface and uses a variable steam bubble diameter.The injection of supersonic steam jet in subcooled water tank was simulated using the Euler-Euler multiphase flow model of Fluent 6.3 code with the condensation model incorporated. The findings of the computational fluid dynamics（CFD） simulations were compared with the published experimental data and fairly good agreement was observed between the two,thus validating the condensation model.The results of CFD simulations for dimensionless penetration length of steam plume varies from 2.73-7.33,while the condensation heat transfer coefficient varies from 0.75-0.917 MW·（m ^2 ·K）^ -1 for water temperature in the range of 293-343 K.

Quadrature Method of Moments for Nanoparticle Coagulation and Diffusion in the Planar Impinging Jet Flow

A computational model combining large .eddy simulation with quadrature moment method was em-ployed to study nanoparticle evolution in a confined impinging jet. The investigated particle size is limited in the transient regime, and the particle collision kernel was obtained by using the theory of flux matching. The simulation was validated by comparing it with the experimental results. The numerical results show coherent structure acts to dominate particle number intensity, size and polydispersity distributions, and it also induce particle-laden iet to be diluted by .the ambient.The evolution of particle dynarnics in.the impinging jet flow are strongly related to the Rey-nolds number and nozzle-to-plate distance, and their relationships were analyzed.

Lattice Boltzmann simulation of a laminar square jet in cross flows

A three-dimensional, nineteen-velocity(D3Q19) Lattice Boltzmann Method(LBM) model was developed to simulate the fluid flow of a laminar square jet in cross flows based on the single relaxation time algorithm. The code was validated by the mathematic solution of the Poiseuille flow in a square channel, and was further validated with a previous well studied empirical correlation for the central trajectory of a jet in cross flows. The developed LBM model was found to be able to capture the dominant vortex, i.e. the Counter-rotating Vortex Pair(CVP) and the upright wake vortex. Results show that the incoming fluid in the cross flow channel was entrained into the leeside of the jet fluid, which contributes to the blending of the jet. That the spread width of the transverse jet decreases with the velocity ratio. A layer-organized entrainment pattern was found indicating that the incoming fluid at the lower position is firstly entrained into the leeside of the jet, and followed by the incoming fluid at the upper position.

Analysis of chamber effect on intermittent pulsation

A numerical study was conducted to seek an optimized dimension of jet chamber in the pulsating impinging flow.The flow and heat transfer effect of the pulsation flow through a jet chamber was investigated.The numerical results indicate that heat transfer effective enhances near the stagnation region for the intermittent pulsed flow with jet chamber compared to that without jet chamber.Simulations of the flow through a jet chamber show that the heat transfer rate on the impingement surface is highly dependent on the velocity at the position which is really close to target surface.Examination of the velocity field suggests that the velocity exists a maximum value as the axis distance increases.In addition,the velocity at the jet hole is enlarged by the jet chamber due to the entrainment effect,and the velocity is amplified even greater as the size of the jet chamber becomes bigger.Nevertheless,the velocity declines quickly while the flow axis distance is more than a certain range,leading to poor heat transfer.Thus,intermittent pulsed flow with jet chamber is suggested as a method of improving heat transfer by employing larger dimensions of jet chamber for appropriate jet-to-surface spacing.

Large-eddy Simulation of Bubble-Liquid Confined Jets

The Large-eddy simulation (LES) with two-way coupling is used to study bubble-liquid two-phase confined multiple jets discharged into a 2D channel.The LES results reveal the large-eddy vortex structures of both liquid flow and bubble motion,the shear-generated and bubble-induced liquid turbulence,and indicate much stronger bubble fluctuation than that of the liquid,the enhancement of liquid turbulence by bubbles.Both shear and bubble-liquid interaction are important for the liquid turbulence generation in the case studied.

Analysis of disagreement between numerically predicted and experimental heat transfer data of impinging jet

The method of numerical simulation was applied to investigate the effects of jet impinging plate thickness and its thermal conductivity on the local heat flux distribution along the impinging plate. The results show that the two factors have great effects on the heat flux distribution. The non-uniformity of the local heat-flux on the impinging plate surface gets more profound as the plate becomes thicker and thermal conductivity gets larger. When Reynolds number is 5000, the ratio of nozzle-to-plate spacing to nozzle diameter is 5 and thermal conductivity is 16W/(m·K), and even for the plate with only 25μm in thickness, the non-uniformity of the heat flux cannot be neglected. When the plate thickness is 50 μm, only when thermal conductivity is as small as 1W/(m·K), the heat flux curve can be approximately treated as an iso-heat-flux boundary. In the experimental research, a real non-iso-heat-flux boundary is treated as an iso-heat-flux boundary, which would result in under-estimated Nusselt number value in the stagnation zone and an over-estimated value outside. Such an experimental Nusselt number distribution is taken to evaluate turbulent model, and the conclusion would be drawn that the turbulent model over-predicts the stagnation heat transfer. This is one of the important reasons why many literatures reported that k-ε turbulent model dramatically over-predicts the impinging jet heat transfer in the stagnation region.

Effect of mesoscale wind stress-SST coupling on the Kuroshio extension jet

Effect of mesoscale wind stress-SST coupling on the Kuroshio extension jet is studied using the Regional Ocean Modeling System. The mesoscale wind stress perturbation( τ_(MS)) is diagnostically determined from modelled mesoscale SST perturbation(SST_(MS)) by using their empirical relationship derived from corresponding observation. From comparing two experiments with and without the τ_(MS) feedback, it is found that the interactively represented τ_(MS)-SST_(MS) coupling can modulate the kinetic energy along the Kuroshio extension jet, with little effect on the Kuroshio pathway. Similar results are also obtained in three additional sensitivity experiments, which consider half strength of the τ_(MS), and the momentum flux and heat flux effect induced by τ_(MS), respectively. That means simply taking into account the τ_(MS)-SST_(MS) coupling has little effect on improving the simulation of the Kuroshio Current system.

Numerical Analysis of Slag Splashing in a Steelmaking Converter

Some variables that influence the slag splashing phenomenon in an oxygen steelmaking converter are numerically analyzed in this work. The effect of lance height, jet velocity, jet exit angle and slag viscosity on the washing and ejection mechanisms of slag splashing is studied employing transient two-dimensional computational fluid dynamics simulations. A parameter here called average slag volume fraction is proposed for the quantitative evaluation of the slag splashing efficiency. Besides, a qualitative comparison is made between the computational fluid dynamics results and physical model results from literature.

Rarefied planar jets into vacuum

This paper presents a fundamental gaskinetic study on high speed rarefied jets expanding into vacuum from a cluster of planar exits.Based on the corresponding exact expressions for one planar jet,this paper straightforwardly derives the combined multiple jet flowfield solutions of density and velocity components,however,for the combined temperature and pressure solutions,extra attention shall be practiced.Several direct simulation Monte Carlo simulation results are provided and they validate these analytical solutions of rarefied planar jet flows.