Droplet Vaporization in a Supercritical Environment. M. F. Trujillo Applied Research Laboratory Pennsylvania State University State College, PA 16803
|
|
- Penelope Day
- 6 years ago
- Views:
Transcription
1 ILASS Americas 19th Annual Conference on Liquid Atomization and Spray Systems, Toronto, Canada, May 2006 Droplet Vaporization in a Supercritical Environment M. F. Trujillo Applied Research Laboratory Pennsylvania State University State College, PA D.J. Torres Theoretical Division Los Alamos National Laboratory Los Alamos, NM Abstract A theoretical framework for the description of the physical processes involved in droplet vaporization leading towards the thermodynamic critical condition is proposed here. As opposed to previous work, an explicit expression for the di cult-to-determine thermal di usion coe cient is presented. This coe cient is responsible for the cross-e ect in heat and di usion fluxes and it is valid for non-ideal mixtures, which are typically encountered in the conditions of interest. Preliminary calculations of droplet vaporization at ambient pressures are shown to agree very well with micro-gravity data. At higher pressures a discrepancy in the heat-up period is noted, although the vaporization rate constants in the D 2 law portion of the droplet lifetime agree fairly well. Corresponding Author mft10@psu.edu
2 Introduction The description and modeling of droplet vaporization in supercritical environments has been the subject of a sizable number of investigations as documented in reviews written by Givler and Abraham [1], Yang [2], and Bellan [3]. While the behavior of vaporization prior to the attainment of a thermodynamic critical condition is well established, the phenomena occurring at the onset and after the interface has reached this point is much less understood. In particular, the physics governing the critical surface regression have received seemingly contradictory interpretations in the literature [4 6]. Our long term objectives are to study the onset critical point attainment taking into account as precisely as possible the anomalies encountered in transport coe cients and to consider the issue of thermodynamics instability for the entire droplet system. Prior to this, however, a complete theoretical framework needs to be put in place to incorporate all the dominant physics, in particular the cross-terms occurring in the heat and di usion fluxes, which have been shown to be important near the critical regime [7]. Among the various investigators that have contributed to the field, only Bellan and coworkers have advanced such a generalization of the physics. The general expressions for both the heat and mass fluxes that resulted from their work have already been known for quite some time in the nonequilibrium thermodynamics field. The new aspect was the application of this theory to the problem of droplet vaporization where it is expected that both the liquid and gas mixtures are non-ideal. The current work distinguishes itself from Bellan s group principally in the manner in which the essential thermal di usion coe cients are determined. These coe cients are responsible for the cross-e ects given by the Dufour and Soret terms. Borrowing from work aimed at hydrocarbon reservoirs [8, 9], an explicit expression for the thermal di usion coe - cient is presented that is valid for non-ideal mixtures. In what follows, the governing equations are described along with the interfacial conditions and the thermophysical properties. The expressions for heat and di usion flux are derived and the expression for the thermal di usion coe cient is developed. Preliminary results consisting of surface temperature and droplet diameter time series are presented for various conditions of temperature and pressure. Governing Equations The system under consideration is multicomponent with species =1,...,N, and with transport coe cients that are functions of pressure, temperature, and chemical composition. The equations are: for Y for species mass for momentum; + r ( u) + r ( Y u)= r j, + r ( v v) = rp + r, DT C p Dt ) ) P,Y Dt + : rv + X r q r j H, for energy. The di usion flux, j, and heat flux, q, originate from formulations found in nonequilibrium thermodynamics and are described below. Additionally the Equation-of-State (EoS) employed in this work is the Peng-Robinson EoS (PR- EoS) given by P = R ut b m a m 2 +2b m b 2, (1) m where the parameters a m and b m are calculated from mixing rules [10] and discussed in more detail in [11]. In order to isolate the core physics from geometrical complexities produced by buoyancy and/or convective forces, we consider vaporization in a micro-gravity environment. This leads to significant simplifications of the governing equations by a reduction to a spherical symmetry. As a consequence every vector and scalar field is only a function of the radial coordinate, i.e. u =(U, 0, 0), q =(q, 0, 0), and j =(j, 0, 0). This reduces the respective gov-
3 erning + r (r2 ( Y )+ r (r2 Y U)= r (r2 ( U)+ r (r2 U + @r r C ( T )+ r (r2 TU) ) ) P,Y Dt r (r2 q)+ X r (r2 j )H " 2 U + µ 2 r 3 r + 2U #. r Interfacial relations (2) The first of interfacial relation is the phaseequilibrium condition. In general, phase equilibrium is expressed by equating the chemical potential of species at the interface. This can be expressed as an equation for fugacities f g i = f l i, (3) where for the PR EoS, fugacities are given by, ln f i = b i (Z 1) ln(z B ) Px i b m P A N 2 p 2 j=1 x! j(1 k ji )(a j a i ) 1/2 b i 2B a m b m ln Z + (1 + p! 2)B Z ( p. (4) 2 1)B Additionally, we solved the following set of jump conditions for mass; g sy g,s(u g s g s( r s u g s) l s( r s u l s)=0, (5) Ṙ) l sy,s(u l l s r s )= j g + j l n, (6) for species; X L [ l sy (Ṙ l ul s) j ]+k l g rts g k l rts l +(P l P g )Ṙ = R utj A k T W avg Y A Y B g s R u Tj A k T W avg Y A Y B for energy. The interfacial constraints introduce the following (6+2N) unknown quantities; they are ( g s, l s, Ṙ, ug s,u l s,t s, Y, g Y )., l assuming that pressure is given, i.e. that the pressure change across the interface is of secondary importance (this assumption becomes even more accurate as the crit- µu r ical point is approached and the surface tension be- l s (7) gins to vanish). The interfacial system reduces to (6+2(N-1)) variables due to P Y g = P Y l = 1. For a binary system the result is 8 unknowns. The number of equations include two phase-equilibrium relations (3), jump condition for mass (5), for species (6), and energy (7), as well as the EoS (for both liquid and gas) giving 7 equations. Due to the excess of unknowns a variable will remain undetermined. This variable is chosen to be the liquid phase velocity which is solved iteratively until both governing equations and interfacial conditions are satisfied. Numerical Solution The numerical solution of the above equations follows an explicit version of an ALE method presented by Torres et al. [12]. We employ a staggered grid where all of the variables are defined at the cell center with the exception of the fluid velocity which is defined at the cell edges. Figure 1 shows this arrangement. A hyperbolic tangent is used to cluster the nodes close to the interface to allow the e cient capture of the thermal and composition boundary layers. Heat and Di usion Fluxes Most of the literature in droplet vaporization expresses the heat and di usion flux in terms of the Fourier heat conduction and Fick s Law of mass diffusion, respectively. Near the critical regime, however, it has been reported by Bellan and co-workers [13], that other modes of heat and mass di usion become important. To account for them consistently, we employ here the Maxwell-Stefan Form of these fluxes derived from irreversible thermodynamics [14, 15]. For a n-component mixture, the heat flux q is given by q = krt + NX =1 NX =1;6= NX =1 j h W + R u TX X D T j C j, (8) where j and j are the di usion fluxes for species and, respectively (to be described shortly). A constraint for these di usion fluxes is that they must satisfy P n =1 j = 0. Moreover, the multicomponent
4 thermal di usion coe cients obey the following relationship P n =1 DT = 0 [14], and the Maxwell-Stefan di usivities C are symmetric [15] (i.e.c = C ). For a binary mixture consisting of species A and B, the heat flux becomes q = krt + h A j A + h B j B W A W B + R utx A X B A + R utx B X A B DA T ja C AB A DB T jb C BA B j B B j A A. (9) Since the mixture is binary, we have j A = j B, DA T = DT B, and C AB = C BA which give ha hb q = krt + j A W A W B D T A 1 + R u TX A X B j A + 1. (10) C AB A B A B where = A + B. Simplifying this further we obtain a working version of the heat flux ha hb q = krt + j A W A W B + R utd T A C AB j A W avg A B W A W B. (11) This expression clearly shows the contributions from the traditional Fourier conduction term, the enthalpy di usion, and the Dufour term. With respect to the di usion flux, j, we begin with the definition of the barycentric or mass averaged velocity v, v = P N =1 v. The di usion flux is then defined as j = (v v). (12) This expression can be written in many cases in terms of the gradient of species concentration, resulting in the well known Fick s di usion form. In general, however, the di usion flux is influenced by other gradients besides mass fraction. A useful general expression is given by Bird et al. [15] relating the di erence in molecular velocities to these gradients. For zero external forces, we have NX =1;6= X X C (v v )= X (r ln a ) T,P 1 R u T (! )rp NX X X D T C =1;6= D T! r ln(t ) (13) For a binary mixture of A and B, this reduces to X A X B (v A v B )= X A (r ln a A ) T,P C AB 1 Y A )rp R u T ( A X A X B D T A C AB A From (12) and j A = v A v B = D T B B j B we have ja A r ln(t ) (14) j B B = j A A B. (15) Additionally, since DA T = DT B, (14) becomes j A = A B C h AB X A r ln a A X A X B + 1 R u T ( A Y A )rp + X AX B C AB D T A r ln T A B i. (16) Typically, we would like to deal with something more tangible than the gradient of the activity a A in the di usion flux equation. First we can rearrange this term to ln aa X A r ln a A = rx ln X A The activity is defined as the ratio a A = f A (T,P,X)/fA o (T,P), where f A(T,P,X) and fa o (T,P) are the fugacities at some arbitrary state (T,P,X) and a reference state, respectively. With this definition of activity and noting that fa o (T,P) is not a function of composition, the above expression ln fa (T,P,X) X A r ln a A = rx ln X A Substituting this into (16) and ignoring the expected weak contribution from the pressure gradient, yields j A = A B C AB X A X B + X AX B D T A C AB ln(fa ) rx ln X A A B r ln T i. (17) This expression can be further simplified by noting that rx A = WavgrY 2 A /(W A W B ) and performing some manipulations to give ln(f A ) j A = C AB ry A + D T rt A. ln X A T The binary di usion coe cient, D AB = C AB (@ ln(f A )/@ ln X A ), which is typically the
5 one provided in measurements, is equal to the Maxwell-Stevan di usivity in the case of ideal mixtures. The second term in this equation represents the Soret e ect. The major di culty in predicting this e ect is that for non-ideal mixtures, which are the type encountered in this work, the determination of the thermal di usion coe cient, DA T, is considerably di cult. Actually even in ideal mixtures the determination of DA T is far from easy [16]. Firoozabadi and co-workers [7 9] have provided a means of predicting the thermal di usion factor that can be related to the thermal di usion coe cient. This is examined below. Thermal Di usion Coe cient Under the same conditions examined above, i.e. no external forces (micro-gravity) and negligible pressure gradient, Ghorayeb and Firoozabadi [7] derive a similar expression to (18) for the di usion flux of a multicomponent mixture. Casting this relationship for a binary mixture yields ln f A j A = H AB X A W avg rx A W avg where and + H AB W avg K T A A i, (19) K T A = M AX A M B X B L 0 Aq W avg R u TL AA (20) H AB = L 0 Aq = L AA Q A W A W avgrl AA M A M B X A X B. (21) KA T is the thermal di usion ratio. The phenomenological coe cients are given respectively by L 0 Aq and L AA (their introduction comes from the theory of non-equilibrium thermodynamics where they appear in the linear relationship between the fluxes, j and q, and the thermodynamic forces rt, rx, rp [17].) An expression exists [8] between L 0 q and L,, which for a binary mixture reduces to. (22) Here Q = Q where [8] Q B W B H is the net heat of transport, q = nx =1 Q W j, (23) i.e. Q is the heat of transport of component. Substituting (22) into the equation for KA T gives KA T = M AX A M B X B Q A Q B. (24) W avg R u T W A W B This clearly shows as documented in [8] that for a binary mixture the thermal di usion ratio is independent of the Onsager phenomenological coe cients. Firoozabadi et al. [8] extends this statement to conclude that the thermal di usion coe cients themselves (D T ) are independent of the molecular di usion coe cients (D AB ) in the binary mixture case. However, as we will see shortly the expression for D T still retains L AA through H AB. Unless L AA can be expressed in terms of something other than a molecular di usion parameter, D T cannot be made independent from it. By comparing (18) with (20), we have D T A = W avg H AB K T A. (25) If we were to use (24) in this expression, L AA remains as previously mentioned. If instead we substitute (20), then L AA cancels, but L 0 Aq remains. To overcome this problem, another comparison between (18) and (20) yields A B C ln(f A ) = ln f A 1 H AB X A W avg. X A X ln X A W ln X A X A (26) This can be solved for H AB giving H AB = W AW B Wavg 2 C AB. (27) This equation along with (24) can be substituted into (25) resulting in a final expression for the thermal di usion coe cient DA T = W A W B W avg Wavg 2 C AB Q A M A X A M B X B W avg R u T W A Q B. (28) W B Relations for the net heat of transport, Q, have been developed by Firoozabadi and co-workers [8, 9] improving on the earlier work of Dougherty and Drickamer [18]. Namely, the net heat of transport for a binary mixture is given as apple XA Ī A A + X B ĪB B Q A = Ī A A + VA X A VA + X B VB, (29) where I is the partial molar internal energy departure and is the ratio of latent heat to the energy of viscous flow [19], for species respectively.
6 Thermophysical Properties For gas di usion, we employ the Chapman- Enskog dilute gas expression [10] modified for high pressure by using the Takahashi correction factor [20]. For liquid di usivity, the Hayduk and Minhas correlation is employed [21]. The method of Chung et al. [22] is used to calculate gas conductivity. Liquid conductivity is obtained from the fuel libraries of KIVA-4; it has been tested against data found in [10] for C 7 H 16. Latent heat of vaporization is calculated directly from the PR-EoS (details in [11]). Liquid specific heats are calculated based on a corresponding-of-states principle [10] and validated against data found in Vargaftik [23]. Gas specific heat are first obtained from ideal predictions [10] and corrected for high pressure using the internal energy departure calculated from the PR- EoS. Results The numerical scheme has been tested to ensure that the results presented here are numerically convergent with respect to time step and grid size. Due to the explicit nature of the numerical scheme, stability imposed a much higher restriction on t than accuracy. This meant that all of the calculations which successfully ran were already converged with respect to t. The calculations presented here are preliminary; the incorporation of the Soret and Dufour e ects previously described are still being incorporated into the code. Furthermore, the viscous dissipation term in the energy equation and the spatial variation of pressure are not yet accounted for. Nevertheless, for the conditions investigated in these initial runs the incorporation of these terms should not cause significant impacts. To validate our predictions, comparisons have been performed against the micro-gravity experimental data of Nomura et al. [24] for C 7 H 16 droplets vaporizing in a pure Nitrogen environment. In their experiments, a droplet in the diameter range of 0.6 to 0.8 mm and having an initial temperature of 300 o K, was generated at the tip of a silica fiber and introduced into an electric furnace. The elapsed time for the introduction was 0.16 seconds. In Figure 2 the surface temperature histories for cases ran at various ambient conditions are shown. As expected the rate of the temperature rise significantly increases with both ambient pressure and temperature. Comparisons of (d/d 2 o) as a function of time for the low pressure case are shown respectively in Figures 3 and 4. The agreement to experiments is particularly good at this pressure. As reported by Harstad and Bellan [25], one of the most unreliable transport coe cient predictions corresponds to the pressure e ect on mass di usivity. They proposed an alternate method for correcting for high pressure; in this work we have opted for the traditional Takahashi correction factor [20]. Comparisons can be made between both of these methods for the system being investigated here. In any case at ambient pressures, the predictions of the thermophysical properties are the most accurate; hence, the agreement shown in Figures 3 and 4 is expected. It should be noted that during the 0.16 seconds of elapsed time between the temperature measurement instant and the subsequent introduction into the chamber, the droplet was exposed to a convective environment caused by this translation which significantly enhances mass and heat transfer. The degree to which the temperature field was a ected can be estimated from the surface temperature history profiles shown in Figure 2; although because of the convective environment this plot actually shows a much milder rise. The key finding is that even though the elapse time is small, 5%, compared to the droplet lifetime, it occurs in a period of time when the rates of temperature increase are the largest. Hence, it is expected that for cases at higher pressures and temperatures, the initial field may deviate substantially from 300 K o. Additionally, other discrepancies are incurred by heat transfer through the supporting fiber [26] or from buoyancy e ects which become important at high pressure, even under micro-gravity conditions [3]. Motivated by this potential significant change in initial conditions, the next set of plots shown in Figures 5 (P=1.0 MPa, Ta=466K o ), 6 (P=2.0 MPa, Ta=452K o ), and 7 (P=2.0 MPa, Ta=656K o ) have an additional calculation performed with an initial droplet temperature of 340 K o. This additional calculation is not performed with the objective of replicating the exact details of the initial time period, but to show the impact on the results. In all of these latter plots, the lower temperature calculations show a considerably greater amount of initial droplet swelling. Once the droplets reaches the D 2 law region, the calculated vaporization rate constants are practically the same in both sets of predictions and agree reasonably well with experiments. Comparing the predictions between the fully resolved physics in both the gas and the liquid to those of only the liquid phase agree with the results of Zhu et al. [27]. It is observed that when fully resolving the gas phase and its unsteadiness, the calculations show a more pronounced vaporization behavior.
7 Summary The governing equations for the physics governing droplet phase change in supercritical environments have been advanced based on the generalization of the di usion and heat fluxes from irreversible thermodynamics [14, 15]. An explicit expression has also been obtained for the thermal diffusion coe cient which is responsible for the crosse ect between fluxes. Calculations for this are currently in progress. Computations for droplet diameter as a function of time have been performed at various pressures and temperatures. Good agreement is found at low pressures. At higher pressures a consistent deviation exists between the predictions and the data. This will be studied more thoroughly in the coming weeks. A potential source of experimental error, however, is the unaccounted convective transport occurring during a short, but significant initial elapsed time. This is almost negligible at lower pressures but becomes important at higher pressures.
8 Nomenclature a m,b m Peng-Robinson parameters C Maxwell-Stefan multicomponent di usivity C p Constant pressure specific heat C v Constant volume specific heat D Mass di usion coe cient D T Thermal di usion coe cient f Fugacity k Thermal conductivity H, h Enthalpy, Enthalpy per unit mass I Internal energy per unit mass j Di usion flux L Latent heat of vaporization M Mass P Pressure Pvap o Pure species vapor pressure q Heat flux r s Droplet radius r s Time rate of change of droplet radius R u Universal gas constant R Gas constant R u /W avg S Entropy per unit mass t Time T Temperature u, u Velocity V Volume v Droplet velocity Specific volume W Molecular weight X, X Mole fraction, mole fraction vector Y Mass fraction Z Compressibility factor Greek Specific heat ratio C p /C v µ Viscosity Density Viscous stress tensor! Acentric factor Subscripts, Corresponding to species, 1 Cell quantity c Critical r Reduced quantity, e.g. P r = P/P c o Reference state (ideal-gas) Superscripts X Specifies a molar quantity for X s Surface value g Gas phase l Liquid phase References [1] S.D. Givler and J. Abraham. Progress in Energy and Combustion Science, 22:1 28, [2] V. Yang. Proceedings of the Combustion Institute, 28: , [3] J. Bellan. Progress in Energy and Combustion Science, 26: , [4] P. Haldenwang, C. Nicoli, and J. Daou. Int. J. Heat Mass Transfer, 39: , [5] V. Yang, K.C. Hsieh, and J.S. Shuen. AIAA , 31st Aerospace Sciences Meeting and Exhibit, [6] W.A. Sirignano. Fluid Dynamics and Transport of Droplets and Sprays. Cambridge University Press, New York, [7] K. Ghorayeb and A. Firoozabadi. AIChE Journal, 46: , [8] A. Firoozabadi, K. Ghorayeb, and K. Shukla. AIChE Journal, 46: , [9] K. Shukla and A. Firoozabadi. Ind. Eng. Chem. Res., 37: , [10] B.E. Poling, J.M. Prausnitz, and J.P. O Connell. The Properties of Gases and Liquids, 5th Edition. Mc Graw Hill, New York, [11] M.F. Trujillo, D.J. Torres, and P.J. O Rourke. Int. J. Engine Res., 5: , [12] D.J. Torres, P.J. O Rourke, and A.A. Amsden. Combustion Theory Modelling, 7:67 86, [13] K. Harstad and J. Bellan. International Journal of Multiphase Flow, 26: , [14] C.F. Curtiss and R.B. Bird. Ind. Eng. Chem. Res., 38: , [15] R.B. Bird, W.E. Stewart, and E.N. Lightfoot. Transport Phenomena. John Wiley & Sons, Inc., 2 edition, [16] J.H. Ferziger and H.G. Kaper. Mathematical Theory of Transport Processes in Gases. North-Holland Publishing Company, Amsterdam, [17] S.R. de Groot and P. Mazur. Non-Equilibrium Thermodynamics. North-Holland Publishing Company, Amsterdam, [18] E.J. Dougherty and H.G. Drickamer. J. Chemical Physics, 23: , 1955.
9 [19] S. Glasstone, K.J. Laidler, and H. Eyring. The Theory of Rate Processes. McGraw-Hill, New York, [20] Takahashi S. J. Chem. Engng. Japan, 6: , [21] W. Hayduk and B.S. Minhas. Canadian Journal of Chemical Engineering, 60, [22] T-H. Chung, M. Ajlan, L.L. Lee, and K. Starling. Ind. Eng. Chem. Res., 27: , [23] N.B. Vargaftik. Tables on the Thermophysical properties of Liquids and Gases. Hemisphere Publishing Corporation, Washington, D.C., [24] H. Nomura, Y. Ujhe, H.J. Rath, J. Sato, and M. Kono. Twenty-Sixth Symposium on Combustion, pp , [25] K. Harstad and J. Bellan. Ind. Eng. Chem. Res., 43: , [26] C. Morin, C. Chauveau, P. Dagaut, I. Gokalp, and M. Cathonnet. Combust. Sci. and Tech., 176: , [27] G.S. Zhu, R.D. Reitz, and S.K. Aggarwal. International Journal of Heat and Mass Transfer, 44: , 2001.
10 Drop boundary Nodes nl 2 nl 1 nl nl+1 nl+2 nl+ng 1 nl+ng Cells nl 1 nl nl+1 nl+2 nl+ng 1 nl+ng Cell quantities: T, Y, Vol, and ρ Figure 1. Staggered grid arrangement used for numerical solution of drop vaporization. 420 Initial elapse period Droplet Surface Temperature (K) P= 0.1 MPa Ta=471 P=1.0 MPa Ta=466 K P=2.0 MPa Ta=452 K P=2.0 MPa, Ta= 656 K t/d 2 o (s/mm 2 ) Figure 2. Droplet surface temperature time series under various conditions.
11 P=0.1MPa, Ta=471 K Calculations Experiments (d/d o ) t/d 2 o (s/mm 2 ) Figure 3. Droplet size as a function of time, P=0.1 MPa, ambient temperature = 471 o K P=0.1 MPa Ta=741 K Calculations Experiments (d/d o ) t/d 2 o (s/mm 2 ) Figure 4. Droplet size as a function of time, P=0.1 MPa, ambient temperature = 741 o K.
12 1 0.8 (d/d o ) P=1 MPa Ta= 466 K 0.4 Calculations T init =340 K Calculations T init =300 K Experiment t/d 2 o (s/mm 2 ) Figure 5. Droplet size as a function of time, P=1.0 MPa, ambient temperature = 466 o K (d/d o ) Pa = 2.0 MPa and Ta=452 K 0.4 T_init=340K T_init=300 K Experiment t/d 2 o (s/mm 2 ) Figure 6. Droplet size as a function of time, P=2.0 MPa, ambient temperature = 452 o K.
13 1 0.8 (d/d o ) Pa= 2 MPa Ta=656 K Calculations T init =340 K Calculations T init =300 K Experiments t/d 2 o (s/mm 2 ) Figure 7. Droplet size as a function of time, P=2.0 MPa, ambient temperature = 656 o K.
Lecture 2 Constitutive Relations
Lecture 2 Constitutive Relations Conservation equations Dv = @ + r v =0 @t rp + r + g DY i + r ( Y iv i )= i i =1, 2,...,N { De = pr v + r q Dh + r q Z Y i p = R h = Y i h o i + c p d o W i We already
More informationSubcritical and Supercritical Droplet Evaporation within a Zero-Gravity Environment: Low Weber Number Relative Motion
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Mechanical & Materials Engineering Faculty Publications Mechanical & Materials Engineering, Department of 4-2008 Subcritical
More informationCFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM EVAPORATOR
Distillation Absorption 2010 A.B. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA2010 copyright notice CFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM
More informationFluid Dynamics and Balance Equations for Reacting Flows
Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and
More informationLecture 12. Droplet Combustion Spray Modeling. Moshe Matalon
Lecture 12 Droplet Combustion Spray Modeling Spray combustion: Many practical applications liquid fuel is injected into the combustion chamber resulting in fuel spray. Spray combustion involves many physical
More informationUnified Model for Nonideal Multicomponent Molecular Diffusion Coefficients
THERMODYNAMICS Unified Model for Nonideal Multicomponent Molecular Diffusion Coefficients Alana Leahy-Dios and Abbas Firoozabadi Chemical Engineering Dept., Mason Lab, Yale University, New Haven, CT 06520
More informationTHERMODYNAMIC CONSISTENCY TESTS FOR PHASE EQUILIBRIUM IN LIQUID SOLUTE+SUPERCRITICAL SOLVENT MIXTURES
THERMODYNAMIC CONSISTENCY TESTS FOR PHASE EQUILIBRIUM IN LIQUID SOLUTE+SUPERCRITICAL SOLVENT MIXTURES José O. Valderrama 1, and Víctor H. Alvarez 1 Fac. de Ingeniería, Depto. Ing. Mecánica, Univ. de la
More informationMulticomponent diffusion in gases and plasma mixtures
High Temperatures ^ High Pressures, 2002, volume 34, pages 109 ^ 116 15 ECTP Proceedings pages 1337 ^ 1344 DOI:10.1068/htwu73 Multicomponent diffusion in gases and plasma mixtures Irina A Sokolova Institute
More informationTransient supercritical droplet evaporation with emphasis on the e ects of equation of state
International Journal of Heat and Mass Transfer 43 (2000) 1157±1171 www.elsevier.com/locate/ijhmt Transient supercritical droplet evaporation with emphasis on the e ects of equation of state G.S. Zhu 1,
More informationReview of Chemical Equilibrium Introduction
Review of Chemical Equilibrium Introduction Copyright c 2015 by Nob Hill Publishing, LLC This chapter is a review of the equilibrium state of a system that can undergo chemical reaction Operating reactors
More informationFUNDAMENTALS OF CHEMISTRY Vol. II - Irreversible Processes: Phenomenological and Statistical Approach - Carlo Cercignani
IRREVERSIBLE PROCESSES: PHENOMENOLOGICAL AND STATISTICAL APPROACH Carlo Dipartimento di Matematica, Politecnico di Milano, Milano, Italy Keywords: Kinetic theory, thermodynamics, Boltzmann equation, Macroscopic
More informationNotes on Entropy Production in Multicomponent Fluids
Notes on Entropy Production in Multicomponent Fluids Robert F. Sekerka Updated January 2, 2001 from July 1993 Version Introduction We calculate the entropy production in a multicomponent fluid, allowing
More informationTHE PROPERTIES OF GASES AND LIQUIDS
THE PROPERTIES OF GASES AND LIQUIDS Bruce E. Poling University of Toledo John M. Prausnitz University of California at Berkeley John P. O'Connell University of Virginia Fifth Edition McGRAW-HILL New York
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 informationSimulating Condensation in a Supercritical Gas Jet
Simulating Condensation in a Supercritical Gas Jet L. Qiu, Y. Wang, H. Wang, Q. Jiao and R. D. Reitz * Engine Research Center University of Wisconsin Madison, WI 53706 USA Abstract The KIVA code was modified
More informationA Shrinking Core Model for Steam Hydration of CaO-Based Sorbents Cycled for CO2 Capture: Supplementary Information
A Shrinking Core Model for Steam Hydration of CaO-Based Sorbents Cycled for CO2 Capture: Supplementary Information John Blamey 1, Ming Zhao 2, Vasilije Manovic 3, Edward J. Anthony 3, Denis R. Dugwell
More informationChemical and Engineering Thermodynamics
Chemical and Engineering Thermodynamics Third Edition Stanley I. Sandler University of Delaware John Wiley & Sons, Inc. New York Chichester Weinheim Brisbane Singapore Toronto Contents NOTATION xv CHAPTER1
More informationarxiv: v2 [physics.flu-dyn] 4 May 2018
Transient Behavior near Liquid-Gas Interface at Supercritical Pressure Jordi Poblador-Ibanez 1,, William A. Sirignano 2 University of California, Irvine, CA 92697-3975, United States arxiv:1802.10149v2
More informationReview of Chemical Equilibrium Introduction
Review of Chemical Equilibrium Introduction Copyright c 2015 by Nob Hill Publishing, LLC This chapter is a review of the equilibrium state of a system that can undergo chemical reaction Operating reactors
More informationCHE 611 Advanced Chemical Reaction Engineering
CHE 611 Advanced Chemical Reaction Engineering Dr. Muhammad Rashid Usman Institute of Chemical Engineering and Technology University of the Punjab, Lahore 54590 mrusman.icet@pu.edu.pk 1 Diffusion and reaction
More informationOverall: 75 ECTS: 7.0
Course: Chemical Engineering Thermodynamics Language: English Lecturer: Prof. dr. sc. Marko Rogošić TEACHING WEEKLY SEMESTER Lectures 3 45 Laboratory 1 15 Seminar 1 15 Overall: 75 ECTS: 7.0 PURPOSE: Within
More informationEffects of the Initial Droplet Temperature on the Vaporization Process at High Pressure
Effects of the Initial Droplet Temperature on the ernando. achini fachini@lcp.inpe.br Inst Nac Pesquisas Espaciais - INPE 630-000 Cachoeira Pta,.SP. Brazil Effects of the Initial Droplet Temperature on
More informationMixing and Combustion in Dense Mixtures by William A. Sirignano and Derek Dunn-Rankin
Mixing and Combustion in Dense Mixtures by William A. Sirignano and Derek Dunn-Rankin At very high pressures and densities, what is different and what is similar about the processes of Injection and Atomization,
More informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationA Fully Resolved UWS Droplet Simulation
ILASS-Americas 24th Annual Conference on Liquid Atomization and Spray Systems, San Antonio, TX, May 2012 A Fully Resolved UWS Droplet Simulation D.T.Ryddner, M.F. Trujillo Department of Mechanical Engineering
More informationA Numerically Stable Method for Integration of the Multicomponent Species Diffusion Equations
Journal of Computational Physics 174, 460 472 (2001) doi:101006/jcph20016930, available online at http://wwwidealibrarycom on A Numerically Stable Method for Integration of the Multicomponent Species Diffusion
More informationOutline. Definition and mechanism Theory of diffusion Molecular diffusion in gases Molecular diffusion in liquid Mass transfer
Diffusion 051333 Unit operation in gro-industry III Department of Biotechnology, Faculty of gro-industry Kasetsart University Lecturer: Kittipong Rattanaporn 1 Outline Definition and mechanism Theory of
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationJ. Chem. Thermodynamics
J. Chem. Thermodynamics 43 (2011) 1192 1203 Contents lists available at ScienceDirect J. Chem. Thermodynamics journal homepage: www.elsevier.com/locate/jct Form of multicomponent Fickian diffusion coefficients
More informationRigorous calculation of LNG flow reliefs using the GERG-2004 equation of state
Rigorous calculation of LNG reliefs using the GERG-2004 equation of state Luigi Raimondi Process Simulation Services www.xpsimworld.com Via Galvani 105, 20025 Legnano (MI) - Italy The design of process
More informationModeling Evaporation and Microexplosion of Water-in-Alkane Emulsion Droplets
Modeling Evaporation and Microexplosion of Water-in-Alkane Emulsion Droplets Patrick C. Le Clercq, Berthold Noll and, Manfred Aigner DLR, German Aerospace Center, Stuttgart, Germany Introduction Water
More informationComputational modeling of the effects of support fibers on evaporation of fiber-supported droplets in reduced gravity
Paper # 070HE-0020 Topic: Heterogeneous combustion, sprays & droplets 8 th US National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University
More informationPHASE EQUILIBRIUM OF MULTICOMPONENT MIXTURES: CONTINUOUS MIXTURE GIBBS FREE ENERGY MINIMIZATION AND PHASE RULE
Chemical Engineering Communications Volume 54, Pages139-148, 1987 PHASE EQUILIBRIUM OF MULTICOMPONENT MIXTURES: CONTINUOUS MIXTURE GIBBS FREE ENERGY MINIMIZATION AND PHASE RULE Pin Chan Du and G.Ali Mansoori
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationChemical Derivatives and Elemental Transport Coefficients in Plasma Flows Near Local Equilibrium
hemical Derivatives and Elemental Transport oefficients in Plasma Flows Near Local Equilibrium Alessio rsini and Elena V. Kustova Avio Aerospace Propulsion, via I Maggio 99, 10040 Rivalta di Torino, Italy
More informationWhy d 2 -law does not hold during droplet vaporization in microgravity conditions? Christian Chauveau a, Madjid Birouk b and Iskender Gökalp a*
Why d 2 -law does not hold during droplet vaporization in microgravity conditions? Christian Chauveau a, Madjid Birouk b and Iskender Gökalp a* a ICARE, CNRS, C, Avenue de la Recherche Scientifique, 457
More informationThermophysical Properties of Ethane from Cubic Equations of State
Thermophysical Properties of Ethane from Cubic Equations of State MIHAELA NOUR, DANIELA DUNA, MIRELA IONITA, VIOREL FEROIU *, DAN GEANA Politehnica University Bucharest, Department of Inorganic Chemistry,
More informationSimulation of evaporation and combustion of droplets using a VOF method
Simulation of evaporation and combustion of droplets using a VOF method P. Keller, P.A. Nikrityuk, B. Meyer June 23, 2010 Motivation Mathematics Physics Validation & Test Cases Conclusions Motivation (1)
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations
A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
More informationComparison of distillation arrangement for the recovery process of dimethyl sulfoxide
Korean J. Chem. Eng., 24(3), 438-444 (2007) SHORT COMMUNICATION Comparison of distillation arrangement for the recovery process of dimethyl sulfoxide Jungho Cho and Dong Min Kim* Department of Chemical
More informationThermodynamics I. Properties of Pure Substances
Thermodynamics I Properties of Pure Substances Dr.-Eng. Zayed Al-Hamamre 1 Content Pure substance Phases of a pure substance Phase-change processes of pure substances o Compressed liquid, Saturated liquid,
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 informationTheoretical Developments in Group Combustion of Droplets and Sprays
Theoretical Developments in Group Combustion of Droplets and Sprays William A. Sirignano University of California, Irvine Collaborations: Guang Wu, current student; Randall Imaoka, former student, US Navy;
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 informationRotational-translational relaxation effects in diatomic-gas flows
Rotational-translational relaxation effects in diatomic-gas flows V.V. Riabov Department of Computer Science, Rivier College, Nashua, New Hampshire 03060 USA 1 Introduction The problem of deriving the
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationA Corresponding State Theory for the Viscosity of Liquids Bull. Korean Chem. Soc. 2008, Vol. 29, No Articles
A Corresponding State Theory for the Viscosity of Liquids Bull. Korean Chem. Soc. 2008, Vol. 29, No. 1 33 Articles A Corresponding State Theory for the Viscosity of Liquids Wonsoo Kim * and Sukbae Lee
More informationNON-EQUILIBRIUM THERMODYNAMICS
NON-EQUILIBRIUM THERMODYNAMICS S. R. DE GROOT Professor of Theoretical Physics University of Amsterdam, The Netherlands E MAZUR Professor of Theoretical Physics University of Leiden, The Netherlands DOVER
More informationConvective Heat and Mass Transfer Prof. A.W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay
Convective Heat and Mass Transfer Prof. A.W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay Module No. # 01 Lecture No. # 32 Stefan Flow Model We are now familiar with
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationIntroduction to Mass Transfer
Introduction to Mass Transfer Introduction Three fundamental transfer processes: i) Momentum transfer ii) iii) Heat transfer Mass transfer Mass transfer may occur in a gas mixture, a liquid solution or
More informationEQUATION OF STATE DEVELOPMENT
EQUATION OF STATE DEVELOPMENT I. Nieuwoudt* & M du Rand Institute for Thermal Separation Technology, Department of Chemical Engineering, University of Stellenbosch, Private bag X1, Matieland, 760, South
More informationCOMPUTATIONAL STUDY OF PARTICLE/LIQUID FLOWS IN CURVED/COILED MEMBRANE SYSTEMS
COMPUTATIONAL STUDY OF PARTICLE/LIQUID FLOWS IN CURVED/COILED MEMBRANE SYSTEMS Prashant Tiwari 1, Steven P. Antal 1,2, Michael Z. Podowski 1,2 * 1 Department of Mechanical, Aerospace and Nuclear Engineering,
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationFriction Coefficient Analysis of Multicomponent Solute Transport Through Polymer Membranes
Friction Coefficient Analysis of Multicomponent Solute Transport Through Polymer Membranes NARASIMHAN SUNDARAM NIKOLAOS A. PEPPAS* School of Chemical Engineering, Purdue University, West Lafayette, Indiana
More informationEnhancement of Heat Transfer by an Electric Field for a Drop Translating at Intermediate Reynolds Number
Rajkumar Subramanian M. A. Jog 1 e-mail: milind.jog@uc.edu Department of Mechanical, Industrial, and Nuclear Engineering, University of Cincinnati, Cincinnati, OH 45221-0072 Enhancement of Heat Transfer
More informationof Nebraska - Lincoln
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Combustion Research at University of Nebraska- Lincoln Mechanical & Materials Engineering, Department of March 2006 Raghavan,
More informationBAE 820 Physical Principles of Environmental Systems
BAE 820 Physical Principles of Environmental Systems Estimation of diffusion Coefficient Dr. Zifei Liu Diffusion mass transfer Diffusion mass transfer refers to mass in transit due to a species concentration
More informationCOEFFICIENTS OF VISCOSITY FOR A FLUID IN A MAGNETIC FIELD OR IN A ROTATING SYSTEM
Hooyman, G.J. Physica XXI Mazur, P. " 355-359 Groot, S. R. de 1955 COEFFICIENTS OF VISCOSITY FOR A FLUID IN A MAGNETIC FIELD OR IN A ROTATING SYSTEM by G. J. HOOYMAN, P. MAZUR and S. R. DE GROOT Instituut
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 informationOn the Validity of the Assumption of Local Equilibrium in Non-Equilibrium Thermodynamics
On the Validity of the Assumption of Local Equilibrium in Non-Equilibrium Thermodynamics Arieh Ben-Naim Department of Physical Chemistry The Hebrew University of Jerusalem Givat Ram, Jerusalem 91904 Israel
More informationBoundary Conditions in Fluid Mechanics
Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial
More informationThermodynamics Qualifying Exam Study Material
Thermodynamics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of undergraduate engineering thermodynamics topics. These topics are listed below for clarification.
More informationModelling of methane gas hydrate incipient conditions via translated Trebble-Bishnoi-Salim equation of state
Modelling of methane gas hydrate incipient conditions via translated Trebble-Bishnoi-Salim equation of state Carlos Giraldo and Matthew Clarke Department of Chemical and Petroleum Engineering, the University
More informationCooling Temperatures of Binary Mixed Refrigerants: Vapor-Liquid-Liquid Equilibrium versus Vapor-Liquid Equilibrium
1 Cooling Temperatures of Binary Mixed Refrigerants: Vapor-Liquid-Liquid Equilibrium versus Vapor-Liquid Equilibrium N. Tzabar, H.J.M. ter Brake Energy Materials and Systems Faculty of Science and Technology
More informationAppendix 4. Appendix 4A Heat Capacity of Ideal Gases
Appendix 4 W-143 Appendix 4A Heat Capacity of Ideal Gases We can determine the heat capacity from the energy content of materials as a function of temperature. The simplest material to model is an ideal
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationPART ONE TWO-PHASE FLOW
PART ONE TWO-PHASE FLOW 1 Thermodynamic and Single-Phase Flow Fundamentals 1.1 States of Matter and Phase Diagrams f Pure Substances 1.1.1 Equilibrium States Recall from thermodynamics that f a system
More informationSpray evaporation model sensitivities
Center for Turbulence Research Annual Research Briefs 20 23 Spray evaporation model sensitivities By Shashank, E. Knudsen AND H. Pitsch. Motivation and objective The energy density of solid- and liquid-phase
More informationTheory of Interfacial Tension of Partially Miscible Liquids
Theory of Interfacial Tension of Partially Miscible Liquids M.-E. BOUDH-HIR and G.A. MANSOORI * University of Illinois at Chicago (M/C 063) Chicago, Illinois USA 60607-7052 Abstract The aim of this work
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 informationPressure-coupled vaporization response of n-pentane fuel droplet at subcritical and supercritical conditions
Available online at www.sciencedirect.com Proceedings of the Combustion Institute 33 (211) 1997 23 Proceedings of the Combustion Institute www.elsevier.com/locate/proci Pressure-coupled vaporization response
More informationEXTENDED SMOKER S EQUATION FOR CALCULATING NUMBER OF STAGES IN DISTILLATION
EXTENDED SMOKER S EQUATION FOR CALCULATING NUMBER OF STAGES IN DISTILLATION Santanu Bandyopadhyay Energy Systems Engineering and Department of Mechanical Engineering, Indian Institute of Technology, Bombay,
More informationAccuracy of vapour ^ liquid critical points computed from cubic equations of state
High Temperatures ^ High Pressures 2000 volume 32 pages 449 ^ 459 15 ECTP Proceedings pages 433 ^ 443 DOI:10.1068/htwu303 Accuracy of vapour ^ liquid critical points computed from cubic equations of state
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationADSORPTION IN MICROPOROUS MATERIALS: ANALYTICAL EQUATIONS FOR TYPE I ISOTHERMS AT HIGH PRESSURE
ADSORPTION IN MICROPOROUS MATERIALS: ANALYTICAL EQUATIONS FOR TYPE I ISOTHERMS AT HIGH PRESSURE A. L. MYERS Department of Chemical and Biomolecular Engineering University of Pennsylvania, Philadelphia
More informationThermodynamic Properties of Refrigerant R116 from Cubic Equations of State
Thermodynamic Properties of Refrigerant R116 from Cubic Equations of State DANIELA DUNA, MIRELA IONITA, VIOREL FEROIU *, DAN GEANA Politehnica University Bucharest, Department of Applied Physical Chemistry
More informationAN EXPERIMENTAL STUDY ON THE DROPLET VAPORIZATION: EFFECTS OF HEAT CONDUCTION THROUGH THE SUPPORT FIBER.
Paper ID ILASS8-4- ILASS8-A2 ILASS 28 Sep. 8-, 28, Como Lake, Italy AN EXPERIMENTAL STUDY ON THE DROPLET VAPORIZATION: EFFECTS OF HEAT CONDUCTION THROUGH THE SUPPORT FIBER. Chauveau, C.*, Halter, F., Lalonde,
More informationThermal Systems. What and How? Physical Mechanisms and Rate Equations Conservation of Energy Requirement Control Volume Surface Energy Balance
Introduction to Heat Transfer What and How? Physical Mechanisms and Rate Equations Conservation of Energy Requirement Control Volume Surface Energy Balance Thermal Resistance Thermal Capacitance Thermal
More informationChemical Potential. Combining the First and Second Laws for a closed system, Considering (extensive properties)
Chemical Potential Combining the First and Second Laws for a closed system, Considering (extensive properties) du = TdS pdv Hence For an open system, that is, one that can gain or lose mass, U will also
More informationSubcritical and Supercritical Nanodroplet Evaporation: A Molecular Dynamics Investigation
ILASS Americas 2th Annual Conference on Liquid Atomization and Spray Systems, Chicago, IL, May 27 Subcritical and Supercritical Nanodroplet Evaporation: A Molecular Dynamics Investigation E. S. Landry,
More informationA NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER
SCIENTIFIC RESEARCH AND EDUCATION IN THE AIR FORCE-AFASES 2016 A NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER Alexandru DUMITRACHE*, Florin FRUNZULICA ** *Institute of
More informationOCN 623: Thermodynamic Laws & Gibbs Free Energy. or how to predict chemical reactions without doing experiments
OCN 623: Thermodynamic Laws & Gibbs Free Energy or how to predict chemical reactions without doing experiments Definitions Extensive properties Depend on the amount of material e.g. # of moles, mass or
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [Ecole Centrale Paris] On: 27 January 29 Access details: Access Details: [subscription number 79383386] Publisher Taylor & Francis Informa Ltd Registered in England and
More informationDroplet evaporation: A molecular dynamics investigation
JOURNAL OF APPLIED PHYSICS 102, 124301 2007 Droplet evaporation: A molecular dynamics investigation E. S. Landry, S. Mikkilineni, M. Paharia, and A. J. H. McGaughey a Department of Mechanical Engineering,
More informationResearch Article. Slip flow and heat transfer through a rarefied nitrogen gas between two coaxial cylinders
Available online wwwjocprcom Journal of Chemical and Pharmaceutical Research, 216, 8(8):495-51 Research Article ISSN : 975-7384 CODEN(USA) : JCPRC5 Slip flow and heat transfer through a rarefied nitrogen
More informationNumerical Study of Laminar Annular Two-Phase Flow in Effervescent Atomizers
ILASS Americas 28th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2016 Numerical Study of Laminar Annular Two-Phase Flow in Effervescent Atomizers C.K. Mohapatra and M.A.
More informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
More information"Energy Applications: Impact of Data and Models"
"Energy Applications: Impact of Data and Models" Energy Applications refers in this particular case to the wide application of equations of state upstream in the Production of Oil and Gas. The petroleum
More informationheat transfer process where a liquid undergoes a phase change into a vapor (gas)
Two-Phase: Overview Two-Phase two-phase heat transfer describes phenomena where a change of phase (liquid/gas) occurs during and/or due to the heat transfer process two-phase heat transfer generally considers
More informationADVANCES IN THERMODIFFUSION AND THERMOPHORESIS (SORET EFFECT) IN LIQUID MIXTURES
Frontiers in Heat and Mass Transfer Available at www.thermalfluidscentral.org ADVANCES IN THERMODIFFUSION AND THERMOPHORESIS (SORET EFFECT) IN LIQUID MIXTURES Morteza Eslamian* School of Engineering and
More informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
More informationConvective Burning of Gaseous Fuel Pockets and Supercritical Droplets
Convective Burning of Gaseous Fuel Pockets and Supercritical Droplets J. DAOU and B. ROGG* Lehrstuhl für Strömungsmechanik, Institut für Thermo- und Fluiddynamik, Ruhr-Universität, D-44780 Bochum, Germany
More informationEquations of State. Equations of State (EoS)
Equations of State (EoS) Equations of State From molecular considerations, identify which intermolecular interactions are significant (including estimating relative strengths of dipole moments, polarizability,
More informationScienceDirect. Modelling CO 2 Water Thermodynamics Using SPUNG Equation of State (EoS) concept with Various Reference Fluids
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 51 (2014 ) 353 362 7th Trondheim CCS Conference, TCCS-7, June 5-6 2013, Trondheim, Norway Modelling CO 2 Water Thermodynamics Using
More informationPart I.
Part I bblee@unimp . Introduction to Mass Transfer and Diffusion 2. Molecular Diffusion in Gasses 3. Molecular Diffusion in Liquids Part I 4. Molecular Diffusion in Biological Solutions and Gels 5. Molecular
More informationNon equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi
Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale Miguel Rubi References S.R. de Groot and P. Mazur, Non equilibrium Thermodynamics, Dover, New York, 1984 J.M. Vilar and
More information