Abstract: The paper presents the results of a detailed assessment of several modern Reynolds Averaged Navier-Stokes (RANS) turbulence models for prediction of C3X vane film cooling at various injection regimes. Three models are considered, namely the Shear Stress Transport (SST) model, the modification of the SST model accounting for the streamlines curvature (SST-CC), and the Explicit Algebraic Reynolds Stress Model (EARSM). It is shown that all the considered models face with a problem in prediction of the adiabatic effectiveness in the vicinity of the cooling holes; however, accounting for the Reynolds stress anisotropy within the EARSM model noticeably increases the solution accuracy. On the other hand, further downstream all the models provide a reasonable agreement with the experimental data for the adiabatic effectiveness and among the considered models the most accurate results are obtained with the use EARMS.
Abstract: This document provides numerical and experimental optimization of the aerodynamic performance of a drone equipped with three types of horizontal stabilizer. To build this optimal configuration, an experimental and numerical study was conducted on three parameters: the geometry of the stabilizer (horizontal form or reverse V form), the position of the horizontal stabilizer (up or down), and the landing gear position (closed or open). The results show that up-stabilizer position with respect to the horizontal plane of the fuselage provides better aerodynamic performance, and that the landing gear increases the lift in the zone of stability, that is to say where the flow is not separated.
Abstract: The wide range of industrial applications involved with boiling flows promotes the necessity of establishing fundamental knowledge in boiling flow phenomena. For this purpose, a number of experimental and numerical researches have been performed to elucidate the underlying physics of this flow. In this paper, the improved wall boiling models, implemented on ANSYS CFX 14.5, were introduced to study subcooled boiling flow at elevated pressure. At the heated wall boundary, the Fractal model, Force balance approach and Mechanistic frequency model are given for predicting the nucleation site density, bubble departure diameter, and bubble departure frequency. The presented wall heat flux partitioning closures were modified to consider the influence of bubble sliding along the wall before the lift-off, which usually happens in the flow boiling. The simulation was performed based on the Two-fluid model, where the standard k-ω SST model was selected for turbulence modelling. Existing experimental data at around 5 bars were chosen to evaluate the accuracy of the presented mechanistic approach. The void fraction and Interfacial Area Concentration (IAC) are in good agreement with the experimental data. However, the predicted bubble velocity and Sauter Mean Diameter (SMD) are over-predicted. This over-prediction may be caused by consideration of only dispersed and spherical bubbles in the simulations. In the future work, the important physical mechanisms of bubbles, such as merging and shrinking during sliding on the heated wall will be incorporated into this mechanistic model to enhance its capability for a wider range of flow prediction.
Abstract: A Computational Fluid Dynamics (CFD) technique has been applied to simulate the gas-liquid phase in double stirred tank of Rushton impeller. Eulerian-Eulerian model was adopted to simulate the multiphase with standard correlation of Schiller and Naumann for drag co-efficient. The turbulence was modeled by using standard k-ε turbulence model. The present CFD model predicts flow pattern, local gas hold-up, and local specific area. It also predicts local kLa (mass transfer rate) for single impeller. The predicted results were compared with experimental and CFD results of published literature. The predicted results are slightly over predicted with the experimental results; however, it is in reasonable agreement with other simulated results of published literature.
Abstract: Numerical simulations of vortex-induced vibration of a three-dimensional flexible tube under uniform turbulent flow are calculated when Reynolds number is 1.35×104. In order to achieve the vortex-induced vibration, the three-dimensional unsteady, viscous, incompressible Navier-Stokes equation and LES turbulence model are solved with the finite volume approach, the tube is discretized according to the finite element theory, and its dynamic equilibrium equations are solved by the Newmark method. The fluid-tube interaction is realized by utilizing the diffusion-based smooth dynamic mesh method. Considering the vortex-induced vibration system, the variety trends of lift coefficient, drag coefficient, displacement, vertex shedding frequency, phase difference angle of tube are analyzed under different frequency ratios. The nonlinear phenomena of locked-in, phase-switch are captured successfully. Meanwhile, the limit cycle and bifurcation of lift coefficient and displacement are analyzed by using trajectory, phase portrait, and Poincaré sections. The results reveal that: when drag coefficient reaches its minimum value, the transverse amplitude reaches its maximum, and the “lock-in” begins simultaneously. In the range of lock-in, amplitude decreases gradually with increasing of frequency ratio. When lift coefficient reaches its minimum value, the phase difference undergoes a suddenly change from the “out-of-phase” to the “in-phase” mode.
Abstract: In recent years, as global warming, the sea-ice extent of North Arctic undergoes an evident decrease and Arctic channel has attracted the attention of shipping industry. Ice crystals existing in the seawater of Arctic channel which enter the seawater system of the ship with the seawater were found blocking the seawater pipe. The appearance of cooler paralysis, auxiliary machine error and even ship power system paralysis may be happened if seriously. In order to reduce the effect of high temperature in auxiliary equipment, seawater system will use external ice-water to participate in the cooling cycle and achieve the state of its flow. The distribution of ice crystals in seawater pipe can be achieved. As the ice slurry system is solid liquid two-phase system, the flow process of ice-water mixture is very complex and diverse. In this paper, the flow process in seawater pipe of ice slurry is simulated with fluid dynamics simulation software based on k-ε turbulence model. As the ice packing fraction is a key factor effecting the distribution of ice crystals, the influence of ice packing fraction on the flowing process of ice slurry is analyzed. In this work, the simulation results show that as the ice packing fraction is relatively large, the distribution of ice crystals is uneven in the flowing process of the seawater which has such disadvantage as increase the possibility of blocking, that will provide scientific forecasting methods for the forming of ice block in seawater piping system. It has important significance for the reliability of the operating of polar ships in the future.
Abstract: A recently developed one-equation turbulence model
has been successfully applied to simulate turbulent flows with
various complexities. The model, which is based on the
transformation of the k-ε closure, is wall-distance free and equipped
with lagging destruction/dissipation terms. Test cases included shockboundary-
layer interaction flows over the NACA 0012 airfoil, an
axisymmetric bump, and the ONERA M6 wing. The capability of the
model to operate in a Scale Resolved Simulation (SRS) mode is
demonstrated through the simulation of a massive flow separation
over a circular cylinder at Re= 1.2 x106. An assessment of the results
against available experiments Menter (k-ε)1Eq and the Spalart-
Allmaras model that belongs to the single equation closure family is
made.
Abstract: There is a gap at combustor-turbine interface where leakage flow comes out to prevent hot gas ingestion into the gas turbine nozzle platform. The leakage flow protects the nozzle endwall surface from the hot gas coming from combustor exit. For controlling flow’s stream, the gap’s geometry is transformed by changing fillet radius size. During the operation, step configuration is occurred that was unintended between combustor-turbine platform interface caused by thermal expansion or mismatched assembly. In this study, CFD simulations were performed to investigate the effect of the fillet and step on heat transfer and film cooling effectiveness on the nozzle platform. The Reynolds-averaged Navier-stokes equation was solved with turbulence model, SST k-omega. With the fillet configuration, predicted film cooling effectiveness results indicated that fillet radius size influences to enhance film cooling effectiveness. Predicted film cooling effectiveness results at forward facing step configuration indicated that step height influences to enhance film cooling effectiveness. We suggested that designer change a combustor-turbine interface configuration which was varied by fillet radius size near endwall gap when there was a step at combustor-turbine interface. Gap shape was modified by increasing fillet radius size near nozzle endwall. Also, fillet radius and step height were interacted with the film cooling effectiveness and heat transfer on endwall surface.
Abstract: It is well known that secondary flow loses account about one third of the total loss in any axial turbine. Modern gas turbine height is smaller and have longer chord length, which might lead to increase in secondary flow. In order to improve the efficiency of the turbine, it is important to understand the behavior of secondary flow and device mechanisms to curtail these losses. The objective of the present work is to understand the effect of a stream wise end-wall fence on the aerodynamics of a linear turbine cascade. The study is carried out computationally by using commercial software ANSYS CFX. The effect of end-wall on the flow field are calculated based on RANS simulation by using SST transition turbulence model. Durham cascade which is similar to high-pressure axial flow turbine for simulation is used. The aim of fencing in blade passage is to get the maximum benefit from flow deviation and destroying the passage vortex in terms of loss reduction. It is observed that, for the present analysis, fence in the blade passage helps reducing the strength of horseshoe vortex and is capable of restraining the flow along the blade passage. Fence in the blade passage helps in reducing the under turning by 70 in comparison with base case. Fence on end-wall is effective in preventing the movement of pressure side leg of horseshoe vortex and helps in breaking the passage vortex. Computations are carried for different fence height whose curvature is different from the blade camber. The optimum fence geometry and location reduces the loss coefficient by 15.6% in comparison with base case.
Abstract: In this paper, results concerning flame propagation of various fuels in a particular combustion chamber with four tilted valves were elucidated. Flame propagation was represented by the evolution of spatial distribution of temperature in various cut-planes within combustion chamber while the flame front location was determined by dint of zones with maximum temperature gradient. The results presented are only a small part of broader on-going scrutinizing activity in the field of multidimensional modeling of reactive flows in combustion chambers with complicated geometries encompassing various models of turbulence, different fuels and combustion models. In the case of turbulence two different models were applied i.e. standard k-ε model of turbulence and k-ξ-f model of turbulence. In this paper flame propagation results were analyzed and presented for two different hydrocarbon fuels, such as CH4 and C8H18. In the case of combustion all differences ensuing from different turbulence models, obvious for non-reactive flows are annihilated entirely. Namely the interplay between fluid flow pattern and flame propagation is invariant as regards turbulence models and fuels applied. Namely the interplay between fluid flow pattern and flame propagation is entirely invariant as regards fuel variation indicating that the flame propagation through unburned mixture of CH4 and C8H18 fuels is not chemically controlled.
Abstract: Automotive designers have been trying to use dimples to reduce drag in vehicles. In this work, a car model has been applied with dimple surface with a parameter called dimple ratio DR, the ratio between the depths of the half dimple over the print diameter of the dimple, has been introduced and numerically simulated via k-ε turbulence model to study the aerodynamics performance with the increasing depth of the dimples The Ahmed body car model with 25 degree slant angle is simulated with the DR of 0.05, 0.2, 0.3 0.4 and 0.5 at Reynolds number of 176387 based on the frontal area of the car model. The geometry of dimple changes the kinematics and dynamics of flow. Complex interaction between the turbulent fluctuating flow and the mean flow escalates the turbulence quantities. The maximum level of turbulent kinetic energy occurs at DR = 0.4. It can be concluded that the dimples have generated extra turbulence energy at the surface and as a result, the application of dimples manages to reduce the drag coefficient of the car model compared to the model with smooth surface.
Abstract: A computational study on bio-inspired NACA634-021 hydrofoils with leading-edge protuberances has been carried out to investigate their hydrodynamic flow control characteristics at a Reynolds number of 14,000 and different angles-of-attack. The numerical simulations were performed using ANSYS FLUENT and based on Reynolds-Averaged Navier-Stokes (RANS) solver mode incorporated with k-ω Shear Stress Transport (SST) turbulence model. The results obtained indicate varying flow phenomenon along the peaks and troughs over the span of the hydrofoils. Compared to the baseline hydrofoil with no leading-edge protuberances, the leading-edge modified hydrofoils tend to reduce flow separation extents along the peak regions. In contrast, there are increased flow separations in the trough regions of the hydrofoil with leading-edge protuberances. Interestingly, it was observed that dissimilar flow separation behaviour is produced along different peak- or trough-planes along the hydrofoil span, even though the troughs or peaks are physically similar at each interval for a particular hydrofoil. Significant interactions between adjacent flow structures produced by the leading-edge protuberances have also been observed. These flow interactions are believed to be responsible for the dissimilar flow separation behaviour along physically similar peak- or trough-planes.
Abstract: This article presents a numerical analysis of a turbulent flow past DTMB 4119 marine propeller by the means of RANS approach; the propeller designed at David Taylor Model Basin in USA. The purpose of this study is to predict the hydrodynamic performance of the marine propeller, it aims also to compare the results obtained with the experiment carried out in open water tests; a periodical computational domain was created to reduce the unstructured mesh size generated. The standard kw turbulence model for the simulation is selected; the results were in a good agreement. Therefore, the errors were estimated respectively to 1.3% and 5.9% for KT and KQ.
Abstract: In this study, the three-dimensional cavitating
turbulent flow in a complete Francis turbine is simulated using
mixture model for cavity/liquid two-phase flows. Numerical analysis
is carried out using ANSYS CFX software release 12, and standard k-ε
turbulence model is adopted for this analysis. The computational
fluid domain consist of spiral casing, stay vanes, guide vanes, runner
and draft tube. The computational domain is discretized with a threedimensional
mesh system of unstructured tetrahedron mesh. The
finite volume method (FVM) is used to solve the governing equations
of the mixture model. Results of cavitation on the runner’s blades
under three different boundary conditions are presented and
discussed. From the numerical results it has been found that the
numerical method was successfully applied to simulate the cavitating
two-phase turbulent flow through a Francis turbine, and also
cavitation is clearly predicted in the form of water vapor formation
inside the turbine. By comparison the numerical prediction results
with a real runner; it’s shown that the region of higher volume
fraction obtained by simulation is consistent with the region of runner
cavitation damage.
Abstract: In this paper numerical studies have been carried out
to examine the pre-ignition flow features of high-performance solid
propellant rocket motors with two different port geometries but with
same propellant loading density. Numerical computations have been
carried out using a validated 3D, unsteady, 2nd-order implicit, SST k-
ω turbulence model. In the numerical study, a fully implicit finite
volume scheme of the compressible, Reynolds-Averaged, Navier-
Stokes equations is employed. We have observed from the numerical
results that in solid rocket motors with highly loaded propellants
having divergent port geometry the hot igniter gases can create preignition
pressure oscillations leading to thrust oscillations due to the
flow unsteadiness and recirculation. We have also observed that the
igniter temperature fluctuations are diminished rapidly thereby
reaching the steady state value faster in the case of solid propellant
rocket motors with convergent port than the divergent port
irrespective of the igniter total pressure. We have concluded that the
prudent selection of the port geometry, without altering the propellant
loading density, for damping the total temperature fluctuations within
the motor is a meaningful objective for the suppression and control of
instability and/or thrust oscillations often observed in solid propellant
rocket motors with non-uniform port geometry.
Abstract: The present study applies the inverse method and
three-dimensional CFD commercial software in conjunction with the
experimental temperature data to investigate the heat transfer and fluid
flow characteristics of the plate-fin heat sink in a rectangular closed
enclosure. The inverse method with the finite difference method and
the experimental temperature data is applied to determine the
approximate heat transfer coefficient. Later, based on the obtained
results, the zero-equation turbulence model is used to obtain the heat
transfer and fluid flow characteristics between two fins. T0 validate
the accuracy of the results obtained, the comparison of the heat transfer
coefficient is made. The obtained temperature at selected
measurement locations of the fin is also compared with experimental
data. The effect of the height of the rectangular enclosure on the
obtained results is discussed.
Abstract: In this paper comprehensive studies have been carried
out for the design optimization of a waste heat recovery system for
effectively utilizing the domestic air conditioner heat energy for
producing hot water. Numerical studies have been carried for the
geometry optimization of a waste heat recovery system for domestic
air conditioners. Numerical computations have been carried out using
a validated 2d pressure based, unsteady, 2nd-order implicit, SST k-ω
turbulence model. In the numerical study, a fully implicit finite
volume scheme of the compressible, Reynolds-Averaged, Navier-
Stokes equations is employed. At identical inflow and boundary
conditions various geometries were tried and effort has been taken for
proposing the best design criteria. Several combinations of pipe line
shapes viz., straight and spiral with different number of coils for the
radiator have been attempted and accordingly the design criteria has
been proposed for the waste heat recovery system design. We have
concluded that, within the given envelope, the geometry optimization
is a meaningful objective for getting better performance of waste heat
recovery system for air conditioners.
Abstract: Numerical studies have been carried out using a
validated two-dimensional standard k-omega turbulence model for
the design optimization of a thrust vector control system using shock
induced self-impinging supersonic secondary double jet. Parametric
analytical studies have been carried out at different secondary
injection locations to identifying the highest unsymmetrical
distribution of the main gas flow due to shock waves, which produces
a desirable side force more lucratively for vectoring. The results from
the parametric studies of the case on hand reveal that the shock
induced self-impinging supersonic secondary double jet is more
efficient in certain locations at the divergent region of a CD nozzle
than a case with supersonic single jet with same mass flow rate. We
observed that the best axial location of the self-impinging supersonic
secondary double jet nozzle with a given jet interaction angle, built-in
to a CD nozzle having area ratio 1.797, is 0.991 times the primary
nozzle throat diameter from the throat location. We also observed
that the flexible steering is possible after invoking ON/OFF facility to
the secondary nozzles for meeting the onboard mission requirements.
Through our case studies we concluded that the supersonic self-impinging
secondary double jet at predesigned jet interaction angle
and location can provide more flexible steering options facilitating
with 8.81% higher thrust vectoring efficiency than the conventional
supersonic single secondary jet without compromising the payload
capability of any supersonic aerospace vehicle.
Abstract: This paper reports the development and application of a 2D1 depth-averaged model. The main goal of this contribution is to apply the depth averaged equations to a wind park model in which the treatment of the geometry, introduced on the mathematical model by the mass and momentum source terms. The depth-averaged model will be used in future to find the optimal position of wind turbines in the wind park. κ − ε and 2D LES turbulence models were consider in this article. 2D CFD2 simulations for one hill was done to check the depth-averaged model in practise.
Abstract: Steady three-dimensional and two free surface waves
generated by moving bodies are presented, the flow problem to be
simulated is rich in complexity and poses many modeling challenges
because of the existence of breaking waves around the ship hull, and
because of the interaction of the two-phase flow with the turbulent
boundary layer. The results of several simulations are reported. The
first study was performed for NACA0012 of hydrofoil with different
meshes, this section is analyzed at h/c= 1, 0345 for 2D. In the second
simulation a mathematically defined Wigley hull form is used to
investigate the application of a commercial CFD code in prediction of
the total resistance and its components from tangential and normal
forces on the hull wetted surface. The computed resistance and wave
profiles are used to estimate the coefficient of the total resistance for
Wigley hull advancing in calm water under steady conditions. The
commercial CFD software FLUENT version 12 is used for the
computations in the present study. The calculated grid is established
using the code computer GAMBIT 2.3.26. The shear stress k-ωSST
model is used for turbulence modeling and the volume of fluid
technique is employed to simulate the free-surface motion. The
second order upwind scheme is used for discretizing the convection
terms in the momentum transport equations, the Modified HRIC
scheme for VOF discretization. The results obtained compare well
with the experimental data.