CHAPTER 6 APPLICATION OF NUMERICAL SIMULATION TO VENTILATED FLOW PROBLEMS
|
|
- Magdalen Black
- 5 years ago
- Views:
Transcription
1 101
2 in the direction of the experiment data from Restivo (1979). The same situation was also found by Davidson and Nielsen (1998) for Large Eddy Simulation. However, additional simulations could be performed in order to rule out problems with grid density, inadequate subgrid scale models and improved boundary conditions at the inlet. But there could also be inaccuracy in the experimental data, so further experiments would be needed to investigate the discrepancies between the Large Eddy Simulation and the experimental data from Restivo (1979). However, the simulations concur well with the findings made by Armaly et al. (1983), for all the four different Reynolds numbers both within the transitional and the fully turbulent regions. This further indicates that the data by Restivo could be subjected to uncertainty. Further results about the backward facing geometry and the modelling of transitional flows is given in Bennetsen (C) (2001). Next, a test from the Air Physic Lab at Research Centre Bygholm, The Danish Institute of Agricultural Sciences is studied. The geometry of the room is changed to a different room width, a different slot height and a different outlet location. This has an effect upon the flow geometry which changes from two-dimensional to fully three-dimensional. 6.4 The SJVF-Room Figure 6.89: Geometry of the SJVF-room; heigth = 3 m; length = 8.5 m; width = 5 m; height of slot inlet is m, length of slot outlet is 0.11m and width of 4.8 m. The last room geometry studied within this thesis is based on an experimental room in the Air Physics Laboratory, at Research Centre Bygholm. The length of the room is 8.5 m, the heigh is 3 m and the width is 5.0 m. The inlet which consists of a slot opening in the full width of the room and is placed just below the ceiling. The height of the inlet slot was 0.1 m when fully open. The actual opening height could be adjusted within the range from 0 to 0.1 m with a 0.24 m bottom-hinged flap. But in the currently studied case the inlet heigh was set at m. The outlet was placed in the floor below the ceiling, but it was only 4.8 m wide and 0.11 m long. 168
3 This, however, is small modification of the original geometry where the outlet is though a slatted floor located at the same place. The room was originally designed as extension of the previously studied geometry the Annex 20 test case and was designed for the creation of a two-dimensional recirculating flow. This room has previously been studied by Bjerg et al. (1999) using the commercial CFD code CFX with the k-i turbulence model with different boundary conditions (see Bjerg et al. for further information about experimental setup etc.). They did find a three-dimensional air flow pattern with the given geometry. The Reynolds number based on the inlet height was 4800, which gives an inlet velocity of 3.8 m/s. At first a steady state solution was computed by using Solve4kI, where the two-layer k-i turbulence model was employed. The iso-surface of velocity depicting the three-dimensional flow pattern is shown in Figure 6.90, where also the vector vectors in the return flow close to the floor is depicted. A hyperbolic tangent function was used to calculate a stretched grid. The minium grid cell dimension was selected as (Fx min /H, Fy min /H, Fz min /H) = (0.002, , 0.004), respectively. The number of grid points was The numerical setup and parameters were similar to the computational setup of the Annex 20 test. The kinetic energy at the inlet was prescribed to 4%. The three-dimensional flow pattern is clearly visible and the simulated direction of the return flow was 215 compared to measurement of 214 (Bjerg et al. (1999)) at (x, y, z) = (6.5, 0.2, 2.5). A return direction of 180 would have indicated a twodimensional return flow direction. Further grid refinement did not change the return flow direction. The inlet profile is assumed uniform and a 4% random fluctuation is superimposed. By mean of the computed k-i solution, a Large Eddy Simulation was started from here by using the implicit Solve4LES with the previously described dynamic one-equation subgrid scale model. The CFLnumber was kept at 1. One Large Eddy Turn Over Time (LETOT) was set to 2 L / U in Z5 sec. Figure 6.90: Iso-surface of velocity computed by Solve4kI by using the two layer k-i turbulence model with grid points. Light Green: 0.16 U in. The slice with velocity vector is located 1/15 H above the floor. Re =
4 A total of 200 LETOT were simulated to collect averaged values. This was done to ensure that any predominant effect from the k-i solution had diminished. In the next four figures different iso-surfaces of instantaneous velocity are depicted, thereby showing the fluctuation and separation of the inlet jet, at first closely to the side walls, and later the evolution of the three-dimensional flow pattern is shown. The flow begins to separate the side wall after only 1/3 of the room length. Figure 6.91: Instantaneous iso-surface of velocity = 0.37 computed by use of Solve4LES and the dynamic one-eqn. model sgs. Figure 6.92: Instantaneous iso-surface of velocity = 0.26 computed by use of Solve4LES and the dynamic one-eqn. model sgs. Figure 6.93: Instantaneous iso-surface of velocity = 0.21 computed by use of Solve4LES and the dynamic one-eqn. model sgs. Figure 6.94: Instantaneous iso-surface of velocity = 0.16 computed by use of Solve4LES and the dynamic one-eqn. model sgs. 170
5 Figure 6.95: Return flow direction of the monitor point in the symmetrical plane of the room (x, y, z) = (6.5, 0.2, 2.5). Direction of 180 indicate two-dimensional flow. In Figure 6.95 the instantaneous direction of the return flow is depicted 0.2 m above the floor. The averaged direction of the return flow was 217, which is close to the experimental data and the previous k-i solution. The fluctuation of the flow in the return direction indicates that the flow is unsteady. Although a switch-over of the inlet jet from the right downstream corner (right side wall) to the left downstream corner (left side wall) was captured during one measurement 1) it could not be reproduced during the Large Eddy Simulation. The change in the experiments from one side wall to the other one occurred without any obvious disturbances within the room. The period of switch-over was measured to sec. and compared to the time spent within the Large Eddy Simulation which represented less than 10%, this was not enough to capture any switch-over to occur. In Figure 6.96, the instantaneous velocity vectors for different planes perpendicular to the streamwise direction are depicted and coloured with the instantaneous streamwise velocity component. Vortical structures are visible just as in the Annex 20 test case. These vortical structures are generated due to the shear layer between the fluid within the inlet jet and the fluid outside the jet. Also, the friction at the wall and at the ceiling contributes to the generation of these secondary motions. Although this has not been identified, these secondary motions could contribute to the generation of the three-dimensional flow pattern within the room, which differs in geometry from the Annex 20 test case in form of a different 1 ) Personal communication with Senior Scientist Guoqiang Zhang, Danish Institute of Agricultural Sciences, Research Centre Bygholm. 171
6 inlet heigh, a larger width and a higher air flow rate. But also the difference in inlet height and the large width could contribute to the generation of the three-dimensional flow pattern. Figure 6.96: Instantaneous velocity vector plot at different locations in the SJVF-room: x/h = 1/3, 2/3, 4/3 and 15/6 computed using the implicit Solve4LES. The velocity vector has been projected on the yz-plane to illustrate secondary flow structures within the flow. Re = The velocity vectors are coloured with the streamwise instantaneous velocity component ( ). Figure 6.97: SJVF-room: Instantaneous velocity vectors coloured with normalized helicity within the xy-plane at z/h = 1/30 from left side wall. 172
7 In the Figures , the instantaneous velocity vectors in the xy-plane at different distances from the side wall are depicted. These velocity vectors are coloured with the normalized helicity. Figure 6.98: SJVF-room: Instantaneous velocity vectors coloured with normalized helicity within the xy-plane at z/h = 0.1 from left side wall. Figure 6.99: SJVF-room: Instantaneous velocity vectors coloured with normalized helicity within the xy-plane at z/h = 0.5 (symmetrical plane) from left side wall. Figure 6.100: SJVF-room: Instantaneous velocity vectors coloured with normalized helicity within the xy-plane at z/h = 0.9 from left side wall. 173
8 Helicity is a scalar quantity defined as an inner (dot) product of velocity and vorticity vectors. The vorticity q is defined as q = i. The normalized helicity is the cosine of the angle between the velocity and the vorticity vector. The extreme values of the normalized helicity, e.g. -1 and 1, indicate where the flow is highly three-dimensional, since there is alignment between the velocity and the vorticity vectors (i.e. the vorticity and the velocity vectors in the streamwise direction). The value zero means that the flow is basically two-dimensional, since the vectors are not aligned (vorticity is normal to the velocity). As a result, the normalized helicity value will be very close to 1.0 or -1.0 at the core of the secondary vortices perpendicular to the streamwise direction of the flow. This can also be utilized to locate the core of those streamwise vortices. When velocity vectors are coloured with normalized helicity, an indication of the three-dimensionality of the flow field can be studied. In the upper right corner a recirculating zone is clearly visible and also within the return flow there is a significantly high level of helicity that will indicate an unsteady flow pattern. Very close to the left side wall the angle between the velocity and vorticity vectors are large. In the symmetrical plane of the room extreme values of the angle can furthermore be observed both within the inlet jet and in the return flow area. Figure 6.101: Instantaneous iso-surface of velocity = 0.3 coloured with the normalized helicity to indicate the three-dimensionality of the flow for the Annex 20 test case. Implicit Solve4LES was used with grid 3 (Table 6.3) and the dynamic one-equation subgrid scale model. Red: indicates high level of helicity ( > 0.9); Green: middle level of helicity ( = 0.5); Blue: low level of helicity (< 0.1). 174
9 The level of helicity is compared to the Annex 20 test case depicted in Figure 6.101, which shows the instantaneous iso-surface of the streamwise velocity, coloured with the absolute value of the normalized helicity. For the inlet jet the level of helicity is at a moderate level, and as the back wall is approached the level of helicity is increased thus indicating that the flow is becoming more and more threedimensional and chaotic. The first part of the Annex 20 room also indicates that the flow pattern is almost two-dimensional. This was also confirmed by comparing the probability density for points in the symmetrical plane and close to the side walls. Finally in Figures show the velocity profile comparing the two-layer k-i turbulence within Solve4kI and the time averaged velocity profile from the Large Eddy Simulation with the dynamic one-equation subgrid scale model to the measured data (Bjerg et al. (1999)). Figure 6.102: Comparison between computed velocity profile by the two-layer k-i turbulence model and the time-averaged velocity from LES by using the dynamic one-equation sgs model close to the left side wall. (x/h=1.5; z/h = 1/6). Red: twolayer k-i turbulence mode; Blue: LES with Dynamic One-eqn. sgs. Markes: Measurement: Bjerg et al.(1999). Figure 6.103: Comparison between computed velocity profile by the two-layer k-i turbulence model and the time-averaged velocity from LES by using the dynamic one-equation sgs model at the symmetry plane. (x/h=1.5; z/h = 5/6). Red: twolayer k-i turbulence mode; Blue: LES with Dynamic One-eqn. sgs. Markes: Measurement: Bjerg et al.(1999). Just below the ceiling at y/h = 0.99 and x/h = 1.5, both types of simulation agree well with the measurement. Although for z/h = 3/2, the LES with the dynamic one-equation subgrid scale model predicts the jet peak velocity a little closer to the measurement. In the return flow the maximum velocity is also predicted quite well, compared to the experiments. The differences between the two-layer k-i turbulence model and the LES is observed in the middle of the room, where the low velocity region exists. Unfortunately, no experimental data exist for that area of the room. A very interesting point remains from these simulations when LES is used: Can LES predict the switch-over of the inlet jet? Since the present LES was only run for 1000 sec. and the measurements were more than 10 times longer, this was not attempted in this project. 175
10 Next, a closer look at the performance of the Solve4LES program is taken. Figure 6.104: Comparison between computed velocity profile by the two-layer k-i turbulence model and the time-averaged velocity from LES by use of the dynamic one-equation sgs model at the symmetry plane. (x/h = 1.5; z/h = 3/2). Red: two-layer k-i turbulence mode; Blue: LES with dynamic one-equation subgrid scale model. Markes: Measurement by Bjerg et al.(1999). 3.5 The performance of the explicit and implicit version of Solve4LES In the previous sections results from the use of different versions of the code Solve4LES are applied to compute transitional and fully turbulent flows inside ventilated enclosures by use of the Large Eddy Simulation. Since LES is a quite time-consuming method of dealing with turbulence modelling, because it computes the evolution of the flow field in time, an efficient numerical solver will be preferable in order to reduce the overall time demand for computing the solution. The basic numerical method in Solve4LES is the same as the Calc-LES code developed by Prof. Lars Davidson and others at Chalmers University, Sweden. But in the current project other numerical methods were tested to evaluate the overall performance of these methods and also in the search for the better numerical methods when using LES. Davidson & Nielsen (1996) and Emvin (1997) mention that around 80% of the total cpu-time required was spent on solving the pressure Poisson equations at each time-step when using an explicit formulation (see Chapter 4). As mentioned earlier in the present codes, the velocity component in the momentum equation can be solved by using either the Gauss-Seidel method, the Strongly Implicit Procedure method (SIP) (Stone, 1968) or the Modified Strongly Implicit Procedure (MSIP) (Schneider and Zedan, 1981) method. Davidson & Nielsen (1996) use the conjugate gradient (CG) method with incomplete Choleski preconditioning when solving the pressure Poisson equations within the explicit method. Preconditioning is a very important task for solving a system of linear equations by an iterative method, like CG and 176
11 other Krylov subspace methods in other to accelerate the convergence rate. During the implementation of the explicit version of Solve4LES the same Sparse Linear Algebra Package (SLAP) 2) as Davidson and Nielsen (1996) were tested briefly, simulating the Annex 20 test case. The CFL number was kept below 0.4, and only the Smagorinsky model was applied. The cpu-time per time step/grid point was found to (sec./grid point). The number of grid points used for this test were only in the x-, y- and z-directions, respectively, and the cpu-time increased dramatically when larger number of grid points were used. This means that 95% of the total cpu-time per time step will be required to solve the pressure Poisson equations, although this explicit method has the advantage that the matrix arising from the discretion of the pressure Poisson equations will have the same coefficients, which implies that the calculation of the preconditioning matrix only has to be carried out once within the first time step, by using the SLAP package for the calculation of the incomplete Choleski preconditioner required approximably the same time as 12 time steps during the Large Eddy Simulation. However, since this calculation will only be needed once, it represents only a fraction of the total computational time when using Large Eddy Simulation. In order to test the performance of other preconditioners, a separate implementation of the Conjugate Gradient method together with an incomplete Choleski (IC) preconditioner, the ILUT (Saad, 1996) which combines incomplete factorization and a threshold strategy for dropping the numerical fill-in values in the preconditioning matrix was implemented and tested. The dropping residual was set at 10-6 and number of fill-in s were set at 7 (see Saad (1996) for further information). Moreover, the multilevel-like incomplete factorization method by Botta &Wubs (1999) named, the Matrix-Renumbering Incomplete LU factorization method (MRILU), which is able to achieve near grid independent convergence rate as also the multigrid method is know for was implemented and tested. Finally, a multigrid method (Llorente and Melson, 1998) was implemented and tested for both the explicit and implicit versions of Solve4LES. Some other important factors to be considered when different numerical methods are implemented are the computer architectural issues. The testing of different storage formats for the sparse banded matrix arising from the pressure Poisson equations like those mentioned in the SPARSKIT 3) package can give performances between 4 and 52 Mflops on the currently used PC workstation. A technique to be used for improving performance on most workstations may be the re-use of the cache-memory, which will often require some sort of re-ordering of the storage-structure or the sparse-matrix bandwidth minimization. This, however, has not been utilized in the current code. Improvements to the code were made by carefull loop unrolling (the creation of fat loop), reduction of procedure, function calls within loops (by in-lining functions and even procedure calls), and the reduction of branches (conditional statements like if-then) within loops. The original matrix is stores as a banded sparse matrix and the preconditioning matrix using compressed row format (Saad, 1996). The same code running on vector supercomputers and cache-based RISC workstations will usually involve severe problems regarding the overall performance. This was also 2 SLAP: available from 3 ) SPARSKIT II available from 177
12 illustrated with the LESROOM code since this was optimized for the Cray vector supercomputer by means of a data structure that favours long vector lengths that could easily be pipelined by this kind of processor. Comparing the performance of LESROOM between the Cray and the PC based on the Mflops values, the ratio was between in favour of the Cray. Even if the Linpack 4) benchmark shows a ratio between the Cray and the currently used PC workstation to be 3.6 and 3, respectively, when the sizes of the dense system of linear equations are and For the explicit version of Solve4LES: in Table 6.8, the cpu-time per time step and grid point is displayed by using the different preconditioning methods, previously described. The Conjugate Gradient method with the MRILU preconditioner which was provided by Ass. Prof. Frank Wubs (see Acknowledgment) are one of the better methods including the multigrid method. The reason why the multigrid method is not the fastest, is that a the low number of grid points are used and that the number of coarsing levels within the multigrid method is limited, due to the number of grid points. Another major advantage of the last two mentioned methods are that they do not provide the same increase in cpu-time as the number of grid points is increased. An comparison between the different preconditioning methods used with the Conjugate Gradient solver and the multigrid method is depicted in Figure This clearly shows the advantage of using the CG with MRILU preconditioning and the multigrid method, because the convergence rate will be nearly grid independent. Furthermore, if the implementation of the currently used conjugate gradient is compared with the incomplete Choleski preconditioning to the same numerical method within the SLAP implementation, the first implementation will turn out to be more than 6 times faster on the same computer. Method of preconditioning for the pressure Poisson equations CPU Time / grid point [sec.] PC[Pentium II Xeon] Incomplete Choleski 4.01 l 10-5 ILUT 3.30 l 10-5 MRILU 2.11 l 10-5 Multigrid 2.14 l 10-5 Table 6.8: Comparison of cpu time per grid point for each time step the explicit Solve4LES when using different versions of preconditioning for the Conjugate Gradient solver and the multigrid solver. The second column represents a pc workstation with a Pentium II Xeon 450 Mhz and Windows NT 4.0 when using Digital Visual Fortran compiler version 6.0 and full compiler optimization. All computations were carried out by use of using grid points, and the CFL-number was kept at 0.4. The residual was set to Linpack benchmark: is used to solve a dense system of linear equations. It reflects the perfor-mance of a computer for solving a dense system of linear equations, and since the problem is very regular, the performance achieved will be quite high, and the performance numbers will give a good representation of the peak performance The sizes of the dense system of linear equations will usually be and See for additional information about the Linpack benchmark program and software. 178
13 Figure 6.105: Comparison between Conjugate Gradient solver with different preconditioning methods and a multigrid method in terms of total cpu-time as a function of number of grid points. Explicit version of SolveLES used on a PC workstation Pentium II 450 Mhz with the Smagorinsky sgs model The ratio between solving the pressure Poisson equation, the velocity components and the Smagorinsky sgs model was 80%, 15% and 5% of the total cpu time, respectively, obtained by using 64x32x32 grid points. Large Eddy Simulation Time / grid point [sec.] CPU requirements Subgrid Scale Model PC[Pentium II Xeon] relative to Smagrinsky sgs model Smagorinsky Dynamic with plane averaging Dynamic no averaging Dynamic One-Eqn Table 6.9: Comparison of cpu time per grid point for the explicit Solve4LES obtained by using different version of subgrid scale models.the second column represents an pc workstation with a Pentium II Xeon 450 Mhz and Windows NT 4.0 and using Digital Visual Fortran compiler version 6.0 and full compiler optimization. All computations were carried out by using grids points and the CFL-number was kept at 0.4. The residual was set to
14 Large Eddy Simulation. Time / grid point [sec.] CPU requirements relative to Smagorinsky Subgrid Scale Model PC[Pentium II Xeon] sgs model Smagorinsky Dynamic with plane averaging Dynamic no averaging Dynamic One-Eqn Table 6.10: Comparison of cpu time per grid point for the implicit Solve4LES by using different versions of subgrid scale models.the second column represents an pc workstation with a Pentium II Xeon 450 Mhz and Windows NT 4.0 and using Digital Visual Fortran compiler version 6.0and full compiler optimization. All computations were carried out by using grid points, and the CFL-number was kept at 1.0. The residual was set to 10-3, defined by the L 2 - normwhich requires 2-3 iterations within each time step. The cpu-time per grid point is taken as an average between 2 and 3 iterations. In Tables 6.9 and 6.10, the different subgrid scale models are compared in terms of cpu-time used in the explicit and implicit code. The reason why the more advanced dynamic subgrid scale models seem less expensive when the implicit code is used, is that more time is spent on solving the pressure and velocity component because of iterations within the implicit method. The subgrid scale model only needs to be solved once every time step. In the next 6 figures the cpu-time obtained by using the explicit and implicit Figure 6.106: CPU-time per time step obtained by using explicit Solve4LES with grid points and the Smagorinsky subgrid model. MSIP is used for velocity component and MG for pressure Poisson equations. Figure 6.107: CPU-time per time step obtained by using explicit Solve4LES with grid points and the dynamic subgrid model with plane averaging. MSIP is used for velocity component and MG for pressure Poisson equations. 180
15 Figure 6.108: CPU-time per time step obtained by Figure 6.109: CPU-time per time step obtained using explicit Solve4LES with grid by using implicit Solve4LES with grid points and the dynamic one-equation subgrid points and the Smagorinsky subgrid model. SIP is model. MSIP is used for velocity component and used for velocity component and MG for pressure MG for pressure Poisson equations. Poisson equations. versions of Solve4LES is depicted for each component, i.e. pressure, velocity and subgrid scale model. In the Figures , the total cpu-time per time step and for the other components only the cpu-time per implicit iteration are depicted. Therefore, variations in the implicit iteration are clearly depicted within the curve for the total cpu-time. Also the cpu-time for the u-velocity component is larger, due to some setup calculation that has to be performed before all the velocity component can be computed. The variations between the three different solvers for the velocity component were small. The SIP and MSIP were 1% and almost 3% faster, respectively, than the Gauss-Seidel method. Figure 6.110: CPU-time per time step obtained by using implicit Solve4LES with grid points and the dynamic subgrid model with plane averaging. SIP is used for velocity component and MG for pressure Poisson equations. Figure 6.111: CPU-time per time step obtained by using implicit Solve4LES with grid points and the dynamic one-eqn. subgrid model. SIP is used for velocity component and MG for pressure Poisson equations. 181
16 From the depicted results on the requirement of cpu-time when using the currently described implementation featuring either an explicit or an implicit method of time advancement, Large Eddy Simulation is still quite time-consuming compared to other methods of turbulence modelling where the flow field is decomposed into a steady and a fluctuating part, the Reynolds Averaging, although the implicit method has the potential to even larger time steps than the currently used limits of one. But it could be questionable to use large time steps in a transitional evolving flow field. This has so far not been tested. But LES is capable of providing much more information about the flow field, although not superior when mean averaging velocity profiles are compared. Furthermore, LES is capable of computing a solution to a flow field that is not fully developed, thus making it very attractive for many ventilation problems where the Reynolds number is often below 5000 and for many cases within the transitional flow range. Finally, the reported cpu measurement on the PC workstation with a Pentium II, Xeon processor, which is approximable one year behind the latest processor trend, could be reduced by upgrading to the latest processor. This modern processor will be nearly capable of reducing the cpu-time by a factor of two and a factor of almost four when the latest UNIX-based workstation with the Alpha processor is used. This assumption is based on the currently available specmark 5) ( Spec_fp95 ) numbers for floating point speeds. The additional use of parallel computer would certainly make Large Eddy Simulations attractive. But the interaction between implementation, the numerical methods, the architectural issues, and the compiler technology remains a challenge. The easy gains in computational speed are over. 5 SPECMARK: see for further information. 182
Assessment of the numerical solver
Chapter 5 Assessment of the numerical solver In this chapter the numerical methods described in the previous chapter are validated and benchmarked by applying them to some relatively simple test cases
More informationModelling of a Wall Inlet in Numerical Simulation of Airflow in Livestock Buildings
1 Modelling of a Wall Inlet in Numerical Simulation of Airflow in Livestock Buildings B. Bjerg 1, K. Svidt 2, S Morsing 3, G. Zhang 3 and J. O. Johnsen 3 1 The Royal Veterinary and Agricultural University,
More informationBackward facing step Homework. Department of Fluid Mechanics. For Personal Use. Budapest University of Technology and Economics. Budapest, 2010 autumn
Backward facing step Homework Department of Fluid Mechanics Budapest University of Technology and Economics Budapest, 2010 autumn Updated: October 26, 2010 CONTENTS i Contents 1 Introduction 1 2 The problem
More informationSimulation of Turbulent Flow over the Ahmed Body
Simulation of Turbulent Flow over the Ahmed Body 58:160 Intermediate Mechanics of Fluids CFD LAB 4 By Timur K. Dogan, Michael Conger, Maysam Mousaviraad, and Fred Stern IIHR-Hydroscience & Engineering
More informationLarge Eddy Simulation of Flow over a Backward Facing Step using Fire Dynamics Simulator (FDS)
The 14 th Asian Congress of Fluid Mechanics - 14ACFM October 15-19, 2013; Hanoi and Halong, Vietnam Large Eddy Simulation of Flow over a Backward Facing Step using Fire Dynamics Simulator (FDS) Md. Mahfuz
More informationVerification and Validation of Turbulent Flow around a Clark-Y Airfoil
Verification and Validation of Turbulent Flow around a Clark-Y Airfoil 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 2 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University
More informationIntroduction to CFX. Workshop 2. Transonic Flow Over a NACA 0012 Airfoil. WS2-1. ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.
Workshop 2 Transonic Flow Over a NACA 0012 Airfoil. Introduction to CFX WS2-1 Goals The purpose of this tutorial is to introduce the user to modelling flow in high speed external aerodynamic applications.
More informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationSimulation of Flow Development in a Pipe
Tutorial 4. Simulation of Flow Development in a Pipe Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D turbulent fluid flow in a pipe. The pipe networks are common
More informationIntroduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
More informationModeling Unsteady Compressible Flow
Tutorial 4. Modeling Unsteady Compressible Flow Introduction In this tutorial, FLUENT s density-based implicit solver is used to predict the timedependent flow through a two-dimensional nozzle. As an initial
More informationCFD MODELING FOR PNEUMATIC CONVEYING
CFD MODELING FOR PNEUMATIC CONVEYING Arvind Kumar 1, D.R. Kaushal 2, Navneet Kumar 3 1 Associate Professor YMCAUST, Faridabad 2 Associate Professor, IIT, Delhi 3 Research Scholar IIT, Delhi e-mail: arvindeem@yahoo.co.in
More informationStream Function-Vorticity CFD Solver MAE 6263
Stream Function-Vorticity CFD Solver MAE 66 Charles O Neill April, 00 Abstract A finite difference CFD solver was developed for transient, two-dimensional Cartesian viscous flows. Flow parameters are solved
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationNumerical Study of Turbulent Flow over Backward-Facing Step with Different Turbulence Models
Numerical Study of Turbulent Flow over Backward-Facing Step with Different Turbulence Models D. G. Jehad *,a, G. A. Hashim b, A. K. Zarzoor c and C. S. Nor Azwadi d Department of Thermo-Fluids, Faculty
More informationANSYS AIM Tutorial Turbulent Flow Over a Backward Facing Step
ANSYS AIM Tutorial Turbulent Flow Over a Backward Facing Step Author(s): Sebastian Vecchi, ANSYS Created using ANSYS AIM 18.1 Problem Specification Pre-Analysis & Start Up Governing Equation Start-Up Geometry
More informationMesh techniques and uncertainty for modelling impulse jetfans. O. A. (Sam) Alshroof. CFD manager Olsson Fire and Risk
Mesh techniques and uncertainty for modelling impulse jetfans O. A. (Sam) Alshroof CFD manager Olsson Fire and Risk Email: Sam.Alshroof@olssonfire.com Abstract This study presents the numerical modelling
More informationTutorial 17. Using the Mixture and Eulerian Multiphase Models
Tutorial 17. Using the Mixture and Eulerian Multiphase Models Introduction: This tutorial examines the flow of water and air in a tee junction. First you will solve the problem using the less computationally-intensive
More informationVerification of Laminar and Validation of Turbulent Pipe Flows
1 Verification of Laminar and Validation of Turbulent Pipe Flows 1. Purpose ME:5160 Intermediate Mechanics of Fluids CFD LAB 1 (ANSYS 18.1; Last Updated: Aug. 1, 2017) By Timur Dogan, Michael Conger, Dong-Hwan
More informationRANS COMPUTATION OF RIBBED DUCT FLOW USING FLUENT AND COMPARING TO LES
RANS COMPUTATION OF RIBBED DUCT FLOW USING FLUENT AND COMPARING TO LES Máté M., Lohász +*& / Ákos Csécs + + Department of Fluid Mechanics, Budapest University of Technology and Economics, Budapest * Von
More informationLab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.
More informationModeling External Compressible Flow
Tutorial 3. Modeling External Compressible Flow Introduction The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a nonzero angle of attack. You will use the Spalart-Allmaras
More informationSimulation of Laminar Pipe Flows
Simulation of Laminar Pipe Flows 57:020 Mechanics of Fluids and Transport Processes CFD PRELAB 1 By Timur Dogan, Michael Conger, Maysam Mousaviraad, Tao Xing and Fred Stern IIHR-Hydroscience & Engineering
More informationReproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software
Reports of Research Institute for Applied Mechanics, Kyushu University No.150 (71 83) March 2016 Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software Report 3: For the Case
More informationCalculate a solution using the pressure-based coupled solver.
Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings
More informationCompressible Flow in a Nozzle
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 1 Compressible Flow in a Nozzle Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification Consider air flowing at high-speed through a
More informationAerodynamic Study of a Realistic Car W. TOUGERON
Aerodynamic Study of a Realistic Car W. TOUGERON Tougeron CFD Engineer 2016 Abstract This document presents an aerodynamic CFD study of a realistic car geometry. The aim is to demonstrate the efficiency
More informationMultiphase flow metrology in oil and gas production: Case study of multiphase flow in horizontal tube
Multiphase flow metrology in oil and gas production: Case study of multiphase flow in horizontal tube Deliverable 5.1.2 of Work Package WP5 (Creating Impact) Authors: Stanislav Knotek Czech Metrology Institute
More informationPossibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics
Possibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics Masanori Hashiguchi 1 1 Keisoku Engineering System Co., Ltd. 1-9-5 Uchikanda, Chiyoda-ku,
More informationTurbulent Premixed Combustion with Flamelet Generated Manifolds in COMSOL Multiphysics
Turbulent Premixed Combustion with Flamelet Generated Manifolds in COMSOL Multiphysics Rob J.M Bastiaans* Eindhoven University of Technology *Corresponding author: PO box 512, 5600 MB, Eindhoven, r.j.m.bastiaans@tue.nl
More informationAdvanced ANSYS FLUENT Acoustics
Workshop Modeling Flow-Induced (Aeroacoustic) Noise 14.5 Release Advanced ANSYS FLUENT Acoustics 2011 ANSYS, Inc. November 7, 2012 1 Introduction This tutorial demonstrates how to model 2D turbulent flow
More informationNUMERICAL SIMULATION OF FLOW FIELD IN AN ANNULAR TURBINE STATOR WITH FILM COOLING
24 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES NUMERICAL SIMULATION OF FLOW FIELD IN AN ANNULAR TURBINE STATOR WITH FILM COOLING Jun Zeng *, Bin Wang *, Yong Kang ** * China Gas Turbine Establishment,
More informationReproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software
Reports of Research Institute for Applied Mechanics, Kyushu University, No.150 (60-70) March 2016 Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software Report 2: For the Case
More informationVerification and Validation of Turbulent Flow around a Clark-Y Airfoil
1 Verification and Validation of Turbulent Flow around a Clark-Y Airfoil 1. Purpose ME:5160 Intermediate Mechanics of Fluids CFD LAB 2 (ANSYS 19.1; Last Updated: Aug. 7, 2018) By Timur Dogan, Michael Conger,
More informationShape optimisation using breakthrough technologies
Shape optimisation using breakthrough technologies Compiled by Mike Slack Ansys Technical Services 2010 ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary Introduction Shape optimisation technologies
More informationIntroduction to C omputational F luid Dynamics. D. Murrin
Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena
More informationON THE NUMERICAL MODELING OF IMPINGING JET HEAT TRANSFER
ON THE NUMERICAL MODELING OF IMPINGING JET HEAT TRANSFER Mirko Bovo 1,2, Sassan Etemad 2 and Lars Davidson 1 1 Dept. of Applied Mechanics, Chalmers University of Technology, Gothenburg, Sweden 2 Powertrain
More informationDES Turbulence Modeling for ICE Flow Simulation in OpenFOAM
2 nd Two-day Meeting on ICE Simulations Using OpenFOAM DES Turbulence Modeling for ICE Flow Simulation in OpenFOAM V. K. Krastev 1, G. Bella 2 and G. Campitelli 1 University of Tuscia, DEIM School of Engineering
More informationA Study of the Development of an Analytical Wall Function for Large Eddy Simulation of Turbulent Channel and Rectangular Duct Flow
University of Wisconsin Milwaukee UWM Digital Commons Theses and Dissertations August 2014 A Study of the Development of an Analytical Wall Function for Large Eddy Simulation of Turbulent Channel and Rectangular
More informationCoupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić
Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions Milovan Perić Contents The need to couple STAR-CCM+ with other theoretical or numerical solutions Coupling approaches: surface and volume
More informationComputational Fluid Dynamics as an advanced module of ESP-r Part 1: The numerical grid - defining resources and accuracy. Jordan A.
Computational Fluid Dynamics as an advanced module of ESP-r Part 1: The numerical grid - defining resources and accuracy Jordan A. Denev Abstract: The present paper is a first one from a series of papers
More informationDetached Eddy Simulation Analysis of a Transonic Rocket Booster for Steady & Unsteady Buffet Loads
Detached Eddy Simulation Analysis of a Transonic Rocket Booster for Steady & Unsteady Buffet Loads Matt Knapp Chief Aerodynamicist TLG Aerospace, LLC Presentation Overview Introduction to TLG Aerospace
More informationNumerical Simulation Study on Aerodynamic Characteristics of the High Speed Train under Crosswind
2017 2nd International Conference on Industrial Aerodynamics (ICIA 2017) ISBN: 978-1-60595-481-3 Numerical Simulation Study on Aerodynamic Characteristics of the High Speed Train under Crosswind Fan Zhao,
More informationNon-Newtonian Transitional Flow in an Eccentric Annulus
Tutorial 8. Non-Newtonian Transitional Flow in an Eccentric Annulus Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D, turbulent flow of a non-newtonian fluid. Turbulent
More informationSimulation of Turbulent Flow in an Asymmetric Diffuser
Simulation of Turbulent Flow in an Asymmetric Diffuser 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 3 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University of Iowa C.
More informationSimulation of Turbulent Flow around an Airfoil
1. Purpose Simulation of Turbulent Flow around an Airfoil ENGR:2510 Mechanics of Fluids and Transfer Processes CFD Lab 2 (ANSYS 17.1; Last Updated: Nov. 7, 2016) By Timur Dogan, Michael Conger, Andrew
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationQNET-CFD WIKI KNOWLEDGE BASE UNDERLYING FLOW REGIME DOCUMENT TEMPLATE
QNET-CFD WIKI KNOWLEDGE BASE UNDERLYING FLOW REGIME DOCUMENT TEMPLATE UNDERLYING FLOW REGIME DOCUMENT GUIDANCE NOTES You are first asked to provide a brief review of the state of the art for this UFR,
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationCOMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS
COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS Dr W. Malalasekera Version 3.0 August 2013 1 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE
More informationSimulation and Validation of Turbulent Pipe Flows
Simulation and Validation of Turbulent Pipe Flows ENGR:2510 Mechanics of Fluids and Transport Processes CFD LAB 1 (ANSYS 17.1; Last Updated: Oct. 10, 2016) By Timur Dogan, Michael Conger, Dong-Hwan Kim,
More informationNumerical Simulation of Coastal Wave Processes with the Use of Smoothed Particle Hydrodynamics (SPH) Method
Aristotle University of Thessaloniki Faculty of Engineering Department of Civil Engineering Division of Hydraulics and Environmental Engineering Laboratory of Maritime Engineering Christos V. Makris Dipl.
More informationInvestigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)
Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,
More informationNIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011
NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011 First-Order Hyperbolic System Method If you have a CFD book for hyperbolic problems, you have a CFD book for all problems.
More informationTeam 194: Aerodynamic Study of Airflow around an Airfoil in the EGI Cloud
Team 194: Aerodynamic Study of Airflow around an Airfoil in the EGI Cloud CFD Support s OpenFOAM and UberCloud Containers enable efficient, effective, and easy access and use of MEET THE TEAM End-User/CFD
More informationShape optimization and active flow control for improved aerodynamic properties Siniša Krajnovic
Shape optimization and active flow control for improved aerodynamic properties Siniša Krajnovic HPC resources used Computer resources at C3SE at Chalmers in Göteborg Computer cluster: Neolith NSC Linköping
More informationCFD Analysis of a Fully Developed Turbulent Flow in a Pipe with a Constriction and an Obstacle
CFD Analysis of a Fully Developed Turbulent Flow in a Pipe with a Constriction and an Obstacle C, Diyoke Mechanical Engineering Department Enugu State University of Science & Tech. Enugu, Nigeria U, Ngwaka
More informationGrid. Apr 09, 1998 FLUENT 5.0 (2d, segregated, lam) Grid. Jul 31, 1998 FLUENT 5.0 (2d, segregated, lam)
Tutorial 2. Around an Airfoil Transonic Turbulent Flow Introduction: The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a non-zero angle of attack. You will use the
More informationAn Introduction to SolidWorks Flow Simulation 2010
An Introduction to SolidWorks Flow Simulation 2010 John E. Matsson, Ph.D. SDC PUBLICATIONS www.sdcpublications.com Schroff Development Corporation Chapter 2 Flat Plate Boundary Layer Objectives Creating
More informationA B C D E. Settings Choose height, H, free stream velocity, U, and fluid (dynamic viscosity and density ) so that: Reynolds number
Individual task Objective To derive the drag coefficient for a 2D object, defined as where D (N/m) is the aerodynamic drag force (per unit length in the third direction) acting on the object. The object
More informationTUTORIAL#3. Marek Jaszczur. Boundary Layer on a Flat Plate W1-1 AGH 2018/2019
TUTORIAL#3 Boundary Layer on a Flat Plate Marek Jaszczur AGH 2018/2019 W1-1 Problem specification TUTORIAL#3 Boundary Layer - on a flat plate Goal: Solution for boudary layer 1. Creating 2D simple geometry
More informationUsing a Single Rotating Reference Frame
Tutorial 9. Using a Single Rotating Reference Frame Introduction This tutorial considers the flow within a 2D, axisymmetric, co-rotating disk cavity system. Understanding the behavior of such flows is
More informationDevelopment of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics
Development of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics I. Pantle Fachgebiet Strömungsmaschinen Karlsruher Institut für Technologie KIT Motivation
More informationThe Spalart Allmaras turbulence model
The Spalart Allmaras turbulence model The main equation The Spallart Allmaras turbulence model is a one equation model designed especially for aerospace applications; it solves a modelled transport equation
More informationCFD Application in Offshore Structures Design at PETROBRAS
CFD Application in Offshore Structures Design at PETROBRAS Marcus Reis ESSS CFD Director Mooring System Design of Floating Production Systems; Current and Wind Loads; Wave Induced Drag Coefficients. Case
More informationAalborg Universitet. Numerical 3-D Modelling of Overflows Larsen, Torben; Nielsen, L.; Jensen, B.; Christensen, E. D.
Aalborg Universitet Numerical 3-D Modelling of Overflows Larsen, Torben; Nielsen, L.; Jensen, B.; Christensen, E. D. Published in: Confernce Proceedings : 11th International Conference on Urban Drainage
More informationThe viscous forces on the cylinder are proportional to the gradient of the velocity field at the
Fluid Dynamics Models : Flow Past a Cylinder Flow Past a Cylinder Introduction The flow of fluid behind a blunt body such as an automobile is difficult to compute due to the unsteady flows. The wake behind
More informationCFD VALIDATION FOR SURFACE COMBATANT 5415 STRAIGHT AHEAD AND STATIC DRIFT 20 DEGREE CONDITIONS USING STAR CCM+
CFD VALIDATION FOR SURFACE COMBATANT 5415 STRAIGHT AHEAD AND STATIC DRIFT 20 DEGREE CONDITIONS USING STAR CCM+ by G. J. Grigoropoulos and I..S. Kefallinou 1. Introduction and setup 1. 1 Introduction The
More informationNUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING
Review of the Air Force Academy No.3 (35)/2017 NUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING Cvetelina VELKOVA Department of Technical Mechanics, Naval Academy Nikola Vaptsarov,Varna, Bulgaria (cvetelina.velkova1985@gmail.com)
More informationPotsdam Propeller Test Case (PPTC)
Second International Symposium on Marine Propulsors smp 11, Hamburg, Germany, June 2011 Workshop: Propeller performance Potsdam Propeller Test Case (PPTC) Olof Klerebrant Klasson 1, Tobias Huuva 2 1 Core
More informationSolidWorks Flow Simulation 2014
An Introduction to SolidWorks Flow Simulation 2014 John E. Matsson, Ph.D. SDC PUBLICATIONS Better Textbooks. Lower Prices. www.sdcpublications.com Powered by TCPDF (www.tcpdf.org) Visit the following websites
More informationCFD Best Practice Guidelines: A process to understand CFD results and establish Simulation versus Reality
CFD Best Practice Guidelines: A process to understand CFD results and establish Simulation versus Reality Judd Kaiser ANSYS Inc. judd.kaiser@ansys.com 2005 ANSYS, Inc. 1 ANSYS, Inc. Proprietary Overview
More informationSimulation of Turbulent Flow over the Ahmed Body
1 Simulation of Turbulent Flow over the Ahmed Body ME:5160 Intermediate Mechanics of Fluids CFD LAB 4 (ANSYS 18.1; Last Updated: Aug. 18, 2016) By Timur Dogan, Michael Conger, Dong-Hwan Kim, Maysam Mousaviraad,
More informationExample 24 Spring-back
Example 24 Spring-back Summary The spring-back simulation of sheet metal bent into a hat-shape is studied. The problem is one of the famous tests from the Numisheet 93. As spring-back is generally a quasi-static
More informationProfile Catalogue for Airfoil Sections Based on 3D Computations
Risø-R-58(EN) Profile Catalogue for Airfoil Sections Based on 3D Computations Franck Bertagnolio, Niels N. Sørensen and Jeppe Johansen Risø National Laboratory Roskilde Denmark December 26 Author: Franck
More informationCOMPUTATIONAL FLUID DYNAMICS USED IN THE DESIGN OF WATERBLAST TOOLING
2015 WJTA-IMCA Conference and Expo November 2-4 New Orleans, Louisiana Paper COMPUTATIONAL FLUID DYNAMICS USED IN THE DESIGN OF WATERBLAST TOOLING J. Schneider StoneAge, Inc. Durango, Colorado, U.S.A.
More informationUsage of CFX for Aeronautical Simulations
Usage of CFX for Aeronautical Simulations Florian Menter Development Manager Scientific Coordination ANSYS Germany GmbH Overview Elements of CFD Technology for aeronautical simulations: Grid generation
More informationComputational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent
MEGR 7090-003, Computational Fluid Dynamics :1 7 Spring 2015 Computational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent Rahul R Upadhyay Master of Science, Dept of Mechanical
More informationIntroduction to Parallel Computing
Introduction to Parallel Computing W. P. Petersen Seminar for Applied Mathematics Department of Mathematics, ETHZ, Zurich wpp@math. ethz.ch P. Arbenz Institute for Scientific Computing Department Informatik,
More informationThis tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following:
Tutorial 22. Modeling Solidification Introduction This tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following: Define a
More informationFLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016
FLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016 Note: These instructions are based on an older version of FLUENT, and some of the instructions
More informationFlow in an Intake Manifold
Tutorial 2. Flow in an Intake Manifold Introduction The purpose of this tutorial is to model turbulent flow in a simple intake manifold geometry. An intake manifold is a system of passages which carry
More informationMultigrid Solvers in CFD. David Emerson. Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK
Multigrid Solvers in CFD David Emerson Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK david.emerson@stfc.ac.uk 1 Outline Multigrid: general comments Incompressible
More informationTHE APPLICATION OF AN ATMOSPHERIC BOUNDARY LAYER TO EVALUATE TRUCK AERODYNAMICS IN CFD
THE APPLICATION OF AN ATMOSPHERIC BOUNDARY LAYER TO EVALUATE TRUCK AERODYNAMICS IN CFD A SOLUTION FOR A REAL-WORLD ENGINEERING PROBLEM Ir. Niek van Dijk DAF Trucks N.V. CONTENTS Scope & Background Theory:
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
More informationProblem description. The FCBI-C element is used in the fluid part of the model.
Problem description This tutorial illustrates the use of ADINA for analyzing the fluid-structure interaction (FSI) behavior of a flexible splitter behind a 2D cylinder and the surrounding fluid in a channel.
More informationDrivAer-Aerodynamic Investigations for a New Realistic Generic Car Model using ANSYS CFD
DrivAer-Aerodynamic Investigations for a New Realistic Generic Car Model using ANSYS CFD Thomas Frank (*), BenediktGerlicher (*), Juan Abanto (**) (*) ANSYS Germany, Otterfing, Germany (**) ANSYS Inc.,
More informationLaser speckle based background oriented schlieren measurements in a fire backlayering front
Laser speckle based background oriented schlieren measurements in a fire backlayering front Philipp Bühlmann 1*, Alexander H. Meier 1, Martin Ehrensperger 1, Thomas Rösgen 1 1: ETH Zürich, Institute of
More informationExperimental and Numerical Analysis of Near Wall Flow at the Intake Valve and its Influence on Large-Scale Fluctuations
Experimental and Numerical Analysis of Near Wall Flow at the Intake Valve and its Influence on Large-Scale Fluctuations Frank Hartmann, Stefan Buhl, Florian Gleiß, Christian Hasse Philipp Barth, Martin
More informationCFD PREDICTION OF WIND PRESSURES ON CONICAL TANK
CFD PREDICTION OF WIND PRESSURES ON CONICAL TANK T.A.Sundaravadivel a, S.Nadaraja Pillai b, K.M.Parammasivam c a Lecturer, Dept of Aeronautical Engg, Satyabama University, Chennai, India, aerovelu@yahoo.com
More informationNUMERICAL ANALYSIS OF WIND EFFECT ON HIGH-DENSITY BUILDING AERAS
NUMERICAL ANALYSIS OF WIND EFFECT ON HIGH-DENSITY BUILDING AERAS Bin ZHAO, Ying LI, Xianting LI and Qisen YAN Department of Thermal Engineering, Tsinghua University Beijing, 100084, P.R. China ABSTRACT
More informationTutorial 2. Modeling Periodic Flow and Heat Transfer
Tutorial 2. Modeling Periodic Flow and Heat Transfer Introduction: Many industrial applications, such as steam generation in a boiler or air cooling in the coil of an air conditioner, can be modeled as
More informationMOMENTUM AND HEAT TRANSPORT INSIDE AND AROUND
MOMENTUM AND HEAT TRANSPORT INSIDE AND AROUND A CYLINDRICAL CAVITY IN CROSS FLOW G. LYDON 1 & H. STAPOUNTZIS 2 1 Informatics Research Unit for Sustainable Engrg., Dept. of Civil Engrg., Univ. College Cork,
More informationLES Applications in Aerodynamics
LES Applications in Aerodynamics Kyle D. Squires Arizona State University Tempe, Arizona, USA 2010 Tutorial School on Fluid Dynamics: Topics in Turbulence Center for Scientific Computation and Mathematical
More informationStrömningslära Fluid Dynamics. Computer laboratories using COMSOL v4.4
UMEÅ UNIVERSITY Department of Physics Claude Dion Olexii Iukhymenko May 15, 2015 Strömningslära Fluid Dynamics (5FY144) Computer laboratories using COMSOL v4.4!! Report requirements Computer labs must
More informationTransition Flow and Aeroacoustic Analysis of NACA0018 Satish Kumar B, Fred Mendonç a, Ghuiyeon Kim, Hogeon Kim
Transition Flow and Aeroacoustic Analysis of NACA0018 Satish Kumar B, Fred Mendonç a, Ghuiyeon Kim, Hogeon Kim Transition Flow and Aeroacoustic Analysis of NACA0018 Satish Kumar B, Fred Mendonç a, Ghuiyeon
More informationDetached-Eddy Simulation of a Linear Compressor Cascade with Tip Gap and Moving Wall *)
FOI, Stockholm, Sweden 14-15 July, 2005 Detached-Eddy Simulation of a Linear Compressor Cascade with Tip Gap and Moving Wall *) A. Garbaruk,, M. Shur, M. Strelets, and A. Travin *) Study is carried out
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY. Analyzing wind flow around the square plate using ADINA Project. Ankur Bajoria
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Analyzing wind flow around the square plate using ADINA 2.094 - Project Ankur Bajoria May 1, 2008 Acknowledgement I would like to thank ADINA R & D, Inc for the full
More informationCFD modelling of thickened tailings Final project report
26.11.2018 RESEM Remote sensing supporting surveillance and operation of mines CFD modelling of thickened tailings Final project report Lic.Sc.(Tech.) Reeta Tolonen and Docent Esa Muurinen University of
More informationSTAR-CCM+: Wind loading on buildings SPRING 2018
STAR-CCM+: Wind loading on buildings SPRING 2018 1. Notes on the software 2. Assigned exercise (submission via Blackboard; deadline: Thursday Week 3, 11 pm) 1. NOTES ON THE SOFTWARE STAR-CCM+ generates
More information