Search NASA⌕ Search

SEARCH · Search NASA

Results for “shape optimization”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Compression dynamics of quasi-spherical wire arrays with different linear mass profiles

Results of experimental studies of the implosion of quasi-spherical wire (or metalized fiber) arrays are presented. The goal of the experiments was to achieve synchronous three-dimensional compression of the plasma produced in different regions of a quasi-spherical array into its geometrical center. To search for optimal synchronization conditions, quasi-spherical arrays with different initial profiles of the linear mass were used. The following dependences of the linear mass on the poloidal angle were used: m{sub l}(θ) ∝ sin{sup –1}θ and m{sub l}(θ) ∝ sin{sup –2}θ. The compression dynamics of such arrays was compared with that of quasi-spherical arrays without linear mass profiling, m{sub l}(θ) = const. To verify the experimental data, the spatiotemporal dynamics of plasma compression in quasi-spherical arrays was studied using various diagnostics. The experiments on three-dimensional implosion of quasi-spherical arrays made it possible to study how the frozen-in magnetic field of the discharge current penetrates into the array. By measuring the magnetic field in the plasma of a quasi-spherical array, information is obtained on the processes of plasma production and formation of plasma flows from the wire/fiber regions with and without an additionally deposited mass. It is found that penetration of the magnetic flux depends on the initial linear mass profile m{sub l}(θ) of the quasi-spherical array. From space-resolved spectral measurements and frame imaging of plasma X-ray emission, information is obtained on the dimensions and shape of the X-ray source formed during the implosion of a quasi-spherical array. The intensity of this source is estimated and compared with that of the Z-pinch formed during the implosion of a cylindrical array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Solidification and crystallographic texture modeling of laser powder bed fusion Ti-6Al-4V using finite difference-monte carlo method

Laser powder bed fusion (LPBF) additive manufacturing makes near-net-shaped parts with reduced material cost and time, rising as a promising technology to fabricate Ti-6Al-4V, a widely used titanium alloy in aerospace and medical industries. However, LPBF Ti-6Al-4V parts produced with 67° rotation between layers, a scan strategy commonly used to reduce microstructure and property inhomogeneity, have varying grain morphologies and weak crystallographic textures that change depending on processing parameters. Here, this study predicts LPBF Ti-6Al-4V solidification at three energy levels using a finite difference-Monte Carlo method and validates the simulations with large-area electron backscatter diffraction (EBSD) scans. The developed model accurately shows that a <001> texture forms at low energy and a <111> texture occurs at higher energies parallel to the build direction but with a lower strength than the textures observed from EBSD. A validated and well-established method of combining spatial correlation and general spherical harmonics representation of texture is developed to calculate a difference score between simulations and experiments. The quantitative comparison enables effective fine-tuning of nucleation density (N 0 ) input, which shows a nonlinear relationship with increasing energy level. Future improvements in texture prediction code and a more comprehensive study of N 0 with different energy levels will further advance the optimization of LPBF Ti-6Al-4V components. These developments contribute a novel understanding of crystallographic texture formation in LPBF Ti-6Al-4V, the development of robust model validation and calibration pipeline methodologies, and provide a platform for mechanical property prediction and process parameter optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Partnership for Edge Physics Simulation

A major goal of our participation in the Edge Physics Simulation project has been to contribute to the understanding of the self-organization of tokamak turbulence fluctuations resulting in the formation of a staircase structure in the ion temperature. A second important goal is to demonstrate how small scale turbulence in plasmas self-organizes with dynamically driven quasi-stationary flow shear. These goals have been accomplished through the analyses of the statistical properties of XGC1 flux driven Gyrokinetic electrostatic ion temperature gradient (ITG) turbulence simulation data in which neutrals are included. The ITG turbulence data, and in particular fluctuation data, were obtained from a massively parallel flux-driven gyrokinetic full-f particle-in-cell simulation of a DIII-D like equilibrium. Below some the findings are summarized. It was observed that the emergence of staircase structure is related to the variations in the normalized temperature gradient length (R/LT) and the poloidal flow shear. Average turbulence intensity is found to be large in the vicinity of minima in R/LTi, where ITG growth is expected to be lower. The distributions of the occurrences of potential fluctuation are found to be Gaussian away from the staircase-step locations, but they are found to be non-Gaussian in the vicinity of staircase-step locations. The results of analytically derived expressions for the distribution of the occurrences of turbulence intensity and intensity flux were compared with the corresponding quantities computed using XGC1 simulation data and good agreement is found. The derived expressions predicts inward and outward propagation of turbulence intensity flux in an intermittent fashion. The outward propagation of turbulence intensity flux occurs at staircase-step locations and is related to the change in poloidal flow velocity shear and to the change in the ion temperature gradient. The standard deviation, skewness and kurtosis for turbulence quantities were computed and found to be large in the vicinity of the staircase-step structures. Large values of skewness and kurtosis can be explained by a temporary opening and closing of the structure which allows turbulence intensity events to propagate. The staircase patterns may reduce the ion heat transport and a manipulation of these patterns may be used to optimize heat transport in tokamaks. An additional objective of the research in support of the Edge Physics Simulation initiative has been to improve the understanding of scrape-off layer thermal transport. In planning experiments and designing future tokamaks, it is important to understand the physical effects that contribute to divertor heat-load fluxes. The research accomplished will contribute to developing new models for the scrape-off layer region. The XGC0 code was used to compute the heat fluxes and the heat-load width in the outer divertor plates of C-Mod and DIII-D tokamaks. It was observed that the width of the XGC0 neoclassical heat-load was approximately inversely proportional to the total plasma current. Anomalous transport in the H-mode pedestal region of five Alcator C-Mod discharges, representing a collisionality scan, was analyzed. The understanding of anomalous transport in the pedestal region is important for the development of a comprehensive model for the H-mode pedestal slope. It was found that the electron thermal anomalous diffusivities at the pedestal top increase with the electron collisionality. This dependence can point to the DRIBM as the modes that drive the anomalous transport in the plasma edge of highly collisional discharges. The effects of plasma shaping on the H-mode pedestal structure was also investigated. The differences in the predicted H-mode pedestal width and height for the DIII-D discharges with different elongation and triangularities were discussed. For the discharges with higher elongation, it was found that the gradients of the plasma profiles in the H-mode pedestal reach semi-steady states. In these simulations, the pedestal slowly continued to evolve to higher pedestal pressures and bootstrap currents until the peeling ballooning stability conditions were satisfied. The discharges with lower elongation do not reach the semi-steady state, and ELM crashes were triggered at earlier times. The plasma elongation was found to have a stronger stabilizing effect than the plasma triangularity. For the discharges with lower elongation and lower triangularity, the ELM frequency was large, and the H-mode pedestal evolves rapidly. It was found that the temperature of neutrals in the scrape-off-layer region can affect the dynamics of the H-mode pedestal buildup. However, the final pedestal profiles were nearly independent of the neutral temperature. The elongation and triangularity affected the pedestal widths of plasma density and electron temperature profiles differently. This study illustrated a new mechanism for controlling the pedestal bootstrap current and the pedestal stability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Fabrication and Characterization of Core–Shell Density-Graded 316L Stainless Steel Porous Structure

In the present work, cylindrical shape 316 L cellular structure with different densities in its core and outer layer (shell) is fabricated by using carbamide as a space holder, via powder metallurgy route and layer-by-layer technique. The arrangement of the created pore is the same as face-centered cubic atomic structure, and different densities are created in two regions by using carbamide particles having two different sizes in the range of 1.7-2 and 2-2.4 mm. The effect of creating a structure with higher porosity (64.5%) in the core and lower porosity (53.8%) in the shell and vice versa and also change in the ratio of the core to the cylinder cross-sectional areas, on the mechanical properties and compaction load bearing of the fabricated foam samples were investigated. In leaching process of the carbamide particles, as an important step of porous structure’s fabrication, it is shown that discontinuous leaching process is more favorable than continuous, by which it would be possible to remove more carbamide particles (around 20%) at the same time. Furthermore, the deformation of the density-graded foam shows the parallel mechanism in the core and shell sections, and the contribution of each part depends on its density and thickness. The energy absorption behavior of the fabricated specimens is evaluated optimally in terms of the energy absorption value associated with the ideal adsorption behavior. The maximum ideal energy absorption efficiency for the samples with more porosity in the shell was approximately equal to 0.91, while for the sample includes the lower amount of porosity in the shell, this value was in the range of 0.85 and 0.89.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Imprinting high-gradient topographical structures onto optical surfaces using magnetorheological finishing: Manufacturing corrective optical elements for high-power laser applications

Corrective optical elements form an important part of high-precision optical systems. We have developed a method to manufacture high-gradient corrective optical elements for high-power laser systems using deterministic magnetorheological finishing (MRF) imprinting technology. Several process factors need to be considered for polishing ultraprecise topographical structures onto optical surfaces using MRF. They include proper selection of MRF removal function and wheel sizes, detailed MRF tool and interferometry alignment, and optimized MRF polishing schedules. Dependable interferometry also is a key factor in high-gradient component manufacture. A wavefront attenuating cell, which enables reliable measurement of gradients beyond what is attainable using conventional interferometry, is discussed. The results of MRF imprinting a 23 μm deep structure containing gradients over 1.6 μm / mm onto a fused-silica window are presented as an example of the technique’s capabilities. As a result, this high-gradient element serves as a thermal correction plate in the high-repetition-rate advanced petawatt laser system currently being built at Lawrence Livermore National Laboratory.

42 ENGINEERING↗
DOE OSTI2016

CFD and X-ray Analysis of Gaseous Direct Injection from an Outward Opening Injector

Using natural gas in an internal combustion engine (ICE) is emerging as a promising way to improve thermal efficiency and reduce exhaust emissions. In the development of such engine platforms, computational fluid dynamics (CFD) plays a fundamental role in the optimization of geometries and operating parameters. One of the most relevant issues in the simulation of direct injection (DI) gaseous processes is the accurate prediction of the gas jet evolution. The simulation of the injection process for a gaseous fuel does not require complex modeling, nevertheless properly describing high-pressure gas jets remains a challenging task. At the exit of the nozzle, the injected gas is under-expanded, the flow becomes supersonic and shocks occur due to compressibility effects. These phenomena lead to challenging computational requirements resulting from high grid resolution and low computational time-steps. In this paper, the simulation of gaseous injection from an outward opening injector was carried out with the commercial CFD software CONVERGE. Argon was injected into a quiescent chamber filled with Nitrogen. Two free-jet cases with injection pressure of 15 and 10 bar, and ambient pressure of 1 and 2 bar respectively, were tested. For the 15 bar injection pressure and 1 bar ambient pressure, a jet impinging on a barrier case was performed as well. The numerical results were validated against experimental mass measurements using X-ray radiography. The results, presented in terms of 2D maps and transverse sections of equivalent path length of argon, showed good agreement between modeling and experiments in terms of jet penetration, shape and gas mass distribution. These emphasized the significant impact of mesh size on the results for the simulation of this type of injector. Due to the small geometrical dimensions a fine mesh is required in order to deliver an accurate description of the mixture formation. Numerical results were discussed with a particular emphasis on the application to engine simulations, where a trade-off between accuracy and computational time is a necessity.

Bartolucci, Lorenzo↗
DOE OSTI2019

Clustering representative days for power systems generation expansion planning: Capturing the effects of variable renewables and energy storage

Decision makers rely on models to make important regulatory, policy, and investment decisions. For power systems, these models must capture (i) the future challenges introduced by the integration of large quantities of variable renewable energy sources and (ii) the role that energy storage technologies should play. In this paper, we explore several different approaches to selecting representative days for generation expansion planning models, focusing on capturing these dynamics. Further, we propose a new methodology for adjusting the outputs of clustering algorithms that provides three advantages: the targeted level of net demand is captured, the underlying net demand shapes that define ramping challenges are accurately represented, and the relationship between annual energy and peak demand is captured. This weighting methodology reduces the magnitude of the error in the representative day based generation expansion planning models estimation of costs by 61% on average. The results also demonstrate the importance of carefully performing the clustering of representative days for both system costs and technology mix. In most cases improvements to the total cost of different representative day based expansion plans are realised where conventional generation capacity is substituted for energy storage. Based on the energy storage technology selected we conclude this capacity is being used to address ramping challenges as opposed to shifting renewable generation from off to on peak periods, reinforcing the importance of capturing detailed intraday dynamics in the representative day selection process.

Clustering↗
DOE OSTI2024

Optimizing DOE Opportunities to Research Land–Atmosphere Interactions in the U.S. Southeast (Workshop Report)

The southeastern United States (Southeast), with its complex and varied environments, is an area of tremendous economic, ecological, and societal importance to the country. The region is characterized by heterogeneous landscapes (i.e., geology and soil type) and a long history of human land use coupled with a warm temperature regime and high precipitation. As a result, soil erosion and deposition are pronounced, vegetation recovery is rapid, and human modification is extensive across the region. To better understand land–atmosphere interactions in this important and complex region, research communities supported by the U.S. Department of Energy’s Biological and Environmental Research (BER) program identified the Southeast as a priority region of interest. In fall 2024, the third Atmospheric Radiation Measurement Mobile Facility (AMF3), one of three mobile monitoring facilities designed to collect atmospheric and climate data from undersampled regions around the world, will begin operations in northwestern Alabama’s Bankhead National Forest (BNF). The AMF3-BNF 5-year deployment, from 2024–2029, will monitor the effects of feedbacks among aerosols, clouds, and precipitation on plant physiology and canopy-scale fluxes. It will also focus on scale aggregation to resolve the role of local forcing on larger-scale processes. To enable broader AMF3 involvement by the science community, the BER Environmental System Science (ESS) program organized the Southeast Land– Atmosphere Research Opportunities (SELARO) workshop in August 2023. The purpose was to identify gaps in scientific understanding of terrestrial processes in the Southeast (defined as states bounded by the Gulf of Mexico to the south, the Atlantic Ocean to east, the Mississippi River to the west, and extending through Tennessee and North Carolina to the north) and explore opportunities to use the AMF3-BNF deployment to coordinate and leverage research efforts across the region. Many parts of the Southeast have experienced repeated anthropogenic forcings. Farming, hunting, burning, and settlement of the region by Indigenous Peoples first shaped the distribution of plant communities, which in turn influenced European colonization patterns. Timber harvesting was common during the expansion of European settlements, and production forestry continues today. Agricultural production was extensive and then waned through the 20th century, creating a period of afforestation following agricultural abandonment. Today, many formerly agricultural landscapes are undergoing rapid urbanization and suburbanization. Overlying these patterns of anthropogenic land use are frequent disturbances from hurricanes, tornadoes, wildfires, drought, flooding, ice storms, and the occasional blizzard. An additional characteristic of the Southeast is its overall landscape complexity. Unlike the western United States, where broad expanses may share similar characteristics, Southeast topography, drainage patterns, vegetation, and development patterns vary widely across relatively small spatial scales (<1 km). This is due to the region’s underlying geology and soil development, species biodiversity patterns, and land ownership and use coupled with strong forces of erosion, weathering, and rapid plant growth in the warm, wet climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Geometric concepts for stellarator permanent magnet arrays

The development of stellarators that use permanent magnet arrays to shape their confining magnetic fields has been a topic of recent interest, but the requirements for how such magnets must be shaped, manufactured, and assembled remain to be determined. To address these open questions, we have performed a study of geometric concepts for magnet arrays with the aid of the newly developed MAGPIE code. A proposed experiment similar to the National Compact Stellarator Experiment (NCSX) is used as a test case. Two classes of magnet geometry are explored: curved bricks that conform to a regular grid in cylindrical coordinates, and hexahedra that conform to the toroidal plasma geometry. In addition, we test constraints on the magnet polarization. While magnet configurations constrained to be polarized normally to a toroidal surface around the plasma are unable to meet the required magnetic field parameters when subject to physical limitations on the strength of present-day magnets, configurations with unconstrained polarizations are shown to satisfy the physics requirements for a targeted plasma.

stellarator, permanent magnet, optimization↗
DOE OSTI2017

Method of predicting the mean lung dose based on a patient's anatomy and dose-volume histograms

The aim of this study was to propose a method to predict the minimum achievable mean lung dose (MLD) and corresponding dosimetric parameters for organs-at-risk (OAR) based on individual patient anatomy. For each patient, the dose for 36 equidistant individual multileaf collimator shaped fields in the treatment planning system (TPS) was calculated. Based on these dose matrices, the MLD for each patient was predicted by the homemade DosePredictor software in which the solution of linear equations was implemented. The software prediction results were validated based on 3D conformal radiotherapy (3D-CRT) and volumetric modulated arc therapy (VMAT) plans previously prepared for 16 patients with stage III non–small-cell lung cancer (NSCLC). For each patient, dosimetric parameters derived from plans and the results calculated by DosePredictor were compared. The MLD, the maximum dose to the spinal cord (D{sub max} {sub cord}) and the mean esophageal dose (MED) were analyzed. There was a strong correlation between the MLD calculated by the DosePredictor and those obtained in treatment plans regardless of the technique used. The correlation coefficient was 0.96 for both 3D-CRT and VMAT techniques. In a similar manner, MED correlations of 0.98 and 0.96 were obtained for 3D-CRT and VMAT plans, respectively. The maximum dose to the spinal cord was not predicted very well. The correlation coefficient was 0.30 and 0.61 for 3D-CRT and VMAT, respectively. The presented method allows us to predict the minimum MLD and corresponding dosimetric parameters to OARs without the necessity of plan preparation. The method can serve as a guide during the treatment planning process, for example, as initial constraints in VMAT optimization. It allows the probability of lung pneumonitis to be predicted.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

Facet Impact of CuMn 2 O 4 Spinel Nanocatalysts on Enhancement of the Oxygen Reduction Reaction in Alkaline Media

Surface-structure engineering represents an attractive strategy to optimize the energy conversion performance of nanocatalysts using their deliberately controlled exposed facets. To further exploit the potential of non-PGM-based spinel catalysts for alkaline oxygen reduction reaction (ORR), a cathodic fuel cell reaction, we hereby report a strategy of ORR improvement by controlling the crystallographic facets of ultra-small CuMn 2 O 4 spinel nanocatalysts through a developed colloidal synthesis approach. The synthesis of CuMn 2 O 4 nanocrystals with morphological control relies on the design and selection of the Cu-/Mn-precursors with striking discrepancies in reaction kinetics. Following carbon loading and an annealing post-treatment of the as-synthesized nanocatalysts, the exclusively {101} facet-exposed CuMn 2 O 4 spinel nano-octahedra exhibit improved electrocatalytic activity toward ORR in 1 M KOH, when compared to their spherical counterparts, exhibiting a mass activity (MA) of 37.6 A/g at 0.85 V. After 10,000 cycles of ORR durability test, the nano-octahedra still retain an MA of 24.5 A/g, which is twice that of the CuMn 2 O 4 spinel nanospheres. Structural characterizations after durability testing indicate that the MA decay is likely associated with a decrease in the Mn 3+ fraction and the emergence of Cu + on CuMn 2 O 4 nano-octahedral surfaces. Finally, as a paradigm, this synthesis approach could be extended to other Mn-based spinel nanocatalysts with precise shape control, enabling us to understand and establish the relationship between the surface lattice/valence state and electrocatalytic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantification of LEU Holdup using gamma ray imaging and inverse transport solver

Holdup is the residual amount of special nuclear material (SNM) remaining in a processing facility after the bulk materials have been cleaned out. In commercial uranium processing facilities, quantification of holdup is a major challenge because of the highly variable shapes and sizes of the deposits. Any method that attempts to generalize and calibrate deposit shapes in order to quantify holdup will be prone to high uncertainties. Uncertainties on the order of ±50% are typical in holdup results. In international safeguards applications, a ±50% uncertainty can result in a large amount of material unaccounted for (MUF) thereby increasing the difficulty of detecting material diversion and facility misuse. An imaging-based methodology has been developed with the objective of significantly reducing this uncertainty by using the true deposit shape, instead of relying on oversimplified geometric assumptions. The project is a collaboration between ORNL, Y-12, and the University of Tennessee, Knoxville, TN. Uranium sources of known masses were measured using the Germanium Gamma-ray Imager (GeGI), a high-resolution imaging spectrometer, creating a pixelated map for each spectral bin. Two different gamma imaging methods are employed in this work: coded aperture imaging and Compton imaging. A validated MonteCarlo model of the detector has been developed using the GEANT4 code for determining the intrinsic response of the detector, its enclosure, and the coded aperture mask. An inverse transport solver based on the Markov Chain Monte-Carlo approach known as Differential Evolution Adaptive Metropolis (DREAM) is employed to use the measurement data from the image pixels (coded aperture or Compton) to solve for the mass of 235 U in the deposit. A reliable method based on the DREAM solver has been developed to flag the infinite thickness condition of a uranium deposit. The project team is working towards improving the image reconstruction for Compton imaging so that a better localization of the source can be achieved. Besides treating the coded aperture and Compton imaging methods independently, the project is also evaluating a combined method that uses the Compton scatter data from a coded aperture measurement. GEANT4 simulations are being performed to evaluate the combined approach. The impact on the DREAM optimization as the source thickness progressively approaches infinite thickness is being evaluated. A number of uranium sources available at ORNL have been measured, and the DREAM results have been tested and validated for the coded aperture imaging. A similar effort will be carried out to validate the Compton based method once the development of algorithms for better localization are complete. The imaging based quantification is very amenable to unattended monitoring of holdup accumulation at key measurement points. A proof of concept measurement has been completed to demonstrate this capability The current work used the high energy resolution imager GeGI. However, the approach and methodologies are applicable to other imagers such as the cadmium zin telluride (CZT) based imager manufactured by H3D, Inc.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Ducted Fuel Injection And Cooled Spray Technologies For Particulate Control In Heavy-duty Diesel Engines (Final Report)

Cooled Spray (CS) and Ducted Fuel Injection (DFI) are in-cylinder technologies for diesel engines that can reduce particulate matter and soot emissions and data has been published showing that these technologies can reduce soot emissions by 75-100% for some engines at some operating conditions. However, little is known about scaling the devices for engine size. Additionally, the performance of either technology over the engine duty cycle has not been explored. This project addresses both of these points through single-cylinder engine investigations. The objectives of this project are to provide details about dimensional scaling of these devices and to demonstrate 75% PM reduction over a range of operating conditions on a single-cylinder engine. Two engines were used for this project: a 125mm bore optically accessible engine at Sandia National Laboratories and a 168mm bore metal engine at Southwest Research Institute. The optical engine was used to study the performance of DFI and CS inserts for a large injector orifice diameter injector that is characteristic of a locomotive engine and to perform scaling studies for DFI. The metal engine was used to perform scaling and alignment studies for CS and to evaluate the technology for both EGR and non-EGR engines over the engine operating map. Modifications were required for both engines to accept the prototype inserts being tested. The optical engine required a new fuel injector, cylinder head and piston so that tests could be run at the pressures and engine speeds required. Additionally, a novel rotating stage was designed for the optical engine to simplify alignment of the modules. The metal engine required a modified cylinder head to accept CS inserts and a modified piston to provide additional space around the fuel injector for the CS inserts. Tests on the optical engine showed that DFI reduces PM emissions for both small injector orifices (0.170mm diameter) and large injector orifices (0.290mm). For high load testing, the DFI modules were not as effective as at low load testing, but it was acknowledged that minimal geometric optimization was performed and more improvements may be possible. Comparing DFI to CS and conventional diesel combustion (CDC), DFI performed better than CS or CDC. The CS geometries used in these studies may not be ideal for that engine and additional modifications likely would improve performance. Tests on the metal engine showed PM reductions as high has 80% at some operating conditions with duty-cycle PM reductions of ~50% for EGR and non-EGR configurations. The CS testing on the metal engine showed that chamfering of the fuel passage inlet either through hydro-erosion or mechanical grinding provided significant improvements in the PM reduction capabilities of the insert. Additionally, alignment sensitivities were explored and the data show that the tolerance to misalignment is approximately 0.05 to 0.1mm for the inserts that were studied here. Air-fuel ratio was shown to be important in the effectiveness of the CS inserts. In several tests, it was shown that the CS inserts are more effective at reducing the PM for high-AFR operating conditions compared to low AFR conditions. In summary, multiple designs were evaluated on both engines. It was found that for the conditions and configurations studied here, a fuel passage diameter of ~2.5mm performed best overall. Significant duty-cycle PM reductions are possible using these technologies and sensitivities to AFR, alignment fuel passage diameter and inlet fuel passage shaping were explored and are reported here. More PM reduction may be possible with improved geometric design and attention to alignment practices.

02 PETROLEUM↗
DOE OSTI2019

Droplet size distributions in turbulent clouds: experimental evaluation of theoretical distributions

Precipitation efficiency and optical properties of clouds, both central to determining Earth's weather and climate, depend on the size distribution of cloud particles. Here in this paper, theoretical expressions for cloud droplet size distribution shape are evaluated using measurements from controlled experiments in a convective-cloud chamber. The experiments are a unique opportunity to constrain theory because they are in steady-state and because the initial and boundary conditions are well characterized compared to typical atmospheric measurements. Three theoretical distributions obtained from a Langevin drift-diffusion approach to cloud formation via stochastic condensation are tested: (a) stochastic condensation with a constant removal time-scale; (b) stochastic condensation with a size-dependent removal time-scale; (c) droplet growth in a fixed supersaturation condition and with size-dependent removal. In addition, a similar Weibull distribution that can be obtained from the drift-diffusion approach, as well as from mechanism-independent probabilistic arguments (e.g., maximum entropy), is tested as a fourth hypothesis. Statistical techniques such as the χ2 test, sum of squared errors of prediction, and residual analysis are employed to judge relative success or failure of the theoretical distributions to describe the experimental data. An extensive set of cloud droplet size distributions are measured under different aerosol injection rates. Five different aerosol injection rates are run both for size-selected aerosol particles, and six aerosol injection rates are run for broad-distribution, polydisperse aerosol particles. In relative comparison, the most favourable comparison to the measurements is the expression for stochastic condensation with size-dependent droplet removal rate. However, even this optimal distribution breaks down for broad aerosol size distributions, primarily due to deviations from the measured large-droplet tail. Lastly, a possible explanation for the deviation is the Ostwald ripening effect coupled with deactivation/activation in polluted cloud conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Plant Diversity and Fertilizer Management Shape the Belowground Microbiome of Native Grass Bioenergy Feedstocks

Plants may actively cultivate microorganisms in their roots and rhizosphere that enhance their nutrition. To develop cropping strategies that substitute mineral fertilizers for beneficial root symbioses, we must first understand how microbial communities associated with plant roots differ among plant taxa and how they respond to fertilization. Arbuscular mycorrhizal (AM) fungi and rhizobacteria are of particular interest because they enhance nutrient availability to plants and perform a suite of nutrient cycling functions. The purpose of this experiment is to examine the root and soil microbiome in a longterm switchgrass ( Panicum virgatum ) biofuel feedstock experiment and determine how AM fungi and rhizobacteria respond to plant diversity and soil fertility. We hypothesize that intra- and interspecific plant diversity, nitrogen fertilization (+N), and their interaction will influence the biomass and community composition of AM fungi and rhizobacteria. We further hypothesize that +N will reduce the abundance of nitrogenase-encoding nifH genes on the rhizoplane. Roots and soils were sampled from three switchgrass cultivars (Cave-in-Rock, Kanlow, Southlow) grown in monoculture, intraspecific mixture, and interspecific planting mixtures with either Andropogon gerardii or diverse native tallgrass prairie species. Molecular sequencing was performed on root and soil samples, fatty acid extractions were assessed to determine microbial biomass, and quantitative polymerase chain reaction (qPCR) was performed on nifH genes from the rhizoplane. Sequence data determined core AM fungal and bacterial microbiomes and indicator taxa for plant diversity and +N treatments. We found that plant diversity and +N influenced AM fungal biomass and community structure. Across all plant diversity treatments, +N reduced the biomass of AM fungi and nifH gene abundance by more than 40%. The AM fungal genus Scutellospora was an indicator for +N, with relative abundance significantly greater under +N and in monoculture treatments. Community composition of rhizobacteria was influenced by plant diversity but not by +N. Verrucomicrobia and Proteobacteria were the dominant bacterial phyla in both roots and soils. Our findings provide evidence that soil fertility and plant diversity structure the root and soil microbiome. Optimization of soil communities for switchgrass production must take into account differences among cultivars and their unique responses to shifts in soil fertility.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Neutronics Analysis Supporting Baseline Fuel-Rodlet Irradiation Design in TREAT Multi-SERTTA Vehicle

The Transient Reactor Test (TREAT) Facility is a thermal-spectrum test reactor designed to evaluate reactor fuels and structural material under simulated nuclear excursion and transient power/cooling mismatch situations in nuclear reactors. TREAT was historically operated from 1959 through 1994. There are ongoing activities to resume its transient testing capabilities to support advanced reactor designs and validate computational predictions of fuel and core behavior. Such validation will be an important factor for review and approval of next generation power reactor design and safety evaluations. Historic methods used in TREAT were very approximate and served to estimate a predicted transient power shape. A number of calibration test experiments were then required prior to the actual transient test to better estimate the power that would have been deposited in the experiment during a transient pulse. Much of the experience in the use of these legacy methods has also been lost over the past two decades. Monte Carlo criticality codes are being utilized for calculation of the initial critical state, control rod reactivity worth, and power coupling between core and experiments. Current efforts primarily support neutronics characterization with safety and design analyses for the Multi-SERTTA (Static Environment Rodlet Transient Test Apparatuses) experiment test vehicle, which can accommodate up to four concurrent rodlets under separate environmental conditions. Ongoing efforts are necessary in preparation of transient testing support for Accident Tolerant Fuel (ATF) characterization and other future fuels testing requirements. The power coupling factor (PCF) is one of two key parameters utilized in designing and characterizing the performance of a TREAT fuel irradiation experiment. The second key parameter is the Transient Correction Factor (TCF) which accounts for time-dependent effects. The PCF is the effective quantity of fission-generated energy per mass of fuel specimen per total core energy. The quantity of energy deposited into a fuel rodlet under various test environments can be used to evaluate overall fuel performance, including failure mechanisms and thresholds. Baseline fuel rodlet calculations have been performed to estimate the effectiveness of the Multi-SERTTA test vehicle and optimization of PCF for various test conditions. The baseline serves as reference fuel pin conditions against which calibration and baseline experiments can be performed prior to actual ATF transient testing. Baseline experiments further serve to validate computational analyses and experimental test methods prior to design planning, experimentation, and evaluation of the ATF candidates. This paper serves to document activities performed to evaluate a reference Pressurized Water Reactor (PWR) rodlet design in the TREAT Multi-SERTTA test vehicle. (authors)

42 ENGINEERING↗
DOE OSTI2016

Parameterization retrieval of trace gas volume mixing ratios from Airborne MAX-DOAS

We present a parameterization retrieval of volume mixing ratios (VMRs) from differential slant column density (dSCD) measurements by Airborne Multi-AXis Differential Optical Absorption Spectroscopy (AMAX-DOAS). The method makes use of the fact that horizontally recorded limb spectra (elevation angle 0°) are strongly sensitive to the atmospheric layer at instrument altitude. These limb spectra are analyzed using reference spectra that largely cancel out column contributions from above and below the instrument, so that the resulting limb dSCDs, i.e., the column integrated concentration with respect to a reference spectrum, are almost exclusively sensitive to the atmospheric layers around instrument altitude. The conversion of limb dSCDs into VMRs is then realized by calculating box air mass factors (Box-AMFs) for a Rayleigh atmosphere and applying a scaling factor constrained by O 4 dSCDs to account for aerosol extinction. An iterative VMR retrieval scheme corrects for trace gas profile shape effects. Benefits of this method are (1) a fast conversion that only requires the computation of Box-AMFs in a Rayleigh atmosphere; (2) neither local aerosol extinction nor the slant column density in the DOAS reference (SCD ref ) needs to be known; and (3) VMRs can be retrieved for every measurement point along a flight track, thus increasing statistics and adding flexibility to capture concentration gradients. Sensitivity studies are performed for bromine monoxide (BrO), iodine monoxide (IO) and nitrogen dioxide (NO 2 ), using (1) simulated dSCD data for different trace gas and aerosol profiles and (2) field measurements from the Tropical Ocean tRoposphere Exchange of Reactive halogen species and Oxygenated VOC (TORERO) field experiment. For simulated data in a Rayleigh atmosphere, the agreement between the VMR from the parameterization method (VMR para ) and the true VMR (VMR true ) is excellent for all trace gases. Offsets, slopes and R 2 values for the linear fit of VMR para over VMR true are, respectively (0.008 ± 0.001) pptv, 0.988 ± 0.001, 0.987 for BrO; (-0.0066 ± 0.0001) pptv, 1.0021 ± 0.0003, 0.9979 for IO; (-0.17 ± 0.03) pptv, 1.0036 ± 0.0001, 0.9997 for NO 2 . The agreement for atmospheres with aerosol shows comparable R 2 values to the Rayleigh case, but slopes deviate a bit more from one: (0.093 ± 0.002) pptv, 0.933 ± 0.002, 0.907 for BrO; (0.0021 ± 0.0004) pptv, 0.887 ± 0.001, 0.973 for IO; (8.5 ± 0.1) pptv, 0.8302 ± 0.0006, 0.9923 for NO 2 . VMR para from field data are further compared with optimal estimation retrievals (VMR OE ). Least orthogonal distance fit of the data give the following equations: BrO para = (0.1 ± 0.2) pptv + (0.95 ± 0.14) × BrO OE ; IO para = (0.01 ± 0.02) pptv + (1.00 ± 0.12) × IO OE ; NO 2para = (3.9 ± 2.5) pptv + (0.87 ± 0.15) × NO 2OE . Overall, we conclude that the parameterization retrieval is accurate with an uncertainty of 20 % for IO, 30 % for BrO and NO 2 , but not better than 0.05 pptv IO, 0.5 pptv BrO and 10 pptv NO 2 . Finally, the retrieval is applicable over a wide range of atmospheric conditions and measurement geometries and not limited to the interpretation of vertical profile measurements in the remote troposphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Thermal and Mechanical Energy Performance Analysis of Closed-loop Systems in Hot-Dry-Rock and Hot-Wet-Rock Reservoirs

To understand the potential and limitations for recovering thermal and mechanical energy from closed-loop geothermal systems a collaborative study is underway that will investigate an array of system configurations, working fluids, geothermal reservoir characteristics, operational periods, and heat transfer enhancements. Closed-loop geothermal systems are distinguished from hydrothermal or enhanced geothermal systems (EGS) in that the working fluid only circulates through drilled boreholes. Principal objectives of this study are to determine upper limits for thermal and mechanical energy recovery and optimal operational and configuration parameters for each scenario. Teams of scientists and engineers are applying a suite of numerical simulation and analytical tools to model the heat recovery from closed-loop geothermal systems, and then optimizing operational and configuration parameters to maximize the thermal and mechanical energy recovery. Results from the suite of numerical simulators and analytical tools, such as outlet and inlet states and temperature profiles in the geothermal reservoir over time are intercompared to increase confidence in the analysis. This paper documents the study findings for closed-loop systems in hot-dry-rock and hot-wet-rock reservoirs, where water is the working fluid. The characteristics of the hot-dry-rock reservoir were based on the U.S. Department of Energy’s Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site, near Milford Utah. Two objective functions are defined to optimize the operational and configuration parameters of the system, one each for the recovery of mechanical and thermal energy over the period of operation. For both objective functions, a surface plant thermal to mechanical energy conversion factor and an energy drilling cost is required. In keeping with the study objectives the surface plant conversion factor is determined from a second-law of thermodynamics analysis of a generic binary plant, and drilling costs are based on those from the Utah FORGE site and current national electrical costs.

Closed-loop geothermal systems, hot-dry-rock, hot-↗