Gas-Kinetic Finite Volume Methods

Total Page:16

File Type:pdf, Size:1020Kb

Gas-Kinetic Finite Volume Methods Gas-Kinetic Finite Volume Methods K. Xu, L. Martinelli, A. Jameson Department of M.A.E, Princeton University, Prificeton N.J. 08544, USA 1 Introduction Gas-kinetic schemes developed from the BGK model have been successfully applied to 1-D and 2-D flows [1,2]. One of the advantages of the gas-kinetic approach over more conventional methods is realized when one considers the multidimensional Euler and Navier-Stokes equations. In fact, by requiring only a single scalar gas distribu- tion function f, the gas-kinetic approach greatly simplifies the calculation of mass, momentum, energy densities and their fluxes. In this paper we construct a novel 3-D method for the system of conservation laws of fluid flow. A Lax-Wendroff type scheme is developed from gas-kinetics and applied to solve time dependent problems. In the limit of infinite collision time (v = ~) our formulation reduces precisely to the kinetic representation of flux vector splitting for the Euler equations. Moreover, we introduce a new time-independent flux formulation which, when applied in conjunction with Jameson's multigrid time stepping scheme, yields efficient and accurate solutions of steady compressible flows. 2 Three Dimensional Finite Volume Gas-Kinetic Schemes In the finite volume method, the discretization is accomplished by dividing the flow into a large number of small subdomains, and applying the conservation laws in the integral form , s=0 f~ to each subdomain f~ with boundary 0fL In this equation U is the macroscopic state vector, defined as V = (p, P~, Py, Pz, e)T, where p,/~, and e are the mass, momentum, and energy densities, and F is the flux vector. In a gas-kinetic finite volume scheme the flux vectors across cell boundaries are constructed by computing the gas distribution function f. In three dimensions, we assume the BGK model as the governing equation for the distribution function f f~ + uf~ + vfy --t- Wfz : (g -- f)/'r (1) Here f is a function of space (x, y, z), time t, particle velocity (u, v, w) and internal variable ~ with K degrees of freedom(i.e. K = 2 for 7 = 1.4 gases). The relations 107 between mass p, momentum fi and energy c densities with the distribution function f are (P, P~, Py, Pz, ()T = JCjdZ, = 1,2, ...5 (2) where ¢, is the vector of moments ¢~ = (1, u,v,w,~(u1 2 + V2 + W 2 +~2)) T and dE = dudvdwd~ is the volume element in the phase space. In the BGK model the equilibrium state is described by a Maxwellian distribution g = Ae -~((~-U)~+(v-v)2+(w-W)2+~2) (3) where U,V and W are macroscopic gas velocities. For an equilibrium gas flow with f = g, the Euler equations in three dimensional space can be obtained by taking the moments of Ca in Eq.(1). On the other hand, to the first order of r, the Navier-Stokes equations, with a dynamic viscosity coefficient of v = rp (where p is the pressure), can be derived from the Chapman-Enskog expansion. Since mass, momentum and energy are conservative quantities in the process of gas evolution, f and g have to satisfy the conservation constraint f(g f)¢~dE 0, 0~ !, 2, ..:5 (4) at any point in space and time. The general solution for f in Eq.(1) in three dimensions at the position of (x, y, z) and t is f(x,y, z,t, u, v, w,~) = -;JogS1vt "x',y , z',t, u, w,~)e-(t-t')/r dt ' (5) + e-t/fro(x- ut,y- vt, z- wt) where z' = x - u(t - t'), y' = y - v(t - t'), z' = z - w(t - t') are the trajectory of a particle motion, and f0 is the initial nonequilibrium gas distribution function f at the beginning of each time step (t = 0). Two unknowns g and f0 must be determined in the above equation to obtain the solution f, so as to evaluate the fluxes across cell boundaries. We will consider the evaluation of fluxes across a boundary separating two cells in the x direction and, to simplify the notation, the point for evaluating fluxes at the cell boundary will be assumed at (x = 0, y = 0, z = 0). Generally, f0 and g can be expanded around the cell boundary as gl(1 + a~x + bly -4- clz), x < 0 fo= [.g~(l+a~x+b~y+c~z), z>0 (6) and g = g0(1 + ax + by+ ez + J]t) (7) where gt, g~ and go are local Maxwellian distribution functions. The dependence of at,b z, ...,A on the particle velocity can be obtained from the Taylor expansion of a Maxwellian such as -A = -~i J¢- -A2 'it + ~i3~ --{--~4 w Jr- .As(~ 2 + v 2 + w 2 --~ ~2) (8) 108 In order to get all the parameters in Eq.(6) and Eq.(7) at time t = 0, we need first interpolate the initial macroscopic variables. The interpolation is carried out by using a Symmetric Limited Positive (SLIP) technique originally developed by Jameson from Local Extremum Diminishing (LED) considerations [3]. All the coefficients in f0 can be obtained directly. Then, go in Eq.(7) at (x = 0, y = 0,z = 0,t = 0) can be evaluated automatically by taking the limit of Eq. (5) as t --* 0 and substituting into Eq.(4) to obtain j..+..=, .:,,.,...o <o, u>O u<O The other terms of ~, b and a in Eq.(8) at t = 0 can be computed from the new mass, momentum and energy interpolations which are continuous across the cell boundary in all three directions. Now, the only unknown term left in Eq.(8) is A. This can be evaluated as follows by substituting Eq.(6) and Eq.(7) into EQ.(5), we get f(O, O, O, t, u, v, w) = (1 - e-'lT)go -F (v(-1 -I- e -tiT) + te-tlr)(ua -t- vb -t- wa)go + r(t/r - 1 + e-q')/.go + e-qTfo(-ut,-vt,-wt) (10) with j" gl(1 :- aZut - b%t - cZwt), u > 0 fo(-ut, ~ wt~ (11) t gr(1 arut - brvt - Cwt), u < O Both f (Eq.(10)) and g (Eq.(7)) contain A. After applying the condition (4) at (x = 0, y = 0, z = 0) and integrating it over the whole time step T, such as L T i(g - f)¢c, dEdt = 0 (12) five moment equations of A can obtained, from which the five constants in A of Eq.(8) can be uniquely determined. Finally, the time-dependent numerical fluxes at a cell boundary can be computed. These fluxes satisfy the consistency condition Jc(U, U) = :P(U) for a homogeneous uniform flow, where ~c(U) is identical to the corresponding Euler fluxes in the 3-D case. Eq.(10) gives explicitly the time-dependent gas distribution function f at the cell boundary. Several limiting cases can be obtained, in particular: (i) in the hydrody- namic limit of r << T and in a smooth region one obtains schemes which are identical to these obtained from the one-step Lax-Wendroff scheme for the Navier-Stokes equa- tions, (ii) in the limit of r = oo corresponding to the collisionless Boltzmann equation, our formulation reduces precisely to the kinetic representation of flux vector splitting [4] for the Euler equations, (iii) for steady state calculations, the relaxation process can be simplified by ignoring all high-order spatial and temporal terms in the expan- sion of f and g to yield f= go + e-tlr (fo - go) = go + e(fo - go) The first term go represents the Euler fluxes, while the second term e(f0-g0) gives rise to a diffusion term which should be large near discontinuities to help fix a nonequi- 109 librium state. Thus, e -t/r can be regarded as a diffusive parameter which should be adapted to the local flow. It can he proved [5] that the entropy of go is always larger than that of f0, which guarantees that the local physical system will approach the state with larger entropy. This, for example, prevents the formation of unphysical rarefaction shocks. 3 Unsteady Flow Calculations A forward-facing step test is carried out on a uniform mesh of 240 × 80 cells and the numerical results are presented in Fig.(1.a) and Fig.(1.b) for the density and entropy distribution. Here, the collision time used is equivalent to a Reynolds number of Re ~_ 50,000 for the upstream gas flow, taking the wind tunnel height as the characteristic length scale. A slip boundary condition is imposed in order to avoid using finer meshes close to the boundary. There is no special treatment around the corner, and we never found any expansion shocks emerging from the corner. 4 Steady Flow Calculations with Multigrid Accel- eration The time independent gas-kinetic discretization scheme formulated in a previous sec- tion has been implemented in the full approximation multigrid time stepping scheme of Jameson [6]. Steady state transonic flow calculations for several airfoils have been performed to validate the simplified relaxation scheme. In these calculations, the selective parameter c is determined by a switching function calculated from local pressure gradients. Using subscripts i and j to label the mesh cells, the switching function for fluxes in the i direction is e = 1 - e -c~max(Pi+l,j'Pi,i) where a is a constant, and P/j is an appropriate pressure switch. The computational domain is an O-mesh with 160 cells in the circumferential direction and 32 cells in the radial direction. This is a fine enough mesh to produce accurate answers with standard high resolution difference schemes.
Recommended publications
  • Schlieren Imaging of Loud Sounds and Weak Shock Waves in Air Near the Limit of Visibility
    Shock Waves DOI 10.1007/s00193-009-0226-6 ORIGINAL ARTICLE Schlieren imaging of loud sounds and weak shock waves in air near the limit of visibility Michael John Hargather · Gary S. Settles · Matthew J. Madalis Received: 27 April 2009 / Accepted: 4 August 2009 © Springer-Verlag 2009 Abstract A large schlieren system with exceptional sensi- Keywords Shock waves · Nonlinear acoustics · Sound · tivity and a high-speed digital camera are used to visualize Noise · Schlieren · High-speed imaging loud sounds and a variety of common phenomena that pro- duce weak shock waves in the atmosphere. Frame rates var- PACS 42.79.Mt · 43.25.+y ied from 10,000 to 30,000frames/s with microsecond frame exposures. Sound waves become visible to this instrumenta- tion at frequencies above 10kHz and sound pressure levels 1 Introduction in the 110dB (6.3Pa) range and above. The density gradient produced by a weak shock wave is examined and found to Visualizing sound-wave motion in air has been a dream of sci- depend upon the profile and thickness of the shock as well as entists and engineers for at least a century and a half, although the density difference across it. Schlieren visualizations of sound waves and weak shock waves were initially confused. weak shock waves from common phenomena include loud August Toepler’s (1836–1912) epitaph, “He was the first to trumpet notes, various impact phenomena that compress a see sound,” actually referred to his schlieren imaging of weak bubble of air, bursting a toy balloon, popping a champagne shock waves from electric sparks [1].
    [Show full text]
  • Gas Kinetics in Galactic Disk. Why We Do Not Need a Dark Matter to Explain Rotation Curves of Spirals? Anton Lipovka
    Gas kinetics in galactic disk. Why we do not need a dark matter to explain rotation curves of spirals? Anton Lipovka To cite this version: Anton Lipovka. Gas kinetics in galactic disk. Why we do not need a dark matter to explain rotation curves of spirals?. 2018. hal-01832309v3 HAL Id: hal-01832309 https://hal.archives-ouvertes.fr/hal-01832309v3 Preprint submitted on 24 Jul 2020 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Gas kinetics in galactic disk. Why we do not need a dark matter to explain rotation curves of spirals? Lipovka A.A. Department of Research for Physics, Sonora University, 83000, Hermosillo, Sonora, México May 27, 2020 Abstract In present paper, we analyze how the physical properties of gas affect the Rotation Curve (RC) of a spiral galaxy. It is shown, that the observed non-Keplerian RCs measured for outer part of disks, and the observed radial gas distribution are closely related by the diffusion equation, which clearly indicates that no additional mass (Dark Matter) is need. It is stressed, that while the inner part of the RC is subject of the Kepler law, for the correct description of the outer part of the RC, the collisional property of the gas should also be taken into account.
    [Show full text]
  • Project RAND Publication: RA-15022
    , , ,-.. ,'. --. - ~.-~~ UNCLASSIFIED (rA) AERODYNAMICS, GAS DYNAMICS AND HEAT TRANSFER PROBLEMS OF A SATELLITE ROCKET RA-15022 February 1, 1947 UNCLASSI FI ED DOUGLAS AIRCRAFT COMPANY, INC. == DECLASSIFIED lAW E012958 BY EO Review Team DATE &' - ) r; --- :J-cJl r} This document contains information affecting the national defense of the United States within the meaning of the Espionage Act, SO U.S.C., 31 and 32, as amended. Its transmission or the revelation of its contents in any manner to an unauthorized person is prohibited by law. SECRET TABLE OF CONTENTS LIST OF FIGURES ...... , ......................................................... v SUMMARy .........................................................................vii LIST OF SyMBOLS ............................................................... ix I. AERODYNAMICS .............................................................. I A. Drag ....................................................................... l 1. The Drag Forces Exerted on a Supersonic Rocket ......................... 1 2. Approximations Used in the Drag Analysis ............................... 4 3. Methods Employed in Evaluating the Drag ................................ 6 4. The Optimum Shape for a Satellite Rocket ......... , ..................... 8 5. Rarefied Gas Dynamics (Superaerodynamics) ............................. 20 B. Aerodynamics of Stability and Control. .. , ..............•.... : ............. 20 L Control Manent ..................................... '............ , ...... 21 2. Body Moment ..........................................................
    [Show full text]
  • Gas Phase Reaction Kinetics of Neutral Oxygen Species
    National Bureau ot Standards Library, £01 Admin. Bldg. FEB 3 1969 A11102 mbQQb NSRDS-NBS 20 All 1021 46006 /NSflOS-MI? N6S-.UB-CW NBS PUBLICATIONS Gas Phase Reaction Kinetics Of Neutral Oxygen Species U.S. DEPARTMENT OF COMMERCE NATIONAL BUREAU OF STANDARDS National Standard Reference Data Series National Bureau of Standards National Standard Reference Data System, Plan of Operation, NSRDS-NBS 1 — 15 cents* Thermal Properties of Aqueous Uni-univalent Electrolytes, NSRDS-NBS 2 — 45 cents* Selected Tables of Atomic Spectra, Atomic Energy Levels and Multiplet Tables — Si II, Si III, Si IV, NSRDS-NBS 3, Section 1 — 35 cents* Selected Tables of Atomic Spectra, Atomic Energy Levels and Multiplet Tables — Si I, NSRDS-NBS 3, Section 2 — 20 cents* Atomic Transition Probabilities, Volume I, Hydrogen Through Neon, NSRDS-NBS 4 — $2.50* The Band Spectrum of Carbon Monoxide, NSRDS-NBS 5 — 70 cents* Tables of Molecular Vibrational Frequencies, Part 1, NSRDS-NBS 6 — 40 cents* High Temperature Properties and Decomposition of Inorganic Salts, Part 1, Sulfates, NSRDS-NBS 7-35 cents* Thermal Conductivity of Selected Materials, NSRDS-NBS 8 — 11.00* Tables of Bimolecular Gas Reactions, NSRDS-NBS 9 — $2.00* Selected Values of Electric Dipole Moments for Molecules in the Gas Phase, NSRDS- NBS 10-40 cents* Tables of Molecular Vibrational Frequencies, Part 2, NSRDS-NBS 11 — 30 cents* Tables for the Rigid Asymmetric Rotor: Transformation Coefficient for Symmetric to Asymmetric Bases and Expectation Values of P§, P|, and Pf, NSRDS-NBS 12 — 60 cents* Hydrogenation
    [Show full text]
  • Gas Kinetics Introduction
    Gas Kinetics Introduction • Most deposition techniques rely on gas flow in a vacuum. • We need a model for understanding; – the speed and energy of the gas molecules as a function of temperature and pressure, – how these molecules interact with each other and their surroundings, – and how mass, heat and momentum is transported by these molecules. The Size of a Molecule (a very rough calculation) • Take water (liquid H2 O) as an example. • Oxygen has 8 protons and 8 neutrons for a total of 16 nucleons and Hydrogen has 1 proton. So water has 16+1+1=18 nucleons −24 m≈PN1 m . ≈ 66 × 10g −23 m18HO= mN 2 = . 99 × g 10 2 d HO2 N1= m 3 = . 35 ×22 molecules 10 / g HO2 HO2 ρ 1= .g 00 cm3 / HO2 1 V = 2= . 99 × −23 10cm 3 molecule / HO2 ρ N HOHO2 2 3 d=2 ×3 V5 . = 64− ×8 10cm HO2 4 HO2 The Distance Between Molecules in a Gas • Now take water vapor (gaseous H2 O). ρ 0= . 8 ×g 10−3 cm3 / HO2 1 −20 3 l V = 3= . 73 × 10cm molecule / HO2 HO2 ρ N HOHO2 2 l=3 V3 . = 34 ×−7 10cm HO2 HO2 l≈ 6 d HO2 HO2 Molecular Velocities • Basic assumptions: – We’ll assume an ideal gas where the gas molecules interact elastically (collisions are similar to the collisions of hard billiard balls). – The distance between molecules are large compared to their sizes. – There are no attractive or repulsive forces between the molecules and each molecule moves independently of the others. Maxwell-Boltzmann Distribution • Under these assumptions, the molecules of a gas have velocities that are distributed according to: 1 dn M ⎛ Mv2 ⎞ where f is the fractional number of 2 ⎜ ⎟ molecules, v is the velocity, M is the f() v= =4πv exp⎜− ⎟ molecular weight, T is the temperature n dv 2πRT ⎝ 2RT ⎠ and R is the universal gas constant.
    [Show full text]
  • Mathematical Modeling of the Shockless Explosion Combustion
    Dissertation Mathematical Modeling of the Shockless Explosion Combustion Phillip Berndt im August 2016 eingereicht am Fachbereich Mathematik und Informatik Freie Universität Berlin zur Erlangung des Grades eines Doktors der Naturwissenschaften Erstgutachter & Betreuer: Prof. Dr.-Ing. Rupert Klein Zweitgutachter: Prof. Dr.-Ing. Christian Oliver Paschereit Drittgutachter: Prof. Dr. rer.nat. Volker John Tag der Disputation: 02. Dezember 2016 Contents 1. Introduction 1 1.1. The Shockless Explosion Combustion . .4 1.2. The Reactive Euler Equations of Fluid Dynamics . .7 1.2.1. The One-Dimensional Reactive Euler Equations . .8 1.2.2. Ignition Delay Time and Excitation Time . 11 1.2.3. Lagrangian Picture . 13 1.2.4. Quasi One-Dimensional Flow . 16 2. Fuel Requirements and Blend Selection 19 2.1. Model Reduction for SEC Fuels . 20 2.2. Model for Idealized SEC Fuels . 24 2.2.1. Single-Stage Model . 25 2.2.2. Two-Stage Extension . 27 2.3. Fuel Selection for the Atmospheric SEC Test Rig . 29 2.4. Numerical Methods for the Simulation of Chemical Kinetics . 32 2.5. Intermediate Conclusion . 35 3. Numerical Methods for Multi-Species Gaseous Flow 37 3.1. Solution of the Riemann Problem for Ideal Gases . 38 3.1.1. Linearization of the Equations . 38 3.1.2. Primitive Form of the Equations . 39 3.1.3. Riemann Invariants . 40 3.1.4. Solution Strategy . 42 3.2. The HLLE Solver in Multi-Species Simulations . 43 3.2.1. Positivity . 44 3.2.2. Pressure Oscillations at Contact Discontinuities . 46 3.2.3. A Note on Larrouturou’s Positivity Correction . 49 3.3.
    [Show full text]
  • On Ionospheric Aerodynamics*
    Prog. Aerospace Sei. 1975. Vo}. 16, No. 3. pp. 273-297. Pergamon Press. Printed in Great Britain ON IONOSPHERIC AERODYNAMICS* V. C. LIU The University of Michigan, Ann Arbor Summary--This paper presents theoretical methods to determine the gas dynamic and the electrostatic effects due to the interaction caused by a rapidly moving body in the ionosphere. The principles of the methods are derived from the kinetic theory of collision-free plasma. It is shown that the collective behavior of the collision-free plasma makes it possible to use the fluid approach to treat the problems of ionospheric aerodynamics. Various solutions to the system of fluid and field equations that have direct bearing on the ionospheric aerodynamics are presented and discussed. Physical significances of the mathematical results are stressed. Some outstanding unsolved problems in ionospheric aerodynamics are elaborated and discussed. LIST OF SYMBOLS AND UNITS (EM) l. INTRODUCTION All symbols are defined in the text. Only those that are repeatedly used are listed here. In the description of The advent of space flight which entails the electromagnetic fields Gaussian cgs units are adopted interaction between a rapidly moving, electrically throughout. These symbols are compiled in a separate charged body, e.g. a spacecraft, and a tenuous group with their corresponding quantities in rationalized ionized gas, e.g. the terrestrial ionosphere or the inks units with conversion factors listed in tandem. interplanetary gas, brings the aerodynamicists to GROUP A (electromagnetic) bear with a new aspect of aerodynamics--the ionospheric aerodynamics.t This interdisciplinary Gaussian Quantity mks science which combines gas dynamics with elec- B magnetic induction (47r/p.o)'r-B H magnetic field (4 ~-p.o)'nH tromagnetic fields poses some distinctive traits of c speed of light (~ot~o)-''~ its own even though it shares the basic principles of E electric field (4rr¢o)'r~E plasma physics with studies on electrical discharges 4, electric potential (4,n.eo)'1:~ and nuclear fusion.
    [Show full text]
  • Biosignature Gases in H2-Dominated Atmospheres on Rocky Exoplanets S
    The Astrophysical Journal, 777:95 (19pp), 2013 November 10 doi:10.1088/0004-637X/777/2/95 C 2013. The American Astronomical Society. All rights reserved. Printed in the U.S.A. BIOSIGNATURE GASES IN H2-DOMINATED ATMOSPHERES ON ROCKY EXOPLANETS S. Seager1,2,W.Bains1,3, and R. Hu1 1 Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA 2 Department of Physics, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA 3 Rufus Scientific Ltd., 37 The Moor, Melbourn, Herts SG8 6ED, UK Received 2013 April 2; accepted 2013 August 13; published 2013 October 18 ABSTRACT Super-Earth exoplanets are being discovered with increasing frequency and some will be able to retain stable H2- dominated atmospheres. We study biosignature gases on exoplanets with thin H2 atmospheres and habitable surface temperatures, using a model atmosphere with photochemistry and a biomass estimate framework for evaluating the plausibility of a range of biosignature gas candidates. We find that photochemically produced H atoms are the most abundant reactive species in H2 atmospheres. In atmospheres with high CO2 levels, atomic O is the major destructive species for some molecules. In Sun-Earth-like UV radiation environments, H (and in some cases O) will rapidly destroy nearly all biosignature gases of interest. The lower UV fluxes from UV-quiet M stars would produce a lower concentration of H (or O) for the same scenario, enabling some biosignature gases to accumulate. The favorability of low-UV radiation environments to accumulate detectable biosignature gases in an H2 atmosphere is closely analogous to the case of oxidized atmospheres, where photochemically produced OH is the major destructive species.
    [Show full text]
  • Nationaladvisorycommitt
    — . -. .!, + ..? .- - . .%- . .-. .. NATIONALADVISORY COMMITTEE FOR AERONAUTICS - TECHNICAL MEMORANDUM 1270 -J m # m OF VERY HIGH FLIGHT SPEEIX By Eugen Sanger 0 -4- Translation of ZWB Forschungsbericht Nr. 972, , May 19-38 i -+ Washington May 1950 . i. .. ----- . .. .. .... .-. , -L .- ● . NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS IIKENICAL WORANDUM 1270 TKE (MS KCNETIC!SOF WRY HI(2EFLIG3T SPEEDS* By Eugen Sanger ABSTRACT Ih ordinary gas dynamics we use assumptions which also agree with kbetic theoq of gases for smell mean free paths of the air molecules. The air forces thus calculated have to be r=xsmined, if the mean free path is comp=able with the dimensions of the moving body or even with its boundsry layer. This case is very clifficult to celculate. The con- ditions, however, are more simple, if the meen free path is large c- pared to the body length, so that the collisias of molecules with each other can be neglected compared to the collisions with We body surface. Ih order to study the influence of the lsrge mean free path, calcu- lations ere first carried out for the case of extreme rarefaction. Furthermore, the calculatims on this “completely ideal” gas will.be carried out under consideration of -Maxwelltsveloci@ distribution and under the assumption of certain experimental.lyestablished reflection laws for the trsnsl.ational end nontrenslational molecular de~ees of freedao The results thus obtained allow us to find, besides the air pressure forces perpendicular to the surface, also the friction stresses parsJ2el to the surface. — Their general results are calculated out for two practically important cases: for the thin smooth plate end for a projectile-shaped body moving axisUy.
    [Show full text]
  • Kinetic Measurements Using Flow Tubes the Journal of Physical Chemistry, Vol
    Kinetic Measurements Using Flow Tubes The Journal of Physical Chemistry, Vol. 83, No. 1, 1979 3 for diatomics substantial progress has been achieved, as bring experimenitalists and theorists together and to show evidenced by some of the present papers. that the field of thermal elementary reaction kinetics is Lastly, the symposium did achieve its major goal: to alive and well. Kinetic Measurements Using Flow Tubes Carleton J. Howard Aeronomy Laboratory, NOAA Environmental Research Laboratory, Boulder, Colorado 80303 (Received October 6, 1978) Publication costs assisted by the NOAA Environmental Research Laboratory This paper is a selective survey of the chemical kinetic literature involving flow tube measurements of elementary reaction rate constants. It describes the origins of the flow tube method, the experimental technique, the measurement of rate constants, and an analysis of the inherent errors. Emphasis is placed on the discussion of the strengths and limitations of the method as a source of kinetic data. I. Introduction TABLE I: Comparison of Flow Tube and Flash-Photo- In recent years there has been an increasing demand for lysis Kinetic Tecliniques gas phase reaction rate data. Laser development,l at- flow tube flash photolysis mospheric chemi~try,~,~and combustion4 are examples of temperature 200-600 K 100-600 K fields of application of elementary reaction studies. range Committees and organizations have been created to collect, pressure 1-10 torr 5 torr-several evaluate, and disseminate kinetic information. In at- range atmospheres mospheric chemistry, for example, there are serious rate constant (1O-lo-lO-l ') cm3 (10~'o-10~'8)cm3 economic and social implications derived from the ap- range molecule-' s-' molecule-' s-l detection excellent requires fast plication of kinetic data in computer models that assess versatility detector the impact of anthropogenic chemicals on stratospheric reactant excellent limited ozone.
    [Show full text]
  • Gas Kinetics, Vacuum Technology 1. A) Write a Compact Sentence That Describes the Molecular Origin of Gas Pressure, Including Its Temperature and Mass Dependence
    3.155J/6.152J Microelectronic Processing Technology Fall Term, 2003 Bob O'Handley Martin Schmidt Quiz # 1 Out: Oct. 22, 2003 Due: Oct. 22, 2003 Gas kinetics, vacuum technology 1. a) Write a compact sentence that describes the molecular origin of gas pressure, including its temperature and mass dependence. [5] A: P is due to change in momentum Dp = m Dv of molecules striking a surface; velocity increases as T so and accounting for the flux of m Dp ~ mT nv atoms hitting the surface Dp µ J mk T µ x mk T , where n is the B 2 B number density of molecules in chamber. b) In a 50 cm diameter vacuum chamber at room temperature and a pressure of 10 † mT, mostly N2 ( assume a nitrogen molecular diameter of 2 Å), what is the mean free path of the N2 molecules and describe their transport as either viscous or ballistic. [6] A: l = 1.75 cm, viscous. c) You want to reduce the pressure in your diffusion-pumped vacuum chamber from 10-5 Torr to 10-6 Torr. Comment on the efficacy of the following options. i) There is a good diffusion pump 5 feet from your chamber; it is known to produce a vacuum of 10-7 Torr in its chamber. You run a 1/2” tygon tube from that diffusion pump to an inlet valve on your chamber. [3] ii) Bake out your system at about 150oC overnight using heating tape. [3] iii) Pump liquid nitrogen through an unused substrate holder in your chamber. [3] A: i) Conductance of a long thin tube is too poor to allow much pumping; 1/2” diameter ≈ l, so flow is viscous.
    [Show full text]
  • The Atmospheric Chemistry of Hydrogen Cyanide (HCN)
    UC Irvine UC Irvine Previously Published Works Title The atmospheric chemistry of hydrogen cyanide (HCN) Permalink https://escholarship.org/uc/item/98d9x8wp Journal Journal of Geophysical Research, 88(C15) ISSN 0148-0227 Authors Cicerone, RJ Zellner, R Publication Date 1983 DOI 10.1029/JC088iC15p10689 License https://creativecommons.org/licenses/by/4.0/ 4.0 Peer reviewed eScholarship.org Powered by the California Digital Library University of California JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 88, NO. C15, PAGES 10,689-10,696, DECEMBER 20, 1983 The AtmosphericChemistry of Hydrogen Cyanide (HCN) R. J. CICERONE National Centerfor AtmosphericResearch R. ZELLNER UniversitdtGSttin•Ten, lnstitut fur PhysikalischeChemie Since 1981,three groupshave reportedspectroscopic detections and measurementsof hydrogency- anidein the atmosphere.HCN concentrations(volume mixing ratios) of (1.5-1.7)x 10-•ø appearto characterizethe stratosphereand the northernhemisphere's nonurban troposphere. In this paper we explorethe atmosphericbehavior of HCN by examiningits chemicaland photochemicalproperties. Its principalsinks are reactionswith atmosphericOH and O(•D); precipitationappears to be a negligible sink. In the stratosphere,vacuum UV photonsalso attack HCN. Atmosphericmodel calculationsshow that HCN shouldbe relativelywell mixedin the troposphereand that its concentrationdecreases slowly with altitudein the stratosphere.Its atmosphericresidence time appearsto be about 2.5 years,although 1-5 yearsis a possiblerange. To maintain the observedatmospheric burden of HCN, an annual source of about 2 x 10TM g nitrogenas HCN is required;we speculateas to the identityof thesesources. Oxidation of HCN by OH, while the major sink for atmosphericHCN, is not simpleor direct.Instead, oxidation proceedsfrom the HCN-OH adduct formed in HCN + OH reactions.These pathways and their uncertainties are outlined here.
    [Show full text]