Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnitude”

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 919 records · Page 51

AFRL Additive Manufacturing Modeling Series: Challenge 1, Characterization of Residual Strain Distribution in Additively-Manufactured Metal Parts Using Energy-Dispersive Diffraction

A machine component fabricated by additive manufacturing (AM) of metal powders often possesses steep residual stress gradients with significant magnitudes due to large temperature gradients that transpire within a localized area during fabrication. These processing-induced residual stresses can cause distortion of the part, and if they are of significant magnitude, they could induce cracking of the component, degrade printability, and/or diminish subsequent mechanical performance. The ability to predict these residual stresses imparted by AM is an important step in permeating AM technology for advanced manufacturing. Calibration and validation of AM process models used for prediction are, therefore, a critical step in understanding the origin and mitigating the challenges associated with residual stresses inherent to the AM process. In the present work, the residual strain distributions in components with simple geometries fabricated by a laser powder bed fusion (LBPF) process were characterized non-destructively using energy-dispersive X-ray diffraction in support of the US Air Force Research Laboratory Additive Manufacturing Modeling Challenge Series Groeber et al. (JOM 70:441-448, 2018). The measurement setup and approach are described in detail so that the data can be used as a benchmark to calibrate and validate models for the prediction of macro-scale residual stresses due to the LPBF process.

36 MATERIALS SCIENCE↗

Assessment of CMIP6 Performance and Projected Temperature and Precipitation Changes Over South America

We evaluate the performance of a large ensemble of Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6) over South America for a recent past reference period and examine their projections of twenty-first century precipitation and temperature changes. The future changes are computed for two time slices (2040–2059 and 2080–2099) relative to the reference period (1995–2014) under four Shared Socioeconomic Pathways (SSPs, SSP1–2.6, SSP2–4.5, SSP3–7.0 and SSP5–8.5). The CMIP6 GCMs successfully capture the main climate characteristics across South America. However, they exhibit varying skill in the spatiotemporal distribution of precipitation and temperature at the sub-regional scale, particularly over high latitudes and altitudes. Future precipitation exhibits a decrease over the east of the northern Andes in tropical South America and the southern Andes in Chile and Amazonia, and an increase over southeastern South America and the northern Andes—a result generally consistent with earlier CMIP (3 and 5) projections. However, most of these changes remain within the range of variability of the reference period. In contrast, temperature increases are robust in terms of magnitude even under the SSP1–2.6. Future changes mostly progress monotonically from the weakest to the strongest forcing scenario, and from the mid-century to late-century projection period. There is an increase in the seasonality of the intra-annual precipitation distribution, as the wetter part of the year contributes relatively more to the annual total. Furthermore, an increasingly heavy-tailed precipitation distribution and a rightward shifted temperature distribution provide strong indications of a more intense hydrological cycle as greenhouse gas emissions increase. The relative distance of an individual GCM from the ensemble mean does not substantially vary across different scenarios. We found no clear systematic linkage between model spread about the mean in the reference period and the magnitude of simulated sub-regional climate change in the future period. Overall, these results could be useful for regional climate change impact assessments across South America.

54 ENVIRONMENTAL SCIENCES↗

Crystal chemistry at high pressure

The chemistry we are taught in school, and we experience in our daily existence occurs at 1 atm. However, pressure spans an astounding 62 orders of magnitude in the Universe in going from the void of interstellar space to the crushing conditions at the center of a neutron star. The way in which pressure affects chemistry is important for Earth and planetary sciences, materials science, in understanding the extreme conditions experienced in nuclear explosions, and it may be key in addressing the future energy needs of our society. In this article we outline how the often neglected pressure variable affects chemistry, beginning from the way in which atomic energy levels are altered. This lays the foundation for understanding the unique crystal and electronic structures that emerge when matter is squeezed, as well as pressure’s effect on chemical reactivity. Finally, we give an overview of the main concepts behind conventional, or phonon-mediated, superconductivity, and describe how the pressure variable may be key in discovering and designing light-element based materials whose superconducting critical temperatures approach room temperature. Here, we discuss some of the main families of superconducting hydrides that have been predicted computationally, and the experimental successes in this exciting and rapidly developing field.

Chemical bonding↗

Distorted dislocation cores and asymmetric glide resistances in titanium

The determination of Critical Resolved Shear Stress (CRSS) in titanium for basal, prismatic, and pyramidal slip-planes without empirical constants is presented by combining Density Functional Theory (DFT) and anisotropic elasticity. A new mechanism leading to tension-compression (T-C) asymmetry in the CRSS levels has been revealed for the first time. The conditions for this asymmetry are established, involving a complex interplay between the dislocation core-structure and core-advance behavior. The three conditions for T-C asymmetry that need to be simultaneously satisfied can be summarized as: (1) a medium stacking fault width, d, (3 < d/b F < 10, where b F is the magnitude of the Burgers vector of the full dislocation), (2) an asymmetric core-structure of the extended dislocation (ξ 1 ≠ ξ 2 , where ξ 1 and ξ 2 are the core-widths of the first and second partials, respectively), and (3) intermittent motion of the partials (Δd/b F ≠ 0, where Δd is the magnitude of fluctuation in stacking fault width during intermittent motion). Pyramidal-slip in titanium satisfies all three conditions, resulting in significant T-C asymmetry. The CRSS theory considers a Wigner-Seitz (W-S) cell based domain area assigned to each lattice site for the calculations of core-energies accurately capturing the slip-plane lattice. The W-S based approach is essential due to the lower symmetry of the HCP crystal circumventing potential errors associated with the one-dimensional atomic-row approximation. Dislocation core structures are obtained for all the slip-systems in titanium showing significant distortion of the disregistry profile governing the core-advance behavior. The CRSS values predicted from the theory show agreement with the experimental CRSS levels reported in the literature.

Dislocation core↗

Strain-rate hardening enhances fatigue resistance of AlSi10Mg alloy at 350°C

The high cycle fatigue behavior of laser powder bed fusion processed AlSi10Mg has been investigated at 350 °C (T/T m ∼ 0.7). The alloy exhibited a fatigue strength of 30 MPa defined by runout after 10 7 cycles at a conventional loading frequency of 20 Hz, corresponding to a notable fatigue strength to ultimate tensile strength ratio of 0.73. The surprising fatigue resistance was attributed to the strain-rate hardening effect at high fatigue loading frequency relative to tensile loading rates at 350 °C. The strain-rate hardening effect was validated by performing ultrasonic fatigue tests (20 kHz loading frequency) with three orders of magnitude higher strain-rates than those at conventional loading frequency. The higher strain-rates in ultrasonic fatigue increased the magnitude of strain-rate hardening resulting in longer AlSi10Mg fatigue lives compared to fatigue at conventional frequency, thus confirming the strain-rate hardening effect. The fatigue crack initiation mechanism was strain-rate dependent. Post-mortem microstructural examination revealed intergranular cavitation inside clusters of fine equiaxed grains. The cavities interlinked with each other to initiate near-surface fatigue cracks at conventional frequency. Cavitation occurred to a lesser extent at the ultrasonic frequency. As a result, fatigue cracks initiated at pre-existing processing defects near the surface in ultrasonic fatigue samples. Finally, this investigation underscores the role of fine grain clusters in promoting high-temperature fatigue crack initiation and indicates a possible trade-off between printability via grain refinement and high-temperature fatigue resistance of additively manufactured alloys.

Additive manufacturing↗

Which way does the dendrite grow? Competition among epitaxy, preferred growth direction, and thermal gradients in powder bed fusion additive manufacturing

The as-processed microstructure of metal alloy parts manufactured through laser powder bed fusion (LPBF) is heavily derived from the cellular dendritic solidification. The growth direction of dendrites within the melt pool is determined through competition among epitaxial growth, preferred growth directions, and maximum thermal gradients. However, the dominant factor and the specific role of each in developing melt pool microstructures remain unknown. Here, in this study, we performed single laser track scans on an SS316L single crystal substrate and combined experimental characterization of microstructure and crystal orientations with Computational Fluid Dynamics simulations of thermal gradients to evaluate the role of each factor in determining dendritic growth direction and evolution. Our results reveal that epitaxial growth dominates microstructure development by preferentially growing along a single 〈100〉 variant of the single crystal substrate adjacent to the melt pool boundary. Under LPBF’s highly curved and rapidly evolving thermal field, this preferential dendrite variant selection and its continued growth from the melt pool boundary to the centerline are governed by the local temperature gradient magnitude at the solid-liquid interface, rather than by the instantaneous maximum temperature gradient direction alone. Using these findings, we successfully predict changes in the dendrite growth direction with changing laser scan direction on a single crystal substrate, and show that the geometric melt pool centerline can deviate from the microstructural centerline because asymmetric local temperature gradient magnitudes transiently limit growth, resulting in different dendrite travel distances on each side of the melt pool.

36 MATERIALS SCIENCE↗

Quantifying carbon budget, crop yields and their responses to environmental variability using the ecosys model for U.S. Midwestern agroecosystems

As one of the major agricultural production areas in the world, the United States (U.S.) Midwest plays a vital role in the global food supply and agricultural ecosystem services. Although significant efforts have been made in modeling the carbon cycle dynamics over this area, large uncertainty still exists in the previous simulations in terms of reproducing individual components of the carbon cycle and their responses to environmental variability. Here we evaluated the performance of an advanced agroecosystem model, ecosys, in simulating carbon budgets over the U.S. Midwest, considering both the magnitude of carbon flux/yield and its response to environmental (climate and soil) variability. We conducted model simulations and evaluations at 7 cropland eddy-covariance sites as well as over 293 counties of Illinois, Indiana, and Iowa in the U.S. Midwest. The site-level simulations showed that ecosys captured both the magnitude and seasonal patterns of carbon fluxes (i.e., net ecosystem carbon exchange (NEE), ecosystem gross primary production (GPP), and ecosystem respiration (Reco)), leaf area index (LAI), and dynamic plant carbon allocation processes, with R2 equal to 0.92, 0.87, 0.87, and 0.78 for GPP, NEE, Reco, and LAI, respectively across all the sites compared with the observations. For regional scale simulations, ecosys reproduced the spatial distribution and interannual variability of corn and soybean yields with the constraints of observed yields and a new remotely sensed GPP product, with R2 of multi-year averaged simulated and observed yield equal 0.83 and 0.80 for corn and soybean, respectively. The simulated responses of carbon cycle dynamics to environmental variability were consistent with that from the empirical observations at both site and regional scales. Furthermore, our results demonstrated the applicability of ecosys in simulating the carbon cycle and soil carbon dynamics of the U.S. Midwestern agroecosystems under different climate and soil conditions.

54 ENVIRONMENTAL SCIENCES↗

Singular value decomposition of adjoint flux distributions for Monte Carlo variance reduction

Monte Carlo (MC) shielding calculations often use weight windows (WWs) and biased sources formed from a deterministic estimate of the adjoint flux to improve the convergence rate of tallies. This requires a significant amount of computer memory, which can limit the memory available for high-resolution tally output. A new method is proposed for reducing these memory requirements by using singular value decomposition (SVD) in linear or logarithmic space to approximate the adjoint flux. This method’s performance is evaluated using the Shift and Denovo codes for streaming and diffusion base case problems, followed by problems using the Westinghouse AP1000 and the Joint European Torus. The log SVD reduced WW memory requirements by an order of magnitude in all cases without a significant performance penalty. Additionally, the linear SVD reduced biased source memory requirements by an order of magnitude, but further investigation is needed to account for observed limitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Reactor and fuel cycle performance of light water reactor fuel with 235 U enrichments above 5%

Recent advances in nuclear fuel materials research, particularly on the topic of accident-tolerant fuels, have brought up potential opportunities for expanding the operating envelope of existing light water reactors. As many of the performance improvements offered by these technologies may be most fully realized by increasing fuel enrichment beyond the standard 5% limit, this paper examines the potential reactor performance and fuel cycle performance of low-enriched uranium oxide fueled light water reactors by generically considering pressurized water reactors with 235 U enrichment from 5 to 7%. Advanced cladding, including accident-tolerant cladding, has the potential to increase fuel burnup limits related to hydrogen in the cladding that coincide with those limits associated with end-of-life reactivity. Therefore, higher enrichment will be necessary in order to realize the higher fuel burnups. This work includes evaluation of the fuel cycle length, discharge burnup, reactivity coefficients, and fuel cycle performance, including radioactive waste and environmental impact metrics per unit energy generated. The analysis was performed using the evaluation metrics from the US Department of Energy Office of Nuclear Energy Fuel Cycle Evaluation and Screening Study. The reactor performance and safety analysis show that enrichments between 5 and 7% would have similar fuel temperature and moderator temperature coefficients. However, the soluble boron coefficient would decrease in magnitude, requiring more corrosive boric acid in the coolant or other methods of reactivity control during the fuel cycle. At these higher enrichments the maximum burnup at the rim of the fuel pellet would increase by almost a factor of two, which is expected to impact the formation of high-burnup structure in the fuel and the corresponding thermo-mechanical fuel properties. The fuel cycle performance assessment shows that increasing enrichment reduces the quantity of high-level waste disposed per unit energy generated, but it increases the natural resource requirements normalized to a gigawatt-electricity-per-year basis. Another impact is the slightly higher discharge burnup, resulting in somewhat different activity levels of the spent nuclear fuel and high-level waste radioactivity at 100 and 100,000 years after fuel discharge. The environmental impacts—including land use, water use, carbon emission, and radiological exposure—are of the same magnitude per unit energy generated. Finally, however, the impacts are distributed differently. Less than 5% enrichment has marginally more impact on the back-end of the fuel cycle, and greater than 5% enrichment has marginally more impact on the front-end of the fuel cycle. Ultimately, no neutronic or reactor safety hindrances to employing light water reactor fuel with enrichments greater than 5% are identified; given the achievable reactor performance benefits with advanced fuels, further practical exploration of increased enrichment fuel is recommended.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Accelerating the solution of the S N equations with highly anisotropic scattering using the Fokker-Planck approximation

The discrete ordinates method can model forward-peaked transport problems accurately. However, convergence of discrete ordinates solution can become arbitrarily slow upon use of standard iterative procedures like source iteration. Standard zeroth and first moment-based acceleration methods like nonlinear diffusion acceleration and diffusion synthetic acceleration are ineffective in accelerating such problems because these methods do not correct higher order Legendre-moments of angular flux. In this paper, we explore the idea of using the Fokker-Planck approximation as a preconditioner to accelerate the convergence of solution of forward-peaked transport problems. Finally, we observe that such preconditioning can successfully reduce the iteration count by up to four orders of magnitude, and solver runtime by up to three orders of magnitude when compared to diffusion synthetic acceleration.

42 ENGINEERING↗

Assessment of tritium effluent from Prototype Generation IV Sodium-cooled Fast Reactor

We report a tritium transport model for a pool-type sodium cooled fast reactor was developed for estimation of tritium effluents from the Prototype Generation IV Sodium-cooled Fast Reactor (PGSFR). For this purpose, the transport model developed for Experimental Breeder Reactor II and Fast Flux Test Facility was updated and validated. The tritium production rate from ternary fission and activations in PGSFR were calculated by the Argonne Reactor Computation (ARC) code suite for fast reactor analysis, and due to the potential uncertainties in source terms and variations of cold trap performance and material permeability during reactor operation, the tritium distributions and tritium effluents were estimated for several bounding cases. The upper bounding gaseous and liquid tritium release rates from PGSFR are similar to ~440.7 Ci/year and 14.7 Ci/year, respectively. The gaseous tritium effluent is about few orders of magnitude smaller than the regulatory constraint of the tritium concentration in the air. The peak tritium concentration in the steam-generator feedwater is similar to ~1.1x10 -3 mu Ci/g, which is required to be diluted prior to releasing as tritiated water. The accumulated tritium in the cold trap after five cycle operations is ~ 52,000 Ci. Because of significant conservatism in these upper bounding estimations, it is expected that the tritium effluent at the nominal case with plausible source terms would be about 2 orders of magnitude smaller than that of the upper bounding case.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An uncertainty quantification method relevant to material test reactors

Within material test reactor calculations, energy dependent flux and reaction rate uncertainties are typically not quantified when performing as-run analyses to determine the neutron field experienced by the experiment. When high fidelity Monte-Carlo codes are used in such analyses, straight forward methods to calculate output uncertainties are not available, instead expert opinion is used to postulate computational uncertainties. New methods to propagate uncertainties through these high fidelity simulations are available when sufficient computational power is available. A tool is developed here for sampling any part of an MCNP input from random distributions to determine output uncertainties based on those inputs. Another tool is developed to sample nuclear data cross-section in ACE format using multi-group nuclear data covariances. The Total Monte-Carlo Method and Gesellschaft für Anlagen-und Reaktorsicherheit method (GRS) are implemented and compared to one another as well as MCNP sensitivity and uncertainty calculations. The methods were applied to the Godiva critical sphere k-eigenvalue, the UAM pin-cell benchmark energy dependent flux and reaction rates, and the Advanced Test Reactor energy dependent flux within an experimental location. Furthermore, the two methods agree well, with GRS allowing for an order of magnitude speedup for reaction rate uncertainty calculations and several orders of magnitude for eigenvalue uncertainty calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Numerical simulation of involute-plate research reactor flow behavior using RANS, LES and DNS

This paper investigates the flow behavior of involute-plate research reactors by performing Reynolds-Averaged Navier Stokes simulation (RANS), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) of the channel flow between fuel plates. By modeling turbulence with different numerical approaches, this study provides data with three levels of fidelity. For the RANS simulation, three widely used turbulence models, i.e., k-ε, k-ω, Reynolds Stress Turbulence model (RST) are applied by using the commercial CFD code STAR-CCM +. For LES and DNS, the open-source CFD code, Nek5000, is used given its outstanding scalability on High Performance Computer (HPC) and high-order technique. The results from RANS simulations are compared with that from LES and DNS for benchmarking. Both macroscale parameters and turbulence statistics, such as velocity magnitude, lateral velocity and turbulence kinetic energy, are presented and analyzed. The results from RANS simulation achieve good agreement with LES and DNS on velocity and turbulence kinetic energy prediction. The RST turbulence model predicts the most similar flow pattern of lateral velocity as compared to LES and DNS. The Lambda-2 (λ2) criterion with a reasonable threshold is used to demonstrate the instantaneous vortices distribution in the involute channel from both LES and DNS calculation. The DNS simulation captures more detailed turbulence especially near the corner, which explains the discrepancy between LES and DNS results near the corner. The normalized RMS error are defined and calculated to assess the performance of those turbulence models. The RST model captures the anisotropic feature of turbulence, which enable it to outperform other turbulence models for predicting the flow behavior in an involute channel. Although some discrepancies are found between LES and DNS results in the corner, the overall deviations between LES and DNS are found to be small. In conclusion, given that the computational cost of DNS calculation is an order of magnitude higher, using LES data for benchmarking RANS model is a cost-effective approach.

DNS↗

Impact of size and thermal gradient on supercooling of phase change materials for thermal energy storage

Phase change material based thermal energy storage has many current and potential applications in the heating and cooling of buildings, battery and electronics thermal management, thermal textiles, and dry cooling of power plants. However, connecting lab scale thermal data obtained with differential scanning calorimetry (DSC) to the performance of large-scale practical systems has been a major challenge primarily due to the dependence of supercooling on the size and temperature gradient of the system. Here we show how a phase change material’s supercooling behavior can be characterized experimentally using common lab scale thermal analysis techniques. We then develop a statistics based theoretical model that uses the lab-scale data on small samples to quantitatively predict the supercooling performance for a general thermal energy storage application of any size, including also allowing for the possibility of temperature gradients. Finally, we validate the modeling methodology by comparing to experimental results for solid-solid phase change in neopentyl glycol, which shows how the model successfully predicts the changes in supercooling temperature across a large range of cooling rates (2 orders of magnitude) and volumes (3 orders of magnitude). By accounting for thermal gradients, the model avoids ~2x error incurred by lumped approximations.

36 MATERIALS SCIENCE↗

Strong temporal variability in methane fluxes from natural gas well pad soils

We measured methane and carbon dioxide fluxes at natural gas well pad soils and undisturbed soils in the Rocky Mountain and Gulf Coast regions of the United States, including producing and gas storage wells. We collected both short-term (15 min) and multi-day (between 3 and 8), continuous measurements at 47 well pads and two undisturbed locations. Methane fluxes varied by more than an order of magnitude over periods as short as 30 min (e.g., 19–593 mg m -2 h -1 in one instance), and diurnal and seasonal variability was also significant (e.g., spring-to-fall change from 509 to 14174 mg m -2 h -1 ). We hypothesize that short-term flux variability was caused by pulsed flow of methane during its migration through the subsurface. Barometric pressure and well conditions likely impacted fluxes, but we found only weak evidence for this. Bacterial methanotrophy appeared to impact methane flux magnitude and variability. We injected methane into the subsurface at one well, and we found that, while fluxes of methane and carbon dioxide, and combustible soil gas concentrations, increased in response to the injection, the response was not uniform, and fluxes exhibited high hourly-scale variability, in spite of a constant injection rate. Methane fluxes tended to be higher at well pad soils compared to background soils (often much higher), and fluxes tended to be higher at well pad locations closer to the well head.

03 NATURAL GAS↗

Photometric redshifts probability density estimation from recurrent neural networks in the DECam local volume exploration survey data release 2

Photometric wide-field surveys are imaging the sky in unprecedented detail. These surveys face a significant challenge in efficiently estimating galactic photometric redshifts while accurately quantifying associated uncertainties. In this work, we address this challenge by exploring the estimation of Probability Density Functions (PDFs) for the photometric redshifts of galaxies across a vast area of 17,000 square degrees, encompassing objects with a median 5 σ point-source depth of g = 24.3, r = 23 . 9 , i = 23.5, and z = 22.8 mag. Our approach uses deep learning, specifically integrating a Recurrent Neural Network architecture with a Mixture Density Network, to leverage magnitudes and colors as input features for constructing photometric redshift PDFs across the whole DECam Local Volume Exploration (DELVE) survey sky footprint. Subsequently, we rigorously evaluate the reliability and robustness of our estimation methodology, gauging its performance against other well-established machine learning methods to ensure the quality of our redshift estimations. Our best results constrain photometric redshifts with the bias of − 0 . 0013 , a scatter of 0.0293, and an outlier fraction of 5.1%. These point estimates are accompanied by well-calibrated PDFs evaluated using diagnostic tools such as Probability Integral Transform and Odds distribution. We also address the problem of the accessibility of PDFs in terms of disk space storage and the time demand required to generate their corresponding parameters.We present a novel Autoencoder model that reduces the size of PDF parameter arrays to one-sixth of their original length, significantly decreasing the time required for PDF generation to one-eighth of the time needed when generating PDFs directly from the magnitudes.

79 ASTRONOMY AND ASTROPHYSICS↗

Survey-wide asteroid discovery with a high-performance computing enabled non-linear digital tracking framework

Modern astronomical surveys detect asteroids by linking together their appearances across multiple images taken over time. This approach faces limitations in detecting faint asteroids and handling the computational complexity of trajectory linking. Here, we present a novel method that adapts “digital tracking” – traditionally used for short-term linear asteroid motion across images – to work with large-scale synoptic surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time (Rubin/LSST). Our approach combines hundreds of sparse observations of individual asteroids across their non-linear orbital paths to enhance detection sensitivity by several magnitudes. To address the computational challenges of processing massive data sets and dense orbital phase spaces, we developed a specialized high-performance computing architecture. We demonstrate the effectiveness of our method through experiments that take advantage of the extensive computational resources at Lawrence Livermore National Laboratory. This work enables the detection of significantly fainter asteroids in existing and future survey data, potentially increasing the observable asteroid population by orders of magnitude across different orbital families, from near-Earth objects (NEOs) to Kuiper belt objects (KBOs).

Asteroid discovery↗

Turbulent transport and reactions of plant-emitted hydrocarbons in an Amazonian rain forest

The processes governing the temporal and spatial patterns of isoprene and monoterpenes emitted by a rainforest in the central Amazon region of Brazil is investigated using a combination of field experiments and numerical simulations. Specifically, Large Eddy Simulations (LES) were used to resolve emissions of isoprene and monoterpenes, turbulent transport, and air chemistry. The coupled chemistry-transport LES included the effects of isoprene and monoterpenes reactivity due to reactions with hydroxyl radical and ozone. The LES results are used to compute vertically resolved budgets of isoprene and monoterpenes in the rainforest canopy in response to emissions, turbulent transport, surface deposition, and air chemistry. Our results indicated that emission and dispersion dominated the isoprene budget as the gases were transported out of the canopy space. In a region limited by nitrogen oxides (with prevailing nitric oxide levels of < 0.5 parts per billion), the in-canopy chemical destruction removed approximately 10% of locally emitted monoterpenes. Hydroxyl radical production rates from the ozonolysis of monoterpenes amounted to ≈ 2 x 10 6 radicals cm -3 s -1 and had similar magnitude to the light-dependent hydroxyl radical formation. One key conclusion was that the Amazonia rainforest abundantly emitted monoterpenes whose in-canopy ozonolysis yielded hydroxyl radicals in amounts similar to the magnitude of light-dependent formation. Reactions of monoterpenes and isoprene with hydroxyl radical and ozone were necessary for the maintenance of the Amazon rainforest canopy as a photochemically active environment suitable to generate oxidants and secondary organic aerosols.

54 ENVIRONMENTAL SCIENCES↗