Search results

A mathematical model has been developed to study incompressible, isothermal, turbulent, confined, swirling flows. The model solves the conservation equations of mass, momentum, and two additional equations for the turbulent kinetic energy and the rate of dissipation of turbulent kinetic energy. The numerical predictions show a recirculation zone in the form of a one‐celled toroidal vortex at the combustor centreline. High levels of turbulence characterize the recirculation zone. The length, diameter and maximum velocity of the recirculation zone first decrease and then increase as the magnitude of the outer swirl number is first decreased from counter‐swirl to zero and then increased to co‐swirl flow conditions. Counter‐swirl produces steeper velocity gradients at the inter‐jet shear layer and promotes faster mixing than co‐swirl. The numerical results also indicate that the mass of the recirculation zone first decreases and then increases as the outer swirl number is first decreased from counter‐swirl to zero and then increased to co‐swirl conditions. The diameter, maximum velocity and mass of the recirculation zone are monotonically increasing functions of the inner jet swirl number. The recirculation zone length, diameter and mass are almost independent of the Reynolds number and outer‐to‐inner jet axial velocity ratio.

An orthogonal array technique is used in the present work to investigate, numerically, the effects of the swirler and the primary jets on the characteristics of the recirculation zone of a can‐type gas turbine combustor. The computer code used for this purpose is first validated with the available experimental data. The effects of change in the percentage flow rate through the swirler, the swirl number, the hub diameter of the swirler and the diameter of the primary injection holes (which influences the velocity of the jets) are estimated first. It is found that the flow rate through the swirler and the size of the primary injection hole have much more influence on the characteristics of the recirculation zone than the swirl number and the hub diameter of the swirler. But the earlier studies show that for a given flow rate through the swirler, the swirl number and swirler geometry have considerable influence on the characteristics of the recirculation zone in the absence of primary jets. Therefore it is inferred that there may be a critical point, based on the ratio of flow rate through the swirler to that of primary holes, beyond which the effects of swirl number and the swirler geometry dominate the effect of primary jets in determining the characteristics of the recirculation zone. This critical point is determined by gradually reducing the flow through the primary holes. It is found that, initially, the recirculation ratio (ratio of the mass of fluid recirculated to that sum of the mass flow rate through the swirler and through that of primary hole) reduces because of weakening of the primary jets but after the critical point it increases because of the swirler effect taking over the role of providing the recirculation. It is also observerd that the length of the recirculation zone increases as the strength of the primary jets reduces.

Over the past few decades, the flow around circular cylinders has been one of the highly researched topics in the field of offshore engineering and fluid-structure interaction (FSI). In the current study, numerical simulations for flow around a fixed circular cylinder are performed at Reynolds number (Re) = 3900 with the LES method using the ICEM-CFD and ANSYS Fluent tool for meshing and analysis, respectively. Previously, similar studies have been conducted at the same Reynolds number, but there have been discrepancies in the results, particularly in calculating the recirculation length and angle of separation. In addition, the purpose of this study is to address the impact of time interval averaging to obtain the fully converged solution.

This study presents the LES method, using the ICEM-CFD and ANSYS fluent tool for meshing and analysis.

In the current study, turbulence statistics are sampled for 25, 50, 75 and 100 vortex-shedding cycles with the CFL value O (1). The recirculation length, angle of separation, hydrodynamic coefficients and the wake behind the cylinder are investigated up to ten diameters. The drag coefficient and Strouhal number are observed to be less sensitive, whereas the recirculation length appeared to be highly dependent on the average time statistics and the non-dimensional time step. Similarly, the mean streamwise and cross-flow velocity are observed to be sensitive to the average time statistics and non-dimensional time step in the wake region near the cylinder.

In the current investigation, turbulence statistics are sampled for 25, 50, 75 and 100 vortex-shedding cycles with the CFL value O (1), using large eddy simulation method at Re = 3900 around a circular cylinder. The impact of time interval averaging to obtain the fully converged mean flow field is addressed. No such consideration is yet published in the literature.

Fluid flows in pipes whose cross-sectional area are increasing in the stream-wise direction are prone to separation of the recirculation region. This paper aims to investigate such fluid flow in expansion pipe systems using direct numerical simulations. The flow in circular diverging pipes with different diverging half angles, namely, 45, 26, 14, 7.2 and 4.7 degrees, are considered. The flow is fed by a fully developed laminar parabolic velocity profile at its inlet and is connected to a long straight circular pipe at its downstream to characterise recirculation zone and skin friction coefficient in the laminar regime. The flow is considered linearly stable for Reynolds numbers sufficiently below natural transition. A perturbation is added to the inlet fully developed laminar velocity profile to test the flow response to finite amplitude disturbances and to characterise sub-critical transition.

Direct numerical simulations of the Navier–Stokes equations have been solved using a spectral element method.

It is found that the onset of disordered motion and the dynamics of the localised turbulence patch are controlled by the Reynolds number, the perturbation amplitude and the half angle of the pipe.

The authors clarify different stages of flow behaviour under the finite amplitude perturbations and shed more light to flow physics such as existence of Kelvin–Helmholtz instabilities as well as mechanism of turbulent puff shedding in diverging pipe flows.

For higher swirling flows (swirl > 0.5), flow confinement significantly impacts fluid flow, flame stability, flame length and heat transfer, especially when the confinement ratio is less than 9. Past numerical studies on helical axial swirler type systems are limited to non-reacting or reacting flows type Reynolds averaged Navier Stokes closure models, mostly are non-parametric studies. Effects of parametric studies like swirl angle and confinement on the unsteady flow field, either numerical or experimental, are very minimal. The purpose of this paper is to document modeling practices for a large eddy simulation (LES) type grid, predict the confinement effects of a single swirler lean direct injection (LDI) system and validate with literature data.

The first part of the paper discusses the approach followed for numerical modelling of LES with the minimum number of cells required across critical sections to capture the spectrum of turbulent energy with good accuracy. The numerical model includes all flow developing sections of the LDI swirler, right from the axial setting chamber to the exit of the flame tube, and its length is effectively modelled to match the experimental data. The computational model predicts unsteady features like vortex breakdown bubble, represented by a strong recirculation zone anchored downstream of the fuel nozzle. It is interesting to note that the LES is effective in predicting the secondary recirculation zones in the divergent section as well as at the corners of the tube wall.

The predictions of a single helical axial swirler with a vane tip angle of 60°, with a duct size of 2 × 2 square inches, are compared with the experimental data at several axial locations as well as with centerline data. Both mean and unsteady turbulent quantities obtained through the numerical simulations are validated with the experimental data (Cai et al., 2005). The methodology is extended to the confinements effect on mean flow characteristics. The time scale and length scale are useful parameters to get the desired results. The results show that with an increase in the confinement ratio, the recirculation length increases proportionally. A sample of three cases has been documented in this paper.

The novelty of the paper is the modelling practices (grid/unsteady models) for a parametric study of LDI are established, and the mean confinement effects are validated with experimental data. The spectrum of turbulent energies is well captured by LES, and trends are aligned with experimental data. The methodology can be extended to reacting flows also to study the effect of swirl angle, fuel injection on aerodynamics, droplet characteristics and emissions.

This paper aims to perform numerical simulations through different shaped double stenoses in a vascular tube for a better understanding of arterial blood flow patterns, and their possible role during the progression of atherosclerosis. The dynamics of flow features have been studied by wall pressure, streamline contour and wall shear stress distributions for all models.

A finite volume method has been employed to solve the governing equations for the two‐dimensional, steady, laminar flow of an incompressible and Newtonian fluid.

The paper finds that impact of pressure drop, reattachment length and peak wall shear stress at each restriction primarily depends upon percentage of restriction, if restriction spacing is sufficient. The quantum of impact of pressure drop, reattachment length and peak wall shear stress is much effected for smaller restriction spacing. If recirculating bubble of first restriction merges with the recirculating bubble formed behind the second restriction in this smaller restriction spacing. The similar effect of smaller restriction spacing is observed, if Reynolds number increases also.

The effect of different shaped stenoses, restriction spacing and Reynolds number on the flow characteristics has been investigated and the role of all the flow characteristics on the progression of the disease, atherosclerosis, is discussed.

When a reflected shock interacts with the boundary layer in a shock tube, the shock bifurcation occurs near the walls. Although the study of the shock bifurcation has been carried out by many researchers for several decades, little attention has been devoted to investigate the instability pattern of the bifurcation. This research work aims to successfully capture the asymmetry of the whole flow field, and attempt to achieve the instability mechanism of the shock bifurcation by a direct numerical simulation of the reflected shock wave/boundary layer interaction at Ma = 1.9. In addition, the reason for the formation of the bifurcated structure is also explored.

The spatial and temporal evolution of the shock bifurcation is obtained by solving the two-dimensional compressible Navier–Stokes equations using a seventh-order accurate weighted essentially non-oscillatory (WENO) scheme and a three-step Runge–Kutta time advancing approach.

The results show that the formation of shock bifurcation is mainly because of the shock/gradient field interaction, and the height of the bifurcated foot increases with the growth of the shock intensity and the gradient field. The unsteady asymmetry of the upper and bottom shock bifurcated structures is because of the vortex shedding with high frequency in the rear recirculation zone, which leads to the fluctuation of the recirculation area. The vortex shedding process behind the bifurcated structure closely resembles the Karman vortex street formed by the flow around the cylinder. The dimensionless vortex shedding frequency varies between 0.01 and 0.02. In comparison to the scenario at Ma = 1.9, the occurring time of instability is delayed and the upper and bottom bifurcated feet intersect in a relatively short time at Ma = 3.5. The region behind the bifurcated shock is a transitional flow field containing obvious cell structures and “isolated islands.”

This paper discovers an unsteady flow pattern of the shock bifurcation, and the mechanism of this instability in the reflected shock/boundary layer interaction is revealed in detail.

One major challenge in turbulent flow applications is to control the recirculation zone behind the backward‐facing step (BFS). One simple idea to do so is to modify the original BFS geometry, of course, without causing adverse or undesirable impacts on the original characteristics of the primary stream. The main objective of this work is to examine the solidity of the recirculation zone behind several different geometries which are slightly to moderately different from the original BFS geometry.

The implemented modifications cause complicated irregularities at the boundaries of the domain. The experience shows that the mesh distribution around these irregularities plays a critical role in the accuracy of the numerical solutions. To achieve the most accurate solutions with the least computational efforts, we use a robust hybrid strategy to distribute the computational grids in the domain. Additionally, a suitable numerical algorithm capable of handling hybrid grid topologies is properly extended to analyze the flow field. The current fully implicit method utilizes a physical pressure‐based upwinding scheme capable of working on hybrid mesh.

The extended algorithm is very robust and obtains very accurate solutions for the complex flow fields despite utilizing very coarse grid resolutions. Additionally, different proposed geometries revealed very similar separated regions behind the step and performed minor differences in the location of the reattachment points.

The current study is fulfilled two‐dimensionally. However, the measurements in testing regular BFS problems have shown that the separated shear layer behind the step is not affected by 3D influences provided that the width of channel is sufficiently wide. A similar conclusion is anticipated here.

The problem occurs in the pipe and channel expansions, combustion chambers, flow over flying objects with abrupt contraction on their external surfaces, etc.

A novel pressure‐based upwinding strategy is properly employed to solve flow on multiblocked hybrid grid topologies. This strategy takes into account the physics associated with all the transports in the flow field. To study the impact of shape improvement, several modified BFS configurations were suggested and examined. These configurations need only little additional manufacturing cost to be fabricated.

The purpose of this paper is to carry out a numerical study on the dynamic and thermal behavior of a fluid with a constant property and flowing turbulently through a two-dimensional horizontal rectangular channel. The upper surface was put in a constant temperature condition, while the lower one was thermally insulated. Two transverse, solid-type obstacles, having different shapes, i.e. flat rectangular and V-shaped, were inserted into the channel and fixed to the top and bottom walls of the channel, in a periodically staggered manner to force vortices to improve the mixing, and consequently the heat transfer. The flat rectangular obstacle was put in the first position and was placed on the hot top wall of the channel. However, the second V-shaped obstacle was placed on the insulated bottom wall, at an attack angle of 45°; its position was varied to find the optimum configuration for optimal heat transfer.

The fluid is considered Newtonian, incompressible with constant properties. The Reynolds averaged Navier–Stokes equations, along with the standard k-epsilon turbulence model and the energy equation, are used to control the channel flow model. The finite volume method is used to integrate all the equations in two-dimensions; the commercial CFD software FLUENT along with the SIMPLE-algorithm is used for pressure-velocity coupling. Various values of the Reynolds number and obstacle spacing were selected to perform the numerical runs, using air as the working medium.

The channel containing the flat fin and the 45° V-shaped baffle with a large Reynolds number gave higher heat transfer and friction loss than the one with a smaller Reynolds number. Also, short separation distances between obstacles provided higher values of the ratios Nu/Nu₀ and f/f₀ and a larger thermal enhancement factor (TEF) than do larger distances.

This is an original work, as it uses a novel method for the improvement of heat transfer in completely new flow geometry.

Stent malapposition is one of the most significant precursors of stent thrombosis and restenosis. Adverse haemodynamics may play a key role in establishing these diseases, although numerical studies have used idealised drug transport models to show that drug transport from malapposed drug-eluting stent struts can be significant. This paper aims to study whether drug transport from malapposed struts is truly significant. Another aim is to see whether a streamlined strut profile geometry – with a 61% smaller coating but a 32% greater coating-tissue contact area – can mitigate the adverse haemodynamics associated with stent malapposition while enhancing drug uptake.

Two- and three-dimensional computational fluid dynamics simulations were used in this study. Unlike past simulations of malapposed drug-eluting stent struts, a qualitatively validated drug-transport model which simulates the non-uniform depletion of drug within the drug coating was implemented.

It was shown that even a 10-µm gap between the strut and tissue dramatically reduces drug uptake after 24 h of simulated drug transport. Furthermore, the streamlined strut profile was shown to minimise the adverse haemodynamics of malapposed and well-apposed stent struts alike and enhance drug uptake.

Unlike prior numerical studies of malapposed stent struts, which did not model the depletion of drug in the drug coating, it was found that stent malapposition yields negligible drug uptake. The proposed semicircular-profiled strut was also shown to be advantageous from a haemodynamic and drug transport perspective.