CONTINUOUS THERMODYNAMICS FINITE DIFFUSION MODEL FOR MULTICOMPONENT FUEL SPRAY EVAPORATION
|
|
- Felix Reynolds
- 5 years ago
- Views:
Transcription
1 CONTINUOUS THERMODYNAMICS FINITE DIFFUSION MODEL FOR MULTICOMPONENT FUEL SPRAY EVAPORATION Dongyao Wang 1 and Chia-fon F. Lee Departent of Theoretical and Applied Mechanics 1 Departent of Mechanical and Industrial Engineering University of Illinois at Urbana-Chapaign Urbana, IL INTRODUCTION The different volatility of the coponents of a fuel ixture will result in different vaporization rates of coponents at the droplet surface, and the fuel becoes non-unifor inside the liquid phase because of the ass diffusion resistance inside the droplet. For practical, ultidiensional siulations of IC engines, zero-diensional odels for the liquid phase are necessary because of their low coputational cost and the fact that thousands of coputational parcels are typically needed to represent the fuel spray in an engine. Zeng and Lee [1] developed a zero-diensional odel capable of including finite diffusion and preferential vaporization. This odel uses discrete coponent representation. Discrete ulticoponent fuel odels use a set of coponents to atch the distillation behavior of a real fuel. A large nuber of coponents are typically required in order to represent a coercial fuel adequately. However, the coputational cost increases quickly as the nuber of the coponents increases, since each coponent adds an additional transport equation. Coercial fuels usually contain hundreds of coponents, so it is usually not practical to represent every coponent using a discrete representation for spray siulation. An alternative approach is continuous representation, which uses a probability density function based on the principles of continuous therodynaics to represent the fuel. For engineering calculations, a typical coercial fuel consisting of hundreds of coponents can be reasonably represented by a continuous distribution function. As the ixture is assued of alkanes and properties of the species can be characterized by the variable of olecular weight, so only a few paraeters are required to describe the characteristics of the ixture, such as the ean and variance of the distribution. Thus, only the paraeters describing the ulticoponent fuel have to be deterined beforehand, rather than having to solve a transport equation for each coponent. The ethod of continuous therodynaics was previously used by Lippert and Reitz [], Yi et al.[3] and Zhu et al. [4] for ultidiensional engine odeling. These studies ade iportant contributions to the application of continuous therodynaics to IC engine siulations. However, these studies are based on the assuption of infinite diffusion in the liquid phase, so non-unifority inside the droplet (i.e., the effects of finite diffusion) is not considered. In this paper, a coprehensive odel considering preferential vaporization of a coplex fuel ixture using continuous distributions is presented. The odel consists of a gas phase sub-odel, which deterines evaporation fluxes, and a liquid phase sub-odel with finite diffusion. This odel was validated with experiental data for the vaporization of isolated JP-4 fuel droplet, and then it was used to analyze the vaporization behavior of single diesel fuel droplet. Finally, gasoline hollow-cone spray siulations were ade using this odel. This odel was ipleented into KIVA-3V for coputations. MATHEMATICAL MODEL Gas Phase Sub-Model A fuel ixture is described by a continuous distribution function, f(w), characterized by the olecular weight, w, of each coponent. In this study, f(w) is a Γ-distribution function, which is often used to represent petroleu fractions (Cotteran et al. [5]): α 1 ( w γ ) w γ f ( w) = exp( ), (1) α β Γ( α) β 1
2 where γ is the origin, α and β are paraeters controlling the shape of the distribution, with the ean olecular weight, θ=αβ+γ, and the variance, σ =αβ. In previous studies, using the second oent was found to be ore convenient for calculation than using the variance. Here, θ is the first oent of f(w), and =σ +θ is defined as the second oent. The evaporation fluxes are given by Equations ()-(4): Sh cd yfr yf N = ln(1 + B), B =, () R 1 yfr ~ D Nθ = yfr + (1 yfr ) θ RN D, (3) Dˆ N = yfr + (1 yfr ) R N, (4) D where the first and second oent fluxes, N θ and N, are associated with ass and variance fluxes. The effect of convection on vaporization is assued to be adequately expressed by a Sherwood nuber, Sh, and a Nusselt nuber, Nu, for olar flux and heat flux, respectively. The heat flux is obtained as: NC p ( Ts T ) q = H fg. (5) RNC p 1 exp( ) Nu k Liquid Phase Sub-Model The transient behavior of the liquid phase significantly affects the vaporization of a spray and thus, the vapor distribution in an IC engine. In the present study, an effort is ade to properly odel the nonunifority of the liquid phase due to preferential vaporization rather than using infinite diffusion odel (i.e., well ixed assuption) as in previous studies. It is assued that the effect of internal circulation can be expressed adequately by an effective diffusivity and conductivity [6,7]. It is reasonable to assue that every coponent in the liquid phase diffuses at the sae rate, though we consider the effect of olecular weight on the diffusion rate in the gas phase. As Abrazon and Sirignano [7] and Runge et al. [8] assued unifor liquid diffusivity for all the coponents of the liquid ixture, unifor liquid phase diffusivity was adopted in the present study. It is possible to find a relation between the surface and average properties of a droplet to approxiate the effect of finite diffusion rates in the liquid phase by using the sae approach as Zeng and Lee [1]. An expression for the difference between the surface and ean values can be obtained: dφ d λ1 E( Pev) Rφ Pev = 0. (1 + ) Φ d, (6) dt ts κ λ1 where, Φ d = Φ s Φ, λ 1 = , t s = R / κ, Pev = RR / κ and E(P ev ) is the coefficient of the effect of droplet surface regression, which is noralized to 1 at P ev = 0 and decreases as P ev increases; and ~ θ D Nθ R Nθ 1 Φ = = =, κ D ˆ, φ N R N. (7) c CpT k / ccp q / Cp The average quantities are deterined by: d dt R 3 g θ Nθ c = 3R N, (8) CpT q
3 where the fluxes are given by the gas phase sub-odel. Thus, the surface value, Φ s = Φ d + Φ can be obtained. The effective diffusivity and conductivity correlations of Abrazon and Sirignano [7] are used. The vapor-liquid equilibriu relations and fuel vapor property correlations given by Cotteran et al. [5] and Tai and Hallet [9] are used here. The present odel involves fuel liquid properties, so the proper correlations are introduced. The liquid viscosity correlation is based on Orrick and Erbar Method [10]. The Sastri Method [10] and the Hayduk and Minas Method [10] are extended in this odel to use for conductivity and diffusivity correlation in the liquid phase, respectively. The above odel was ipleented into a odified version of KIVA-3V [11]. With a continuous fuel representation, the vapor-liquid equilibriu and liquid phase quantities are olar based; however, original KIVA is ass based. Thus, olar and ass basis variables are converted in each tie step. RESULTS AND DISCUSSIONS Three cases of siulations were ade for validation of the present odel, analysis of single droplet vaporization behavior, and siulation of spray. The paraeters of distribution by Hallett [1] were used for our coputations. For diesel, gasoline and JP-4 fuels, the initial eans and variances are θ = 185 and σ = 43, θ = 85.5 and σ = 35.8 and θ = 113 and σ = 7.5, respectively, and with γ = 0. Coparisons With Experiental Data The present continuous distribution vaporization odel was used to siulate the evaporation of a JP-4 droplet. Runge et al. [8] easured the teperature of a JP-4 droplet in 94K air and oving at 3 /s. This experient showed that the droplet teperature initially dropped rapidly and then increased to about 93K. This behavior is due to the lighter coponents vaporizing rapidly and absorbing heat fro the liquid. Siulation using the present odel was ade at the sae conditions as those of the experient. Figure 1 shows that the coputational results agree well with the experiental data. Runge et al. [8] also showed coputational results with a discrete 33-coponent odel. Copared with their results, the present odel not only reduced the coputational cost draatically, but also iproved the agreeent with experient rearkably. Single Stagnant Droplet Vaporization Diesel droplet vaporization using continuous therodynaics was studied with a finite diffusion odel in the liquid phase. The coputations were ade using a single, 100 µ droplet with an initial teperature of 300 K, in abient air at 1 at and 1000 K. Figure shows the difference between the coposition at the droplet surface and the average coposition of the diesel fuel droplet when 50% of the fuel ass has evaporated. Since lighter coponents at the droplet surface vaporize faster than heavier coponents, the coposition at the surface is different fro that inside the droplet. If the infinite diffusion odel is used, the coposition of the droplet is considered unifor (i.e, the two curves will overlap). The infinite diffusion odel siplifies the proble enorously, but it neglects soe iportant physical processes. In the present odel, diffusion in the liquid phase is at a finite rate. The siulation predicts that the surface distribution has a higher ean olecular weight and a saller variance than the average distribution. Since the coposition of the droplet is varying, as the ore volatile coponents vaporize, the overall boiling point increases and the droplet never reaches an equilibriu teperature throughout its whole lifetie, which is an iportant feature of the distillation of practical fuels. Gasoline Spray The present odel was used to spray siulations under typical engine operating conditions suarized as follows. Abient conditions are 450 K and 4.4 bars, the total fuel ass injected is 0.01 gra, and the duration of injection is 1.3 s. A two-diensional cylindrical sector is used for coputational efficiency, with a grid of by 3 in the radial and axial directions, respectively. The initial spray cone angle is set as 45 degrees with a cone angle thickness of 10 degrees. Two cases 3
4 were investigated: one with iso-octane as the fuel with the vaporization odel using discrete representation [1], and the other with a ulti-coponent gasoline fuel using continuous representation and the present vaporization odel. The spray evaporation history is shown in Figure 3. Initially, the ulti-coponent fuel vaporizes slightly faster, but later on, it vaporizes uch slower than the single coponent case. This phenoenon is due to the lighter coponents vaporizing rapidly at first, while the heavier coponents with higher noral boiling points vaporize ore slowly. Since heavier coponents vaporize later, the overall fuel vapor ean olecular weight and average variance increase. The teporal variation of fuel vapor coposition is shown in Figure 4. The difference in vaporization behavior will naturally affect the global spray behavior. Droplet positions predicted in the two cases are shown in Figure 5. It is noticed that there are ore sall droplets in the gasoline spray than in the iso-octane spray, especially in the center area of the spray. In the ulticoponent case, a partially vaporized droplet consists of heavier coponents and thus can last longer. In the center of the spray, with a lower pressure, sall single coponent droplets can be depleted quickly, however, ulti-coponent droplets with heavier coponents ay persist. Another iportant difference is that the single coponent iso-octane spray tips see to collapse inward, whereas, the ulti-coponent spray tips stretch outward due to the heavier coposition of the droplets, which results in slower vaporization and higher penetration. SUMMARY AND CONCLUSIONS A preferential vaporization odel for coplex coercial fuels has been developed using continuous therodynaics. The fuel coposition is odeled using continuous distribution function, the effect of finite diffusion in liquid phase is odeled by the differences between surface and average values. Siulation of single JP-4 droplet vaporization agrees with the experient very well. Siulation of single diesel droplet shows that the present odel can capture the ajor features of coercial fuels. This odel is zero-diensional with a very low coputational cost, and so it could be eployed for siulations of a hollow-cone gasoline spray. Copared with an iso-octane spray, the spray structure of a ulti-coponent fuel is substantially different, which will lead to considerable effects on cobustion characteristics. ACKNOWLEDGEMENTS This work was supported in part by the National Science Foundation under grant No. CTS , with Farley Fisher as technical onitor, and by the Departent of Energy under grant No. 99- DOE SBC-L/P-0030DOE. We also thank Dar-Lon Chang and Yu-Yu Lee for their assistance in preparing the anuscript. REFERENCES [1] Zeng, Y. and Lee, C.F., ILASS Aericas 11th Annual Conference on Liquid Atoization and Spray Systes, Sacraento, CA, pp , (1998). [] Lippert, A.M. and Reitz, R.D., SAE Paper [3] Yi, Y., Zhu, G.S., and Reitz, R., SAE Paper [4] Zhu, G.S., Reitz, R., Xin, J., and Takabayashi, T., SAE Paper [5] Cotteran, R.L. and Prausnitz, J.M., Ind. Eng. Che. Process Des. Dev. 3:434 (1985). [6] Jin, D. and Boran, G.L., SAE Paper [7] Abrazon, B., and Sirignano, W.A., Int. J. Heat and Mass Transfer 3:1605 (1989). [8] Runge, T., Teske, M., and Polyeropoulos, C.E., Atoization & Sprays, 8:5-44(1998). [9] Tai, J., and Hallett, W.L.H., Che. Engng Sci. 50:933 (1995). [10] Reid, R.C., Prausnitz, J.M., and Poling, B.E., 1987, The Properties of Gases and Liquids (4 th edition), McGraw-Hill, New York. [11] Asden, A., KIVA-3V: A Block-Structured KIVA Progra for Engines with Vertical or Canted Valves, Los Alaos Natl. Lab. Rep. LA MS,
5 [1] Hallett, W.L.H., Cobust. Flaes 11:334 (000). NOMENCLATURE B Transfer nuber C p heat capacity D diffusion coefficient c Molar density f distribution function k theral conductivity H fg Enthalpy of vaporization, overall q Heat flux N olar vaporization flux, overall R droplet radius T teperature t Tie y Molar fraction w Molecular weight α,β,γ Paraeters of distribution function κ effective diffusivity θ Mean of distribution function Φ diffusion variable σ Variance of distribution function φ Boundary condition ψ nd oent of distribution function Subscript d difference between the surface and ean values F Fuel g gas phase liquid phase average value s surface value R Surface of droplet Abient Overhead -, ~, Average and weighted average over distribution function Average Surface 90 Teperature K 85 Coputation Experient Molar Fraction Tie s Figure 1. Coparison of droplet teperature for the coputation and the experient by Runge et al [8]. A JP-4 fuel droplet with initial diaeter of 646 µ oves at 3 /s in air at 94K and 1 at Molecular Weight Figure. Surface and average distributions of the droplet when 50% of the droplet ass has evaporated. 5
6 Gasoline Isooctane Fuel vapor ass (g) Coposition Mean Molecular Weight Variance Tie s Figure 3. Overall fuel vapor ass history of the spray for single coponent and continuous representation fuel odels Tie s Figure 4. Variation of fuel vapor coposition as the spray vaporizes using continuous fuel representation. Gasoline fuel with continuous representation 1 s 3 s 5 s Gasoline fuel with single coponent representation ( iso-octane) Figure 5. Positions of droplets of sprays predicted with single coponent and continuous representation fuel odels. 1
Developed Correlations for Prediction of The Enthalpies of Saturated Vapor Liquid Coexisting Phases
Nahrain University, College of Engineering Journal (NUCEJ) Vol.13 No.2, 2010 pp.116-128 Developed Correlations for Prediction of he Enthalpies of Saturated Vapor Liquid Coexisting Phases Mahoud Oar bdullah
More informationFuel Droplet Evaporation Modeling using Continuous Thermodynamics Method with Multi-Distribution Functions
International Multidimensional Engine Modeling User s Group Meeting at the SAE Congress, April 23, 2012, Detroit, MI Fuel Droplet Evaporation Modeling using Continuous Thermodynamics Method with Multi-Distribution
More informationChapter 4: Hypothesis of Diffusion-Limited Growth
Suary This section derives a useful equation to predict quantu dot size evolution under typical organoetallic synthesis conditions that are used to achieve narrow size distributions. Assuing diffusion-controlled
More informationUrbana-Champaign, USA. Tsinghua University, China; and
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 Finite Diffusion Multi-Components Fuel Droplet Vaporization Modeling
More informationAn Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period
An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction
More informationTHERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS = (1) ln = COMPOSITES & POLYCON 2009
Aerican Coposites Manufacturers Association January 15-17, 29 Tapa, FL USA Abstract THERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS by Thore M. Klaveness, Reichhold AS In
More informationKinetic Theory of Gases: Elementary Ideas
Kinetic Theory of Gases: Eleentary Ideas 17th February 2010 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion
More informationMODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION
MODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION A. L. Venezuela a, R. F. Cantão a, R. S. Ongaratto b, and R. N. Haneda c a Universidade Federal de São
More informationMeasuring Temperature with a Silicon Diode
Measuring Teperature with a Silicon Diode Due to the high sensitivity, nearly linear response, and easy availability, we will use a 1N4148 diode for the teperature transducer in our easureents 10 Analysis
More informationSIMULATION OF THE HEATING STEP WITHIN THE THERMOFORMING PROCESS USING THE FINITE DIFFERENCE METHOD
SIMULATION OF THE HEATING STEP WITHIN THE THERMOFORMING PROCESS USING THE FINITE DIFFERENCE METHOD A. Fertschej 1 *, G.R. Langecker 1 University of Leoben artur.fertschej@u-leoben.at; Franz-Josef Strasse
More informationSIMULATION OF SOIL TEMPERATURE VARIATION FOR GEOTHERMAL APPLICATIONS
International Journal of Mechanical Engineering and Technology (IJMET) Volue 8, Issue 2, February 2017, pp. 167 175 Article ID: IJMET_08_02_020 Available online at http://www.iaee.co/ijmet/issues.asp?jtype=ijmet&vtype=8&itype=2
More information1 (40) Gravitational Systems Two heavy spherical (radius 0.05R) objects are located at fixed positions along
(40) Gravitational Systes Two heavy spherical (radius 0.05) objects are located at fixed positions along 2M 2M 0 an axis in space. The first ass is centered at r = 0 and has a ass of 2M. The second ass
More informationKinetic Theory of Gases: Elementary Ideas
Kinetic Theory of Gases: Eleentary Ideas 9th February 011 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion of
More informationOptimum Design of Assembled Cavity Dies for Precision Forging Process
International Syposiu on Material, Energy and Environent Enginee (ISM3E 2015) Optiu Design of Assebled Cavity Dies for Precision Forging Process Jun-song Jin1 and Xin-yun Wang1,* 1 State Key Laboratory
More informationLecture #8-3 Oscillations, Simple Harmonic Motion
Lecture #8-3 Oscillations Siple Haronic Motion So far we have considered two basic types of otion: translation and rotation. But these are not the only two types of otion we can observe in every day life.
More informationThe Thermal Dependence and Urea Concentration Dependence of Rnase A Denaturant Transition
The Theral Dependence and Urea Concentration Dependence of Rnase A Denaturant Transition Bin LI Departent of Physics & Astronoy, University of Pittsburgh, Pittsburgh, PA 15260, U.S.A Feb.20 th, 2001 Abstract:
More informationREPRODUCING KERNEL PARTICLE METHOD FOR CONDUCTION RADIATION INTERACTION IN 3D PARTICIPATING MEDIA
Proceedings of the Asian Conference on Theral Sciences 2017, 1st ACTS March 26-30, 2017, Jeu Island, Korea ACTS-P00337 REPRODUCING KERNEL PARTICLE METHOD FOR CONDUCTION RADIATION INTERACTION IN 3D PARTICIPATING
More informationAn analytical relation between relaxation time spectrum and molecular weight distribution
An analytical relation between relaxation tie spectru and olecular weight distribution Wolfgang Thi, Christian Friedrich, a) Michael Marth, and Josef Honerkap b) Freiburger Materialforschungszentru, Stefan-Meier-Straße
More informationTitle. Author(s)Izumida, Yuki; Okuda, Koji. CitationPhysical review E, 80(2): Issue Date Doc URL. Rights. Type.
Title Onsager coefficients of a finite-tie Carnot cycle Author(s)Izuida, Yuki; Okuda, Koji CitationPhysical review E, 80(2): 021121 Issue Date 2009-08 Doc URL http://hdl.handle.net/2115/39348 Rights 2009
More informationCFD SIMULATION OF PRESSURE LOSS IN HVDC TRANSFORMER WINDING
Journal of Energy VOLUME 63 2014 journal hoepage: http://journalofenergy.co/ Ralf Wittaack Sieens AG, Nürnberg, Gerany Ralf.Wittaack@Sieens.co CFD SIMULATION OF PRESSURE LOSS IN HVDC TRANSFORMER WINDING
More informationThe Thermal Conductivity Theory of Non-uniform Granular Flow and the Mechanism Analysis
Coun. Theor. Phys. Beijing, China) 40 00) pp. 49 498 c International Acadeic Publishers Vol. 40, No. 4, October 5, 00 The Theral Conductivity Theory of Non-unifor Granular Flow and the Mechanis Analysis
More informationScattering and bound states
Chapter Scattering and bound states In this chapter we give a review of quantu-echanical scattering theory. We focus on the relation between the scattering aplitude of a potential and its bound states
More informationDensity and structure of undercooled liquid titanium
Article Condensed Matter Physics March 2012 Vol.57 No.7: 719 723 doi: 10.1007/s11434-011-4945-6 Density and structure of undercooled liquid titaniu WANG HaiPeng, YANG ShangJing & WEI BingBo * Departent
More informationANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER
IEPC 003-0034 ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER A. Bober, M. Guelan Asher Space Research Institute, Technion-Israel Institute of Technology, 3000 Haifa, Israel
More informationEFFECT OF THE NON-CONDENSABLE GAS TYPE DURING CONDENSATION OF WATER VAPOR
THERMAL SCIENCE: Year 217, Vol. 21, No. 6A, pp. 2457-2468 2457 EFFECT OF THE NON-CONDENSABLE GAS TYPE DURING CONDENSATION OF WATER VAPOR by Kaoutar ZINE-DINE *, Youness EL HAMMAMI, Rachid MIR, Touria MEDIOUNI,
More informationNumerical simulation of VOC Emissions from Dry Materials
Yang, X., Chen, Q., Zhang, J.S, Magee, R., Zeng, J. and Shaw,C.Y. 2001. Nuerical siulation of VOC eissions fro dry aterials, Building and Environent, 36(10), 1099-1107. Nuerical siulation of VOC Eissions
More informationAnalyzing Simulation Results
Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient
More informationExtension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels
Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique
More informationThe Development of a Generalised Jet Mixing Model for the Derivation of Exclusion Zones in Hazardous Area Classification
The Developent of a Generalised Jet Mixing Model for the Derivation of Exclusion Zones in Hazardous Area Classification Katherine McCobe, Peter Tait, Jason Whitley LogiCas Ltd. 33 Boundary Street, Spring
More informationUfuk Demirci* and Feza Kerestecioglu**
1 INDIRECT ADAPTIVE CONTROL OF MISSILES Ufuk Deirci* and Feza Kerestecioglu** *Turkish Navy Guided Missile Test Station, Beykoz, Istanbul, TURKEY **Departent of Electrical and Electronics Engineering,
More informationEVAPORATION EFFECT IN NONLINEAR PENETRATION OF HIGH ENERGY BEAM DRILLING
1 Journal of Marine Science and echnology, Vol. 17, No., pp. 1-17 (9) EVAPORAION EFFEC IN NONLINEAR PENERAION OF HIGH ENERGY BEAM DRILLING Je-Ee Ho* and Chen-Lung Yen** Key words: enthalpy ethod. ABSRAC
More informationPhysics 139B Solutions to Homework Set 3 Fall 2009
Physics 139B Solutions to Hoework Set 3 Fall 009 1. Consider a particle of ass attached to a rigid assless rod of fixed length R whose other end is fixed at the origin. The rod is free to rotate about
More informationEffect of Channel Geometries on Flow Reversal in Microchannel Evaporators
Purdue University Purdue e-pubs International Refrigeration and Air Conditioning Conference School of Mechanical Engineering 216 Effect of Channel Geoetries on Flow Reversal in Microchannel Evaporators
More informationChapter 1 Introduction and Kinetics of Particles
Chapter 1 Introduction and Kinetics of Particles 1.1 Introduction There are two ain approaches in siulating the transport equations (heat, ass, and oentu), continuu and discrete. In continuu approach,
More informationAssessment of wind-induced structural fatigue based on joint probability density function of wind speed and direction
The 1 World Congress on Advances in Civil, Environental, and Materials Research (ACEM 1) eoul, Korea, August 6-3, 1 Assessent of wind-induced structural fatigue based on oint probability density function
More informationDESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *
IJST, Transactions of Mechanical Engineering, Vol. 39, No. M1, pp 89-100 Printed in The Islaic Republic of Iran, 2015 Shira University DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS
More informationEVALUATION OF THERMAL CONDUCTIVITY IN PITCH- BASED CARBON FIBER REINFORCED PLASTICS
16 TH INTERNATIONA CONFERENCE ON COMPOSITE MATERIAS EVAUATION OF THERMA CONDUCTIVITY IN PITCH- BASED CARBON FIBER REINFORCED PASTICS Shinji Ogihara*, Makoto Yaaguchi**, Takahito Chiba**, Junichi Shiizu****,
More informationSpine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area
Proceedings of the 006 WSEAS/IASME International Conference on Heat and Mass Transfer, Miai, Florida, USA, January 18-0, 006 (pp13-18) Spine Fin Efficiency A Three Sided Pyraidal Fin of Equilateral Triangular
More informationProjectile Motion with Air Resistance (Numerical Modeling, Euler s Method)
Projectile Motion with Air Resistance (Nuerical Modeling, Euler s Method) Theory Euler s ethod is a siple way to approxiate the solution of ordinary differential equations (ode s) nuerically. Specifically,
More informationI affirm that I have never given nor received aid on this examination. I understand that cheating in the exam will result in a grade F for the class.
Che340 hysical Cheistry for Biocheists Exa 3 Apr 5, 0 Your Nae _ I affir that I have never given nor received aid on this exaination. I understand that cheating in the exa will result in a grade F for
More informationCFD SIMULATION OF A MEMBRANE DISTILLATION MODULE CHANNEL
CFD SIMULATION OF A MEMBRANE DISTILLATION MODULE CHANNEL A. Cipollina *, A. Di Miceli *, J. Koschikowski, G. Micale *, L. Rizzuti * * Dipartiento di Ingegneria Chiica dei Processi e dei Materiali, Università
More informationDistillation. The Continuous Column. Learning Outcomes. Recap - VLE for Meth H 2 O. Gavin Duffy School of Electrical Engineering DIT Kevin Street
Distillation The Continuous Colun Gavin Duffy School of Electrical Engineering DIT Kevin Street Learning Outcoes After this lecture you should be able to.. Describe how continuous distillation works List
More informationSolidification of Porous Material under Natural Convection by Three Phases Modeling
Solidification of Porous Material under Natural Convection by Three Phases Modeling Hassan Basirat Tabrizi, Meber, IAENG and F. Sadeghpour Abstract The perforance of natural convective flow over a rectangular
More informationKINETIC THEORY. Contents
KINETIC THEORY This brief paper on inetic theory deals with three topics: the hypotheses on which the theory is founded, the calculation of pressure and absolute teperature of an ideal gas and the principal
More informationMolecular dynamics algorithm for multiple time scales: Systems with long range forces
Molecular dynaics algorith for ultiple tie scales: Systes with long range forces Mark E. Tuckeran@ Bruce J. Berne Departent of Cheistry, Colubia University, New York, New York 10027 Glenn J. Martyna Departent
More informationAdsorption of 2,4-Xylidine on Granulated Activated Carbon in a Packed Bed Absorber. f f
Proceedings o the 6 th Nordic Syposiu on Petrophysics, 5-6 May, Trondhei, Norway Nordic Energy Research Prograe Norwegian U. o Science and Technology Adsorption o,4-xylidine on Granulated Activated Carbon
More informationHALF-LIFE TIME FOR VOC EMISSION AND SORPTION OF POROUS BUILDING MATERIALS
HAF-IFE TIME FOR VOC EMISSION AND SORPTION OF POROUS BUIDING MATERIAS C.-S. ee is Ph.D. Candidate, Departent of Building, Civil, and Environental Engineering; F. Haghighat is Professor, Departent of Building,
More informationModelling of the Through-air Bonding Process
Modelling of the Through-air Bonding Process M. Hossain 1, M. Acar, Ph.D. 2, W. Malalasekera 2 1 School of Engineering, The Robert Gordon University, Aberdeen, UNITED KINDOM 2 Mechanical and Manufacturing
More informationANALYSIS OF THE EFFECT OF THE CHEMICAL SPECIES CONCENTRATIONS ON THE RADIATION HEAT TRANSFER IN PARTICIPATING GASES USING A MONTE CARLO METHODOLOGY
Proceedings of the 11 th Brazilian Congress of Theral Sciences and Engineering -- ECIT 2006 Braz. Soc. of Mechanical Sciences and Engineering -- ABCM, Curitiba, Brazil,- Dec. 5-8, 2006 Paper CIT06-0366
More informationln P 1 saturation = T ln P 2 saturation = T
More Tutorial at www.littledubdoctor.co Physical Cheistry Answer each question in the space provided; use back of page if extra space is needed. Answer questions so the grader can READILY understand your
More informationMolecular Speeds. Real Gasses. Ideal Gas Law. Reasonable. Why the breakdown? P-V Diagram. Using moles. Using molecules
Kinetic Theory of Gases Connect icroscopic properties (kinetic energy and oentu) of olecules to acroscopic state properties of a gas (teperature and pressure). P v v 3 3 3 But K v and P kt K v kt Teperature
More informationOptical Properties of Plasmas of High-Z Elements
Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer
More informationPhysically Based Modeling CS Notes Spring 1997 Particle Collision and Contact
Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which
More informationA Simple Continuous Mixture Droplet Evaporation Model with Multiple Distribution Functions
Introduction A Simple Continuous Mixture Droplet Evaporation Model with Multiple Distribution Functions William Hallett and Claire Grimwood Dept. of Mechanical Engineering, University of Ottawa, Ottawa,
More informationTHE EFFECT OF SOLID PARTICLE SIZE UPON TIME AND SEDIMENTATION RATE
Bulletin of the Transilvania University of Braşov Series II: Forestry Wood Industry Agricultural Food Engineering Vol. 5 (54) No. 1-1 THE EFFECT OF SOLID PARTICLE SIZE UPON TIME AND SEDIMENTATION RATE
More informationSolutions to the problems in Chapter 6 and 7
Solutions to the probles in Chapter 6 and 7 6.3 Pressure of a Feri gas at zero teperature The nuber of electrons N and the internal energy U, inthevoluev,are N = V D(ε)f(ε)dε, U = V εd(ε)f(ε)dε, () The
More information2002 University of Porto, Faculty of Engineering (FEUP)
Holberg H, Ahtila P. Drying phenoenon in a fixed bed under the bio fuel ulti stage drying. In: Oliveira A, Afonso C, Riffat S, editors. Proceedings of the st International Conference on Sustainable Energy
More information1 Brownian motion and the Langevin equation
Figure 1: The robust appearance of Robert Brown (1773 1858) 1 Brownian otion and the Langevin equation In 1827, while exaining pollen grains and the spores of osses suspended in water under a icroscope,
More information1 Bounding the Margin
COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #12 Scribe: Jian Min Si March 14, 2013 1 Bounding the Margin We are continuing the proof of a bound on the generalization error of AdaBoost
More informationREDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION
ISSN 139 14X INFORMATION TECHNOLOGY AND CONTROL, 008, Vol.37, No.3 REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION Riantas Barauskas, Vidantas Riavičius Departent of Syste Analysis, Kaunas
More informationI. Concepts and Definitions. I. Concepts and Definitions
F. Properties of a syste (we use the to calculate changes in energy) 1. A property is a characteristic of a syste that can be given a nuerical value without considering the history of the syste. Exaples
More informationKinetic Molecular Theory of Ideal Gases
Lecture -3. Kinetic Molecular Theory of Ideal Gases Last Lecture. IGL is a purely epirical law - solely the consequence of experiental obserations Explains the behaior of gases oer a liited range of conditions.
More informationCrystallization of Supercooled Liquid Elements Induced by Superclusters Containing Magic Atom Numbers Abstract: Keywords: 1.
Crystallization of Supercooled Liquid Eleents Induced by Superclusters Containing Magic Ato Nubers Robert F. Tournier, CRETA /CNRS, Université Joseph Fourier, B.P. 166, 804 Grenoble cedex 09, France. E-ail:
More information3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme
P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University
More informationSF Chemical Kinetics.
SF Cheical Kinetics. Lecture 5. Microscopic theory of cheical reaction inetics. Microscopic theories of cheical reaction inetics. basic ai is to calculate the rate constant for a cheical reaction fro first
More information4.5 Evaporation and Diffusion Evaporation and Diffusion through Quiescent Air (page 286) bulk motion of air and j. y a,2, y j,2 or P a,2, P j,2
4.5 Evaporation and Diffuion 4.5.4 Evaporation and Diffuion through Quiecent Air (page 86) z bul otion of air and j z diffuion of air (a) diffuion of containant (j) y a,, y j, or P a,, P j, z 1 volatile
More informationCFD-Parametric Study in Stator Heat Transfer of an Axial Flux Permanent Magnet Machine
World Acadey of Science, Engineering and Technology CFD-Paraetric Study in Stator Heat Transfer of an Axial Flux Peranent Magnet Machine Alireza Rasekh, Peter Sergeant, Jan Vierendeels Abstract This paper
More informationSolvent Extraction Research and Development, Japan, Vol. 22, No 1, (2015)
Solvent Extraction Research and Developent, Japan, Vol. 22, No 1, 95 11 (215) Notes Separation of Nitrogen Heterocyclic Copounds fro Model Oils by an Eulsion Liquid Mebrane Hiroaki HABAKI, Yoshiyuki SHIMADA
More informationHumidity parameters. Saturation (equilibrium) vapor pressure Condensation balances evaporation
uidity paraeters Saturation (equilibriu) vapor pressure Condensation balances evaporation Miing ratio & specific huidity Mass ratio of water vapor and air and water content and wet air. Dew point & frost
More informationChapter 2: Introduction to Damping in Free and Forced Vibrations
Chapter 2: Introduction to Daping in Free and Forced Vibrations This chapter ainly deals with the effect of daping in two conditions like free and forced excitation of echanical systes. Daping plays an
More informationHyperbolic Horn Helical Mass Spectrometer (3HMS) James G. Hagerman Hagerman Technology LLC & Pacific Environmental Technologies April 2005
Hyperbolic Horn Helical Mass Spectroeter (3HMS) Jaes G Hageran Hageran Technology LLC & Pacific Environental Technologies April 5 ABSTRACT This paper describes a new type of ass filter based on the REFIMS
More informationSchool of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, United Kingdom
pubs.acs.org/jpcb Calculation of Partition Functions and Free Energies of a Binary Mixture Using the Energy Partitioning Method: Application to Carbon Dioxide and Methane Haina Do,* Jonathan D. Hirst,
More informationKeywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution
Testing approxiate norality of an estiator using the estiated MSE and bias with an application to the shape paraeter of the generalized Pareto distribution J. Martin van Zyl Abstract In this work the norality
More informationSimplified Thermal Model of a Radial Flux Motor for Electric Vehicle Application
Journal of Asian Electric Vehicles, Volue 8, Nuber 2, Deceber 2010 Siplified Theral Model of a adial Flux Motor for Electric Vehicle Application Moez Ayadi 1, Hana Kebaili 2, Mohaed Aine Fakhfakh 3, and
More informationExperimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis
City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationCouncil for Innovative Research
ISSN 347-3487 Turbulent fil condensation in a vertical tube in presence of non condensable gas Y. Belkassi, K. Gueraoui, N. Hassanain, A. Elbouzidi Tea of vapour gas in fluid echanics and environent, LPT,
More informationDISTRIBUTION OF THE HYDRAULIC PARAMETERS AT RIVER BENDS
DISTRIBUTION OF THE HYDRAULIC PARAMETERS AT RIVER BENDS Isa Issa Oran *, Riyad Hassan Al-Anbari ** and Walaa Khalil Ali *** * Assist. Professor, Foundation of Technical Education ** Assist. Professor,
More informationMathematical Models of the Fluid Flowing for Geothermal and Hydrocarbon Wells
Proceedings World Geotheral Congress 005 Antalya, Turkey, 4-9 April 005 atheatical odels of the Fluid Flowing for Geotheral and Hydrocarbon Wells Olga B. Vereina Geological Institute of Russian Acadey
More informationUNCERTAINTIES IN THE APPLICATION OF ATMOSPHERIC AND ALTITUDE CORRECTIONS AS RECOMMENDED IN IEC STANDARDS
Paper Published on the16th International Syposiu on High Voltage Engineering, Cape Town, South Africa, 2009 UNCERTAINTIES IN THE APPLICATION OF ATMOSPHERIC AND ALTITUDE CORRECTIONS AS RECOMMENDED IN IEC
More informationMultiscale Entropy Analysis: A New Method to Detect Determinism in a Time. Series. A. Sarkar and P. Barat. Variable Energy Cyclotron Centre
Multiscale Entropy Analysis: A New Method to Detect Deterinis in a Tie Series A. Sarkar and P. Barat Variable Energy Cyclotron Centre /AF Bidhan Nagar, Kolkata 700064, India PACS nubers: 05.45.Tp, 89.75.-k,
More informationStatistical associating fluid theory-dimer of the solid phase of the pearl-necklace model
Indian Journal of Pure & Applied Physics Vol., epteber 00, pp 688-69 tatistical associating fluid theory-dier of the solid phase of the pearl-necklace odel Vivek Kuar ingh & K N Khanna Departent of Physics,
More informationBayes Decision Rule and Naïve Bayes Classifier
Bayes Decision Rule and Naïve Bayes Classifier Le Song Machine Learning I CSE 6740, Fall 2013 Gaussian Mixture odel A density odel p(x) ay be ulti-odal: odel it as a ixture of uni-odal distributions (e.g.
More informationProc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES
Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co
More informationSimulation of Geomechanical Behavior during SAGD Process using COMSOL Multiphysics
Siulation of Geoechanical Behavior during SAGD Process using COMSOL Multiphysics X. Gong 1, R.Wan *2 Departent of Civil Engineering, University of Calgary, Calgary, Alberta, Canada. *Corresponding author:
More informationJ. R. Barbosa, Jr. Departamento de Engenharia Mecânica Universidade Federal de Santa Catarina Florianópolis, SC. Brazil
J. R. Barbosa, Jr. Departaento de Engenharia Mecânica Universidade Federal de Santa Catarina 884-9 Florianópolis, SC. Brazil jrb@nrva.ufsc.br Two-Phase Non-Equilibriu Models: the Challenge of Iproving
More informationComparison of Stability of Selected Numerical Methods for Solving Stiff Semi- Linear Differential Equations
International Journal of Applied Science and Technology Vol. 7, No. 3, Septeber 217 Coparison of Stability of Selected Nuerical Methods for Solving Stiff Sei- Linear Differential Equations Kwaku Darkwah
More informationCelal S. Konor Release 1.1 (identical to 1.0) 3/21/08. 1-Hybrid isentropic-sigma vertical coordinate and governing equations in the free atmosphere
Celal S. Konor Release. (identical to.0) 3/2/08 -Hybrid isentropic-siga vertical coordinate governing equations in the free atosphere This section describes the equations in the free atosphere of the odel.
More informationDepartment of Building Science, Tsinghua University, Beijing , China
Supporting Inforation A SPME based C a -history ethod for easuring SVOC diffusion coefficients in clothing aterial Jianping Cao 1,, Ningrui Liu 1,, Yinping Zhang 1,,* 1 Departent of Building Science, Tsinghua
More informationA Numerical Investigation of Turbulent Magnetic Nanofluid Flow inside Square Straight Channel
A Nuerical Investigation of Turbulent Magnetic Nanofluid Flow inside Square Straight Channel M. R. Abdulwahab Technical College of Mosul, Mosul, Iraq ohaedalsafar2009@yahoo.co Abstract A nuerical study
More informationME 300 Thermodynamics II Exam 2 November 13, :00 p.m. 9:00 p.m.
ME 300 Therodynaics II Exa 2 Noveber 3, 202 8:00 p.. 9:00 p.. Nae: Solution Section (Circle One): Sojka Naik :30 a.. :30 p.. Instructions: This is a closed book/notes exa. You ay use a calculator. You
More informationMAKE SURE TA & TI STAMPS EVERY PAGE BEFORE YOU START
Laboratory Section: Last Revised on Deceber 15, 2014 Partners Naes: Grade: EXPERIMENT 8 Electron Beas 0. Pre-Laboratory Work [2 pts] 1. Nae the 2 forces that are equated in order to derive the charge to
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer
More informationTwo Dimensional Consolidations for Clay Soil of Non-Homogeneous and Anisotropic Permeability
Two Diensional Consolidations for Clay Soil of Non-Hoogeneous and Anisotropic Pereability Ressol R. Shakir, Muhaed Majeed Thiqar University, College of Engineering, Thiqar, Iraq University of Technology,
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 informationPremixed turbulent flames are of increasing practical importance
Nuerical siulation of a laboratory-scale turbulent V-flae J. B. Bell*, M. S. Day*, I. G. Shepherd, M. R. Johnson, R. K. Cheng, J. F. Grcar*, V. E. Beckner*, and M. J. Lijewski* *Center for Coputational
More information12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015
18.354J Nonlinear Dynaics II: Continuu Systes Lecture 12 Spring 2015 12 Towards hydrodynaic equations The previous classes focussed on the continuu description of static (tie-independent) elastic systes.
More informationMonitoring and system identification of suspension bridges: An alternative approach
Monitoring and syste identification of suspension bridges: An alternative approach Erdal Şafak Boğaziçi University, Kandilli Observatory and Earthquake Reseach Institute, Istanbul, Turkey Abstract This
More informationSupplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators
Suppleentary Inforation for Design of Bending Multi-Layer Electroactive Polyer Actuators Bavani Balakrisnan, Alek Nacev, and Elisabeth Sela University of Maryland, College Park, Maryland 074 1 Analytical
More informationFirst of all, because the base kets evolve according to the "wrong sign" Schrödinger equation (see pp ),
HW7.nb HW #7. Free particle path integral a) Propagator To siplify the notation, we write t t t, x x x and work in D. Since x i, p j i i j, we can just construct the 3D solution. First of all, because
More information