Home About us Contact
|
Home About us Contact | ||
|
|
||
Turbulence Models (turbulence + models)
Kinds of Turbulence Models Selected AbstractsTwo-dimensional prediction of time dependent, turbulent flow around a square cylinder confined in a channelINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2010M. 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] On the construction of manufactured solutions for one and two-equation eddy-viscosity modelsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 2 2007L. Eça Abstract This paper presents manufactured solutions (MSs) for some well-known eddy-viscosity turbulence models, viz. the Spalart & Allmaras one-equation model and the TNT and BSL versions of the two-equation k,, model. The manufactured flow solutions apply to two-dimensional, steady, wall-bounded, incompressible, turbulent flows. The two velocity components and the pressure are identical for all MSs, but various alternatives are considered for specifying the eddy-viscosity and other turbulence quantities in the turbulence models. The results obtained for the proposed MSs with a second-order accurate numerical method show that the MSs for turbulence quantities must be constructed carefully to avoid instabilities in the numerical solutions. This behaviour is model dependent: the performance of the Spalart & Allmaras and k,, models is significantly affected by the type of MS. In one of the MSs tested, even the two versions of the k,, model exhibit significant differences in the convergence properties. Copyright © 2006 John Wiley & Sons, Ltd. [source] Using vorticity to define conditions at multiple open boundaries for simulating flow in a simplified vortex settling basinINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 1 2007A. N. Ziaei Abstract In this paper a method is developed to define multiple open boundary (OB) conditions in a simplified vortex settling basin (VSB). In this method, the normal component of the momentum equation is solved at the OBs, and tangential components of vorticity are calculated by solving vorticity transport equations only at the OBs. Then the tangential vorticity components are used to construct Neumann boundary conditions for tangential velocity components. Pressure is set to its ambient value, and the divergence-free condition is satisfied at these boundaries by employing the divergence as the Neumann condition for the normal-direction momentum equation. The 3-D incompressible Navier,Stokes equations in a primitive-variable form are solved using the SIMPLE algorithm. Grid-function convergence tests are utilized to verify the numerical results. The complicated laminar flow structure in the VSB is investigated, and preliminary assessment of two popular turbulence models, k,, and k,,, is conducted. Copyright © 2006 John Wiley & Sons, Ltd. [source] Comparative study of the continuous phase flow in a cyclone separator using different turbulence models,INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2005H. Shalaby Abstract Numerical calculations were carried out at the apex cone and various axial positions of a gas cyclone separator for industrial applications. Two different NS-solvers (a commercial one (CFX 4.4 ANSYS GmbH, Munich, Germany, CFX Solver Documentation, 1998), and a research code (Post-doctoral Thesis, Technical University of Chemnitz, Germany, September, 2002)) based on a pressure correction algorithm of the SIMPLE method have been applied to predict the flow behaviour. The flow was assumed as unsteady, incompressible and isothermal. A k,, turbulence model has been applied first using the commercial code to investigate the gas flow. Due to the nature of cyclone flows, which exhibit highly curved streamlines and anisotropic turbulence, advanced turbulence models such as Reynolds stress model (RSM) and large eddy simulation (LES) have been used as well. The RSM simulation was performed using the commercial package activating the Launder et al.'s (J. Fluid. Mech. 1975; 68(3):537,566) approach, while for the LES calculations the research code has been applied utilizing the Smagorinsky model. It was found that the k,, model cannot predict flow phenomena inside the cyclone properly due to the strong curvature of the streamlines. The RSM results are comparable with LES results in the area of the apex cone plane. However, the application of the LES reveals qualitative agreement with the experimental data, but requires higher computer capacity and longer running times than RSM. This paper is organized into five sections. The first section consists of an introduction and a summary of previous work. Section 2 deals with turbulence modelling including the governing equations and the three turbulence models used. In Section 3, computational parameters are discussed such as computational grids, boundary conditions and the solution algorithm with respect to the use of MISTRAL/PartFlow-3D. In Section 4, prediction profiles of the gas flow at axial and apex cone positions are presented and discussed. Section 5 summarizes and concludes the paper. Copyright © 2005 John Wiley & Sons, Ltd. [source] Anisotropic, isothermal, turbulent swirling flow in a complex combustor geometryINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10-11 2005L. N. Jones Abstract The performance of popular second moment closure (LRR, SSG) turbulence models is assessed and compared against experimental data for anisotropic swirling flow in a cylindrical combustion chamber. In contrast to previous studies, where the dissipation anisotropy is correlated with the stress anisotropy, the benefit of approximating the former for swirling flows in terms of the mean strain and vorticity is investigated. Second moment closure models are found to predict mean and turbulent flow quantities reasonably well everywhere except near the wall. The anisotropic dissipation model is found to improve prediction of mean flow quantities near the chamber axis and acts to preserve turbulence further downstream. Copyright © 2005 John Wiley & Sons, Ltd. [source] Turbulence model and numerical scheme assessment for buffet computationsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2004Eric Goncalves Abstract The prediction of shock-induced oscillations over transonic rigid airfoils is important for a better understanding of the buffeting phenomenon. The unsteady resolution of the Navier,Stokes equations is performed with various transport-equation turbulence models in which corrections are added for non-equilibrium flows. The lack of numerical efficiency due to the CFL stability condition is circumvented by the use of a wall law approach and a dual time stepping method. Moreover, various numerical schemes are used to try and be independent of the numerical discretization. Comparisons are made with the experimental results obtained for the supercritical RA16SC1 airfoil. They show the interest in using the SST correction or realizability conditions to get correct predictions of the frequency, amplitude and pressure fluctuations over the airfoil. Copyright © 2004 John Wiley & Sons, Ltd. [source] Assessment of two-equation turbulence modelling for high Reynolds number hydrofoil flowsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2004N. Mulvany Abstract This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two-equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k,,model, the realizable k,,model, the standard k,,model and the shear-stress-transport (SST) k,,model. It has found that the realizable k,,turbulence model used with enhanced wall functions and near-wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. Copyright © 2004 John Wiley & Sons, Ltd. [source] Computation of strongly swirling confined flows with cubic eddy-viscosity turbulence modelsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2003Xiaodong Yang Graduate Student Abstract An investigation on the predictive performance of four cubic eddy-viscosity turbulence models for two strongly swirling confined flows is presented. Comparisons of the prediction with the experiments show clearly the superiority of cubic models over the linear k,,model. The linear k,,model does not contain any mechanism to describe the stabilizing effects of swirling motion and as a consequence it performs poorly. Cubic models return a lower level of Reynolds stresses and the combined forced-free vortex profiles of tangential velocity close to the measurements in response to the interaction between swirl-induced curvature and stresses. However, a fully developed rotating pipe flow is too simple to contain enough flow physics, so the calibration of cubic terms is still a topic of investigation. It is shown that explicit algebraic stress models require fewer calibrations and contain more flow physics. Copyright © 2003 John Wiley & Sons, Ltd. [source] Computation of an unsteady complex geometry flow using novel non-linear turbulence modelsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2003Paul G. Tucker Abstract Non-linear zonal turbulence models are applied to an unsteady complex geometry flow. These are generally found to marginally improve predicted turbulence intensities. However, relative to linear models, convergence is mostly difficult to achieve. Clipping of some non-linear Reynolds stress components is required along with velocity field smoothing or alternative measures. Smoothing is naturally achieved through multilevel convergence restriction operators. As a result of convergence difficulties, generally, non-linear model computational costs detract from accuracy gains. For standard Reynolds stress model results, again computational costs are prohibitive. Also, mean velocity profile data accuracies are found worse than for a simple mixing length model. Of the non-linear models considered, the explicit algebraic stress showed greatest promise with respect to accuracy and stability. However, even this shows around a 30% error in total (the sum of turbulence and unsteadiness) intensity. In strong contradiction to measurements the non-linear and Reynolds models predict quasi-steady flows. This is probably a key reason for the total intensity under-predictions. Use of LES in a non-linear model context might help remedy this modelling aspect. Copyright © 2003 John Wiley & Sons, Ltd. [source] URANS computations for an oscillatory non-isothermal triple-jet using the k,, and second moment closure turbulence modelsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2003M. Nishimura Abstract Low Reynolds number turbulence stress and heat flux equation models (LRSFM) have been developed to enhance predictive capabilities. A new method is proposed for providing the wall boundary condition for dissipation rate of turbulent kinetic energy, ,, to improve the model capability upon application of coarse meshes for practical use. The proposed method shows good agreement with accepted correlations and experimental data for flows with various Reynolds and Prandtl numbers including transitional regimes. Also, a mesh width about 5 times or larger than that used in existing models is applicable by using the present boundary condition. The present method thus enhanced computational efficiency in applying the complex turbulence model, LRSFM, to predictions of complicated flows. Unsteady Reynolds averaged Navier,Stokes (URANS) computations are conducted for an oscillatory non-isothermal quasi-planar triple-jet. Comparisons are made between an experiment and predictions with the LRSFM and the standard k,, model. A water test facility with three vertical jets, the cold in between two hot jets, simulates temperature fluctuations anticipated at the outlet of a liquid metal fast reactor core. The LRSFM shows good agreement with the experiment, with respect to mean profiles and the oscillatory motion of the flow, while the k,, model under-predicts the mixing due to the oscillation, such that a transverse mean temperature difference remains far downstream. Copyright © 2003 John Wiley & Sons, Ltd. [source] Computation of turbulent free-surface flows around modern shipsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 4 2003Tingqiu Li Abstract This paper presents the calculated results for three classes of typical modern ships in modelling of ship-generated waves. Simulations of turbulent free-surface flows around ships are performed in a numerical water tank, based on the FINFLO-RANS SHIP solver developed at Helsinki University of Technology. The Reynolds-averaged Navier,Stokes (RANS) equations with the artificial compressibility and the non-linear free-surface boundary conditions are discretized by means of a cell-centred finite-volume scheme. The convergence performance is improved with the multigrid method. A free surface is tracked using a moving mesh technology, in which the non-linear free-surface boundary conditions are given on the actual location of the free surface. Test cases recommended are a container ship, a US Navy combatant and a tanker. The calculated results are compared with the experimental data available in the literature in terms of the wave profiles, wave pattern, and turbulent flow fields for two turbulence models, Chien's low Reynolds number k,,model and Baldwin,Lomax's model. Furthermore, the convergence performance, the grid refinement study and the effect of turbulence models on the waves have been investigated. Additionally, comparison of two types of the dynamic free-surface boundary conditions is made. Copyright © 2003 John Wiley& Sons, Ltd. [source] Evaluation of one- and two-equation low- Re turbulence models.INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2003Axisymmetric separating, Part I, swirling flows Abstract This first segment of the two-part paper systematically examines several turbulence models in the context of three flows, namely a simple flat-plate turbulent boundary layer, an axisymmetric separating flow, and a swirling flow. The test cases are chosen on the basis of availability of high-quality and detailed experimental data. The tested turbulence models are integrated to solid surfaces and consist of: Rodi's two-layer k,, model, Chien's low-Reynolds number k,, model, Wilcox's k,, model, Menter's two-equation shear-stress-transport model, and the one-equation model of Spalart and Allmaras. The objective of the study is to establish the prediction accuracy of these turbulence models with respect to axisymmetric separating flows, and flows of high streamline curvature. At the same time, the study establishes the minimum spatial resolution requirements for each of these turbulence closures, and identifies the proper low-Mach-number preconditioning and artificial diffusion settings of a Reynolds-averaged Navier,Stokes algorithm for optimum rate of convergence and minimum adverse impact on prediction accuracy. Copyright © 2003 John Wiley & Sons, Ltd. [source] Evaluation of one- and two-equation low- Re turbulence models.INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2003Part II, Vortex-generator jet, diffusing S-duct flows Abstract This second segment of the two-part paper systematically examines several turbulence models in the context of two flows, namely a vortex flow created by an inclined jet in crossflow, and the flow field in a diffusing S-shaped duct. The test cases are chosen on the basis of availability of high-quality and detailed experimental data. The tested turbulence models are integrated to solid surfaces and consist of: Rodi's two-layer k,, model, Wilcox's k,, model, Menter's two-equation shear,stress-transport model, and the one-equation model of Spalart and Allmaras. The objective of the study is to establish the prediction accuracy of these turbulence models with respect to three-dimensional separated flows with streamline curvature. At the same time, the study establishes the minimum spatial resolution requirements for each of these turbulence closures, and identifies the proper low-Mach-number preconditioning and artificial diffusion settings of a Reynolds-averaged Navier,Stokes algorithm for optimum rate of convergence and minimum adverse impact on prediction accuracy. Copyright © 2003 John Wiley & Sons, Ltd. [source] Calculation of turbulent fluid flow and heat transfer in ducts by a full Reynolds stress modelINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 2 2003Masoud Rokni Abstract A computational method has been developed to predict the turbulent Reynolds stresses and turbulent heat fluxes in ducts by different turbulence models. The turbulent Reynolds stresses and other turbulent flow quantities are predicted with a full Reynolds stress model (RSM). The turbulent heat fluxes are modelled by a SED concept, the GGDH and the WET methods. Two wall functions are used, one for the velocity field and one for the temperature field. All the models are implemented for an arbitrary three-dimensional channel. Fully developed condition is achieved by imposing cyclic boundary conditions in the main flow direction. The numerical approach is based on the finite volume technique with a non-staggered grid arrangement. The pressure,velocity coupling is handled by using the SIMPLEC-algorithm. The convective terms are treated by the van Leer scheme while the diffusive terms are handled by the central-difference scheme. The hybrid scheme is used for solving the , equation. The secondary flow generation using the RSM model is compared with a non-linear k,, model (non-linear eddy viscosity model). The overall comparison between the models is presented in terms of the friction factor and Nusselt number. Copyright © 2003 John Wiley & Sons, Ltd. [source] Convergence acceleration by self-adjusted time stepsize using Bi-CGSTAB method for turbulent separated flow computationINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 2 2002W. B. Tsai Abstract Poor convergence behavior is usually encountered when numerical computations on turbulent separated flow are performed. A design of self-adjusted stepsize concept both in time span and spatial coordinate systems to achieve faster convergence is demonstrated in this study. The determination of the time stepsize based on the concept of minimization of residuals using the Bi-CGSTAB algorithm is proposed. The numerical results show that the time stepsize adjusted by the proposed method indeed improves the convergence rate for turbulent separated flow computations using advanced turbulence models in low-Reynolds number forms. Copyright © 2002 John Wiley & Sons, Ltd. [source] Hydrogen as burner fuel: modelling of hydrogen,hydrocarbon composite fuel combustion and NOx formation in a small burnerINTERNATIONAL JOURNAL OF ENERGY RESEARCH, Issue 11 2005Mustafa Ilbas Abstract The objective of this work is to investigate numerically the turbulent non-premixed hydrogen (H2) and hydrogen,hydrocarbon flames in a small burner. Numerical studies using Fluent code were carried out for air-staged and non-staged cases. The effects of fuel composition from pure hydrogen to natural gas (100%H2, 70%H2+30%CH4, 10%H2+90%CH4, and 100%CH4) were also investigated. The predictions are validated and compared against the experimental results previously obtained and results from the literature. Turbulent diffusion flames are investigated numerically using a finite volume method for the solution of the conservation equations and reaction equations governing the problem. Although, three different turbulence models were tested, the standard k,, model was used for the modelling of the turbulence phenomena in the burner. The temperature and major pollutant concentrations (CO and NOx) distributions are in good agreement with the existing experimental results. Air staging causes rich and lean combustion regions thus lower NOx emissions through the combustor exit. Blending hydrogen with methane causes considerable reduction in temperature levels and thus NO emissions. Increasing the mixture ratio from stoichiometric to leaner mixtures also decreases the temperature and thus NO emissions. Hydrogen may be considered a good alternative fuel for burners, as its use reduces the emission of pollutants, and as it is a renewable synthetic fuel. Copyright © 2005 John Wiley & Sons, Ltd. [source] CFD modeling of heat transfer in turbulent pipe flowsAICHE JOURNAL, Issue 9 2000S. S. Thakre Twelve versions of low Reynolds number k-, and two low Reynolds number Reynolds stress turbulence models for heat transfer were analyzed comparatively. Predictions of the mean axial temperature, the radial and axial turbulent heat fluxes, and the effect of Prandtl number on Nusselt number were compared with the experimental data. The model by Lai and So from the k-, group and Lai and So from the Reynolds stress group had the best overall predictive ability for heat transfer in turbulent pipe flow. The Lai and So model was attributed to its success in the predictions of flow parameters such as mean axial velocity, turbulent kinetic energy, eddy diffusivity, and the overall energy dissipation rate. The k-, models performed relatively better than the Reynolds stress models for predicting the mean axial temperature and the Nusselt number. This qualitative and quantitative study found the need for more sophisticated near-wall experimental measurements and the accuracy of the dissipation (of turbulent energy) and the pressure-scrambling models. [source] An enhanced-physics-based scheme for the NS-, turbulence modelNUMERICAL METHODS FOR PARTIAL DIFFERENTIAL EQUATIONS, Issue 6 2010William W. Miles Abstract We study a new enhanced-physics-based numerical scheme for the NS-alpha turbulence model that conserves both energy and helicity. Although most turbulence models (in the continuous case) conserve only energy, NS-alpha is one of only a very few that also conserve helicity. This is one reason why it is becoming accepted as the most physically accurate turbulence model. However, no numerical scheme for NS-alpha, until now, conserved both energy and helicity, and thus the advantage gained in physical accuracy by modeling with NS-alpha could be lost in a computation. This report presents a finite element numerical scheme, and gives a rigorous analysis of its conservation properties, stability, solution existence, and convergence. A key feature of the analysis is the identification of the discrete energy and energy dissipation norms, and proofs that these norms are equivalent (provided a careful choice of filtering radius) in the discrete space to the usual energy and energy dissipation norms. Numerical experiments are given to demonstrate the effectiveness of the scheme over usual (helicity-ignoring) schemes. A generalization of this scheme to a family of high-order NS-alpha-deconvolution models, which combine the attractive physical properties of NS-alpha with the high accuracy gained by combining ,-filtering with van Cittert approximate deconvolution. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010 [source] Convergence acceleration by varying time-step size using Bi-CGSTAB method for turbulent flow computationNUMERICAL METHODS FOR PARTIAL DIFFERENTIAL EQUATIONS, Issue 5 2001W. B. Tsai Abstract A design of varying step size approach both in time span and spatial coordinate systems to achieve fast convergence is demonstrated in this study. This method is based on the concept of minimization of residuals by the Bi-CGSTAB algorithm, so that the convergence can be enforced by varying the time-step size. The numerical results show that the time-step size determined by the proposed method improves the convergence rate for turbulent computations using advanced turbulence models in low Reynolds-number form, and the degree of improvement increases with the degree of the complexity of the turbulence models. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 454,474, 2001. [source] A low reynolds number k-, modelling of turbulent pipe flow: Flow pattern and energy balanceTHE CANADIAN JOURNAL OF CHEMICAL ENGINEERING, Issue 2 2001Shirish S. Thakre Abstract The present paper addresses a comparative analysis of twelve different versions of low Reynolds number k -, turbulence models. The predictive capability of the models have been tested on the basis of the flow patterns and energy balance. Numerical simulations were performed at the Reynolds numbers of 7400, 22 000 and 500 000. The predicted mean axial velocity and turbulent kinetic energy were compared with the experimental data of Durst et al. (1995) and Schildknecht et al.(1979) for the Reynolds number of 7400 and 22000 respectively. The overall energy balance was established at three Reynolds numbers of 7400, 22 000 and 500000. A comparison of all the models has been predicted. On décrit dans cet article une analyse comparative de douze versions différentes de modèles de turbulence à faibles nombres de Reynolds k -,. La capacité de prédiction de ces modèles a été testée d'après les profils d'écoulement et le bilan énergétique. Des simulations numériques ont été réalisées à des nombres de Reynolds de 7400, 22 000 et 500 000. La vitesse axiale et l'énergie cinétique turbulente moyennes prédites ont été comparées aux données expérimentales de Durst et al. (1995) et Schildknecht et al. (1979) pour les nombres de Reynolds de 7400 et 22 000, respectivement. Le bilan énergétique global a été établi pour les trois nombres de Reynolds. Une comparaison de tous les modèles a été effectuée. [source] Experimental and Numerical Simulation Study of Heat Transfer Due to Confined Impinging Circular JetCHEMICAL ENGINEERING & TECHNOLOGY (CET), Issue 10 2007L. Chang-geng Abstract An experimental and numerical simulation study of heat transfer due to a confined impinging circular jet is presented. In this research, a stainless steel foil heated disk was used as the heat transfer surface of a simulated chip, and the thermocouples were mounted symmetrically along the diameter of the foil to measure the temperature distribution on the surface. Driven by a small pump, a circular air jet (1.5,mm and 1,mm in diameter) impinged on the heat-transfer surface with middle and low Reynolds numbers. The parameters, such as Reynolds number and ratio of height-to-diameter, were changed to investigate the radial distribution of the Nusselt number and the characteristics of heat transfer in the stagnation region. Numerical computations were performed by using several different turbulence models. In wall bounded turbulent flows, near-wall modeling is crucial. Therefore, the turbulence models enhanced wall treatment, such as the RNG ,-, model, may be superior for modeling impingement flows. The numerical results showed reasonable agreement with the experimental data for local heat transfer coefficient distributions. The impinging jet may be an effective method to solve the cooling problem of high power density electronic packaging. [source] Simulation of Turbulent Flow in a Packed BedCHEMICAL ENGINEERING & TECHNOLOGY (CET), Issue 5 2006B. 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] | |||||||||||||