Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical equilibrium”

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.

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

At least 343 records · Page 19

New Reduced Two-Time Step Method for Calculating Combustion and Emission Rates of Jet-A and Methane Fuel With and Without Water Injection

A simplified kinetic scheme for Jet-A, and methane fuels with water injection was developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) or even simple FORTRAN codes that are being developed at Glenn. The two time step method is either an initial time averaged value (step one) or an instantaneous value (step two). The switch is based on the water concentration in moles/cc of 1x10(exp -20). The results presented here results in a correlation that gives the chemical kinetic time as two separate functions. This two step method is used as opposed to a one step time averaged method previously developed to determine the chemical kinetic time with increased accuracy. The first time averaged step is used at the initial times for smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used with higher water concentrations, gives the chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure and temperature (T4). The simple correlations would then be compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates were then used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide and NOx were obtained for Jet-A fuel and methane with and without water injection to water mass loadings of 2/1 water to fuel. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, and pressure.

Molnar, Melissa↗

A Pseudo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

A multi-grid, flux-difference-split, finite-volume code, VULCAN, is presented for solving the elliptic and parabolized form of the equations governing three-dimensional, turbulent, calorically perfect and non-equilibrium chemically reacting flows. The space marching algorithms developed to improve convergence rate and or reduce computational cost are emphasized. The algorithms presented are extensions to the class of implicit pseudo-time iterative, upwind space-marching schemes. A full approximate storage, full multi-grid scheme is also described which is used to accelerate the convergence of a Gauss-Seidel relaxation method. The multi-grid algorithm is shown to significantly improve convergence on high aspect ratio grids.

White, J. A.↗

Beneficial effects of humic acid on micronutrient availability to wheat

Humic acid (HA) is a relatively stable product of organic matter decomposition and thus accumulates in environmental systems. Humic acid might benefit plant growth by chelating unavailable nutrients and buffering pH. We examined the effect of HA on growth and micronutrient uptake in wheat (Triticum aestivum L.) grown hydroponically. Four root-zone treatments were compared: (i) 25 micromoles synthetic chelate N-(4-hydroxyethyl)ethylenediaminetriacetic acid (C10H18N2O7) (HEDTA at 0.25 mM C); (ii) 25 micromoles synthetic chelate with 4-morpholineethanesulfonic acid (C6H13N4S) (MES at 5 mM C) pH buffer; (iii) HA at 1 mM C without synthetic chelate or buffer; and (iv) no synthetic chelate or buffer. Ample inorganic Fe (35 micromoles Fe3+) was supplied in all treatments. There was no statistically significant difference in total biomass or seed yield among treatments, but HA was effective at ameliorating the leaf interveinal chlorosis that occurred during early growth of the nonchelated treatment. Leaf-tissue Cu and Zn concentrations were lower in the HEDTA treatment relative to no chelate (NC), indicating HEDTA strongly complexed these nutrients, thus reducing their free ion activities and hence, bioavailability. Humic acid did not complex Zn as strongly and chemical equilibrium modeling supported these results. Titration tests indicated that HA was not an effective pH buffer at 1 mM C, and higher levels resulted in HA-Ca and HA-Mg flocculation in the nutrient solution.

Non-NASA Center↗

Summary of Simplified Two Time Step Method for Calculating Combustion Rates and Nitrogen Oxide Emissions for Hydrogen/Air and Hydrogen/Oxygen

A simplified single rate expression for hydrogen combustion and nitrogen oxide production was developed. Detailed kinetics are predicted for the chemical kinetic times using the complete chemical mechanism over the entire operating space. These times are then correlated to the reactor conditions using an exponential fit. Simple first order reaction expressions are then used to find the conversion in the reactor. The method uses a two time step kinetic scheme. The first time averaged step is used at the initial times with smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, temperature, and pressure. The second instantaneous step is used at higher water concentrations (greater than l x 10(exp -20)) moles per cc) in the mixture which gives the chemical kinetic time as a function of the instantaneous fuel and water mole concentrations, pressure and temperature (T(sub 4)). The simple correlations are then compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. This time is regressed over the complete initial conditions using the Excel regression routine. Chemical kinetic time equations for H2 and NOx are obtained for H2/Air fuel and for H2/O2. A similar correlation is also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium temperature (T(sub 4)) as a function of overall fuel/air ratio, pressure and initial temperature (T(sub 3)). High values of the regression coefficient R squared are obtained.

Marek, C. John↗

Simplified Two-Time Step Method for Calculating Combustion Rates and Nitrogen Oxide Emissions for Hydrogen/Air and Hydorgen/Oxygen

A simplified single rate expression for hydrogen combustion and nitrogen oxide production was developed. Detailed kinetics are predicted for the chemical kinetic times using the complete chemical mechanism over the entire operating space. These times are then correlated to the reactor conditions using an exponential fit. Simple first order reaction expressions are then used to find the conversion in the reactor. The method uses a two-time step kinetic scheme. The first time averaged step is used at the initial times with smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, temperature, and pressure. The second instantaneous step is used at higher water concentrations (> 1 x 10(exp -20) moles/cc) in the mixture which gives the chemical kinetic time as a function of the instantaneous fuel and water mole concentrations, pressure and temperature (T4). The simple correlations are then compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. This time is regressed over the complete initial conditions using the Excel regression routine. Chemical kinetic time equations for H2 and NOx are obtained for H2/air fuel and for the H2/O2. A similar correlation is also developed using data from NASA s Chemical Equilibrium Applications (CEA) code to determine the equilibrium temperature (T4) as a function of overall fuel/air ratio, pressure and initial temperature (T3). High values of the regression coefficient R2 are obtained.

Molnar, Melissa↗

Simplified Two-Time Step Method for Calculating Combustion and Emission Rates of Jet-A and Methane Fuel With and Without Water Injection

A simplified kinetic scheme for Jet-A, and methane fuels with water injection was developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) or even simple FORTRAN codes. The two time step method is either an initial time averaged value (step one) or an instantaneous value (step two). The switch is based on the water concentration in moles/cc of 1x10(exp -20). The results presented here results in a correlation that gives the chemical kinetic time as two separate functions. This two time step method is used as opposed to a one step time averaged method previously developed to determine the chemical kinetic time with increased accuracy. The first time averaged step is used at the initial times for smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used with higher water concentrations, gives the chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure and temperature (T4). The simple correlations would then be compared to the turbulent mixing times to determine the limiting rates of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide and NOx are obtained for Jet-A fuel and methane with and without water injection to water mass loadings of 2/1 water to fuel. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, and pressure.

Molnar, Melissa↗

ThermoBuild: Online Method Made Available for Accessing NASA Glenn Thermodynamic Data

The new Web site program "ThermoBuild" allows users to easily access and use the NASA Glenn Thermodynamic Database of over 2000 solid, liquid, and gaseous species. A convenient periodic table allows users to "build" the molecules of interest and designate the temperature range over which thermodynamic functions are to be displayed. ThermoBuild also allows users to build custom databases for use with NASA's Chemical Equilibrium with Applications (CEA) program or other programs that require the NASA format for thermodynamic properties. The NASA Glenn Research Center has long been a leader in the compilation and dissemination of up-to-date thermodynamic data, primarily for use with the NASA CEA program, but increasingly for use with other computer programs.

McBride, Bonnie↗

Diode Laser Sensors for Arc-Jet Characterization

The development and application of tunable diode laser (TDL) absorption sensors to monitor the health and operating conditions in the large-scale 60 MW arc-heated- plasma wind-tunnel at NASA Ames Research Center is reported. The interactive heating facility (THF) produces re-entry flow conditions by expanding the gas heated in a constricted plasma arc-heater to flow at high velocity over a model located in a test cabin. This facility provides the conditions needed to test thermal protective systems for spacecraft re-entering the earth s atmosphere. TDL sensors are developed to monitor gas flows in both the high-temperature constricted flow and the supersonic expansion flow into test cabin. These sensors utilize wavelength-tuned diode lasers to measure absorption transitions of atomic oxygen near 777.2 nm, atomic nitrogen near 856.8 nm, and atomic copper near 793.3 nm. The oxygen and nitrogen sensors measure the population density in exited electronic states of these atoms. The measurements combined with the assumption of local thermal and chemical equilibrium yield gas temperature (typically near 7,000K). The nitrogen and oxygen population temperatures are redundant, and their close agreement provides an important test of the local thermal equilibrium assumption. These temperature sensors provide time-resolved monitors of the operating conditions of the arc-heater and can be used to verify and control the test conditions. An additional TDL sensor was developed to monitor the copper concentration in the arc-heater flow yielding values as high as 13 ppm. Measurements of copper in the flow can identify flow conditions with unacceptably rapid electrode erosion, and hence this sensor provides valuable information needed to schedule maintenance to avoid costly arc-heater failure. TDL sensors were also developed for measurements in the test cabin, where absorption measurements of the populations of argon and molecular nitrogen in excited metastable electronic states established that the number density of these excited species is much lower than estimated using frozen-chemistry approximations. This key finding suggests that in the post-expansion region there is not a significant energy sequestration in electronically excited species. Finally, TDL measurements of atomic potassium seeded into the test cabin flow were used to directly measure the static temperature of the test gas. The results of this study illustrate the high potential of time-resolved TDL measurements for routine and economical sensing of arc-heater health (gas temperature and electrode erosion) as well as the time-resolved test-cabin-flow conditions in front of the model.

Hanson, Ronald K.↗

Calculation of Free-Atom Fractions in Hydrocarbon-Fueled Rocket Engine Plume

Free atom fractions (Beta) of nine elements are calculated in the exhaust plume of CH4- oxygen and RP-1-oxygen fueled rocket engines using free energy minimization method. The Chemical Equilibrium and Applications (CEA) computer program developed by the Glenn Research Center, NASA is used for this purpose. Data on variation of Beta in both fuels as a function of temperature (1600 K - 3100 K) and oxygen to fuel ratios (1.75 to 2.25 by weight) is presented in both tabular and graphical forms. Recommendation is made for the Beta value for a tenth element, Palladium. The CEA computer code was also run to compare with experimentally determined Beta values reported in literature for some of these elements. A reasonable agreement, within a factor of three, between the calculated and reported values is observed. Values reported in this work will be used as a first approximation for pilot rocket engine testing studies at the Stennis Space Center for at least six elements Al, Ca, Cr, Cu, Fe and Ni - until experimental values are generated. The current estimates will be improved when more complete thermodynamic data on the remaining four elements Ag, Co, Mn and Pd are added to the database. A critique of the CEA code is also included.

Verma, Satyajit↗

Thermal Analysis on Plume Heating of the Main Engine on the Crew Exploration Vehicle Service Module

The crew exploration vehicle (CEV) service module (SM) main engine plume heating is analyzed using multiple numerical tools. The chemical equilibrium compositions and applications (CEA) code is used to compute the flow field inside the engine nozzle. The plume expansion into ambient atmosphere is simulated using an axisymmetric space-time conservation element and solution element (CE/SE) Euler code, a computational fluid dynamics (CFD) software. The thermal analysis including both convection and radiation heat transfers from the hot gas inside the engine nozzle and gas radiation from the plume is performed using Thermal Desktop. Three SM configurations, Lockheed Martin (LM) designed 604, 605, and 606 configurations, are considered. Design of multilayer insulation (MLI) for the stowed solar arrays, which is subject to plume heating from the main engine, among the passive thermal control system (PTCS), are proposed and validated.

Wang, Xiao-Yen J.↗

Oxygen Fugacity at High Pressure: Equations of State of Metal-Oxide Pairs

Oxygen fugacity (fO2) varies by orders of magnitude in nature, and can induce profound changes in the chemical state of a substance, and also in the chemical equilibrium of multicomponent systems. One prominent area in high pressure geochemistry, in which fO2 is widely recognized as a principal controlling factor, is that of metal-silicate partitioning of siderophile trace elements (e.g., [1]). Numerous experiments have shown that high pressures and temperatures can significantly affect metal/silicate partitioning of siderophile and moderately siderophile elements. Parameterization of these experimental results over P, T, X, and fO2 can allow the observed siderophile element composition of the mantle to be associated with particular thermodynamic conditions [2]. However, this is best done only if quantitative control exists over each thermodynamic variable relevant to the experiments. The fO2 values for many of these partitioning experiments were determined relative to a particular metal-oxide buffer (e.g., Fe-FeO (IW), Ni-NiO (NNO), Co-CoO, Re-ReO2 (RRO)), but the parameterization of all experimental results is weakened by the fact that the pressure-induced relative changes between these buffer systems are imprecisely known.

Campbell A. J.↗

Implementation of Radiation, Ablation, and Free Energy Minimization Modules for Coupled Simulations of Hypersonic Flow

A description of models and boundary conditions required for coupling radiation and ablation physics to a hypersonic flow simulation is provided. Chemical equilibrium routines for varying elemental mass fraction are required in the flow solver to integrate with the equilibrium chemistry assumption employed in the ablation models. The capability also enables an equilibrium catalytic wall boundary condition in the non-ablating case. The paper focuses on numerical implementation issues using FIRE II, Mars return, and Apollo 4 applications to provide context for discussion. Variable relaxation factors applied to the Jacobian elements of partial equilibrium relations required for convergence are defined. Challenges of strong radiation coupling in a shock capturing algorithm are addressed. Results are presented to show how the current suite of models responds to a wide variety of conditions involving coupled radiation and ablation.

Gnoffo, Peter A.↗

NASA Tech Briefs, November 1995

The contents include: 1) Mission Accomplished; 2) Resource Report: Marshall Space Flight Center; 3) NASA 1995 Software of the Year Award; 4) Microbolometers Based on Epitaxial YBa2Cu3O(sub 7-x) Thin Films; 5) Garnet Random-Access Memory; 6) Fabrication of SNS Weak Links on SOS Substrates; 7) High-Voltage MOSFET Switching Circuit; 8) Asymmetric Switching for a PWM H-Bridge Power Circuit; 9) Better Ohmic Contacts for InP Semiconductor Devices; 10) Low-Bandgap Thermovoltaic Materials and Devices; 11) Digital Frequency-Differencing Circuit; 12) Imaging Magnetometer; 13) Computer-Assisted Monitoring of a Complex System; 14) Buffered Telemetry Demodulator; 15) Compact Multifunction Inspection Head; 16) Optical Detection of Fractures in Ceramic Diaphragms; 17) Eddy-Current Detection of Cracks in Reinforced Carbon/Carbon; 18) Apparent Thermal Conductivity of Multilayer Insulation; 19) Optimizing Misch-Metal Compositions in Metal Hydride Anodes; 20) Device for Sampling Surface Contamination; 21) Probabilistic Failure Assessment for Fatigue; 22) Probabilistic Fatigue and Flaw-Propagation Analysis; 23) Windows Program for Driving the TDU-850 Printer; 24) Subband/Transform MATLAB Functions for Processing Images; 25) Computing Equilibrium Chemical Compositions; 26) Program Processes Thermocouple Readings; 27) ICAN-Second-Generation Integrated Composite Analyzer; 28) Integrated Composite Analyzer with Damping Capabilities; 29) Computing Efficiency of Transfer of Microwave Power; 30) Program Calculates Power Demands of Electronic Designs; 31) Cost-Estimation Program; 32) Program Estimates Areas Required by Electronic Designs; 33) Program to Balance Mapped Turbopump Assemblies; 34) BiblioTech; 35) Controlling Mirror Tilt With a Bimorph Actuator; 36) Burst-Disk Device Simulates Effect of Pyrotechnic Device; 37) Bearing-Mounting Concept Accommodates Thermal Expansion; 38) Parallel-Plate Acoustic Absorbers for Hot Environments; 39) Adjustable-Length Strut Withstands Large Cyclic Loads; 40) Tool Indicates Contact Angles in Bearing Raceways; 41) Gravity Slides With Magnetic Braking; 42) High-Torque, Lightweight, Pneumatically Driven Wrench for Small Spaces; 43) Device for Testing Compatibility of an O-Ring; 44) Magnetic Heat Pump Containing Flow Diverters; 45) Variable-Tilt Helicopter Rotor Mast; 46) "Beach-Ball" Robotic Rovers; 47) Apparatus Would Measure Temperatures of Ball Bearings; 48) Flexible Borescope for Inspecting Ducts; 49) Texturing Copper To Reduce Secondary Emission of Electrons; 50) Automated Laser Cutting in Three Dimensions; 51) Algorithm Helps Monitor Engine Operation; 52) Flexible Revision of Data-Processing Communications; 53) Software for Managing the Use of Land; 54) Thermal Strap Increases Cryocooling Efficiency; 55) Reversible Nut With Engagement Indication; 56) Control Algorithms for Kinematically Redundant Manipulators; 57) Computed Hydrogen-Flow Splits in a Rocket Engine; 58) Pressure and Thermal Modeling of Rocket Launches; 59) Field of View of a Spacecraft Antenna: Analysis and Software; 60) Digital Controller for Laser-Beam-Steering Subsystem; 61) More About Beam-Steering Subsystem for Laser Communication; 62) Digital Controller for Laser-Beam-Steering Subsystem: Part 2; 63) Interface Circuit Board for Space-Shuttle Communications; 64) Automated Planning of Spacecraft Telecommunications; 65) Artifacts of Spectral Analysis of Instrument Readings; 66) Neural-Network Controller for Vibration Suppression; 67) Adaptive Finite-Element Computation in Fracture Mechanics; 68) Attitude Control for the Cassini Spacecraft; 69) Analytical Model for Fluid Dynamics in a Microgravity Environment; 70) Study of Rocket-Engine Joints Bonded by NVCU/NARloy-Z; 71) Improved Silicon Nitride for Advanced Heat Engines; 72) Parameters for Welding Aluminum/Lithium Alloys; 73) Lightweight Composite Intertank Structure; 74) Foil Patches Seal Small Vacuum Leaks; 75) Data Base on Cables and Connectors; 76) Effect of Clock Mode on Radiation Hardnessf an ADC; and 77) Fault-Tolerant Control for a Robotic Inspection System.

Source record↗

High C/O Chemistry and Weak Thermal Inversion in the Extremely Irradiated Atmosphere of Exoplanet WASP-12b

The carbon-to-oxygen ratio (C/O) in a planet provides critical information about its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 causes a carbide-dominated interior as opposed to the silicate-dominated composition as found on Earth; the solar C/O is 0.54. Theory, shows that high C/O leads to a diversity of carbon-rich planets that can have very different interiors and atmospheres from those in the solar system. Here we report the detection of C/O greater than or equal to 1 in a planetary atmosphere. The transiting hot Jupiter WASP-12b has a dayside atmosphere depleted in water vapour and enhanced in methane by over two orders of magnitude compared to a solar-abundance chemical equilibrium model at the expected temperatures. The observed concentrations of the prominent molecules CO, CH4, and H2O are consistent with theoretical expectations for an atmosphere with the observed C/O = 1. The C/O ratios are not known for giant planets in the solar system, although they are expected to equal the solar value. If high C/O ratios are common, then extrasolar planets are likely very different in interior composition, and formed very differently, from expectations based on solar composition, potentially explaining the large diversity in observed radii. We also find that the extremely irradiated atmosphere (greater than 2500 K) of WASP-12b lacks a prominent thermal inversion, or a stratosphere, and has very efficient day-night energy circulation. The absence of a strong thermal inversion is in stark contrast to theoretical predictions for the most highly irradiated hot-Jupiter atmospheres.

Madhusudhan, Nikku↗

Implicit Coupling Approach for Simulation of Charring Carbon Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver with nonequilibrium gas/surface interaction for simulation of charring carbon ablators can be performed using an implicit approach. The material thermal response code used in this study is the three-dimensional version of Fully Implicit Ablation and Thermal response program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation method. Coupling between the material response and flow codes is performed by solving the surface mass balance in flow solver and the surface energy balance in material response code. Thus, the material surface recession is predicted in flow code, and the surface temperature and pyrolysis gas injection rate are computed in material response code. It is demonstrated that the time-lagged explicit approach is sufficient for simulations at low surface heating conditions, in which the surface ablation rate is not a strong function of the surface temperature. At elevated surface heating conditions, the implicit approach has to be taken, because the carbon ablation rate becomes a stiff function of the surface temperature, and thus the explicit approach appears to be inappropriate resulting in severe numerical oscillations of predicted surface temperature. Implicit coupling for simulation of arc-jet models is performed, and the predictions are compared with measured data. Implicit coupling for trajectory based simulation of Stardust fore-body heat shield is also conducted. The predicted stagnation point total recession is compared with that predicted using the chemical equilibrium surface assumption

Ablation↗

On Ion Clusters in the Interstellar Gas

In a recent paper V.I. Krassovsky (1958) predicts the occurrence of clusters of large numbers of atoms and molecules around ions in the interstellar gas. He then proposes a number of physicochemical processes that would be considerably enhanced by the high particle density in such clusters. In particular, he suggests that absorption by negative ions formed in the clusters would account for the interstellar extinction without any necessity for the presence of grains. Because of the important consequences that ion clusters could have, it is necessary to examine their occurrence more fully. This note re-examines the formation of ion clusters in space and shows that even ion-molecule pairs are essentially non-existent. Ion clusters have been considered by Bloom and Margenau (1952) from the same point of view as that used by Krassovsky, whose basic reference (Joffe and Semenov 1933) unfortunately is not available. A different approach has been used by Eyring, Hirschfelder, and Taylor (1936) following the methods of chemical equilibrium. Both the references cited here enable one to conclude that clustering is negligible. Therefore, the treatment of Eyring et al. is more appropriate than the method of Bloom and Margenau, which depends on the statistical equilibrium of an atmosphere in a force field.

Donn, Bertram↗

Treatment of Overlapping Gaseous Absorption with the Correlated-k Method in Hot Jupiter and Brown Dwarf Atmosphere Models

The correlated-k method is frequently used to speed up radiation calculations in both one-dimensional and three-dimensional atmosphere models. An inherent difficulty with this method is how to treat overlapping absorption, i.e. absorption by more than one gas in a given spectral region. We have evaluated the applicability of three different methods in hot Jupiter and brown dwarf atmosphere models, all of which have been previously applied within models in the literature: (i) Random overlap, both with and without resorting and rebinning, (ii) equivalent extinction and (iii) pre-mixing of opacities, where (i) and (ii) combine k-coefficients for different gases to obtain k-coefficients for a mixture of gases, while (iii) calculates k-coefficients for a given mixture from the corresponding mixed line-by-line opacities. We find that the random overlap method is the most accurate and flexible of these treatments, and is fast enough to be used in one-dimensional models with resorting and rebinning. In three-dimensional models such as GCMs it is too slow, however, and equivalent extinction can provide a speed-up of at least a factor of three with only a minor loss of accuracy while at the same time retaining the flexibility gained by combining k-coefficients computed for each gas individually. Pre-mixed opacities are significantly less flexible, and we also find that particular care must be taken when using this method in order to to adequately resolve steep variations in composition at important chemical equilibrium boundaries. We use the random overlap method with resorting and rebinning in our one-dimensional atmosphere model and equivalent extinction in our GCM, which allows us to e.g. consistently treat the feedback of non-equilibrium chemistry on the total opacity and therefore the calculated P-T profiles in our models.

atmospheres↗

Wildfire-driven Thunderstorms Cause a Volcano-Like Stratospheric Injection of Smoke

Intense heating by wildfires can generate deep, smoke-infused thunderstorms, known as pyrocumulonimbus (pyroCb), which can release a large quantity of smoke particles above jet aircraft cruising altitudes. Injections of pyroCb smoke into the lower stratosphere have gained increasing attention over the past 15 years due to the rapid proliferation of satellite remote sensing tools. Impacts from volcanic eruptions and other troposphere-to-stratosphere exchange processes on stratospheric radiative and chemical equilibrium are well recognized and monitored. However, the role of pyroCb smoke in the climate system has yet to be acknowledged. Here, we show that the mass of smoke aerosol particles injected into the lower stratosphere from five near-simultaneous intense pyroCbs occurring in western North America on 12 August 2017 was comparable to that of a moderate volcanic eruption, and an order of magnitude larger than previous benchmarks for extreme pyroCb activity. The resulting stratospheric plume encircled the Northern Hemisphere over several months. By characterizing this event, we conclude that pyroCb activity, considered as either large singular events, or a full fire season inventory, significantly perturb the lower stratosphere in a manner comparable with infrequent volcanic intrusions.

Pyrocumulonimbus (pyroCb)↗