Numerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
|
|
- Horace Whitehead
- 5 years ago
- Views:
Transcription
1 Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin Chu, aiwan * yylung@chu.edu.tw Abstract. he periodic vortex shedding past a cylinder is presented. he program was validated by the solution of one-dimensional Stokes second problem first. hen, a two-dimensional steady flow through a cylinder at Reynolds number.54 was used to compare with the streamlines with experiment and also at Reynolds number 28.4, to check the separation point and the length of wake. At Reynolds number 280 where the wake is laminar, the solution becomes un-symmetric after 25 cycles of iterations automatically and becomes fully periodic flow after 56 cycles of iterations. For Reynolds number 000 with a turbulent wake, the numerical solution starts to un-symmetric after 2 cycles of iterations, and becomes fully periodic after 3 cycles of iterations. he streamlines of the periodic solution are compared very well with Arnone s result. Using 3500 steps in a cycle, the error on the Strouhal number is 9%. However, increasing to steps in a cycle, the error on the Strouhal number is less than %. Wilcox s high and low Reynolds number κ-ω models, Jones-Launder κ-ε model and Launder-Sharma κ-ε model were used to study the structure of turbulent vortex street at Reynolds number he high Reynolds number κ-ω model gives the largest eddy viscosity among these turbulent models. On the other hands, the solution using Jones-Launder κ-ε model was similar to laminar one. he spreading angle and the distance between trailing vortex are better simulated by using high Reynolds number κ-ω model, Besides, based on the accuracy of the predicted Strouhal number, the Wilcox s high Reynolds number κ-ω models is suggested. Keywords:Stokes second problem, Strouhal number, wall-distance-free turbulence models. Introduction he unsteady vortex shedding behind a body causes hysteric forces on the structure. If the frequency of the unsteady loading closes to the natural frequency of the structure, the vortex could induce severe structure vibration. he catastrophe of acoma Bridge is a famous example. However, the relationship between the natural unsteadiness and the evolution of a symmetric numerical solution into a cyclic solution is not well posed in the past. Hence, in this study a numerical program to reproduce the natural unsteadiness is an important demonstration of a robust scheme. Lienhard [] has summarized the vortex pattern behind a cylinder. For Reynolds number below 5, no separation is occurred. For Reynolds number between 5 and 40, there is a fixed pair of vortices in the wake. For Reynolds number between 40 and 300, the Karman vortex street is laminar. For Reynolds number between 300 and 3x0 5, the vortex street is fully turbulent. For Reynolds number between 3x0 5 and 3.5x0 6, the laminar boundary layer has undergone turbulent transition. he wake is narrower and disorganized. No vortex street is apparent. For Reynolds number above 3.5x0 6, the boundary layer is turbulent and the wake is thinner. Kimura and sutahara [2] have arranged a groove on the cylinder to control the cylinder separation. Placing the groove at 84 degree from the leading edge can minimize the drag coefficient of the cylinder. Bearman and Harvey [3] also studied the effect of dimples on the cylinder to trickle the boundary from laminar to turbulent. An effective reduction on the drag coefficient of the cylinder between Reynolds number 4x0 4 and 3x0 5 has been achieved. he numerical simulation of unsteady flow field has been addressed extensively. Jameson [4] extended the multigrid methodology to allow practical time steps with the use of an explicit code. Arnone et al. [5] expanded the Jameson s subiterative technique into a dual time stepping. A good prediction on the Strouhal number was shown in their work. Rumsey et al. [6] have explored the accuracy and efficiency of two different types of subiterations in both explicit and implicit codes for unsteady flows. he pseudo time subiterations removes the time-step limitations in practice. A pseudo time version of the code may be used to speed up the solution. Because the primary focus of this paper is on the turbulent modeling. he efficiency of the subiterative time-stepping method is not performed for the circular cylinder flow. he solution is marched by local time-stepping at the beginning till the flow field becomes unsymmetric. hen a prescribed time step is used for every cell, several cycles of calculations is proceeded to ensure the flow field is periodic [0].
2 he large region of separated flow behind the cylinder gives a challenge on the turbulent modeling. he wall-distance-free turbulent models may be more promising for calculations in complex flow field. In this paper, results with the Jones and Launder [7] κ-ε model, Launder and Sharma [8] κ-ε model, and Wilcox s high and low Reynolds number [9] κ-ω models are compared for the unsteady wake behind the cylinder. he unsteady hysteresis effect is reproduced computationally using these turbulence models. 2. Mathematical modeling he two-dimensional unsteady Reynolds-Averaged Navier-Stokes equations in conservative form are written as: U Ec Fc Ev Fv + + = + () t x y x y where ρ ρu ρv ρu 2 U = u p E ρ + ρuv ρv c = F ρuv c = 2 ρv + p ρ E ρue ( + p/ ρ) ρve ( + p/ ρ ) 0 0 τ xx τ xy E v = τ xy F v = τ yy uτ xx + vτ xy + K x uτ xy vτ yy K + + y he pressure is obtained by the equation of state for a perfect gas: 2 2 p = ( γ ) ρ[ E ( u + v ) / 2] (2) he stress tensor has been modeled using the Boussinesq approximation: τ = ( μ+ μ )[ u x + u x 2/3 u x δ ] 2/ 3ρk δ ij t ( i j j i) k k ij ij (3) and μt is the eddy-viscosity. A time-marching solver has been employed for the solution of the equations. he solver uses a fourth-order Runge-Kutta scheme. he cell-centered averaged scheme is used for the discretisation of the inviscid and viscous fluxes. A fouth-order artificial viscosity is added to reduce numerical oscillation. he solution is marched by local time-stepping at the beginning till the flow field becomes un-symmetric. hen a prescribed time step is used for every cell, several cycles of calculations is proceeded to ensure the flow field is periodic. he accuracy of the numerical solution can be checked by its Strouhal mumber nd St = (4) U where n is the frequency of the vortex shedding D is the diameter of the cylinder and U is the free stream velocity. he structure of the unsteady wake is more easily to present through the entropy distribution, the definition of the entropy is P s = C p ln Rln (5) P where C is the specific heat at constant pressure. p 3. urbulence models In this paper, four different wall-distance-free turbulent models of κ-ε and κ-ω have been assessed. hese are: () Jones-Launder (JL) κ-ε (2) Launder-Sharma (LS) κ-ε (3) Wilcox high-reynolds number (κ-ωh) κ-ω (4) Wilcox low-reynolds number (κ-ωl) κ-ω All models employ the boundary condition for the turbulent kinetic energy at solid boundaries k w = 0 (6) For the JL and LS models the boundary condition for the turbulent dissipative rate is ε w = 0 (7) In the case of κ-ωh and κ-ωl models, a zero gradient condition for specific dissipative rate on the solid wall was used. At the far field of the computational domain, the turbulent kinetic energy, turbulent dissipative rate, and specific dissipative rate are propagated outward. 4. Results he numerical program was validated by a one-dimensional Stokes second problem first. Figure (a) shows the velocity distribution at the beginning of the sine wave (t/=0) using three different numbers of grids with 40 intervals in a cycle. he maximum error from the analytical solution is less than 0.04% by using 2 grid points and only 0.023% by using 8 grid points. Figure (b) gives the velocity distribution at the quarter of the sine wave (t/=0.25). he agreement is very well again. able and able 2 give the maximum numerical error at t/=0 and t/=0.25 respectively. From this study, the accuracy of the solution is more sensitive to grid resolution. Figure 2(a) shows the computational O-grid of the circular cylinder flow and Figure 2(b) plots the two symmetrical monitor points on the back of the cylinder. here are 20 grids around the cylinder and2 grids in the normal direction. he grid point extends 5 diameters upstream and 0 diameters downstream. 2
3 he first calculation is performed for a Reynolds number of.54. Figure 3 shows the streamlines between the calculation and the experiment []. here is no separation in this case. he agreement in the streamlines is very good. Another validation of the program has been performed for the Reynolds number at Figure 4 plots the instantaneous pressure variation of the monitor points. No un-symmetric solution was found for this Reynolds number. Figure 5 shows the streamlines of the re-circulating region between numerical solution and experiment []. he agreement on the separation points and the length of laminar wake are also very close. he natural vortex shedding past a cylinder is reproduced in the numerical solution. he successive pressure variation behind the cylinder is plotted in Fig. 6 for Reynolds number 280. he solution becomes un-symmetric after 25 cycles of iterations and becomes fully periodic vortex shedding after 56 cycles of iterations. he solution is marched by local time-stepping at the beginning till the flow field becomes un-symmetric. hen a constant time step is used to march the periodic solution. he predicted error in Strouhal number is 0% by using steps in a cycle and reduces to 7% by using steps in a cycle. For Reynolds number 000, the vortex street is fully turbulent. Figure 7 plots the horizontal and vertical velocity cyclic variation at the monitor point. he numerical solution is well periodic as shown by this plot. Figure 8 compares the instantaneous streamlines in nine instants over a cycle between the current solution and Arnone s results for laminar calculation. he agreement of the shedding of the vortex between these two solutions is very evident. he mechanism of vortex formation and merging is clearly shown in these plots. here are two different grid points used to study the grid resolution on the solution. he Reynolds number 000 was used. he calculation of 20x8 meshes requires a much more iterations than 20x2 meshes before the solution become un-symmetric as shown in Fig. 9. he instantaneous entropy contours between these two meshes are shown in Fig. 0. he finer the mesh, the better resolution of the wake is obtained. he predicted Strouhal number is quite similar in both cases. he evolution of pressure using four different time steps is shown in Fig.. he computed Strouhal number is normalized by the experimental value of Using 3500 iterations ( Δ t = 2 x0 ) in a cycle, the accuracy of Strouhal number is only 8%. 7 A time step of Δ t = x0 improves the accuracy of Strouhal number to 99%. he boundary layer is laminar over the front, separate, and breaks up into a turbulent wake above Reynolds number 300. A Reynolds number of 2000 was used for turbulent calculation in this study. he separation points moving forward as the Reynolds number reaches 2000, the boundary layer has now attained their upstream limit ahead of maximum thickness. he instantaneous entropy contour is presented for different turbulence models in Fig. 2. he calculation using κ-ωh model gives a much stronger distribution of eddy viscosity. he spreading of the wake structure is smallest among these four cases. he distance between the downstream vortices by using the κ-ωh model is also largest among these models. he low Reynolds number κ-ωl model gives a slighter less eddy viscosity. Hence, the spreading of the wake becomes wider. he JL and LS models give an even larger spreading angle. Also, the speed of vortex sheet is large so the trailing vortices are narrowly packed. he evolution of pressure using the four different turbulence models is shown in Fig. 3. he predicted Strouhal number is better by using the high Reynolds number κ-ω model. 5. Conclusions A Rung-Kutta integration method coupled with wall-distance-free turbulent models was used to study the periodic vortex shedding past a cylinder. he solution was validated by the streamlines for Reynolds number.54 and Numerical comparison between the current results and Arnone s results at Reynolds 000 was made for in nine instants over a cycle. he finer the mesh, the more iteration is required to form vortex sheet. he accuracy of Strouhal number requires a much smaller time step in the calculation. he use of different turbulence model in the codes has been explored for unsteady rowing up of vortex sheet. Only the wall-distance-free turbulence models were used to calculate the rowing up of vortex sheet. heκ-ωh gives a much larger eddy viscosity in the result. he distance between the vortices is increased. he resolution of vortex formation and rowing up is very similar to laminar solution by using JL turbulence model. he prediction of the Strouhal number is better by using the High Reynolds κ-ω turbulence model. More numerical test cases and maybe more advanced wall-distance-free turbulence models should be included to gain confidence in its predictive capability. 6. Acknowledgment he authors would like to thank the Chung Hua University (Contract No. CHU-94-M-05) for financially supporting this research. 7. References [] Lienhard, J. H., Synopsis of Lift, Drag and Vortex Frequency Data for Rigid Circular Cylinders, Washington State University, College of Engineering, Research Division Bulletin, pp. 300,
4 [2] Kimura,., and sutahara,., Fluid Dynamic Effects of Grooves on Circular Cylinder Surface, AIAA Journal, Vol. 29, No. 2, pp , 99. [3] Bearman, P. W., and Harvey, J. K., Control of Circular Cylinder Flow by the Use of Dimples, AIAA Journal, Vol. 3, No. 0, pp , 993. [4] Jameson, A., ime Dependent Calculations Using Multigrid with Applications to Unsteady Flow Past Airfoils and Wings, AIAA-9-596, 99. [5] Arnone, A., Liou, M. S. and Povinelli, L. A., Integration of Navier-Stokes Equations Using Dual time Stepping and a Multigrid Method, AIAA Journal, Vol. 33, No. 6, pp , 995. [6] Rumsey, C. L., Sanetrik, M. D., Biedron, R.., Melson, N. D., and Parlette, E. B., Efficiency and Accuracy of ime-accurate urbulent Navier-Stokes Computations, Computers & Fluids, Vol. 25, No. 2, pp , 996. [7] Jones, W. P., and Launder, B. E., he Prediction of Laminarization with a wo-equation Model of urbulence, International Journal of Heat and Mass ransfer, Vol. 5, pp , Feb [8] Launder, B. E., and Sharma, B. I., Application of the Energy Dissipation Model of urbulence to the Calculation of Flow Near a Spinning Disc, Letters in Heat and Mass ransfer, Vol., No. 2, pp. 3-38, 974. [9] Wilcox, D. C., Comparison of wo-equation urbulence Models for Boundary Layers with Pressure Gradient, AIAA Journal, Vol. 3, No. 8, pp , 993. [0] Wang, G. L., Numerical Calculation of Natural Unsteadiness in Laminar and urbulent Flows, M.S. Dissertation, Department of Mechanical Engineering, Chung Hua University, Hsin Chu, aiwan, 2006 [] Van Dyke, M., An Album of Fluid Motion, Department of Mechanical Engineering, Stanford University, he Parabolic Press, Stanford, California, Fourth printing, pp , 988. u/uo (40 X 20) (40 X 40) (40 X 80) analytical Fig.(a) Solution of Stokes second problem at t/=0 u/uo (20 X 40) (40 X 40) (80 X 40) 0.3 analytical η Fig.(b) Solution of Stokes second problem at t/=0.25 η 4
5 Fig. 3 he streamlines comparison at Re=.54. P P P, P2 Fig. 2(a) he computational grid for the cylinder Fig. 4 he evolution of pressure at the two monitor points (Re=28.4) Fig. 2(b) he monitor points behind the cylinder Fig. 5 he streamlines comparison at Re=28.4 5
6 .08 P P time = 4 P, P Fig. 6 he evolution of pressure at the two monitor points (Re=28.4) time = 5 time = Y-COMPONEN OF VELOCIY time = 7 time = X-COMPONENOFVELOCIY Fig. 7 Instantaneous u and v velocity variation time = 9 Fig. 8 Comparison of evolution of streamlines between the current solution (left) and Arnone s solution (right).06 ( 20 x 8 ) ( 20 x 2 ).04 time =.02 P 0.9 time = Fig. 9 he evolution of pressure at the monitor point using different meshes (Re=000) time = 3 6
7 20x2 κ-ωh 20x8 κ-ωl Fig. 0 Entropy plots by two different meshes on the wake structure (Re=000) D = 2.D-6 D =.D-6 D = 5.D-7 D =.D-7 P Fig. he instantaneous pressure variation at the monitor point using four different time steps (Re=000) LS 7
8 JL Fig. 2 Instantaneous Entropy plots by different turbulence models (Re=2000) Low Reynolds Wilcox High Reynolds Wilcox Jones-Launder Launder-Sharma Laminar P Fig. 3 he instantaneous pressure variation at the monitor point using four different time steps (Re=2000) able Maximum error for Stokes second problem at t/=0 ime Space % 0.04% 0.097% % 0.050% 0.049% % 0.023% 0.02% able2 Maximum error for Stokes second problem at t/= ime Space % 0.056% 0.053% % 0.03% 0.029% % 0.09% 0.08% 8
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University
More informationNumerical Simulation of a Blunt Airfoil Wake
6th Australasian Fluid Mechanics Conference Crown Plaza, Gold Coast, Australia 2-7 December 7 Numerical Simulation of a Blunt Airfoil Wake C.J. Doolan School of Mechanical Engineering University of Adelaide,
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 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 informationActive Control of Separated Cascade Flow
Chapter 5 Active Control of Separated Cascade Flow In this chapter, the possibility of active control using a synthetic jet applied to an unconventional axial stator-rotor arrangement is investigated.
More information1. Fluid Dynamics Around Airfoils
1. Fluid Dynamics Around Airfoils Two-dimensional flow around a streamlined shape Foces on an airfoil Distribution of pressue coefficient over an airfoil The variation of the lift coefficient with the
More informationNumerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder
Numerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder Ali Kianifar, Edris Yousefi Rad Abstract In many applications the flow that past bluff bodies have frequency nature (oscillated)
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationAerodynamic force analysis in high Reynolds number flows by Lamb vector integration
Aerodynamic force analysis in high Reynolds number flows by Lamb vector integration Claudio Marongiu, Renato Tognaccini 2 CIRA, Italian Center for Aerospace Research, Capua (CE), Italy E-mail: c.marongiu@cira.it
More informationNonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations p.1/50
Nonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations Matthew McMullen Advisor: Antony Jameson Co-advisor: Juan Alonso Sponsor: Accelerated Strategic Computing Initiative
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
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 informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Contours of Pressure.8.6.4...4.6.8 Matthew McMullen Antony Jameson Juan J. Alonso Aero Astro Dept. Stanford
More informationON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
Conference Applications of Mathematics 212 in honor of the 6th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 212 ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
More informationTURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS
HEFAT2014 10 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 2014 Orlando, Florida TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS Everts, M.,
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 informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More informationQ. Xiao 1 and H. M. Tsai 2 Temasek Laboratories National University of Singapore Singapore
AIAA Paper 24-56 42nd AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada 5-8 Jan 24 A NUMERICAL SUDY OF RANSONIC BUFFE ON A SUPERCRIICAL AIRFOIL Q. Xiao and H. M. sai 2 emasek Laboratories National
More informationComparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2
Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2 Department of Mechanical and Materials Engineering, Military Engineering Institute,
More informationSuppression of the unsteady vortex wakes of a circular cylinder pair by a doublet-like counter-rotation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids (2009) Published online in Wiley InterScience (www.interscience.wiley.com)..2075 Suppression of the unsteady vortex wakes
More informationMath 575-Lecture Failure of ideal fluid; Vanishing viscosity. 1.1 Drawbacks of ideal fluids. 1.2 vanishing viscosity
Math 575-Lecture 12 In this lecture, we investigate why the ideal fluid is not suitable sometimes; try to explain why the negative circulation appears in the airfoil and introduce the vortical wake to
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More information[N175] Development of Combined CAA-CFD Algorithm for the Efficient Simulation of Aerodynamic Noise Generation and Propagation
The 32nd International Congress and Exposition on Noise Control Engineering Jeju International Convention Center, Seogwipo, Korea, August 25-28, 2003 [N175] Development of Combined CAA-CFD Algorithm for
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 informationChapter 9 Flow over Immersed Bodies
57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 009 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,
More informationComputation of NACA0012 Airfoil Transonic Buffet Phenomenon with Unsteady Navier-Stokes Equations
5th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 9-2 January 22, Nashville, Tennessee AIAA 22-699 Computation of NACA2 Airfoil Transonic Buffet Phenomenon with
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 informationResolving the dependence on free-stream values for the k-omega turbulence model
Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok This report is based on
More informationVORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1395-1399 1395 VORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS by Li-Zhong HUANG a* and De-Ming NIE b a State Key Laboratory
More informationA fundamental study of the flow past a circular cylinder using Abaqus/CFD
A fundamental study of the flow past a circular cylinder using Abaqus/CFD Masami Sato, and Takaya Kobayashi Mechanical Design & Analysis Corporation Abstract: The latest release of Abaqus version 6.10
More informationFourth-Order Implicit Runge-Kutta Time Marching Using A Newton-Krylov Algorithm. Sammy Isono
Fourth-Order Implicit Runge-Kutta Time Marching Using A Newton-Krylov Algorithm by Sammy Isono A thesis submitted in conformity with the requirements for the degree of Master of Applied Science Graduate
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 informationNumerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation
, pp.49-58 http://dx.doi.org/10.1457/ijast.016.9.06 Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation Mojtaba Daneshi Department of Mechanical Engineering,
More informationNUMERICAL INVESTIGATION OF THE FLOW OVER A GOLF BALL IN THE SUBCRITICAL AND SUPERCRITICAL REGIMES
NUMERICAL INVESTIGATION OF THE FLOW OVER A GOLF BALL IN THE SUBCRITICAL AND SUPERCRITICAL REGIMES Clinton Smith 1, Nikolaos Beratlis 2, Elias Balaras 2, Kyle Squires 1, and Masaya Tsunoda 3 ABSTRACT Direct
More informationInvestigation of Co-Flow Jet Airfoil Mixing Mechanism Using Large Eddy Simulation
4st AIAA Fluid Dynamics Conference and Exhibit 27-3 June 2, Honolulu, Hawaii AIAA 2-398 Investigation of Co-Flow Jet Airfoil Mixing Mechanism Using Large Eddy Simulation Hong-Sik IM, Ge-Cheng Zha, Bertrand
More informationNumerical study of the effects of trailing-edge bluntness on highly turbulent hydro-foil flows
Numerical study of the effects of trailing-edge bluntness on highly turbulent hydro-foil flows T. Do L. Chen J. Tu B. Anderson 7 November 2005 Abstract Flow-induced noise from fully submerged lifting bodies
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationρ t + (ρu j ) = 0 (2.1) x j +U j = 0 (2.3) ρ +ρ U j ρ
Chapter 2 Mathematical Models The following sections present the equations which are used in the numerical simulations documented in this thesis. For clarity, equations have been presented in Cartesian
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationNUMERICAL SIMULATION AND MODELING OF UNSTEADY FLOW AROUND AN AIRFOIL. (AERODYNAMIC FORM)
Journal of Fundamental and Applied Sciences ISSN 1112-9867 Available online at http://www.jfas.info NUMERICAL SIMULATION AND MODELING OF UNSTEADY FLOW AROUND AN AIRFOIL. (AERODYNAMIC FORM) M. Y. Habib
More informationSIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS
ICAS 2000 CONGRESS SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS H Yao, R K Cooper, and S Raghunathan School of Aeronautical Engineering The Queen s University of Belfast, Belfast BT7 1NN,
More informationNumerical Study of Transonic Buffet on a Supercritical Airfoil
AIAA JOURNAL Vol. 44, No. 3, March 2006 Numerical Study of Transonic Buffet on a Supercritical Airfoil Q. Xiao and H. M. Tsai National University of Singapore, Singapore 119260, Republic of Singapore and
More informationVortex Induced Vibrations
Vortex Induced Vibrations By: Abhiroop Jayanthi Indian Institute of Technology, Delhi Some Questions! What is VIV? What are the details of a steady approach flow past a stationary cylinder? How and why
More informationAalto University School of Science and Technology CFD-group/ Department of Applied Mechanics. MEMO No CFD/MECHA DATE: March 15, 2012
Aalto University School of Science and Technology CFD-group/ Department of Applied Mechanics MEMO No CFD/MECHA-2-212 DATE: March 15, 212 TITLE The Effect of Free-Stream Turbulence Parameters on the SST
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
More informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationDirect Numerical Simulations of Plunging Airfoils
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 010, Orlando, Florida AIAA 010-78 Direct Numerical Simulations of Plunging Airfoils Yves Allaneau
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 informationNUMERICAL SIMULATION OF UNSTEADY CAVITATING FLOWS
NUMERICAL SIMULATION OF UNSTEADY CAVITATING FLOWS Charles C.S. Song and Qiao Qin St. Anthony Falls Laboratory, University of Minnesota Mississippi River at 3 rd Ave. SE, Minneapolis, MN 55414, USA ABSTRACT
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 informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationNumerical Simulation of Flow Around An Elliptical Cylinder at High Reynolds Numbers
International Journal of Fluids Engineering. ISSN 0974-3138 Volume 5, Number 1 (2013), pp. 29-37 International Research Publication House http://www.irphouse.com Numerical Simulation of Flow Around An
More informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
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 informationModule 2: External Flows Lecture 12: Flow Over Curved Surfaces. The Lecture Contains: Description of Flow past a Circular Cylinder
The Lecture Contains: Description of Flow past a Circular Cylinder Experimental Results for Circular Cylinder Flow file:///d /Web%20Course%20(Ganesh%20Rana)/Dr.%20gautam%20biswas/Final/convective_heat_and_mass_transfer/lecture12/12_1.htm[12/24/2014
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
More informationApplication of a Virtual-Boundary Method for the Numerical Study of Oscillations Developing Behind a Cylinder Near A Plane Wall
Fluid Dynamics, Vol. 39, No. 1, 2004, pp. 61 68. Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2004, pp. 69 77. Original Russian Text Copyright 2004 by Kit, Nikitin,
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 informationNumerical Investigation of Vortex Induced Vibration of Two Cylinders in Side by Side Arrangement
Numerical Investigation of Vortex Induced Vibration of Two Cylinders in Side by Side Arrangement Sourav Kumar Kar a, 1,, Harshit Mishra a, 2, Rishitosh Ranjan b, 3 Undergraduate Student a, Assitant Proffessor
More informationOpenFOAM selected solver
OpenFOAM selected solver Roberto Pieri - SCS Italy 16-18 June 2014 Introduction to Navier-Stokes equations and RANS Turbulence modelling Numeric discretization Navier-Stokes equations Convective term {}}{
More informationImproved numerical simulation of bridge deck aeroelasticity by model validation
Improved numerical simulation of bridge deck aeroelasticity by model validation A.Šarkić, R. Höffer Building Aerodynamics Laboratory, Bochum, Germany anina.sarkic@rub.de Abstract In this study, the results
More information41st AIAA Aerospace Sciences Meeting and Exhibit January 6 9, 2003/Reno, NV
AIAA 3 267 Calculation of Unsteady Transonic Flow by an Euler Method with Small Perturbation Boundary Conditions C. Gao, S. Luo, F. Liu Department of Mechanical and Aerospace Engineering University of
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More information2-D and 3-D Assessment of Cambered and Symmetric Airfoils: A CFD Study
Clemson University TigerPrints All Theses Theses 12-2009 2-D and 3-D Assessment of Cambered and Symmetric Airfoils: A CFD Study Tuncay Kamas Clemson University, tkamas@clemson.edu Follow this and additional
More informationMasters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16
Masters in Mechanical Engineering Aerodynamics st Semester 05/6 Exam st season, 8 January 06 Name : Time : 8:30 Number: Duration : 3 hours st Part : No textbooks/notes allowed nd Part : Textbooks allowed
More informationMestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13
Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 212/13 Exam 2ª época, 2 February 213 Name : Time : 8: Number: Duration : 3 hours 1 st Part : No textbooks/notes allowed 2 nd Part :
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 informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More informationFLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationHEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE
HEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE A.K. Pozarlik 1, D. Panara, J.B.W. Kok 1, T.H. van der Meer 1 1 Laboratory of Thermal Engineering,
More informationAvailable online at ScienceDirect. Procedia Engineering 79 (2014 ) 49 54
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 79 (2014 ) 49 54 37th National Conference on Theoretical and Applied Mechanics (37th NCTAM 2013) & The 1st International Conference
More informationChapter 9 Flow over Immersed Bodies
57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 014 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,
More informationAerodynamics. High-Lift Devices
High-Lift Devices Devices to increase the lift coefficient by geometry changes (camber and/or chord) and/or boundary-layer control (avoid flow separation - Flaps, trailing edge devices - Slats, leading
More informationSelf-Excited Vibration in Hydraulic Ball Check Valve
Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationUNIT IV BOUNDARY LAYER AND FLOW THROUGH PIPES Definition of boundary layer Thickness and classification Displacement and momentum thickness Development of laminar and turbulent flows in circular pipes
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationSENSITIVITY ANALYSIS OF THE FACTORS AFFECTING FORCE GENERATION BY WING FLAPPING MOTION
Proceedings of the ASME 2013 International Mechanical Engineering Congress and Exposition IMECE2013 November 15-21, 2013, San Diego, California, USA IMECE2013-65472 SENSITIVITY ANALYSIS OF THE FACTORS
More informationPeriodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain.
Flow Characterization of Inclined Jet in Cross Flow for Thin Film Cooling via Large Eddy Simulation Naqavi, I.Z. 1, Savory, E. 2 and Martinuzzi, R. J. 3 1,2 The Univ. of Western Ontario, Dept. of Mech.
More informationModule 3: Velocity Measurement Lecture 15: Processing velocity vectors. The Lecture Contains: Data Analysis from Velocity Vectors
The Lecture Contains: Data Analysis from Velocity Vectors Velocity Differentials Vorticity and Circulation RMS Velocity Drag Coefficient Streamlines Turbulent Kinetic Energy Budget file:///g /optical_measurement/lecture15/15_1.htm[5/7/2012
More informationPublication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada
ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature
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 informationAerodynamic Rotor Model for Unsteady Flow and Wake Impact
Aerodynamic Rotor Model for Unsteady Flow and Wake Impact N. Bampalas, J. M. R. Graham Department of Aeronautics, Imperial College London, Prince Consort Road, London, SW7 2AZ June 28 1 (Steady Kutta condition)
More informationAPPLICATION OF SPACE-TIME MAPPING ANALYSIS METHOD TO UNSTEADY NONLINEAR GUST-AIRFOIL INTERACTION PROBLEM
AIAA 2003-3693 APPLICATION OF SPACE-TIME MAPPING ANALYSIS METHOD TO UNSTEADY NONLINEAR GUST-AIRFOIL INTERACTION PROBLEM Vladimir V. Golubev* and Axel Rohde Embry-Riddle Aeronautical University Daytona
More informationTURBULENCE MODELING VALIDATION FOR AFTERBODY FLOWS A. Hadjadj 1 and A.N. Kudryavtsev 2
TURBULENCE MODELING VALIDATION FOR AFTERBODY FLOWS A. Hadjadj 1 and A.N. Kudryavtsev 2 1 LMFN-CORIA, UMR 6614, INSA de Rouen, B.P. 8, 76801 Saint-Etienne du Rouvray, France 2 Institute of Theoretical and
More informationLecture-4. Flow Past Immersed Bodies
Lecture-4 Flow Past Immersed Bodies Learning objectives After completing this lecture, you should be able to: Identify and discuss the features of external flow Explain the fundamental characteristics
More informationThe current issue and full text archive of this journal is available at
The current issue and full text archive of this journal is available at www.emeraldinsight.com/0961-5539.htm HFF 16,6 660 Received February 2005 Revised December 2005 Accepted December 2005 3D unsteady
More informationINTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 1, No 3, 2010
CFD analysis of 2D unsteady flow around a square cylinder Gera.B, Pavan K. Sharma, Singh R.K Reactor Safety Division, Bhabha Atomic Research Centre, Trombay, Mumbai, India 400085 pa1.sharma@gmail.com ABSTRACT
More informationTurbulence modelling. Sørensen, Niels N. Publication date: Link back to DTU Orbit
Downloaded from orbit.dtu.dk on: Dec 19, 2017 Turbulence modelling Sørensen, Niels N. Publication date: 2010 Link back to DTU Orbit Citation (APA): Sørensen, N. N. (2010). Turbulence modelling. Paper presented
More informationcompression corner flows with high deflection angle, for example, the method cannot predict the location
4nd AIAA Aerospace Sciences Meeting and Exhibit 5-8 January 4, Reno, Nevada Modeling the effect of shock unsteadiness in shock-wave/ turbulent boundary layer interactions AIAA 4-9 Krishnendu Sinha*, Krishnan
More informationNon-Synchronous Vibrations of Turbomachinery Airfoils
Non-Synchronous Vibrations of Turbomachinery Airfoils 600 500 NSV Frequency,!, hz 400 300 200 F.R. Flutter 100 SFV 0 0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 Rotor Speed,!, RPM Kenneth C. Hall,
More informationStabilization of Explicit Flow Solvers Using a Proper Orthogonal Decomposition Technique
50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 09-12 January 2012, Nashville, Tennessee AIAA 2012-1096 Stabilization of Explicit Flow Solvers Using a Proper
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 informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationSIMULATION OF UNSTEADY TURBULENT FLOW OVER A STALLED AIRFOIL
Computational Fluid Dynamics JOURNAL 14(4):42 January 2006 (pp.359 378) SIMULATION OF UNSTEADY TURBULENT FLOW OVER A STALLED AIRFOIL B. STOCKDILL G. PEDRO A. SULEMAN F. MAGAGNATO N. DJILALI Abstract Unsteady,
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationCalculation of the Flow around the KVLCC2M Tanker
Calculation of the Flow around the KVLCC2M Tanker L. Eça 1, M. Hoekstra 2 and S.L. Toxopeus 2 1 Instituto Superior Técnico, Portugal 2 Maritime Research Institute, Netherlands SUMMARY The flow around the
More information