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.     

Numerical investigation of the variable nozzle effect on the mixed flow turbine performance characteristics

There are various matching ways between turbocharger and engine, the variable nozzle turbine is the most significant method. The turbine design must be economic with high efficiency and large capacity over a wide range of operational conditions. These design intents are used in order to decrease thermal load and improve thermal efficiency of the engine. This paper presents an original design method of a variable nozzle vane for mixed flow turbines developed from previous experimental and numerical studies. The new device is evaluated with a numerical simulation over a wide range of rotational speeds, pressure ratios, and different vane angles. The compressible turbulent steady flow is solved using the ANSYS CFX software. The numerical results agree well with experimental data in the nozzleless configuration. In the variable nozzle case, the results show that the turbine performance characteristics are well accepted in different open positions and improved significantly in low speed regime and at low pressure ratio.     

Numerical simulation of pulsation processes in hydraulic turbine based on 3D model of cavitating flow

A new approach was proposed for simulation of unsteady cavitating flow in the flow passage of a hydraulic power plant. 1D hydro-acoustics equations are solved in the penstock domain. 3D equations of turbulent flow of isothermal compressible liquid-vapor mixture are solved in the turbine domain. Cavitation is described by a transfer equation for liquid phase with a source term which is responsible for evaporation and condensation. The developed method was applied for simulation of pulsations in pressure, discharge, and total energy propagating along the flow conduit of the hydraulic power plant. Simulation results are in qualitative and quantitative agreement with experiment. The influence of key physical and numerical parameters like discharge, cavitation number, penstock length, time step, and vapor density on simulation results was studied.     

Experimental investigation of the magnetohydrodynamic parachute effect in a hypersonic air flow

New data on experimental implementation of the magnetohydrodynamic (MHD) parachute configuration in an air flow with Mach number M = 6 about a flat plate are considered. It is shown that MHD interaction near a flat plate may transform an attached oblique shock wave into a normal detached one, which considerably extends the area of body-incoming flow interaction. This effect can be employed in optimizing return space vehicle deceleration conditions in the upper atmosphere.     

Numerical simulation of cavitating flow around a slender body with slip boundary condition

In this paper, we perform a numerical simulation of the cavitating flow around an underwater hemispherical-head slender body running at a high speed. For the first time, the slip boundary condition is introduced into this problem, and we find that the slip boundary condition has a big influence on the cavitation in the flow-separation zone. By simulating the cavitating flow under different cavitation numbers, we demonstrate that the slip boundary condition can effectively reduce the intensity of cavitation, as represented by the length of cavitation bubbles. The present paper provides a new method for utilization of new surface materials to control the cavitation on the underwater moving objects.     

Experimental study of pressure pulsations in the flow duct of a medium-size model hydroelectric generator with Francis turbine

In the present study, we report on the results of an experimental study of pressure pulsations in the flow duct of a medium-scale hydrodynamic bench with Francis turbine. In various regimes, integral and pulsation characteristics of the turbine were measured. With the help of high-speed filming, the structure of the flow behind the turbine runner was analyzed, and the influence of this structure on the intensity and frequency of pressure pulsations in the flow duct was demonstrated.     

Numerical investigation of separation-induced transition in a lowpressure turbine cascade in a low-disturbance environment

In this paper,the separation-induced transition in an LPT(low-pressure turbine)cascade is investigated at low Reynolds number with DNS(direct numerical simulation).The transition process is accurately predicted giving good agreements between the DNS and experimental results.To illustrate the secondary instability of separation-induced transition in a low-disturbance environment,the results are comprehensively analyzed in both Fourier space and physical space.It is illustrated that the effect of hyperbolic instability dominates around the saddle point of hyperbolic streamlines.This instability mechanism is responsible for the emergence of the streamwise vortices in the braid region.Elongated and intensified because of the“stretching”effect of the background flow,these vortices become the most noticeable characteristic of the flow field.Fundamental modes of small spanwise wavelength are excited in the braid region,so as some low-frequency modes.The elliptical instability plays a minor role than hyperbolic instability.It is also observed that the fundamental mode with a larger spanwise wavelength is unstable in the vortex core which is associated with the deformation of the vortex core via elliptical instability.There is no convincing evidence for the existence of subharmonic instability.     

Numerical investigation of the flow in the cathode melt droplet of electric arc

The role of electromagnetic forces and the forces of viscous friction with the arc plasma flow in the flow formation within the cathode melt is considered within the framework of numerical modelling; a comparative estimation of the separate influence of each of the above forces is carried out. The melt flow pattern is found to form mainly by electromagnetic forces. The character of the effect of electromagnetic forces is determined to a considerable extent by the ratio of the radius of cathode attachment of the arc on the melt droplet and the rod cathode radius.     

Numerical simulation of transitional flow on a wind turbine airfoil with RANS-based transition model

This paper presents a numerical investigation of transitional flow on the wind turbine airfoil DU91-W2-250 with chord-based Reynolds number Re_c = 1.0 × 10⁶. The Reynolds-averaged Navier–Stokes based transition model using laminar kinetic energy concept, namely the k ? k_L ? ω model, is employed to resolve the boundary layer transition. Some ambiguities for this model are discussed and it is further implemented into OpenFOAM-2.1.1. The k ? k_L ? ω model is first validated through the chosen wind turbine airfoil at the angle of attack (AoA) of 6.24° against wind tunnel measurement, where lift and drag coefficients, surface pressure distribution and transition location are compared. In order to reveal the transitional flow on the airfoil, the mean boundary layer profiles in three zones, namely the laminar, transitional and fully turbulent regimes, are investigated. Observation of flow at the transition location identifies the laminar separation bubble. The AoA effect on boundary layer transition over wind turbine airfoil is also studied. Increasing the AoA from ?3° to 10°, the laminar separation bubble moves upstream and reduces in size, which is in close agreement with wind tunnel measurement.     

Numerical investigation of supersonic flow breakdown at the inlet duct throttling

The work presents the results of investigating the process of supersonic flow deceleration in a duct of the two-dimensional inlet throttled by variation of the outlet cross-sectional area. An inlet with three external compression shock waves designed for the freestream Mach number M_d = 7 was considered as an example for the investigation. A one-dimensional analysis of the conditions for realization of the supersonic flow deceleration regimes in the inlet duct with two throats — in the inlet entrance and at the inlet duct outlet, has been carried out. The parametric numerical computations of two-dimensional inviscid or turbulent flows in the inlet were performed with the use of the Euler and Navier—Stokes codes of the program package FLUENT. The critical conditions for the nonuniform flow in the outlet throat bringing to choking the inlet duct were determined.     

Numerical investigation of air free convection in a room with heat source

Mathematical modelling of three-dimensional steady free convective incompressible viscous gas flows in the rooms with a heat source is carried out within the framework of the Navier — Stokes equations with the effective viscosity determined on the basis of the κ-ɛ turbulence model. The investigation is carried out for the case of model rooms and heat sources having the form of rectangular parallelepipeds with square bases. The influence of the heat source power and the sizes of the room base on local and averaged values of the air velocity and temperature in the rooms is analysed. The flow pattern in the room is shown to have a torus-like shape. It is found that the variation of sizes of the room base rather than the capacity variation of the heat source is of determining importance for the gas motion character in a closed volume.     

�������ЧӦ�Ե�������Ӱ�����ֵ�о�

在考虑到捕获电子效应的情况下,对求解二维Fokker-Planck方程的编程进行了反弹平均的修改,使用了交替方向隐式法来求解方程。分析和计算了在不同扩散系数和不同共振区间的情况下,捕获电子效应对驱动电流的影响。结果显示：随着逆纵横比的增加驱动电流密度有明显的下降,在磁轴附近捕获电子效应对电流驱动影响很小;提高波功率并不能很好的改善捕获电子效应对电流驱动的影响;右移共振区间提高共振电子的速度,也不能很好的改善捕获电子效应对电流驱动的影响。所得结果与理论分析基本一致。     

捕获电子效应对电流驱动影响的数值研究

在考虑到捕获电子效应的情况下,对求解二维Fokker-Planck方程的编程进行了反弹平均的修改,使用了交替方向隐式法来求解方程。分析和计算了在不同扩散系数和不同共振区间的情况下,捕获电子效应对驱动电流的影响。结果显示：随着逆纵横比的增加驱动电流密度有明显的下降,在磁轴附近捕获电子效应对电流驱动影响很小;提高波功率并不能很好的改善捕获电子效应对电流驱动的影响;右移共振区间提高共振电子的速度,也不能很好的改善捕获电子效应对电流驱动的影响。所得结果与理论分析基本一致。     

Numerical investigation of the effect of insoluble surfactant on drop formation in microfluidic device

Surfactants are commonly used in droplet-based microfluidics to stabilize the droplet interface. In this study, we investigate the effect of insoluble surfactant on drop formation in a capillary microfluidic device. We use a diffuse-interface method to describe the evolution of interface involving insoluble surfactant. The Navier-Stokes/Cahn-Hilliard equations and the surfactant conservation equation are solved by a finite element method along with a grid deformation method. As the surfactant has a non-uniform distribution during the drop formation in general, the surface tension has a gradient on the interface, which affects the flow field and interface evolution. The surfactant effect is discussed for dripping and jetting regimes.     

超声速流场条件下激光辐照耦合效应数值模拟

激光辐照结构物包含复杂的多物理场耦合问题,其存在流、热、固多种机制的耦合效应。结合计算流体力学(CFD)和有限元方法,对超声速条件下的激光辐照平板问题进行了热流固耦合分析。采用CFD方法得到平板附近流场分布,利用有限元方法计算平板的温度分布,并将二者结合起来实现流体和固体间的数据交互。理论分析确定了流场效应的最主要影响参数为来流马赫数与攻角。对于不同马赫数,激光区域在6 Ma条件下存在温度的谷值,小于等于6 Ma条件下主要体现为冷却效应,而6 Ma以上主要体现为气动加热效应。攻角增大会导致激光区流体质量流量的增加,使冷却效应更加明显。最后综合分析了流场气动加热和冷却两种效应的产生机制。     

超声速流场条件下激光辐照耦合效应数值模拟

激光辐照结构物包含复杂的多物理场耦合问题,其存在流、热、固多种机制的耦合效应。结合计算流体力学（CFD）和有限元方法,对超声速条件下的激光辐照平板问题进行了热流固耦合分析。采用CFD方法得到平板附近流场分布,利用有限元方法计算平板的温度分布,并将二者结合起来实现流体和固体间的数据交互。理论分析确定了流场效应的最主要影响参数为来流马赫数与攻角。对于不同马赫数,激光区域在6Ma条件下存在温度的谷值,小于等于6Ma条件下主要体现为冷却效应,而6Ma以上主要体现为气动加热效应。攻角增大会导致激光区流体质量流量的增加,使冷却效应更加明显。最后综合分析了流场气动加热和冷却两种效应的产生机制。     

Experimental investigation of the effect of small-obstacle-induced vortex sheet on the separated flow in cavity

In the present paper, we report results of an experimental study of the influence which a vortex-generating element installed upstream of the main obstacle has on the separated flow and heat transfer in a cross-flow cavitytrench. The element was a small cross-flow rib whose height was an order of magnitude smaller than the depth of the cavity. In the experiments, the variable parameters were the angle of inclination of the frontal and rear walls of the cavity, the rib height, and the rib-to-cavity distance. It is shown that the introduction of additional vortical perturbations into the recirculation zone leads to a substantial modification of both the vortex production process and the distributions of pressure and heat-transfer coefficients. Optimal height of the mini-turbulizer and its optimal location are defined by the fall of the re-attachment point of mini-rib-generated flow onto the rear wall of cavity. In the latter situation, the maximal value of the heat-transfer coefficient increases as compared to the case with no vortex generator used, the increase amounting to 30 %.