Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “root productivity”

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 721 records · Page 40

Small-x structure of oxygen and neon isotopes as seen by the Large Hadron Collider

Results on collisions of O-16 nuclei performed at the Relativistic Heavy Ion Collider (RHIC) have been presented for the first time at Quark Matter 2023 by the STAR collaboration. O-16+O-16 collisions are also expected to take place in the near future at the Large Hadron Collider (LHC) at much higher beam energies. We explore the potential of beam-energy-dependent studies for this system to probe small-x dynamics and QCD evolution. We perform 3+1D IP-Glasma simulations to predict the rapidity dependence of the initial geometry of light-ion collisions, focusing on O-16+O-16 and Ne-20+Ne-20 collisions at root s(NN) = 70 GeV and 7 TeV. The choice of Ne-20 is motivated by its strongly elongated geometry, which may respond differently to the effect of the high-energy evolution compared to the more spherical O-16. We find that smearing induced by soft gluon production at high energy causes mild variations in the initial-state eccentricities as a function of the collision energy. These effects could be resolved in future experiments and deserve further investigation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Microthruster Propulsion for the Space Technology 7 (ST7) Technology Demonstration Mission

For future applications to precision formation flying missions, NASA's New Millennium Program is scheduled to test colloid micro-Newton thrusters (CMNTs) on the ST7 technology demonstration mission. These CMNTs are part of a disturbance reduction system (DRS) on the ESA SMART-2 Spacecraft or LISA Pathfinder. The goal of the ST7 DRS is to demonstrate technologies necessary to meet the nanometer precision positioning control requirements of the LISA mission. In order to achieve these goals, the CMNTs are required to demonstrate a thrust resolution of less than 0.1 micro-N and a thrust noise of less than 0.1 micro-N/[square root]Hz for thrust levels between 5 and 30 micro-N. Developed by Busek Co. with support from JPL in testing an design, the CMNT has been developed over the last four years into a flight-ready microthrust system. The development, validation testing, and flight unit production of the CMNTs are described. Development tests and analysis include preliminary wear tests, propellant loading process verification, flow testing, and performance verification. Validation and flight unit verification includes thermal and structural analysis, life testing, thermal and dynamic load testing, and performance verification. Final delivery of the units is planned in 2007 with and planned launch and flight demonstration 2009.

microthruster↗

Using mini-rockwool blocks as growing media for limited-cluster tomato production

Rockwool is an excellent growing medium for the hydroponic production of tomato; however, the standard size rockwool blocks [4 x 4 x 2.5 inches (10 x 10 x 6.3 cm) or 3 x 3 x 2.5 inches (7.5 x 7.5 x 6.3 cm)] are expensive. The following experiments were conducted with less expensive minirock wool blocks (MRBs), on rayon polyester material (RPM) as a bench top liner, to reduce the production cost of tomatoes (Lycopersicon esculentum) grown in a limited-cluster, ebb and flood hydroponic cultivation system. Fruit yield for single-cluster plants growing in MRBs [2 x 2 x 1.6 inches (5 x 5 x 4 cm) and 1.6 x 1.6 x 1.6 inches (4 x 4 x 4 cm)] was not significantly different from plants grown in larger sized blocks (3 x 3 x 2.5 inches). When the bench top was lined with RPM, roots penetrated the RPM, and an extensive root mat developed between the RPM and the bench top. The fruit yield from plants on RPM was significantly increased compared to plants without RPM due to increases in fruit size and fruit number. RPM also significantly reduced the incidence of blossom-end rot. In a second experiment, single- and double-cluster plants were grown on RPM. Fruit yield for double-cluster plants was 40% greater than for single-cluster plants due to an increase in fruit number, although the fruit were smaller in size. As in the first experiment, fruit yield for all plants grown in MRBs was not significantly different from plants grown in the larger sized blocks. MRBs and a RPM bench liner are an effective combination in the production of limited-cluster hydroponic tomatoes.

Non-NASA Center↗

Impacts of Legacy and Contemporary Nitrogen Inputs on N 2 O and CO 2 Emissions in Miscanthus and Maize Cultivated Soils

ABSTRACT Nutrient inputs influence the sustainability of bioenergy crop production through contemporary (shortly after addition) and legacy effects (persisting over years) on microbial nitrogen (N) and carbon cycling, which contribute to greenhouse gas emissions. However, the relative importance of contemporary and legacy effects and how that could vary by crop functional types is poorly understood. Considering its rhizomatous roots and perennial growth, we hypothesized that Miscanthus × giganteus (M×g) would be more sensitive to legacy N fertilization and the historical context of its environment than an annual crop like maize. To test this hypothesis, we examined the effects of legacy and contemporary N inputs on nitrous oxide (N 2 O) and carbon dioxide (CO 2 ) emissions, as well as key N cycling genes in soils where M×g and maize were grown. A 150‐day soil incubation experiment was conducted using soils from a long‐term M×g and maize fertility experiment with three historic N fertilization rates (0, 112, and 336 kg N ha −1 year −1 ) and a contemporary amendment (60 mg N kg −1 ) with negative control (0 mg N kg −1 ). We observed significant increases in cumulative N 2 O emissions in Mxg soils relative to maize soils, particularly at higher legacy fertilization rates, while contemporary N had no significant effect. Bacterial amo A gene abundance, which plays a significant role in nitrification in nutrient‐rich soils, also increased with higher legacy fertilization rates in M×g soils but was unaffected by the contemporary N. In maize soils, legacy and contemporary N did not significantly affect N 2 O emissions, but cumulative CO 2 emissions and amo A gene abundance significantly increased. The abundances of nor B genes were not significantly influenced by either legacy fertilization or contemporary N amendments in either soil. Our findings demonstrate the greater importance of fertilization history over contemporary N in mediating soil N 2 O emissions, particularly for perennial bioenergy crops.

09 BIOMASS FUELS↗

germs-lab/PAPER_miscanthus_soil_response

Nutrient inputs influence the sustainability of bioenergy crop production through contemporary (shortly after addition) and legacy effects (persisting over years) on microbial nitrogen (N) and carbon cycling, which contribute to greenhouse gas emissions. However, the relative importance of contemporary and legacy effects and how that could vary by crop functional types is poorly understood. Considering its rhizomatous roots and perennial growth, we hypothesized that Miscanthus × giganteus (M×g) would be more sensitive to legacy N fertilization and the historical context of its environment than an annual crop like maize. To test this hypothesis, we examined the effects of legacy and contemporary N inputs on nitrous oxide (N2O) and carbon dioxide (CO2) emissions, as well as key N cycling genes in soils where M×g and maize were grown. A 150-day soil incubation experiment was conducted using soils from a long-term M×g and maize fertility experiment with three historic N fertilization rates (0, 112, and 336 kg N ha−1 year−1) and a contemporary amendment (60 mg N kg−1) with negative control (0 mg N kg−1). We observed significant increases in cumulative N2O emissions in Mxg soils relative to maize soils, particularly at higher legacy fertilization rates, while contemporary N had no significant effect. Bacterial amoA gene abundance, which plays a significant role in nitrification in nutrient-rich soils, also increased with higher legacy fertilization rates in M×g soils but was unaffected by the contemporary N. In maize soils, legacy and contemporary N did not significantly affect N2O emissions, but cumulative CO2 emissions and amoA gene abundance significantly increased. The abundances of norB genes were not significantly influenced by either legacy fertilization or contemporary N amendments in either soil. Our findings demonstrate the greater importance of fertilization history over contemporary N in mediating soil N2O emissions, particularly for perennial bioenergy crops.

Lee, Jaejin↗

An evaluation of the National Meteorological Center weekly hindcast of upper-ocean temperature along the eastern Pacific equator in January 1992

The upper-ocean temperature distribution along the Pacific equator from 139 to 103 deg W was observed in January 1992 with temperature profiles recorded from a ship and inferred from an ocean general circulation model calculation involving data assimilation (i.e., hindcast). An El Nino episode was in progress. The l00-m-thick mixed layer depth, the mixed-layer temperature, and the depth-averaged temperature below the thermocline were similar in both data products. Considerable differences occurred in the representation of the 15-25 C thermocline, such as the depth-averaged temperatures above and below the 20 C isotherm, the east-west slope of the 20 C isotherm, and a 1000-km-wide depression. The longitudinal-averaged root-mean-square difference between the hindcast and observed depths of the center of the thermocline was 17 m. Most of the disparities could be attributed to a high wavenumber transient event that the model-based assimilation system was not intended to resolve.

Halpern, David↗

Production of muonic kaon atoms at high-energy colliders

We develop a framework for the formation of exotic muonic kaon atoms ($K\mu$) in semileptonic $D^{0}$ decays, using the effective weak Hamiltonian, a helicity-based treatment of the leptonic current, and a nonrelativistic bound-state projection. The resulting branching ratio, $\mathrm{BR}(D^{0}\!\to(K\mu )\nu_{\mu})=2.29\times10^{-10}$, is implemented in a ROOT-based code to estimate yields at RHIC, LHC, and STCF. We show quantitatively that $K\mu$ atoms—also produced through coalescence in the quark--gluon plasma (QGP)—provide a sensitive probe of low-momentum primordial muons and early-time electromagnetic radiation, offering complementary constraints in an otherwise unexplored phase space for thermal dilepton and photon emission. Newly estimated dissociation cross sections in detector material indicate that secondary‑vertex reconstruction should be experimentally feasible, allowing clean experimental identification of the atoms. Projected yields from QGP coalescence in LHC and RHIC heavy-ion collisions, and from $D^{0}$ decays in LHC high-luminosity $p+p$ collisions indicate that the first observation of $K\mu$ atoms is within reach. Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Article funded by SCOAP3 and published under licence by Chinese Physical Society and the Institute of High Energy Physics of the Chinese Academy of Science and the Institute of Modern Physics of the Chinese Academy of Sciences and IOP Publishing Ltd.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simulation model for plant growth in controlled environment systems

The role of the mathematical model is to relate the individual processes to environmental conditions and the behavior of the whole plant. Using the controlled-environment facilities of the phytotron at North Carolina State University for experimentation at the whole-plant level and methods for handling complex models, researchers developed a plant growth model to describe the relationships between hierarchial levels of the crop production system. The fundamental processes that are considered are: (1) interception of photosynthetically active radiation by leaves, (2) absorption of photosynthetically active radiation, (3) photosynthetic transformation of absorbed radiation into chemical energy of carbon bonding in solube carbohydrates in the leaves, (4) translocation between carbohydrate pools in leaves, stems, and roots, (5) flow of energy from carbohydrate pools for respiration, (6) flow from carbohydrate pools for growth, and (7) aging of tissues. These processes are described at the level of organ structure and of elementary function processes. The driving variables of incident photosynthetically active radiation and ambient temperature as inputs pertain to characterization at the whole-plant level. The output of the model is accumulated dry matter partitioned among leaves, stems, and roots; thus, the elementary processes clearly operate under the constraints of the plant structure which is itself the output of the model.

Raper, C. D., Jr.↗

Carbon Dynamics in Vegetation and Soils

The overall goals of CD-08 team in Phase I were to quantify the contributions of different components of the carbon cycle to overall ecosystem carbon balance in Amazonian tropical forests and to undertake process studies at a number of sites along the eastern LBA transect to understand how and why these fluxes vary with site, season, and year. We divided this work into a number of specific tasks: (1) determining the average rate (and variability) of tree growth over the past 3 decades; (2) determining age demographics of tree populations, using radiocarbon to determine tree age; (3) assessing the rate of production and decomposition of dead wood debris; (4) determining turnover rates for organic matter in soils and the mean age of C respired from soil using radiocarbon measurements; and (5) comparing our results with models and constructing models to predict the potential of tropical forests to function as sources or sinks of C. This report summarizes the considerable progress made towards our original goals, which have led to increased understanding of the potential for central Amazon forests to act as sources or sinks of carbon with altered productivity. The overall picture of tropical forest C dynamics emerging from our Phase I studies suggests that the fraction of gross primary production allocated to growth in these forests is only 25-30%, as opposed to the 50% assumed by many ecosystem models. Consequent slow tree growth rates mean greater mean tree age for a given diameter, as reflected in our measurements and models of tree age. Radiocarbon measurements in leaf and root litter suggest that carbon stays in living tree biomass for several years up to a decade before being added to soils, where decomposition is rapid. The time lags predicted from 14C, when coupled with climate variation on similar time scales, can lead to significant interannual variation in net ecosystem C exchange.

Trumbore, Susan↗

Refining water and carbon fluxes modeling in terrestrial ecosystems via plant hydraulics integration

Plant hydraulics substantially affects terrestrial water and carbon cycles by modulating water transport and carbon assimilation. Despite improved drought simulations in certain ecosystems through their integration into land surface models (LSMs), the broader application of plant hydraulics in diverse ecosystems and hydroclimates is still underexplored. Here, in this study, we implemented the recently developed Noah-Multiparameterization Land Surface Model (Noah-MP LSM) equipped with a plant hydraulics scheme (Noah-MP-PHS) across 40 FLUXNET sites globally. Employing the Shuffled Complex Evolution-University of Arizona (SCE-UA) auto-calibration algorithm, we optimized key plant hydraulics parameters for these sites spanning eight vegetation types in both arid and humid climates. Noah-MP-PHS significantly improves the simulation of evapotranspiration (ET) and gross primary production (GPP) by better representing atmospheric and soil water stress compared to traditional soil hydraulic schemes (SHSs, such as Noah and CLM). The augmented Noah-MP-PHS models reduce surface flux overestimation and underestimation, exhibiting an average increase of 0.14 and 0.15 in Kling-Gupta Efficiency (KGE) compared to Noah and CLM, respectively. The explicit consideration of plant capacitance in PHS reveals substantial deep-layer and nocturnal root water uptake especially under dry conditions. We employed eXplainable Machine learning (XML) to quantify the model’s relative sensitivity to newly introduced leaf-, stem and root-related parameters in PHS. The sensitivity analysis reveals a rise in root parameter importance and a decline in leaf and stem parameters as conditions shift from humid to arid. These findings indicate that as aridity states vary, the most influential parameters affecting surface fluxes variation may change in parameter calibration for PHS applications. Our findings underscore the importance of incorporating plant hydraulics into LSMs to enhance simulations of terrestrial water and carbon dynamics. These findings are crucial for understanding ecosystem responses to global climate changes and guide the broader application of PHS at larger scales.

54 ENVIRONMENTAL SCIENCES↗

Effect of microgravity on stress ethylene and carbon dioxide production in sweet clover (Melilotus alba L.)

The study of higher plant growth and development in the microgravity (micro-g) environment continues to be a challenge. This is in part a result of the available flight qualified hardware with restrictive closed gas environments. This point is underscored by considering that gas exchange of seedlings grown in microgravity may be further limited owing to a thicker layer of water wicked onto the roots and to the absence of convective mixing. We hypothesized that seedlings grown under such conditions will experience greater hypoxia in microgravity than at Earth gravity, and thus produce greater stress ethylene. We compared flight and ground samples of sweet clover seedlings grown in the Fluid Processing Apparatus (FPA) during STS-57 and found them to contain extremely high levels of carbon dioxide (CO2) and stress ethylene. There were time dependent increases for both gases, and seedling growth was greatly inhibited. We repeated these experiments aboard STS-60 using modified chambers which increased, by fifty fold, the air available to the developing seedlings. Sweet clover seed germination and subsequent seedling growth to eight days within the FPA modified with a gas permeable membrane is not compromised by the microgravity environment.

Gallegos, Gregory L.↗

Winter wheat: A model for the simulation of growth and yield in winter wheat

The basic ideas and constructs for a general physical/physiological process level winter wheat simulation model are documented. It is a materials balance model which calculates daily increments of photosynthate production and respiratory losses in the crop canopy. The partitioning of the resulting dry matter to the active growing tissues in the plant each day, transpiration and the uptake of nitrogen from the soil profile are simulated. It incorporates the RHIZOS model which simulates, in two dimensions, the movement of water, roots, and soluble nutrients through the soil profile. It records the time of initiation of each of the plant organs. These phenological events are calculated from temperature functions with delays resulting from physiological stress. Stress is defined mathematically as an imbalance in the metabolite supply; demand ratio. Physiological stress is also the basis for the calculation of rates of tiller and floret abortion. Thus, tillering and head differentiation are modeled as the resulants of the two processes, morphogenesis and abortion, which may be occurring simulaneously.

Baker, D. N.↗

Height and Biomass of Mangroves in Africa from ICEsat/GLAS and SRTM

The accurate quantification of forest 3-D structure is of great importance for studies of the global carbon cycle and biodiversity. These studies are especially relevant in Africa, where deforestation rates are high and the lack of background data is great. Mangrove forests are ecologically significant and it is important to measure mangrove canopy heights and biomass. The objectives of this study are to estimate: 1. The total area, 2. Canopy height distributions and 3. Aboveground biomass of mangrove forests in Africa. To derive mangrove 3-D structure and biomass maps, we used a combination of mangrove maps derived from Landsat ETM+, LiDAR canopy height estimates from ICEsat/GLAS (Ice, Cloud, and land Elevation Satellite/Geoscience Laser Altimeter System) and elevation data from SRTM (Shuttle Radar Topography Mission) for the African continent. More specifically, we extracted mangrove forest areas on the SRTM DEM using Landsat based landcover maps. The LiDAR (Light Detection and Ranging) measurements from the large footprint GLAS sensor were used to derive local estimates of canopy height and calibrate the Interferometric Synthetic Aperture Radar (InSAR) data from SRTM. We then applied allometric equations relating canopy height to biomass in order to estimate above ground biomass (AGB) from the canopy height product. The total mangrove area of Africa was estimated to be 25 960 square kilometers with 83% accuracy. The largest mangrove areas and greatest total biomass was 29 found in Nigeria covering 8 573 km2 with 132 x10(exp 6) Mg AGB. Canopy height across Africa was estimated with an overall root mean square error of 3.55 m. This error also includes the impact of using sensors with different resolutions and geolocation error which make comparison between measurements sensitive to canopy heterogeneities. This study provides the first systematic estimates of mangrove area, height and biomass in Africa. Our results showed that the combination of ICEsat/GLAS and SRTM data is well suited for vegetation 3-D mapping on a continental scale.

Fatoyinbo, Temilola E.↗

The Ocean Colour Climate Change Initiative: I. A Methodology for Assessing Atmospheric Correction Processors Based on In-Situ Measurements

The Ocean Colour Climate Change Initiative intends to provide a long-term time series of ocean colour data and investigate the detectable climate impact. A reliable and stable atmospheric correction procedure is the basis for ocean colour products of the necessary high quality. In order to guarantee an objective selection from a set of four atmospheric correction processors, the common validation strategy of comparisons between in-situ and satellite derived water leaving reflectance spectra, is extended by a ranking system. In principle, the statistical parameters such as root mean square error, bias, etc. and measures of goodness of fit, are transformed into relative scores, which evaluate the relationship of quality dependent on the algorithms under study. The sensitivity of these scores to the selected database has been assessed by a bootstrapping exercise, which allows identification of the uncertainty in the scoring results. Although the presented methodology is intended to be used in an algorithm selection process, this paper focusses on the scope of the methodology rather than the properties of the individual processors.

Climate Change↗

Triacylglycerol, total fatty acid, and biomass accumulation of metabolically engineered energycane grown under field conditions confirms its potential as feedstock for drop‐in fuel production

Abstract Metabolic engineering for hyperaccumulation of lipids in vegetative tissues of high biomass crops promises a step change in oil yields for the production of advanced biofuels. Energycane is the ideal feedstock for this approach due to its exceptional biomass production and persistence under marginal conditions. Here, we evaluated metabolically engineered energycane with constitutive expression of the lipogenic factors WRINKLED 1 ( WRI 1), DIACYLGLYCEROL ACYLTRANSFERASE 1 ( DGAT 1), and OLEOSIN 1 ( OLE 1) for the accumulation of triacylglycerol (TAG), total fatty acid (TFA), and biomass under field conditions at the University of Florida‐IFAS experiment station near Citra, Florida. TAG and TFA accumulation were highest in leaves (up to 9.9% and 12.9% of DW, respectively), followed by juice from crushed stems, stems, and roots. TAG and TFA accumulation increased up to harvest time and correlated highest with OLE 1 and DGAT 1 expression. Biomass dry weight, TAG, and TFA content differed greatly depending on DGAT 1 and OLE 1 expression in transgenic lines with similar WRI 1 expression. Biomass did not significantly differ between WT and line L2 with DAGT 1 and OLE 1 expressed at low levels and TAG and TFA accumulating to 12‐ and 1.6‐fold that of WT leaves, respectively. In contrast, line L13, with intron‐mediated enhancement of DGAT 1 expression, displayed a 245‐ to 330‐fold increase in TAG and a 4.75‐ to 6.45‐fold increase in TFA content compared with WT leaves and a biomass reduction of 52%. These results provide the basis for developing novel feedstocks for expanding plant lipid production and point to new prospects for advanced biofuels.

59 BASIC BIOLOGICAL SCIENCES↗

A methodology for decay heat characterization in molten salt reactors

Accurate decay heat prediction in molten salt reactors (MSRs) faces dual challenges: complex operational uncertainties and the need for interpretable models compatible with engineering workflows. This work presents a hybrid machine learning and segmented polynomial methodology that addresses both requirements through three key innovations. First, a modular data architecture encodes MSR-specific operational parameters (power density: 1-100 W cm -3 , humidity: 0-0.1 wt %, air ingress: 0-0.1 mol %) with uncertainty-aware temporal discretization spanning 15 orders of magnitude. Second, region-optimized machine learning models achieve 92.3 % root mean square error (RMSE) reduction over conventional polynomials while maintaining physical interpretability through automated piecewise equation generation. Third, dual front-end interfaces accelerate safety analyses — a Jupyter environment enables researchers to explore 10,000+ parameter combinations via interactive widgets, while a Streamlit web application reduces design iteration cycles through production-grade visualization tools. Operational deployment demonstrates prediction times of only a couple hundred milliseconds for 10 4 years decay profiles, enabling real-time optimization of spent fuel container designs.

42 - ENGINEERING↗

On the Relationship Between Methane Production in Anaerobic Incubations of Peat Material and In Situ Methane Emissions

Anaerobic incubations of peat have been widely used to explore soil processes, but this in vitro technique raises many questions as to how well it reproduces in situ conditions. To investigate this, we conducted 60–100 days (+25 days pre‐incubation) anaerobic, temperature‐controlled incubation experiments across a temperature range of 1–26°C on samples from bog and fen habitats, at two different depths (9–19 and 25–35 cm). Here, we observed exponential increases in CO 2 and methane production with temperature in all conditions. We then compared field‐based measurements of methane emission with modeled expectations by extrapolating incubation‐determined methane production rates based on (a) soil temperature profiles, (b) the observed incubation temperature‐methane production relationship, and (c) seasonal thaw depth from each site. The resulting incubation‐extrapolated methane production agreed with measured emission rates within a factor of two at both sites and corresponded to 182 ± 54% and 59 ± 14% of the measured average yearly fluxes from the field for the bog and fen, respectively. The underestimation of fen methane fluxes may be due to the lack of living plant root‐derived dissolved organic carbon inputs in incubations, a key process in fens. Conversely, the overestimation in bogs could be attributed to methane oxidation in the field, which is absent in anaerobic incubation conditions. Nonetheless incubations predicted greenhouse gas emissions from a northern peatland within a factor of two.

Biological and medical sciences,Geosciences,Enviro↗

A critical analysis of the accuracy of several numerical techniques for combustion kinetic rate equations

A detailed analysis of the accuracy of several techniques recently developed for integrating stiff ordinary differential equations is presented. The techniques include two general-purpose codes EPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREK1D, and GCKP4 developed specifically to solve chemical kinetic rate equations. The accuracy study is made by application of these codes to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas-phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. To illustrate the error variation in the different combustion regimes the species are divided into three types (reactants, intermediates, and products), and error versus time plots are presented for each species type and the temperature. These plots show that CHEMEQ is the most accurate code during induction and early heat release. During late heat release and equilibration, however, the other codes are more accurate. A single global quantity, a mean integrated root-mean-square error, that measures the average error incurred in solving the complete problem is used to compare the accuracy of the codes. Among the codes examined, LSODE is the most accurate for solving chemical kinetics problems. It is also the most efficient code, in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that use of the algebraic enthalpy conservation equation to compute the temperature can be more accurate and efficient than integrating the temperature differential equation.

Radhadrishnan, Krishnan↗