Numerical simulation of transient processes in hydroturbines

A method for calculation of unsteady 3D turbulent flows in hydro turbines of power plants developed for simulation of the transient processes is presented herein. The method is based on joint solution to the Reynolds-averaged Navier—Stokes equations for incompressible fluid flow written for moving mesh, the equation of runner rotation and the system of 1D equations describing propagation of elastic hydraulic shock in the flow domain. The exchange in flow parameters between the penstock and hydro turbine regions is considered in this approach. Results of transient processes simulation are presented for several modes: start-up to the turbine regime, instantaneous load shedding, and output power decrease. Comparison with experimental data is performed.     

管道后传声数值模拟的不稳定波抑制

在管道后传声的数值模拟中,必须考虑平均流剪切层的散射效应,然而在非均匀剪切流动下时域求解线化欧拉方程会面临Kelvin-Helmholtz不稳定波产生和放大的难题。已有的不稳定波抑制技术通常很难获得令人满意的结果。本文采用一种混合方法,首先引入有限时段的宽频声源波包将声波和不稳定波分离,进而采用声源滤波器技术对不稳定波进行抑制。数值验证算例选择半无限长轴对称环形硬壁直管道,采用计算气动声学方法时域求解2.5维线化欧拉方程,无背景流动的数值解与解析解符合很好,验证了程序的精度与可靠性,非均匀流动算例则表明所采用波包加声源滤波器混合方法对不稳定波抑制效果明显,对声场影响很小,充分显示了该方法的精度与可行性。     

叶栅内定常及非定常粘性流动数值模拟

本文给出了一个模拟叶栅内准三维定常和非定常粘性流动的数值方法。对于定常流动,采用ＴＶＤ Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式和代数湍流模型求解雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,使用当地时间步长和多网格技术使计算加速收敛到定常状态;对于非定常流动,使用双时间步长和全隐式离散,采用与求解定常流动相似的多网格方法求解隐式离散方程。文中给出了ＶＫＩ透平叶栅内的定常流结果和１．５级透平叶栅内的非定常数值结果。     

Computational Model of Thermal‐Fluid Flow in a Radio‐Frequency Plasma Torch

The paper aims to clarify the modelling results concerning the heat transfer and fluid flow in a radio‐frequency plasma torch with argon at atmospheric pressure. Fluid numerical simulation requires the coupling of magnetohydrodynamics (MHD) and thermal phenomena. This model combines Navier–Stokes equations with the Maxwell's equations for compressible fluid and electromagnetic phenomena successively. A numerical formulation based on the finite element method is used. In this study, fluid flow and temperature equations are simultaneously solved (direct method, instead of using the indirect method) using a finite elements method (FEM) for optically thin argon plasmas under the assumptions of local thermodynamic equilibrium (LTE) and laminar flow. Appropriate boundary conditions are given, and nonlinear parameters such as the thermal and electrical conductivity of the gas and input power used in the simulation are detailed. We have found that the source of power is located on the torch wall in this type of inductive discharge. The center can be heated by conduction and convection via electromagnetic phenomena (power loss and Lorentz force). (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of steady cavitating flow of viscous fluid in a Francis hydroturbine

Numerical technique was developed for simulation of cavitating flows through the flow passage of a hydraulic turbine. The technique is based on solution of steady 3D Navier—Stokes equations with a liquid phase transfer equation. The approch for setting boundary conditions meeting the requirements of cavitation testing standard was suggested. Four different models of evaporation and condensation were compared. Numerical simulations for turbines of different specific speed were compared with experiment.     

Regularities of heat- and mass transfer and evaporation of activators of thermoemissive cathodes

The problem of heat- and mass transfer and evaporation of emission-activating elements of refractory rod cathodes is formulated and solved in the 2D approximation with nonlinear boundary conditions. Time-dependent heat conduction equations, current continuity equations, as well as the diffusion equations and equations describing evaporation of activators are formulated simultaneously and solved numerically taking into account the temperature dependence of the electrode material. The thermal state of cathode assemblies depending on extrinsic parameters of plasma units is analyzed. The results of simulation of the thermal state and heat- and mass transfer in cathode assemblies with different geometries are considered.     

高超声速层流绕三维压缩拐角的数值模拟

通过数值模拟研究了高超声速来流绕过压缩拐角的层流分离三维流动特性.数值方法采用三维N-S方程,结合2阶精度Roe格式以及分区结构网格有限体积法进行离散.数值模拟的空间激波结构与实验纹影结果符合较好;激波/边界层干扰区内3条纵向线上的计算压力分布与实验结果进行了对比分析,计算获得在三维楔侧面存在低压力区,与实验结果反映的规律一致,计算结果表明低压力区是由楔体侧缘尖端发起的二次涡的抽吸作用造成的.此外,在楔体后端尾流区的低压沿边界层内的亚声速区往上游传递了一定距离.      

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

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

The MAST FV/FE scheme for the simulation of two-dimensional thermohaline processes in variable-density saturated porous media

A novel methodology for the simulation of 2D thermohaline double diffusive processes, driven by heterogeneous temperature and concentration fields in variable-density saturated porous media, is presented. The stream function is used to describe the flow field and it is defined in terms of mass flux. The partial differential equations governing system is given by the mass conservation equation of the fluid phase written in terms of the mass-based stream function, as well as by the advection–diffusion transport equations of the contaminant concentration and of the heat. The unknown variables are the stream function, the contaminant concentration and the temperature. The governing equations system is solved using a fractional time step procedure, splitting the convective components from the diffusive ones. In the case of existing scalar potential of the flow field, the convective components are solved using a finite volume marching in space and time (MAST) procedure; this solves a sequence of small systems of ordinary differential equations, one for each computational cell, according to the decreasing value of the scalar potential. In the case of variable-density groundwater transport problem, where a scalar potential of the flow field does not exist, a second MAST procedure has to be applied to solve again the ODEs according to the increasing value of a new function, called approximated potential. The diffusive components are solved using a standard Galerkin finite element method. The numerical scheme is validated using literature tests.     

Kinetic modeling of a high power fast-axial-flow CO2 laser with computational fluid dynamics method

A new computational fluid dynamics (CFD) method for the simulation of fast-axial-flow CO2 laser is developed.The model which is solved by CFD software uses a set of dynamic differential equations to describe the dynamic process in one discharge tube.The velocity,temperature,pressure and turbulence energy distributions in discharge passage are presented.There is a good agreement between the theoretical prediction and the experimental results.This result indicates that the parameters of the laser have significant effect on the flow distribution in the discharge passage.It is helpful to optimize the output of high power CO2 laser by mastering its kinetic characteristics.     

On the accuracy of direct forcing immersed boundary methods with projection methods

Direct forcing methods are a class of methods for solving the Navier–Stokes equations on nonrectangular domains. The physical domain is embedded into a larger, rectangular domain, and the equations of motion are solved on this extended domain. The boundary conditions are enforced by applying forces near the embedded boundaries. This raises the question of how the flow outside the physical domain influences the flow inside the physical domain. This question is particularly relevant when using a projection method for incompressible flow. In this paper, analysis and computational tests are presented that explore the performance of projection methods when used with direct forcing methods. Sufficient conditions for the success of projection methods on extended domains are derived, and it is shown how forcing methods meet these conditions. Bounds on the error due to projecting on the extended domain are derived, and it is shown that direct forcing methods are, in general, first-order accurate in the max-norm. Numerical tests of the projection alone confirm the analysis and show that this error is concentrated near the embedded boundaries, leading to higher-order accuracy in integral norms. Generically, forcing methods generate a solution that is not smooth across the embedded boundaries, and it is this lack of smoothness which limits the accuracy of the methods. Additional computational tests of the Navier–Stokes equations involving a direct forcing method and a projection method are presented, and the results are compared with the predictions of the analysis. These results confirm that the lack of smoothness in the solution produces a lower-order error. The rate of convergence attained in practice depends on the type of forcing method used.     

二维槽道湍流拟序结构的大涡模拟

本文采用大涡模拟的方法,对二维槽道湍流流动进行了数值模拟。采用Chorin的分步投影法求解大尺度涡运动的Navier-Stokes方程,小尺度涡采用三种亚格子（SGS）模式分别模拟,给出了不同亚格子涡粘性模式下的模拟结果。对固壁面采用了壁函数。模拟结果再现了二维槽道流动拟序结构的发展演变过程。通过对不同入口速度下的瞬态流场的比较,揭示了入口速度分布对流场的影响。     

剪切流动条件下液滴变形和断裂的数值模拟

本文采用扩散界面法研究了剪切流动条件下悬浮液滴变形和断裂的动力学机制。控制方程采用考虑表面张力影响的Navier-Stokes-Cahn-Hilliard方程描述。计算网格采用均匀矩形交错网格。采用基于压力增量的近似投影法计算 Navier-Stokes方程,采用完全近似多重网格法计算Cahn-Hilliard方程。稳态液滴变形规律及液滴拉伸断裂过程的计算结果与试验结果符合较好,表明本文模型能够很好的研究液滴变形及断裂机理。     

A diffuse-interface method for two-phase flows with soluble surfactants

A method is presented to solve two-phase problems involving soluble surfactants. The incompressible Navier-Stokes equations are solved along with equations for the bulk and interfacial surfactant concentrations. A non-linear equation of state is used to relate the surface tension to the interfacial surfactant concentration. The method is based on the use of a diffuse interface, which allows a simple implementation using standard finite difference or finite element techniques. Here, finite difference methods on a block-structured adaptive grid are used, and the resulting equations are solved using a non-linear multigrid method. Results are presented for a drop in shear flow in both 2D and 3D, and the effect of solubility is discussed.     

一种新型宽带光纤拉曼放大器的理论研究

主要比较分析了光子晶体光纤在设计宽带光纤拉曼放大器上的优越性。简化了多波长宽带光纤拉曼放大器功率耦合方程,分两步来确定各泵浦波的频率及输入功率的大小。首先通过模拟煺火算法迭代出满足条件的泵浦波频率,再利用平均功率分析方法,采用四阶阿当迭代方法计算出各泵浦波输入功率的大小,设计出了具有较宽平坦增益带宽的光纤拉曼放大器。分别对光子晶体光纤和色散位移光纤进行了模拟计算,仿真结果表明光子晶体光纤可以用来设计成短长度、高效的光纤拉曼放大器。     

高功率微波与等离子体相互作用理论和数值研究

研究高功率微波与等离子体的相互作用,对于微波放电和电磁兼容研究均具有重要意义.基于波动方程、等离子体的流体力学方程以及波尔兹曼方程,建立高功率微波脉冲与等离子体相互作用的理论模型,并结合等离子体的特征参数,采用时域有限差分方法分析了等离子体电子密度和高功率微波传输特性的变化.结果表明,由于高功率微波的电子加热作用,等离子体中的非线性效应明显,发生击穿使得等离子体电子密度增大,从而导致微波的反射增强,透过率降低.所提出的模型和相关结果对于高功率微波和电磁脉冲防护具有指导意义.     

Ghost Fluid方法与双介质可压缩流动计算

应用带有Isobaric修正的GhostFluid方法配合LevelSet方法计算可压缩双介质无粘流动.该方法可以消除计算流体界面时所产生的数值跳动和耗散,且编程上比界面跟踪法简单.应用WENO格式数值求解欧拉方程和LevelSet方程,对由刚性气体状态方程所支配的一二维双介质流动进行数值计算,得到了分辨率较高的计算结果.     

Heat Transfer and Thermocapillary Convection during the Laser Deposition of Metal Powders Implemented in Additive Technologies

Heat-transfer- and thermocapillary-convection macroprocesses observed during direct laser metal deposition (DLMD) with coaxial powder injection are examined. The study is performed using the 3D mathematical model incorporating self-consistent equations for free surface evolution, heat transfer, and hydrodynamics, which allow for powder-particle embedding into the thermocapillary convection zone under DLMD. The processes under consideration refer to the main ones underlying additive laser technologies, which determine the microstructural properties and quality of synthesized parts. The convection-diffusion equations are numerically solved using the final volume method. Calculations are carried out for the thermocapillary convection of H13 steel powder. The influence of laser-radiation characteristics (power, scanning rate, intensity distribution in the beam) and the powder-mass flow velocity on temperature fields, the structure of convective melt flow (including a maximum melt velocity), and the geometric characteristics (height and width) of the object formed is investigated.     

Modeling nonlinear ultrasound propagation in heterogeneous media with power law absorption using a k-space pseudospectral method

The simulation of nonlinear ultrasound propagation through tissue realistic media has a wide range of practical applications. However, this is a computationally difficult problem due to the large size of the computational domain compared to the acoustic wavelength. Here, the k-space pseudospectral method is used to reduce the number of grid points required per wavelength for accurate simulations. The model is based on coupled first-order acoustic equations valid for nonlinear wave propagation in heterogeneous media with power law absorption. These are derived from the equations of fluid mechanics and include a pressure-density relation that incorporates the effects of nonlinearity, power law absorption, and medium heterogeneities. The additional terms accounting for convective nonlinearity and power law absorption are expressed as spatial gradients making them efficient to numerically encode. The governing equations are then discretized using a k-space pseudospectral technique in which the spatial gradients are computed using the Fourier-collocation method. This increases the accuracy of the gradient calculation and thus relaxes the requirement for dense computational grids compared to conventional finite difference methods. The accuracy and utility of the developed model is demonstrated via several numerical experiments, including the 3D simulation of the beam pattern from a clinical ultrasound probe.     

Phase-space moment-equation model of highly relativistic electron-beams in plasma-wakefield accelerators

We formulate a new procedure for modelling the transverse dynamics of relativistic electron beams with significant energy spread when injected into plasma-based accelerators operated in the blow-out regime. Quantities of physical interest, such as the emittance, are furnished directly from solution of phase space moment equations formed from the relativistic Vlasov equation. The moment equations are closed by an Ansatz, and solved analytically for prescribed wakefields. The accuracy of the analytic formulas is established by benchmarking against the results of a semi-analytic/numerical procedure which is described within the scope of this work, and results from a simulation with the 3D quasi-static PIC code HiPACE.