Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “volume distribution”

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,063 records · Page 59

THE POTENTIAL FOR REUSE OF PRODUCED WATER IN THE CRITICAL MINERALS SUPPLY CHAIN IN THE POWDER RIVER BASIN OF WYOMING AND MONTANA

The Powder River Basin (PRB) Carbon Ore, Rare Earth, and Critical Minerals (CORE-CM) project is part of the Department of Energy sponsored CORE-CM initiative focused on domestic production of critical minerals and promoting economic development in traditionally fossil fuel producing basins. An important part of this project is regional assessment of waste streams and development of basinal strategies for waste stream reuse within the critical minerals (CM) supply chain. The availability of water for all parts of the CM supply chain is of particular concern in the arid mountain west. In the Wyoming PRB, average annual production of produced water (PW) by the oil and gas industry is more than 16 billion gallons (2015-2022). Through published data and engagement with industry partners, we conducted a preliminary investigation of PW in the PRB to understand volumes, quality, and the potential for use in the CM supply chain. PW was assessed for potential as a direct feedstock and for use in CM processing of a variety of conventional and unconventional feedstocks, including ore from the Bear Lodge Alkaline Complex, coal, and coal byproducts. Preliminary findings suggest that PRB PW contains significantly lower concentrations of rare earth elements (REE) and lithium than would be considered economic at this time. For example, concentrations of REE+Sc measured in PRB PW are <1 ppb (DE-EE0007603), and only seven PRB PW samples reported in the USGS PW database had lithium concentrations ≥50 ppm. However, enrichment of CM in residual brine and sludge from PW evaporation and solidification treatment ponds may yield concentrations of interest over a pond lifetime and should be considered in future studies. Coal bed methane (CBM) PW accounted for 42% of all PW in the Wyoming PRB in 2022. Compared to non-CBM PW, PRB CBM water quality generally meets standards allowing surface discharge and reuse in local industries such as livestock. CBM PW may help to fill water needs within the CM supply chain. Understanding the distribution of CBM PW that is discharged versus reused will help develop basinal reuse strategies.

Jackson, Lily↗

PV Inverter Systems Enabled by Monolithically Integrated SiC based Four Quadrant Power Switch (4-QPS) [BiDFET]

The purpose of this project was to develop a new breed of Power Conversion Systems (PCS) for PV integration that is enabled by the newly developed 4-Quadrant Single Die SiC Power Semiconductor Switches (4-QPS) or also referred to as “Bidirectional FET (BIDFET)”. This work includes semiconductor die development, advanced packaging, converter design, development, and testing of 4-QPS enabled hardware prototypes at 1 kW (for single phase residential application) and 10 kW (for three phase commercial application). The BiDirectional Field-Effect Transistor (BiDFET) can enable circuit topologies requiring four-quadrant switches, that were earlier designed using discrete combinations of MOSFETs, IGBTs, GaN HEMTs, and PiN diodes. The monolithic nature of the BiDFET allows lower device count, smaller switch volume, lower inductance, and simpler packaging, and hence more reliable and commercially viable implementation in power electronics converters. The matrix converter topologies, now feasible using BiDFETs, can eliminate the bulky and unreliable dc link capacitors or inductors required for conventional voltage-source or current-source converters in ac–ac and ac–dc applications. The 1.2 kV BiDFET has the potential to disrupt all the applications utilizing 1.2 kV switches, including electric vehicle (EV) drivetrain, bidirectional EV chargers, industrial motor drives, solid-state transformers, datacenter power supplies, elevator drives, dc microgrids, energy storage grid integration, solid-state breakers, etc.

24 POWER TRANSMISSION AND DISTRIBUTION↗

H i asymmetries in LVHIS, VIVA, and HALOGAS galaxies

ABSTRACT We present an analysis of morphological, kinematic, and spectral asymmetries in observations of atomic neutral hydrogen (H i) gas from the Local Volume H i Survey (LVHIS), the VLA Imaging of Virgo in Atomic Gas (VIVA) survey, and the Hydrogen Accretion in Local Galaxies Survey. With the aim of investigating the impact of the local environment density and stellar mass on the measured H i asymmetries in future large H i surveys, we provide recommendations for the most meaningful measures of asymmetry for use in future analysis. After controlling for stellar mass, we find signs of statistically significant trends of increasing asymmetries with local density. The most significant trend we measure is for the normalized flipped spectrum residual (Aspec), with mean LVHIS and VIVA values of 0.204 ± 0.011 and 0.615 ± 0.068 at average weighted 10th nearest-neighbour galaxy number densities of log (ρ10/Mpc−3) = −1.64 and 0.88, respectively. Looking ahead to the Widefield ASKAP L-band Legacy All-sky Blind survey on the Australian Square Kilometre Array Pathfinder, we estimate that the number of detections will be sufficient to provide coverage over 5 orders of magnitude in both local density and stellar mass increasing the dynamic range and accuracy with which we can probe the effect of these properties on the asymmetry in the distribution of atomic gas in galaxies.

Reynolds, T. N.↗

Moments of nucleon isovector structure functions in 2+1+1-flavor QCD

We present results on the isovector momentum fraction, $\langle x \rangle_{u–d}$ , helicity moment, $\langle x \rangle_{Δ u–Δd}$ , and the transversity moment, $\langle x \rangle_{δu–δd}$, of the nucleon obtained using nine ensembles of gauge configurations generated by the MILC Collaboration using 2 + 1 + 1 -flavors of dynamical highly improved staggered quarks. The correlation functions are calculated using the Wilson-Clover action, and the renormalization of the three operators is carried out nonperturbatively on the lattice in the RI'–MOM scheme. The data have been collected at lattice spacings a ≈ 0.15 , 0.12, 0.09, and 0.06 fm and $M_π$ ≈ 310 , 220, and 135 MeV, which are used to obtain the physical values using a simultaneous chiral-continuum-finite-volume fit. The final results, in the $\overline{\text{MS}}$ scheme at 2 GeV, are $\langle$ x $\rangle$ u – d = 0.173 ( 14 ) ( 07 ) , $\langle x \rangle_{Δu–Δd}$ = 0.213 ( 15 ) ( 22 ) , and $\langle x \rangle_{δu–δd}$ = 0.208 ( 19 ) ( 24 ) , where the first error is the overall analysis uncertainty and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with other recent lattice calculations and phenomenological global fit values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The SAMI Galaxy Survey: stellar population and structural trends across the Fundamental Plane

ABSTRACT We study the Fundamental Plane (FP) for a volume- and luminosity-limited sample of 560 early-type galaxies from the SAMI survey. Using r-band sizes and luminosities from new multi-Gaussian expansion photometric measurements, and treating luminosity as the dependent variable, the FP has coefficients a = 1.294 ± 0.039, b = 0.912 ± 0.025, and zero-point c = 7.067 ± 0.078. We leverage the high signal-to-noise ratio of SAMI integral field spectroscopy, to determine how structural and stellar population observables affect the scatter about the FP. The FP residuals correlate most strongly (8σ significance) with luminosity-weighted simple stellar population (SSP) age. In contrast, the structural observables surface mass density, rotation-to-dispersion ratio, Sérsic index, and projected shape all show little or no significant correlation. We connect the FP residuals to the empirical relation between age (or stellar mass-to-light ratio Υ⋆ ) and surface mass density, the best predictor of SSP age amongst parameters based on FP observables. We show that the FP residuals (anti)correlate with the residuals of the relation between surface density and Υ⋆ . This correlation implies that part of the FP scatter is due to the broad age and Υ⋆ distribution at any given surface mass density. Using virial mass and Υ⋆, we construct a simulated FP and compare it to the observed FP. We find that, while the empirical relations between observed stellar population relations and FP observables are responsible for most (75 per cent) of the FP scatter, on their own they do not explain the observed tilt of the FP away from the virial plane.

D’Eugenio, Francesco↗

Electrolyte and Cutoff Potential Effects on Cycle Life of Li4Ti5O12/LiNi0.9Mn0.1O2 Batteries for Behind-the-Meter Storage Applications

Behind-the-Meter Storage (BTMS) is a stationary battery energy storage system that is connected to the electrical distribution system on the customer's side of the utility's service meter. BTMS systems are used to store electrical energy from the grid as well as inconstant, renewable energy, such as local solar and wind generation. A successful BTMS system will allow the customer to pair their energy generation and storage to optimize electrical consumption from the grid, improving reliability and minimizing cost. For BTMS applications, batteries must be designed and optimized with different set of criteria from other leading segments of the Li-ion battery market, like transportation, due the system being stationary and proximal to the residential or commercial building it's benefitting. BTMS applications prioritize safety, cost (low/no-critical materials), reliability (20-year calendar life), and durability (10,000 cycle life), while having the ability to (minimally) compromise energy density and rate capability. Lithium titanate (Li4Ti5O12-, LTO) is a promising anode candidate for BTMS applications due to its high safety and capacity retention, while maintaining a reasonable 160 mAhg-1 reversable capacity and composition of relatively abundant materials. (1) Specifically, LTO has a high working voltage which helps to prevent Li dendrite formation, improving safety. Furthermore, LTO also has negligible lithiation-based volume change, leading to less mechanical pulverization, or loss of active material, upon cycling. For the cathode, materials with little or no Co are of high interest due to the high cost and low abundance of Co. LiMn2O4 (LMO) has been paired with LTO for BTMS applications in the past due to its safety, low cost (abundancy), and reasonably high operating voltage. (2-4) However, the low capacity of LMO limits energy density and specific energy. While not the highest priority for BTMS applications, increasing energy density will enable deployment in space constrained BTMS applications and decrease total cost. LiNi0.9Mn0.1O2 (LN-MO) is a recently developed material with promise due to its high operating voltage and relatively low price. (5) However, Ni-rich layered oxides, including LNMO, tend to struggle with capacity retention during high-voltage cycling due to mechanical pulverization, irreversible phase transitions, and unstable solid-electrolyte interphase. The study presented here focuses on building an understanding of how electrolyte solvent and varied cutoff potentials will impact the cycle life of LTO/LN-MO cells. Specifically, a comparison is provided between ethylene carbonate (EC), ethyl methyl carbonate (EMC), fluoroethylene carbonate (FEC), and Gen2 electrolyte solvents with 1M Lithium hexafluorophosphate (LiPF6) salt, cycling to two upper termination potentials, 2.6V and 2.7V. Electrochemical testing and diagnostics (e.g., differential capacity analysis, area specific impedance, constant voltage hold, and rate capability) and post-mortem characterization will be used to understand the aging behavior and failure mechanisms of the 8 cell combinations (four electrolytes and two voltage cutoffs). Cells with FEC electrolyte showed a lower initial capacity compared to cells with Gen2, EMC, and EC cycling at both voltages; however, the cells with FEC showed consistent trends in capacity retention with 2.6V and 2.7V termination potentials, while the cells with the other electrolytes showed much higher rates of capacity loss when cycling to the higher voltage. These results indicate that FEC may play a role in improving durability of high-voltage, Ni-rich electrode systems for use in high-cycle applications, such as BTMS.

electrolyte↗

Tree damage and mortality measurements across seven ForestGEO plots in the tropics between Oct 2016 and Mar 2023

Annual records (raw data) on tree survival and structural completeness of 36,524 trees (2,467 species) collected across 29 censuses in seven tropical forests distributed across the Neotropics (Amacayacu, Colombia; Barro Colorado Island (BCI), Panamá; Yasuní, Ecuador) and Asia (Fushan, Taiwan; Huai Kha Khaeng (HKK), Thailand; Khao Chong (KC), Thailand; Pasoh, Malaysia). This dataset was used to compare aboveground biomass loss via damage to living trees relative to total AGB loss (mortality + damage). Variable definitions: site: name of the ForestGEO plot stemID.ams: Unique ID for the stem in the annual mortality surveys (ams) treeID.ams: Unique ID for the tree in the ams date.full.census: date of the previous full census of the plot, format: YYYY-MM-DD dbh.full.census: diameter at the breast height (dbh in mm) measured during the previous full census of the plot home: height of measurement of dbh (in m) meanWD: species-level wood density (g cm-3) date.ams: date of the ams, format: YYYY-MM-DD status: survival status of the tree (A: alive; D: dead; NF: not found; "?": unknown) H_considering_damage: living length of the main axis in meters; provides an estimate of the amount of remaining living tissues along the main axis of the stem (e.g., the height of breakage or the height discounting wood decay) b: the remaining proportion of branch volume within the living length (b∈[0,1]). weights.ind: frequency of the [size class x species] bins within the forest plot relative to their frequency in the sample. Necessary to extrapolate estimates to the whole plot.

54 ENVIRONMENTAL SCIENCES↗

A multi-criteria land suitability assessment of field allocation decisions for switchgrass

With the goal of evaluating biomass resource availability in certain regions in the United States, the Department of Energy, Department of Agriculture, Environmental Protection Agency, national laboratories, and Forest Service research laboratories, together with academic and industry collaborators, performed a resource analysis study, 2016 Billion-Ton Report: Advancing Domestic Resources for a Thriving Bioeconomy, Volume 1: Economic Availability of Feedstocks, using POLYSYS model to estimate biomass potential. The study provided county- and state-level biomass potential projections, with biomass feedstock production estimates and suitability assessment of lands acreage likely to transition to energy crops. For this analysis, we estimate biomass supply by year 2028, more specifically, lands acreage transition to switchgrass in the Midwest of the United States, across 50 counties in 3 states (Nebraska, Kansas, and Colorado). We propose a multi-criteria site suitability framework, to determine the potential acreage and production for switchgrass, at the field level. Site selection criteria span agronomic and ecosystem services such as previous land cover, crop rotation, soil water holding capacity, crop productivity, field efficiency, topography, soil organic matter, soil leaching propensity, distance to surface water bodies, and soil erosivity. Efficient agricultural field selection for dedicated energy crops field allocation is critical to evaluate biomass resource availability, and it is our position that estimation at the field level would provide a more accurate and useful estimation to support an emerging bioeconomy. Using field suitability scores and county-level assessments, we model land allocation decision for energy crop in the future, estimating the spatial distribution of dynamics of energy crop adoption with priority to fields scoring the highest. Additionally, compared with POLYSYS results from the report mentioned above, our method shows good performance overall, with about 2% and 1.2% of difference in switchgrass acreage allocation and crop production, respectively.

54 ENVIRONMENTAL SCIENCES↗

Field intercomparison of ice nucleation measurements: the Fifth International Workshop on Ice Nucleation Phase 3 (FIN-03)

Abstract. The third phase of the Fifth International Ice Nucleation Workshop (FIN-03) was conducted at the Storm Peak Laboratory in Steamboat Springs, Colorado, in September 2015 to facilitate the intercomparison of instruments measuring ice-nucleating particles (INPs) in the field. Instruments included two online and four offline measurement systems for INPs, which are a subset of those utilized in the laboratory study that comprised the second phase of FIN (FIN-02). The composition of the total aerosols was characterized using the Particle Analysis by Laser Mass Spectrometry (PALMS) and Wideband Integrated Bioaerosol Sensor (WIBS) instruments, and aerosol size distributions were measured by a laser aerosol spectrometer (LAS). The dominant total particle compositions present during FIN-03 were composed of sulfates, organic compounds, and nitrates, as well as particles derived from biomass burning. Mineral-dust-containing particles were ubiquitous throughout and represented 67 % of supermicron particles. Total WIBS fluorescing particle concentrations for particles with diameters of > 0.5 µm were 0.04 ± 0.02 cm−3 (0.1 cm−3 highest; 0.02 cm−3 lowest), typical of the warm season in this region and representing ≈ 9 % of all particles in this size range as a campaign average. The primary focus of FIN-03 was the measurement of INP concentrations via immersion freezing at temperatures > −33 °C. Additionally, some measurements were made in the deposition nucleation regime at these same temperatures, representing one of the first efforts to include both mechanisms within a field campaign. INP concentrations via immersion freezing agreed within factors ranging from nearly 1 to 5 times on average between matched (time and temperature) measurements, and disagreements only rarely exceeded 1 order of magnitude for sampling times coordinated to within 3 h. Comparisons were restricted to temperatures lower than −15 °C due to the limits of detection related to sample volumes and very low INP concentrations. Outliers of up to 2 orders of magnitude occurred between −25 and −18 °C; a better agreement was seen at higher and lower temperatures. Although the 5–10 factor agreement of INP measurements found in FIN-03 aligned with the results of the FIN-02 laboratory comparison phase, giving confidence in progress of this measurement field, this level of agreement still equates to temperature uncertainties of 3.5 to 5 °C that may not be sufficient for numerical cloud modeling applications that utilize INP information. INP activity in the immersion-freezing mode was generally found to be an order of magnitude or more, making it more efficient than in the deposition regime at 95 %–99 % water relative humidity, although this limited data set should be augmented in future efforts. To contextualize the study results, an assessment was made of the composition of INPs during the late-summer to early-fall period of this study inferred through comparison to existing ice nucleation parameterizations and through measurement of the influence of thermal and organic carbon digestion treatments on immersion-freezing ice nucleation activity. Consistent with other studies in continental regions, biological INPs dominated at temperatures of > −20 °C and sometimes colder, while arable dust-like or other organic-influenced INPs were inferred to dominate below −20 °C.

54 ENVIRONMENTAL SCIENCES↗

Investigation of theoretical scaling laws using large eddy simulations for airborne spreading of viral contagion from sneezing and coughing

Using a set of large eddy point-particle simulations, we explore the fluid dynamics of an ejected puff resulting from a cough/sneeze. The ejection contains over 61,000 potentially virus-laden droplets at an injection Reynolds number of about 46,000, comparable to an actual cough/ sneeze. We observe that global puff properties, such as centroid, puff volume, momentum, and buoyancy vary little across realizations. Other properties, such as maximum extent, shape, and edge velocity of the puff, may exhibit substantial variation. In many realizations, a portion of the puff splits off and advances along a random direction, while keeping airborne droplet nuclei afloat. This peeled-off portion provides a mechanism for virus-laden droplets to travel over large distances in a short amount of time. We also observe that the vast majority of droplets remain suspended within the puff after all liquid has evaporated. The main objectives of the study are to (i) evaluate assumptions of Balachandar’s et al. theory [Int. J. Multiphase Flow 132, 103439 (2020)], which include buoyancy effects, shape of the puff, and droplet evaporation rate, (ii) obtain values of closure parameters, which include location and time of the virtual origin, and puff entrainment and drag coefficients, and (iii) evaluate the accuracy of the theory in predicting the shape, size, and location of the puff, as well as droplet number density long after ejection. The theory adequately predicts global puff properties including size, velocity, and distance traveled, the largest size of droplets that exit the puff due to settling, and the droplet size distribution within the puff long after ejection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN

Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.

58 GEOSCIENCES↗

Multiport Control with Partial Power Processing in Solid-State Transformer for PV, Storage, and Fast-Charging Electric Vehicle Integration

This article proposes a multiport control method to enable partial power processing (PPP) in a medium-voltage (MV) multiport solid-state transformer (SST). MV multiport SSTs are promising in integrating low-voltage DC sources or loads such as solar photovoltaic, energy storage, and electric vehicles into smart grids without bulky line-frequency transformers. Compared to voltage-source SST, current-source (CS) SST features single-stage isolated bidirectional AC/AC, AC/DC, or DC/DC conversion using an inductive DC link. For a multiport CS SST, it is revealed in this article that the PPP capability can be enabled through the proposed control without extra hardware, different from the case of voltage-source converters where special hardware architecture is required for the PPP. With the PPP, most power exchange between LV ports is processed by only a fraction of the entire conversion stage, leading to reduced DC-link current, volume, loss, and improved efficiency. The proposed multiport PPP control scheme is analyzed to verify the advantages across a wide voltage and power range against conventional full power processing (FPP) multiport control, using the soft-switching solid-state transformer (S4T) with reduced conduction loss as an example. Comparative experimental results based on a SiC three-port S4T prototype verify the effectiveness of the proposed PPP scheme against the FPP scheme under both steady state and dynamic conditions. Here, the DC-link current reduction is measured to be more than 36%. Significantly, the proposed multiport PPP control scheme is generic and applicable to any hard-switching or soft-switching CS SSTs without extra hardware.

14 SOLAR ENERGY↗

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗

In situ high-temperature 3D imaging of the damage evolution in a SiC nuclear fuel cladding material

Silicon carbide (SiC)-based nuclear fission fuel rod cladding has been considered as one of the possible designs for accident tolerant fuels. It is in the form of a SiC fibre reinforced SiC matrix composite tube (SiC f -SiC m ) with monolithic SiC outer and/or inner coating layers. This study focuses on the deformation and fracture processes in this material using in situ X-ray micro-computed tomography (XCT) at room temperature (RT) and 1200 °C in an inert gas environment in a C-ring compression loading configuration. Prior to testing, local properties and residual stresses were characterised using nanoindentation and Raman spectroscopy since they can impact the mechanical behaviour of the material. The 3D strain distribution, crack formation and propagation processes including the toughening mechanisms (e.g., crack deflection, micro-cracking, crack bridging and bifurcation) are investigated in the coating and underlying composites at RT and 1200 °C. There is no particular sequence which toughening mechanism occurs first – this is very different from the conventional toughening theory in ceramic-matrix composites under uniaxial tension loading. Indeed, no evidence of fibre pull-out or fibre fracture was observed in this SiC f -SiC m nuclear cladding material in the current C-ring compression configuration. The correlation between local measurements and bulk mechanical behaviour are discussed.

3D imaging↗

Dissolution zone model of the oxide structure in additively manufactured dispersion-strengthened alloys

The structural evolution of oxides in dispersion-strengthened superalloys during laser-powder bed fusion is considered in detail. Alloy chemistry and process parameter effects on oxide structure are assessed through a parameter study on the model alloy Ni-20Cr, doped with varying concentrations of Y 2 O 3 and Al. Small angle neutron scattering measurements of the dispersoid size distribution show the dispersoid size increases with higher laser power, slower scan speed, and increasing Y 2 O 3 and Al content. Complementary electron microscopy measurements reveal reactions between Y 2 O 3 and Al, even in nanoscale dispersoids, and the presence of micron-scale oxide slag inclusions in select specimens. A scaling analysis of mass and momentum transport within the melt pool, presented here, establishes that diffusional structural evolution mechanisms dominate for nanoscale dispersoids, while fluid forces and advection become significant for larger slag inclusions. These findings are developed into a theory of dispersoid structural evolution, integrating quantitative models of diffusional processes – dispersoid dissolution, nucleation, growth, coarsening – with a reduced order model of time-temperature trajectories of fluid parcels within the melt pool. Calculations of the dispersoid size in single-pass melting reveal a zone in the center of the melt track in which the oxide feedstock fully dissolves. Within this zone the final Y 2 O 3 size is independent of feedstock size and determined by nucleation and growth kinetics. If the dissolution zones of adjacent melt tracks overlap sufficiently with each other to dissolve large oxides, formed during printing or present in the powder feedstock, then the dispersoid structure throughout the build volume is homogeneous and matches that from a single pass within the dissolution zone. Gaps between adjacent dissolution zones result in oxide accumulation into larger slag inclusions. Predictions of final dispersoid size and slag formation using this dissolution zone model match the present experimental data and explain process-structure linkages speculated in the open literature.

36 MATERIALS SCIENCE↗

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↗

Accumulator Sizing for the Hybrid Hydraulic Electric Architecture (HHEA) Using Dynamic Programming

Abstract The Hybrid Hydraulic Electric Architecture (HHEA) is an approach to electrifying and improving that does not require very large electric machines. The majority power is provided hydraulically via a set of common pressure rails while small electric machines are used to modulate that power to meet exact demand. Each common pressure rail has an accumulator, which keeps the pressure of each rail near constant as long as the volume in the accumulator does not change much. A main pump adds or removes fluid from each accumulator as needed to reduce volume changes in the accumulators. It is the operation of this main pump which is investigated here. In previous analysis it was assumed that the size of the accumulators was large enough to account for any difference in fluid volume without significantly changing the pressure. In this paper, the required accumulator sizes are studied using a dynamic programming approach. At each time step the inlet and outlet of the main pump could be connected to any of the common pressure rails. These decisions for the main pump were optimized to minimize the required accumulator size. Case studies are conducted using drive cycles from a 22-ton excavator, a 20-ton wheel loader and a 5-ton excavator. The machines could be operated with total accumulator sizes of less than 20, 10, and 5 liters respectively. Required accumulator size was found to vary with pump flow rate, the frequency at which the pump/motor was able to change operating conditions, the tolerance level of pressure deviation and the accumulator’s polytropic index.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Fission gas retention of densely packed uranium carbonitride tristructural-isotropic fuel particles in a 3D printed SiC matrix

The Transformational Challenge Reactor (TCR) fuel form was designed to contain large, densely packed uranium carbonitride (UCN) tristructural-isotropic (TRISO) fuel particles within a 3D printed SiC matrix, increasing the uranium density compared to conventional TRISO fuel forms and offering full geometric freedom for core design. Here, this work summarizes initial low-burnup, high-power irradiation testing of TCR fuel materials, including loose UCN TRISO particles and integral fuel compacts with ~55% TRISO particles by volume, to evaluate fission gas retention. Fission gasses were fully retained in all loose particle tests and in integral compacts irradiated at low (<250 °C) surface temperatures. Initial testing at higher (~700–750 °C) fuel surface temperatures showed fission gas release (FGR) and complete fracture of three compacts, but no FGR was observed in later high temperature tests (~300–750 °C) of both fueled compacts and loose TRISO particles. Calculated thermal stresses in the failed compacts were far less than the measured strength of the SiC matrix and the stresses in some failed compacts were less than those in compacts that did not show FGR. Thermal stress-induced matrix cracks also would not cause complete fracture because the tensile stresses transition to compression in the higher temperature regions. Therefore, fuel failure was likely not caused by thermal stresses and may have been related to leakage currents from the electrical heaters and erratic fuel surface temperatures that were only observed in the test for which failure was observed. In any case, the matrix cracks propagated through the coatings of TRISO particles located in the high-density matrix regions on the peripheries of the compacts, resulting in measurable fission gas release. The discussion focuses on the importance of understanding matrix density distributions and the particle-matrix interface properties to prevent matrix cracks from causing TRISO particle failures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗