Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HOMOGENIZATION METHODS”

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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Microfluidic Uranium Microspheres Production for TRISO and Advanced Fuel Concepts

The purpose of this project was to expand existing internal gelation sol-gel capabilities at PNNL to explore producing uranium dioxide spheres for potential use as fuel kernels in next generation Tri-isotropic (TRISO) particle fuel. This project expands on previous sol-gel efforts at PNNL by (1) increasing the size regime of sphere production from the micro-fluidic range to the milli-fluidic range, and (2) producing uranium spheres. The approach involved first scaling up the channel size of the fluidic system to the millimeter range, with radiation safety considerations in mind; testing and demonstration on non-radioactive surrogate material, cerium oxide; then transitioning to uranium production and finally optimizing system parameters. Commercially available fluidic chips in the desired size range were could not be found, therefore the project designed and fabricated a T-junction with 1mm channels for droplet production. Because the production process is temperature sensitive, prior efforts have involved performing droplet production in a lab freezer. To reduce radiological waste and footprint, two alternative chilling methods were explored using aluminum thermal beads as a chill bath and a custom aluminum block fit to reagent reservoir sizes. Both were successful in the cerium tests, however the aluminum block design outperformed the thermal bead bath and was further adapted for the radiological test and production run in the Radiochemical Processing Laboratory (RPL). Gelation trials were performed to determine an acceptable range of feed solution parameters for the uranium dioxide gels, characterized by R-values, which is determined by the ratio of uranium nitrate to Hexamethylenetetramine (HMTA)/urea in the feed solution. R-values ranging from 1.6-2 were examined in the gelation trials, with only the 1.8 condition being tested in production. The project was successful in demonstrating a proof of concept design for producing uranium dioxide spheres, however further optimization is needed to dial in production parameters and improve sphere quality and homogeneity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1994

Ice surface temperature retrieval from AVHRR, ATSR, and passive microwave satellite data: Algorithm development and application

One essential parameter used in the estimation of radiative and turbulent heat fluxes from satellite data is surface temperature. Sea and land surface temperature (SST and LST) retrieval algorithms that utilize the thermal infrared portion of the spectrum have been developed, with the degree of success dependent primarily upon the variability of the surface and atmospheric characteristics. However, little effort has been directed to the retrieval of the sea ice surface temperature (IST) in the Arctic and Antarctic pack ice or the ice sheet surface temperature over Antarctica and Greenland. The reason is not one of methodology, but rather our limited knowledge of atmospheric temperature, humidity, and aerosol vertical, spatial and temporal distributions, the microphysical properties of polar clouds, and the spectral characteristics of snow, ice, and water surfaces. Over the open ocean the surface is warm, dark, and relatively homogeneous. This makes SST retrieval, including cloud clearing, a fairly straightforward task. Over the ice, however, the surface within a single satellite pixel is likely to be highly heterogeneous, a mixture of ice of various thicknesses, open water, and snow cover in the case of sea ice. Additionally, the Arctic is cloudy - very cloudy - with typical cloud cover amounts ranging from 60-90 percent. There are few observations of cloud cover amounts over Antarctica. The goal of this research is to increase our knowledge of surface temperature patterns and magnitudes in both polar regions, by examining existing data and improving our ability to use satellite data as a monitoring tool. Four instruments are of interest in this study: the AVHRR, ATSR, SMMR, and SSM/I. Our objectives are as follows. Refine the existing AVHRR retrieval algorithm defined in Key and Haefliger (1992; hereafter KH92) and applied elsewhere. Develop a method for IST retrieval from ATSR data similar to the one used for SST. Further investigate the possibility of estimating surface temperature from passive microwave data (in conjunction with AVHRR clear sky samples) through the use of 'effective emissivities' and physical relationships between skin temperature and subsurface temperature. Use the general method outlined in MK93 to calculate a 12-year record of clear sky equivalent surface temperatures, or possibly all-sky snow-ice interface physical temperatures, from SMMR and SSM/I, compare these temperatures to climatologies, ECMWF modeled surface temperatures, and surface temperatures predicted by a 2-D ice model. And intercompare several ice surface retrieval methods and validate them against ground measurements from the Swiss Camp on the Greenland ice sheet. Additionally, we intend to develop a surface temperature product based on AVHRR data and possibly blended with drifting buoy and meteorological station temperatures.

Key, Jeff↗
DOE OSTI2021

Emerging contact force heterogeneity in ordered soft granular media

Under external perturbations, inter-particle forces in disordered granular media are well known to form a heterogeneous distribution with filamentary patterns. Better understanding these forces and the distribution is important for predicting the collective behavior of granular media, the media second only to water as the most manipulated material in global industry. However, studies in this regard so far have been largely confined to granular media exhibiting only geometric heterogeneity, leaving the dimension of mechanical heterogeneity a rather uncharted area. In this work, through a FEM contact mechanics model, we show that a heterogeneous inter-particle force distribution can also emerge from the dimension of mechanical heterogeneity alone. Specifically, we numerically study inter-particle forces in packing of mechanically heterogeneous disks arranged over either a square or a hexagonal lattice and under quasi-static isotropic compression. Our results show that, at the system scale, a hexagonal packing exhibit a more heterogeneous inter-particle force distribution than a square packing does; At the particle scale, for both packing lattices, preliminary analysis shows the consistent coexistence of outliers (i.e., softer disks sustaining larger forces while stiffer disks sustaining smaller forces) in comparison to their homogeneous counterparts, which implies the existence of nonlocal effect. Further analysis on the portion of outliers and on spatial contact force correlations suggest that the hexagonal packing shows more pronounced nonlocal effect over the square packing under small mechanical heterogeneity. However, such trend is reversed when assemblies becomes more mechanically heterogeneous. Lastly, we confirm that, in the absence of particle reorganization events, contact friction merely plays the role of packing stabilization while its variation has little effect on inter-particle forces and their distribution.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Burning Velocity Measurements in Aluminum-Air Suspensions using Bunsen Type Dust Flames

Laminar burning velocity (sometimes also referred in literature as fundamental or normal flame propagation speed) is probably the most important combustion characteristic of the premixed combustible mixture. The majority of experimental data on burning velocities in gaseous mixtures was obtained with the help of the Bunsen conical flame. The Bunsen cone method was found to be sufficiently accurate for gaseous mixtures with burning velocities higher than 10-15 cm/s at normal pressure. Hans Cassel was the first to demonstrate that suspensions of micron-size solid fuel particles in a gaseous oxidizer can also form self-sustained Bunsen flames. He was able to stabilize Bunsen flames in a number of suspensions of different nonvolatile solid fuels (aluminum, carbon, and boron). Using the Bunsen cone method he estimated burning velocities in the premixed aluminum-air mixtures (particle size less than 10 microns) to be in the range of 30-40 cm/s. Cassel also found, that the burning velocity in dust clouds is a function of the burner diameter. In our recent work, we have used the Bunsen cone method to investigate dependence of burning velocity on dust concentration in fuel-rich aluminum dust clouds. Burning velocities in stoichiometric and fuel-rich aluminum dust suspensions with average particle sizes of about 5 microns were found to be in the range of 20-25 cm/s and largely independent on dust concentration. These results raise the question to what degree burning velocities derived from Bunsen flame specifically and other dust flame configurations in general, are indeed fundamental characteristics of the mixture and to what degree are they apparatus dependent. Dust flames in comparison to gas combustion, are thicker, may be influenced by radiation heat transfer in the flame front, respond differently to heat losses, and are fundamentally influenced by the particular flow configuration due to the particles inertia. Since characteristic spatial scales of dust flames are larger, one can expect that they will also be more sensitive than homogeneous combustion to a particular experimental geometric configuration of the flame and the flow. With such sensitivity the introduction of the very concept of the fundamental flame speed may be problematic for dust combustion. With this in mind, the objective of the present work is to further investigate Bunsen dust flames and evaluate to what degree burning velocities derived from Bunsen cone depend on experimental conditions (i.e. flow rate and nozzle diameter).

Lee, John↗
DOE OSTI2021

Ionic Dopant-Induced Ordering Enhances the Thermoelectric Properties of a Polythiophene-Based Block Copolymer

Conjugated polymer-based block copolymers (CP-BCPs) are an unexplored class of materials for organic thermoelectrics. In this study, the authors report on the electronic conductivity (σ) and Seebeck coefficient (α) of a newly synthesized CP-BCP, poly(3-hexylthiophene)-block-poly (oligo-oxyethylene methacrylate) (P3HT-b-POEM), upon solution co-processing with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI), and subsequently vapor-doping with a molecular dopant, 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ). It is found that the addition of the hydrophilic block POEM greatly enhances the processability of P3HT, enabling homogeneous solution-mixing with LiTFSI. Notably, interactions between P3HT-b-POEM with ionic species significantly improve molecular order and unexpectedly cause electrical oxidizing doping of P3HT block both in solution and solid-states, a phenomenon that has not been previously observed in Li-salt containing P3HT. Vapor doping of P3HT-b-POEM-LiTFSI thin films with F4TCNQ further enhances σ and yields a thermoelectric power factor PF = α 2 σ of 13.0 μW m -1 K -2 , which is more than 20 times higher than salt-free P3HT-b-POEM sample. Through modeling thermoelectric behaviors of P3HT-b-POEM with the Kang-Snyder transport model, the improvement in PF is attributed to higher electronic charge mobility originating from the enhanced molecular ordering of P3HT. The results demonstrate that solution co-processing CP-BCPs with a salt is a powerful method to control structure and performance of organic thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Leak-tight crack repair for 304L stainless steel using friction surfacing

Here, friction surfacing (FS) was investigated as a method for sealing 50 μm wide thru-cracks in 3-mm-thick 304L stainless steel (304L) plate. Friction surfacing was performed using two sizes of 304L consumable rod, with diameters of 9.52 mm and 12.7 mm, on two substrate conditions (clean and oxidized). Friction surfacing was able to deposit a 600 μm thick coating and repair cracks to a depth of 100 μm - 200 μm below the original surface, when using the 12.7-mm-diameter consumable rod for both substrate conditions. Helium leak rates of 10 -10 atm-cc/s were achieved on crack repairs, designating them as leak-tight by the ANSI N14.5 standard. Optical and scanning electron microscopy were used to investigate the coating microstructure and bond interfaces. In this study, the 12.7-mm-diameter rod performed better with disruption of the oxide layer on the oxidized substrate and providing a more homogenous coating. Mechanical properties of the coated samples were evaluated by performing micro-indentation, tensile, bending, and adhesion testing. Higher hardness in the coatings was observed due to the fine equiaxed grains resulting from dynamic recrystallization. Complete consolidation of plastically deformed material on the substrate and a strong diffusion bond across the interface was observed. This was reflected in tensile and bending properties of the FS coatings being comparable to those of uncoated specimen. No delamination of the FS coating was observed during adhesion testing. The results demonstrate that friction surfacing is a viable option to seal cracks in stainless steels.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Control of Meridional Flow in Circular Cylinders by a Travelling Axial Magnetic Field

Convective flow in a Bridgman or float zone configuration significantly affects the interface shape and segregation phenomena. While the primary causative factor for this flow is buoyancy induced convection in an enclosed Bridgman melt, the presence of a free surface gives rise to surface tension driven flows in the floating zone processing of melts. It is of interest to curtail these flows in order to realize near quiescent growth conditions that have shown to result in crystals with good longitudinal and radial homogeneity and thereby of better overall quality. While buoyancy effects can be reduced by careful processing in a low gravity (space) environment, the reduction of Marangoni flows due to surface tension variations is not that straight forward. Attempts have been made with some limited success with the use of external fields to affect the melt thermo-fluid behavior. The use of a static magnetic field that reduces convective contamination through the effects of a non-intrusively induced, dissipative Lorentz force in an electrically conducting melt is one such approach. Experiments have shown that axial fields of the order of 5 Tesla can significantly eliminate convection and yield close to diffusion limited crystal growth conditions. The generation and use of such high magnetic fields require substantial hardware and incur significant costs for its operation. Lately, the use of rotating magnetic fields has been tested in semiconductor crystal growth. The method is fairly well known and commonly used in metal processing but its adaptation to crystal growth of semiconductors is fairly recent. The elegance of the technique rests in its low power requirement (typically 10-20 milli-Tesla at 50-400 Hz) and its efficacy in curtailing deleterious temperature fluctuations in the melt. A rotating magnetic field imposes a rotational force and thereby induces a circulation within the melt that tends to dominate other sporadic convective effects. Thus a known low level of convective flow is introduced into the system. A new novel variation of the Lorentz force mechanism is proposed and investigated in this study. Since one of the desired process conditions in melt crystal growth is the minimization of convective effects, this investigation examines the use of an external field of magnetic origin to counteract existing convective flow within the melt. This is accomplished by utilizing a running or traveling axial magnetic wave in the system. The concept is similar to the use of vibrational means in order to induce streaming flows that oppose buoyant or surface tension driven convection in the system. The rotation direction as well as the magnitude (strength) of this circulation can be easily controlled by external inputs thus affording a direct means of controlling the developing shape of the crystallizing front (interface). The theoretical model of this technique is fully developed and presented in this paper. Results from the solution of the developed governing equations and boundary conditions are also presented. An experimental demonstration of the concept is presented through the suppression of natural convective flow in a mercury column. Implications to crystal growth systems will be fully explored in the final manuscript.

Mazuruk, K.↗
DOE OSTI2020

Development of Mischmetal–Fe–Co–B Permanent Magnet Alloys via High-Throughput Methods

Additive manufacturing synthesis using laser engineered net shaping (LENS) is utilized in this work to rapidly print libraries of mischmetal (MM = La, Ce, Nd, and Pr) containing R 2 TM 14 B alloys (R = MM + separated Nd and TM = Fe and Co) enabling robust evaluation of physical properties over a wide composition range. High-throughput characterization of the magnetic and thermal properties are used to identify compositions for potential high-temperature, high-performance permanent magnets with reduced critical rare-earth elements. Improved Curie temperature ( T c ~ 450 °C) is obtained with substitution of Fe by Co in pseudoternary R 2 TM 14 B alloys. Furthermore, a 4-fold decrease in the Nd content can be achieved through substitution with less critical Ce- and La-rich MM, while retaining high T c . Guided by the properties of the LENS printed samples, selected compositions with and without TiC additions are synthesized via melt-spinning techniques to produce nanostructured ribbons. The maximum room temperature coercivity ( H c ) and energy product ((BH) max ) without TiC are found to be 5.8 kOe, 8.5 MGOe, respectively, while TiC additions as a grain refiner gave H c and (BH) max of 4.9 kOe, 9.8 MGOe, respectively. Structural characterization of the melt-spun ribbons shows homogeneous grain refinement with TiC additions, which leads to an increase in the energy product.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Refining Methods to Determine the Isotopic Composition of Uranium Particles by Laser Ablation MC-ICP-MS

We report on efforts to mitigate the generation of isotopic anomalies during the ablation of micrometer-sized uranium oxide particles and analysis by MC-ICP-MS. The results of testing on particles of U200 indicate that laser fluence and frequency can affect isotopic data produced by laser ablation, but no settings were tested that could eradicate the signal spiking effect and generation of anomalous isotopic data for 234U/238U and 236U/238U. These anomalies are more frequent in samples with higher 235U enrichments, which is expected given their higher abundances of 234U and 236U. Of the standards tested here, only U005-A was anomaly-free. Efforts to compare the laser traces for particles with normal and anomalous isotopic compositions showed subtle differences in the behavior of the 234U/238U and 236U/238U ratios. However, it would be challenging to identify anomalous data points from a population of unknowns using this distinction, i.e., this is unlikely to be diagnostic. Thus, using current analytical hardware, we cannot eradicate the signal spiking phenomenon and cannot unambiguously identify isotopic anomalies from laser ablation traces. Ultimately, laser ablation ICP-MS would either require dramatic improvement to the ablation process and generation of more homogenous aerosols and/or improvements to detector electronics to identify and correct for the spiking phenomenon. Even if the reliability of the technique could be improved, a broader question to address is whether current data quality is of a high enough standard for laser ablation MC-ICP-MS to be used as a complementary technique to LG-SIMS for Safeguards.

organic↗
DOE OSTI2024

Key Experimental Considerations When Evaluating Functional Ionic Liquids for Combined Capture and Electrochemical Conversion of CO 2

Ionic liquids (ILs) are considered functional electrolytes for the electrocatalytic reduction of CO 2 (ECO 2 R) due to their role in the double-layer structure formation and increased CO 2 availability at the electrode surface, which reduces the voltage requirement. However, not all ILs are the same, considering the purity and degree of the functionality of the IL. Further, there are critical experimental factors that impact the evaluation of ILs for ECO 2 R including the reference electrode, working electrode construction, cosolvent selection, cell geometry, and whether the electrochemical cell is a single compartment or a divided cell. Here, we describe improved synthesis methods of imidazolium cyanopyrrolide IL for electrochemical studies in consideration of precursor composition and reaction time. We explored how IL with cosolvents (i.e. acetonitrile, dimethylformamide, dimethyl sulfoxide, propylene carbonate, and n-methyl-2-pyrrolidone) affects conductivity, CO 2 mass transport, and ECO 2 R activation overpotential together with the effects of electrode materials (Sn, Ag, Au, and glassy carbon). Acetonitrile was found to be the best solvent for lowering the onset potential and increasing the catalytic current density for the production of CO owing to the enhanced ion mobility in combination with the silver electrode. Further, the ECO 2 R activity of molecular catalysts Ni(cyclam)Cl 2 and iron tetraphenylsulfonato porphyrin (FeTPPS) on the carbon cloth electrode maintained high Faradaic efficiencies for CO in the presence of the IL. This study presents best practices for examining nontraditional multifunctional electrolytes amenable to integrated CO 2 capture and conversion technologies for homogeneous and heterogeneous ECO 2 R.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Source scaling comparison and validation in Central Italy: data intensive direct S waves versus the sparse data coda envelope methodology

SUMMARY Robustness of source parameter estimates is a fundamental issue in understanding the relationships between small and large events; however, it is difficult to assess how much of the variability of the source parameters can be attributed to the physical source characteristics or to the uncertainties of the methods and data used to estimate the values. In this study, we apply the coda method by Mayeda et al. using the coda calibration tool (CCT), a freely available Java-based code (https://github.com/LLNL/coda-calibration-tool) to obtain a regional calibration for Central Italy for estimating stable source parameters. We demonstrate the power of the coda technique in this region and show that it provides the same robustness in source parameter estimation as a data-driven methodology [generalized inversion technique (GIT)], but with much fewer calibration events and stations. The Central Italy region is ideal for both GIT and coda approaches as it is characterized by high-quality data, including recent well-recorded seismic sequences such as L'Aquila (2009) and Amatrice–Norcia–Visso (2016–2017). This allows us to apply data-driven methods such as GIT and coda-based methods that require few, but high-quality data. The data set for GIT analysis includes ∼5000 earthquakes and more than 600 stations, while for coda analysis we used a small subset of 39 events spanning 3.5 < Mw < 6.33 and 14 well-distributed broad-band stations. For the common calibration events, as well as an additional 247 events (∼1.7 < Mw < ∼5.0) not used in either calibration, we find excellent agreement between GIT-derived and CCT-derived source spectra. This confirms the ability of the coda approach to obtain stable source parameters even with few calibration events and stations. Even reducing the coda calibration data set by 75 per cent, we found no appreciable degradation in performance. This validation of the coda calibration approach over a broad range of event size demonstrates that this procedure, once extended to other regions, represents a powerful tool for future routine applications to homogeneously evaluate robust source parameters on a national scale. Furthermore, the coda calibration procedure can homogenize the Mw estimates for small and large events without the necessity of introducing any conversion scale between narrow-band measures such as local magnitude (ML) and Mw, which has been shown to introduce significant bias.

Morasca, Paola (ORCID:0000000265254867)↗
DOE OSTI2023

Metaoptics Enabled Multifunctional Imaging

Despite the fact that most microorganisms occupy two- and three-dimensional space in heterogeneous arrangements and in proximity to other microorganisms of different species, much of our current knowledge about their metabolic processes is derived from homogenizing, then analyzing, the chemical content of cultures containing a single microbial species. This project addressed this fundamental dichotomy by developing a set of optical imaging principles combining multiple new technologies. The imaging strategies developed in this project combine metaoptics structures with active in situ nanoscale control of the chemical environment and applied them to a microbial system, Myxococcus xanthus, with particular relevance to the DoE mission. The combination of metaoptical architectures and nanoscale control over the molecular environment enables: (a) precise control over the electromagnetic (EM) field at length scales smaller than the wavelength of light; (b) control of the interaction of the EM field with critical molecular systems in DOE relevant microbes; (c) control over the chemical environment – especially the presence and quantity of reactive oxygen species (ROS) that can affect redox homeostasis; and (d) the ability to ask new kinds of questions not accessible to ‘omics’ approaches or standard methods of biological imaging. These capabilities are applicable to detailed studies of metabolic pathways in microbes and to lignocellulosic biomass deconstruction. To accomplish these objectives, we pursued two over-arching technical goals: (1) the development of new metaoptics-enabled approaches to imaging and spectroscopic characterization; and (2) the development of tools to control the chemical environment of a microbial sample with nanometer-scale precision. Goal 1 was addressed through the design, fabrication, and characterization of new metasurfaces capable of super-resolution imaging through extreme confinement of the optical field. Goal 2 was addressed by controlling the redox potential on the nanoscale in microbial communities and characterizing their effect on intrinsic bacterial fluorophores which act as molecular sentinels and through characterization of soluble factors secreted by Myxococcus xanthus by confocal Raman imaging. The optical imaging/sensing approaches developed here make it possible to use these powerful new imaging and sensing modalities in metabolic studies by making it possible to visualize and track the spatial and temporal expression patterns of natural or engineered pathways in microorganisms.

09 BIOMASS FUELS↗
DOE OSTI2020

Microkinetic Modeling: A Tool for Rational Catalyst Design

The design of heterogeneous catalysts relies on understanding the fundamental surface kinetics that controls catalyst performance, and microkinetic modeling is a tool that can help the researcher in streamlining the process of catalyst design. Microkinetic modeling is used to identify critical reaction intermediates and rate-determining elementary reactions, thereby providing vital information for designing an improved catalyst. In this review, we summarize general procedures for developing microkinetic models using reaction kinetics parameters obtained from experimental data, theoretical correlations, and quantum chemical calculations. We examine the methods required to ensure the thermodynamic consistency of the microkinetic model. We describe procedures required for parameter adjustments to account for the heterogeneity of the catalyst and the inherent errors in parameter estimation. We discuss the analysis of microkinetic models to determine the rate determining reactions using the degree of rate control and reversibility of each elementary reaction. We introduce incorporation of Brønsted–Evans–Polanyi relations and scaling relations in microkinetic models and the effects of these relations on catalytic performance and formation of volcano curves. We review the analysis of reaction schemes in terms of the maximum rate of elementary reactions, and we outline a procedure to identify kinetically significant transition states and adsorbed intermediates. Here, we explore the application of generalized rate expressions for the prediction of optimal binding energies of important surface intermediates and to estimate the extent of potential rate improvement. We also explore the application of microkinetic modeling in homogeneous catalysis, electro-catalysis, and transient reaction kinetics. We conclude by highlighting the challenges and opportunities in the application of microkinetic modeling for catalyst design.

42 ENGINEERING↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
NASA NTRS1987

Study of one-dimensional spectral dynamic equations of the Reynolds stresses in homogeneous anisotropic turbulence: Application to split-spectrum modeling

The CTR numerical data base generated by direct simulation of homogeneous anisotropic turbulence was used to calculate all of the terms in the spectral balance equations for the turbulent Reynolds stresses. The aim in not only to test the main closure assumptions used in the split-spectrum models, but also to try to devise improved hypotheses deduced from the statistical information. Numerical simulations of turbulent flows provide a large amount of data, a thought provoking wealth of information. The main advantage of this type of comparison is that a great variety of flows can be considered, and this is necessary to test closure hypotheses. Moreover various initial conditions can be introduced in the calculation, even if they are not experimentally feasible. All the terms in the spectral equations can be calculated. The limited Reynolds numbers of the simulations and the statistical noise caused by a small sample, particularly at the large scales, causes some difficulty in the interpretation of the results, but the method of approach proved to be a powerful tool for testing and improving spectral closures.

Schiestel, R.↗
NASA NTRS1993

Maximum a posteriori classification of multifrequency, multilook, synthetic aperture radar intensity data

We present a maximum a posteriori (MAP) classifier for classifying multifrequency, multilook, single polarization SAR intensity data into regions or ensembles of pixels of homogeneous and similar radar backscatter characteristics. A model for the prior joint distribution of the multifrequency SAR intensity data is combined with a Markov random field for representing the interactions between region labels to obtain an expression for the posterior distribution of the region labels given the multifrequency SAR observations. The maximization of the posterior distribution yields Bayes's optimum region labeling or classification of the SAR data or its MAP estimate. The performance of the MAP classifier is evaluated by using computer-simulated multilook SAR intensity data as a function of the parameters in the classification process. Multilook SAR intensity data are shown to yield higher classification accuracies than one-look SAR complex amplitude data. The MAP classifier is extended to the case in which the radar backscatter from the remotely sensed surface varies within the SAR image because of incidence angle effects. The results obtained illustrate the practicality of the method for combining SAR intensity observations acquired at two different frequencies and for improving classification accuracy of SAR data.

Rignot, E.↗
DOE OSTI2023

Effect of Ligand Chemistry on the Electronic Properties and Reactivity of Cobalt Acetate Autoxidation Catalysts

Autoxidation chemistry catalyzed by cobalt(II) acetate is an industrial pathway for converting hydrocarbons into oxygenated compounds and has promising potential applications in plastic waste deconstruction. However, the chemical properties of the Co acetate-based catalysts and their roles in the reaction mechanism are poorly understood, and as a result, attempts to redesign the autoxidation chemical process for novel applications lack clear direction. In this work, we investigate the structure and electronic properties of a series of multinuclear Co(III) compounds that have been proposed as active participants in catalytic autoxidation. Through a combination of X-ray spectroscopic measurements [Co K-edge X-ray absorption spectroscopy (XAS), extended X-ray absorption fine structure (EXAFS), Kβ X-ray emission spectroscopy (XES), and high-energy resolution fluorescence detection (HERFD) XAS] and theoretical methods, we characterize the interactions of Co with acetate and hydroxyl ligands in these compounds. Here we show that the substitution of acetate ligands with hydroxyls in these compounds causes an increase in the electron density on Co, driven by the loss of a π back-bonding interaction between Co and acetate and an increase in electron donation to Co from the hydroxyl ligands. We further classify the bonding between Co and acetate based on the orbital overlap and show that the experimental absorption and emission spectra are well-described by the resulting molecular orbitals. Finally, we predict the impacts of acetate/hydroxyl ligand exchange on autoxidation catalytic properties, showing that the reaction free energies for key Co oxidation and reduction steps are extremely sensitive to metal-ligand bonding interactions and thus have meaningful impacts on overall reactivity. The results of this study constitute an important set of design criteria for developing novel homogeneous autoxidation catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Partially polaron-transformed quantum master equation for exciton and charge transport dynamics

Polaron-transformed quantum master equation (PQME) offers a unified framework to describe the dynamics of quantum systems in both limits of weak and strong couplings to environmental degrees of freedom. Thus, the PQME serves as an efficient method to describe charge and exciton transfer/transport dynamics for a broad range of parameters in condensed or complex environments. However, in some cases, the polaron transformation (PT) being employed in the formulation invokes an over-relaxation of slow modes and results in premature suppression of important coherence terms. A formal framework to address this issue is developed in the present work by employing a partial PT that has smaller weights for low frequency bath modes. In this work, it is shown that a closed form expression of a second order time-local PQME including all the inhomogeneous terms can be derived for a general form of partial PT, although more complicated than that for the full PT. All the expressions needed for numerical calculation are derived in detail. Applications to a model of a two-level system coupled to a bath of harmonic oscillators, with test calculations focused on those due to homogeneous relaxation terms, demonstrate the feasibility and the utility of the present approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗