Instability Analysis of the Separated Boundary Layer in Shock Tubes
|
|
- Lisa Shepherd
- 5 years ago
- Views:
Transcription
1 AIAA SciTech Forum 9 - January 7, Grapevine, Texas 55th AIAA Aerospace Sciences Meeting.54/ Instability Analysis of the Separated Boundary Layer in Shock Tubes Kevin Grogan and Matthias Ihme Department of Mechanical Engineering, Stanford University, Stanford, CA 945 The instabilities manifesting from the shock-separated boundary layer in shock tube reaction-kinetic experiments are examined using linear stability analysis. The shockseparated boundary layer is modeled using a polynomial profile with matching conditions across the shock, and the compressible Rayleigh equation is employed to determine the most hydrodynamically unstable wavelengths and growth rates. An instability time scale is identified to characterize the mixing of the boundary layer fluid with the test gas using the length scales and time scales determined from the stability analysis. From this analysis, the shock-separated boundary layer is found to be unstable for all Mach numbers examined. Furthermore, an argon-diluted mixture is found to yield a more stable boundary layer when compared to a nitrogen-diluted mixture. From this, it is concluded that the dilution of the test gas with argon would yield less mixing of the boundary layer fluid with the core flow and therefore, would reduce the likelihood of a non-ideal ignition in a shock tube reaction-kinetic experiment for the low-temperature regime analyzed. Introduction Shock tubes are a primary tool for chemical kinetic measurements, providing ignition delay times, species time-histories, and elementary reaction rates. They replicate zero-dimensional reactors to high accuracy for test times up to the millisecond-range allowing for idealized modeling of the combustion process at well-defined temperature and pressure conditions. A shock tube consists of a driver and a driven section, and an x-t diagram illustrating the relevant regions is shown in Fig.. A shock tube experiment is initiated by the rupture of a diaphragm separating the driver and driven sections, which are denoted as 4 and, respectively, in Fig.. The incident shock propagates towards the end wall with the shocked gas in region immediately following. The incident shock reflects off the end wall and intersects the shocked gas in region to produce the twice-shocked test gas in region 5; the test gas in region 5 is ideally stagnant due to the no-penetration boundary condition at the wall. A contact surface, separating the driver and driven gas and traveling at the velocity of the fluid, trails the incident shock and interacts with the reflected shock wave. Figure illustrates a condition where the driver and driven gas are the same; however, experiments with a tailored driver gas extend the test time, τ TEST, from that shown by modifying how the contact surface interacts with the reflected shock. The expansion fan initially propagates into region 4 as a procession of Mach waves traveling at the local speed of sound; the expansion fan reflects off the driver section s wall in a complex fashion with interaction between the reflected part of the fan and that still propagating towards the wall. The ideal combustion mode for shock tube experiments is strong ignition, which is characterized by the uniform ignition of the test gas beginning at the end wall. However, inhomogeneous or weak ignition has been found in shock tube systems under certain conditions; 4 8 weak ignition can occur when pockets of reactive gas ignite prior to the bulk gas, and this is primarily due to a confluence of chemical sensitivities and gasdynamic fluctuations. There has been increased interest in recent years for chemical kinetics measurements in shock tubes within the negative temperature coefficient (NTC) region due to advancements in extending shock tube test times., 9 In the NTC regime, the ignition time decreases with temperature, which could lead to the preignition of the cool boundary layer fluid near the shock tube walls; furthermore, the interaction of Graduate Student, AIAA member Assistant Professor, AIAA member of Copyright 7 by the, American Inc. All Institute rights reserved. of Aeronautics and Astronautics
2 7 6 t 5 TEST 4 x l DRIVER l DRIVEN Figure. Diagram illustrating a typical shock tube experiment. Regions are labeled in accordance with standard convention. Shock waves are shown by solid RED lines, the contact surface is depicted by a dashed GREEN line, and the expansion fan is represented by thin BLUE lines Figure. Detailed simulations of SBLI in a shock tube. Snapshots of the temperature field are taken subsequent to shock-reflection from the end wall. The nominal conditions behind the reflected shock are constant while the incident shock Mach number, M s, is varied. A 78.% nitrogen-diluted n-heptane/oxygen mixture is simulated. The specific heat ratio is γ =.. the reflected shock wave with the boundary layer could yield enhanced mixing of the boundary layer fluid with the core gas leading to weak ignition. Shock boundary layer interaction (SBLI) demonstrates a regime-dependence based on the incident shock Mach number and the specific heat ratio; canonical SBLI regimes include incipient separation and shock bifurcation. Previous detailed simulations of SBLI in shock tube experiments are shown in Fig.. As shown in this figure, the same nominal test conditions for a shock tube experiment (i.e., T 5, p 5, and composition) can produce drastically differing flow fields when the incident shock Mach number, M s, is varied. Of interest to this work is the intermediate case (M s =.47) between the quiescent simulation (M s =.5) and the simulation showing bifurcation (M s =.79). This intermediate case shows the separation of the incident boundary layer by the reflected shock wave and the ejection of boundary layer fluid into the bulk gas. Hence, it is proposed that the instabilities resultant from the SBLI can produce of
3 u, + U s x s 5 U s y p, p 5, x 5 Figure. Illustration showing the development of instabilities in a separated shear layer immediately subsequent to the reflected shock. Velocities are given in the frame of the reflected shock. inhomogeneities in temperature that could lead to non-ideal ignition phenomenon for mixtures in the NTC regime. Therefore, it is the objective of this paper to analyze the intermediate SBLI shown in Fig. and to develop a relevant time scale for comparison with shock tube experiments. From this, insight will be gleaned about a source of inhomogeneity that could give rise to non-ideal ignition in shock tubes. This paper proceeds with the mathematical model in Sec. I, where the theory and numerical method are described; subsequently, Sec. II includes the results from the application of the model; finally, a discussion of the conclusions drawn from this work is contained in Sec. III. I. Mathematical Model When the reflected shock wave is sufficiently strong, the flow reverses direction in the laboratory frame, and a separated shear layer forms. An illustration of this process is shown in Fig.. The separated shear layer is proposed to engender instabilities. After a certain time scale related to the growth of the linear instability, instabilities in the shear layer become non-linear, and vortices are formed. Furthermore, the instabilities formed by the SBLI could transition the flow into turbulence. Hence, a characteristic time scale for the ejection of nominal vortices into the test gas, τ INST, is developed using the length scales and time scales of the linear instability as well as physical dimensions of the boundary layer. The stability of the separated shear layer is inferred from the eigenvalues of the compressible Rayleigh equation: ˆp + [ T 5 ũ 5 T 5 (Ũs + ũ 5 ) c ] ˆp α T 5 { T 5 M [(Ũs + ũ 5 ) c] }ˆp =, () where for an arbitrary variable, χ, the tilde, χ, indicates a non-dimensional quantity (excluding α, c, and ω); the hat, ˆχ, is used to denote a quantity for which the normal modes hypothesis has been applied (i.e., χ = ˆχ(ỹ) exp ( iα x iω t ) ); and the prime, χ, indicates a derivative in y. The pressure fluctuation is normalized by the pressure in the free stream: p = p/p 5, ; the coordinates, x, y, are non-dimensionalized by the post-shock boundary layer thickness (e.g., x = x/δ 5 ); the velocities (i.e., Ũs, and ũ 5 ) are normalized by the speed of the reflected shock, U s ; the temperature is normalized with respect to the free stream value, T 5 = T 5 /T 5, ; M = U s /a 5, is the free stream Mach number; and the phase speed is related to the growth rate and the wave number by c = ω/α C. In the following, a temporal stability analysis is undergone using the wave speed, c, as the eigenvalue. From the application of Eq., it is assumed that the flow is parallel and that the instability mechanism for this flow is primarily inviscid. A cubic profile is assumed for the incident momentum and thermal boundary layers: 4 u = T = [ (y/δ ) (y/δ ) ] u,, for y < δ (a) u,, for y δ [ T + (T, T ) (y/δ ) (y/δ ) ], for y < δ (b) T,, for y δ, where the thicknesses of the boundary layers have been presumed to be similar due to the Prantl number of
4 of gases analyzed being close to unity. Additionally, the wall is assumed to be isothermal and at a value equal to the initial gas temperature (i.e., K). Other functions may be used to describe the intercepted boundary layer such as the similarity solution of Mirels; 5 however, a simplistic profile is thought to be sufficient to capture the shape of the separated profile while making the exploration of the parameter space more computationally efficient. A quartic profile is assumed for the separated boundary layer: 4 i= u 5 = a i (y/δ 5 ) i, for y < δ 5 (a), for y δ 5 [ T + (T 5, T ) T 5 = (y/δ 5) (y/δ 5) ], for y < δ 5 (b) T 5,, for y δ, wherein the thermal boundary layer is proposed to adjust with the momentum boundary layer. As before, the wall is assumed to be isothermal. The quartic equation has six parameters that need to be accounted for (i.e., a i for i =,..., 4, and δ 5 ). Hence, the six proposed conditions on Eq. are. No slip condition on wall (u 5 y= = ).. Inflection point on the wall (d u 5 /dy y= = ).. Matching condition in free stream (u 5 y=δ5 = ). 4. Stress free condition in free stream (du 5 /dy y=δ5 = ). 5. Conservation of mass across the reflected shock. 6. Balance of momentum across the reflected shock. Mass conservation and momentum balance across the shock yield δ δ ρ (u + U s )dy + ρ, (u, + U s )(δ 5 δ ) = ρ (u + U s ) dy + ρ, (u, + U s ) (δ 5 δ ) + p δ 5 = δ5 δ5 ρ 5 (u 5 + U s )dy, ρ 5 (u 5 + U s ) dy + p 5 δ 5, where the free-stream component of the mass conservation and momentum balance cancel through the Rankine-Hugoniot relations, and the post-shock boundary layer is assumed to thicken. The density profile is found via the assumed thermal boundary layer and the ideal gas equation of state. Since the model utilizes invicid, parallel flows, the pressure is constant across the shear layer. Also, the shear stress at the wall is neglected over the interaction region since its sign necessarily changes for separated flows, which likely makes the integral value of the shear stress on the wall small. Using Eq., the cubic eigenvalue problem is rewritten in the following form: (4a) (4b) where the operators, L k, are given by c L ˆp + c L ˆp + cl ˆp + L ˆp =, (5) [ L = α M (Ũs + ũ 5 ) α T5 (Ũs + ũ 5 ) (Ũs + ũ 5 ) + T ] 5 ũ T 5 T d 5 5 dỹ + (Ũs + ũ 5 ) d dỹ, L = α M (Ũs + ũ 5 ) + α T5 T 5 d T 5 T 5 dỹ d dỹ, L = α M (Ũs + ũ 5 ) T 5, (6c) L = α M T 5. (6a) (6b) (6d) 4 of
5 Equation 5 is discretized using a second-order, centered scheme on a non-uniform grid (i.e., ˆp p and L k L k, where p C N, L k R N N, and N is the number of grid points). The non-uniform grid points are distributed using an arctangent profile, and points are found to be sufficient for convergence in the eigenvalue of the most unstable mode. Finally, the cubic eigenvalue problem is solved by rewriting Eq. 5 as a generalized eigenvalue problem, Aπ = cbπ, where in block-matrix form, A = L L L I I, B = L I I, π = p cp. (7) c p The identity and zero matrices in Eq. 7 are given by I and, respectively. The solution of the discrete eigenvalue problem yields the eigenvalue, c, and the eigenvector, π. Additionally, a Neumann boundary condition is imposed on the pressure fluctuation in the free stream and at the wall. The effect of the instability of the separated boundary layer is the introduction of inhomogeneity into the test gas. A quantification of the inhomogeneity added to the flow within a given timespan can be made using the length scales and time scales obtained from the solution of Eq.. A metric for the inhomogeneity in the flow is postulated to be V BL /V, where V BL is the volume of the fluid ejected into the test gas from the boundary layer, and V is the total volume of fluid behind the reflected shock at a given time; the total volume is taken to be V = U s τ INST πd /4 where τ INST refers to the characteristic time scale of the instability, and D is the diameter of the shock tube. Thus, the volume of fluid ejected by the boundary layer is given by V BL = N VORT j= V VORT,j, (8) where the volume of the j th nominal ejected vortex is denoted by V VORT,j. The volume of this nominal vortex is assumed to have a toroidal shape: V VORT,j = πd(πl INST,j /4), where l INST,j is the characteristic diameter of the j th nominal ejected vortex. The characteristic length of the instability taken to be twice the wave number at the most unstable mode: l INST,j = πδ 5,j /α max. Hence, the volume of the fluid ejected into the test gas is approximated by V VORT,j = π4 Dδ5,j αmax. (9) Similarly, the frequency of the ejection of nominal vortices into the test gas is assumed to be related to the linear growth rate by a constant: f INST,j = b f ω i,max U s /(πδ 5,j ), where b f is a proportionality constant. Hence, the period of the ejected vortex is given by τ j = f INST,j. Equation 8 can now be rewritten as V BL = N VORT j= τinst V VORT,j f INST,j τ j, V VORT f INST dτ, where the approximation of the sum by an integral should hold for large N VORT (i.e., the regime of interest with this model). The value of δ 5 is known from the separated shear layer model and Mirel s similarity solution, where the boundary layer thickness has the functional form ν w (U s + U s )τ δ = (U s u, ) g(m s, γ, Pr), () () where the Prandtl number is denoted by Pr, and the function g is found from the similarity solution of Mirels. 5 By positing that there exists a critical threshold for the ratio of the injected boundary layer fluid to test gas, (V BL /V ) crit, where the quality of ignition changes (i.e., strong to weak), the time scale of the instability can be obtained from Eq. 8 as ( ) α τ INST = β ( max δ ω i,max δ 5 ) [ ] ( D D δ (τ = D/U s ) U s ), () 5 of
6 .5 γ =. (n-heptane/o/n) γ =.47 (n-heptane/o/ar) M s =..5 γ =. (n-heptane/o/n) γ =.47 (n-heptane/o/ar) y/δ.5 y/δ.5 M s =..5 M s =..5 M s = u 5 /U s T 5 /T 5, Figure 4. Comparison of the modeled separated boundary layer profiles for a nitrogen-diluted and an argondiluted n-heptane/o mixture. Both mixtures are diluted at 78.% and have an equivalence ratio of φ =.5. Normalized velocities are given in the laboratory frame. where the non-dimensional constant is given by β = (V BL /V ) crit /(4π b f ), and both α max and ω i,max are found from the solution of Eq.. It is supposed that β is mostly a fuel-dependent parameter that accounts for the sensitivity of the mixture to fluctuations. II. Results The modeled profiles of the separated shear layers are shown in Fig 4. The boundary layer profiles of two n-heptane mixtures are shown for comparison due to the fuel s prominent NTC characteristics. As shown in the figure, the monatomic diluent, argon, exhibits a much thinner profile relative to the incident boundary layer. This is likely due to monatomic gasses being stiffer than diatomic gasses and therefore, less compliant to pressure gradients. Additionally, the variation of the profiles with respect to Mach number is shown to decrease at higher Mach numbers. The most unstable wavelengths and growth rates are shown in Fig. 5; these are found by applying Eq. to the model boundary layers shown in Fig. 4. For all the incident Mach numbers examined, the boundary layer is found to be unstable. This indicates that there are only regimes in which the separated boundary layer is quasi-stable; that is, the manifestation of the instability is much slower than the chemistry. A jump in the most unstable wavelength, α max, is shown for both mixtures between an incident shock Mach number of.5 to.7. This occurs because the growth rate of the lower wavelength mode becomes negligible at higher Mach numbers. Interestingly, the comparative size of the maximum growth rate, ω i,max, between the two gases switches as the most unstable mode changes. The growth rates of the most unstable mode as a function of the wavelength is shown in Fig. 6. At the lower Mach number, two local maxima are demonstrated, and at the higher Mach number, a single global maximum exists. The two maxima at the lower Mach number correspond to distinct modes, and a single mode is tracked at the higher Mach number until approximately α = 4.. Figure 7 shows the eigenfunctions of the separated boundary layer at the most unstable mode for two different Mach numbers. While the pressure fluctuation, ˆp, is found from the solution of Eq., the normal velocity fluctuation, ˆv, is determined using the y-momentum equation: T 5 ˆp ˆv = () iα[(ũs + ũ 5 ) c]γm As shown in Fig. 7, the eigenfunctions demonstrate the highest variability within the boundary layer (i.e., y/δ 5 < ). The normal velocity fluctuations are shown to decay in the free stream and at the wall for all modes; due to the Neumann boundary condition applied to Eq. 7, the derivative of the pressure fluctuation is zero at the wall and in the free stream. Also, the eigenfunctions at the higher Mach number 6 of
7 .5.5 ω i,max α max M s Figure 5. Comparison of the most unstable wavenumber, α max, and the corresponding growth rate, ω i,max. RED lines denote an argon-diluted n-heptane mixture, and BLUE lines denote a nitrogen-diluted n-heptane mixture. Both mixtures are diluted at 78.% and have an equivalence ratio of φ =.5. ωi,max M s =. M s =.5 M s =.8 M s =. 4 5 α (a) N -Dilution, γ =. ωi,max M s =. M s =.5 M s =.8 M s =. 4 5 α (b) Ar-Dilution, γ =.47 Figure 6. Comparison of the growth rates of the most unstable modes for two selected Mach numbers. Both mixtures are diluted at 78.% and have an equivalence ratio of φ =.5. and, correspondingly, higher wavelength are shown to have increased variation. For both modes and gases, the normal velocity fluctuation demonstrates a peak near the outer edge of the separated shear layer; hence, these mode shapes could give rise to vortex formation at this location and entrainment of boundary layer fluid. A comparison of the instability time scale with the ignition delay of n-heptane is shown in Fig. 8. The parameters in Eq. are set to be b f = and (V BL /V ) CRIT =.. The initial temperature is presumed to be T = K. The ignition delay time is taken from an adiabatic, isobaric (HP) reactor using the time at which the concentration of the OH radical has peaked. The chemistry is modeled with a 94-species reduced n-heptane mechanism derived from the detailed mechanism due to Mehl et al. 6 using the procedure of Stagni et al. 7 It is shown in Fig. 8 that the chemical time scale, τ IGN, is much larger than the instability time scale at low temperatures (i.e, T 5 < K for argon dilution and T 5 < 8 K for nitrogen dilution); this indicates that at low temperatures, inhomogeneities emanating from the separated boundary layer may become prevalent and lead to weak ignition. The dimensional time scale indicates that the argon mixture, while less stable than the nitrogen mixture, 7 of
8 .5 M s =. M s =..5 M s =. M s =. y/δ5.5 y/δ y/δ ˆv (a) N -Dilution, γ =. M s =. M s = ˆp (c) N -Dilution, γ =. y/δ ˆv (b) Ar-Dilution, γ =.47 M s =. M s = ˆp (d) Ar-Dilution, γ =.47 Figure 7. Comparison of the eigenfunctions at the most unstable wavelength for two diluted n-heptane mixtures. Both mixtures are diluted at 78.% and have an equivalence ratio of φ =.5. would likely induce less mixing of the cool boundary layer fluid with the bulk gas. Since the dynamic viscosities of the argon-diluted and nitrogen-diluted mixtures are similar, it is concluded that the discrepancy between the two time scales is due to the differing specific heat ratios of the mixtures. Disparate specific heat ratios can yield different shock velocities, kinematic viscosities, and sound speeds, which as shown by Eq., can cause the observed discrepancy in the instability time scale. The lower wavelength mode is shown to produce a smaller instability timescale at low temperatures, which indicates that more boundary layer fluid may become entrained within the test gas in a shorter period of time. Hence, the instability mode may change for low temperature shock tube experiments in a fashion that has a higher propensity for a non-ideal ignition. Also, preheating the test gas would yield a leftward shift in the instability timescale curves due to a lower requisite T 5 /T temperature ratio and correspondingly, a lower Mach number. Hence, the model implies that preheating the test gas would reduce the entrainment of boundary layer fluid and therefore, would lead to a decreased likelihood of weak ignition. III. Conclusions A low-order model is developed to describe the temperature and velocity profiles of a shock-separated boundary layer, and the hydrodynamic stability of the model boundary layer is analyzed using the compressible Rayleigh equation. The profile is found to be unstable for all incident Mach numbers examined. 8 of
9 τ IGN τ INST, Ar-diluted τ INST, N -diluted τ [s] /T 5 [K ] Figure 8. Comparison of the instability time scale with that of the ignition delay curve of n-heptane using a HP reactor model with a test pressure of p 5 =.7 bar. The mixtures are diluted at 78.% and have an equivalence ratio of φ =.5. Additionally, this boundary layer instability is proposed to cause inhomogeneities in the core gas, which could potentially yield weak ignition. An instability time scale is developed for comparison to the ignition time of a test mixture. It is found that while an argon-diluted mixture is more unstable due to a thinner shear layer profile, a nitrogen-diluted mixture has increased mixing of the boundary layer fluid with the core gas. Hence, our analysis suggests that argon dilution is preferable to combat inhomogeneous ignition in shock tube experiments. Acknowledgments The authors gratefully acknowledge financial support through the Air Force Office of Scientific Research under Award No. FA with Dr. Chiping Li as program manager. References R. K. Hanson and D. F. Davidson. Recent advances in laser absorption and shock tube methods for studies of combustion chemistry. Prog. Energy Combust. Sci., 44: 4, 4. O. Trass and D. Mackay. Contact surface tailoring in a chemical shock tube. AIAA J., (9):6 6, 96. G. Ben-Dor, O. Igra, and T. Elpherin, editors. Handbook of Shock Waves. Elsevier, New York,. 4 V. V. Voevodsky and R. I. Soloukhin. On the mechanism and explosion limits of hydrogen-oxygen chain self-ignition in shock waves. Proc. Combust. Inst., :79 8, J. W. Meyer and A. K. Oppenheim. On the shock-induced ignition of explosive gases. Proc. Combust. Inst., :5 64, D. J. Vermeer, J. W. Meyer, and A. K. Oppenheim. Auto-ignition of hydrocarbons behind reflected shock waves. Combust. Flame, 8:7 6, H. Yamashita, J. Kasahara, Y. Sugiyama, and A. Matsuo. Visualization study of ignition modes behind bifurcatedreflected shock waves. Combust. Flame, 59: ,. 8 K. P. Grogan and M. Ihme. Weak and strong ignition of hydrogen/oxygen mixtures in shock-tube systems. Proc. Combust. Inst., 5():8 89, 5. 9 M. F. Campbell, S. Wang, C.S. Goldenstein, R. M. Spearrin, A. M. Tulgestke, L. T. Zaczek, D. F. Davidson, and R. K. Hanson. Constrained reaction volume shock tube study of n-heptane oxidation: Ignition delay times and time-histories of multiple species and temperature. Proc. Combust. Inst., 5: 9, 5. K. P. Grogan and M. Ihme. Regimes describing shock boundary layer interaction and ignition in shock tubes. Proc. Combust. Inst., 6, 7. In Press. D. R. Chapman, D. M. Kuehn, and H. K. Larson. Investigation of separated flows in supersonic and subsonic streams with emphasis on the effect of transition. NACA-TN-56, Ames Aeronautical Lab, Moffett Field, CA, of
10 H. Mark. The interaction of a reflected shock wave with the boundary layer in a shock tube. J. Aeronaut. Sci., 4:4 6, 957. W. Criminale, T. Jackson, and R. Joslin. Theory and Computation of Hydrodynamic Stability. Cambridge University Press, Cambridge,. 4 W. Kays, M. Crawford, and B. Weigand. Convective Heat and Mass Transfer. McGraw-Hill, 5. 5 H. Mirels. Laminar boundary layer behind shock advancing into stationary fluid. NACA TN 4, Lewis Flight Propulsion Laboratory, Cleveland, OH, M. Mehl, W. J. Pitz, C. K. Westbrook, and H. J. Curran. Kinetic modeling of gasoline surrogate components and mixtures under engine conditions. Proc. Combust. Inst., :9,. 7 A. Stagni, A. Frassoldati, A. Cuoci, T. Faravelli, and E. Ranzi. Skeletal mechanism reduction through species-targeted sensitivity analysis. Combust. Flame, 6:8 9, 6. of
Mild Ignition Phenomena in Rapid Compression Machines
25 th ICDERS August 2 7, 2015 Leeds, UK Kevin P. Grogan a, S. Scott Goldsborough b, Matthias Ihme a a Stanford University, Stanford, CA 94305 b Argonne National Laboratory, Argonne, IL 60439 1 Introduction
More informationIgnition delay-time study of fuel-rich CH 4 /air and CH 4 /additive/air mixtures over a wide temperature range at high pressure
25 th ICDERS August 2 7, 2015 Leeds, UK Ignition delay-time study of fuel-rich CH 4 /air and CH 4 /additive/air mixtures over a wide temperature range at high pressure Jürgen Herzler, Mustapha Fikri, Oliver
More informationEffects of Variation of the Flame Area and Natural Damping on Primary Acoustic Instability of Downward Propagating Flames in a Tube
5 th ICDERS August 7, 015 Leeds, UK Effects of Variation of the Flame Area and Natural Damping on Primary Acoustic Instability of Downward Propagating Flames in a Tube Sung Hwan Yoon and Osamu Fujita Division
More informationTemperature time-history measurements in a shock tube using diode laser absorption of CO 2 near 2.7 µm
23 rd ICDERS July 24-29, 2011 Irvine, USA Temperature time-history measurements in a shock tube using diode laser absorption of CO 2 near 2.7 µm Wei Ren, Sijie Li, David F Davidson, and Ronald K Hanson
More informationEffect of orientation on the ignition of stoichiometric ethylene mixtures by stationary hot surfaces
Josue Melguizo-Gavilanes and Joseph E.Shepherd. Effect of orientation on the ignition of stoichiometric ethylene mixtures by stationary hot surfaces. Paper No. 981, 26th International Colloquium on the
More informationPostdoctoral Scholar Hanson Research Group, Mechanical Engineering, Stanford University, CA
SHENGKAI WANG EDUCATION Phone: +1 (650) 391-6853 Email: sk.wang@stanford.edu Address: 418 Panama Mall, Rm. 106 Stanford, CA, 94305 Ph.D. in Mechanical Engineering, Stanford University, Stanford, CA 01/2017
More informationAA210A Fundamentals of Compressible Flow. Chapter 13 - Unsteady Waves in Compressible Flow
AA210A Fundamentals of Compressible Flow Chapter 13 - Unsteady Waves in Compressible Flow The Shock Tube - Wave Diagram 13.1 Equations for irrotational, homentropic, unsteady flow ρ t + x k ρ U i t (
More informationS. Kadowaki, S.H. Kim AND H. Pitsch. 1. Motivation and objectives
Center for Turbulence Research Annual Research Briefs 2005 325 The dynamics of premixed flames propagating in non-uniform velocity fields: Assessment of the significance of intrinsic instabilities in turbulent
More informationHierarchical approach
Chemical mechanisms Examine (i) ways in which mechanisms are constructed, (ii)their dependence on rate and thermodynamic data and (iii) their evaluation using experimental targets Copyright 2011 by Michael
More informationLaminar Premixed Flames: Flame Structure
Laminar Premixed Flames: Flame Structure Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Course Overview Part I: Fundamentals and Laminar Flames Introduction Fundamentals and mass balances of
More informationANALYSIS OF THE PREDICTION ABILITY OF REACTION MECHANISMS FOR CFD MODELING OF THE COMBUSTION IN HIGH VELOCITY ENGINES
ANALYSIS OF THE PREDICTION ABILITY OF REACTION MECHANISMS FOR CFD MODELING OF THE COMBUSTION IN HIGH VELOCITY ENGINES V. Kopchenov*, S. Batura*, L. Bezgin*, N. Titova*, A. Starik* *CIAM, Russia Keywords:
More informationHigh-pressure shock-tube study of the ignition of fuel-rich CH 4 /air and CH 4 /additive/air mixtures over a wide temperature range
High-pressure shock-tube study of the ignition of fuel-rich CH 4 /air and CH 4 /additive/air mixtures over a wide temperature range J. Herzler, M. Fikri, O. Welz, C. Schulz Institute for Combustion and
More informationAsymptotic Structure of Rich Methane-Air Flames
Asymptotic Structure of Rich Methane-Air Flames K. SESHADRI* Center for Energy and Combustion Research, Department of Mechanical and Aerospace Engineering, University of California at San Diego, La Jolla,
More informationDriver-gas Tailoring For Test-time Extension Using Unconventional Driver Mixtures
University of Central Florida Electronic Theses and Dissertations Masters Thesis (Open Access) Driver-gas Tailoring For Test-time Extension Using Unconventional Driver Mixtures 006 Anthony Amadio University
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationThe evolution and cellular structure of a detonation subsequent to a head-on interaction with a shock wave
Combustion and Flame 151 (2007) 573 580 www.elsevier.com/locate/combustflame The evolution and cellular structure of a detonation subsequent to a head-on interaction with a shock wave Barbara B. Botros
More informationFluid Mechanics. Chapter 9 Surface Resistance. Dr. Amer Khalil Ababneh
Fluid Mechanics Chapter 9 Surface Resistance Dr. Amer Khalil Ababneh Wind tunnel used for testing flow over models. Introduction Resistances exerted by surfaces are a result of viscous stresses which create
More informationTHERMODYNAMIC ANALYSIS OF COMBUSTION PROCESSES FOR PROPULSION SYSTEMS
2nd AIAA Aerospace Sciences Paper 2-33 Meeting and Exhibit January -8, 2, Reno, NV THERMODYNAMIC ANALYSIS OF COMBUSTION PROCESSES FOR PROPULSION SYSTEMS E. Wintenberger and J. E. Shepherd Graduate Aeronautical
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 informationAnalysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data
45th AIAA Aerospace Science Meeting and Exhibit, January 8 11, 25/Reno, Nevada Analysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data M. Wu and M.P. Martin Mechanical
More informationModeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet Model with Heat-Loss Effects
9 th U. S. National Combustion Meeting Organized by the Central States Section of the Combustion Institute May 17-20, 2015 Cincinnati, Ohio Modeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet
More informationLecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory
Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory 8.-1 Systems, where fuel and oxidizer enter separately into the combustion chamber. Mixing takes place by convection and diffusion. Only where
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationNUMERICAL INVESTIGATION OF IGNITION DELAY TIMES IN A PSR OF GASOLINE FUEL
NUMERICAL INVESTIGATION OF IGNITION DELAY TIMES IN A PSR OF GASOLINE FUEL F. S. Marra*, L. Acampora**, E. Martelli*** marra@irc.cnr.it *Istituto di Ricerche sulla Combustione CNR, Napoli, ITALY *Università
More informationSupersonic Combustion Simulation of Cavity-Stabilized Hydrocarbon Flames using Ethylene Reduced Kinetic Mechanism
42nd AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit 9-12 July 2006, Sacramento, California AIAA 2006-5092 Supersonic Combustion Simulation of Cavity-Stabilized Hydrocarbon Flames using Ethylene
More informationNumerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure
21 st ICDERS July 23-27, 27 Poitiers, France Numerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure Ahmad S. El Sayed, Cécile B. Devaud Department of Mechanical
More informationThe role of diffusion at shear layers in irregular detonations
The role of diffusion at shear layers in irregular detonations Marco Arienti 1 Joseph E. Shepherd 2 1 United Technologies Research Center, 411 Silver Lane, East Hartford, CT 06108 2 California Institute
More information1. Introduction Some Basic Concepts
1. Introduction Some Basic Concepts 1.What is a fluid? A substance that will go on deforming in the presence of a deforming force, however small 2. What Properties Do Fluids Have? Density ( ) Pressure
More informationTURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows
TURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows Presented by William A. Sirignano Mechanical and Aerospace Engineering University of California
More informationAME 513. " Lecture 8 Premixed flames I: Propagation rates
AME 53 Principles of Combustion " Lecture 8 Premixed flames I: Propagation rates Outline" Rankine-Hugoniot relations Hugoniot curves Rayleigh lines Families of solutions Detonations Chapman-Jouget Others
More informationImpact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions
Journal of Physics: Conference Series PAPER OPEN ACCESS Impact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions To cite this article: A Rosiak and
More informationWell Stirred Reactor Stabilization of flames
Well Stirred Reactor Stabilization of flames Well Stirred Reactor (see books on Combustion ) Stabilization of flames in high speed flows (see books on Combustion ) Stabilization of flames Although the
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 informationApplication of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures
Supplemental material for paper published in the International J of Hydrogen Energy, Vol. 30, 6044-6060, 2014. http://dx.doi.org/10.1016/j.ijhydene.2014.01.182 Application of a Laser Induced Fluorescence
More informationCFD Simulation of Internal Flowfield of Dual-mode Scramjet
CFD Simulation of Internal Flowfield of Dual-mode Scramjet C. Butcher, K. Yu Department of Aerospace Engineering, University of Maryland, College Park, MD, USA Abstract: The internal flowfield of a hypersonic
More informationPreliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data
35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander
More informationAvailable online at ScienceDirect. Procedia Engineering 99 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 99 (5 ) 6 APISA4, 4 Asia-Pacific International Symposium on Aerospace echnology, APISA4 Effect of Combustion Heat Release on
More informationShock Wave Boundary Layer Interaction from Reflecting Detonations
Shock Wave Boundary Layer Interaction from Reflecting Detonations J. Damazo 1, J. Ziegler 1, J. Karnesky 2, and J. E. Shepherd 1 1 Introduction The present work is concerned with the differences in how
More informationComputer Aided Detailed Mechanism Generation for Large Hydrocarbons: n-decane
23 rd ICDERS July 24 29, 2011 Irvine, USA Computer Aided Detailed Mechanism Generation for Large Hydrocarbons: n-decane Martin Hilbig 1, Lars Seidel 1, Xiaoxiao Wang 1, Fabian Mauss 1 and Thomas Zeuch
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 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 informationOxidation of C 3 and n-c 4 aldehydes at low temperatures
Oxidation of C 3 and n-c 4 aldehydes at low temperatures M. Pelucchi*, A. Frassoldati*, E. Ranzi*, T. Faravelli* matteo.pelucchi@polimi.it * CRECK-Department of Chemistry, Materials and Chemical Engineering
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationA Priori Model for the Effective Lewis Numbers in Premixed Turbulent Flames
Paper # 070LT-0267 Topic: Turbulent Flames 8 th US National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013.
More informationNumerical investigation of flame propagation in small thermally-participating, adiabatic tubes
Paper # 7MI-181 Topic: Microcombustion and New Combustion Devices 8 th US National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationInvestigation of ignition dynamics in a H2/air mixing layer with an embedded vortex
Paper # 070LT-0211 The 8th US National Meeting of the Combustion Institute, Park City, UT, May 19-22, 2013 Investigation of ignition dynamics in a H2/air mixing layer with an embedded vortex S.K. Menon
More informationIgnition Delay Time of Small Hydrocarbons-Nitrous Oxide(-Oxygen) Mixtures
24 th ICDERS July 28 - August 2, 2013 Taipei, Taiwan Ignition Delay Time of Small Hydrocarbons-Nitrous Oxide(-Oxygen) Mixtures Rémy Mével and Joseph Shepherd Graduate Aerospace Laboratories, California
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationCombustion Behind Shock Waves
Paper 3F-29 Fall 23 Western States Section/Combustion Institute 1 Abstract Combustion Behind Shock Waves Sandeep Singh, Daniel Lieberman, and Joseph E. Shepherd 1 Graduate Aeronautical Laboratories, California
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
More informationEffect of Filter Choice on OH* Chemiluminescence Kinetics at Low and Elevated Pressures
7 th US National Technical Meeting of the Combustion Institute Hosted by the Georgia Institute of Technology, Atlanta, GA March 20-23, 2011 Effect of Filter Choice on OH* Chemiluminescence Kinetics at
More informationDNS of Reacting H 2 /Air Laminar Vortex Rings
46th AIAA Aerospace Sciences Meeting and Exhibit 7-10 January 2008, Reno, Nevada AIAA 2008-508 DNS of Reacting H 2 /Air Laminar Vortex Rings Jeff Doom and Krishnan Mahesh University of Minnesota, Minneapolis,
More informationFINAL REPORT. Office of Naval Research. entitled. Anatol Roshko Theodore Von Karman Professor of Aeronautics
to F11 F rnpv FINAL REPORT 4to Office of Naval Research on Contract No. N00014-85-C-0646 Work Unit No. 4328-434 entitled STRUCTURE AND MIXING IN TURBULENT SHEAR FLOWS 1 July 1985-31 October 1988 Anatol
More information41st Aerospace Sciences Meeting and Exhibit 6-9 January 2003, Reno, Nevada Modeling shock unsteadiness in shock/turbulence interaction
4st Aerospace Sciences Meeting and Exhibit 6-9 January 003, Reno, Nevada Modeling shock unsteadiness in shock/turbulence interaction AIAA 003-65 Krishnendu Sinha*, Krishnan Mahesh and Graham V. Candler
More informationA validation study of the flamelet approach s ability to predict flame structure when fluid mechanics are fully resolved
Center for Turbulence Research Annual Research Briefs 2009 185 A validation study of the flamelet approach s ability to predict flame structure when fluid mechanics are fully resolved By E. Knudsen AND
More informationHIGH PRESSURE METHANE-OXYGEN COMBUSTION KINETIC ANALYSIS
HIGH PRESSURE METHANE-OXYGEN COMBUSTION KINETIC ANALYSIS G. Saccone*, P. Natale*, F. Battista* g.saccone@cira.it p.natale@cira.it f.battista@cira.it *CIRA Italian Aerospace Research Centre, Capua Italy,
More informationinter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE
Copyright SFA - InterNoise 2000 1 inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering 27-30 August 2000, Nice, FRANCE I-INCE Classification: 7.0 VIBRATIONS OF FLAT
More informationNUMERICAL SIMULATION OF HYDROGEN COMBUSTION. Jan-patrice SIMONEAU, FRAMATOME - FRANCE
FR0200503 NUMERICAL SIMULATION OF HYDROGEN COMBUSTION Jan-patrice SIMONEAU, FRAMATOME - FRANCE Novatome - 10, rue Juliette Recamier- F 69456 LYON cedexo6 - France Ph : +33 4 72 74 73 75 - Facs : +33 4
More informationAn Experimental Study of Butene Isomers and 1,3-Butadiene Ignition Delay Times at Elevated Pressures
An Experimental Study of Butene Isomers and 1,3-Butadiene Ignition Delay Times at Elevated Pressures Yang Li, 1, Eoin O Connor 1, Chong-Wen Zhou 1, Henry J. Curran 1 1 Combustion Chemistry Centre, National
More informationOn fully developed mixed convection with viscous dissipation in a vertical channel and its stability
ZAMM Z. Angew. Math. Mech. 96, No. 12, 1457 1466 (2016) / DOI 10.1002/zamm.201500266 On fully developed mixed convection with viscous dissipation in a vertical channel and its stability A. Barletta 1,
More informationFluid Mechanics Theory I
Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to
More informationShock tunnel operation and correlation of boundary layer transition on a cone in hypervelocity flow
Shock tunnel operation and correlation of boundary layer transition on a cone in hypervelocity flow J.S. Jewell 1, J.E. Shepherd 1, and I.A. Leyva 2 1 Introduction The Caltech T reflected shock is used
More informationSummary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer
1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic
More informationUnsteady Wave Motion Shock Tube Problem - Shock Reflection
Unsteady Wave Motion Shock Tube Problem - Shock Reflection Niklas Andersson Division of Fluid Dynamics Department of Applied Mechanics Chalmers University of Tecnology 8 februari 09 Shock Reflection When
More informationCOMBUSTION OF THE BUTANOL ISOMERS: REACTION PATHWAYS AT ELEVATED PRESSURES FROM LOW-TO-HIGH TEMPERATURES
COMBUSTION OF THE BUTANOL ISOMERS: REACTION PATHWAYS AT ELEVATED PRESSURES FROM LOW-TO-HIGH TEMPERATURES Michael R. Harper, William H. Green* Massachusetts Institute of Technology, Department of Chemical
More informationA Novel FEM Method for Predicting Thermoacoustic Combustion Instability
Excerpt from the Proceedings of the COMSOL Conference 009 Milan A Novel FEM Method for Predicting Thermoacoustic Combustion Instability G. Campa *, S.M. Camporeale Politecnico di Bari * campa@imedado.poliba.it,
More informationFlame Propagation in Poiseuille Flow under Adiabatic Conditions
Flame Propagation in Poiseuille Flow under Adiabatic Conditions J. DAOU and M. MATALON* Department of Engineering Sciences and Applied Mathematics, Northwestern University, Evanston, IL 60208-3125, USA
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationInitiation of stabilized detonations by projectiles
Initiation of stabilized detonations by projectiles P. Hung and J. E. Shepherd, Graduate Aeronautical Laboratory, California Institute of Technology, Pasadena, CA 91125 USA Abstract. A high-speed projectile
More informationChemical Kinetics of Ethane Oxidation and Methane Oxidation with Platinum
Abstract Chemical Kinetics of Ethane Oxidation and Methane Oxidation with Platinum Jerry J. Zhang University of Southern California Professor Kenneth Brezinsky University of Illinois at Chicago Aleksandr
More informationLOW TEMPERATURE MODEL FOR PREMIXED METHANE FLAME COMBUSTION
ISTP-16, 2005, PRAGUE 16TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA LOW TEMPERATURE MODEL FOR PREMIXED METHANE FLAME MBUSTION M. Forman, J.B.W.Kok,M. Jicha Department of Thermodynamics and Environmental
More informationComputation of Shock Waves in Inert Binary Gas Mixtures in Nonequilibrium Using the Generalized Boltzmann Equation
50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 09-12 January 2012, Nashville, Tennessee AIAA 2012-0361 Computation of Shock Waves in Inert Binary Gas Mixtures
More informationLES of an auto-igniting C 2 H 4 flame DNS
Center for Turbulence Research Annual Research Briefs 2011 237 LES of an auto-igniting C 2 H 4 flame DNS By E. Knudsen, E. S. Richardson, J. H. Chen AND H. Pitsch 1. Motivation and objectives Large eddy
More informationFigure 11.1: A fluid jet extruded where we define the dimensionless groups
11. Fluid Jets 11.1 The shape of a falling fluid jet Consider a circular orifice of a radius a ejecting a flux Q of fluid density ρ and kinematic viscosity ν (see Fig. 11.1). The resulting jet accelerates
More informationLarge-eddy simulation of an industrial furnace with a cross-flow-jet combustion system
Center for Turbulence Research Annual Research Briefs 2007 231 Large-eddy simulation of an industrial furnace with a cross-flow-jet combustion system By L. Wang AND H. Pitsch 1. Motivation and objectives
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationThermoacoustic Instabilities Research
Chapter 3 Thermoacoustic Instabilities Research In this chapter, relevant literature survey of thermoacoustic instabilities research is included. An introduction to the phenomena of thermoacoustic instability
More information1. (20 pts total 2pts each) - Circle the most correct answer for the following questions.
ME 50 Gas Dynamics Spring 009 Final Exam NME:. (0 pts total pts each) - Circle the most correct answer for the following questions. i. normal shock propagated into still air travels with a speed (a) equal
More informationDynamics of Zonal Shear Collapse in Hydrodynamic Electron Limit. Transport Physics of the Density Limit
Dynamics of Zonal Shear Collapse in Hydrodynamic Electron Limit Transport Physics of the Density Limit R. Hajjar, P. H. Diamond, M. Malkov This research was supported by the U.S. Department of Energy,
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationarxiv: v1 [physics.ins-det] 9 Jun 2017
On the Uncertainty of Temperature Estimation in a Rapid Compression Machine arxiv:1706.04243v1 [physics.ins-det] 9 Jun 2017 Bryan W. Weber, Chih-Jen Sung, Michael Renfro Department of Mechanical Engineering,
More informationInsights into Model Assumptions and Road to Model Validation for Turbulent Combustion
Insights into Model Assumptions and Road to Model Validation for Turbulent Combustion Venke Sankaran AFRL/RQR 2015 AFRL/RQR Basic Research Review UCLA Jan 20, 2015 AFTC PA Release# 15011, 16 Jan 2015 Goals
More informationShock and Expansion Waves
Chapter For the solution of the Euler equations to represent adequately a given large-reynolds-number flow, we need to consider in general the existence of discontinuity surfaces, across which the fluid
More informationFLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam
More informationChapter 5 Phenomena of laminar-turbulent boundary layer transition (including free shear layers)
Chapter 5 Phenomena of laminar-turbulent boundary layer transition (including free shear layers) T-S Leu May. 3, 2018 Chapter 5: Phenomena of laminar-turbulent boundary layer transition (including free
More informationEVALUATION OF HYDROGEN, PROPANE AND METHANE-AIR DETONATIONS INSTABILITY AND DETONABILITY
EVALUATION OF HYDROGEN, PROPANE AND METHANE-AIR DETONATIONS INSTABILITY AND DETONABILITY Borzou, B. 1 and Radulescu, M.I. 2 1 Mechanical Engineering Department, University of Ottawa, 161 Louis Pasteur,
More informationCombustion Chemistry
Combustion Chemistry Hai Wang Stanford University 2015 Princeton-CEFRC Summer School On Combustion Course Length: 3 hrs June 22 26, 2015 Copyright 2015 by Hai Wang This material is not to be sold, reproduced
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationLecture 7 Detonation Waves
Lecture 7 etonation Waves p strong detonation weak detonation weak deflagration strong deflagration / 0 v =/ University of Illinois at Urbana- Champaign eflagrations produce heat Thermal di usivity th
More informationURANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver
URANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver *Yohan Choe 1), Hyeongjun Kim 1), Chongam Kim 2) 1), 2) Department of Aerospace Engineering, Seoul National
More informationApplication of Computational Fluid Dynamics in Discontinuous Unsteady Flow with Large Amplitude Changes; the Shock Tube Problem
Application of Computational Fluid Dynamics in Discontinuous Unsteady Flow with Large Amplitude Changes; the Shock Tube Problem KHALED ALHUSSAN Research Assistant Professor, Space Research Institute, King
More informationTHERMOCHEMICAL INSTABILITY OF HIGHLY DILUTED METHANE MILD COMBUSTION
THERMOCHEMICAL INSTABILITY OF HIGHLY DILUTED METHANE MILD COMBUSTION G. Bagheri*, E. Ranzi*, M. Lubrano Lavadera**, M. Pelucchi*, P. Sabia**, A. Parente***, M. de Joannon**, T. Faravelli* tiziano.faravelli@polimi.it
More informationShear Turbulence. Fabian Waleffe. Depts. of Mathematics and Engineering Physics. Wisconsin
Shear Turbulence Fabian Waleffe Depts. of Mathematics and Engineering Physics, Madison University of Wisconsin Mini-Symposium on Subcritical Flow Instability for training of early stage researchers Mini-Symposium
More informationMAE 210C FLUID MECHANICS III SPRING 2015 Rayleigh s equation for axisymmetric flows Forthestabilityanalysisofparallelaxisymmetricflowswithbasicvelocity profileū = (U x,u r,u θ ) = (U(r),0,0) weintroducemodalperturbationsoftheform(v
More informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
More informationAn analytical model for direct initiation of gaseous detonations
Issw21 An analytical model for direct initiation of gaseous detonations C.A. Eckett, J.J. Quirk, J.E. Shepherd Graduate Aeronautical Laboratories, California Institute of Technology, Pasadena CA 91125,
More informationUQ in Reacting Flows
UQ in Reacting Flows Planetary Entry Simulations High-Temperature Reactive Flow During descent in the atmosphere vehicles experience extreme heating loads The design of the thermal protection system (TPS)
More informationThis paper is part of the following report: UNCLASSIFIED
UNCLASSIFIED Defense Technical Information Center Compilation Part Notice ADP023624 TITLE: Ignition Kinetics in Fuels Oxidation DISTRIBUTION: Approved for public release, distribution unlimited This paper
More information