Turbulent flow around single concentric long capsule in a pipe

A numerical solution was developed for the equations governing the turbulent flow around single concentric long capsule in a pipe. First, a turbulence model was established for the concentric annulus between the capsule and the pipe to simulate the flow as axi-symmetric, two dimensional, steady flow without edge effect. Second, the same case was considered taking into account the edge effect. Finally, turbulence modelling was established to simulate the case as a three dimensional steady flow, with a view of investigating the validity of axi-symmetric flow assumption. Three different turbulence models were used: an algebraic model (Baldwin–Lomax model) and two types of two-equation models (k–ε and k–ω). Obtained results of pressure gradient along the capsule were compared with available experimental data to verify the used models. In addition, experimental data of the velocity profiles of other investigators were also used in this concern. The results predicted by the three different turbulence models were shown to agree well with the experimental data, though precision differed from one to another.     

Numerical evaluation of the suppression effect of a free-to-rotate triangular fairing on the vortex-induced vibration of a circular cylinder

The suppression of vortex-induced vibration (VIV) of a circular cylinder with a free-to-rotate triangular fairing in the Reynolds number range of Re = 1100–6100 is numerically investigated using computational fluid dynamics. The unsteady Reynolds-averaged Navier–Stokes equations and the shear stress transport k–ω turbulence model coupled with an improved fourth-order Runge–Kutta method are used to solve the wake flow, the structure's vibration, and the fairing's rotation. The computational model is validated with the available experimental results for a cylinder with an attached short-tail fairing. The numerical results indicate that the triangular fairing has a positive role in suppressing vibration when it achieves a stable position deflected from the flow direction. The suppression effect is sensitive to the incoming flow velocity. The fairing shifts from a stable state to an unstable one when the flow velocity varies. Therefore, maintaining the hydrodynamic stability of the fairing is the key to achieving success in vibration suppression, and the stability is dependent on the characteristic length and the rotational friction. Although the strong flapping of the 70° triangular fairing excites a more vigorous vibration, it may be used as an amplifier of VIV for energy harvesting.     

Investigation on the effect of trailing edge ejection on a turbine cascade

This paper presents a detailed experimental and numerical investigation for a turbine cascade with different trailing edge ejection. The numerical simulation is based on Three-Dimensional Navier–Stokes equations coupled with an effective ejection model, where a high resolution non-oscillatory scheme, LU-SGS implicit algorithm and Baldwin-Lomax turbulence model are employed. The experiments presented in this paper focused on a transonic turbine cascade performance with different ejection to validate the numerical simulation results. The results show that the blowing ratio has a small effect on the Mach number distribution and exit flow angle with two slot types. However the energy loss coefficient increases initially, and subsequently has a decrease tendency with the increasing of blowing ratio. The ejection from the symmetry slot blows away the vortex at the blade trailing edge and strengthens the mixing between the wake and main flow. The ejection from the pressure side cutback only clears up the vortex near the slot surface, and has small effect on the flow field near the trailing edge.     

Assessment of turbulence modeling for gas flow in two-dimensional convergent–divergent rocket nozzle

In the present study, the turbulent gas flow dynamics in a two-dimensional convergent–divergent rocket nozzle is numerically predicted and the associated physical phenomena are investigated for various operating conditions. The nozzle is assumed to have impermeable and adiabatic walls with a flow straightener in the upstream side and is connected to a plenum surrounding the nozzle geometry and extended in the downstream direction. In this integrated component model, the inlet flow is assumed a two-dimensional, steady, compressible, turbulent and subsonic. The physics based mathematical model of the considered flow consists of conservation of mass, momentum and energy equations subject to appropriate boundary conditions as defined by the physical problem stated above. The system of the governing equations with turbulent effects is solved numerically using different turbulence models to demonstrate their numerical accuracy in predicting the characteristics of turbulent gas flow in such complex geometry. The performance of the different turbulence models adopted has been assessed by comparing the obtained results of the static wall pressure and the shock position with the available experimental and numerical data. The dimensionless shear stress at the nozzle wall and the separation point are also computed and the flow field is illustrated. The various implemented turbulence models have shown different behavior of the turbulent characteristics. However, the shear-stress transport (SST) k–ω model exhibits the best overall agreement with the experimental measurements. In general, the proposed numerical procedure applied in the present paper shows good capability in predicting the physical phenomena and the flow characteristics encountered in such kinds of complex turbulent flow.     

Three‐dimensional boundary layer on wind turbine blades

This paper describes a three‐dimensional analysis of the laminar boundary layer that develops on the blades of an horizontalaxis wind turbine. The main aim was to investigate a fundamental phenomenon: the effect of rotation on the blade boundary layer of a wind turbine in conjunction with the widely observed phenomenon of stall‐delay. The separation position in retarded flows with pressure gradients is calculated and compared for the rotation and non‐rotation cases. It is concluded that the stall is linearly postponed due to the Coriolis force and the separation point is delayed as a result of increasing rotation speed or decreasing blade spanwise position. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of turbulence models in simulating swirling pipe flows

Swirling flow is a common phenomenon in engineering applications. A numerical study of the swirling flow inside a straight pipe was carried out in the present work with the aid of the commercial CFD code fluent. Two-dimensional simulations were performed, and two turbulence models were used, namely, the RNG k–ε model and the Reynolds stress model. Results at various swirl numbers were obtained and compared with available experimental data to determine if the numerical method is valid when modeling swirling flows. It has been shown that the RNG k–ε model is in better agreement with experimental velocity profiles for low swirl, while the Reynolds stress model becomes more appropriate as the swirl is increased. However, both turbulence models predict an unrealistic decay of the turbulence quantities for the flows considered here, indicating the inadequacy of such models in simulating developing pipe flows with swirl.     

An assessment of turbulence models applied to the simulation of a two-dimensional submerged jet

This paper is concerned with the investigation of the performance of different turbulence models in the numerical prediction of transient flow caused by a confined submerged jet. Several widely used models, i.e., the standard k–ε, RNG k–ε, low Reynolds number k–ε models and the differential Reynolds stress model, as included in CFD codes, were compared with each other for a two-dimensional, incompressible, turbulent jet flow and with reported experimental data. A flapping oscillation was predicted regardless of the model used. A chosen Strouhal (St) number definition brought the fundamental frequencies from both the experiments and computations into close proximity. However, different turbulence models have exhibited quite different behaviours in terms of the frequency and regularity of the oscillation and in terms of the scale and duration of the vortices generated. All versions of the k–ε model yielded regular oscillations, which agree with experimental observations. On the other hand, the Reynolds stress (RS) model produced a complex pattern but a slower dissipation of vortices. In addition, some aspects of gridding and inflow representation are also discussed.     

Berechnung der reibungslosen Strömung für ein vorgegebenes ebenes Schaufelgitter bei hohen Unterschallgeschwindigkeiten

Summary A simple approximation method based onPrandtl-Glauert's rule is given for the calculation of inviscous compressible flow through a two-dimensional cascade. It is applicable to cascades of any solidity and stagger angle with blade sections of small thickness and camber. The calculation of compressible flow for a predetermined cascade at a given Mach number is reduced to the calculation of an ‘associated incompressible flow’ through a cascade having blades of the same blade section, a higher solidity and larger angle of stagger. This associated incompressible cascade flow can be conveniently calculated on the basis of a method disclosed byH. Schlichting [2]. Comparison of some computed pressure distributions along the blade contour reveals satisfactory agreement with measurements in the high-speed cascade wind tunnel of the Deutsche Forschungsanstalt für Luftfahrt (DFL; German Research Centre for Aviation), Brunswick.

---

---

Eine ausführliche Ver?ffentlichung dieser Untersuchungen ist in ?Forschung auf dem Gebiete des Ingenieurwesens?24, 19–28, (1958), erschienen.     

Numerical modelling of subcritical open channel flow using the K-ε turbulence model and the penalty function finite element technique

A numerical model has been developed that employs the penalty function finite element technique to solve the vertically averaged hydrodynamic and turbulence model equations for a water body using isoparametric elements. The full elliptic forms of the equations are solved, thereby allowing recirculating flows to be calculated. Alternative momentum dispersion and turbulence closure models are proposed and evaluated by comparing model predictions with experimental data for strongly curved subcritical open channel flow. The results of these simulations indicate that the depth-averaged two-equation k-ε turbulence model yields excellent agreement with experimental observations. In addition, it appears that neither the streamline curvature modification of the depth-averaged k-ε model, nor the momentum dispersion models based on the assumption of helicoidal flow in a curved channel, yield significant improvement in the present model predictions. Overall model predictions are found to be as good as those of a more complex and restricted three-dimensional model.     

Numerical study on the fluid-structure interaction in a model aquatic canopy flow

This work presents numerical simulations and selected results of the flow over aquatic canopies, consisting of artificial flexible rectangular blades, arranged in a well-defined order. The results obtained with three different Reynolds and Cauchy numbers are compared with experimental data achieving good agreement. The considered range of Cauchy numbers represents three different types of canopies ranging from rigid up to highly flexible plants. The transient flow data and blade positions are statistically analyzed to gain deeper understanding of the complex physical processes for this kind of fluid structure interaction. For example, the correlation of role of large scale motion of the flexible blades in conjunction with coherent vortex structures of the flow is addressed. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

RBF Neural Network Based Prediction on Blade Surface Pressure Fields in Compressors北大核心CSCD

The airflow characteristics of the internal flow path of an aero-engine compressor are complex, and the vortex flow field around the blade is characterized by high pressure, high speed, rotation, and unsteadiness. Therefore, there is an urgent need to calculate and predict the aerodynamic characteristics of the complex flow field around the compressor blade efficiently and accurately. The computational fluid dynamics (CFD) method was used to generate the aerodynamic load distribution on the blade surface under different operating conditions for the study of the complex flow fields around aero-engine blades. The radial based function (RBF) neural network was applied to establish the pressure surface aerodynamic load prediction model, and the neural network modeling method was combined with the flow field calculation. The neural network method can learn and train the CFD-based data set to properly compensate the errors from the CFD, which provides a reference for the effective prediction of the complex flow fields around aero-engine compressor blades. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Numerical modeling of turbulence-induced interfacial instability in two-phase flow with moving interface

The present paper introduces a new interfacial marker-level set method (IMLS) which is coupled with the Reynolds averaged Navier–Stokes (RANS) equations to predict the turbulence-induced interfacial instability of two-phase flow with moving interface. The governing RANS equations for time-dependent, axisymmetric and incompressible two-phase flow are described in both phases and solved separately using the control volume approach on structured cell-centered collocated grids. The transition from one phase to another is performed through a consistent balance of kinematic and dynamic conditions on the interface separating the two phases. The topological changes of the interface are predicted by applying the level set approach. By fitting a number of interfacial markers on the intersection points of the computational grids with the interface, the interfacial stresses and consequently, the interfacial driving forces are easily estimated. Moreover, the normal interface velocity, calculated at the interfacial markers positions, can be extended to the higher dimensional level set function and used for the interface advection process. The performance of linear and non-linear two-equation k–ε turbulence models is investigated in the context of the considered two-phase flow impinging problem, where a turbulent gas jet impinging on a free liquid surface. The numerical results obtained are evaluated through the comparison with the available experimental and analytical data. The nonlinear turbulence model showed superiority in predicting the interface deformation resulting from turbulent normal stresses. However, both linear and nonlinear turbulence models showed a similar behavior in predicting the interface deformation due to turbulent tangential stresses. In general, the developed IMLS numerical method showed a remarkable capability in predicting the dynamics of the considered two-phase immiscible flow problems and therefore it can be applied to quite a number of interface stability problems.     

Aus- und Eintrittsstösse an Schaufelgittern

Summary Though the flow in the cylindrical surface of turbines or compressors may be treated as flow in a plane by development of the cylinder into a plane, the double connection once around the cylinder introduces difficulties in the law of forbidden signals for supersonic flow with a subsonic axial component. In cases of large gaps between the stages and blades of zero thickness and straight exit and entry zones, the methods of steady one-dimensional flow can shed light on the adaptability of the flow to the transient conditions by means of special kinds of shock polars for the exit and entry flow at given blade anglesβ (Figures 7 and 8). The shocks may sometimes be a combination ofHugoniot's discontinuities andCarnot's transient zones, but are in all cases governed by mass, impulse and energy balances with an increase in entropy.      

Numerical simulation of particulate-flow in spiral separators: Part I. Low solids concentration (0.3% & 3% solids)

The aim of the present study is the simulation of the particulate flow in spiral separators. The study is based on Eulerian approach and turbulence modeling. The results focus on particulate-flow characteristics such as the velocity, the distribution, and concentration of particulates on the spiral trough. The predicted results are compared with the experimental findings from LD9 coal spiral. The comparison shows good agreement and indicates that the most accurate turbulence model is RNG K–ε.     

Numerical simulation of UV disinfection reactors: Evaluation of alternative turbulence models

Six turbulence models, including standard k–ε, k–ε RNG, k–ω (88), revised k–ω (98), Reynolds stress transport model (RSTM), and two-fluid model (TFM), were applied to the simulation of a closed conduit polychromatic UV reactor. Predicted flow field and turbulent kinetic energy were compared with the experimental data from a digital particle image velocimetry (DPIV). All of the predicted flow fields were combined with a multiple segment source summation (MSSS) fluence rate model and three different microbial response kinetic models to simulate the disinfection process at two UV lamp power conditions. Microbial transport was simulated using the Lagrangian particle tracking method. The results show that the fluence distributions and the effluent inactivation levels were sensitive to the turbulence model selection. The level of sensitivity was a function of the operating conditions and the UV response kinetics of the microorganisms. Simulations with operating conditions that produced higher log inactivation or utilized microorganisms with higher UV sensitivity showed greater sensitivity to the turbulence model selection. In addition, a broader fluence distribution was found with turbulence models that predicted a larger wake region behind the lamps.     

Influence of turbulence on the performance of a laboratory scale HAWT

The oncoming wind to horizontal axis wind turbines (HAWT) may change its speed and direction stochastically in time. Hence, turbine blades are exposed to flows both with fluctuating angle of attack and fluctuating yaw angles. The modern wind turbines are reacting to those changes by pitch angle and torque control not only to exploit as much power as possible but also stabilize energy production and prevent any damage of the turbine. However, time scales of wind fluctuations and sudden changes of wind properties can be very short and with very high in amplitude. In the present study we focus on the influence of turbulence on the performance of a HAWT. Main motivation of the investigations is to figure out best strategies for the aerodynamic design of the blades operating under turbulent conditions. A laboratory scale HAWT and a performance measurement set-up are employed to measure the influence of the oncoming wind. The tests are conducted in the closed loop wind tunnel of our institute. The test section of the tunnel is 1.87 m in width, 1.4 m in height and 2 m in length. The rotor blades are specially designed and optimized for this wind tunnel and the generator used. The turbulence is generated by two static squared mesh grids; fine and coarse one. Hence, two mainly different turbulence scales are obtained. In addition, the distance between the wind-turbine and the grid is adjusted to have additional sub-turbulence scales for each grid. The turbulence is nearly isotropic and decays in the flow direction. The developments of Taylor's micro scale (λ_g) and integral scale (L_g) of the turbulence in the flow direction at various incoming wind velocities (8−16 m/s) are measured. Hence, the facility allows to expose the wind-turbine to turbulence with various energy and length scale content. Those measurements are conducted with hot-wire anemometry in the absence of the wind-turbine. Upstream and downstream turbulence intensities (TI) distributions are measured to give insight on the surrounding free stream and turbine wake interaction and how can different turbulence eddies scales contribute in the influence of the performance of the turbine. Performance measurements are conducted with and without turbulence and the results are compared. The study shows that the higher the turbulence, the more the power extracted by the turbine. This is due to the higher interaction of large eddies with the turbine wake and with the boundary layer, which helps to keeping it attached. Furthermore, higher TI's help in suppressing the tip vortex, thus, reduce turbine tip losses. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical study of the effect of changing the cyclone cone length on the gas flow field

Models of six different cyclone separator designs have been investigated numerically to study the effect of changing cone length on the gas flow field. The cone length was changed in two ways, one by fixing the cone diameter and the other by fixing the cone angle. The Reynolds stress model was applied to predict the turbulence of gas flow. The model was validated by comparing the numerical results to published experimental data. The results show that, for the cyclone without cylindrical part, the pressure drop and the value of tangential and axial velocity components increase with increasing cone length with fixed cone angle until they reach maximum value. It was found that increasing the cone length with the same cone diameter increases the back flow and the lip leakage below the vortex finder. The back flow and the lip leakage phenomena was observed to be reduced by increasing cone length with fixed cone angle. The effect of contraction ratio D_x/D, which is defined as the ratio between the vortex finder diameter and the cone diameter, on the turbulent kinetic energy and the velocity components also was investigated. As the contraction ratio decreases, the maximum value of the turbulent kinetic energy and the tangential and axial velocity components increases. The pressure drop was fitted into a correlation as a function of dimensionless quantities (Re, D_x/D, and l_c/D) for each design, so that the designers could select the appropriate dimensions according to their design limitations.     

A Lagrangian simulation of particle deposition in a turbulent boundary layer in the presence of thermophoresis

Knowledge of particle deposition in turbulent flows is often required in engineering situations. Examples include fouling of turbine blades, plate-out in nuclear reactors and soot deposition. Thus it is important for numerical simulations to be able to predict particle deposition. Particle deposition is often principally determined by the forces acting on the particles in the boundary layer. The particle tracking facility in the CFD code uses the eddy lifetime model to simulate turbulent particle dispersion, no specific boundary layer being modelled. The particle tracking code has been modified to include a boundary layer. The non-dimensional yplus, y⁺, distance of the particle from the wall is determined and then values for the fluid velocity, fluctuating fluid velocity and eddy lifetime appropriate for a turbulent boundary layer used. Predictions including the boundary layer have been compared against experimental data for particle deposition in turbulent pipe flow. The results giving much better agreement. Many engineering problems also involve heat transfer and hence temperature gradients. Thermophoresis is a phenomena by which small particles experience a force in the opposite direction to the temperature gradient. Thus particles will tend to deposit on cold walls and be repulsed by hot walls. The effect of thermophoresis on the deposition of particles can be significant. The modifications of the particle tracking facility have been extended to include the effect of thermophoresis. A preliminary test case involving the deposition of particles in a heated pipe has been simulated. Comparison with experimental data from an extensive experimental programme undertaken at ISPRA, known as STORM (Simplified Tests on Resuspension Mechanisms), has been made.     

Direct numerical simulation of turbulent non-Newtonian flow using OpenFOAM

Understanding transition and turbulence in the flow of shear-thinning non-Newtonian fluids remains substantially unresolved and additional research is required to develop better computational methods for wall-bounded turbulent flows of these fluids. Previous DNS studies of shear-thinning fluids mainly use purpose-built codes and simple geometries such as pipes and channels. However in practical application, the geometry of mixing vessels, pumps and other process equipment is far more complex, and more flexible computational methods are required. In this paper a general-purpose DNS approach for shear-thinning fluids is undertaken using the OpenFOAM CFD library. DNS of turbulent Newtonian and non-Newtonian flow in a pipe flow are conducted and the accuracy and efficiency of OpenFOAM are assessed against a validated high-order spectral element-Fourier DNS code – Semtex. The results show that OpenFOAM predicts the flow of shear-thinning fluids to be a little more transitional than the predictions from Semtex, with lower radial and azimuthal turbulence intensities and higher axial intensity. Despite this, the first and second order turbulence statistics differ by at most 16%, and usually much less. An assessment of the parallel scaling of OpenFOAM indicates that OpenFOAM scales very well for the CPUs from 8 to 512, but the intranode scalability is poor for less than 8CPUs. The present work shows that OpenFOAM can be used for DNS of shear-thinning fluids in the simple case of pipe flow, and suggests that more complex flows, where flow separation is often important, are likely to be simulated with accuracies that are acceptably good for engineering application.