PREDICTION OF ICE CRYSTAL ACCRETION WITH IN-HOUSE TOOL TAICE
|
|
- Scott Palmer
- 6 years ago
- Views:
Transcription
1 PREDICTION OF ICE CRYSTAL ACCRETION WITH IN-HOUSE TOOL TAICE Erdem Ayan, Serkan Özgen, Erhan Tarhan, Murat Canıbek TAI Turkish Aerospace Industries Inc. SAE 2015 International Conference on Icing of Aircraft, Engines and Structures, June 22-25, 2015, Prague, Czech Republic
2 Outline Introduction TAI ice accretion prediction in house developed tool TAICE Modifications towards ice crystal accretion prediction Ice Crystal Accretion Drag Models Heat Transfer & Phase Change Model Impingement Model Secondary Particle Behaviour Conclusions & Future Work
3 Introduction Ice crystal ingestion to aircraft engines may cause ice to accrete on internal components, leading to flameout, mechanical damage, rollback, etc. Many in-flight incidents have occurred in the last decades due to engine failures especially at high altitude convective weather conditions. Being a solid substance, ice crystals may not cause any problem for airframe components at low temperatures, since they only bounce off when they hit the airframe. However, when they face warmer flow conditions, they can change phase and pose a threat for the aircraft. These warmer flow conditions generally occur around or within the engine, on heated components (i.e. pitot and aoa probes) where local temperatures are significantly higher than the freezing temperature
4 Introduction Numerous in-flight incidents, which have occurred generally in tropical regions and at high altitudes made researchers to investigate ice accretion due to the ice crystals. In recent years, there have been important enhancements in simulating SLD. However, ice accretion due to ice crystals has not yet been studied extensively. Thus, in the framework of HAIC FP7 European Project, the physical mechanisms of ice accretion on surfaces exposed to ice-crystals and mixedphase conditions are investigated. TAI will implement models related to the ice crystal accretion calculation to the existing ice accumulation prediction program for droplets, namely TAICE. The aim of the study is to present the current status of the work performed by TAI in this field within the context of drag coefficient models, heat transfer and phase change models and impingement models
5 TAI Ice Accretion Prediction Tool TAICE Hess-Smith Panel Method for flow field computation around wings & wing sections. Lagrangian approach for particle trajectory calculations. Integral Boundary Layer Method for convective heat transfer coefficient calculations. Extended Messinger Model for computation of ice accretion rates
6 Droplet to Ice Crystal In order to use the existing ice accretion prediction tool for ice crystals, following modifications are required. Non-spherical drag model implementation, Heat-transfer and phase change models, Impact & re-emission model, Secondary particle trajectory calculations, Modifications required for solving mixed phase glaciated icing conditions
7 Ice Crystal Particle Geometric Properties The calculation of sphericity differs for an oblate (E < 1) or prolate (E > 1) particle [1]. sphericity Φ = πdp2 A crosswise sphericity d p = 3 6v p π Φ = πrp2 A A denotes the surface area of the particle, d p refers to the volume equivalent diameter and v p is the volume of the particle The orientation of the particle with respect to the flow can be taken into consideration by crosswise sphericity. A is the projected area of the particle with respect to the flow. Aspect ratio (E) and eccentricity (e) E = d 1 d 2, e = 1 min d 1 d 2, d 2 d 1 2 where d 1 and d 2 are diameters of the particle parallel and perpendicular to the revolution axis, respectively. Thus, for oblate particles Φ = 4E2/3 2+ E2 1+e ln e 1 e for prolate cases Φ = 2E 2/3 1+E arcsin(e) e Finally if the particle is assumed to be in the most stable orientation during its trajectory, crosswise sphericity can be defined as follows for oblate and prolate cases, respectively: Φ = E 2/3 Φ = E 1/
8 Drag Coefficient Models Several models exist for drag coefficient calculation in the literature. Some of these models are applicable to non-spherical particles and take orientation change into consideration [1]: Haider & Levenspiel suggest a model in terms of Φ and Re p [2]: C d = 24 Re p (1 + a(re p b ) + cre p Re p +d Re p= ρ ad p v rel μ a = exp(2.4486φ Φ ) b = (0.5565Φ ) c = exp( Φ Φ Φ ) d = exp( φ Φ Φ ) Ganser model also takes crosswise sphericity into consideration [3]: C d = 24 Re ( (Re p ) Re p p Re p Re [ log Φ ] p = [ log Φ ]
9 Drag Coefficient Models Additional drag model is suggested by Hölzer and Sommerfeld. Their model also includes the effect of orientation of the particle [4]: C d = Re p Φ Re p Φ log Φ Re p Φ3/4 Φ A study is performed to evaluate the effect of drag coefficient modeling on collection efficiency of the icing cases. Aspect ratio and diameter of the particle are varied. Following running conditions are analyzed with these drag models. Run Diameter a (μm) Aspect ratio E Equivalent diameter d (μm) Φ Φ Reynolds number E E E E E E
10 Run # 1 Run # 2 Run # 3 When sphericity differs from 1, models predictions start to deviate from each other. When Φ < 0.6 this difference is much bigger. However, sphericity is not the only parameter that effects the drag coefficient. Generally, the higher the Stokes number, the lower the difference between the models. Run Diameter a (μm) Aspect ratio E Equivalent diameter d (μm) Φ Φ Reynolds number Stokes Number % Difference in Max. β Value E E E E E E
11 Run # 4 Run # 5 Run # 6 For lower sphericity values, Haider & Levenspiel model collection efficiency prediction is generally higher than the predictions of other drag models. Collection efficiency curves of the Ganser and Hölzer & Sommerfeld drag models behave similar to each other in the run conditions which are considered. Run Diameter a (μm) Aspect ratio E Equivalent diameter d (μm) Φ Φ Reynolds number Stokes Number % Difference in Max. β Value E E E E E E
12 Heat Transfer and Phase Change Models The general form of the heat equation is: m p c p dt p,m dt = πd p Nu Φ k a(t a T p,s ) For non-spherical particles, Nusselt number can be defined in terms of Re, Pr, Φ and Φ using the Richter correlation: 7.2 Nu = ΦRe 1/2 Pr 1/3 Φ Re 2/3 Pr 1/3 Φ Theoretically, this expression can be applied to all particles. However, for oblate particles and when aspect ratio is lower than 0.8, this relation yields unrealistic results. As it is suggested by Villedieu et al [1], Frössling correlation [5] for non-spherical particles can be extended by using the relationship between C d and Nu. Thus, the final equation is as follows: Nu = 2 Φ Pr 1/3 Φ 1/4 Re p 1/2 One obtains the Frössling correlation by setting Φ = 1: Nu = Pr 1/3 Re p 1/2 Φ
13 Heat Transfer and Phase Change Models According to the Mason s study, for ice crystal melting and evaporation, following assumptions can be made [6]: The ice core temperature is uniform and the particle begins to melt when its temperature has reached the melting temperature, T f. The ice core temperature remains constant and equal to T f during all the melting phase. The ice core and the surrounding liquid film are spherical and concentric. The heat transfer between the liquid film and the ice core is purely conductive. The ice particles do not shed their liquid water film during the melting phase. In this study, the specified model is used which is the extended version of Mason s model The model consists of three steps: In the first step when T p < T f, pure ice particle temperature increases up to melting temperature and mass exchange may occur at the particle surface by sublimation. Total mass and the final heat transfer equation is as follows: dm p Sh = m dt sub = πd p ρ Φ ad v,a y v,s y v,a dt p m p c i dt = πd Nu p Φ k a T a T p m sub L f + L v Sh = 2 Φ Sc 1/3 Φ 1/4 1/2 Re p ; Sherwood number, Sc = ν D; Schmidt number
14 Heat Transfer and Phase Change Models In the second step T p = T f, particle starts to melt and is surrounded by a concentric water layer. Mass exchange may occur at the particle surface by evaporation. This step continues until ice particle is completely melted. At the end of this step, particle becomes spherical. Evaporation rate m ev can be calculated with the following equation dm p dt = m Sh ev = πd p Φ ρ ad v,a (y v,s y v,a ) Then the final heat transfer equation is as follows: Nu πd p Φ k a T a T f = m ev L v + m f L f In the third step, when T p > T f droplet temperature increases from melting temperature and mass exchange may occur by evaporation. dt p m p c w dt = πd pnuk a T a T p m ev L v In all these steps, temperature inside the particle is supposed to be uniform
15 Case # Flow Velocity m/s, RH % Sphericity Particle init. diam. (μm) Flow Temp., o C Initial Temp., o C 1, , , 6, , 0.70, , 551, , , , , Implemented heat transfer and phase change model is validated with the Hauk s experiments which were performed within the HAIC project [7]. The maximum difference is around 15%. Experimental melting time for irregular shaped particles are smaller than the calculated melting time due to their higher surface area. Final liquid droplet diameters also seem compatible with each other. Maximum difference between the experimental and calculated droplet diameter is around 4%. The success of the models seems similar for both spherical and irregular shapes
16 Impingement Model Ice crystal impact on a surface may result in three scenarios depending on particle and surface properties [1]: Sticking regime: the particle totally adheres to the wall. Bouncing regime: the particle bounces off the wall upon impact. In this case the size and shape of the particle does not change but its velocity is altered. Shattering regime: The particle breaks down into smaller fragments. Depending on the presence of a liquid film on the wall or water inside the impinging particle, some fragments may be re-emitted into the flow, while others may adhere to the surface. The threshold between bouncing and shattering regimes depends on the particle kinetic energy before impact. A dimensionless number L, which is the ratio of the normal kinetic energy to the surface energy is defined: L = πρpd p 3 v pn 2 12 πe σ T p d p 2 = 1 12 ρ p d p v2 pn where v pn is the component of the impact velocity vector perpendicular to the wall and e σ is the surface energy per unit area. e σ T p
17 Impingement Model The current model takes three possible scenarios into account. Accordingly [1]: Case L 0.5 : quasi-elastic bouncing with no fracturing, Case 0.5 L 90 : inelastic bouncing with particle fragmentation and loss of kinetic energy. Case L 90 : highly inelastic impact with high particle fragmentation and high loss of kinetic energy. The most important effect of the presence of liquid water layer on the wall or within the particle is the probability of particle adhering partially or completely to the wall. Water layer thickness within the particle h p, can be calculated using the particle liquid water content: 2h p = d p d p,i = 6 π m p m p,i ρ w + m p,i ρ p,i The model is based on the following expression for P D : P D = 1 P B π m p,i ρ p,i 1 3 with P B = min 1, K B ξ n δ p + δ w The dimensionless film thicknesses δ p and δ w are δ p = h p d p and δ w = h w d p
18 Impingement Model K B is an adjustable parameter and ξ n is the normal restitution coefficient defined as: ξ n = ξ n B = 1 L 0.5 ξ n B = 0.5 L F ξ nn = 0.5 L L 90 L 90 In the Lagrangian approach, the model is implemented as follows: For each particle impact, P D is calculated. If P D = 1, the particle is assumed to totally adhere to the wall and there is no reemission of secondary particles. If P D < 1, a fraction of the impinging mass flow rate equal to P D is assumed to adhere to the wall and the remaining part equal to 1 P D is re-emitted to the air flow
19 Secondary Particle Behaviour If L 90, the impinging particle either bounces of the wall with a probability P B = 1 P D or adhere to the wall with a probability P D. If P B = 0, the particle is assumed to totally adhere to the wall. If P B > 0, the particle partially adheres to the wall and a secondary particle is re-emitted. If L 90, the impinging particle partially adheres to the wall or fragments into smaller particles, where P D is the ratio of the deposited mass to the incident mass. If P D = 1, the particle is assumed to totally adhere to the wall. If P D < 1, the particle partially adheres to the wall and a secondary particle is re-emitted to the air. Following secondary particle properties are evaluated depending on the calculated threshold value. o o o o o o Mass flow rate, Particle diameter, Solid and liquid mass fractions, Ice core density, Sphericity, Particle velocity
20 Diameter a (μm) Aspect ratio E Equivalent diameter d (μm) Φ Φ Reynolds number Stokes Number Mach Number E E Circular X-Section The model is based on the following expression for P D : P D = 1 P B with P B = min 1, K B ξ n δ p + δ w NACA 0012 K B strongly affects the collection efficiency curves. When the number of K B increases, the collection efficiency decreases. The effect of K B on collection efficiency curves differs depending on the test case even if all other flow parameters are same
21 Run Diameter a (μm) Aspect ratio E Eq. Dia. d (μm) Run Reynolds Number Stokes Number E E E E E E E E-01 For the analyzed cases, with increasing Reynolds Number, sphericity and crosswise sphericity particles spreads in a larger area. All of the drag models give similar results. However, Heider&Levenspiel results differs more from others. Φ Φ
22 Conclusion Ganser, Haider & Levenspiel, Hölzer & Sommerfeld Drag coefficient models, Heat Transfer & Phase Change and Impingement models, secondary particle trajectory calculations were successfully implemented to TAICE. Currently used models will be enhanced when the results of HAIC experiments become available. Collection efficiency curves which are obtained with different drag models for the particles having various diameter and aspect ratio are compared. Implemented heat transfer and phase change model results are compared with Hauk s experimental results. Melting time and final liquid diameter values are similar both for calculated and experimental results. Influence of the adjustment variable K B in impingement model is discussed. Thus, further calibration of K B is required for accurate ice accretion determination for each icing test case
23 References 1. Villedieu P., Trontin P., Chauvin R., Glaciated and mixed phase ice accretion modeling using ONERA 2D icing suite, 6 th AIAA Atmospheric and Space Environments Conference, Atlanta, AIAA , Haider A., Levenspiel O., Drag coefficient and terminal velocity of spherical and nonspherical particles, Powder Technology, Vol. 58, pp , Ganser, G.H., A rational approach to drag prediction of spherical and nonspherical particles, Powder Technology, Vol. 77, pp , Hölzer A., Sommerfeld M., New simple correlation formula for the drag coefficient of non-spherical particles, Powder Technology, Vol. 184, pp R-M, Frössling, N., The evaporation of falling drops (in German), Gerlands Beiträge zur Geophysik, Vol. 52, pp , Mason, B.J., On the melting of hailstones, Quart. J. Roy. Meteor. Soc., Vol. 82, pp , Hauk P., Roisman I., Tropea C., Investigation of the Melting Behaviour of Ice particles in an Acoustic Levitator, 6 th AIAA Atmospheric and Space Environments Conference, AIAA, Reston, VA,
In-Flight Mixed Phase Ice Accretion Prediction on Finite Wings with TAICE-3D
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND AEROSPACE SCIENCES (EUCASS) DOI: 339 In-Flight Mixed Phase Ice Accretion Prediction on Finite Wings with TAICE-3D Erdem Ayan and Serkan Özgen Turkish Aerospace
More informationAircraft Icing Icing Physics
Aircraft Icing Icing Physics Prof. Dr. Dept. Aerospace Engineering, METU Fall 2015 Outline Formation of ice in the atmosphere Supercooled water droplets Mechanism of aircraft icing Icing variations Ice
More informationPhysics of Aircraft Icing: A Predictive Challenge
Physics of Aircraft Icing: A Predictive Challenge Cameron Tropea Institute for Fluid Mechanics and Aerodynamics (SLA), Technische Universität Darmstadt, Germany MUSAF III ONERA, CERFACS 27-29 September
More informationEffect of Mixed Icing Conditions on Thermal Ice Protection Systems
Effect of Mixed Icing Conditions on Thermal Ice Protection Systems By Kamel Al-Khalil, Ph.D. Manager, LeClerc Icing Research Laboratory New York, NY 114 FAA Specialists Workshop I. INTRODUCTION The purpose
More informationAircraft Icing FAR/CS-25, Appendix O and P charts
Aircraft Icing FAR/CS-25, Appendix O and P charts Prof. Dr. Serkan ÖZGEN Dept. Aerospace Engineering, METU Fall 2015 Outline Appendix O and P - Background Existing CS-25 certification specifications for
More informationThe Ice Crystal Weather Threat to Engines
Jeanne Mason Boeing Commercial Airplanes The Ice Crystal Weather Threat to Engines BOEING is a trademark of Boeing Management Company. Filename.ppt 1 Agenda Introduction Recognition of engine power-loss
More informationIN-FLIGHT ICING SIMULATION WITH SUPERCOOLED LARGE DROPLET EFFECTS
HEFAT010 7 th International Conference on Heat Transfer Fluid Mechanics and Thermodynamics 19-1 July 010 Antalya Turkey IN-FLIGHT ICING SIMULATION WITH SUPERCOOLED LARGE DROPLET EFFECTS Özgen S.* and Canıbek
More informationAircraft Icing. FAR 25, Appendix C charts. Prof. Dr. Serkan ÖZGEN. Dept. Aerospace Engineering, METU Fall 2015
Aircraft Icing FAR 25, Appendix C charts Prof. Dr. Serkan ÖZGEN Dept. Aerospace Engineering, METU Fall 2015 Outline FAR 25 and FAR 29 Appendix C charts Using FAR 25 Appendix C charts Liquid water content
More informationGraupel and Hail Growth
Graupel and Hail Growth I. Growth of large ice particles In this section we look at some basics of graupeln and hail growth. Important components of graupeln and hail growth models include production of
More informationA Description of Convective Weather Containing Ice Crystals Associated with Engine Powerloss and Damage
A Description of Convective Weather Containing Ice Crystals Associated with Engine Powerloss and Damage The Boeing Company 1 Photo: courtesy of Ian McPherson The Boeing Company acknowledges the contributions
More informationTotal Temperature Measurements Using a Rearward Facing Probe in Supercooled Liquid Droplet and Ice Crystal Clouds
https://ntrs.nasa.gov/search.jsp?r=186883 18-11-9T16:4:47+:Z Total Temperature Measurements Using a Rearward Facing Probe in Supercooled Liquid Droplet and Ice Crystal Clouds Juan H. Agui, 1 and Peter
More informationEvaporator Design for an Isokinetic Total Water Content Probe in a Naturally Aspirating Configuration
16 th Australasian Fluid Mechanics Conference Crowne Plaza, Gold Coast, Australia -7 December 007 Evaporator Design for an Isokinetic Total Water Content Probe in a Naturally Aspirating Configuration D.R.
More informationBIOLOGICAL WARFARE AGENT SLURRY DROPLET METHODOLOGY
BIOLOGICAL WARFARE AGENT SLURRY DROPLET METHODOLOGY 30 th Annual AAAR Conference Rosen Shingle Creek, Orlando, Florida 3-7 October 2011 by Timothy Bauer, Z24 Naval Surface Warfare Center Dahlgren Division
More informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
More informationThe Ice Crystal Icing Hazard & Risk Mitigation: Delta Air Lines Perspective
The Ice Crystal Icing Hazard & Risk Mitigation: Delta Air Lines Perspective Bob Culver, Chief Line Check Pilot & Tom Fahey, Mgr. Meteorology Friends/Partners Aviation Weather (FPAW) 19 November 2015 Las
More informationNUMERICAL SIMULATION OF AIRFOIL ICE ACCRETION AND THERMAL ANTI-ICING SYSTEMS
24th INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES NUMERICAL SIMULATION OF AIRFOIL ICE ACCRETION AND THERMAL ANTI-ICING SYSTEMS Jeroen E Dillingh, Harry WM Hoeijmakers Group Engineering Fluid Dynamics,
More informationIn-flight Ice Accretion Prediction Code
In-flight Ice Accretion Prediction Code Vladimír HORÁK*, Zdeněk CHÁRA** *Corresponding author University of Defence in Brno, Kounicova 65, 612 00 Brno, Czech Republic vladimir.horak@unob.cz **Institute
More informationLeveraging STAR-CCM+ for Aircraft Applications. Durrell Rittenberg, Ph.D.
Leveraging STAR-CCM+ for Aircraft Applications Durrell Rittenberg, Ph.D. Overview of Icing with STAR-CCM+ Icing in aerospace Common applications Impact of icing on Aircraft safety Common icing conditions
More information2 D. Terminal velocity can be solved for by equating Fd and Fg Fg = 1/6πd 3 g ρ LIQ = 1/8 Cd π d 2 ρ air u
ecture 8 Collection Growth of iquid Hydrometeors. I. Terminal velocity The terminal velocity of an object is the maximum speed an object will accelerate to by gravity. At terminal velocity, a balance of
More informationDroplet behaviour in a Ranque-Hilsch vortex tube
Journal of Physics: Conference Series Droplet behaviour in a Ranque-Hilsch vortex tube To cite this article: R Liew et al 2 J. Phys.: Conf. Ser. 38 523 View the article online for updates and enhancements.
More informationRogers and Yau Chapter 10: Drop breakup, snow, precip rate, and bulk models
Rogers and Yau Chapter 10: Drop breakup, snow, precip rate, and bulk models One explanation for the negative exponential (M-P) distribution of raindrops is drop breakup. Drop size is limited because increased
More informationKessler warm rain microphysics scheme. Haifan Yan
Kessler warm rain microphysics scheme Haifan Yan INTRODUCTION One moment scheme,and many available bulk schemes have followed the approach of Kessler The purpose of the scheme is to increase understanding
More informationChapter 8 - Precipitation. Rain Drops, Cloud Droplets, and CCN
Chapter 8 - Precipitation Rain Drops, Cloud Droplets, and CCN Recall the relative sizes of rain drops, cloud drops, and CCN: raindrops - 2000 μ m = 2 mm fall at a speed of 4-5 ms -1 cloud drops - 20 μ
More informationChapter 17. Fundamentals of Atmospheric Modeling
Overhead Slides for Chapter 17 of Fundamentals of Atmospheric Modeling by Mark Z. Jacobson Department of Civil & Environmental Engineering Stanford University Stanford, CA 94305-4020 August 21, 1998 Mass
More informationStudy on the Heat Transfer Characteristics in aircraft icing
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (2014) 000 000 www.elsevier.com/locate/procedia APISAT2014, 2014 Asia-Pacific International Symposium on Aerospace Technology,
More informationStrategy in modelling irregular shaped particle behaviour in confined turbulent flows
Title Strategy in modelling irregular shaped particle behaviour in confined turbulent flows M. Sommerfeld F L Mechanische Verfahrenstechnik Zentrum Ingenieurwissenschaften 699 Halle (Saale), Germany www-mvt.iw.uni-halle.de
More informationChapter 7: External Forced Convection. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University
Chapter 7: External Forced Convection Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University Objectives When you finish studying this chapter, you should be able to: Distinguish between
More informationSummary of Recent Modelling Work at Cambridge. J. Dean, C. Taltavull, P. Earp & T. W. Clyne
Summary of Recent Modelling Work at Cambridge J. Dean, C. Taltavull, P. Earp & T. W. Clyne Experimental Overview Powder Feed Thermocouples Plasma Gun Pitot Tube Target Substrate Experimental Data Thermal
More informationPhysical Science Exam 3 Study Guide. Dr. Karoline Rostamiani. Chapter 3
Chapter 3 Section 1 States of Matter What is matter made of? What are the three most common states of matter? How do particles behave in each state of matter? Solids, Liquids, and Gases Materials can be
More informationSummary of riming onset conditions for different crystal habits. Semi-dimension: width / lateral dimension (perpendicular to c-axis)
Summary of riming onset conditions for different crystal habits Semi-dimension: width / lateral dimension (perpendicular to c-axis) HEAT BALANCE FOR GRAUPEL PARTICLES Consider a graupel particle growing
More informationTRANSIENT BEHAVIOR IN AXIAL COMPRESSORS IN EVENT OF ICE SHED
Proceedings of ASME Turbo Expo 2015: Turbine Technical Conference and Exposition GT2015 June 15 19, 2015, Montréal, Canada GT2015-42413 TRANSIENT BEHAVIOR IN AXIAL COMPRESSORS IN EVENT OF ICE SHED Swati
More informationAssessing the Likelihood of Hail Impact Damage on Wind Turbine Blades
Assessing the Likelihood of Hail Impact Damage on Wind Turbine Blades Hamish Macdonald Prof. Margaret Stack, Prof. David Nash Environmental Conditions Temperature Salinity in the air UV radiation from
More informationChapter 7: Precipitation Processes. ESS5 Prof. Jin-Yi Yu
Chapter 7: Precipitation Processes From: Introduction to Tropical Meteorology, 1st Edition, Version 1.1.2, Produced by the COMET Program Copyright 2007-2008, 2008, University Corporation for Atmospheric
More informationPrecipitations. Terminal Velocity. Chapter 7: Precipitation Processes. Growth of Cloud Droplet Forms of Precipitations Cloud Seeding
Chapter 7: Precipitation Processes Precipitations Water Vapor Saturated Need cloud nuclei Cloud Droplet formed around Cloud Nuclei Growth of Cloud Droplet Forms of Precipitations Cloud Seeding Precipitation
More informationSupporting Information for Nonlinear Response of Iceberg Side Melting to Ocean Currents
GEOPHYSICAL RESEARCH LETTERS Supporting Information for Nonlinear Response of Iceberg Side Melting to Ocean Currents A. FitzMaurice 1, C. Cenedese 2, and F. Straneo 2 Contents of this file 1. Text S1 to
More informationPrecipitation Processes METR σ is the surface tension, ρ l is the water density, R v is the Gas constant for water vapor, T is the air
Precipitation Processes METR 2011 Introduction In order to grow things on earth, they need water. The way that the earth naturally irrigates is through snowfall and rainfall. Therefore, it is important
More informationAnalysis of hydrodynamic forces on non-spherical particles (Spherocylinder)
Coimbra, 6-8 March 2013. International workshop Fibre Suspension Flow Modelling French program ANR PLAYER Analysis of hydrodynamic forces on non-spherical particles (Spherocylinder) Rafik OUCHENE (LEMTA,
More informationIntroduction to Heat and Mass Transfer. Week 12
Introduction to Heat and Mass Transfer Week 12 Next Topic Convective Heat Transfer» Heat and Mass Transfer Analogy» Evaporative Cooling» Types of Flows Heat and Mass Transfer Analogy Equations governing
More informationChapter 3- Energy Balance and Temperature
Chapter 3- Energy Balance and Temperature Understanding Weather and Climate Aguado and Burt Influences on Insolation Absorption Reflection/Scattering Transmission 1 Absorption An absorber gains energy
More informationIn this chapter we explain the processes by which nonprecipitating cloud droplets and ice crystals grow large enough to fall as precipitation
Goals for today: 19 Oct., 2011 Ch 7, Precipitation Processes In this chapter we explain the processes by which nonprecipitating cloud droplets and ice crystals grow large enough to fall as precipitation
More informationICING SIMULATION ON JET ENGINE WITH TEMPERATURE CHANGE OF SUPER-COOLED DROPLET
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) Ryosuke Hayashi and
More informationConvective Vaporization and Burning of Fuel Droplet Arrays
Convective Vaporization and Burning of Fuel Droplet Arrays Guang Wu and William A. Sirignano Mechanical & Aerospace Engineering University of California, Irvine May 19, 2009 Background and Task In many
More informationCollision and Coalescence 3/3/2010. ATS 351 Lab 7 Precipitation. Droplet Growth by Collision and Coalescence. March 7, 2006
ATS 351 Lab 7 Precipitation March 7, 2006 Droplet Growth by Collision and Coalescence Growth by condensation alone takes too long ( 15 C -) Occurs in clouds with tops warmer than 5 F Greater the speed
More informationSummary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer
1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic
More informationFORMULA SHEET. General formulas:
FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to
More information1. describe the two methods by which cloud droplets can grow to produce precipitation (pp );
10 Precipitation Learning Goals After studying this chapter, students should be able to: 1. describe the two methods by which cloud droplets can grow to produce precipitation (pp. 232 236); 2. distinguish
More informationChapter 13 - States of Matter. Section 13.1 The nature of Gases
Chapter 13 - States of Matter Section 13.1 The nature of Gases Kinetic energy and gases Kinetic energy: the energy an object has because of its motion Kinetic theory: all matter is made if particles in
More informationHeat and Mass Transfer Unit-1 Conduction
1. State Fourier s Law of conduction. Heat and Mass Transfer Unit-1 Conduction Part-A The rate of heat conduction is proportional to the area measured normal to the direction of heat flow and to the temperature
More information3.1and 3.2 Thermal. Rise in temperature in deg C Final temperature in C A B C D
Name: Date: 3.1and 3.2 Thermal 1. During an experiment, a solid is heated from 285 K to 298 K. Which one of the following gives the rise in temperature, in deg C, and the final temperature, in C, of the
More informationModule 9: Mass Transfer Lecture 40: Analysis of Concentration Boundary Layer. The Lecture Contains: The concentration boundary layer
The Lecture Contains: The concentration boundary layer Heat and Mass Transfer Analogy Evaporate Cooling file:///d /Web%20Course%20(Ganesh%20Rana)/Dr.%20gautam%20biswas/Final/convective_heat_and_mass_transfer/lecture40/40_1.html[12/24/2014
More informationList of symbols. Latin symbols. Symbol Property Unit
Abstract Aircraft icing continues to be a threat for modern day aircraft. Icing occurs when supercooled large droplets (SLD s) impinge on the body of the aircraft. These droplets can bounce off, freeze
More informationModeling of dispersed phase by Lagrangian approach in Fluent
Lappeenranta University of Technology From the SelectedWorks of Kari Myöhänen 2008 Modeling of dispersed phase by Lagrangian approach in Fluent Kari Myöhänen Available at: https://works.bepress.com/kari_myohanen/5/
More informationModeling of Humidification in Comsol Multiphysics 4.4
Modeling of Humidification in Comsol Multiphysics 4.4 Indrajit Wadgaonkar *1 and Suresh Arikapudi 1 1 Tata Motors Ltd. Pimpri, Pune, India, 411018. *Corresponding author: Indrajit Wadgaonkar, Tata Motors
More informationObservation of Falling Motion of Columnar Snow Crystals
276 Journal of the Meteorological Society of Japan Vol. 54, No. 5 Observation of Falling Motion of Columnar Snow Crystals By Masahiro Kajikawa Department of Earth Science, Akita University, Akita, Japan
More informationChapter 7: External Forced Convection
Chapter 7: External Forced Convection Yoav Peles Department of Mechanical, Aerospace and Nuclear Engineering Rensselaer Polytechnic Institute Copyright The McGraw-Hill Companies, Inc. Permission required
More informationWMO Aeronautical Meteorology Scientific Conference 2017
Session 1 Science underpinning meteorological observations, forecasts, advisories and warnings 1.1 En route phenomena 1.1.1 Ice crystal icing, and airframe icing research Signatures of supercooled liquid
More informationSTATES OF MATTER. Chapter 3
STATES OF MATTER Chapter 3 Labs done so far for ch. 3 sections 1 and 2: 1. Distilled wood and related read of temperatures with plateaus for substances produced 2. Distilling solution X (BP/CP evaporation/condensation)
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II External Flows 1 Introduction In this chapter we will consider several fundamental flows, namely: the flat plate, the cylinder, the sphere, several other body shapes, and banks
More informationExam 2: Cloud Physics April 16, 2008 Physical Meteorology Questions 1-10 are worth 5 points each. Questions are worth 10 points each.
Exam : Cloud Physics April, 8 Physical Meteorology 344 Name Questions - are worth 5 points each. Questions -5 are worth points each.. Rank the concentrations of the following from lowest () to highest
More informationEverything in the universe can be classified as either matter or energy. Kinetic Energy Theory: All particles of matter are in constant motion.
Physical Science Everything in the universe can be classified as either matter or energy. Kinetic Energy Theory: All particles of matter are in constant motion. State of Matter Bose- Einstein Condensate
More information5. SPRAY/WALL IMPINGEMENT
5. SPRAY/WALL IMPINGEMENT 5.1 Wall Interaction Regimes Wachters and Westerling (1966), Akao et al. (1980), Senda et al. (1994) and Nagaoka et al. (1994) describe in detail the phenomena observed when drops
More informationMelting of ice particles:
Melting of ice particles: When ice particles fall below 0 C they begin to melt, but the process takes some time since heat transfer needs to occur (heat from ambient environment has to supply the latent
More informationExternal Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
External Forced Convection Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Drag and Heat Transfer in External flow Fluid flow over solid bodies is responsible
More informationFREEZING CONTAMINATION : AIRCRAFT ICING
FREEZING CONTAMINATION : AIRCRAFT ICING EFFECTS ON AIRCRAFT Different types of accretion Intensity of ice accretion Consequences of accretion Vulnerability factors examples Specific vulnerabilities Detection
More informationMECHANISMS OF CLOUD ELECTRIFICATION
8th AGILE MiniWorkshop The Third Birthday 28 April 2010, Bologna (Italy) MECHANISMS OF CLOUD ELECTRIFICATION F. Prodi [1,2]; F. Porcù [1,2]; M. Buiat [2] [1] ISACCNR Bologna [2] University of Ferrara,
More informationPrecipitation Processes. Precipitation Processes 2/24/11. Two Mechanisms that produce raindrops:
Precipitation is any form of water that falls from a cloud and reaches the ground. How do cloud drops grow? Chapter 7 When air is saturated with respect to a flat surface it is unsaturated with respect
More informationSolutions to questions from chapter 9 in GEF Cloud Physics
Solutions to questions from chapter 9 in GEF4310 - Cloud Physics i.h.h.karset@geo.uio.no Problem 1 What is the difference between collision, coalescence and collection? Problem 2 Name and very briefly
More informationParallel Plate Heat Exchanger
Parallel Plate Heat Exchanger Parallel Plate Heat Exchangers are use in a number of thermal processing applications. The characteristics are that the fluids flow in the narrow gap, between two parallel
More informationNumerical Study of Atmospheric Ice Accretion on Various Geometric Cross-sections
Numerical Study of Atmospheric Ice Accretion on Various Geometric Cross-sections by Muhammad S. Virk REPRINTED FROM WIND ENGINEERING VOLUME 35, NO. 5, 2011 MULTI-SCIENCE PUBLISHING COMPANY 5 WATES WAY
More informationA MODEL STUDY OF POWDER PARTICLE SIZE EFFECTS IN COLD SPRAY DEPOSITION
A ODEL STUDY OF POWDER PARTICLE SIZE EFFECTS IN COLD SPRAY DEPOSITION D. Helfritch and V. Champagne U.S. Army Research Laboratory, Aberdeen Proving Ground, D ABSTRACT The importance of high velocity for
More informationIntroduction. Effect of aerosols on precipitation: - challenging problem - no agreement between the results (quantitative and qualitative)
Introduction Atmospheric aerosols affect the cloud mycrophysical structure & formation (observations, numerical studies) An increase of the aerosol particles: - increases CCN concentrations - decreases
More informationLecture 28. Key words: Heat transfer, conduction, convection, radiation, furnace, heat transfer coefficient
Lecture 28 Contents Heat transfer importance Conduction Convection Free Convection Forced convection Radiation Radiation coefficient Illustration on heat transfer coefficient 1 Illustration on heat transfer
More informationCutting Fluid Effects on Machine Shop Air Quality
Cutting Fluid Effects on Machine Shop Air Quality For IAB000 Zhong Chen Advisor: Dr. Steven Y. Liang Oct. 18, 000 Problem Description Occupational Safety and Health Administration (OSHA) requirements 0.5
More informationA Lagrangian parcel based mixing plane method for calculating water based mixed phase particle flows in turbo-machinery
Comp. Part. Mech. (2015) 2:39 50 DOI 10.1007/s40571-015-0033-z A Lagrangian parcel based mixing plane method for calculating water based mixed phase particle flows in turbo-machinery Colin S. Bidwell Received:
More informationCHM Solids, Liquids, and Phase Changes (r15) Charles Taylor 1/9
CHM 111 - Solids, Liquids, and Phase Changes (r15) - 2015 Charles Taylor 1/9 Introduction In CHM 110, we used kinetic theory to explain the behavior of gases. Now, we will discuss solids and liquids. While
More informationFluid-Particles Interaction Models Asymptotics, Theory and Numerics I
Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I J. A. Carrillo collaborators: T. Goudon (Lille), P. Lafitte (Lille) and F. Vecil (UAB) (CPDE 2005), (JCP, 2008), (JSC, 2008) ICREA
More informationSeparating Cloud Forming Nuclei from Interstitial Aerosol
Chapter 16 Separating Cloud Forming Nuclei from Interstitial Aerosol Gourihar Kulkarni Additional information is available at the end of the chapter http://dx.doi.org/10.5772/50589 1. Introduction Our
More informationCombustion MATHEMATICAL MODEL FOR TRANSIENT. S. M. Frolov Λ,F.S.Frolov Λ, and B. Basara y
Combustion MATHEMATICAL MODEL FOR TRANSIENT DROPLET VAPORIZATION S. M. Frolov Λ,F.S.Frolov Λ, and B. Basara y Λ N. N. Semenov Institute of Chemical Physics Russian Academy of Sciences Moscow, Russia y
More informationCondensation, Evaporation, Weather, and the Water Cycle
Last name: First name: Date: Period: COM INC COM* Condensation, Evaporation, Weather, and the Water Cycle Evaporation and boiling are both types of vaporization, in which a liquid changes into a gas. The
More informationQuiz Review Topical Questions
Quiz Review Topical Questions Kinetic Theory of Matter Expansion and Contraction Solids, Liquids, Gases States of Matter Phase Changes Distillation Water Properties Kinetic Theory 1. The kinetic theory
More informationTutorial 1. Where Nu=(hl/k); Reynolds number Re=(Vlρ/µ) and Prandtl number Pr=(µCp/k)
Tutorial 1 1. Explain in detail the mechanism of forced convection. Show by dimensional analysis (Rayleigh method) that data for forced convection may be correlated by an equation of the form Nu = φ (Re,
More informationIntroduction to Cloud Microphysics
Introduction to Cloud Microphysics Mountain Weather and Climate ATM 619: Atmospheric Science Seminar Series Department of Earth and Atmospheric Sciences University at Albany W. James Steenburgh Department
More information8.6 Drag Forces in Fluids
86 Drag Forces in Fluids When a solid object moves through a fluid it will experience a resistive force, called the drag force, opposing its motion The fluid may be a liquid or a gas This force is a very
More informationTHE blades and vanes of a gas turbine engine are susceptible to
JOURNAL OF AIRCRAFT Analysis Method for Inertial Particle Separator Farooq Saeed and Ahmad Z. Al-Garni King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia DOI: 10.2514/1.20245 This
More informationChapter 3. States of Matter
Chapter 3 States of Matter 1. Solid 2. Liquid 3. Gas States of Matter Two More (discuss later) Plasma Bose-Einstein condensate States of Matter Solid (definite shape and volume) Particles are tightly packed
More informationPrecipitation. AT350: Ahrens Chapter 8
Precipitation AT350: Ahrens Chapter 8 Precipitation Formation How does precipitation form from tiny cloud drops? Warm rain process The Bergeron (ice crystal) process Most important at mid and northern
More informationScale-up problems are often perceived as difficult. Here the reaction calorimetry has proven to be
APPLICATION OF REACTION CALORIMETRY FOR THE SOLUTION OF SCALE-UP PROBLEMS A paper from the RC User Forum Europe, Interlaken, 1995 Francis Stoessel, Ciba AG, Basel, Switzerland. Scale-up problems are often
More informationPHASE CHANGE. Freezing Sublimation
Melting Graphic Organizer Deposition PHASE CHANGE Freezing Sublimation Boiling Evaporation Condensation PHASE CHANGE Phase change happens as the temperature changes. All matter can move from one state
More informationEmpirical Co - Relations approach for solving problems of convection 10:06:43
Empirical Co - Relations approach for solving problems of convection 10:06:43 10:06:44 Empirical Corelations for Free Convection Use T f or T b for getting various properties like Re = VL c / ν β = thermal
More informationChemistry A: States of Matter Packet Name: Hour: Page!1. Chemistry A States of Matter Packet
Chemistry A: States of Matter Packet Name: Hour: Page!1 Chemistry A States of Matter Packet Chemistry A: States of Matter Packet Name: Hour: Page!2 Worksheet #1: States of Matter In this packet we will
More informationMECHANICAL PROPERTIES OF FLUIDS:
Important Definitions: MECHANICAL PROPERTIES OF FLUIDS: Fluid: A substance that can flow is called Fluid Both liquids and gases are fluids Pressure: The normal force acting per unit area of a surface is
More informationFederal Aviation Administration
High Ice Water Content Research International Collaboration and Field Campaign Presented to: FPAW By: Tom Bond Date: Statement of Need Over the past 10+ years, it has been recognized that jet engine powerloss
More informationComputational Study of Sprays for the Development of a Monte Carlo Model
38th Dayton-Cincinnati Aerospace Sciences Symposium Computational Study of Sprays for the Development of a Monte Carlo Model Presenter: Murat Dinc West Virginia University Donald D. Gray West Virginia
More informationWhy Should You Consider a Freezing Point Depressant Ice Protection System? Icing Certification Present and Future. CAV Aerospace Limited
Why Should You Consider a Freezing Point Depressant Ice Protection System? Icing Certification Present and Future OVERVIEW Presentation Aircraft Review of Current Icing Environments Pending Changes to
More informationWater in the Atmosphere
Water in the Atmosphere Characteristics of Water solid state at 0 o C or below (appearing as ice, snow, hail and ice crystals) liquid state between 0 o C and 100 o C (appearing as rain and cloud droplets)
More informationLagrangian Particle Tracking
FD with OpenFOAM software Lagrangian article Tracking Jelena Andric December 2009, Gothenburg About multiphase flows great importance can occur even more frequently than single phase flows correct formulation
More informationChemistry B11 Chapter 6 Gases, Liquids, and Solids
Chapter 6 Gases, Liquids, and Solids States of matter: the physical state of matter depends on a balance between the kinetic energy of particles, which tends to keep them apart, and the attractive forces
More informationOn the Coefficients in Meteor Physics Equations
On the Coefficients in Meteor Physics Equations aria Yu. Khanukaeva epartment of Applied Mathematics, Moscow Institute of Physics and Technology (State University), Institutsky 9, 141700 olgoprudny, Moscow
More informationDimensionless Numbers
1 06.10.2017, 09:49 Dimensionless Numbers A. Salih Dept. of Aerospace Engineering IIST, Thiruvananthapuram The nondimensionalization of the governing equations of fluid flow is important for both theoretical
More information* Defining Temperature * Temperature is proportional to the kinetic energy of atoms and molecules. * Temperature * Internal energy
* Defining Temperature * We associate temperature with how hot or cold an object feels. * Our sense of touch serves as a qualitative indicator of temperature. * Energy must be either added or removed from
More information