scholarly journals The Use of a Distributed Body Force to Simulate Viscous Effects in 3D Flow Calculations

Author(s):  
J. D. Denton

Three dimensional viscous flow calculations methods for turbomachinery are starting to become available but are not yet sufficiently well developed to be used for design purposes. Three dimensional inviscid calculations on the other hand are now well developed and are widely used for design purposes. This paper describes a method intermediate between fully viscous methods and inviscid methods. The viscous effects are approximated by a very simple model which can be tuned empirically to get the correct overall level of loss and which reproduces many of the details of real viscous flow, such as boundary layers and secondary flows. The method is a simple extension to a widely used inviscid method and enables viscous effects to be simulated with little extra computational cost compared to a 3D inviscid calculation.

Energies ◽  
2021 ◽  
Vol 14 (23) ◽  
pp. 8080
Author(s):  
Iván Castro-Fernández ◽  
Ricardo Borobia-Moreno ◽  
Rauno Cavallaro ◽  
Gonzalo Sánchez-Arriaga

The validity of using a low-computational-cost model for the aerodynamic characterization of Airborne Wind Energy Systems was studied by benchmarking a three-dimensional Unsteady Panel Method (UnPaM) with experimental data from a flight test campaign of a two-line Rigid-Framed Delta kite. The latter, and a subsequent analysis of the experimental data, provided the evolution of the tether tensions, the full kinematic state of the kite (aerodynamic velocity and angular velocity vectors, among others), and its aerodynamic coefficients. The history of the kinematic state was used as input for UnPaM that provided a set of theoretical aerodynamic coefficients. Disparate conclusions were found when comparing the experimental and theoretical aerodynamic coefficients. For a wide range of angles of attack and sideslip angles, the agreement in the lift and lateral force coefficients was good and moderate, respectively, considering UnPaM is a potential flow tool. As expected, UnPaM predicts a much lower drag because it ignores viscous effects. The comparison of the aerodynamic torque coefficients is more delicate due to uncertainties on the experimental data. Besides fully non-stationary simulations, the lift coefficient was also studied with UnPaM by assuming quasi-steady and steady conditions. It was found that for a typical figure-of-eight trajectory there are no significant differences between unsteady and quasi-steady approaches allowing for fast simulations.


2000 ◽  
Author(s):  
L. J. Lenke ◽  
H. Simon

Abstract From the wide range of applications in which centrifugal compressors are used, two different return channels will be investigated, to demonstrate the influence of flow coefficients onto the flow structure. The investigated return channels are typical to join the exit from one stage of a centrifugal machine to the inlet of the next stage and cover the range of very small and large flow coefficients. Starting with a comparison between measurements and numerical results to demonstrate the performance of the calculations, the comparison of both return channels show that three-dimensional phenomena and viscous effects such as secondary flow and wakes have appreciable effect on the fluid dynamics and performance of centrifugal compressors. Especially the deceleration of the flow introduces large separations and recirculations which will decrease the efficiency. Furthermore, a variation of the 180°-bend demonstrates the influence of streamline curvature onto the separation behaviour within the whole return channel. Due to the strongly three-dimensional flow structure with high streamline curvature and secondary flows on hub and shroud of the return channel vanes, a modified explicit algebraic Reynolds stress model will be used for all calculations.


2021 ◽  
Vol 931 ◽  
Author(s):  
Osamu Sano ◽  
Timir Karmakar ◽  
G.P. Raja Sekhar

Viscous flow around spherical macroscopic cavities in a granular material is investigated. The Stokes equation inside and the Darcy–Brinkman equation outside the cavities are considered. In particular, the interaction of two equally sized cavities positioned in tandem is examined in detail, where the asymptotic effect of the other cavity is taken into account. The present analysis gives a reasonable estimate on the volume flow into the cavity and the local enhancement of stresses. This is applicable to predict the microscale waterway formation in that material, onset of landslides, collapse of cliffs and river banks, etc.


Author(s):  
Marek Pátý ◽  
Sergio Lavagnoli

Abstract The efficiency of modern axial turbomachinery is strongly driven by the secondary flows within the vane or blade passages. The secondary flows are characterized by a complex pattern of vortical structures that origin, interact and dissipate along the turbine gas path. The endwall flows are responsible for the generation of a significant part of the overall turbine loss because of the dissipation of secondary kinetic energy and mixing-out of non-uniform momentum flows. The understanding and analysis of secondary flows requires a reliable vortex identification technique to predict and analyse the impact of specific turbine designs on the turbine performance. However, literature shows a remarkable lack of general methods to detect vortices and to determine the location of their cores and to quantify their strength. This paper presents a novel technique for the identification of vortical structures in a general 3D flow field. The method operates on the local flow field and it is based on a triple decomposition of motion proposed by Kolář. In contrast to a decomposition of velocity gradient into the strain and vorticity tensors, this method considers a third, pure shear component. The subtraction of the pure shear tensor from the velocity gradient remedies the inherent flaw of vorticity-based techniques which cannot distinguish between rigid rotation and shear. The triple decomposition of motion serves to obtain a 3D field of residual vorticity whose magnitude is used to define vortex regions. The present method allows to locate automatically the core of each vortex, quantify its strength and determine the vortex bounding surface. The output may be used to visualize the turbine vortical structures for the purpose of interpreting the complex three-dimensional viscous flow field, as well as to highlight any case-to-case variations by quantifying the vortex strength and location. The vortex identification method is applied to a high-pressure turbine with three optimized blade tip geometries. The 3D flow-field is obtained by CFD computations performed with Numeca FINE/Open. The computational model uses steady-state RANS equations closed by the Spalart-Allmaras turbulence model. Although developed for turbomachinery applications, the vortex identification method proposed in this work is of general applicability to any three-dimensional flow-field.


Author(s):  
Zhiheng Li ◽  
Jiawei Yu ◽  
Dakui Feng ◽  
Kaijun Jiang ◽  
Yujie Zhou

Abstract The virtual propeller model can achieve the rapid numerical prediction of the ship self-propulsion performance through viscous flow, which used the improved body-force method. The two-dimensional lift coefficient CL and the drag coefficient CD are very important parameters in this method, which are generally obtained by the potential flow methods and cannot incorporate viscous effects. This study will perform a fully nonlinear unsteady RANS (Reynolds Average Navier-Stokes) simulation to get the KP505 open-water characteristics and then divide its blade into several parts to get the lift coefficient CL and the drag coefficient CD on each one. Then fitting by multivariate regression method, the relationship between CL, CD and propeller parameters is obtained. The Unsteady Blade Element Theory (UBET) is coupled with RANS in house CFD code HUST-Ship (Hydrodynamic Unsteady Simulation Technology for Ship) to calculate the flow around the propeller. RANS equations are solved by the finite difference method and PISO arithmetic. have been made using structured grid with overset technology. The results show that comparing with the EFD data, the maximum differences of the result of the improved body-force method are 4.32% and 2.7% for the thrust coefficient and the torque coefficient respectively near the propeller operating point.


1996 ◽  
Vol 118 (3) ◽  
pp. 492-502 ◽  
Author(s):  
S. Kang ◽  
C. Hirsch

A Navier–Stokes solver is applied to investigate the three-dimensional viscous flow in a low-speed linear compressor cascade with tip clearance at design and off-design conditions with two different meshes. The algebraic turbulence model of Baldwin–Lomax is used for closure. Relative motion between the blades and wall is simulated for one flow coefficient. Comparisons with experimental data, including flow structure, static and total pressures, velocity profiles, secondary flows and vorticity, are presented for the stationary wall case. It is shown that the code predicts well the flow structure observed in experiments and shows the details of the tip leakage flow and the leading edge horseshoe vortex.


2020 ◽  
Vol 31 (08) ◽  
pp. 2050113
Author(s):  
H. Montaseri ◽  
K. Tavakoli ◽  
S. Evangelista ◽  
P. Omidvar

Lateral intakes are hydraulic structures used for domestic, agricultural and industrial water conveyance, characterized by a very complex three-dimensional morphodynamic behavior: since streamlines near the lateral intake are deflected, some vortices form, pressure gradient, shear and centrifugal forces at the intake generate flow separation and a secondary movement, responsible for local scour and sediment deposition. On the other side, the modeling of flows, besides the sediment transport, in curved channels implies some more complications in comparison with straight channels. In this research, this complex process has been investigated experimentally and numerically, with the mechanism of sediment transport, bed topography evolution, flow pattern and their interactions. Experiments were performed in the Laboratory of Tarbiat Modares University, Iran, where a U-shaped channel with a lateral intake was installed and dry sediment was injected at constant rate into a steady flow. Due to the spiral flow, the bed topography changes significantly and the bed forms in turn affect the sediment entering the intake. Different from the previous works on this topic which were mainly based on laboratory experiments, here, Computational Fluid Dynamics (CFD) numerical simulations with FLUENT software were also performed, specifically with the two-phase Eulerian Model (EM) and Discrete Phase Model (DPM), at the aim of evaluating their performance in reproducing the observed physical processes. This software is used for a large variety of CFD problems, but not much for simulating sediment transport phenomena and bed topography evolution. The comparison of the results obtained through the two models against the laboratory experimental data proved a good performance of both the models in reproducing the main features of the flow, for example, the longitudinal and vertical streamlines and the mechanism of particles movement. However, the EM reveals a better performance than DPM in the prediction of the secondary flows and, consequently, of the bed topography evolution, whereas the DPM well depicts the particles pattern, predicts the location of trapped particles and determines the percentage of sediment entering the intake. The numerical models so calibrated and validated were applied to other cases with different positions of the intake in the bend. The results show that mechanism of sediment entrance into the intake varies in different position. If the intake is installed in the second half of the bend, the sediment accumulates along the inner bank of the bend and enters the intake from downstream edge of intake; on the other side, if it is placed in the first half of the bend, the sediment accumulates along both the inner and the outer bends and, therefore, more sediment enters the intake. Also the results of the simulations performed with the DPM model for different positions of the lateral intake show that for all discharge ratios, the position of 120∘ is the one which guarantees the minimum ratio of sediment diverted to the intake (Gr).


2015 ◽  
Vol 119 (1215) ◽  
pp. 613-630 ◽  
Author(s):  
M. Bouriga ◽  
R. Taher ◽  
F. Morency ◽  
J. Weiss

AbstractThe flow inside a constant-width wind-tunnel contraction is simulated by solving the Reynolds-Averaged Navier-Stokes equations with an eddy-viscosity turbulence model. The results show the presence of longitudinal vortices near the sidewalls centreline. This confirms a former hypothesis involving the generation of skew-induced longitudinal vorticity within the sidewalls boundary layers. Detailed analysis reveals that the flow structure is influenced by viscous effects in the boundary layers and streamline curvature in the potential flow. Three-dimensional boundary-layer profiles on the contraction sidewall are analysed in the framework of the streamline co-ordinate system and its associated hodographic diagram. The resulting profiles help understand the generation of secondary flows and the associated longitudinal vorticity.


2020 ◽  
Vol 124 (1281) ◽  
pp. 1683-1701 ◽  
Author(s):  
M. Righi ◽  
L.E. Ferrer-Vidal ◽  
V. Pachidis

ABSTRACTThis paper describes the application of low-order models to the prediction of the steady performance of axial compressors at sub-idle conditions. An Euler body-force method employing sub-idle performance correlations is developed and presented alongside a mean-line approach employing the same set of correlations. The low-order tools are used to generate the characteristic lines of the compressor in the locked-rotor and zero-torque windmilling conditions. The results are compared against steady-state operating points from three-dimensional (3D) Reynolds-averaged Navier–Stokes (RANS) computational fluid dynamics (CFD) simulations. The accuracy of the low-order tools in reproducing the results from high-fidelity CFD is analysed, and the trade-off with the computational cost of each method is discussed. The low-order tools presented are shown to offer a fast alternative to traditional CFD which can be used to predict the performance in sub-idle conditions of a new compressor design during early development stages.


Author(s):  
Ch. Hirsch ◽  
S. Kang ◽  
G. Pointel

The three-dimensional flow in centrifugal impellers is investigated on the basis of a detailed analysis of the results of numerical simulations. An in-depth validation has been performed, based on the computations of Krain’s centrifugal compressor and a radial pump impeller, both with vaneless diffusers and detailed comparisons with available experimental data, discussed in Part I, provide high confidence in the numerical tools and results. The low energy, high loss ‘wake’ region results from a balance between various contributions to the secondary flows influenced by tip leakage flows and is not necessarily connected to 3D boundary layer separation. A quantitative evaluation of the different contributions to the streamwise vorticity is performed, identifying the main features influencing their intensity. The main contributions are: the passage vortices along the end walls due to the flow turning; a passage vortex generated by the Coriolis forces proportional to the local loading and mainly active in the radial parts of the impeller; blade surface vortices due to the meridional curvature. The analysis provides an explanation for the differences in wake position under different geometries and flow conditions. A secondary flow representation is derived from the calculated 3D flow field for the two geometries validated in Part I, and the identified flow features largely confirm the theoretical analysis.


Sign in / Sign up

Export Citation Format

Share Document