SEARCH · Engineering Papers

Results for “constant”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,711 records · Page 95

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Validity of the two-level model for Viterbi decoder gap-cycle performance

A two-level model has previously been proposed for approximating the performance of a Viterbi decoder which encounters data received with periodically varying signal-to-noise ratio. Such cyclically gapped data is obtained from the Very Large Array (VLA), either operating as a stand-alone system or arrayed with Goldstone. This approximate model predicts that the decoder error rate will vary periodically between two discrete levels with the same period as the gap cycle. It further predicts that the length of the gapped portion of the decoder error cycle for a constraint length K decoder will be about K-1 bits shorter than the actual duration of the gap. The two-level model for Viterbi decoder performance with gapped data is subjected to detailed validation tests. Curves showing the cyclical behavior of the decoder error burst statistics are compared with the simple square-wave cycles predicted by the model. The validity of the model depends on a parameter often considered irrelevant in the analysis of Viterbi decoder performance, the overall scaling of the received signal or the decoder's branch-metrics. Three scaling alternatives are examined: optimum branch-metric scaling and constant branch-metric scaling combined with either constant noise-level scaling or constant signal-level scaling. The simulated decoder error cycle curves roughly verify the accuracy of the two-level model for both the case of optimum branch-metric scaling and the case of constant branch-metric scaling combined with constant noise-level scaling. However, the model is not accurate for the case of constant branch-metric scaling combined with constant signal-level scaling.

Dolinar, S.↗

NASA NTRS1951

Investigation of dynamic characteristics of a turbine-propeller engine

Time constants that characterize engine speed response of a turbine-propeller engine over the cruising speed range for various values of constant fuel flow and constant blade angle were obtained both from steady-state characteristics and from transient operation. Magnitude of speed response to changes in fuel flow and blade angle was investigated and is presented in the form of gain factors. Results indicate that at any given value of speed in the engine cruising speed range, time constants obtained both from steady-state characteristics and from transient operation agree satisfactorily for any given constant fuel flow, whereas time constants obtained from transient operation exceed time constants obtained from steady-state characteristics by approximately 14 percent for any given blade angle.

ENGINES, TURBO-PROPELLER↗

NASA NTRS1928

The Gaseous Explosive Reaction : The Effect of Inert Gases

Attention is called in this report to previous investigations of gaseous explosive reactions carried out under constant volume conditions, where the effect of inert gases on the thermodynamic equilibrium was determined. The advantage of constant pressure methods over those of constant volume as applied to studies of the gaseous explosive reaction is pointed out and the possibility of realizing for this purpose a constant pressure bomb mentioned. The application of constant pressure methods to the study of gaseous explosive reactions, made possible by the use of a constant pressure bomb, led to the discovery of an important kinetic relation connecting the rate of propagation of the zone of explosive reaction within the active gases, with the initial concentrations of those gases: s = K(sub 1)(A)(sup n1)(B)(sup n2)(C)(sup n3)------. By a method analogous to that followed in determining the effect of inert gases on the equilibrium constant K, the present paper records an attempt to determine their kinetic effect upon the expression given above.

Stevens, F W↗

NASA NTRS1993

Microwave Scattering Model for Grass Blade Structures

In this paper, the electromagnetic scattering solution for a grass blade with complex cross-section geometry is considered. It is assumed that the blade cross section is electrically small, but its length is large compared to the incident wavelength. In a recent study it has been shown that the scattering solution for such problems, in the form of a polarizability tensor, can be obtained using the low-frequency approximation in conjunction with the method of moments. In addition, the study shows that the relationship between the polarizability tensor of a dielectric cylinder and its dielectric constant can be approximated by a simple algebraic expression. The results of this study are used to show that this algebraic approximation is valid also for cylinders with cross sections the shape of grass blades, providing that proper values am selected for each of three constants appearing in the expression. These constants are dependent on cylinder shape, and if the relationship between the constants and the three parameters describing a grass blade shape can be determined, an algebraic approximation relating polarizability tensor to blade shape, as well as dielectric constant, can be formed. Since the elements of the polarizability tensor are dependent on only these parameters, this algebraic approximation can replace the cumbersome method of moments model. A conjugate gradient method is then implemented to correctly determine the three constants of the algebraic approximation for each blade shape. A third-order polynomial fit to the data is then determined for each constant, thus providing a complete analytic replacement to the numerical (moment method) scattering model. Comparisons of this approximation to the numerical model show an average error of less than 3%.

Stiles, James M.↗

NASA NTRS1961

Flight-Path Characteristics for Decelerating From Supercircular Speed

Characteristics of the following six flight paths for decelerating from a supercircular speed are developed in closed form: constant angle of attack, constant net acceleration, constant altitude" constant free-stream Reynolds number, and "modulated roll." The vehicles were required to remain in or near the atmosphere, and to stay within the aerodynamic capabilities of a vehicle with a maximum lift-drag ratio of 1.0 and within a maximum net acceleration G of 10 g's. The local Reynolds number for all the flight paths for a vehicle with a gross weight of 10,000 pounds and a 600 swept wing was found to be about 0.7 x 10(exp 6). With the assumption of a laminar boundary layer, the heating of the vehicle is studied as a function of type of flight path, initial G load, and initial velocity. The following heating parameters were considered: the distribution of the heating rate over the vehicle, the distribution of the heat per square foot over the vehicle, and the total heat input to the vehicle. The constant G load path at limiting G was found to give the lowest total heat input for a given initial velocity. For a vehicle with a maximum lift-drag ratio of 1.0 and a flight path with a maximum G of 10 g's, entry velocities of twice circular appear thermo- dynamically feasible, and entries at velocities of 2.8 times circular are aerodynamically possible. The predominant heating (about 85 percent) occurs at the leading edge of the vehicle. The total ablated weight for a 10,000-pound-gross-weight vehicle decelerating from an initial velocity of twice circular velocity is estimated to be 5 percent of gross weight. Modifying the constant G load flight path by a constant-angle-of-attack segment through a flight- to circular-velocity ratio of 1.0 gives essentially a "point landing" capability but also results in an increased total heat input to the vehicle.

Luidens, Roger W.↗

NASA NTRS1999

Analyses of Fatigue Crack Growth and Closure Near Threshold Conditions for Large-Crack Behavior

A plasticity-induced crack-closure model was used to study fatigue crack growth and closure in thin 2024-T3 aluminum alloy under constant-R and constant-K(sub max) threshold testing procedures. Two methods of calculating crack-opening stresses were compared. One method was based on a contact-K analyses and the other on crack-opening-displacement (COD) analyses. These methods gave nearly identical results under constant-amplitude loading, but under threshold simulations the contact-K analyses gave lower opening stresses than the contact COD method. Crack-growth predictions tend to support the use of contact-K analyses. Crack-growth simulations showed that remote closure can cause a rapid rise in opening stresses in the near threshold regime for low-constraint and high applied stress levels. Under low applied stress levels and high constraint, a rise in opening stresses was not observed near threshold conditions. But crack-tip-opening displacement (CTOD) were of the order of measured oxide thicknesses in the 2024 alloy under constant-R simulations. In contrast, under constant-K(sub max) testing the CTOD near threshold conditions were an order-of-magnitude larger than measured oxide thicknesses. Residual-plastic deformations under both constant-R and constant-K(sub max) threshold simulations were several times larger than the expected oxide thicknesses. Thus, residual-plastic deformations, in addition to oxide and roughness, play an integral part in threshold development.

Newman, J. C., Jr.↗

NASA NTRS1999

Clay catalysis of oligonucleotide formation: kinetics of the reaction of the 5'-phosphorimidazolides of nucleotides with the non-basic heterocycles uracil and hypoxanthine

The montmorillonite clay catalyzed condensation of activated monocleotides to oligomers of RNA is a possible first step in the formation of the proposed RNA world. The rate constants for the condensation of the phosphorimidazolide of adenosine were measured previously and these studies have been extended to the phosphorimidazolides of inosine and uridine in the present work to determine of substitution of neutral heterocycles for the basic adenine ring changes the reaction rate or regioselectivity. The oligomerization reactions of the 5'-phosphoromidazolides of uridine (ImpU) and inosine (ImpI) on montmorillonite yield oligo(U)s and oligo(I)s as long as heptamers. The rate constants for oligonucleotide formation were determined by measuring the rates of formation of the oligomers by HPLC. Both the apparent rate constants in the reaction mixture and the rate constants on the clay surface were calculated using the partition coefficients of the oligomers between the aqueous and clay phases. The rate constants for trimer formation are much greater than those dimer synthesis but there was little difference in the rate constants for the formation of trimers and higher oligomers. The overall rates of oligomerization of the phosphorimidazolides of purine and pyrimidine nucleosides in the presence of montmorillonite clay are the same suggesting that RNA formed on the primitive Earth could have contained a variety of heterocyclic bases. The rate constants for oligomerization of pyrimidine nucleotides on the clay surface are significantly higher than those of purine nucleotides since the pyrimidine nucleotides bind less strongly to the clay than do the purine nucleotides. The differences in the binding is probably due to Van der Waals interactions between the purine bases and the clay surface. Differences in the basicity of the heterocyclic ring in the nucleotide have little effect on the oligomerization process.

NASA Discipline Exobiology↗

NASA NTRS1947

Equations for Adiabatic but Rotational Steady Gas Flows without Friction

This paper makes the following assumptions: 1) The flowing gases are assumed to have uniform energy distribution. ("Isoenergetic gas flows," that is valid with the same constants for the the energy equation entire flow.) This is correct, for example, for gas flows issuing from a region of constant pressure, density, temperature, end velocity. This property is not destroyed by compression shocks because of the universal validity of the energy law. 2) The gas behaves adiabatically, not during the compression shock itself but both before and after the shock. However, the adiabatic equation (p/rho(sup kappa) = C) is not valid for the entire gas flow with the same constant C but rather with an appropriate individual constant for each portion of the gas. For steady flows, this means that the constant C of the adiabatic equation is a function of the stream function. Consequently, a gas that has been flowing "isentropically",that is, with the same constant C of the adiabatic equation throughout (for example, in origination from a region of constant density, temperature, and velocity) no longer remains isentropic after a compression shock if the compression shock is not extremely simple (wedge shaped in a two-dimensional flow or cone shaped in a rotationally symmetrical flow). The solution of nonisentropic flows is therefore an urgent necessity.

Schaefer, Manfred↗

NASA NTRSDate not supplied

A Misperception in Reliability Growth Modelling

The Duane reliability growth model is n(t)/t = k t^-alpha (1) The reliability growth rate is alpha, the downward slope of n(t)/t versus t. It usually varies from 0.2 to 0.6. k is a constant. Crow used a 56-failure data set to illustrate reliability growth.1 A graphical Duane model fit to this data gives n(t)/t = 0.640 t^-0.283 (2) A problem in using the Duane-Crow reliability growth model is that it assumes that reliability growth continues and the failure rate decreases throughout the test period. It is more usual that reliability growth stops when the failure rated is low enough. Growth testing is often followed by testing with a low constant failure rate due to rare or uncorrectable failure modes. As more and more low constant rate acceptable failures accumulate after the period of reliability growth, the reliability growth time exponent alpha decreases toward zero. This occurs if constant rate failures are treated as occurring during the reliability growth period. It is more accurate to model a period of initial reliability growth followed by testing without repair to more accurately determine the final constant failure rate. This is done in the abcd model. n(t)/t = a t^-b + c from t = 0 to td (3) = c + d after td, where d = a td^-b (4) The term a t^-b describes the continuous reliability growth that continues out to time td and c is the constant uncorrected failure rate. The parameter d represents an additional constant failure rate due to correctable but uncorrected failure modes. After the reliability growth process is terminated, the failure rate n(t)/t = c + d.

Harry W Jones↗

DOE OSTI2020

General Method for Determining Redox Potentials without Electrolyte

A novel method to determine redox potentials without electrolyte is presented here. The method is based on a new ability to determine the dissociation constant, $K°_d$, for ion pairs formed between any radical anion and any inert electrolyte counterion. These dissociation constants can be used to determine relative shifts of redox potential as a function of electrolyte concentration, connecting referenced potentials determined with electrochemistry (with 0.1 M electrolyte) to electrolyte-free values. Pulse radiolysis created radical anions enabling determination of equilibrium constants for electron transfer between anions of donor and acceptor molecules as a function of electrolyte concentration in THF. The measurements determined “composite equilibrium constants”, $K_{eqC}$, which contain information about the dissociation constant for the electrolyte cations, X + , with the radical anions of both the donor, $K°_d$(D –• ,X + ) and the acceptor, $K°_d$(A –• ,X + ). Dissociation constants were obtained for a selection of radical anions with tetrabutylammonium (TBA + ). The electrolyte was found to shift the reduction potentials of small molecules 1-methylpyrene and trans-stilbene by close to +130 mV whereas oligo-fluorenes and polyfluorenes experienced shifts of only (+25 ± 6) mV due to charge delocalization weakening the ion pair. These shifts for reduction of aromatic hydrocarbon molecules are smaller than shifts of +232 and +451 mV seen previously for benzophenone radical anion with TBA + and Na + respectively where the charge on the radical anion is localized largely on one C=O bond, thus forming a more tightly bound ion pair.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Surface Dynamics of Charge Transport in LaAlO 3 /SrTiO 3 with Time-Resolved Kelvin Probe Force Microscopy

Time-resolved Kelvin probe force microscopy (tr-KPFM) was used to map surface potential dynamics in LAO/STO under its critical thickness, where it acts as an insulator, as polarizing and depolarizing electric fields were applied to planar electrodes in ambient conditions. Two time constants, for polarization and relaxation (i.e., depolarization) each, were derived and were shown to be independent of the direction of electrical field, revealing the reproducibility and recovery of the surfaces. Ambient aging or recent electrical history did not change the surface kinetics as determined by tr-KPFM. During polarization, the faster time constant (~3–6 s) is likely associated with the transport of H+ ions through the full length of the channel. The slower time constant (~40–50 s), on the other hand, is observed across only half the channel width near the higher-potential–more positive–electrode, which is likely associated with an electrochemical reaction such as the dissociation of water. This slower time constant drops to ~10–20 s toward the middle of the channel which could represent the transport of OH– ions. Furthermore, during relaxation, the fast time constant is consistent with polarization and is correlated to the same mechanism while the slower time constant is higher than it was during polarization and may represent the reabsorption of moisture from the atmosphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Bayesian Parameter Estimation of the k-ω Shear Stress Transport Model for Accurate Simulations of Impinging-Jet Heat Transfer

Given the lack of fusion-relevant component test facilities, current estimates of the thermo-fluid performance of plasma-facing components are based for the most part on numerical simulations. A major source of uncertainty in these simulations is the semiempirical turbulence (closure) models for the Reynolds stresses appearing in the governing Reynolds-averaged Navier-Stokes equations, which involve a set of constants that depend upon the flow. The objective of this study is to evaluate Bayesian parameter estimation of turbulence closure constants in ANSYS Fluent to model heat transfer in impinging jets. The Bayesian statistical calibration produces a probability distribution for these constants from experimental data; the maximum a posteriori estimates are then taken to be the calibrated constants, or parameters. The turbulence model constants are calibrated using an experimental study of a submerged jet of air impinging on a flat heated surface at Reynolds numbers Re = O(10 4 ) and impingement distance in jet diameters H/d = 2. Numerical predictions using the calibrated model parameters are then compared with those generated using the default constants. Predictions obtained with model parameters calibrated on datasets of two different sizes are compared to evaluate the effect of the number of calibration samples. Lastly, the extrapolative ability of the calibrated model is examined by predictions at a Re beyond the calibration values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Theoretical Understanding of Effects of Operating Modes on the Performance Durability of Solid Oxide Cells: A Comparison between Potentiostatic and Galvanostatic Operations

In solid oxide cells (SOCs), the choice between galvanostatic (constant current) and potentiostatic (constant voltage) modes does not significantly affect the performance of SOCs as long as the cell components remain unchanged and intact. However, the degradation of cell components, which leads to changes in electrochemical and physical properties of the cell, elevates the importance of the selected operating mode. This paper aims to investigate the effects of galvanostatic and potentiostatic operating modes on the evolving properties of SOCs and their subsequent influence on performance durability. Employing non-equilibrium thermodynamic analysis, a crucial approach for understanding the degradation phenomena within an active electrochemical system, this study aims to provide in-depth insights into how these operating modes affect the longevity and efficacy of SOCs. Key findings include: In cases where oxygen electrode (OE) degradation is accelerated by higher partial pressure of oxygen ( p O 2 ), operating under constant voltage electrolysis can mitigate the high p O 2 at the OE|electrolyte (OE|EL) interface. Conversely, if OE degradation occurs more rapidly under a lower p O 2 , constant current electrolysis is more effective in suppressing degradation by achieving a high p O 2 at the OE|EL interface. For degradation of the fuel electrode (FE) due to higher p O 2 , constant current electrolysis is beneficial for more stable performance, which helps maintain low p O 2 at the FE|EL interface. When FE degradation is accelerated by lower p O 2 , constant voltage electrolysis can avert low p O 2 at the FE|EL interface. In practical scenarios, more complex degradation mechanisms come into play, especially when p O 2 significantly deviates from initial conditions. Degradation in one electrode can influence p O 2 in the other electrode, a phenomenon more pronounced in potentiostatic than in galvanostatic electrolysis.

Electrochemistry↗

DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Effects of Critical Compression Ratio on Rating Gasoline Knock Propensity

It is common practice in the automotive industry to explore the knock limits of fuels on an engine by a comparison of the knock limited spark advance (KLSA) at threshold knock intensity. However, the knock propensity of gasolines can be rated by changing one of three metrics on a variable compression ratio Cooperative Fuels Research (CFR) octane rating engine while holding the other two variables constant: knock intensity, spark timing, and critical compression ratio. The operational differences between the standard research octane number (RON) rating and modern engine operation have been explored in three parts. The first part focused on the effects of lambda and knock characterization. The second part studied the effects of spark timing. This third part explores the knock ratings of several gasolines by comparing the critical compression ratios at constant combustion phasing and knock intensity. The threshold knock intensity was based on the standard octane rating D1 pickup or by maximum amplitude of pressure oscillations (MAPO) measured by a piezoelectric cylinder pressure transducer. Several Fuels for Advanced Combustion Engines (FACE) gasolines, primary reference fuels (PRFs), and toluene standardization fuels (TSFs) were tested on a CFR octane rating engine with advanced data acquisition equipment and a piezoelectric cylinder pressure transducer. These tests deviated from the ASTM D2699 standard octane rating procedure. For each test fuel, the CFR engine was operated at stoichiometry at a constant combustion phasing (CA50) and the compression ratio was modified until a threshold knock intensity was realized. It was found that the chemical composition of the fuels affected the relationship of critical compression ratios between the D1 knockmeter and piezoelectric pressure transducer knock intensity thresholds, as well as the measured combustion maximum pressure rise rate and spark timing setting for constant CA50. For highly aromatic fuels tested at a constant MAPO knock intensity threshold, it was found that the maximum pressure rise rate was two to three times higher than that of highly paraffinic fuels with similar RON and the spark advance was several crank angle degrees less for constant combustion phasing.

Kolodziej, Christopher P↗

DOE OSTI2024

Energy-conserving physics for nonhydrostatic dynamics in mass coordinate models

Abstract. Motivated by reducing errors in the energy budget related to enthalpy fluxes within the Energy Exascale Earth System Model (E3SM), we study several physics–dynamics coupling approaches. Using idealized physics, a moist rising bubble test case, and the E3SM's nonhydrostatic dynamical core, we consider unapproximated and approximated thermodynamics applied at constant pressure or constant volume. With the standard dynamics and physics time-split implementation, we describe how the constant-pressure and constant-volume approaches use different mechanisms to transform physics tendencies into dynamical motion and show that only the constant-volume approach is consistent with the underlying equations. Using time step convergence studies, we show that the two approaches both converge but to slightly different solutions. We reproduce the large inconsistencies between the energy flux internal to the model and the energy flux of precipitation when using approximate thermodynamics, which can only be removed by considering variable latent heats, both when computing the latent heating from phase change and when applying this heating to update the temperature. Finally, we show that in the nonhydrostatic case, for physics applied at constant pressure, the general relation that enthalpy is locally conserved no longer holds. In this case, the conserved quantity is enthalpy plus an additional term proportional to the difference between hydrostatic pressure and full pressure.

58 GEOSCIENCES↗

NASA NTRS1980

A laser flash photolysis-resonance fluorescence kinetics study of the reaction Cl/2P/ + CH4 yields CH3 + HCl

The technique of laser flash photolysis-resonance fluorescence is employed to study the kinetics of the reaction Cl(2P) + CH4 yields CH3 + HCl over the temperature range 221-375 K. At temperatures less than or equal to 241 K the apparent bimolecular rate constant is found to be dependent upon the identity of the chemically inert gases in the reaction mixture. For Cl2/CH4/He reaction mixtures (total pressure = 50 torr) different bimolecular rate constants are measured at low and high methane concentrations. For Cl2/CH4/CCl/He and Cl2/CH4/Ar reaction mixtures, the bimolecular rate constant is independent of methane concentration, being approximately equal to the rate constant measured at low methane concentrations for Cl2/CH4/He mixtures. These rate constants are in good agreement with previous results obtained using the discharge flow-resonance fluorescence and competitive chlorination techniques. At 298 K the measured bimolecular rate constant is independent of the identity of the chemically inert gases in the reaction mixture and in good agreement with all previous investigations. The low-temperature results obtained in this investigation and all previous investigations can be rationalized in terms of a model which assumes that the Cl(2P 1/2) state reacts with CH4 much faster than the Cl(2P 3/2) state. Extrapolation of this model to higher temperatures, however, is not straightforward.

Ravishankara, A. R.↗