Atmospheric Environment
|
|
- Debra Lawrence
- 6 years ago
- Views:
Transcription
1 Atmospheric Environment 43 (2009) Contents lists available at ScienceDirect Atmospheric Environment journal homepage: CFD modelling of small particle dispersion: The influence of the turbulence kinetic energy in the atmospheric boundary layer C. Gorlé a,b, *, J. van Beeck b, P. Rambaud b, G. Van Tendeloo a a University of Antwerp, Department of Physics, EMAT, Groenenborgerlaan 171, 2020 Antwerp, Belgium b Von Karman Institute for Fluid Dynamics, Waterloosesteenweg 72, 1640 Sint-Genesius-Rode, Belgium article info abstract Article history: Received 2 July 2008 Received in revised form 26 September 2008 Accepted 26 September 2008 Keywords: Computational Fluid Dynamics (CFD) Numerical simulation Atmospheric Boundary Layer (ABL) Turbulence kinetic energy Particle dispersion When considering the modelling of small particle dispersion in the lower part of the Atmospheric Boundary Layer (ABL) using Reynolds Averaged Navier Stokes simulations, the particle paths depend on the velocity profile and on the turbulence kinetic energy, from which the fluctuating velocity components are derived to predict turbulent dispersion. It is therefore important to correctly reproduce the ABL, both for the velocity profile and the turbulence kinetic energy profile. For RANS simulations with the standard k e model, Richards and Hoxey (1993. Appropriate boundary conditions for computational wind engineering models using the k e turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 46 47, ) proposed a set of boundary conditions which result in horizontally homogeneous profiles. The drawback of this method is that it assumes a constant profile of turbulence kinetic energy, which is not always consistent with field or wind tunnel measurements. Therefore, a method was developed which allows the modelling of a horizontally homogeneous turbulence kinetic energy profile that is varying with height. By comparing simulations performed with the proposed method to simulations performed with the boundary conditions described by Richards and Hoxey (1993. Appropriate boundary conditions for computational wind engineering models using the k e turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 46 47, ), the influence of the turbulence kinetic energy on the dispersion of small particles over flat terrain is quantified. Ó 2008 Elsevier Ltd. All rights reserved. 1. Introduction Because of the increasing interest in urban air quality modelling, the application of Computational Fluid Dynamics (CFD) to study atmospheric dispersion processes in the lower part of the Atmospheric Boundary Layer (ABL) has become an important research subject. Validation is an essential aspect of this research and several comparative studies between CFD and wind tunnel or field measurements have been performed, e.g. Leitl et al. (1997), Meroney et al. (1999), Garcia Sagrado et al. (2002), Chang and Meroney (2003), Dixon et al. (2006), and Wang and McNamara (2007). In all of these publications the intermittent nature of the dispersion process in the wind tunnel and field measurements as opposed to the Reynolds Averaged Navier Stokes (RANS) solution of the CFD simulations is indicated as a reason for the observed discrepancies. These discrepancies could be reduced by either averaging measured values over a sufficiently long period or by performing Large Eddy * Corresponding author. University of Antwerp, Department of Physics, EMAT, Groenenborgerlaan 171, 2020 Antwerp, Belgium. Tel.: þ ; fax: þ address: gorle@vki.ac.be (C. Gorlé). Simulations (LES), which resolve the larger scales of turbulence Schatzmann and Leitl (2002). Other significant sources of error reported are inaccuracies in the boundary conditions for the flow or the pollutant source and the underestimation of the turbulence kinetic energy (Garcia Sagrado et al., 2002; Dixon et al., 2006). It is obvious that a correct reproduction of the ABL from the wind tunnel or field experiments is essential to obtain accurate simulation results for atmospheric dispersion. This is usually achieved by imposing fully developed inlet profiles for the velocity and turbulence characteristics that represent the influence of upstream roughness elements, which are not included in the domain. These inlet profiles should be identical to the measurement conditions and should also be established at the locations of interest further downstream in the computational domain. However, a frequently reported problem is that the inlet profiles are not horizontally homogeneous throughout the domain (Franke et al., 2004, 2007; Blocken et al., 2007a,b; Riddle et al., 2004; Hargreaves and Wright, 2007). The reason for horizontal variations in the profiles is an incompatibility of the inlet profiles with the applied wall functions. As a result the profiles will quickly adapt to the implemented wall functions, leading to unacceptable changes further downstream in the domain. It has been shown that both the velocity and /$ see front matter Ó 2008 Elsevier Ltd. All rights reserved. doi: /j.atmosenv
2 674 C. Gorlé et al. / Atmospheric Environment 43 (2009) turbulence profiles of the approach flow affect the results for the downstream velocity, turbulence and concentration (Castro and Robins, 1977; Saathoff et al., 1998; Gao and Chow, 2005). Consequently, horizontal variations in the inlet ABL profiles will influence the downstream prediction of these quantities. For RANS simulations with the standard k e model several solutions have been proposed which result in a horizontally homogeneous velocity profile. The most frequently used method is based on the boundary conditions proposed by Richards and Hoxey (1993). They suggested profiles for velocity, turbulence kinetic energy and turbulence dissipation rate in combination with a logarithmic wall function. When using their inlet profiles without the proposed wall function, one will not obtain a horizontally homogeneous velocity profile, so either the correct wall function has to be implemented in the code (Hargreaves and Wright, 2007), or the constants in the existing rough wall functions have to be adapted to fit with the inlet conditions (Blocken et al., 2007b). A drawback of these methods is that the proposed inlet boundary condition for the turbulence kinetic energy is a constant value over the full height of the lower ABL, which might not correspond to the profiles measured in the atmosphere or generated in a wind tunnel. Although the influence hereof on the velocity field is probably limited, the influence on the concentration field could be more important, since turbulence will enhance dispersion. The goal of this study is to propose a method that allows the correct reproduction of the ABL profiles of velocity and turbulence kinetic energy as measured in wind tunnel tests, using RANS simulations with the standard k e turbulence model. The simulations will be performed using Fluent The influence of the proposed correction on the dispersion of tracer particles in the ABL over flat terrain will be determined. 2. Characteristics of modelled ABL The ABL considered in this study was taken from the CEDVAL database. The profile was generated in the BLASIUS wind tunnel at the Meteorological Institute of the University of Hamburg (Donat, 1995). Using modified Standen Spires (Standen, 1972) and a uniform LEGO-roughness an ABL was modelled at a scale of 1:200. Profiles of mean velocity and turbulence were obtained from two-component LDV measurements in one horizontal and one vertical plane. They are reported in terms of three mean velocity components U (longitudinal), V (lateral), W (vertical) [m s 1 ], three turbulence intensities TI u,ti v,ti w [%] and two Reynolds stresses Cu 0 v 0 D and Cu 0 w 0 D [m 2 s 2 ]. All profiles are shown in Fig. 1. For simulations with the k e turbulence model, the boundary conditions for the turbulence should preferably be described in Fig. 2. k profile calculated for ABL from CEDVAL database. terms of the turbulence kinetic energy k and the turbulence dissipation rate e. The reported turbulence intensity profiles should thus be converted in profiles for k and e. In the present study only the u- component of the measured turbulence, Cu 0 D ¼ TI u U=100, is used for determination of the k-profile: k ¼ 3 2 hu0 i 2 (1) The resulting profile for k is shown in Fig. 2. Using Eq. (1), a rather high value of turbulence kinetic energy is achieved, which will clearly demonstrate the problem related to the horizontal homogeneity of the profile. As the k e model is an isotropic turbulence model, the resulting velocity fluctuations extracted by the model for the determination of the turbulent dispersion will be correct for the u-component, but will be overpredicted for the v- and w-components. The profile for e will be determined according to Eq. (8). 3. ABL modelling using RANS with the standard k e turbulence model The simulations presented in this section were performed on a 2D domain of 4 m length and 1 m height (model scale), with a grid of cells. A gradient was applied in the vertical direction to have the centre point of the wall adjacent cell at a height of m. Richards and Hoxey (1993) outline the following conditions for modelling the ABL in steady incompressible two-dimensional flow using the standard k e model: Fig. 1. ABL profiles from CEDVAL database.
3 C. Gorlé et al. / Atmospheric Environment 43 (2009) The vertical velocity is zero. The pressure is constant. The shear stress, s 0 is constant: ðm l þ m t Þ du dz ¼ s 0 ¼ ru 2 * (2) where m l and m t are the laminar and turbulent viscosities, respectively, r is the density and u * is the boundary layer friction velocity. m t is related to k and e using r and a constant C m as follows: m t ¼ rc m k 2 e The turbulence kinetic energy and dissipation rate satisfy their conservation equations which are given by: v mt vk vu 2 þ m vz s k vz t re ¼ 0 (4) vz v mt ve vu 2 e þ C vz s e vz 1 m t vz k C 2r e2 k ¼ 0 (5) where s k is the turbulence kinetic energy Prandtl number, s e the turbulence dissipation Prandtl number and C 1 and C 2 are constants. If the profiles defined for U, k and e satisfy Eqs. (2 5), a horizontally homogeneous ABL should be obtained. In the following sections different combinations of inlet profiles and their effect on the horizontal homogeneity of the ABL are discussed Constant inlet value for k The approach described by Richards and Hoxey (1993) consists of defining the following inlet boundary conditions for U: UðzÞ ¼ u * z þ z0 ln k z 0 where z 0 is the ABL roughness height and k the von Karman constant. A good log law fit for the velocity profile shown in Fig. 1 is found for a model scale roughness of m and a friction velocity of m s 1 for k ¼ (3) (6) For k and e the inlet conditions should be given by: kðzþ ¼ eðzþ ¼ p u2* ffiffiffiffiffi (7) C m u 3 * kðz þ z 0 Þ These inlet conditions imply equilibrium of turbulence production and dissipation: qffiffiffiffiffi eðzþ ¼ C m kðzþ du dz At the ground boundary, the retarding shear stress s g ¼ ru 2 *g should be determined with the local friction velocity equal to: u * g ¼ (8) (9) ku g ln z g þ z 0 =z0 (10) Combining these boundary conditions with the following relation for the constants of the turbulence model provides a solution to the standard k e model: qffiffiffiffiffi k 2 ¼ðC 2 C 1 Þs e C m (11) When using a commercial CFD code, however, the wall boundary condition is usually not as prescribed above. Blocken et al. (2007b) provide a solution for this by deriving a relationship which brings the rough wall functions in equilibrium with the inlet profiles. For Fluent this relation is given by: k s;abl ¼ 9:793z 0 C s (12) where k s is the roughness height as defined in Fluent for the rough wall function and C s is a constant required for the wall function. k s should be smaller than the height of the centre point of the wall adjacent cell and was consequently set to m. The resulting value for C s is The Fluent Manual states that C s should be smaller than or equal to 1. This requirement, however, results in a first cell height of minimum z 0, which is too large to obtain a good resolution. Therefore, the requirement of limiting C s to 1 is ignored and a higher value is defined through a User Defined Function. The inlet conditions for the CFD simulations for k and e as described by Eqs. (7) and (8) are plotted in Fig. 3. For comparison the value of k derived from the wind tunnel measurements is also Fig. 3. Comparison of inlet profiles from wind tunnel measurements and those from Richards and Hoxey (1993) and the profiles at the domain outlet.
4 676 C. Gorlé et al. / Atmospheric Environment 43 (2009) shown. The plot shows that the turbulence kinetic energy is underpredicted when using Eq. (7) with the standard value for C m of The profiles obtained for U, k and e at the outlet are almost identical to the inlet conditions. Identical results were obtained, with a grid refined by 2 in all directions. This refinement does require to decrease k s and increase C s by a factor of 2. These results show that the approach described by Blocken et al. (2007b) results in homogeneous profiles throughout the domain. The inlet condition for k as described by Eq. (7) does however underestimate the values derived from the wind tunnel tests Measured inlet profile for k In Franke et al. (2007) it is recommended to use measured values for the turbulence quantities, if available, as inlet boundary conditions. Therefore a simulation was performed to verify the homogeneity of the ABL when using the measured values for k, combined with the boundary conditions for the velocity, the turbulence dissipation and the wall as described before. The resulting profiles at the in- and outlet are shown in Fig. 4.It is observed that the new inlet values for k are not maintained through the domain. In addition the higher value of k seems to slightly disturb the horizontal homogeneity of the velocity profile. The development of the inlet profiles through the domain can be explained by the fact that the new profile for k in combination with the given profiles for U and e does not satisfy Eqs. (2) (5). The profiles will then develop to values that do provide a solution, which in this case results in underestimation of the turbulence kinetic energy Correction of profile for e and turbulence model constants The results from Fig. 4 show that, if the inlet profiles should be horizontally homogeneous through the domain, profiles for k and e have to be found that result in an approximate solution for the turbulence model when combined with the logarithmic profile for U. If equilibrium between turbulence dissipation and production is imposed, the profile for e can be determined from Eq. (9). Yang et al. (2007) already derived that substituting this relation for e in the conservation equation for k, results in a solution for k of the following form: qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi kðzþ ¼ Alnðz þ z 0 ÞþB (13) where A and B are constants that can be determined by fitting the equations to the measured profile of k. For the profile under consideration A ¼ 0.11 and B ¼ According to Eq. (9), e will then be given by: pffiffiffiffiffi C mu* qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi eðzþ ¼ Alnðz þ z kðz þ z 0 Þ 0 ÞþB (14) Using these relations, the conservation equation for k is satisfied, but Eqs. (2) and (5) are not. Therefore, simply imposing Eqs. (13) and (14) as inlet boundary condition will not provide a correct solution and will most likely still result in a development of the ABL profiles inside the domain. The first remaining equation to be satisfied is Eq. (2). Substituting Eqs. (13) and (14), the relation to be fulfilled is derived as: C m ¼ u 4 * Alnðz þ z 0 ÞþB (15) As it is not possible to impose a value for C m that varies with height, the value for C m in the wall adjacent cell will be used. For the given profile this results in C m ¼ The second remaining equation is the conservation equation for e. Again substituting Eqs. (13) and (14), it can be derived that the following relation should be satisfied for s e : k 2 A=2 þ kðzþ 2 s e ðzþ ¼ u 2 * ðc 2 C 1 ÞkðzÞ½1 Eðz þ z 0 ÞŠ where E is a function given by: E ¼ 1 e B=A p ffiffiffi p qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 þ Erf kðzþ 2 =A q ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2 kðzþ 2 =A (16) (16a) For the present values of k(z), B and A, an approximate solution is obtained when assuming E ¼ 0, which considerably simplifies the equation. The resulting value for s e also varies with height and can be imposed in Fluent through a User Defined Function. Fig. 5 shows the resulting profiles for k, e and s e with a value of C m equal to The profiles were used for a simulation in combination with the logarithmic profile for U (Eq. (6)). k s and C s were identical to the values used before (Eq. (12)). The profiles obtained for U, k and e at the outlet are shown in Fig. 6. The profile for U is perfectly homogeneous, while the profiles Fig. 4. Comparison of inlet profiles from wind tunnel measurements with the profiles at the domain outlet.
5 C. Gorlé et al. / Atmospheric Environment 43 (2009) for k and e still show a small deviation. This is explained by the fact that the provided inlet profiles are not an exact solution for the constant value of C m. When comparing the results to Fig. 4 it is however obvious that a much better reproduction of the turbulence kinetic energy is obtained. 4. Influence of the turbulence kinetic energy profile on small particle dispersion The influence of the different profiles of turbulence kinetic energy on the dispersion of spherical particles with a 1 mm diameter was investigated by performing simulations for dispersion in an ABL over flat terrain. The same problem has been solved with a simple Gaussian model for comparison with the CFD results Computational domain and boundary conditions The domain represents a rectangular wind tunnel test section with dimensions m as shown in Fig. 7. A circular source with a diameter of 1 cm is located at the ground, 1.2 m behind the domain inlet. Air is entering the domain through this source with a velocity of 1 m s 1 and 20,000 particles are released from random positions within the source at the same velocity. The reference wind speed of the incoming boundary layer is 5.61 m s 1 at 0.4 m height. A symmetric boundary condition is applied at half the width of the domain. The grid resolution is 0.5 mm at the source and increases gradually at locations further away from the source. The total amount of cells is 340,000. The centrepoint of the wall adjacent cell, z 1, is at a height of m, resulting in z þ values just above 30. The value of k s was set to m. The grid dependency of the solution was checked by obtaining a solution with a grid refined by 2 in all directions within the region of interest. The results presented did not show significant changes when performing this refinement. The two sets of boundary conditions used will be referred to as Type 1 (low k) and Type 2 (high k) and are presented in Table 1. The values used for fitting the measured profiles of U and k to the equations are u * ¼ m s 1, z 0 ¼ m, A ¼ 0.11 and B ¼ Particle dispersion modelling The particle dispersion is modelled using a discrete phase model, which determines the particle trajectory by integrating the force balance of the particle written in a Lagrangian reference frame. In the present analysis the only forces considered are the inertia and the drag force, hence the force balance can be written as: Fig. 5. Inlet profiles for k and e and corresponding values for s e. Fig. 6. Comparison of inlet profiles from Eqs. (6, 13, 14) with the profiles at the domain outlet.
6 678 C. Gorlé et al. / Atmospheric Environment 43 (2009) crossing time is defined as t cross ¼ sln½1 ðl e =sju U p jþš where s is the particle relaxation time and L e is the eddy length scale. Concentrations at specified points inside the domain are determined by calculating the time that each particle spends inside a small cell volume surrounding the point of interest. Using the total time spent by all particles in the cell, the non-dimensional concentration C * ¼ C local /C source is given by: C cell ¼ Q P t i;cell NV cell (19) Fig. 7. Sketch of simulation domain. where Q is the source strength (source flow rate), N is the total number of particles tracked through the domain and V cell is the volume of the cell. The resulting values for the concentrations presented in the following sections have been verified to be independent of the choice of N and V cell. du p ¼ 18m dt d 2 p r U Up pc c (17) where d p is the particle diameter, r p the particle density (1550 kg m 3 ) and C c is the Cunningham correction factor, which is equal to 1.17 for particles of 1 mm in air at standard conditions. The turbulent dispersion is accounted for by integrating the trajectory equations for individual particles using the instantaneous fluid velocity U þ u 0 (t) along the particle path. The instantaneous fluid velocity is determined using a stochastic method, i.e. the random walk model (Fluent Inc., 2006). In this method the fluctuating velocity components are discrete piecewise constant functions of time given by: qffiffiffiffiffiffiffiffiffiffi u 0 ¼ z 2k=3 (18) z is a normally distributed random number between 1 and 1, kept constant over an interval of time which indicates the time that a particle is assumed to interact with an eddy. This value is given by the smaller value of the characteristic lifetime of the eddies and the particle eddy crossing time. The characteristic life time is given by s e ¼ T L log(r) where r is a random number between 0 and 1 and T L is the integral time scale given by T L zcðk=eþ. C is a constant which was determined from C ¼ 4=3C 0 with C 0 ¼ C 0N = ð1 þ 7:5 C0N 2 ÞRe 1:64 l. In this equation C 0N is approximately equal to 6 (Oesterlé, 2006). The p Taylor-scale Reynolds number was determined from Re l ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ð20=3þre L with Re L ¼ k 2 =en (Pope, 2008). The resulting value for T L is 0.28 for boundary conditions Type 1 and 0.23 for boundary conditions Type 2. The particle eddy Table 1 Boundary conditions and model coefficients for Type 1 and Type 2 (u * and z 0 from fitting the velocity profile; A and B from fitting the k profile; k, C 1 and C 2 are constants (0.4187, 1.44 and 1.92, respectively); z 1 is the centre point of the first cell). Inlet profiles u k e Type 1 (low k) Type 2 (high k) u * ln z þ z 0 u * ln z þ z 0 k z 0 k z 0 u 2* pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi pffiffiffiffiffi Alnðz þ z 0 ÞþB C m u 3 * kðz þ z 0 Þ Turbulence model C m 0.09 Coefficients s e 1.3 Wall function k s z 1 z 1 pffiffiffiffiffi C mu* kðz þ z 0 Þ qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi Alnðz þ z 0 ÞþB u 4 * Alnðz 1 þ z 0 ÞþB A=2 þ kðzþ 2 ðc 2 C 1 Þ u2 * k 2kðzÞ Coefficients C s 9:793z 0 k s 9:793z 0 k s 4.3. Gaussian dispersion model The Gaussian dispersion model is governed by the following equation when considering the dispersion along the plume centerline (y ¼ 0): " " # " ## Q ðh z Cðx;z;HÞ¼ exp d Þ 2 ðhþz 2ps y s z u 2s 2 þexp d Þ 2 z 2s 2 z (20) where Q is the source strength, u is the velocity and should be a constant value instead of a profile as in the CFD calculations. H is the effective source height and should include the influence of the plume rise due to buoyancy and momentum. The value for H was determined from the results of the CFD simulations as m. The value for u was set to 2 m s 1, which is the free stream velocity at this height. s y and s z are the dispersion coefficients, whose values increase with increasing turbulence. They may be determined using s y ¼ ax 0:093 and s d z ¼ bx c. Three different stability classes are d considered: C (slightly unstable), D (neutral) and E (slightly stable). The values used for the constants a, b and c depend on the stability class and are given by: a C ¼ 0.2, b C ¼ 0.111, c C ¼ 0.911, a D ¼ 0.130, b D ¼ 0.105, c D ¼ 0.827, a E ¼ 0.098, b E ¼ and c E ¼ (Owen Harrop, 2002) Results The contours for flow velocity and turbulence kinetic energy on the symmetry plane obtained from the CFD simulations are shown in Fig. 8 for both types of boundary conditions. The plots show that there is no significant difference in the velocity, while the difference in turbulence kinetic energy is around 40% at the location of interest. The resulting differences in the particle tracks can therefore be attributed to the difference in the fluctuating part of the flow velocity, determined from the turbulence kinetic energy. The influence of this difference in turbulence kinetic energy on the particle tracks is shown in Fig. 9. This plot shows the y- and z- position of the particles when they cross a plane perpendicular to the flow direction at several distances behind the source. A distinctly larger dispersion of the particles is observed when using the boundary conditions and model constants that result in a higher turbulence kinetic energy. In Fig. 10 the resulting non-dimensional concentration values on planes perpendicular to the flow direction at 0.02 and 0.05 m behind the source are plotted. The plots show that a 40% lower turbulence kinetic energy results in a 90% higher value for the maximum value of C *. Fig. 11 shows a plot of the values on the symmetry plane (y ¼ 0.65 m) at 0.05 and 0.1 m behind the source. For comparison
7 C. Gorlé et al. / Atmospheric Environment 43 (2009) Fig. 8. Comparison of velocity and turbulence kinetic energy contours at symmetry plane for boundary conditions Type 1 and Type 2. Fig. 9. Particle positions on three planes downstream of the source. Fig. 10. Non-dimensional concentrations on two planes located downstream of the source, perpendicular to the flow direction.
8 680 C. Gorlé et al. / Atmospheric Environment 43 (2009) Fig. 11. Non-dimensional concentrations on the symmetry plane at two locations downstream of the source. the results from the Gaussian dispersion model are included. The plot shows that the difference between the two CFD simulations is approximately equivalent to the difference between stability classes C (slightly unstable) and E (slightly stable). The CFD simulations show a lower concentration for the high turbulence kinetic energy below m at 0.05 m behind the source and below m at 0.1 m behind the source. Above these heights, the concentration obtained for the higher turbulence is higher due to the increase in the dispersion of the particles. This is the major difference with the results from the Gaussian model, which show a lower concentration for the lower stability classes over the full height. This difference can be explained by considering the formula used for s y and s z. s y depends largely on the stability class, while the variation in s z is limited. Consequently, the dispersion is mainly increasing in the y-direction, leading to lower concentrations over the full height at the centerline, while the distribution over the height is less affected. 5. Discussion The method proposed involves the modification of two constants in the standard k e turbulence model. The value of C m is changed from 0.09 to and the value of s e is changed from a constant value of 1.3 to values varying between 2.2 and 3.0 for the case studied. The significant decrease of C m was already proposed by Richards and Hoxey (1993). They obtained a value of for the modelling of the atmospheric boundary layer at Silsoe. According to Eq. (11), they then increased the turbulent dissipation Prandtl number s e to 3.2. They suggested that a possible cause for the deviation of C m is the significant low frequency contribution to the velocity fluctuations, which contribute little or nothing to the Reynolds stress. The Reynolds stresses are related to the turbulent viscosity by ru 0 w 0 ¼ m t ðdu=dzþ. Substituting Eq. (3) in this relation one obtains u 0 w 0 ¼ C m ðk 2 =eþðdu=dzþ. Consequently, if u 0 w 0 is small compared to k, the value of C m should be decreased. The increase in the value of s e implies that the diffusion of e is decreased. By imposing higher values for s e at the wall the diffusion away from the wall of the value of e defined near the wall is limited. The observation that C m should be decreased has some implications for the modelling of dispersion with the convection diffusion equation. When using this method, the dispersion will be governed by the parameter m t /Sc t, where Sc t is the turbulent Schmidt number. By decreasing C m and imposing higher values for k at the inlet boundary, there will be no net increase in m t. Consequently, when using the proposed boundary conditions and modifications of the constants in the turbulence model, no enhanced dispersion will be observed when using the convection diffusion equation. In several studies (e.g. Di Sabatino et al., 2007; Tominaga and Stathopoulos, 2007; Riddle et al., 2004) it is suggested that for the modelling of dispersion in urban areas the turbulent Schmidt number should be lower than the standard value of 0.7. If one considers that the dispersion will be governed by all spatial and temporal scales of the velocity fluctuations, this observation is consistent with the explanation given by Richards and Hoxey (1993) for the lower value of C m : if the velocity fluctuations have significant low frequency contributions, which do not contribute to the Reynolds stresses and therefore not to the value of m t, the value of Sc t should be reduced to obtain a higher value for the dispersion. The effect of changing the constants on the simulation of flows around building configurations should be further investigated. When performing such validation studies it is suggested to determine the boundary conditions using a value for k based on the three fluctuating velocity components: k ¼ 1/2(u 02 þ v 02 þ w 02 ). This should result in a more physically correct value for k. The resulting velocity fluctuations will then be underpredicted in the u-direction and slightly overpredicted in the v- and w-directions. In an attempt to compensate for the isotropy of the turbulence model one could consider the introduction of constant factors in the calculation of the different components of the velocity fluctuations. These could be determined according to the ratio between the magnitude of the three components of the velocity fluctuations that would be expected in the ABL under consideration. The constant C used in the determination of the Lagrangian time scale T L has an important effect on the resulting values for the concentration. The values for C that were determined with the formulas from Oesterlé (2006) and Pope (2008) are between the range of and 0.6 that is given by Oesterlé (2006). They can be considered an educated guess, but are by no means exact values. The default value used in Fluent is When using this value, there is an increase in the maximum concentration of 12 18% compared to, respectively, boundary conditions Type 2 with T L ¼ 0.23 and boundary conditions Type 1 with T L ¼ This demonstrates the importance of verifying the sensitivity of the results to the value of C, since the constant cannot be determined exactly. 6. Conclusions The study presented aimed at proposing a method to correctly model the turbulence kinetic energy in the atmospheric boundary layer and to determine the influence of the turbulence kinetic energy on the dispersion of small particles. The proposed method is based on the boundary conditions that were defined by Richards and Hoxey (1993) combined with the correction of the rough wall function implemented in the CFD code as described by Blocken et al. (2007b). The drawback of this method is the inlet boundary condition for the turbulence kinetic energy, which has to be constant over the height of the ABL.
9 C. Gorlé et al. / Atmospheric Environment 43 (2009) Therefore, the inlet boundary conditions for the turbulence kinetic energy and the turbulence dissipation rate were allowed to vary with height according to the relations derived by Yang et al. (2007). In addition, it was derived that the constants C m and s e, used by the turbulence model, should be determined according to Eq. (15) (using the value in the wall adjacent cell) and Eqs. (16 and 16a). When using these relations, the proposed boundary conditions approximately satisfy the conditions that should be fulfilled to model a horizontally homogeneous boundary layer, both for the velocity and the turbulence kinetic energy. It was shown that the influence of the turbulence kinetic energy on the dispersion of small particles is significant by considering the dispersion of 20,000 particles from a circular source over flat terrain. A comparison was made between a simulation using the correction of the rough wall function to obtain a horizontally homogeneous velocity profile, combined with the constant value for the turbulence kinetic energy and the standard values for the constants in the turbulence model and a simulation that also included the correction of the rough wall function, but combined with a profile for the turbulence kinetic energy varying with height and the modified constants for the turbulence model. This comparison showed that a 40% lower turbulence kinetic energy results in a 90% higher prediction of the maximum concentration close to the source. This difference is comparable to the difference that is obtained between the results obtained from a Gaussian dispersion model for stability classes C and E. Acknowledgements This work was supported by IWT Vlaanderen, the Institute for the Promotion of Innovation by Science and Technology in Flanders, through the SBO project NanoSoc. References Blocken, B., Carmeliet, J., Stathopoulos, T., 2007a. Cfd evaluation of wind speed conditions in passages between parallel buildings effect of wall-function roughness modifications for the atmospheric boundary layer flow. Journal of Wind Engineering and Industrial Aerodynamics 95, Blocken, B., Stathopoulos, T., Carmeliet, J., 2007b. Cfd simulation of the atmospheric boundary layer: wall function problems. Atmospheric Environment 41, Castro, I., Robins, A., The flow around a surface-mounted cube in uniform and turbulent streams. Journal of Fluid Mechanics 79 (2), Chang, C., Meroney, R., Concentration and flow distributions in urban street canyons: wind tunnel and computational data. Journal of Wind Engineering and Industrial Aerodynamics 91, Dixon, N., Boddy, J., Smalley, R., Tomlin, A.S., Evaluation of a turbulent flow and dispersion model in a typical street canyon in York, UK. Atmospheric Environment 40, Di Sabatino, S., Buccolieri, R., Pulvirenti, B., Britter, R., Simulations of pollutant dispersion within idealised urban-type geometries with cfd and integral models. Atmospheric Environment 41 (37), Donat, J., Windkanalexperimente zur ausbreitung von schwergasstrahlen. Ph.D. thesis, University of Hamburg, Meteorological Institute. FLuent 6.3 Users s Guide, Fluent Inc. Franke, J., Hellsten, A., Schlunzen, H., Carissimo, B. (Eds.), Best Practice Guideline for the CFD Simulation of Flows in the Urban Environment COST Action 732. Franke, J., Hirsch, C., Jensen, A., Krus, H., Schatzmann, M., Westbury, P., Miles, S., Wisse, J., Wright, N., Recommendations on the use of cfd in wind engineering. In: COST Action C14, Impact of Wind and Storm on City Life Built Environment, Proceedings of the International Conference on Urban Wind Engineering and Building Aerodynamics. von Karman Institute. Gao, Y., Chow, W., Numerical studies on air flow around a cube. Journal of Wind Engineering and Industrial Aerodynamics 93, Garcia Sagrado, A., van Beeck, J., Rambaud, P., Olivari, D., Numerical and experimental modelling of pollutant dispersion in a street canyon. Journal of Wind Engineering and Industrial Aerodynamics 90, Hargreaves, D., Wright, N., On the use of the k e model in commercial cfd software to model the neutral atmospheric boundary layer. Journal of Wind Engineering and Industrial Aerodynamics 95, Leitl, B., Kastner-Klein, P., Rau, M., Meroney, R.N., Concentration and flow distributions in the vincinity of u-shaped buildings: wind-tunnel and computational data. Journal of Wind Engineering and Industrial Aerodynamics 67 68, Meroney, R., Leitl, B., Rafailidis, S., Schatzmann, M., Wind-tunnel and numerical modelling of flow and dispersion about several building shapes. Journal of Wind Engineering and Industrial Aerodynamics 81, Oesterlé, B. (Ed.), Ecoulements multiphasiques des fondements aux methodes d ingénierie. Owen Harrop, D., Air Quality Assessment and Management A Practical Guide. Clay s Library of Health and the Environment. Pope, S. (Ed.), Turbulent Flows. Richards, P., Hoxey, R., Appropriate boundary conditions for computational wind engineering models using the k e turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 46-47, Riddle, A., Carruthers, D., Sharpe, A., McHugh, C., Stocker, J., Comparisons between fluent and adms for atmospheric dispersion modelling. Atmospheric Environment 38, Saathoff, P., Stathopoulos, T., Wu, H., The influence of freestream turbulence on nearfield dilution of exhaust from building vents. Journal of Wind Engineering and Industrial Aerodynamics 77 78, Schatzmann, M., Leitl, B., Validation and application of obstacle resolving dispersion models. Atmospheric Environment 36, Standen, N., A Spire for Generating Thick Turbulent Shear Layers for Neutral Wind Simulation in Wind Tunnels. Technical Report LTR-LA-94. National Aeronautical Establishment, Ottawa, Canada. Tominaga, Y., Stathopoulos, T., Turbulent schmidt numbers for cfd analysis with various types of flowfield. Atmospheric Environment 41, Wang, X., McNamara, K., Effects of street orientation on dispersion at or near urban street intersections. Journal of Wind Engineering and Industrial Aerodynamics 95, Yang, W., Jin, X., Jin, H., Gu, M., Chen, S., Application of new inflow boundary conditions for modeling equilibrium atmosphere boundary layer in rans-based turbulence models. In: Proceedings of International Conference on Wind Engineering 2007.
PRELIMINARY STUDY OF COMPUTATIONAL SETUP FOR URBAN STREET CANYONS. by MUHAMMAD NOOR AFIQ WITRI, M.Eng
PRELIMINARY STUDY OF COMPUTATIONAL SETUP FOR URBAN STREET CANYONS by MUHAMMAD NOOR AFIQ WITRI, M.Eng 1 CONTENTS 1.Introduction 2.Building Configuration 3.Boundary Condition 4.Previous Works 5.Summary 2
More informationNumerical evaluation of pollutant dispersion in the built environment: comparisons between models and experiments
PREPRINT: Blocken B, Stathopoulos T, Saathoff P, Wang X. 2008. Numerical evaluation of pollutant dispersion in the built environment: comparisons between models and experiments. Journal of Wind Engineering
More informationNumerical simulation of dispersion around an isolated cubic building: Model evaluation of RANS and LES. Yoshihide Tominaga a and Ted Stathopoulos b
Accepted on 3 April for publication in the Building and Environment Numerical simulation of dispersion around an isolated cubic building: Model evaluation of RANS and Yoshihide Tominaga a and Ted Stathopoulos
More informationMODELING THE NEUTRAL ATMOSPHERIC BOUNDARY LAYER BASED ON THE STANDARD k-ε TURBULENT MODEL: MODIFIED WALL FUNCTION
The Seventh Asia-Pacific Conference on Wind Engineering, November 8-12, 29, Taipei, Taiwan MODELING THE NEUTRAL ATMOSPHERIC BOUNDARY LAYER BASED ON THE STANDARD k-ε TURBULENT MODEL: MODIFIED WALL FUNCTION
More informationZhang Yue 1, Guo Wei 2, Wang Xin 3, Li Jiawu 4 1 School of High way, Chang an University, Xi an, Shanxi, China,
The Eighth Asia-Pacific Conference on Wind Engineering, December 10 14, 2013, Chennai, India Analysis of Two Kinds of Boundary Conditions for Simulating Horizontally Homogenous Atmosphere Boundary Layer
More informationPart I: Overview of modeling concepts and techniques Part II: Modeling neutrally stratified boundary layer flows
Physical modeling of atmospheric boundary layer flows Part I: Overview of modeling concepts and techniques Part II: Modeling neutrally stratified boundary layer flows Outline Evgeni Fedorovich School of
More informationANALYZING WIND-DRIVEN RAIN ON A BUILDING FACADE USING THE LASER PRECIPITATION MONITOR (LPM)
ANALYZING WIND-DRIVEN RAIN ON A BUILDING FACADE USING THE LASER PRECIPITATION MONITOR (LPM) Ayyapan Kumaraperumal, Chris H. Sanders, Paul H. Baker, Graham H. Galbraith, and Don McGlinchey ABSTRACT Centre
More informationOn the validation study devoted to stratified atmospheric flow over an isolated hill
On the validation study devoted to stratified atmospheric flow over an isolated hill Sládek I. 2/, Kozel K. 1/, Jaňour Z. 2/ 1/ U1211, Faculty of Mechanical Engineering, Czech Technical University in Prague.
More information17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary
17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary INVESTIGATION OF VENTILATION AND AIR QUALITY IN URBAN SQUARES
More informationEnvironmental Pollution
Gousseau P, Blocken B, van Heijst GJF. 2012. Large-Eddy Simulation of pollutant dispersion around a cubical building: Analysis of the turbulent mass transport mechanism by unsteady concentration and velocity
More informationUncertainty quantification for RANS simulation of flow over a wavy wall
Uncertainty quantification for RANS simulation of flow over a wavy wall Catherine Gorlé 1,2,3, Riccardo Rossi 1,4, and Gianluca Iaccarino 1 1 Center for Turbulence Research, Stanford University, Stanford,
More informationADAPTATION OF THE REYNOLDS STRESS TURBULENCE MODEL FOR ATMOSPHERIC SIMULATIONS
ADAPTATION OF THE REYNOLDS STRESS TURBULENCE MODEL FOR ATMOSPHERIC SIMULATIONS Radi Sadek 1, Lionel Soulhac 1, Fabien Brocheton 2 and Emmanuel Buisson 2 1 Laboratoire de Mécanique des Fluides et d Acoustique,
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationarxiv: v1 [physics.flu-dyn] 11 Oct 2012
Low-Order Modelling of Blade-Induced Turbulence for RANS Actuator Disk Computations of Wind and Tidal Turbines Takafumi Nishino and Richard H. J. Willden ariv:20.373v [physics.flu-dyn] Oct 202 Abstract
More informationCFD analysis of convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer
Accepted for publication in International Journal of Heat and Mass Transfer CFD analysis of convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer Thijs Defraeye a, *,
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationAIJ COOPERATIVE PROJECT FOR PRACTICAL APPLICATIONS OF CFD TO URBAN VENTILATION
The Seventh Asia-Pacific Conference on Wind Engineering, November 8-2, 29, Taipei, Taiwan AIJ COOPERATIVE PROJECT FOR PRACTICAL APPLICATIONS OF CFD TO URBAN VENTILATION Ryuichiro Yoshie, Akashi Mochida
More information17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary
17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary SIMULATIONS OF TRAFFIC RELATED POLLUTANTS IN A MAIN STREET
More informationCFD Study of Flow Over Parallel Ridges with Varying Height and Spacing
Proceedings of the World Congress on Engineering 21 Vol II WCE 21, June 3 - July 2, 21, London, U.K. CFD Study of Flow Over Parallel Ridges with Varying Height and Spacing Lee Chin Yik, Salim Mohamed Salim,
More informationINTER-COMPARISON AND VALIDATION OF RANS AND LES COMPUTATIONAL APPROACHES FOR ATMOSPHERIC DISPERSION AROUND A CUBIC OBSTACLE. Resources, Kozani, Greece
INTER-COMPARISON AND VALIDATION OF AND LES COMPUTATIONAL APPROACHES FOR ATMOSPHERIC DISPERSION AROUND A CUBIC OBSTACLE S. Andronopoulos 1, D.G.E. Grigoriadis 1, I. Mavroidis 2, R.F. Griffiths 3 and J.G.
More informationA combined application of the integral wall model and the rough wall rescaling-recycling method
AIAA 25-299 A combined application of the integral wall model and the rough wall rescaling-recycling method X.I.A. Yang J. Sadique R. Mittal C. Meneveau Johns Hopkins University, Baltimore, MD, 228, USA
More informationBOUNDARY LAYER STRUCTURE SPECIFICATION
August 2017 P09/01X/17 BOUNDARY LAYER STRUCTURE SPECIFICATION CERC In this document ADMS refers to ADMS 5.2, ADMS-Roads 4.1, ADMS-Urban 4.1 and ADMS-Airport 4.1. Where information refers to a subset of
More informationSTREET CANYON MODEL MISKAM AND
STREET CANYON MODEL MISKAM AND CONCENTRATION ESTIMATION COUPLING THE MICROMIXING LAGRANGIAN MODEL Giovanni Leuzzi 1, Márton Balczó 2, Andrea Amicarelli 3, Paolo Monti 1, Joachim Eichhorn 3 and David J.
More informationInvestigation of Flow Profile in Open Channels using CFD
Investigation of Flow Profile in Open Channels using CFD B. K. Gandhi 1, H.K. Verma 2 and Boby Abraham 3 Abstract Accuracy of the efficiency measurement of a hydro-electric generating unit depends on the
More informationFluid Mechanics. Chapter 9 Surface Resistance. Dr. Amer Khalil Ababneh
Fluid Mechanics Chapter 9 Surface Resistance Dr. Amer Khalil Ababneh Wind tunnel used for testing flow over models. Introduction Resistances exerted by surfaces are a result of viscous stresses which create
More informationNumerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders
Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute
More informationPerformance characteristics of turbo blower in a refuse collecting system according to operation conditions
Journal of Mechanical Science and Technology 22 (2008) 1896~1901 Journal of Mechanical Science and Technology www.springerlink.com/content/1738-494x DOI 10.1007/s12206-008-0729-6 Performance characteristics
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationCFD calculations of the test 2-4 experiments. Author: G. de With
CFD calculations of the test 2-4 experiments Author: G. de With 34. Model setup and boundary conditions Dimensions CFD model: x=1000m / y=100m / z=2000m. CFD Model: Transient simulation, with steady-state
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationInvestigation of Opening Positions on the Natural Ventilation in a Low-Rise Building by CFD Analysis
Proceedings of the 3 rd International Conference on Fluid Flow, Heat and Mass Transfer (FFHMT 16) Ottawa, Canada May 2 3, 2016 Paper No. 151 Investigation of Opening Positions on the Natural Ventilation
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationSimulation of Flow around a Surface-mounted Square-section Cylinder of Aspect Ratio Four
Simulation of Flow around a Surface-mounted Square-section Cylinder of Aspect Ratio Four You Qin Wang 1, Peter L. Jackson 2 and Jueyi Sui 2 1 High Performance Computing Laboratory, College of Science and
More informationINFLUENCE OF THE GEOMETRICAL PARAMETERS OF URBAN CANYONS ON THE CONVECTIVE HEAT TRANSFER COEFFICIENT
INFLUENCE OF THE GEOMETRICAL PARAMETERS OF URBAN CANYONS ON THE CONVECTIVE HEAT TRANSFER COEFFICIENT Andrea VALLATI a, Giorgio GALLI a, Chiara COLUCCi a, Pawel OCLON b,* a Department of DIAEE, Sapienza
More informationTurbulence Laboratory
Objective: CE 319F Elementary Mechanics of Fluids Department of Civil, Architectural and Environmental Engineering The University of Texas at Austin Turbulence Laboratory The objective of this laboratory
More informationExperimental Verification of CFD Modeling of Turbulent Flow over Circular Cavities using FLUENT
Experimental Verification of CFD Modeling of Turbulent Flow over Circular Cavities using FLUENT T Hering, J Dybenko, E Savory Mech. & Material Engineering Dept., University of Western Ontario, London,
More informationIntroduction to Turbulence AEEM Why study turbulent flows?
Introduction to Turbulence AEEM 7063-003 Dr. Peter J. Disimile UC-FEST Department of Aerospace Engineering Peter.disimile@uc.edu Intro to Turbulence: C1A Why 1 Most flows encountered in engineering and
More informationCFD modeling of dust dispersion through Najaf historic city centre
INTERNATIONAL JOURNAL OF ENERGY AND ENVIRONMENT Volume 5, Issue 6, 2014 pp.723-728 Journal homepage: www.ijee.ieefoundation.org TECHNICAL PAPER CFD modeling of dust dispersion through Najaf historic city
More informationSymmetry of Turbulent Characteristics Inside Urban Intersection
Colloquium FLUID DYNAMICS 2007 Institute of Thermomechanics AS CR, v. v. i., Prague, October 24-26, 2007 p.1 Symmetry of Turbulent Characteristics Inside Urban Intersection Radka Kellnerová 1,2 Zbyněk
More informationA Comparative Analysis of Turbulent Pipe Flow Using k And k Models
A Comparative Analysis of Turbulent Pipe Flow Using k And k Models 1 B. K. Menge, 2 M. Kinyanjui, 3 J. K. Sigey 1 Department of Mathematics and Physics. Technical University of Mombasa, P.O BOX 90420-80100,Mombasa,
More informationKeywords - Gas Turbine, Exhaust Diffuser, Annular Diffuser, CFD, Numerical Simulations.
Numerical Investigations of PGT10 Gas Turbine Exhaust Diffuser Using Hexahedral Dominant Grid Vaddin Chetan, D V Satish, Dr. Prakash S Kulkarni Department of Mechanical Engineering, VVCE, Mysore, Department
More informationStorm Water Best Management Practice: Development of Debris Filtering Structure for Supercritical Flow
Storm Water Best Management Practice: Development of Debris Filtering Structure for Supercritical Flow Jungseok Ho 1, Todd Marti 2, and Julie Coonrod 3 1 Department of Civil Engineering, University of
More informationGENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION
25 th National Symposium on Wind Engineering, Tokyo, Japan, 3-5 December 2018 第 25 回風工学シンポジウム (2018) GENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION Takafumi
More informationAnalysis of Heat Transfer in Pipe with Twisted Tape Inserts
Proceedings of the 2 nd International Conference on Fluid Flow, Heat and Mass Transfer Ottawa, Ontario, Canada, April 30 May 1, 2015 Paper No. 143 Analysis of Heat Transfer in Pipe with Twisted Tape Inserts
More informationCalculations on a heated cylinder case
Calculations on a heated cylinder case J. C. Uribe and D. Laurence 1 Introduction In order to evaluate the wall functions in version 1.3 of Code Saturne, a heated cylinder case has been chosen. The case
More informationA Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM
206 9th International Conference on Developments in esystems Engineering A Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM Hayder Al-Jelawy, Stefan Kaczmarczyk
More informationDivergence free synthetic eddy method for embedded LES inflow boundary condition
R. Poletto*, A. Revell, T. Craft, N. Jarrin for embedded LES inflow boundary condition University TSFP Ottawa 28-31/07/2011 *email: ruggero.poletto@postgrad.manchester.ac.uk 1 / 19 SLIDES OVERVIEW 1 Introduction
More informationCOMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE
COMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE L. Velázquez-Araque 1 and J. Nožička 2 1 Division of Thermal fluids, Department of Mechanical Engineering, National University
More informationBOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW
Proceedings of,, BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW Yunho Jang Department of Mechanical and Industrial Engineering University of Massachusetts Amherst, MA 01002 Email:
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationABSTRACT INTRODUCTION
Numerical simulation of pulsating flow around a cube C. Dargent, D. Dartus, J. George Institut de Mecanique des Fluides de Toulouse, Avenue du Professeur Camille Soula, 31400 Toulouse, France ABSTRACT
More informationDepartment of Mechanical Engineering
Department of Mechanical Engineering AMEE401 / AUTO400 Aerodynamics Instructor: Marios M. Fyrillas Email: eng.fm@fit.ac.cy HOMEWORK ASSIGNMENT #2 QUESTION 1 Clearly there are two mechanisms responsible
More informationModeling of dispersed phase by Lagrangian approach in Fluent
Lappeenranta University of Technology From the SelectedWorks of Kari Myöhänen 2008 Modeling of dispersed phase by Lagrangian approach in Fluent Kari Myöhänen Available at: https://works.bepress.com/kari_myohanen/5/
More informationSimulation of Aeroelastic System with Aerodynamic Nonlinearity
Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More informationValidation 3. Laminar Flow Around a Circular Cylinder
Validation 3. Laminar Flow Around a Circular Cylinder 3.1 Introduction Steady and unsteady laminar flow behind a circular cylinder, representing flow around bluff bodies, has been subjected to numerous
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
More informationChristophe DUCHENNE 1, Patrick ARMAND 1, Maxime NIBART 2, Virginie HERGAULT 3. Harmo 17 Budapest (Hungary) 9-12 May 2016
Validation of a LPDM against the CUTE experiments of the COST ES1006 Action Comparison of the results obtained with the diagnostic and RANS versions of the flow model Christophe DUCHENNE 1, Patrick ARMAND
More informationLaminar Flow. Chapter ZERO PRESSURE GRADIENT
Chapter 2 Laminar Flow 2.1 ZERO PRESSRE GRADIENT Problem 2.1.1 Consider a uniform flow of velocity over a flat plate of length L of a fluid of kinematic viscosity ν. Assume that the fluid is incompressible
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationMATHEMATICAL MODELING AND NUMERICAL SOLUTION OF 3D ATMOSPHERIC BOUNDARY LAYER
, Vol, Pt, Special Issue Proceedings of International Conference RDAMM 585 MATHEMATICAL MODELING AND NUMERICAL SOLUTION OF D ATMOSPHERIC BOUNDARY LAYER L. Beneš, K. Kozel Department of Technical Mathematics,
More informationWind environmental conditions in passages between two long narrow perpendicular buildings
PRE-PRINT of the article Blocken B, Stathopoulos T, Carmeliet J. 2008. Wind environmental conditions in passages between two long narrow perpendicular buildings. Journal of Aerospace Engineering ASCE 21(4):
More informationNumerical simulation of high pressure gas quenching of H13 steel
journal of materials processing technology 202 (2008) 188 194 journal homepage: www.elsevier.com/locate/jmatprotec Numerical simulation of high pressure gas quenching of H13 steel Jing Wang a, Jianfeng
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 4 HEAT TRANSFER IN CHANNEL FLOW BASIC CONCEPTS BASIC CONCEPTS Laminar
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationDEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M.
DEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M. Sahu 1, Kishanjit Kumar Khatua and Kanhu Charan Patra 3, T. Naik 4 1, &3 Department of Civil Engineering, National Institute of technology,
More informationTurbulence Instability
Turbulence Instability 1) All flows become unstable above a certain Reynolds number. 2) At low Reynolds numbers flows are laminar. 3) For high Reynolds numbers flows are turbulent. 4) The transition occurs
More informationLarge-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall
Large-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall Y. Kim, I.P. Castro, and Z.T. Xie Introduction Wind turbines operate in the atmospheric boundary layer and their
More informationLES evaluation of wind pressures on a standard tall building with and without a neighboring building
LES evaluation of wind pressures on a standard tall building with and without a neighboring building Agerneh K. Dagnew a, Girma T. Bitsuamlak b a PhD candidate, Civil and Environmental Eng. (CEE) Dept./
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationNUMERICAL SIMULATION OF THREE DIMENSIONAL GAS-PARTICLE FLOW IN A SPIRAL CYCLONE
Applied Mathematics and Mechanics (English Edition), 2006, 27(2):247 253 c Editorial Committee of Appl. Math. Mech., ISSN 0253-4827 NUMERICAL SIMULATION OF THREE DIMENSIONAL GAS-PARTICLE FLOW IN A SPIRAL
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
More informationNUMERICAL AND EXPERIMENTAL INVESTIGATIONS OF AIR FLOW AND TEMPERATURE PATTERNS OF A LOW VELOCITY DIFFUSER
NUMERICAL AND EXPERIMENTAL INVESTIGATIONS OF AIR FLOW AND TEMPERATURE PATTERNS OF A LOW VELOCITY DIFFUSER M Cehlin and B Moshfegh Division of Energy and Mechanical Engineering, Department of Technology,
More informationComparison of Numerical Prediction of Pressure Coefficient on Rectangular Small Building
May 26, Volume 3, Issue 5 Comparison of Numerical Prediction of Pressure Coefficient on Rectangular Small Building Neel M. Patel, 2 Satyen D. Ramani PG Student, Department of Civil Engineering, 2 Asst.
More informationPARTICLE DISPERSION IN ENCLOSED SPACES USING A LAGRANGIAN MODEL
IV Journeys in Multiphase Flows (JEM 217) March 27-31, 217, São Paulo, SP, Brazil Copyright 217 by ABCM Paper ID: JEM-217-4 PARTICLE DISPERSION IN ENCLOSED SPACES USING A LAGRANGIAN MODEL Ana María Mosquera
More informationTable of Contents. Foreword... xiii. Preface... xv
Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...
More informationModelling diabatic atmospheric boundary layer using a RANS-CFD code with a k-ε turbulence closure F. VENDEL
Modelling diabatic atospheric boundary layer using a RANS-CFD code with a k-ε turbulence closure F. VENDEL Florian Vendel 1, Guillevic Laaison 1, Lionel Soulhac 1, Ludovic Donnat 2, Olivier Duclaux 2,
More informationComparison of turbulence models for the computational fluid dynamics simulation of wind turbine wakes in the atmospheric boundary layer
Comparison of turbulence models for the computational fluid dynamics simulation of wind turbine wakes in the atmospheric boundary layer D. Cabezón E. Migoya and A. Crespo ABSTRACT An elliptic computational
More informationPRESSURE COEFFICIENTS ACTING UPON THE CYLINDER OBTAINED BY NUMERICAL AND EXPERIMENTAL ANALYSIS
Vol. 13, Issue 2/2017, 149-155, DOI: 10.1515/cee-2017-0020 PRESSURE COEFFICIENTS ACTING UPON THE CYLINDER OBTAINED BY NUMERICAL AND EXPERIMENTAL ANALYSIS Soňa MEDVECKÁ 1,*, Oľga IVÁNKOVÁ 1, Marek MACÁK
More informationBoundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity
Boundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity Tobias Knopp D 23. November 28 Reynolds averaged Navier-Stokes equations Consider the RANS equations with
More informationCHAM Case Study CFD Modelling of Gas Dispersion from a Ruptured Supercritical CO 2 Pipeline
CHAM Limited Pioneering CFD Software for Education & Industry CHAM Case Study CFD Modelling of Gas Dispersion from a Ruptured Supercritical CO 2 Pipeline 1. INTRODUCTION This demonstration calculation
More informationNumerical modeling of the flow conditions in a closed-circuit low-speed wind tunnel
PREPRINT: P. Moonen, B. Blocken, S. Roels and J. Carmeliet, Numerical modeling of the flow conditions in a closed-circuit lowspeed wind tunnel, International Journal of Wind Engineering and Industrial
More informationelements remain in high frequency region and sometimes very large spike-shaped peaks appear. So we corrected the PIV time histories by peak cutting an
The Seventh International Colloquium on Bluff Body Aerodynamics and Applications (BBAA7) Shanghai, China; September 2-6, 2012 LES of fluctuating wind pressure on a 3D square cylinder for PIV-based inflow
More informationWQMAP (Water Quality Mapping and Analysis Program) is a proprietary. modeling system developed by Applied Science Associates, Inc.
Appendix A. ASA s WQMAP WQMAP (Water Quality Mapping and Analysis Program) is a proprietary modeling system developed by Applied Science Associates, Inc. and the University of Rhode Island for water quality
More informationCFD ANALYSIS OF NATURAL VENTILATION IN LARGE SEMI-ENCLOSED BUILDINGS CASE STUDY: AMSTERDAM ARENA FOOTBALL STADIUM Leuven, Belgium
Eleventh International IBPSA Conference Glasgow, Scotland July 27-30, 2009 CFD ANALYSIS OF NATURAL VENTILATION IN LARGE SEMI-ENCLOSED BUILDINGS CASE STUDY: AMSTERDAM ARENA FOOTBALL STADIUM Twan van Hooff
More informationCalculation of the separation streamlines of barchans and transverse dunes
Calculation of the separation streamlines of barchans and transverse dunes H. J. Herrmann a,b, J. S. Andrade Jr. b, V. Schatz a, G. Sauermann a and E. J. R. Parteli a a ICP, University of Stuttgart, Pfaffenwaldring
More informationINTERNATIONAL JOURNAL OF SCIENTIFIC & TECHNOLOGY RESEARCH VOLUME 5, ISSUE 09, SEPTEMBER 2016 ISSN
Numerical Analysis Of Heat Transfer And Fluid Flow Characteristics In Different V-Shaped Roughness Elements On The Absorber Plate Of Solar Air Heater Duct Jitesh Rana, Anshuman Silori, Rohan Ramola Abstract:
More informationNumerical simulations of heat transfer in plane channel flow
Numerical simulations of heat transfer in plane channel flow Najla EL GHARBI 1, 3, a, Rafik ABSI 2, b and Ahmed BENZAOUI 3, c 1 Renewable Energy Development Center, BP 62 Bouzareah 163 Algiers, Algeria
More informationExplicit algebraic Reynolds stress models for internal flows
5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence
More information350 Int. J. Environment and Pollution Vol. 5, Nos. 3 6, 1995
350 Int. J. Environment and Pollution Vol. 5, Nos. 3 6, 1995 A puff-particle dispersion model P. de Haan and M. W. Rotach Swiss Federal Institute of Technology, GGIETH, Winterthurerstrasse 190, 8057 Zürich,
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 19 Turbulent Flows Fausto Arpino f.arpino@unicas.it Introduction All the flows encountered in the engineering practice become unstable
More informationKeywords: Large-eddy simulation, Turbulent coherent structure, Four quadrant analysis, Integral scale
The Eighth Asia-Pacific Conference on Wind Engineering, December 4, 3, Chennai, India NUMERICAL ANALYSIS OF THE MOMENTUM TRANSPORT AND TEMPORAL AND SPATIAL SCALES OF TURBULENT COHERENT STRUCTURES IN THE
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationWind interference effect on an octagonal plan shaped tall building due to square plan shaped tall buildings
Int J Adv Struct Eng (2016) 8:73 86 DOI 10.1007/s40091-016-0115-z ORIGINAL RESEARCH Wind interference effect on an octagonal plan shaped tall building due to square plan shaped tall buildings Rony Kar
More informationRECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD
RECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD N. Jarrin 1, A. Revell 1, R. Prosser 1 and D. Laurence 1,2 1 School of MACE, the University of Manchester,
More informationESTIMATION OF THE FLOW CHARACTERISTICS BETWEEN THE TRAIN UNDERBODY AND THE BALLAST TRACK
BBAA VI International Colloquium on: Bluff Bodies Aerodynamics & Applications Milano, Italy, July, -4 8 ESTIMATION OF THE FLOW CHARACTERISTICS BETWEEN THE TRAIN UNDERBODY AND THE BALLAST TRACK Javier García*,
More informationUniversity of Huddersfield Repository
University of Huddersfield Repository Malviya, Vihar, Gundala, Naresh and Mishra, Rakesh Effect of cross wind on aerodynamic coefficients of ground vehicles. Original Citation Malviya, Vihar, Gundala,
More information