Search NASA⌕ Search

SEARCH · Search NASA

Results for “CAPTURE EFFECT”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The role of cellular structure, non-equilibrium eutectic phases and precipitates on quasi-static strengthening mechanisms of as-built AlSi10Mg parts 3D printed via laser powder bed fusion

The quasi-static loading strengthening mechanism in the as-built state of an AlSi10Mg alloy 3D-printed via Laser Powder Bed Fusion (LPBF) was thoroughly identified and quantified using state-of-the-art electron microscopy and synchrotron X-ray diffraction techniques. The yield strength was comprehensively modelled through an in-depth characterization and quantification of the microstructural features as well as their effective volumes for strengthening. In particular, the non-equilibrium eutectic network was characterized with a two-fold structure: cell boundary particles as well as intracellular lamellar/fibrous networks, each consisting of discrete phases exhibiting a thru-thickness compositional gradient, a semi-coherent interface with the matrix and an abundance of crystal defects such as nano-sized sub-grains, microstrains and stacking faults. Further, these altogether made the eutectic phase the most potent contributor to the yield strength accounting for ~30–40% of the estimated value (i.e., cell boundary and eutectic network strengthening combined). Precipitates strengthening was identified as the second most potent mechanism via a shearing process only. Moreover, the presented methodology was able to capture the effect of LPBF processing variables on the individual strengthening contributions, e.g., the effect of a lower laser scanning speed on increasing the cell size as well as the mean precipitate size, which are shown to exhibit opposing impacts on strengthening.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Ground-motions site and event specificity: Insights from assessing a suite of simulated ground motions in the San Francisco Bay Area

This article presents the results of a research that is part of a larger collaborative effort between the Lawrence Berkeley National Laboratory and the Pacific Earthquake Engineering Research Center, funded by the US Department of Energy Office of Cybersecurity, Energy Security and Emergency Response. The main objective of this study is to assess a suite of near and far-field simulated ground motions obtained from 20 realizations of an M7 Hayward Fault earthquake in the San Francisco Bay Area, California USA, and inform the selection of rupture simulation parameters leading to strong motions. To this aim, comparisons are conducted with NGA-W2 and directivity ground-motion models and a selected population of records. An archetypal steel moment-resisting frame is utilized to assess infrastructure response distributions. The analyses carried out for each simulated event and subdomain with consistent properties in terms of shallow shear-wave velocity proved to be instrumental for better interpreting the differences between simulated motions and empirical models. The main reasons identified for variances between simulations and empirical relationships included (1) directivity effects fully captured by the simulations across the full breadth of rupture models; (2) site vicinity to ruptures that incorporate large-slip patches, particularly if these are in the forward-directivity direction; and (3) presence of geologic structures that can “trap” seismic waves and produce ground motions with large amplitude and long signal duration. The analyses carried out in this work provide a path for interpreting ground-motion site and event specificity obtained from a suite of physics-based simulations, differing only in the rupture model characterization, to inform the selection of simulation scenarios for site-specific engineering analyses under strong excitations. Evidence from this work points to the possibility that current hazard models may underestimate ground-motion intensities in areas where the combined effect of directivity and site conditions results in large ground-motion amplitudes.

58 GEOSCIENCES↗
DOE OSTI2021

Baryonic effects in the Effective Field Theory of Large-Scale Structure and an analytic recipe for lensing in CMB-S4

Upcoming Large-Scale Structure surveys will likely become the next leading sources of cosmological information, making it crucial to have a precise understanding of the influence of baryons on cosmological observables. The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a consistent way to predict the clustering of dark matter and baryons on large scales, where their leading corrections in perturbation theory are given by a simple and calculable functional form even after the onset of baryonic processes. In this paper, we extend the two-fluid-like system up to two-loop order in perturbation theory. Along the way, we show that a new linear counterterm proportional to the relative velocity of the fluids could generically be present, but we show that its effects are expected to be small in our universe. Regardless, we show how to consistently perform perturbation theory in the presence of this new term. We find that the EFTofLSS at two-loop order can accurately account for the details of baryonic processes on large scales. Additionally, we compare our results to a hydrodynamical N-body simulation at many redshifts and find that the counterterms associated with the leading corrections to dark matter and baryons start to differ between redshifts z ≈ 3 and z ≈ 2, signaling the onset of star-formation physics. We then use these fits to compute the lensing power spectrum, show that the understanding of baryonic processes will be important for analyzing CMB-S4 data, and show that the two-loop EFTofLSS accurately captures these effects for ℓ ≲ 2000. Our results are also potentially of interest for current and future weak lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Non-linear relationships between daily temperature extremes and US agricultural yields uncovered by global gridded meteorological datasets

Global agricultural commodity markets are highly integrated among major producers. Prices are driven by aggregate supply rather than what happens in individual countries in isolation. Furthermore, estimating the effects of weather-induced shocks on production, trade patterns and prices hence requires a globally representative weather data set. Recently, two data sets that provide daily or hourly records, GMFD and ERA5-Land, became available. Starting with the US, a data rich region, we formally test whether these global data sets are as good as more fine-scaled country-specific data in explaining yields and whether they estimate similar response functions. While GMFD and ERA5-Land have lower predictive skill for US corn and soybeans yields than the fine-scaled PRISM data, they still correctly uncover the underlying non-linear temperature relationship. All specifications using daily temperature extremes under any of the weather data sets outperform models that use a quadratic in average temperature. Correctly capturing the effect of daily extremes has a larger effect than the choice of weather data. In a second step, focusing on Sub Saharan Africa, a data sparse region, we confirm that GMFD and ERA5-Land have superior predictive power to CRU, a global weather data set previously employed for modeling climate effects in the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Comparison of Fission Product Release Predictions using PARFUME and BISON with Results from the AGR-3/4 Irradiation Experiment

The PARFUME (PARticle Fuel ModEl) fuel performance modeling code and the BISON nuclear fuel performance application built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library were used to predict the fission product release from tristructural isotropic (TRISO) coated fuel particles and compacts during the third and fourth irradiation experiment of the Advanced Gas Reactor (AGR-3/4) Fuel Development and Qualification Program. The fuel performance modeling codes PARFUME and BISON modeled the AGR-3/4 irradiation experiment using the fuel compact time-averaged volume averaged (TAVA) daily temperatures for a total irradiation duration of 369.1 effective full power days (EFPD) to predict the release fraction of the fission product silver (Ag-110m) from a representative TRISO-coated fuel particle from AGR-3/4 compacts. Post-irradiation examination (PIE) measurements provided data on the release of these fission products in the compacts outside of the silicon carbide (SIC) layer. The PARFUME and BISON results were then compared to the silver release measured from compact gamma scanning. The results showed good agreement between PARFUME and BISON but both codes under-predicted the silver release fraction for all the compacts. In addition, BISON was used to model and predict the fission product concentration radial profile outside of the compacts in capsules’ inner and outer rings. These rings were either comprised of matrix and/or structural graphite. To obtain the concentration profiles of silver, cesium, and strontium, a sorption isotherm model was developed in BISON to capture the effects of fission product transport across the gaps between the concentric rings. The general shape of the concentration radial profiles as calculated by BISON were similar in the inner ring (IR) but varied in the outer ring (OR) depending on the fission product of interest or capsule temperature. Using this methodology and model, BISON now has the capability to aid in developing new fission product diffusion coefficients for matrix or structural graphite materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Carbon dioxide capture with aqueous amino acids: Mechanistic study of amino acid regeneration by guanidine crystallization and process intensification

CO 2 capture from powerplant-generated flue gas via a phase-changing process involving absorption with aqueous amino acids (e.g., glycine or sarcosine) and bicarbonate crystallization with bis-iminoguanidines (e.g., glyoxal-bis-iminoguanidine or GBIG) is investigated in this paper. This process is of high interest due to its potential to decrease the energy penalty for CO 2 capture by significantly reducing the solvent regeneration energy typically associated with aqueous amine solvents. A critical step in the proposed CO 2 capture mechanism is the regeneration of the amino acid by removal of protons and bicarbonate ions from solution through crystallization of GBIGH 2 2+ bicarbonate salt. Here, we investigated the thermodynamics and kinetics of glycine regeneration by crystallization of GBIGH 2 2+ (HCO 3 – ) 2 (H 2 O) 2 . A theoretical model was developed and compared to experimental data to simulate and predict the glycine regeneration and determine its reaction mechanism. This combined experimental and theoretical study led to the conclusion that, while the GBIGH 2 2+ bicarbonate crystallization step provides most of the thermodynamic driving force for the glycine regeneration, the rate-limiting step is the protonation of GBIG prior to crystallization. The CO 2 loading and amino acid regeneration steps were combined into a single, intensified process using a bubble column reactor. The CO 2 loading capacity of GBIG was experimentally determined to be roughly 1.36 mol CO 2 per mol GBIG. These results provide the fundamental basis for developing an effective carbon capture technology with phase-changing amino acid/guanidine absorbents.

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

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Asymptotic-preserving gyrokinetic implicit particle-orbit integrator for arbitrary electromagnetic fields

We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in Ricketson and Chacón (2020) to include finite Larmor-radius (FLR) effects in the presence of electric-field length-scales comparable to the particle gyro-radius (the gyro-kinetic limit). We introduce two modifications to the earlier scheme. The first is the explicit gyro-averaging of the electric field at the half time-step, along with an analogous modification to the current deposition, which we show preserves total energy conservation in implicit PIC schemes. The number of gyrophase samples is chosen adaptively, ensuring proper averaging for large timesteps and the recovery of full-orbit dynamics in the small time-step limit. The second modification is an alternating large and small time-step strategy that ensures the particle trajectory samples gyrophases evenly. We show that this strategy relaxes the time-step restrictions on the scheme, allowing even larger speed-ups than previously achievable. We demonstrate the new method with several single-particle motion tests in a variety of electromagnetic field configurations featuring gyro-scale variation in the electric field. Finally, the results demonstrate the advertised ability to capture FLR effects accurately even when significantly stepping over the gyration time-scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Zeolite Membrane Reactor for Pre-Combustion Carbon Dioxide Capture (Final Scientific/Technical Report)

Water-gas shift (WGS) reaction followed by carbon dioxide (CO 2 ) separation is a critical step in the integrated gasification combined cycle (IGCC) process for fossil-fuel-fired electrical power generation with CO 2 capture. To intensify the IGCC process hydrogen-permselective zeolite membrane reactor offers promise to replace the conventional energy-intensive fixed-bed reactors and solvent-based CO 2 capture units. The objectives of this project were to develop a bench-scale zeolite membrane reactor (total membrane area: 932 cm 2 for a 21-tube membrane bundle) for the water-gas-shift reaction of raw syngas from an oxygen-blown coal-gasifier for H 2 production with simultaneous CO 2 separation at the capability of about 2 kilograms H 2 per day (equivalent to 2 kW IGCC power plant) and to demonstrate significant progress toward achieving overall performance goal of 90% CO 2 capture rate with 95% CO 2 purity at the cost of electricity 30% less than the baseline carbon capture approaches. This report summarizes results obtained in this project on scaling up the zeolite membrane reactor by a factor of 200 in membrane area, tests of the bench-scale zeolite membrane reactor for the water-gas-shift reaction at high temperature and high-pressure, and techno-economic analysis of the integration of the zeolite membrane reactor in IGCC power plant for the electrical generation with CO 2 capture. With effective pore modification by catalytic cracking deposition of MDES (methyldiethoxysilane), fabrication of MFI-type zeolite membranes was successfully scaled-up from a lab-scale disk type to bench-scale multiple-tube bundles on low-cost alumina supports. A Co-Mo based sour shift catalyst was evaluated and used in the zeolite membrane reactors for water-gas-shift reaction. The reaction kinetic and gas-permeation equations were developed and employed in mathematic models to guide/predict experiments/performance of zeolite membrane reactors for water-gas shift reaction. Multiple-tube zeolite membrane bundles and reactors were designed, fabricated and tested for gas separation and water-gas-shift reactions with real raw syngas from a coal-fired gasifier operated at high temperature and pressure. The zeolite membrane reactors demonstrated good long-term thermal and chemical stability and constant H 2 permeance (>300 GPU) together with a Knudsen selectivity (~4.7) for H 2 over CO 2 in the field test (cumulative time >28 hours) with high-sulfur coal-derived syngas. The zeolite-membrane-reactor integrated IGCC process was designed using the performance experimentally measured by the University of Cincinnati team with a single-tube zeolite membrane reactor that offers a CO conversion >98% with more than 90% of CO 2 and H 2 captured. With the above integrative approaches, a techno-economic analysis for a cost-benefit comparison was performed to uncover features that determine the power output, capital expenditure, operating expenditure, cost of electricity and cost of CO 2 capture in a 550-MW zeolite-membrane-reactor integrated IGCC process. The integration of the zeolite membrane reactor in IGCC could provide a significant reduction of 80% and 27% in the power consumption for Selexol™ Acid Gas Removal and CO 2 compression, respectively, which lowers the total auxiliary power consumption by 12.5%. However, the low pressure required at permeate stream for maintaining the driving force of hydrogen permeation through the zeolite membrane costs a huge power in permeate compressor, compensating the power consumption reduction achieved with the membrane reactor. Thus, for coal-fired IGCC for electricity generation with 90% CO 2 captured, the integration of the membrane reactor could provide a CO conversion ~99% and a significant drop in cost-of-electricity using zeolite membrane with H 2 permeance >600 GPU and the H 2 /CO 2 selectivity over 70.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Extended Modeling of DOE Sealed Canisters with Updated Chemistry Models

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. These models have coupled the temperature conditions to both the gas phase radiolysis chemistry and the surface chemistry associated with the oxyhydroxide layer. The previous iteration of the model utilized constant G-values for the hydrogen generation a 50-year period. This new iteration of the model utilizes new experimental data to update the hydrogen generation rate, as well as increase the simulated time to a 200-year period. Given the half-life assumed for the primary Cs-137 isotope responsible for gamma radiation in the ATR spent fuel, the 200-year period decreases the decay heat and dose rate of the spent fuel by a factor of 100, which combined with the updated chemistry substantially decreases new hydrogen generation. Continued experimental work has identified trends for aluminum surrogate samples with oxyhydroxide layers for tests done at higher dose rates. At low initial doses a fast generation rate of hydrogen occurs which starts to roll over to a lower generation as the total dose applied increase. A small-scale chemical model was built to replicate as mini-canister surrogate system at SRNL as well as for the smaller capsule tests performed at INL. A variety of chemical models to capture this effect were tested, and a back reaction of H radical absorbing onto the surface, or inhibition of the reaction by significant H 2 cover gas were both able to fit both the mini-canister and the small capsule test data. Both kinetic fits were used to generate data from the new 200-year simulation. As additional experiments continue, the kinetic fits may be adjusted to adapt to new data. However, updated data shows that hydrogen atmosphere has little effect on actual generation data, so the model was reverted to use a star-stepped G-value for low-dose and high-dose regions. For the undried fuel case, the three models differ significantly with an end concentration of 1% for back reaction, 4.2% for inhibition, and 13.1% for constant G-value for the nominal scenario. For the dried fuel case, the nominal cases showed end concentrations of 0.23% for back reaction, 1.9% for inhibition, and 4.2% for const G-value for the nominal scenario. For the constant G-value case that is less conservative than the other two, the pressure for undried fuel increases to 1.94 atm over 200 years for the nominal case and 2.11 atm for the high decay heat case. The total hydrogen concentration after 200 years is 4.92% for the low decay heat case and 27.4% for the high decay heat case for undried fuel and is 1.4% and 9.6% for low and high decay heat fuel for the dried fuel case. If small amount of residual air is present, the potential for nitric acid formation of 595, 1568, and 2816 ppm for the lower, nominal, and upper fuel decay heat can occur. At long-term timeframes no significant shift in major species present occurs, so no appreciable amount of oxygen is present in the system. The continued rate of hydrogen generation in the canister occurs at an increase of 0.04% by mole over the final 10 years. Increasing the model range out to 1000 years continues to drop the hydrogen generation rate through decreasing dose rate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

An experimental investigation of flame and autoignition behavior of propane

We report autoignition delay time data are one important means to develop, quantify, and validate fundamental understanding of combustion chemistry at low temperatures (T<1200 K). However, low-temperature chemistry often has higher uncertainties and scatter in the experimental data compared with high-temperature ignition data (T>1200 K). In this study, autoignition properties of propane and oxygen mixtures were investigated using the University of Michigan rapid compression facility in order to understand the effects of ignition regimes on low-temperature ignition data. For the first time for propane, autoignition delay times were determined from pressure histories, and autoignition characteristics were simultaneously recorded using highspeed imaging of the test section through a transparent end-wall. Propane mixtures with fuel-to-O 2 equivalence ratios of φ = 0.25 and φ = 0.5 and O 2 -to-inert gas molar ratios of 1:3.76 were studied over the pressure range of 8.9 to 11.3 atm and the temperature range of 930 – 1070 K. The results showed homogeneous or strong autoignition occurred for all φ = 0.25 experiments, and inhomogeneous or mixed autoignition occurred for all φ = 0.5 experiments. While a limited temperature range is covered in the study, importantly the data span predicted transitions in autoignition behavior, allowing validation of autoignition regime hypotheses. Specifically, the results agree well with strong-autoignition limits proposed based on the Sankaran Criterion. The autoignition delay time data at the strong-ignition conditions are in excellent agreement with predictions using a well-validated detailed reaction mechanism from the literature and a zerodimensional modeling assumption. However, the experimental data at the mixed autoignition conditions were systematically faster than the model predictions, particularly at lower temperatures (T< ~970 K). The results are an important addition to the growing body of data in the literature that show mixed autoignition phenomena are important sources of the higher scatter observed in the low-temperature autoignition data for propane and other fuels. The results are discussed in terms of different methods to capture the effects of pre-autoignition heat release associated with mixed autoignition conditions and thereby address some of the discrepancies between kinetic modeling and experimental measurements.

42 ENGINEERING↗
DOE OSTI2025

Emulation of the calculations of final r -process abundance patterns with a neural network

This work explores the construction of a fast emulator for the calculation of the final pattern of nucleosynthesis in the rapid neutron capture process (the r-process). An emulator is built using a feed-forward artificial neural network (ANN). We train the ANN with nuclear data and relative abundance patterns. We take as input the β-decay half-lives and the one-neutron separation energy of the nuclei in the rare-earth region. The output is the final isotopic abundance pattern. In this work, we focus on the nuclear data and abundance patterns in the rare-earth region to reduce the dimension of the input and output space. We show that the ANN can capture the effect of the changes in the nuclear physics inputs on the final r-process abundance pattern in the adopted astrophysical conditions. We employ the deep ensemble method to quantify the prediction uncertainty of the neural network emulator. The emulator achieves a speed-up by a factor of about 20 000 in obtaining a final abundance pattern in the rare-earth region. The emulator may be utilized in statistical analyses such as uncertainty quantification, inverse problems, and sensitivity analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Review and construction of interatomic potentials for molecular dynamics studies of hydrogen embrittlement in Fe─C based steels

Abstract Reducing hydrogen embrittlement in the low‐cost Fe─C based steels have the potential to significantly impact the development of hydrogen energy technologies. Molecular dynamics studies of hydrogen interactions with Fe─C steels provide fundamental information about the behavior of hydrogen at microstructural length scales, although such studies have not been performed due to the lack of an Fe─C─H ternary interatomic potential. In this work, the literature on interatomic potentials related to the Fe─C─H systems are reviewed with the aim of constructing an Fe─C─H potential from the published binary potentials. We found that Fe─C, Fe─H, and C─H bond order potentials exist and can be combined to construct an Fe─C─H ternary potential. Therefore, we constructed two such Fe─C─H potentials and demonstrate that these ternary potentials can reasonably capture hydrogen effects on deformation characteristics and deformation mechanisms for a variety of microstructural variations of the Fe─C steels, including martensite that results from γ to α phase transformation, and pearlite that results from the eutectic formation of the Fe 3 C cementite compound.

Zhou, Xiaowang↗
DOE OSTI2021

Quantification of energy and cost reduction from decreasing dissolved oxygen levels in full‐scale water resource recovery facilities

Abstract Aeration systems often lack the efficiency to maintain a desired residual dissolved oxygen (DO) concentration in the tank in part because little consideration is given to the dynamic daily and seasonal loading conditions. Although advanced aeration controllers exist, the majority of plants have DO set points typically based on common practice and literature values rather than site‐specific conditions, which can result in DO set points higher than those necessary to meet treatment objectives. DO set point reduction strategies have primarily been proposed through either static or dynamic simulations. In this study, the substantial improvements associated with DO set point reduction are demonstrated at full scale. A yearlong characterization of full‐scale aeration dynamics captured the effect of diurnal and seasonal fluctuations on oxygen transfer and energy demand and so facilitated the estimation of the potential savings of DO reduction strategies. Full‐scale validation provided direct evidence of DO reduction strategies inducing an overall enhancement of oxygen transfer efficiency along the different bioreactors, while confirming that energy savings as high as 20% were feasible. This study quantifies the influence of oxygen transfer efficiency on operating choices and site‐specific conditions (control strategy, loading conditions, and influent flow variability). Practitioner points We quantified the energy reduction and cost savings associated with a DO reduction in an aeration tank. For each 0.2 mg/L of DO decreased, the average power demand reduction per unit water treated exceeded 17%. Field measurements of dynamic alpha values eliminate the uncertainty in estimating aeration energy and cost savings from DO variations.

Pasini, Federico↗
DOE OSTI2021

Modeling creep behavior in ceramic matrix composites

Here, in this work, a three-dimensional viscoplasticity formulation with progressive damage is developed and used to investigate the complex time-dependent constituent load transfer and progressive damage behavior in ceramic matrix composites (CMCs) subjected to creep. The viscoplasticity formulation is based on Hill's orthotropic plastic potential, an associative flow rule, and the Norton-Bailey creep power law with Arrhenius temperature dependence. A fracture mechanics-informed isotropic matrix damage model is used to account for CMC brittle matrix damage initiation and propagation, in which two scalar damage variables capture the effects of matrix porosity as well as matrix property degradation due to matrix crack initiation and propagation. The Curtin progressive fiber damage model is utilized to simulate progressive fiber failure. The creep-damage formulation is subsequently implemented as a constitutive model in the generalized method of cells (GMC) micromechanics formulation to simulate time-dependent deformation and material damage under creep loading conditions. The developed framework is used to simulate creep of single fiber SiC/SiC microcomposites. Simulation results are in excellent agreement with experimental and numerical data available in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗