On the Turbulence Modelling for an Air Cavity Interface
|
|
- Irma Pitts
- 5 years ago
- Views:
Transcription
1 1 Introduction On the Turbulence Modelling for an Air Cavity Interface Gem Rotte, Maarten Kerkvliet, and Tom van Terwisga TU Delft, Delft/The Netherlands MARIN, Wageningen/The Netherlands The use of air lubrication techniques can significantly reduce a ship s fuel consumption. One of the most promising techniques applicable to ships is the external air cavity technique. An external cavity is created beneath the ship s hull with the help of a cavitator, which is located directly upstream of an air injection point (figure 1a). The cavitator is extended in the span-wise direction and typically has a rectangular or triangular cross section. It is used to separate the mean water flow from the hull, thereby providing a stable air layer. Air cavity applications are claimed to lead to propulsive power reduction percentages of 10-20% by reducing the ship s frictional drag (Gorbachev and Amromin (2014), Zverkhovskyi (2014)). However, a complete understanding of the influence that the ship s hull form has on the relevant twophase flow physics and thereby also on the length and stability of the air cavities is still lacking. This inability to predict the air cavity characteristics hampers the application of air cavity techniques in the shipping industry. Multiphase CFD methods can help us to gain a better understanding of the relevant physics. Fig. 1: The external air cavity technique (a) and two possible mechanisms for air discharge: re-entrant jet (b) and wave pinch-off (c) The largest challenge in predicting the air cavity characteristics lies in correctly modelling their closure region (Zverkhovskyi et al. (2015), Shiri et al. (2012)). Here, the closure region is defined as the region where the separated air-water flow transforms into a more dispersed flow. Both the re-entrant jet mechanism and the wave pinch-off mechanism are cited as possible mechanisms for air discharge from the cavity in the closure region (e.g. Shiri et al. (2012), Zverkhovskyi (2014), Mäkiharju (2012)). The reentrant jet mechanism is provisionally assumed to be similar to the re-entrant jet mechanism responsible for the break-up of natural sheet cavities. Callenaere et al. (2001) describe the mechanism as follows. In the closure region of a sheet cavity, a region with high adverse pressure gradient is formed. This increase in local pressure forces a thin stream of liquid into the cavity. This thin stream is called the re-entrant jet and it is illustrated in figure 1b. Impingement of this jet with the gas-liquid interface results in a disturbance leading to a portion of the attached cavity being pinched off and transported by the mean flow. The wave pinch-off mechanism, as illustrated in figure 1c, is governed by waves on the air-water interface. These waves are believed to be formed by turbulence structures interfering with the interface. If the resulting wave amplitudes are of the same magnitude as the cavity thickness, the free surface interface will hit the ship s bottom and a pocket of air will be shed from the cavity. This article aims to link the physical modelling of the relevant phenomena to their numerical modelling, with an emphasis on the modelling of the re-entrant jet mechanism with a RaNS turbulence model. The article is based on the available literature in the public domain and on knowledge gained from research projects carried out at Delft University of Technology and at the Maritime Research Institute Netherlands (MARIN).
2 The numerical solver used for all simulations is ReFRESCO (MARIN (2017)). It is a viscous-flow CFD code that solves multiphase (unsteady) incompressible flows using the RaNS equations. It is complemented with turbulence models, cavitation models and volume-fraction transport equations for different phases. The equations are discretised using a finite-volume approach. Time integration is performed implicitly with first or second-order backward schemes. The implementation is face-based, which facilitates grids with elements consisting of an arbitrary number of faces (hexahedrals, tetrahedrals, prisms, pyramids, etc.) and, if needed, h-refinement (hanging nodes). For turbulence modelling, RaNS/URaNS, SAS (Scale Adaptive Simulation), DES (Detached Eddy Simulation), PANS (Partially Averaged Navier-Stokes) and LES (Large Eddy Simulation) approaches can be used. 2 Base case calculations The (2D) base case calculations are based on the measurements performed in the Delft Cavitation Tunnel by Zverkhovskyi (2014). The cavitation tunnel has a test section length of 2m and a m crosssectional area. The length of the computational domain upstream of the cavity is chosen such that the boundary layer profile in de simulations approaches the one from the measurements as close as possible. A schematic of the domain including the applied boundary conditions is shown in figure 2. The close-up shows the grid topology in the region including the cavity, with refinement zones in the region around the cavitator, the cavity closure and the expected cavity interface. Close to the wall grid clustering was applied to ensure y+ values well below 1. Unsteady simulations were carried out with a water tunnel inlet velocity of 1 m/s. Based on current experimental observations in the Delft Cavitation tunnel, this is the free-stream velocity at which the re-entrant jet mechanism is dominant over the wave pinch-off mechanism. All computations were carried out using the k kl (KSKL) turbulence model as described by Menter et al. (2006). A second order upwind scheme is used for the time discretisation and the time step used in the calculations was 0.001s to ensure Courant numbers in the order of 1. The discritisation of both the momentum and turbulence equations was done using the QUICK scheme and the volume fraction transport equation was discretised using a Fromm scheme with a SUPERBEE limiter. Fig. 2: Computational domain and boundary conditions Zverkhovskyi et al. (2015) found that prescribing the air flux at the air inflow did not yet give physical solutions. Figure 3 shows the instantaneous cavity profiles for simulations with an air inlet pressure of 40, 45, 50, 55 and 60 Pa respectively. The profiles are plotted using the contour where the air volume fraction equals 0.5. Based on linear theory and the experiments from Zverkhovskyi (2014) a cavity length 0.32 m and a cavity thickness of around 0.01 m is expected. The pressure at the air inlet for the base case calculation is chosen such that the the cavity profile from the simulations matches the one from the experiment as close as possible, which resulted in an air inlet pressure of 50 Pa. Fig. 3: Cavity profile sensitivity to pressure at the air inflow boundary
3 The base case simulation did not show correct cavity closure behaviour. No break-up of the air cavity and thus no shedding was observed. A thin layer of air is constantly present at the closure and downstream of the cavity, influencing the computed wall shear stress. Figure 4 shows that the overall computed velocity profiles compare reasonably well with the measured profiles. However, the computed upstream profile already differs from the experimental one. This could be the reason for the differences in the middle and downstream of the cavity. The simulated profile is somewhat thicker and seems more turbulent. As for the profile at the middle of the cavity, the simulations predict higher velocities outside the cavity. This seems a reasonable solution due to the narrowing of the tunnel as a result of the presence of the cavity. Unfortunately the results of the experimental measurements close to the wall are of insufficient resolution to quantitatively compare the viscous sub-layer of the boundary layer, which determines the wall shear stress. Fig. 4: Time-averaged boundary layer profiles at three different streamwise locations: before in the middle and after the cavity. The 2-phase base case simulations from the current project and the experiments from Zverkhovskyi (2014) are shown 3 Eddy-viscosity correction For the numerical modelling of the re-entrant jet in the air cavity cases, it was investigated to see whether RaNS methods can be used, since these are the most interesting option for application in practical ship design. For the simulation of unsteady natural sheet cavity dynamics, and in particular the re-entrant jet mechanism occurring in those cases, RaNS methods can already be applied. However, the commonly used two-equation turbulence models, such as the k-ɛ and k-ω models, were originally developed for RaNS simulations of single-phase flows. A frequently encountered problem with these turbulence models that are based on the Boussinesq assumption is the over-prediction of the eddy-viscosity around the liquid-gas interface and thus in and downstream of the cavity closure region. According to Reboud et al. (1998), this over-predicted eddy-viscosity reduces the development of the re-entrant jet and can thus prevent the occurrence of shedding. The above mentioned tendency is not just associated with one specific RaNS solver or cavitation model. There are several additional correction models that artificially lower the over-predicted eddy-viscosity, of which the Reboud-correction is just one. Here, the formulation for the turbulence viscosity is modified. The originally computed turbulence viscosity µ t is multiplied by a correction factor based on a function of the local density f (ρ). Computations were carried out using two different correction functions: the function as described in Coutier-Delgosha et al. (2001) and Reboud et al. (1998) and a correction based on a Gaussian function which is proposed here. The same numerical settings as given in section 2 are applied. The function as described by Coutier-Delgosha et al. (2001) is given in Eq. (1). Here, ρ g is the density of the gas, ρ l is the density of the liquid, ρ is the local cell density and n influences the amount of reduction. For natural cavitation problems, n usually has a value of 10. The eddy viscosity correction factor is obtained by dividing f (ρ) by the local cell density, resulting in a value between 0 and 1. ( ) ρg ρ n f (ρ) = ρ g + (ρ l ρ g ) (1) ρ g ρ l
4 The correction based on a Gaussian curve is proposed in Eq. (2). Here, A and W are the amplitude and width of the peak and have a value between 0 and 1, and p controls the order and thus the width of the peak. Figure 5 shows the different correction functions for different values of n for cases of Reboud correction and different values of A, W and p for Gaussian correction cases. It can be observed that the Reboud correction function is asymmetrical and has a bias towards the lower density fluid. It also has large gradients when the mixture approaches pure air or pure water, which can possibly introduce difficulties when numerically solving the problem. For the proposed Gaussian functions the correction only is active when the cell is marked as an interface cell and the function is symmetric. ( ρ 1 2 (3ρ l ρg) ) 2 g(ρ) = Ae 2(W(ρ l ρg)/10) p (2) Fig. 5: Eddy-viscosity correction factor as a function of the local density ρ. The computed eddy viscosity is multiplied by this correction factor. Figure 6 shows the computed cavity profiles for three representative cases: the base case, one case with Reboud correction and one case with Gaussian correction. Most probably, the asymmetric nature of the Reboud correction function is the reason that the cavity profiles as computed including this function are very unstable. This instability is present for all n > 1, thus for all cases were the correction function is applied. For all simulations with n > 1, the cavity grows to its stable length which is similar to to profile of the base case simulation and then starts to get unstable. The air volume fraction inside the cavity is equal to or larger than 0.99, but due to the bias of the Reboud correction function towards the low density fluid, all eddy-viscosity ν t inside the cavity is reduced such that there is no eddy-viscosity at all. This results in large difference between the viscosity inside and outside the cavity, causing instabilities on the air-water interface. The total viscosity difference between inside and outside the cavity is in the order of ten times lower for the base case when compared to the the Reboud cases. The total eddy-viscosity outside the cavity is rather high, due to the developed boundary layer upstream of the cavity (Re x > 10 6 ). An exploratory sensitivity study confirmed that the onset of these instabilities is also the case for smaller time steps and with finer grids around the cavity. In contrast to the unstable behaviour which is found for the simulations including the Reboud correction function, all three simulations as carried out using the Gaussian correction produced very stable cavity profiles (figure 6). All three correction functions with an amplitude of 1 which are shown in figure 5 were applied and produced more or less the same result as found for the base case calculation, showing incorrect closure behaviour. This lead to the conclusion that the reduction of the eddy-viscosity on the air-water interface is probably not the only requirement for the re-entrant jet to form and the production of eddy viscosity at the cavity interface as well as inside the cavity have a large influence on the interface stability. Figure 7 shows the regions where the correction function was active
5 for both the Reboud function and Gaussian function. The Reboud function is active in a large region inside the cavity as well. whereas the Gaussian correction function only acts on the cavity interface. For the boundary layer profiles, similar profiles were found for all correction cases. Only the simulations including the Reboud correction did show some irregular behaviour behind the cavity closure, most probably due to passing pockets of air-water mixture bubbles. Fig. 6: Cavity profiles for different eddy viscosity correction functions Fig. 7: Regions where the eddy viscosity correction function was active; 1: no reduction, 0: maximum reduction 4 Discussion Most of this article focuses on the numerical simulation of the re-entrant jet mechanism for ventilated air cavity flows with a RaNS model. Solely correcting the eddy-viscosity on the air-water interface appears to be insufficient to initiate the re-entrant jet shedding mechanism. However, the computed eddy-viscosity at the cavity interface as well as inside the cavity have a large influence on the interface stability. In contrast to the RaNS modelling of natural cavitation problems, there is a developed boundary layer upstream of the cavity in the ventilated cavity cases. The eddy-viscosity produced in this region is transported with the flow, causing a thick band of viscous flow around the cavity. This could also hamper the formation of the re-entrant jet in the simulations. This hypothesis is to be strengthened by a simulation with a very short no-slip wall upstream of the cavity, thus almost no turbulent boundary layer formation. Maybe a re-entrant jet can then exist (simulation is running now). Another difference with respect to natural sheet cavitation is the lack of a condensation term. There is no destruction of gas at the cavity closure. Since the re-entrant jet mechanism is a 3D phenomenon, simulations also have to be carried out in 3D, this will most probably be added to the final version of this article. In cavitation problems you need very high convergence levels and small time steps. An exploratory study of the influence of the time step size did not show any differences for smaller timesteps. Also, the RaNS modelling of natural cavitation problems usually ask for a low-order and thus diffuse discretisation scheme for the volume fraction transport equation. Here a higher order scheme is used, since it is hypothesised that sharp interfaces are required in order to compute the cavity surface correctly, especially when one is also interested in the wave pinch-off mechanism and the total air flux leaving the cavity.
6 5 Conclusions and further work To conclude, the simulated velocity profiles at different streamwise locations in the boundary layer around the cavity compare reasonably well with the experimental profiles. However, RaNS turbulence models show serious deficiencies for the modelling of the re-entrant jet mechanism for air cavities. Also not after including an ad-hoc eddy-viscosity correction model as proposed by Reboud and in this paper. Maybe a Reynolds Stress Model (RSM) could help, since the isotropic turbulence modelling as used in classical 2-eq. RaNS turbulence models is a violation of the anisotropic nature of the turbulence structures near the interface. The Reynolds stresses across the interface should be close to zero as no significant momentum flux should cross the interface. Next to the modelling of the re-entrant jet, RaNS methods are especially not suitable when one also wants to include the simulation of the wave pinch-off mechanism. This mechanism is hypothesised to be governed by waves formed by turbulence structures disturbing the interface. These structures are not solved using a RaNS method and one would need a scale resolving or hybrid model such LES, RaNS- LES or a PANS method. The use of DES requires caution, since the cavities are located inside a boundary layer. In these regions a regular DES is still in its RaNS-mode, but e.g. a iddes could be an option. From a practical point of view however, these scale resolving models are not believed to be viable for practical ship design purposes yet, due to the significant amount of added computational time. These models can however first be used to identify the physical conditions in the closure region and to study and identify scaling influences. Thereafter, one could develop a RaNS model for the numerical modelling of air cavitites. The dispersion (or diffusion) of the gas phase can then for example be modelled by an extra convection-diffusion term to break up the air layer. This extra term can then be activated when certain conditions characterised by the physics in the closure region are met. Further work will first focus on a performance study for different RaNS and hybrid RaNS-LES turbulence models around air-water interfaces using a benchmark case such as a Kelvin-Helmholtz stability analysis in the vicinity of a wall with a known disturbance. Thereafter hybrid RaNS-LES,LES or PANS simulations could be carried out to study the physical mechanisms responsible for air discharge in the closure region. Acknowledgements We thank the Netherlands Organisation for Scientific Research NWO, by whom this research is financially supported. References Callenaere, M., Franc, J.-P., Michel, J.-M. and Riondet, M. (2001). The cavitation instability induced by the development of a re-entrant jet. J. Fluid Mech 444, Coutier-Delgosha, O., Fortess-Patella, R. and Reboud, J.L. (2003). Evaluation of the turbulence model influence on the numerical simulations of unsteady cavitation. J. of Fluids Eng. 125(1), Gorbachev, Y. and Amromin, E. (2012). Ship Drag Reduction by Hull Ventilation from Laval to Near Future: Challenges and Successes. ATMA 2012 Mäkiharju, S.A. (2012). The Dynamics of Ventilated Partial Cavities over a Wide Range of Reynolds Numbers ans Quantitative 2D X-Ray Densitometry for Multiphase Flow. Ph.D. Thesis The University of Michigan ReFRESCO (2017) Online; accessed 23 March 2017 Menter, F.R., Egorov, Y. and Rusch, D. (2006). Steady and unsteady flow modelling using the k kl model. Ichmt Digital Library Online, Begel House Inc. Reboud, J., Stutz, B. and Coutier, O. (1998). Two-phase Flow Structure of Cavitation: Experiment and modelling of Unsteady Effects. Third International Symposium on Cavitation, Grenoble, France Shiri, A., Leer-Andersen, M.,Bensow, R.E. and Norrby, J. (2012). Hydrodynamics of a Displacement Air Cavity Ship. 29 th Symposium on Naval Hydrodynamics, Gothenburg, Sweden Zverkhovskyi, O. (2014). Ship Drag Reduction by Air Cavities. Ph.D. Thesis, TU Delft Zverkhovskyi, O., Kerkvliet, M., Vaz, G. and van Terwisga, T. (2015). Numerical Study on Air Cavity Flows. Proceedings of NuTTS 2015, Cortona, Italy
Turbulent 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 informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More informationHybrid RANS/LES simulations of a cavitating flow in Venturi
Hybrid RANS/LES simulations of a cavitating flow in Venturi TSFP-8, 28-31 August 2013 - Poitiers - France Jean Decaix jean.decaix@hevs.ch University of Applied Sciences, Western Switzerland Eric Goncalves
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationNumerical investigation of cavitation-regimes in a converging-diverging nozzle
Numerical investigation of cavitation-regimes in a converging-diverging nozzle 1 Polina Gorkh, 1 Steffen J. Schmidt, and 1 Nikolaus A. Adams 1 Institute of Aerodynamics and Fluid Mechanics, Department
More informationNUMERICAL MODELLING OF UNSTEADY PARTIAL CAVITIES BEHIND A
FOURTH INTERNATIONAL SYMPOSIUM ON CAVITATION, JUNE 20-23, 2001, PASADENA, CALIFORNIA, USA NUMERICAL MODELLING OF UNSTEADY PARTIAL CAVITIES BEHIND A BACKWARD FACING STEP THIERRY MAÎTRE, CHRISTIAN PELLONE
More informationComputation of Unsteady Flows With Moving Grids
Computation of Unsteady Flows With Moving Grids Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de Unsteady Flows With Moving Boundaries, I Unsteady flows
More informationNumerical investigation of swirl flow inside a supersonic nozzle
Advances in Fluid Mechanics IX 131 Numerical investigation of swirl flow inside a supersonic nozzle E. Eslamian, H. Shirvani & A. Shirvani Faculty of Science and Technology, Anglia Ruskin University, UK
More informationINTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING---
ILASS-Europe 2002 Zaragoza 9 11 September 2002 INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING--- Z. Tapia, A. Chávez e-mail: ztapia@imp.mx Instituto Mexicano del Petróleo Blvd. Adolfo Ruiz Cortines
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 informationThere are no simple turbulent flows
Turbulence 1 There are no simple turbulent flows Turbulent boundary layer: Instantaneous velocity field (snapshot) Ref: Prof. M. Gad-el-Hak, University of Notre Dame Prediction of turbulent flows standard
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 informationFlow Induced Noise Generation By Partial Cavities
Flow Induced Noise Generation By Partial Cavities Harish Ganesh, Juliana Wu and Steven L. Ceccio University of Michigan, Ann Arbor USA 2nd FLINOVIA Symposium Penn-State University, State College, PA Sponsor:
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 informationComparison of two equations closure turbulence models for the prediction of heat and mass transfer in a mechanically ventilated enclosure
Proceedings of 4 th ICCHMT May 17-0, 005, Paris-Cachan, FRANCE 381 Comparison of two equations closure turbulence models for the prediction of heat and mass transfer in a mechanically ventilated enclosure
More informationNUMERICAL SIMULATION OF CAVITATING FLOW AROUND A 2D HYDROFOIL: A BAROTROPIC APPROACH
Cav03-OS--005 Fifth International Symposium on Cavitation (CAV003) Osaka, Japan, November -4, 003 NUMERICAL SIMULATION OF CAVITATING FLOW AROUND A D HYDROFOIL: A BAROTROPIC APPROACH B. Pouffary Laboratoire
More informationComputations of unsteady cavitating flow on wing profiles using a volume fraction method and mass transfer models
Computations of unsteady cavitating flow on wing profiles using a volume fraction method and mass transfer models Tobias Huuva *, Aurelia Cure, Göran Bark * and Håkan Nilsson * Chalmers University of Technology,
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 informationZonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach
Zonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach S. Arvidson 1,2, L. Davidson 1, S.-H. Peng 1,3 1 Chalmers University of Technology 2 SAAB AB, Aeronautics 3 FOI, Swedish Defence
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More 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 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 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 informationINSTITUTTET FOR BYGNINGSTEKNIK DEPT. OF BUILDING TECHNOLOGY AND STRUCTURAL ENGINEERING AALBORG UNIVERSITET AAU AALBORG DANMARK
INSTITUTTET FOR BYGNINGSTEKNIK DEPT. OF BUILDING TECHNOLOGY AND STRUCTURAL ENGINEERING AALBORG UNIVERSITET AAU AALBORG DANMARK Lars Davidson and Peter V. Nielsen A Study of Laminar Backward-Facing Step
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 informationsmp 11 Workshop - Case 1: DelftFoil
Second International Symposium on Marine Propulsors smp 11, Hamburg, Germany, June 2011 Workshop: Cavitation on foil smp 11 Workshop - Case 1: DelftFoil M. Hoekstra 1, T. van Terwisga 1,2, E.J. Foeth 1
More informationSELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE
Computational Thermal Sciences, 3 (1): 63 72 (2011) SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE M. Zhao, 1, M. Yang, 1 M. Lu, 1 & Y. W. Zhang 2 1
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More 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 informationLarge eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions
June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering
More informationENGINEERING MECHANICS 2012 pp Svratka, Czech Republic, May 14 17, 2012 Paper #195
. 18 m 2012 th International Conference ENGINEERING MECHANICS 2012 pp. 309 315 Svratka, Czech Republic, May 14 17, 2012 Paper #195 NUMERICAL SIMULATION OF TRANSITIONAL FLOWS WITH LAMINAR KINETIC ENERGY
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 informationNumerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel
Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel Arunanshu Chakravarty 1* 1 CTU in Prague, Faculty of Mechanical Engineering, Department of Process Engineering,Technická
More informationCFD as a Tool for Thermal Comfort Assessment
CFD as a Tool for Thermal Comfort Assessment Dimitrios Koubogiannis dkoubog@teiath.gr G. Tsimperoudis, E. Karvelas Department of Energy Technology Engineering Technological Educational Institute of Athens
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 informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationTable of Contents. Foreword... xiii. Preface... xv
Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...
More 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 informationOn the transient modelling of impinging jets heat transfer. A practical approach
Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,
More informationCavitation Control on Hydrofoils
Proceedings of the International Conference on Heat Transfer and Fluid Flow Prague, Czech Republic, August 11-12, 2014 Paper No. 181 Cavitation Control on Hydrofoils Mohammad Mortezazadeh, Ali Katal, Khodayar
More informationNumerical Studies of Droplet Deformation and Break-up
ILASS Americas 14th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2001 Numerical Studies of Droplet Deformation and Break-up B. T. Helenbrook Department of Mechanical and
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 informationCFD SIMULATION OF A SINGLE PHASE FLOW IN A PIPE SEPARATOR USING REYNOLDS STRESS METHOD
CFD SIMULATION OF A SINGLE PHASE FLOW IN A PIPE SEPARATOR USING REYNOLDS STRESS METHOD Eyitayo A. Afolabi 1 and J. G. M. Lee 2 1 Department of Chemical Engineering, Federal University of Technology, Minna,
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 informationExperimental Verification of CFD Modeling of Turbulent Flow over Circular Cavities using FLUENT
Experimental Verification of CFD Modeling of Turbulent Flow over Circular Cavities using FLUENT T Hering, J Dybenko, E Savory Mech. & Material Engineering Dept., University of Western Ontario, London,
More 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 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 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 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 informationLES of synthetic jets in boundary layer with laminar separation caused by adverse pressure gradient
LES of synthetic jets in boundary layer with laminar separation caused by adverse pressure gradient Tetsuya Ozawa, Samuel Lesbros and Guang Hong* University of Technology, Sydney (UTS), Australia *Corresponding
More informationINVESTIGATION OF THE FLOW OVER AN OSCILLATING CYLINDER WITH THE VERY LARGE EDDY SIMULATION MODEL
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
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 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 informationChapter 5 Control Volume Approach and Continuity Equation
Chapter 5 Control Volume Approach and Continuity Equation Lagrangian and Eulerian Approach To evaluate the pressure and velocities at arbitrary locations in a flow field. The flow into a sudden contraction,
More informationBOUNDARY LAYER FLOWS HINCHEY
BOUNDARY LAYER FLOWS HINCHEY BOUNDARY LAYER PHENOMENA When a body moves through a viscous fluid, the fluid at its surface moves with it. It does not slip over the surface. When a body moves at high speed,
More informationNumerical investigation of the flow instabilities in centrifugal fan
Proceedings of the 4th WSEAS International Conference on Fluid Mechanics and Aerodynamics, Elounda, Greece, August 21-23, 26 (pp282-288) Numerical investigation of the flow instabilities in centrifugal
More informationMODELLING TURBULENT HEAT FLUXES USING THE ELLIPTIC BLENDING APPROACH FOR NATURAL CONVECTION
MODELLING TURBULENT HEAT FLUXES USING THE ELLIPTIC BLENDING APPROACH FOR NATURAL CONVECTION F. Dehoux Fluid Mechanics, Power generation and Environment Department MFEE Dept.) EDF R&D Chatou, France frederic.dehoux@edf.fr
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationA combined application of the integral wall model and the rough wall rescaling-recycling method
AIAA 25-299 A combined application of the integral wall model and the rough wall rescaling-recycling method X.I.A. Yang J. Sadique R. Mittal C. Meneveau Johns Hopkins University, Baltimore, MD, 228, USA
More informationSimulation analysis using CFD on vibration behaviors of circular cylinders subjected to free jets through narrow gaps in the vicinity of walls
Fluid Structure Interaction V 85 Simulation analysis using CFD on vibration behaviors of circular cylinders subjected to free jets through narrow gaps in the vicinity of walls K. Fujita Osaka City University,
More informationFlow structure and modeling issues in the closure region of attached cavitation
PHYSICS OF FLUIDS VOLUME 12, NUMBER 4 APRIL 2000 Flow structure and modeling issues in the closure region of attached cavitation Shridhar Gopalan and Joseph Katz Department of Mechanical Engineering, The
More informationA Discussion of Low Reynolds Number Flow for the Two-Dimensional Benchmark Test Case
A Discussion of Low Reynolds Number Flow for the Two-Dimensional Benchmark Test Case M. Weng, P. V. Nielsen and L. Liu Aalborg University Introduction. The use of CFD in ventilation research has arrived
More informationIntroduction to CFD modelling of source terms and local-scale atmospheric dispersion (Part 1 of 2)
1 Introduction to CFD modelling of source terms and local-scale atmospheric dispersion (Part 1 of 2) Atmospheric Dispersion Modelling Liaison Committee (ADMLC) meeting 15 February 2018 Simon Gant, Fluid
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 informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationarxiv: v1 [physics.flu-dyn] 4 Aug 2014
A hybrid RANS/LES framework to investigate spatially developing turbulent boundary layers arxiv:1408.1060v1 [physics.flu-dyn] 4 Aug 2014 Sunil K. Arolla a,1, a Sibley School of Mechanical and Aerospace
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
More informationDetailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition
ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and
More informationExplicit algebraic Reynolds stress models for internal flows
5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence
More informationSG Turbulence models for CFD
SG2218 2012 Turbulence models for CFD Stefan Wallin Linné FLOW Centre Dept of Mechanics, KTH Dept. of Aeronautics and Systems Integration, FOI There are no simple turbulent flows Turbulent boundary layer:
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
More informationHybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model
Hybrid RANS Method Based on an Explicit Algebraic Reynolds Stress Model Benoit Jaffrézic, Michael Breuer and Antonio Delgado Institute of Fluid Mechanics, LSTM University of Nürnberg bjaffrez/breuer@lstm.uni-erlangen.de
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 informationASSESSMENT OF RANS AND DES METHODS FOR THE AHMED BODY
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More informationNumerical Study of Jet Plume Instability from an Overexpanded Nozzle
45th AIAA Aerospace Sciences Meeting and Exhibit 8 - January 27, Reno, Nevada AIAA 27-39 Numerical Study of Jet Plume Instability from an Overexpanded Nozzle Q. Xiao * and H.M. Tsai Temasek Laboratories,
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 informationInvestigation of three-dimensional effects on a cavitating Venturi flow
Investigation of three-dimensional effects on a cavitating Venturi flow Jean Decaix, Eric Goncalvès da Silva To cite this version: Jean Decaix, Eric Goncalvès da Silva. Investigation of three-dimensional
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 informationDirect Numerical Simulations of Transitional Flow in Turbomachinery
Direct Numerical Simulations of Transitional Flow in Turbomachinery J.G. Wissink and W. Rodi Institute for Hydromechanics University of Karlsruhe Unsteady transitional flow over turbine blades Periodic
More informationNumerical modeling of complex turbulent flows
ISTP-16, 5, PRAGUE 16 TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA Numerical modeling of complex turbulent flows Karel Kozel Petr Louda Jaromír Příhoda Dept. of Technical Mathematics CTU Prague, Karlovo
More informationNumerical Analysis of Unsteady Open Water Characteristics of Surface Piercing Propeller
Third International Symposium on Marine Propulsors smp 13, Launceston, Tasmania, Australia, May 2013 Numerical Analysis of Unsteady Open Water Characteristics of Surface Piercing Propeller Kohei Himei
More informationWALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION
WALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION Joongnyon Kim, Kyoungyoun Kim, Hyung Jin Sung Department of Mechanical Engineering, Korea Advanced Institute of Science
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 informationDirect numerical simulation of interfacial instability in gas-liquid flows
Direct numerical simulation of interfacial instability in gas-liquid flows Iain Bethune 1, Lennon Ó Náraigh 2, David Scott 1, Peter Spelt 3,4, Prashant Valluri 5, Zlatko Solomenko 3 1 Edinburgh Parallel
More informationTwo-scale momentum theory for very large wind farms
Submitted to TORQUE 2016, Munich, Germany, 5-7 October 2016 Draft (2 nd May 2016) Two-scale momentum theory for very large wind farms Takafumi Nishino Centre for Offshore Renewable Energy Engineering,
More informationMeysam ATASHAFROOZ, Seyyed Abdolreza GANDJALIKHAN NASSAB, and Amir Babak ANSARI
THERMAL SCIENCE: Year 014, Vol. 18, No., pp. 479-49 479 NUMERICAL INVESTIGATION OF ENTROPY GENERATION IN LAMINAR FORCED CONVECTION FLOW OVER INCLINED BACKWARD AND FORWARD FACING STEPS IN A DUCT UNDER BLEEDING
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
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 informationURANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver
URANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver *Yohan Choe 1), Hyeongjun Kim 1), Chongam Kim 2) 1), 2) Department of Aerospace Engineering, Seoul National
More informationLarge-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall
Large-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall Y. Kim, I.P. Castro, and Z.T. Xie Introduction Wind turbines operate in the atmospheric boundary layer and their
More informationCOMPARATIVE STUDY BETWEEN FLOWS AROUND SPHERE AND POD USING FINITE VOLUME METHOD
Journal of Naval Architecture and Marine Engineering June, 2011 DOI: 10.3329/jname.v8i1.7388 http://www.bangljol.info COMPARATIVE STUDY BETWEEN FLOWS AROUND SPHERE AND POD USING FINITE VOLUME METHOD M.
More informationNonlinear shape evolution of immiscible two-phase interface
Nonlinear shape evolution of immiscible two-phase interface Francesco Capuano 1,2,*, Gennaro Coppola 1, Luigi de Luca 1 1 Dipartimento di Ingegneria Industriale (DII), Università di Napoli Federico II,
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
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 informationCFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 31 (2011) 33-41 CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE Sharmina Hussain Department of Mathematics and Natural Science BRAC University,
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 information