Home  About us  Contact
 

Turbulent Flow (turbulent + flow)

Distribution by Scientific Domains
Engineering55%
Chemistry28%
Medical Sciences7%
Earth and Environmental Science2%
4 Other Domains8%

Kinds of Turbulent Flow

  • incompressible turbulent flow
    		•

    Terms modified by Turbulent Flow

  • turbulent flow condition
    		•

    Selected Abstracts

    Turbulent Flow Through Porous Media

    GROUND WATER, Issue 5 2001
    Douglas W. Barr
    The pressure driving flow through porous media must be equal to the viscous resistance plus the inertial resistance. Formulas are developed for both the viscous resistance and the inertial resistance. The expression for the coefficient of permeability consists of parameters which describe the characteristics of the porous medium and the permeating fluid and which, for unconsolidated isotropic granular media, are all measurable. A procedure is proposed for testing for the occurrence of turbulence and calculating the effective permeability when it occurs. The formulas are applied to a set of data from 588 permeameter runs ranging from laminar to highly turbulent. The equations fit the data from the permeameter closely through the laminar flow conditions and quite closely through the turbulent conditions. In the turbulent range, the plotting of the data separates into three distinct lines for each of the three shapes of particles used in the tests. For the porous medium and fluid of these tests, turbulence begins at a head gradient of about 0.1. [source]

    Combined PIV/PTV-Measurements for the Analysis of Bubble Interactions and Coalescence in a Turbulent Flow

    THE CANADIAN JOURNAL OF CHEMICAL ENGINEERING, Issue 3-4 2003
    Dirk Bröder
    Abstract In order to allow more reliable modeling of coalescence processes in turbulent bubbly flows, detailed experiments in a double loop reactor were performed. Narrow and essentially monomodal bubble size distributions in the range of 2 to 4 mm were created. For simultaneous measurements of bubble size, bubble velocity and liquid velocity a combined system of PIV and PTV was developed and applied. It was possible to determine bubble size distributions and mean, as well as fluctuating velocities for both phases. The spatial changes of the bubble size distribution, due to the influence of bubble coalescence, was analyzed and coalescence rates were calculated. Afin d'améliorer la fiabilité de la modélisation des procédés de coalescence dans les écoulements bouillonnants turbulents, des expériences poussées ont été menées dans un réacteur à double boucle. Des distributions de taille des bulles étroites, essentiellement monodisperses, dans la gamme de 2 à 4 mm, ont été créées. Pour les mesures simultanées de la taille des bulles, de la vitesse des bulles et de la vitesse du liquide, on a mis au point et appliqué un système combiné de PIV et PTV. Il a été possible de déterminer les distributions de taille des bulles et la moyenne mais également les vitesses fluctuantes pour les deux phases. Les changements spatiaux de la distribution de la taille des bulles dus à l'influence de la coalescence des bulles ont été analysés et les vitesses de coalescence ont été calculées. [source]

    Numerical Investigation of Turbulent Flow around a Rotating Stepped Cylinder for Corrosion Study

    THE CANADIAN JOURNAL OF CHEMICAL ENGINEERING, Issue 1 2003
    Kyung-Soo Yang
    Abstract Direct numerical simulation has been carried out for turbulent flow set up by a rotating cylinder with two backward-facing steps axisymmetrically mounted in the circumferential direction. This flow geometry creates a qualitatively similar flow pattern as observed near a sudden pipe expansion or a plane backward-facing step, characterized by flow separation and reattachment. A region of intense turbulence intensity and high wall-shear-stress fluctuations is formed in the recirculating region downstream of the step, where high mass-transfer capacity was also experimentally observed. Since corrosion is frequently mass-transfer controlled, our findings put forward this apparatus as a useful tool for future corrosion research. On a effectué une simulation numérique directe de l'écoulement turbulent créé par un cylindre rotatif ayant deux contractions axisymétriques dans la direction circonférentielle. Cette géométrie crée un profil d'écoulement qualitativement similaire à celui qu'on observe près d'une expansion de conduite soudaine ou d'une contraction planaire, caractérisés par la séparation et le ré-attachement de l'écoulement. Une région d'intense turbulence et de fortes fluctuations de contraintes de cisaillement pariétal se forment dans la région en recirculation en aval de la contraction, où une grande capacité de transfert de matière a également été observée expérimentalement. Étant donné que la corrosion dépend souvent du transfert de matière, nos résultats font la promotion de cet appareillage en tant qu'outil utile pour la recherche future sur la corrosion. [source]

    Simulation of Turbulent Flow in a Packed Bed

    CHEMICAL ENGINEERING & TECHNOLOGY (CET), Issue 5 2006
    B. Guo
    Abstract Numerous models for simulating the flow and transport in packed beds have been proposed in the literature with few reported applications. In this paper, several turbulence models for porous media are applied to the gas flow through a randomly packed bed and are examined by means of a parametric study against some published experimental data. These models predict widely different turbulent eddy viscosity. The analysis also indicates that deficiencies exist in the formulation of some model equations and selection of a suitable turbulence model is important. With this realization, residence time distribution and velocity distribution are then simulated by considering a radial profile of porosity and turbulence induced dispersion, and the results are in good agreement with the available experimental data. [source]

    An Improved LES on Dense Particle-Liquid Turbulent Flows Using Integrated Boltzmann Equations

    THE CANADIAN JOURNAL OF CHEMICAL ENGINEERING, Issue 2 2007
    Xuelin Tang
    Abstract An improved large eddy simulation (LES) using a dynamic second-order subgrid stress (SGS) model has been developed for simulating dense particle-liquid two-phase turbulent flows. The governing equations of each phase are obtained from a microscopic point of view, using the kinetic theory of molecular gas. They are derived by multiplying the Boltzmann equation of each phase by property parameters and integrating over the velocity space. An inter-particle collision term is included in the governing equation of the particle phase. Assuming a Maxwellian distribution of the velocity for particle-phase, an inter-particle collision term is derived. On a amélioré la simulation des grands tourbillons (LES) à l'aide d'un modèle de contraintes dans les mailles inférieures (SGS) de second ordre dynamique afin de simuler des écoulements turbulents diphasiques particules-liquide. Les équations gouvernantes de chaque phase sont obtenues d'un point de vue microscopique, au moyen de la théorie cinétique du gaz moléculaire. Elles sont calculées en multipliant l'équation de Boltzmann de chaque phase par des paramètres de propriétés et en intégrant sur l'intervalle des vitesses. Un terme de collision inter-particulaire est inclus dans l'équation gouvernante de la phase des particules. [source]

    Numerical analysis of turbulent flow separation in a rectangular duct with a sharp 180-degree turn by algebraic Reynolds stress model

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2008
    Hitoshi Sugiyama
    Abstract Turbulent flow in a rectangular duct with a sharp 180-degree turn is difficult to predict numerically because the flow behavior is influenced by several types of forces, including centrifugal force, pressure-driven force, and shear stress generated by anisotropic turbulence. In particular, this type of flow is characterized by a large-scale separated flow, and it is difficult to predict the reattachment point of a separated flow. Numerical analysis has been performed for a turbulent flow in a rectangular duct with a sharp 180-degree turn using the algebraic Reynolds stress model. A boundary-fitted coordinate system is introduced as a method for coordinate transformation to set the boundary conditions next to complicated shapes. The calculated results are compared with the experimental data, as measured by a laser-Doppler anemometer, in order to examine the validity of the proposed numerical method and turbulent model. In addition, the possibility of improving the wall function method in the separated flow region is examined by replacing the log-law velocity profile for a smooth wall with that for a rough wall. The analysis results indicated that the proposed algebraic Reynolds stress model can be used to reasonably predict the turbulent flow in a rectangular duct with a sharp 180-degree turn. In particular, the calculated reattachment point of a separated flow, which is difficult to predict in a turbulent flow, agrees well with the experimental results. In addition, the calculation results suggest that the wall function method using the log-law velocity profile for a rough wall over a separated flow region has some potential for improving the prediction accuracy. Copyright © 2007 John Wiley & Sons, Ltd. [source]

    Experimental study of a turbulent cross-flow near a two-dimensional rough wall with narrow apertures

    AICHE JOURNAL, Issue 10 2008
    Satya Mokamati
    Abstract Turbulent flow over a rough wall with suction or blowing is a common fluid mechanics problem that has many practical applications including pulp screening. To better understand, the complex hydrodynamics at the critical region near the surface of the wall, the streamwise mean and velocity fluctuations were determined experimentally using laser Doppler velocietry. The near-wall streamwise velocity fluctuations and local mean streamwise velocity were shown to be a strong function of the surface roughness, and the aperture and cross-flow velocities. A correlation for the mean velocity and the wall shear stress acting near the wall was determined. © 2008 American Institute of Chemical Engineers AIChE J, 2008 [source]

    Turbulent flows on forested hilly terrain: the recirculation region

    THE QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, Issue 625 2007
    D. Poggi
    Abstract A number of analytical and numerical studies employing first-order closure principles have suggested that canopy flows on gentle sinusoidal hills feature a recirculation region, situated on the lee side, that can dramatically affect scalar transfer between the biosphere and the atmosphere. To date, the onset of this region, and its effects on bulk flow properties, have not been experimentally investigated. We study the applicability of first-order closure schemes jointly with the properties of this recirculation region, using detailed laser Doppler anemometry (LDA) measurements. These experiments are conducted in a neutrally stratified boundary-layer flow within a large flume over a train of gentle and narrow hills. The canopy is composed of an array of vertical cylinders with a frontal-area index concentrated in the upper third, to resemble a tall hardwood forest at maximum leaf area. The LDA measurements are recorded for both sparse and dense canopies. We find that, while the onset of a recirculation region is ambiguous in the sparse-canopy case, it is well delineated in the dense-canopy case. This finding constitutes the first experimental evidence confirming the analytical and numerical model predictions concerning this region in dense canopies on gentle hills. Moreover, we show that the presence of the recirculation region can explain the anomalous pressure variation across the hill (first reported in numerical simulations) using an ,effective hill shape' function. Detailed momentum-flux measurements show, surprisingly, that the effective mixing length leff within the canopy and in the inner layer is not significantly affected by the recirculation region. We expected leff to be comparable to the size of the vortex responsible for the recirculation zone, but the measurements show that leff maintains its canonical canopy turbulence shape. Using laser-induced fluorescence (LIF) measurements, we find that the recirculation region is not characterized by a classical ,rotor', but by a highly intermittent zone with alternating positive and negative velocity values in the lower layers of the canopy. These LIF measurements may explain why leff maintains its canonical canopy turbulence shape in the recirculation region. The LIF measurements also show that the main mechanism for scalar transfer within the recirculation region is a sequence of accumulation,ejection episodes that are quasi-periodic in nature. Copyright © 2007 Royal Meteorological Society [source]

    Turbulent flow over a dune: Green River, Colorado

    EARTH SURFACE PROCESSES AND LANDFORMS, Issue 3 2005
    Jeremy G. Venditti
    Abstract Detailed echo-sounder and acoustic Doppler velocimeter measurements are used to assess the temporal and spatial structure of turbulent flow over a mobile dune in a wide, low-gradient, alluvial reach of the Green River. Based on the geometric position of the sensor over the bedforms, measurements were taken in the wake, in transitional flow at the bedform crest, and in the internal boundary layer. Spatial distributions of Reynolds shear stress, turbulent kinetic energy, turbulence intensity, and correlation coefficient are qualitatively consistent with those over fixed, two-dimensional bedforms in laboratory flows. Spectral and cospectral analysis demonstrates that energy levels in the lee of the crest (i.e. wake) are two to four times greater than over the crest itself, with minima over the stoss slope (within the developing internal boundary layer). The frequency structure in the wake is sharply defined with single, dominant peaks. Peak and total spectral and cross-spectral energies vary over the bedform in a manner consistent with wave-like perturbations that ,break' or ,roll up' into vortices that amalgamate, grow in size, and eventually diffuse as they are advected downstream. Fluid oscillations in the lee of the dune demonstrate Strouhal similarity between laboratory and field environments, and correspondence between the peak frequencies of these oscillations and the periodicity of surface boils was observed in the field. Copyright © 2005 John Wiley & Sons, Ltd. [source]

    Implications of different degrees of arytenoid cartilage abduction on equine upper airway characteristics

    EQUINE VETERINARY JOURNAL, Issue 7 2008
    V. RAKESH
    Summary Reason for performing study: The necessary degree of arytenoid cartilage abduction (ACA) to restore airway patency at maximal exercise has not been determined. Objectives: Use computational fluid dynamics modelling to measure the effects of different degrees of ACA on upper airway characteristics of horses during exercise. Hypothesis: Maximal ACA by laryngoplasty is necessary to restore normal peak airflow and pressure in Thoroughbred racehorses with laryngeal hemiplegia. Methods: The upper airway was modeled with the left arytenoid in 3 different positions: maximal abduction; 88% cross-sectional area of the rima glottis; and 75% cross-sectional area of the rima glottis. The right arytenoid cartilage was maximally abducted. Two models were assumed: Model 1: no compensation of airway pressures; and Model 2: airway pressure compensation occurs to maintain peak airflow. The cross-sectional pressure and velocity distributions for turbulent flow were studied at peak flow and at different positions along the airway. Results: Model 1: In the absence of a change in driving pressure, 12 and 25% reductions in cross-sectional area of the larynx resulted in 4.11 and 5.65% reductions in peak airflow and 3.68 and 5.64% in tidal volume, respectively, with mild changes in wall pressure. Model 2: To maintain peak flow, a 6.27% increase in driving tracheal pressure was required to compensate for a cross-sectional reduction of 12% and a 13.63% increase in driving tracheal pressure was needed for a cross-sectional area reduction of 25%. This increase in negative driving pressure resulted in regions with low intraluminal and wall pressures, depending on the degree of airway diameter reduction. Conclusion: Assuming no increase in driving pressure, the decrease in left ACA reduced airflow and tidal volume. With increasing driving pressure, a decrease in left ACA changed the wall pressure profile, subjecting the submaximally abducted arytenoid cartilage and adjacent areas to airway collapse. Clinical relevance: The surgical target of ACA resulting in 88% of maximal cross-sectional area seems to be appropriate. [source]

    Development of equine upper airway fluid mechanics model for Thoroughbred racehorses

    EQUINE VETERINARY JOURNAL, Issue 3 2008
    V. RAKESH
    Summary Reason for performing study: Computational fluid dynamics (CFD) models provide the means to evaluate airflow in the upper airways without requiring in vivo experiments. Hypothesis: The physiological conditions of a Thoroughbred racehorse's upper airway during exercise could be simulated. Methods: Computed tomography scanned images of a 3-year-old intact male Thoroughbred racehorse cadaver were used to simulate in vivo geometry. Airway pressure traces from a live Thoroughbred horse, during exercise was used to set the boundary condition. Fluid-flow equations were solved for turbulent flow in the airway during inspiratory and expiratory phases. The wall pressure turbulent kinetic energy and velocity distributions were studied at different cross-sections along the airway. This provided insight into the general flow pattern and helped identify regions susceptible to dynamic collapse. Results: The airflow velocity and static tracheal pressure were comparable to data of horses exercising on a high-speed treadmill reported in recent literature. The cross-sectional area of the fully dilated rima glottidis was 7% greater than the trachea. During inspiration, the area of highest turbulence (i.e. kinetic energy) was in the larynx, the rostral aspect of the nasopharynx was subjected to the most negative wall pressure and the highest airflow velocity is more caudal on the ventral aspect of the nasopharynx (i.e. the soft palate). During exhalation, the area of highest turbulence was in the rostral and mid-nasopharynx, the maximum positive pressure was observed at the caudal aspect of the soft palate and the highest airflow velocity at the front of the nasopharynx. Conclusions and clinical relevance: In the equine upper airway collapsible area, the floor of the rostral aspect of the nasopharynx is subjected to the most significant collapsing pressure with high average turbulent kinetic during inhalation, which may lead to palatal instability and explain the high prevalence of dorsal displacement of the soft palate (DDSP) in racehorses. Maximal abduction of the arytenoid cartilage may not be needed for optimal performance, since the trachea cross-sectional area is 7% smaller than the rima glottidis. [source]

    A spatially advancing turbulent flow and heat transfer in a curved channel

    HEAT TRANSFER - ASIAN RESEARCH (FORMERLY HEAT TRANSFER-JAPANESE RESEARCH), Issue 1 2010
    Koji Matsubara
    Abstract Direct numerical simulation was performed for a spatially advancing turbulent flow and heat transfer in a two-dimensional curved channel, where one wall was heated to a constant temperature and the other wall was cooled to a different constant temperature. In the simulation, fully developed flow and temperature from the straight-channel driver was passed through the inlet of the curved-channel domain. The frictional Reynolds number was assigned 150, and the Prandtl number was given 0.71. Since the flow field was examined in the previous paper, the thermal features are mainly targeted in this paper. The turbulent heat flux showed trends consistent with a growing process of large-scale vortices. In the curved part, the wall-normal component of the turbulent heat flux was twice as large as the counterpart in the straight part, suggesting active heat transport of large-scale vortices. In the inner side of the same section, temperature fluctuation was abnormally large compared with the modest fluctuation of the wall-normal velocity. This was caused by the combined effect of the large-scale motion of the vortices and the wide variation of the mean temperature; in such a temperature distribution, large-scale ejection of the hot fluid near the outer wall, which is transported into the near inner-wall region, should have a large impact on the thermal boundary layer near the inner wall. Wave number decomposition was conducted for various statistics, which showed that the contribution of the large-scale vortex to the total turbulent heat flux normal to the wall reached roughly 80% inside the channel 135° downstream from the curved-channel inlet. © 2009 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/htj.20275 [source]

    Prediction of turbulent flow and heat transfer within rotating U-shaped passages

    HEAT TRANSFER - ASIAN RESEARCH (FORMERLY HEAT TRANSFER-JAPANESE RESEARCH), Issue 6 2006
    Liu Chuan Kai
    Abstract Numerical predictions of three-dimensional flow and heat transfer are presented for rotating serpentine passages with and without rib turbulators. The coolant air is pressurized and its operating conditions are selected closely to match actual turbine operating parameters. Two different arrangements of rib turbulators were studied: (1) transverse ribs on the leading and trailing walls and (2) transverse ribs on all four walls. The rib height-to-hydraulic diameter ratio (e/Dh) is 0.143; the rib pitch-to-height ratio (s/e) is 7. Results for the rib-roughened serpentine passages were compared with those of smooth ones calculated in the literature. It was shown that a significant enhancement is achieved by means of rib turbulators in a serpentine passage at a stationary state as well as in a rotating state. In the radially-outward flow passages, the effect of rotation on heat transfer is relatively prominent. The secondary flows induced by the Coriolis forces are most intensive in the channel with four ribbed surfaces. The heat transfer after a 180° sharp turn in the smooth channel is influenced more by the sharp-turn-induced flow than the rib-roughened ones. © 2006 Wiley Periodicals, Inc. Heat Trans Asian Res, 35(6): 410,420, 2006; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/htj.20125 [source]

    Cover Picture: Fabrication of Multicomponent Microsystems by Directed Three-Dimensional Self-Assembly (Adv. Funct.

    ADVANCED FUNCTIONAL MATERIALS, Issue 5 2005
    Mater.
    Abstract Directed three-dimensional self-assembly to assemble and package integrated semiconductor devices is demonstrated by Jacobs and Zheng on p.,732. The self-assembly process uses geometrical shape recognition to identify different components and surface-tension between liquid solder and metal-coated areas to form mechanical and electrical connections. The components (top left) self-assemble in a turbulent flow (center) and form functional multi-component microsystems (bottom right) by sequentially adding parts to the assembly solution. The technique provides, for the first time, a route to enable the realization of three-dimensional heterogeneous microsystems that contain non-identical parts, and connecting them electrically. We have developed a directed self-assembly process for the fabrication of three-dimensional (3D) microsystems that contain non-identical parts and a statistical model that relates the process yield to the process parameters. The self-assembly process uses geometric-shape recognition to identify different components, and surface tension between liquid solder and metal-coated areas to form mechanical and electrical connections. The concept is used to realize self-packaging microsystems that contain non-identical subunits. To enable the realization of microsystems that contain more than two non-identical subunits, sequential self-assembly is introduced, a process that is similar to the formation of heterodimers, heterotrimers, and higher aggregates found in nature, chemistry, and chemical biology. The self-assembly of three-component assemblies is demonstrated by sequentially adding device segments to the assembly solution including two hundred micrometer-sized light-emitting diodes (LEDs) and complementary metal oxide semiconductor (CMOS) integrated circuits. Six hundred AlGaInP/GaAs LED segments self-assembled onto device carriers in two minutes, without defects, and encapsulation units self-assembled onto the LED-carrier assemblies to form a 3D circuit path to operate the final device. The self-assembly process is a well-defined statistical process. The process follows a first-order, non-linear differential equation. The presented model relates the progression of the self-assembly and yield with the process parameters,component population and capture probability,that are defined by the agitation and the component design. [source]

    Simulation of coherent structures in turbulent boundary layer using Gao,Yong equations of turbulence

    HEAT TRANSFER - ASIAN RESEARCH (FORMERLY HEAT TRANSFER-JAPANESE RESEARCH), Issue 5 2004
    Bo Liu
    Abstract The equations of incompressible turbulent flow developed by the Gao,Yong turbulence model have two important features. First, they do not contain any empirical coefficients or wall functions. Second, the series representation of turbulence energy equation reflects multi-scale structures of the nonlinearity of turbulence, and, therefore, is capable of describing both statistical mean flows and the coherent structures. This paper presents some simulation results of a two-dimensional turbulent boundary layer with zero pressure gradient based on Gao,Yong equations of turbulence. With a staggered grid arrangement, an incompressible SIMPLE code was used in the simulations. The simulated coherent structures have verified the adaptability of the newly derived equations to real intermittent turbulent flows. The effect of the orders of the energy equation and computational grid scales on the detection of coherent structures is also investigated. © 2004 Wiley Periodicals, Inc. Heat Trans Asian Res, 33(5): 287,298, 2004; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/htj.20019 [source]

    Variability of shallow overland flow velocity and soil aggregate transport observed with digital videography

    HYDROLOGICAL PROCESSES, Issue 20 2008
    A. Sidorchuk
    Abstract Field experiments at Tiramoana station 30 km north of Christchurch, New Zealand using an erosion plot 16·5 m long, 0·6 m wide, and with a slope of 14,14·5° on rendzina soil aimed to measure the variability of flow velocity and of soil aggregates transport rate in shallow overland flow. Discharge/cross-section area ratio was used to estimate mean velocity, and high-speed digital video camera and image analysis provided information about flow and sediment transport variability. Six flow runs with 0·5,3·0 L s,1 discharges were supercritical with Froude numbers close to or more than 1. Mean flow velocity followed Poiseuille law, float numbers were more than 1·5 and hydraulic resistance was an inverse proportional function of the Reynolds number, which is typical for laminar flows. Hence actual velocity varied through time significantly and the power spectrum was of ,red-noise', which is typical for turbulent flow. Sediment transport rates had even higher variability, and soil aggregates transport was a compound Poisson process. Copyright © 2008 John Wiley & Sons, Ltd. [source]

    An integrated, finite element-based process model for the analysis of flow, heat transfer, and solidification in a continuous slab caster

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 3 2003
    C. H. Moon
    Abstract An integrated, finite element-based process model is presented for the prediction of full three-dimensional flow, heat transfer, and solidification occurring in a continuous caster. Described in detail are the basic models for the analysis of turbulent flow and heat transfer in the liquid steel zone, in the zone of mixture of the liquid steel and solidified steel, and in the solidified zone. Then, the models are integrated to form a process model which can take into account the strong interdependence between the heat transfer behaviour and the flow behaviour. The capability of the process model to reveal the detailed aspects of turbulent flow, heat transfer, and solidification occurring in a continuous caster is demonstrated through a series of process simulations. Copyright © 2003 John Wiley & Sons, Ltd. [source]

    Implicit algebraic model for predicting turbulent heat flux in film cooling flow

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 5 2010
    Mehran Rajabi-Zargarabadi
    Abstract The present study addresses a new effort to improve the prediction of turbulent heat flux in the film cooling flow by applying the implicit algebraic flux (IAF) model of Rogers et al. A three-dimensional symmetry case is investigated using a film hole length-to-diameter ratio of 1.75 and an injection angle of 35,. The low Reynolds number second moment closure (SMC) model with a wall-reflection term is employed for simulating the turbulent flow field right up to the wall. Results obtained from the IAF model are compared with two other algebraic turbulent heat flux models, namely, the simple eddy diffusivity (SED) with a constant turbulent Prandtl number and the generalized gradient diffusion hypothesis (GGDH). Comparisons of the turbulent heat flux components calculated by these models show that the major difference appears in the streamwise turbulent heat flux. These models demonstrate a significant effect on the prediction of film cooling effectiveness. The SED model with a constant prescribed value for the turbulent Prandtl number fails to predict the cooling air spreading in the lateral direction while by employing the GGDH and IAF models, the spreading of the cooling air and the decay of the effectiveness in the core region are reasonably predicted. A combination of the SMC and IAF models for simulating the turbulent flow and heat transfer is capable of predicting the streamwise and lateral film cooling effectiveness in very good agreement with the available experimental data. Copyright © 2009 John Wiley & Sons, Ltd. [source]

    An evolutionary optimization of diffuser shapes based on CFD simulations

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10 2010
    S. Ghosh
    Abstract An efficient and robust algorithm is presented for the optimum design of plane symmetric diffusers handling incompressible turbulent flow. The indigenously developed algorithm uses the CFD software: Fluent for the hydrodynamic analysis and employs a genetic algorithm (GA) for optimization. For a prescribed inlet velocity and outlet pressure, pressure recovery coefficient C (the objective function) is estimated computationally for various design options. The CFD software and the GA have been combined in a monolithic platform for a fully automated operation using some special control commands. Based on the developed algorithm, an extensive exercise has been made to optimize the diffuser shape. Different methodologies have been adopted to create a large number of design options. Interestingly, not much difference has been noted in the optimum C values obtained through different approaches. However, in all the approaches, a better design has been obtained through a proper selection of the number of design variables. Finally, the effect of diffuser length on the optimum shape has also been studied. Copyright © 2009 John Wiley & Sons, Ltd. [source]

    Two-dimensional prediction of time dependent, turbulent flow around a square cylinder confined in a channel

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2010
    M. Raisee
    Abstract This paper presents two-dimensional and unsteady RANS computations of time dependent, periodic, turbulent flow around a square block. Two turbulence models are used: the Launder,Sharma low-Reynolds number k,, model and a non-linear extension sensitive to the anisotropy of turbulence. The Reynolds number based on the free stream velocity and obstacle side is Re=2.2×104. The present numerical results have been obtained using a finite volume code that solves the governing equations in a vertical plane, located at the lateral mid-point of the channel. The pressure field is obtained with the SIMPLE algorithm. A bounded version of the third-order QUICK scheme is used for the convective terms. Comparisons of the numerical results with the experimental data indicate that a preliminary steady solution of the governing equations using the linear k,, does not lead to correct flow field predictions in the wake region downstream of the square cylinder. Consequently, the time derivatives of dependent variables are included in the transport equations and are discretized using the second-order Crank,Nicolson scheme. The unsteady computations using the linear and non-linear k,, models significantly improve the velocity field predictions. However, the linear k,, shows a number of predictive deficiencies, even in unsteady flow computations, especially in the prediction of the turbulence field. The introduction of a non-linear k,, model brings the two-dimensional unsteady predictions of the time-averaged velocity and turbulence fields and also the predicted values of the global parameters such as the Strouhal number and the drag coefficient to close agreement with the data. Copyright © 2009 John Wiley & Sons, Ltd. [source]

    Turbulence in a three-dimensional wall-bounded shear flow

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 8 2010
    A. Holstad
    Abstract A new turbulent flow with distinct three-dimensional characteristics has been designed in order to study the impact of mean-flow skewing on the turbulent coherent vortices and Reynolds-averaged statistics. The skewing of a unidirectional plane Couette flow was achieved by means of a spanwise pressure gradient. Direct numerical simulations of the statistically steady Couette,Poiseuille flow enabled in-depth explorations of the turbulence field in the skewed flow. The imposition of a modest spanwise gradient turned the mean flow about 8° away from the original Couette flow direction and this turning angle remained nearly the same over the entire cross section. Nevertheless, a substantial non-alignment between the turbulent shear stress angle and the mean velocity gradient angle was observed. The structure parameter turned out to slightly exceed that in the pure Couette flow, contrary to the observations made in some other three-dimensional shear flows. Coherent flow structures, which are known to be associated with the Reynolds shear stress in near-wall regions, were identified by the ,2 -criterion. Instantaneous and ensemble-averaged vortices resembled those found in the unidirectional Couette flow. In the skewed flow, however, the vortex structures were turned to align with the local mean-flow direction. The conventional symmetry between Case 1 and Case 2 vortices was broken due to the mean-flow three-dimensionality. The turning of the coherent vortices and the accompanying symmetry-breaking gave rise to secondary and tertiary turbulent shear stress components. By averaging the already ensemble-averaged shear stresses associated with Case 1 and Case 2 vortices in the homogeneous directions, a direct link between the educed near-wall structures and the Reynolds-averaged turbulent stresses was established. These observations provide evidence in support of the hypothesis that the structural model proposed for two-dimensional turbulent boundary layers remains valid also in flows with moderate mean three-dimensionality. Copyright © 2009 John Wiley & Sons, Ltd. [source]

    FEM simulation of turbulent flow in a turbine blade passage with dynamical fluid,structure interaction

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2009
    Lixiang Zhang
    Abstract Results are described from a combined mathematical modeling and numerical iteration schemes of flow and vibration. We consider the coupling numerical simulations of both turbulent flow and structure vibration induced by flow. The methodology used is based on the stabilized finite element formulations with time integration. A fully coupled model of flow and flow-induced structure vibration was established using a hydride generalized variational principle of fluid and solid dynamics. The spatial discretization of this coupling model is based on the finite element interpolating formulations for the fluid and solid structure, while the different time integration schemes are respectively used for fluid and solid structure to obtain a stabilized algorithm. For fluid and solid dynamics, Hughes' predictor multi-corrector algorithm and the Newmark method are monolithically used to realize a monolithic solution of the fully coupled model. The numerical convergence is ensured for small deformation vibrating problems of the structure by using different time steps for fluid and solid, respectively. The established model and the associated numerical methodology developed in the paper were then applied to simulate two different flows. The first one is the lid-driven square cavity flow with different Reynolds numbers of 1000, 400 and 100 and the second is the turbulent flows in a 3-D turbine blade passage with dynamical fluid,structure interaction. Good agreement between numerical simulations and measurements of pressure and vibration acceleration indicates that the finite element method formulations developed in this paper are appropriate to deal with the flow under investigation. Copyright © 2009 John Wiley & Sons, Ltd. [source]

    Simulating turbulent Dean flow in Cartesian coordinates

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2009
    George K. El Khoury
    Abstract A simplified approach to simulate turbulent flows in curved channels is proposed. A set of governing equations of motion in Cartesian coordinates is derived from the full Navier,Stokes equations in cylindrical coordinates. Terms to first order in the dimensionless curvature parameter are retained, whereas higher-order terms are neglected. The curvature terms are implemented in a conventional Navier,Stokes code using Cartesian coordinates. Direct numerical simulations (DNS) of turbulent flow in weakly curved channels are performed. The pronounced asymmetries in the mean flow and the turbulence statistics observed in earlier DNS studies are faithfully reproduced by the present simplified Navier,Stokes model. It is particularly rewarding that also distinct pairs of counter-rotating streamwise-oriented vortices are embedded in the simulated flow field. Copyright © 2008 John Wiley & Sons, Ltd. [source]

    A gradient smoothing method (GSM) for fluid dynamics problems

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10 2008
    G. R. Liu
    Abstract A novel gradient smoothing method (GSM) based on irregular cells and strong form of governing equations is presented for fluid dynamics problems with arbitrary geometries. Upon the analyses about the compactness and the positivity of coefficients of influence of their stencils for approximating a derivative, four favorable schemes (II, VI, VII and VIII) with second-order accuracy are selected among the total eight proposed discretization schemes. These four schemes are successively verified and carefully examined in solving Poisson's equations, subjected to changes in the number of nodes, the shapes of cells and the irregularity of triangular cells, respectively. Numerical results imply us that all the four schemes give very good results: Schemes VI and VIII produce a slightly better accuracy than the other two schemes on irregular cells, but at a higher cost in computation. Schemes VII and VIII that consistently rely on gradient smoothing operations are more accurate than Schemes II and VI in which directional correction is imposed. It is interestingly found that GSM is insensitive to the irregularity of meshes, indicating the robustness of the presented GSM. Among the four schemes of GSM, Scheme VII outperforms the other three schemes, for its outstanding overall performance in terms of numerical accuracy, stability and efficiency. Finally, GSM solutions with Scheme VII to some benchmarked compressible flows including inviscid flow over NACA0012 airfoil, laminar flow over flat plate and turbulent flow over an RAE2822 airfoil are presented, respectively. Copyright © 2008 John Wiley & Sons, Ltd. [source]

    Numerical model for the prediction of dilute, three-dimensional, turbulent fluid,particle flows, using a Lagrangian approach for particle tracking and a CVFEM for the carrier phase

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 5 2008
    L. A. Oliveira
    Abstract A numerical model for dilute, three-dimensional, turbulent, incompressible fluid,solid particle flows and its application to a demonstration problem are presented. An Eulerian description is used to model the flow of the fluid (carrier) phase, and the governing equations are solved using a control-volume finite element method (CVFEM). The motion of the solid (particulate) phase is simulated using a Lagrangian approach. An efficient algorithm is proposed for locating the particles in the finite element mesh. In the demonstration problem, which involves a particle-laden axisymmetric jet, a modified k,, turbulence model is used to characterize the velocity and length scales of the turbulent flow of the fluid phase. The effect of turbulence on the particle trajectories is accounted for through a stochastic model. The effect of the particles on the fluid time,mean velocity and turbulence (two-way coupling) is also addressed. Comparisons between predictions and available experimental data for the demonstration problem are presented. Satisfactory agreement is obtained. Copyright © 2008 John Wiley & Sons, Ltd. [source]

    Numerical analysis of turbulent flow separation in a rectangular duct with a sharp 180-degree turn by algebraic Reynolds stress model

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2008
    Hitoshi Sugiyama
    Abstract Turbulent flow in a rectangular duct with a sharp 180-degree turn is difficult to predict numerically because the flow behavior is influenced by several types of forces, including centrifugal force, pressure-driven force, and shear stress generated by anisotropic turbulence. In particular, this type of flow is characterized by a large-scale separated flow, and it is difficult to predict the reattachment point of a separated flow. Numerical analysis has been performed for a turbulent flow in a rectangular duct with a sharp 180-degree turn using the algebraic Reynolds stress model. A boundary-fitted coordinate system is introduced as a method for coordinate transformation to set the boundary conditions next to complicated shapes. The calculated results are compared with the experimental data, as measured by a laser-Doppler anemometer, in order to examine the validity of the proposed numerical method and turbulent model. In addition, the possibility of improving the wall function method in the separated flow region is examined by replacing the log-law velocity profile for a smooth wall with that for a rough wall. The analysis results indicated that the proposed algebraic Reynolds stress model can be used to reasonably predict the turbulent flow in a rectangular duct with a sharp 180-degree turn. In particular, the calculated reattachment point of a separated flow, which is difficult to predict in a turbulent flow, agrees well with the experimental results. In addition, the calculation results suggest that the wall function method using the log-law velocity profile for a rough wall over a separated flow region has some potential for improving the prediction accuracy. Copyright © 2007 John Wiley & Sons, Ltd. [source]

    Coupled ghost fluid/two-phase level set method for curvilinear body-fitted grids

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2007
    Juntao Huang
    Abstract A coupled ghost fluid/two-phase level set method to simulate air/water turbulent flow for complex geometries using curvilinear body-fitted grids is presented. The proposed method is intended to treat ship hydrodynamics problems. The original level set method for moving interface flows was based on Heaviside functions to smooth all fluid properties across the interface. We call this the Heaviside function method (HFM). The HFM requires fine grids across the interface. The ghost fluid method (GFM) has been designed to explicitly enforce the interfacial jump conditions, but the implementation of the jump conditions in curvilinear grids is intricate. To overcome these difficulties a coupled GFM/HFM method was developed in which approximate jump conditions are derived for piezometric pressure and velocity and pressure gradients based on exact continuous velocity and stress and jump in momentum conditions with the jump in density maintained but continuity of the molecular and turbulent viscosities imposed. The implementation of the ghost points is such that no duplication of memory storage is necessary. The level set method is adopted to locate the air/water interface, and a fast marching method was implemented in curvilinear grids to reinitialize the level set function. Validations are performed for three tests: super- and sub-critical flow without wave breaking and an impulsive plunging wave breaking over 2D submerged bumps, and the flow around surface combatant model DTMB 5512. Comparisons are made against experimental data, HFM and single-phase level set computations. The proposed method performed very well and shows great potential to treat complicated turbulent flows related to ship flows. Copyright © 2007 John Wiley & Sons, Ltd. [source]

    A modification of the artificial compressibility algorithm with improved convergence characteristics

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 4 2007
    Frank Muldoon
    Abstract The artificial compressibility algorithm has a significant drawback in the difficulty of choosing the artificial compressibility parameter, improper choice of which leads either to slow convergence or divergence. A simple modification of the equation for pressure in the artificial compressibility algorithm which removes the difficulty of choosing the artificial compressibility parameter is proposed. It is shown that the choice of the relaxation parameters for the new algorithm is relatively straightforward, and that the same values can be used to provide robust convergence for a range of application problems. This new algorithm is easily parallelized making it suitable for computations such as direct numerical simulation (DNS) which require the use of distributed memory machines. Two key benchmark problems are studied in evaluating the new algorithm: DNS of a fully developed turbulent channel flow, and DNS of a driven-cavity flow, using both explicit and implicit time integration schemes. The new algorithm is also validated for a more complex flow configuration of turbulent flow over a backward-facing step, and the computed results are shown to be in good agreement with experimental data and previous DNS work. Copyright © 2007 John Wiley & Sons, Ltd. [source]

    Numerical calculations of erosion in an abrupt pipe contraction of different contraction ratios

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 1 2004
    M. A. Habib
    Abstract Erosion predictions in a pipe with abrupt contraction of different contraction ratios for the special case of two-phase (liquid and solid) turbulent flow with low particle concentration are presented. A mathematical model based on the time-averaged governing equations of 2-D axi-symmetric turbulent flow is used for the calculations of the fluid velocity field (continuous phase). The particle-tracking model of the solid particles is based on the solution of the governing equation of each particle motion taking into consideration the effect of particle rebound behaviour. Models of erosion were used to predict the erosion rate in mg/g. The effect of Reynolds number and flow direction with respect to the gravity was investigated for three contraction geometries considering water flow in a carbon steel pipe. The results show that the influence of the contraction ratio on local erosion is very significant. However, this influence becomes insignificant when the average erosion rates over the sudden contraction area are considered. The results also indicate the significant influence of inlet velocity variations. The influence of buoyancy is significant for the cases of low velocity of the continuous flow. A threshold velocity below which erosion may be neglected was indicated. Copyright © 2004 John Wiley & Sons, Ltd. [source]

    Large eddy simulation of turbulent concentric annular channel flows

    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2004
    Nan-Sheng Liu
    Abstract Fully developed turbulent concentric annular channel flow has been investigated numerically by use of large eddy simulation (LES) technique coupled with a localized one-equation dynamic subgrid-scale (SGS) model. The objective of this study is to deal with the behaviour of turbulent flow near the inner and outer walls of the concentric annular channel and to examine the effectiveness of LES technique for predicting the turbulent flow influenced by the transverse curvature effect. The computations are performed for the Reynolds number Re,=180, 395 and 640, based on an averaged friction velocity and the annular channel width with the inner and outer cylinder radius being Ri=1 and Ro=2. To validate the present approach, calculated results for turbulent pipe flow and concentric annular channel flow are compared with available experimental data and direct numerical simulation results, which confirms that the present approach can be used to study turbulent concentric annular channel flow satisfactorily. To elucidate turbulence characteristics in the concentric annular channel, some typical quantities, including the resolved velocity, turbulence intensity, turbulent eddy viscosity, SGS kinetic energy, SGS dissipation rate, Reynolds stress budgets, and turbulence structures based on the velocity fluctuations, are analysed. Copyright © 2004 John Wiley & Sons, Ltd. [source]