Search results

The purpose of this paper is to numerically investigate the three-dimensional (3D) reacting turbulent two-phase flow field of a scaled swirl-stabilized gas turbine combustor using the commercial computational fluid dynamic (CFD) software ANSYS FLUENT. The first scope of the study aims to explicitly compare the predictive capabilities of two turbulence models namely Unsteady Reynolds Averaged Navier-Stokes and Scale Adaptive Simulation for a reasonable trade-off between accuracy of results and global computational cost when applied to simulate swirl-stabilized spray combustion. The second scope of the study is to couple chemical reactions to the turbulent flow using a realistic chemistry model and also to model the local chemical non-equilibrium(NEQ) effects caused by turbulent strain such as flame stretching.

Standard Eulerian and Lagrangian formulations are used to describe both gaseous and liquid phases, respectively. The computing method includes a two-way coupling in which phase properties and spray source terms are interchanging between the two phases within each coupling time step. The fuel used is liquid jet-A1 which is injected in the form of a polydisperse spray and the droplet evaporation rate is calculated using the infinite conductivity model. One-component (n-decane) and two-component fuels (n-decane+toluene) are used as jet-A1 surrogates. The combustion model is based on the mean mixture fraction and its variance, and a presumed-probability density function is used to model turbulent-chemistry interactions. The instantaneous thermochemical state necessary for the chemistry tabulation is determined by using initially the equilibrium (EQ) assumption and thereafter, detailed NEQ calculations through the steady flamelets concept. The combustion chemistry of these surrogates is represented through a reduced chemical kinetic mechanism (CKM) comprising 1,045 reactions among 139 species, derived from the detailed jet-A1 surrogate model, JetSurf 2.0 using a sensitivity based method, Alternate Species Elimination.

Numerical results of the gas velocity, the gas temperature and the species molar fractions are compared with their experimental counterparts obtained from a steady state flame available in the literature. It is observed that, SAS coupled to the tabulated flamelet-based chemistry, predicts reasonably the main flame trends, while URANS even provided with the same combustion model and computing resources, leads to a poor prediction of the global flame trends, emphasizing the asset of a proper resolution when simulating spray flames.

The steady flamelet model even coupled with a robust turbulence model does not reproduce accurately the trend of species with slow oxidation kinetics such as CO and H2, because of the restrictiveness of the solutions space of flamelet equations and the assumption of unity Lewis for all species.

This work is adding a contribution for spray flame modeling and can be seen as an extension to the significant efforts for the modeling of gaseous flames using robust turbulence models coupled with the tabulated flamelet-based chemistry approach to considerably reduce computing cost. The exclusive use of a commercial CFD code widely used in the industry allows a direct application of this simulation approach to industrial configurations while keeping computing cost reasonable.

This study is useful to engineers interested in designing combustors of gas turbines and others combustion systems fed with liquid fuels.

This paper aims to investigate numerical solution of transonic flow around NACA0012 airfoil under sinusoidal pitch oscillation. Accordingly, effects of the amplitude and frequency of oscillations on aerodynamic coefficients are evaluated and the efficiency of the turbulent models, K-ω shear-stress transport (SST), scale adaptive simulation (SAS) and delayed detached eddy simulation (DDES), in simulation of the nonlinear phenomena – i.e. the interaction between shock and boundary layer and the shock oscillations – is studied.

K-ω SST, SAS and DDES models are used as turbulence approaches. The numerical results are compared with available experimental and numerical information.

According to the results inside the buffet boundaries, the DDES turbulent model expresses results that are more appropriate; however, SAS and SST models are not efficient enough in evaluating the characteristics of nonlinear flow.

In this research study, hybrid RANS-LES turbulence model is engaged to simulate transonic flow around pitching NACA0012 airfoil, and results are compared to the SAS and Reynolds Average Navier–Stocks simulations as well as available numerical and experimental data. In addition, effects of the amplitude and frequency of oscillations on aerodynamic coefficients are evaluated in buffet regions.

The purpose of this study is to propose a precise and standardized strategy for numerically simulating vehicle aerodynamics.

Error sources in computational fluid dynamics were analyzed. Additionally, controllable experiential and discretization errors, which significantly influence the calculated results, are expounded upon. Considering the airflow mechanism around a vehicle, the computational efficiency and accuracy of each solution strategy were compared and analyzed through numerous computational cases. Finally, the most suitable numerical strategy, including the turbulence model, simplified vehicle model, calculation domain, boundary conditions, grids and discretization scheme, was identified. Two simplified vehicle models were introduced, and relevant wind tunnel tests were performed to validate the selected strategy.

Errors in vehicle computational aerodynamics mainly stem from the unreasonable simplification of the vehicle model, calculation domain, definite solution conditions, grid strategy and discretization schemes. Using the proposed standardized numerical strategy, the simulated steady and transient aerodynamic characteristics agreed well with the experimental results.

Building upon the modified Low-Reynolds Number k-e model and Scale Adaptive Simulation model, to the best of the authors’ knowledge, a precise and standardized numerical simulation strategy for vehicle aerodynamics is proposed for the first time, which can be integrated into vehicle research and design.

Convective heat transfer features of a turbulent slot jet impingement are comprehensively studied using two different computational approaches, namely, URANS (unsteady Reynolds-averaged Navier–Stokes equations) and SAS (scale-adaptive simulation). Turbulent slot jet impingement heat transfer is used where a considerable heat transfer enhancement is required, and computationally, it is a quite challenging flow configuration.

Customized OpenFOAM 4.1, an open-access computational fluid dynamics (CFD) code, is used for SAS (SST-SAS k-ω) and URANS (standard k-ε and SST k-ω) computations. A low-Re version of the standard k-ε model is used, and other models are formulated for good wall-refined calculations. Three turbulence models are formulated in OpenFOAM 4.1 with second-order accurate discretization schemes.

It is observed that the profiles of the streamwise turbulence are under-predicted at all the streamwise locations by SST k-ω and SST SAS k-ω models, but follow similar trends as in the reported results. The standard k-ε model shows improvements in the predictions of the streamwise turbulence and mean streamwise velocity profiles in the zone of outer wall jet. Computed profiles of Nusselt number by SST k-ω and SST-SAS k-ω models are nearly identical and match well with the reported experimental results. However, the standard k-ε model does not provide a reasonable profile or quantification of the local Nusselt number.

Hybrid turbulence model is suitable for efficient CFD computations for the complex flow problems. This paper deals with a detailed comparison of the SAS model with URANS and LES for the first time in the literature. A thorough assessment of the computations is performed against the results reported using experimental and large eddy simulations techniques followed by a detailed discussion on flow physics. The present results are beneficial for scientists working with hybrid turbulence models and in industries working with high-efficiency cooling/heating system computations.

The purpose of this paper is to outline the extensive multi-scale and multi-physics challenges when simulating future aircraft and offer strategies to help deal with some of these challenges.

To help with the multi-scale challenges, in a hierarchical, zonal fashion both the handling of turbulence and geometry is considered.

Such modelling of geometry is necessary to help deal with the increasingly coupled nature of many aerodynamic problems more economically and the drive towards considering ever increasing levels of geometrical complexity/scale.

The proposed unified framework could be exploited all the way, through initial fast preliminary design to final numerical test involving various bespoke combinations of hierarchical components.

The purpose of this paper is to study the pressure fluctuation characteristics in the sidewall gaps of a centrifugal dredging pump in detail and discover the excitation sources.

An unsteady numerical simulation with shear–stress transport–scale-adaptive simulation (SAS-SST) model was conducted for a centrifugal pump considering the sidewall gaps. The numerical codes were validated by a model test carried out in China Water Resources Beifang Investigation, Design and Research Co., Ltd. Fast Fourier transform was used to obtain the frequency components of the pressure fluctuation.

Pressure fluctuation characteristics inside the pump were analyzed for a condition near the design point. In the sidewall gaps, the circumferential, radial and axial distribution of the pressure fluctuation amplitude follow different laws. The non-axisymmetrical distribution of pressure fluctuation in the sidewall gaps shows that the unsteady flow in the volute casing which has a non-axisymmetrical geometry imposes an evident effect on the flow field in the sidewall gaps and the interaction between the main flow and the clearance flow cannot be neglected. There are several frequency components appearing as the dominant frequencies at different locations in the sidewall gaps, but the relatively stronger pressure fluctuations are all dominated by the rotating frequency. It indicates that the rotating impeller, which originally makes the shrouds rotate, is the primarily excitation source of the pressure fluctuations in the sidewall gaps.

The pressure fluctuation characteristics in the sidewall gaps of centrifugal pumps were first comprehensively analyzed. Unsteady flows in the sidewall gaps should be considered during the design and operation of centrifugal pumps.

The purpose of this study is to present and demonstrate a numerical method for solving chemically reacting flows. These are important for energy conversion devices, which rely on chemical reactions as their operational mechanism, with heat generated from the combustion of the fuel, often gases, being converted to work.

The numerical study of such flows requires the set of Navier-Stokes equations to be extended to include multiple species and the chemical reactions between them. The numerical method implemented in this study also accounts for changes in the material properties because of temperature variations and the process to handle steep spatial fronts and stiff source terms without incurring any numerical instabilities. An all-speed numerical framework is used through simple low-dissipation advection upwind splitting (SLAU) convective scheme, and it has been extended in a multi-component species framework on the in-house density-based flow solver. The capability of solving turbulent combustion is also implemented using the Eddy Dissipation Concept (EDC) framework and the recent k-kl turbulence model.

The numerical implementation has been demonstrated for several stiff problems in laminar and turbulent combustion. The laminar combustion results are compared from the corresponding results from the Cantera library, and the turbulent combustion computations are found to be consistent with the experimental results.

This paper has extended the single gas density-based framework to handle multi-component gaseous mixtures. This paper has demonstrated the capability of the numerical framework for solving non-reacting/reacting laminar and turbulent flow problems. The all-speed SLAU convective scheme has been extended in the multi-component species framework, and the turbulent model k-kl is used for turbulent combustion, which has not been done previously. While the former method provides the capability of solving for low-speed flows using the density-based method, the later is a length-scale-based method that includes scale-adaptive simulation characteristics in the turbulence modeling. The SLAU scheme has proven to work well for unsteady flows while the k-kL model works well in non-stationary turbulent flows. As both these flow features are commonly found in industrially important reacting flows, the convection scheme and the turbulence model together will enhance the numerical predictions of such flows.

This study aims to deal with the large eddy simulation (LES) of an ignition sequence and the resulting steady combustion in a swirl-stabilized liquid-fueled combustor. Particular attention is paid to the ease of handling the numerical tool, the accuracy of the results and the reasonable computational cost involved. The primary aim of the study is to appraise the ability of the newly developed computational fluid dynamics (CFD) methodology to retrieve the spark-based flame kernel initiation, its propagation until the full ignition of the combustion chamber, the flame stabilization and the combustion processes governing the steady combustion regime.

The CFD model consists of an LES-based spray module coupled to a subgrid-scale ignition model to capture the flame kernel initiation and the early stage of the flame kernel growth, and a combustion model based on the mixture fraction-progress variable formulation in the line of the flamelet generated manifold (FGM) method to retrieve the subsequent flame propagation and combustion properties. The LES-spray module is based on an Eulerian-Lagrangian approach and includes a fully two-way coupling at each time step to account for the interactions between the liquid and the gaseous phases. The Wall-Adapting Local Eddy-viscosity (WALE) model is used for the flow field while the eddy diffusivity model is used for the scalar fluxes. The fuel is liquid kerosene, injected in the form of a polydisperse spray of droplets. The spray dynamics are tracked using the Lagrangian procedure, and the phase transition of droplets is calculated using a non-equilibrium evaporation model. The oxidation mechanism of the Jet A-1 surrogate is described through a reduced reaction mechanism derived from a detailed mechanism using a species sensitivity method.

By comparing the numerical results with a set of published data for a swirl-stabilized spray flame, the proposed CFD methodology is found capable of capturing the whole spark-based ignition sequence in a liquid-fueled combustion chamber and the main flame characteristics in the steady combustion regime with reasonable computing costs.

The proposed CFD methodology simulates the whole ignition sequence, namely, the flame kernel initiation, its propagation to fully ignite the combustion chamber, and the global flame stabilization. Due to the lack of experimental ignition data on this liquid-fueled configuration, the ability of the proposed CFD methodology to accurately predict ignition timing was not quantitatively assessed. It would, therefore, be interesting to apply this CFD methodology to other configurations that have experimental ignition data, to quantitatively assess its ability to predict the ignition timing and the flame characteristics during the ignition sequence. Such further investigations will not only provide further validation of the proposed methodology but also will potentially identify its shortfalls for better improvement.

This CFD methodology is developed by customizing a commercial CFD code widely used in the industry. It is, therefore, directly applicable to practical configurations, and provides not only a relatively straightforward approach to predict an ignition sequence in liquid-fueled combustion chambers but also a robust way to predict the flame characteristics in the steady combustion regime as significant improvements are noticed on the prediction of slow species.

The incorporation of the subgrid ignition model paired with a combustion model based on tabulated chemistry allows reducing computational costs involved in the simulation of the ignition phase. The incorporation of the FGM-based tabulated chemistry provides a drastic reduction of computing resources with reasonable accuracy. The CFD methodology is developed using the platform of a commercial CFD code widely used in the industry for relatively straightforward applicability.

This paper aims to improve performance prediction and to acquire more detailed flow structures so as to analyze the turbulence in complex rotor-stator flow.

Hydraulic retarder as typical fluid machinery was numerically investigated by using hybrid Reynolds-averaged Navier–Stokes (RANS)/large eddy simulation (LES) models CIDDES Algebraic Wall-Modeled Large Eddy Simulation (LES) (WMLES) S-Ω and dynamic hybrid RANS/LES (DHRL). The prediction results were compared and analyzed with a RANS model shear stress transport (SST) k-omega which was a recommended choice in engineering.

The numerical results were verified by experiment and indicated that the predicted values for three hybrid turbulence models were more accurate. Then, the transient flow field was further analyzed visually in terms of turbulence statistics, Reynolds number, pressure-streamline, vortex structure and eddy viscosity ratio. The results indicated that HRL approaches could capture unsteady flow phenomena.

This study achieves both in performance prediction improvement and better flow mechanism understanding. The computational fluid dynamics (CFD) could be used instead of flow visualization to a certain extent. The improved CFD method, the fine computational grid and the reasonable simulation settings jointly enhance the application of CFD in the rotor-stator flow.

The improvement was quite encouraging compared with the reported literatures, contributing to the CFD playing a more important role in the flow machinery. DHRL provided the detailed explanation of flow transport between rotor and stator, which was not reported before. Through it, the flow mechanism can be better understood.

This paper aims to conduct, a detailed investigation of various Reynolds averaged Navier–Stokes (RANS) models to study their performance in attached and separated flows. The turbulent flow over two airfoils, namely, National Advisory Committee for Aeronautics (NACA)-0012 and National Aeronautics and Space Administration (NASA) MS(1)-0317 with a static stall setup at a Reynolds number of 6 million, is chosen to investigate these models. The pre-stall and post-stall regions, which are in the range of angles of attack 0°–20°, are simulated.

RANS turbulence models with the Boussinesq approximation are the most commonly used cost-effective models for engineering flows. Four RANS models are considered to predict the static stall of two airfoils: Spalart–Allmaras (SA), Menter’s k – ω shear stress transport (SST), k – kL and SA-Bas Cakmakcioglu modified (BCM) transition model. All the simulations are performed on an in-house unstructured-grid compressible flow solver.

All the turbulence models considered predicted the lift and drag coefficients in good agreement with experimental data for both airfoils in the attached pre-stall region. For the NACA-0012 airfoil, all models except the SA-BCM over-predicted the stall angle by 2°, whereas SA-BCM failed to predict stall. For the NASA MS(1)-0317 airfoil, all models predicted the lift and drag coefficients accurately for attached flow. But the first three models showed even further delayed stall, whereas SA-BCM again did not predict stall.

The numerical results at high Re obtained from this work, especially that of the NASA MS(1)-0317, are new to the literature in the knowledge of the authors. This paper highlights the inability of RANS models to predict the stall phenomenon and suggests a need for improvement in modeling flow physics in near- and post-stall flows.