Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduce order modeling”

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

Impact of point spread function higher moments error on weak gravitational lensing – II. A comprehensive study

ABSTRACT Weak lensing is one of the most powerful probes for dark matter and dark energy science, although it faces increasing challenges in controlling systematic uncertainties as the statistical errors become smaller. The point spread function (PSF) needs to be precisely modeled to avoid systematic error on the weak lensing measurements. The weak lensing biases induced by errors in the PSF model second moments, i.e. its size and shape, are well-studied. However, Zhang et al. showed that errors in the higher moments of the PSF may also be a significant source of systematics for upcoming weak lensing surveys. Therefore, this work comprehensively investigate the modelling quality of PSF moments from the 3rd to 6th order, and propagate the PSFEx higher moments modelling error in the HSC survey data set to the weak lensing shear–shear correlation functions and their cosmological analyses. The overall multiplicative shear bias associated with errors in PSF higher moments can cause an ∼0.1σ shift on the cosmological parameters for LSST Y10, while the associated additive biases can induce 1σ uncertainties in cosmology parameter inference for LSST Y10, if not accounted. We compare the PSFEx model with PSF in Full FOV, and find similar performance in modelling the PSF higher moments. We conclude that PSF higher moment errors of the future PSF models should be reduced from those in current methods, otherwise needed to be explicitly modeled in the weak lensing analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On-the-fly response function generation method for composite coarse mesh

The hybrid stochastic deterministic transport code COMET, based on the incident response expansion theory, is used to model reactor cores with high fidelity and formidable computational speed. COMET models a reactor core using a library of incident flux response expansion coefficients that are pre computed for all the unique lattice cells (e.g., fuel assemblies, reflector blocks, etc.) in the core. In order to further improve its computational efficiency in pre-calculating the response library a new response function generation method is developed to compute the response functions for the composite coarse meshes made of a smaller set of unique lattices on the fly within the COMET's deterministic transport core sweep. The efficiency is achieved by eliminating a number of unique lattices that can be made up from the reduced set of unique meshes on the fly. The numerical process consists of the following steps. First, the boundary condition on composite coarse mesh boundaries is projected onto the expansion basis to compute the incident flux moments on external surfaces of all the basic (reduced set of unique) coarse meshes. Secondly, the deterministic sweeping solver in COMET is used to converge on the outgoing/incoming flux expansion moments crossing interfaces between the basic coarse meshes. Thirdly, the response functions for the composite coarse meshes are constructed as a superposition on the fly. The new response function generation method was tested on 88 composite coarse meshes consisting of CANDU fuel bundles and moderator blocks. It was found that response functions generated by the new method agree very well with those generated by direct Monte Carlo calculations. The average and maximum relative differences in the surface-to-surface response coefficients computed by the two methods are 0.10% and 0.20%, respectively. Similarly, the average and maximum relative differences in the response fission densities are 0.13% and 0.43%, respectively. These discrepancies are within one standard deviation of the stochastic uncertainties. The new method is five times faster than the original direct Monte Carlo method. The size of the response function library for the new method is five times smaller than that for the original method, leading to significantly less requirement for the computer hard drive space and memory. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Total kinetic energy release in the fast neutron induced fission of actinide nuclei

The total kinetic energy release and fission mass distributions for the fast neutron (En = 3–100 MeV) induced fission of 232 Th, 233 U, 235 U, 237 Np, 239 Pu, 240 Pu, and 242 Pu have been measured using the LANSCE facility. The neutron energies were deduced from time-of- flight measurements. The fission fragments were detected using Si PIN diode detectors, giving us the fragment energies. The actinide targets were made by vapor deposition leading to high-quality targets, that were thin and uniform with reduced impurities. Corrections were made to the data for pulse height defect and the fragment energy loss in the target and its backing. The TKE distributions were Gaussian in shape and their mean value as a function of incoming neutron energy could be fitted with second order polynomials. In the case of 233 U and 235 U, our measurements agree with prior work. Our measurements for 232 Th are unique. Our data agree with Viola scaling. The constant position of the heavy mass peak is interpreted as being due to the influence of the N = 88 and Z = 50 shells. The GEF model predictions agree with the data in general as do the CGMF model predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

New modelling capabilities in IDT

This work concerns the enhanced modelling capabilities of the discrete ordinates transport solver IDT. The novelties introduced allow for modelling unstructured geometries composed by a collection of X/Y segments and circles, and the use of reciprocity and conservation relations reduce the memory imprint as well as the computational cost of the method. IDT decomposes geometries in modular Cartesian patterns, which are the so-called Heterogeneous Cartesian Cells (HCCs), containing a chunk of the original unstructured geometries. Each HCC can be then discretized by superimposing a XY grid to refine locally the HCC. Unlike the most popular MOC, IDT performs the spatial sweeping by directional collision probabilities instead of trajectories. The sources and interface angular fluxes are expanded up to linear order. The accuracy of ray-tracing, the memory imprints together with the novel mesh refinement capabilities have been verified. A first set of preliminary results on PWR lattice problems will be presented in this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Experimental Characterization and Modeling of High Hole Mobility GeSn Quantum Wells: The Role of Alloy Disorder Scattering

Understanding mechanisms influencing electrical transport in material systems not only provides a scientific explanation for observed behavior but also offers insight into ways to enhance transport in devices. This study reports experimental hole mobility of 8 x 10 4 cm 2 V -1 s -1 in a Ge 0.92 Sn 0.08 , the highest recorded mobility for the GeSn system. A study of the material's quality is presented using structural and electrical characterization techniques, with transport data being supported by simulations using an extensive modeling framework. Quantum Hall measurements further indicate the material's high quality and potential spintronic applications, with extracted values of 0.0689$m$ 0 and 13.6 for the effective mass and effective g‐factor, respectively. It is observed that transport is limited by alloy disorder scattering at cryogenic temperatures. A comparative study between the presented structure and similar quantum well heterostructures revealed that the difference in hole mobilities is captured by a disparity in the reduced nominal alloy disorder scattering potential (Δ U alloy = 0.8 eV), that is lower than the value of a fully random alloy (Δ U alloy = 1.4–1.7 eV) potential. The difference in Δ U alloy suggests that heterostructures with similar geometries and alloy compositions can have different alloy disorder scattering, implying that an underlying mechanism, such as short‐range order, may be responsible and warrants further investigation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

MPEX AI Digital Twins

All magnetically confined plasma fusion power plant concepts (Tokamak, Spherical Tokamak, Stellarator, Mirror, ...) must exhaust the heat and plasma from the core confinement region to the material walls. The primary channel for this exhaust is through a plasma divertor which directs plasma along open magnetic field lines to a material target. The Material Plasma Exposure eXperiment (MPEX) illustrated in Figure 1, is a high-power, steady-state linear plasma device designed to produce the plasma material interaction (PMI) conditions of the divertor of future magnetic confinement fusion power plants: energy flux 20MW/m 2 , ion fluence 1031/m 2 , pulse duration 106 sec. These goals of plasma exposure in MPEX are well beyond those achieved in magnetic fusion experimental devices. Successfully achieving these high power steady state conditions for long pulses requires operational control of the heating and particle sources and the plasma flux to the walls and target. The MPEX AI Hot Spot Controller, proposed in this project, will help achieve the operational milestones of MPEX. The MPEX device will begin commissioning at the end of FY26. A smaller proto-MPEX was operated for 14,666 plasma discharges and will resume operation in September of 2025 as proto-MPEX-lite, with reduced capability, to test a new window for the Helicon plasma source. The proto-MPEX data has undergone surrogate modeling with machine learning methods (R. Archibald, 2022 IEEE International Conference on Big Data). This proto-MPEX data will be used to begin development of the AI digital twins described in this white paper. The scientific mission of MPEX is to qualify materials of different composition for use in the high energy and plasma flux conditions of a fusion power plant. The materials exposed in MPEX will in some cases be exposed to high neutron fluxes at other ORNL facilities to measure the changes to their PMI properties. The targets exposed in MPEX will be transported under vacuum to a Surface Analysis Station (SAS). The SAS will be equipped with the following diagnostics: Focused Ion Beam (FIB) for trench milling, 100-400 angstrom resolution scanning electron microscope (SEM), surface mapping x-ray spectrometer, high resolution camera, and a future upgrade to a laser induced breakdown spectroscopy quadruple mass spectrometer (LIBS-QMS). The MPEX experiments will generate diverse pre- and post-exposure measurement data of detailed material properties down to the crystal grain level in 3D for post-exposure assessment of PMI damage (e.g. cracking, melting, erosion and redeposition of the material). Physics models for the PMI, and how the material composition and manufacturing impact its performance under high energy plasma exposure, need to be validated with MPEX data to guide the selection of new candidate materials. Our vision for the MPEX AI Digital Twins project is to supply experimental and physics model simulation data to train Artificial Intelligence (AI) models for data processing, analysis, operational control, PMI and materials simulation to maximize the scientific output of the MPEX device. Ultimately, an AI digital twin of MPEX material assessment metrics for tested and synthetic material types with simulated PMI will be trained by the AI Modeling Teams on the experimental and physics simulation data submitted to the American Science Cloud by this project. A purely empirical search for the best material is inefficient given the finite number of samples that can be tested on MPEX. In order to expand the material properties database for training the MPEX Material Assessment AI Digital Twin, and to gain physics understanding of the PMI processes, physics models of the material properties and PMI processes are required. The physics simulations provide detailed simulation data, like impact angles for plasma ions, sputtering yields, transport of the ionized sputtered target material in the plasma, and redeposition locations. This simulation data expands the measurement data for deeper physics understanding. The experimental data is essential to validate the PMI and material structure simulation models. The validated models can then be used to generate new simulation data of MPEX material assessments for synthetic material compositions that have not been exposed in MPEX. These predictive simulations, plus the whole experimental dataset, will be used to train the MPEX Material Assessment AI Digital Twin allowing a rapid generative AI search for new materials with reduced PMI damage by interpolating the domain of the training set. These new optimum materials can be simulated with the physics codes and/or tested in MPEX. The ability of AI neural networks to interpolate multi-dimensional parameter spaces and generate virtual data is exploited for a more efficient search for optimum materials. The advent of the Transformational AI Models Consortium (TAIMC) is an opportunity to engage with state of the art private and public AI developers to achieve the goals of the AI digital twins and AI accelerated physics models proposed in this project. Our partners at ORNL from the Advance Scientific Computing Research (ASCR) organization will collaborate in accelerating the integrated plasma material interaction simulation framework. This simulation framework will provide a platform for generating simulation data across a range of physical fidelities, including hybrid methods that produce multi-fidelity results. This data will be leveraged for AI model development, both for generation of surrogates and the automation of simulation campaigns. A part of the research below will include collaborative efforts with the TAIMC to (i) adapt data storage approaches to ensure AI-readiness, (ii) provide a protypical exemplar to inform and exercise constructed workflows, and (iii) generate and share data, using the TAIMC unified AI data standard, for foundational models that will be trained from multiple sources across the DOE complex. We will also collaborate with the TAIMC, as well as the planned AI modeling teams, to develop approaches for reducing the cost of data generation. These include tailored multi-fidelity approaches as well as fine-tuning strategies to augment general, large-scale foundational models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Unraveling the CO Oxidation Mechanism over Highly Dispersed Pt Single Atom on Anatase TiO 2 (101)

Catalysts with noble metals deposited as single atoms on metal oxide supports have recently been studied extensively due to their maximized metal utilization and potential for performing difficult chemical conversions owing to their unique electronic properties. Understanding of the reaction mechanisms on supported single-metal atoms is still limited but is highly important for designing more efficient catalysts. In this study, we report the complexity of the CO oxidation reaction mechanism on Pt single atoms supported on anatase TiO 2 (Pt SA /a-TiO 2 ) by coupling density functional theory (DFT) calculations and microkinetic analysis with kinetic measurements, in situ/operando infrared, and X-ray absorption spectroscopies. Starting from the adsorbed Pt SA occupying an O vacancy induced by reductive pretreatment, we show that CO oxidation follows a complex mechanism consisting of initiation steps to reorganize the active site and multibranch reactive cycles, with the Pt SA /a-TiO 2 catalyst not returning to its initial configuration. The initiation step consists of CO and O 2 adsorption healing the O vacancy, followed by CO oxidation using gas-phase CO to form Pt(CO). The reactive cycle alternates O 2 adsorption and dissociation to oxidize the catalyst to Pt(O)(O)(CO) and branching pathways of competing Langmuir–Hinshelwood (LH)- or Eley–Rideal (ER)-type CO oxidation steps to reduce it again to Pt(CO). In situ/operando infrared experiments, including cryogenic CO adsorption and isotopic CO exchange, confirm the combined involvement of strongly adsorbed CO and gas-phase CO in an Eley–Rideal step along the reaction cycle. Microkinetic modeling shows that Pt single atoms are present in a mixture of Pt(CO), Pt(CO)(O 2 ), Pt(O)(CO)(O 2 ), and Pt(CO)(CO 3 ) structures as the main intermediates during steady-state CO oxidation, all having the C–O vibrational stretch close to the experimentally observed value of 2115 cm –1 . Microkinetic modeling also shows that the fractional orders of CO and O 2 measured experimentally originate from multiple steps with a high degree of rate control and not from a simple competitive adsorption. The results demonstrate the complex reaction pathways that even CO oxidation on a simple single-atom system can follow, providing mechanistic insights for designing efficient Pt-based single-atom catalysts. Here, we further show that microkinetic modeling results are sensitive to changes in energies of intermediate and transition states within errors of density functional theory, which can ultimately lead to incorrect conclusions regarding the reaction pathways and most abundant reaction intermediates if not accounted for by experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Local Chemical Clustering Enabled Ultrahigh Capacitive Energy Storage in Pb-Free Relaxors

Designing Pb-free relaxors with both a high capacitive energy density (W rec ) and high storage efficiency (η) remains a remarkable challenge for cutting-edge pulsed power technologies. Local compositional heterogeneity is crucial for achieving complex polar structure in solid solution relaxors, but its role in optimizing energy storage properties is often overlooked. Here, in this paper, we report that an exceptionally high W rec of 15.2 J cm –3 along with an ultrahigh η of 91% can be achieved through designing local chemical clustering in Bi 0.5 Na 0.5 TiO 3 –BaTiO 3 -based relaxors. A three-dimensional atomistic model derived from neutron/X-ray total scattering combined with reverse Monte Carlo method reveals the presence of subnanometer scale clustering of Bi, Na, and Ba, which host differentiated polar displacements, and confirming the prediction by density functional theory calculations. This leads to a polar state with small polar clusters and strong length and direction fluctuations in unit-cell polar vectors, thus manifesting improved high-field polarizability, steadily reduced hysteresis, and high breakdown strength macroscopically. The favorable polar structure features also result in a unique field-increased η, excellent stability, and superior discharge capacity. Our work demonstrates that the hidden local chemical order exerts a significant impact on the polarization characteristic of relaxors, and can be exploited for accessing superior energy storage performance.

25 ENERGY STORAGE↗
DOE OSTI2022

Geometric and doping effects on radiative recombination in thin-film near-field energy converters

Modeling radiative recombination is crucial to the analysis of radiative energy converters. In this work, a local radiative recombination coefficient is defined and derived based on fluctuational electrodynamics that is applicable to thin-film cells in both the near field and far field. The predicted radiative recombination coefficient of an InAs cell deviates from the van Roosbroeck–Shockley relation when the thickness is less than 10 µm, and the difference exceeds fourfold with a 10 nm film. The local radiative recombination coefficient is orders of magnitude higher when an InAs cell is configured in the near field. The local radiative recombination coefficient reduces as the doping level approaches that of a degenerate semiconductor. The maximum output power and efficiency of a thermoradiative cell would be apparently overpredicted if the electroluminescence coefficient defined in this paper were taken as unity for heavily doped semiconductors.

Feng, Dudong↗
DOE OSTI2024

End-Use Savings Shapes Measure Documentation: Window Film

This documentation focuses on a single end-use savings shape measure - window film. The window film studied in this analysis, called solar control film, is a passive retrofit solution for windows that does not involve window replacement. This type of film is composed of transparent, tinted, or metalized laminated polyester layers and can be attached to an existing window surface (either on the exterior or interior side of the window). The properties of the window film are designed to shift the thermal and optical performances of the overall glazing system in order to serve various needs the customer would have (e.g., heat, glare). While the practical goal of purchasing and installing a window film varies widely in the real market, this study only focuses on the goal of energy savings, and highlighting the corresponding emissions. Other important aspects that customers typically consider include visual comfort, privacy, aesthetics, ultraviolet protection, etc. Thus, in practice, customers often choose a window film product not only to save energy (or cost) but also to mitigate issues around glare, excessive light, daytime privacy, or inconsistent appearance of the building. Window film products that were modeled in this analysis significantly reduced the solar heat gain coefficient of the overall glazing system, resulting in better energy savings for buildings in hot climate regions. However, significantly reducing the solar heat gain coefficient that can block unfavorable heat during the summer can actually harm blocking favorable heat during the winter. By applying window films on a stock of buildings covering various load and weather conditions, this analysis highlights when (e.g., time of day) and where (e.g., geospacial location) we can save energy with window films.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

CUPID: The Next-Generation Neutrinoless Double Beta Decay Experiment

We report CUPID is a next-generation tonne-scale bolometric neutrinoless double beta decay experiment that will probe the Majorana nature of neutrinos and discover lepton number violation in case of observation of this singular process. CUPID will be built on experience, expertise and lessons learned in CUORE and will be installed in the current CUORE infra-structure in the Gran Sasso underground laboratory. The CUPID detector technology, successfully tested in the CUPID-Mo experiment, is based on scintillating bolometers of Li 2 MoO 4 enriched in the isotope of interest 100 Mo. In order to achieve its ambitious science goals, the CUPID collaboration aims to reduce the backgrounds in the region of interest by a factor 100 with respect to CUORE. This performance will be achieved by introducing the high efficient α/β discrimination demonstrated by the CUPID-0 and CUPID-Mo experiments, and using a high transition energy double beta decay nucleus such as 100 Mo to minimize the impact of the gamma background. CUPID will consist of about 1500 hybrid heat-light detectors for a total isotope mass of 250 kg. The CUPID scientific reach is supported by a detailed and safe background model based on CUORE, CUPID-Mo and CUPID-0 results. The required performances have already been demonstrated and will be presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

High-Efficiency Inductive Output Tubes Using a Third Harmonic Drive on the Grid

In this article, we discuss the use of a third harmonic component to the drive voltage on the grid of an inductive output tube (IOT). High-efficiency IOTs are typically characterized by efficiencies up to 70%–75%. However, the achievement of efficiencies greater than 80% would substantially reduce the operating costs of next-generation accelerators. In order to achieve this goal, we consider the addition of a third harmonic component to the drive signal on the grid of an IOT. The use of a third harmonic drive component in IOT guns has been considered in order to apply such a gun as the injector of radio frequency linear accelerators (RF linacs). Furthermore, we consider that the IOT will be used to provide the rf power to drive RF linacs and apply the third harmonic with the intention of increasing the efficiency of the RF output of the IOT. We consider a model IOT with a 700-MHz resonant cavity and using an annular beam with a voltage of 30 kV, an average current of 6.67 A yielding a perveance of about 1.3 μP. We simulate this IOT using the NEMESIS simulation code which has been successfully validated by comparison with the K5H90W-2 IOT developed by Communications and Power Industries. It is found that the effect of the third harmonic on the efficiency is greatest when the phase of the third harmonic is shifted by radians with respect to the fundamental drive signal and with third harmonic powers greater than about 50% that of the fundamental drive power. For the present example, we show that efficiencies approaching 86% are possible by this means.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

The Mu2e experiment — Searching for charged lepton flavor violation

The Mu2e experiment will search for a Standard Model violating rate of neutrinoless conversion of a muon into an electron in the presence of an aluminum nucleus. Observation of this charged lepton flavor violating process would be an unambiguous sign of new physics. Mu2e will improve upon previous searches for this process by four orders of magnitude. This requires the world’s highest-intensity muon beam, a detector system capable of efficiently reconstructing the 105 MeV/c conversion electron signal, and minimizing sensitivity to background events. A pulsed 8 GeV proton beam strikes a target, producing pions that decay into muons. Beam outside the pulse must be suppressed to < 1 0 – 10 to reduce beam-related backgrounds. The muon beam is guided from the production target along the transport system and onto the aluminum stopping target. Conversion electrons leave the stopping target and propagate inside a solenoidal magnetic field to the tracker and electromagnetic calorimeter. The tracker is a system of straw tube panels filled with Ar/CO 2 at 1 atm that tracks particles inside of a solenoidal B-field and measures their momenta with ~ 100 keV/c resolution to resolve signal events from decay-in-orbit backgrounds. The CsI calorimeter provides E / p and is used to seed the track reconstruction algorithm with σ E / E ~ 10 % and σ t < 500 ps . Additionally, a novel cosmic ray veto with greater than 99.99% efficiency brings the expected number of background events to fewer than one over three years of running. Finally, to normalize the experiment, the stopping target monitor measures the rate of capture photons from muons incident on the stopping target by using a system of high-purity germanium and lanthanum bromide scintillators.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Aeroservoelastic Stability of a Floating Wind Turbine

The problem of negative damping undermines the aeroservoelastic stability of floating offshore wind turbines. The negative damping problem is most prevalent around rated wind speed, where the sensitivity of thrust to wind speed is the largest. This paper investigates the implementation of peak shaving, a controller feature that limits the rated thrust by pitching the blades before rated wind speed is reached. Two controller designs are investigated: a de-tuned controller and a nacelle-feedback controller. A time-domain metric is defined, inspired by Lyapunov theory, in order to compute and assess the stability of floating offshore wind turbines. The model of the International Energy Agency 15-MW reference wind turbine mounted on the University of Maine VolturnUS-S floater is simulated in HAWC2 with the National Renewable Energy Laboratory reference open-source controller. Peak shaving is applied to the two controller designs and stability is assessed. According to the chosen metric, peak shaving does not improve the stability of the system. This is due to the trade-off between loads and error tracking: although the loads and displacements in the fore-aft direction of the turbine are reduced, the rotor-speed tracking is poorer, which increases the shaft torsion fatigue load.

17 WIND ENERGY↗
DOE OSTI2024

Once-in-a-lifetime encounter models for neutrino media: From coherent oscillations to flavor equilibration

Collective neutrino oscillations are typically studied using the lowest-order quantum kinetic equation, also known as the mean-field approximation. However, some recent quantum many-body simulations suggest that quantum entanglement among neutrinos may be important and may result in flavor equilibration of the neutrino gas. In this work, we develop new quantum models for neutrino gases in which any pair of neutrinos can interact at most once in their lifetimes. A key parameter of our models is γ = μ Δ z , where μ is the neutrino coupling strength, which is proportional to the neutrino density, and Δ z is the duration over which a pair of neutrinos can interact each time. Our models reduce to the mean-field approach in the limit γ → 0 and achieve flavor equilibration in time t ≫ ( γ μ ) − 1 . These models demonstrate the emergence of coherent flavor oscillations from the particle perspective and may help elucidate the role of quantum entanglement in collective neutrino oscillations. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗
DOE OSTI2020

Effect of Cr(III) Adsorption on the Dissolution of Boehmite Nanoparticles in Caustic Solution

Incorporation of relatively minor impurity metals onto metal (oxy)hydroxides can strongly impact solubility. In complex highly alkaline multicomponent tank waste such as that at the Hanford Nuclear Reservation, tests indicate that the surface area-normalized dissolution rate of boehmite ($γ$-AlOOH) nanomaterials is at least an order of magnitude lower than that predicted for the pure phase. In this work we examine the dissolution kinetics of boehmite coated by adsorbed Cr(III), which adheres at saturation coverages as sparse chemisorbed monolayer clusters. Using 40 nm boehmite nanoplates as a model system, temperature-dependent dissolution rates of pure versus Cr(III)-adsorbed boehmite showed that the initial rate for the latter is consistently several times lower, with an apparent activation energy 16 kJ·mol -1 higher. Although the surface coverage is only around 50%, solution analytics coupled to multimethod solids characterization reveals a phyicochemical armoring effect by adsorbed Cr(III) that substantially reduces the number of dissolution-active sites on particle surfaces. Such findings could help improve kinetics models of boehmite and/or metal ion adsorbed boehmite nanomaterials, ultimately providing a stronger foundation for the development of more robust complex radioactive liquid waste processing strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Revealing the Crystalline Architecture of Semicrystalline Ion Exchange Membranes for the Design of Conductive and Durable Alkaline Anion Exchange Membranes

Alkaline anion exchange membrane (AAEM) fuel cells offer a cost-effective alternative to proton exchange membrane (PEM) fuel cells by eliminating the need for expensive precious metal catalysts. In both PEMs and AAEMs, semicrystalline polymers are a common choice, as the crystalline domains can act as mechanical reinforcements that limit swelling and promote mechanical durability in the material. However, spatially resolved characterization of crystalline organization in ion exchange membranes beyond ensemble-averaged X-ray scattering is underrepresented, likely in part due to ionization damage limitations in soft materials. Here, in this study, we resolve the nanometer-size crystallites in semicrystalline ion exchange membranes by applying cryogenic four-dimensional scanning transmission electron microscopy (cryo-4D-STEM) along with data-processing algorithms designed to optimize signals at a low dose to minimize radiation damage. We investigate the effects of synthesis components, including molecular weight and thermal treatment, on a model system of AAEMs in comparison to Nafion, the most commonly used and commercially successful PEM today. We find that excess water uptake in polymer membranes, a property directly associated with weak mechanical durability and with possible negative impacts on ion conductivity, can be reduced by over 30% by varying the polymer's crystalline morphology through changes in synthesis parameters such as molecular weight and thermal history. Our results indicate that this improvement is correlated with smaller crystalline domains with a more homogeneous distribution. More broadly, these results demonstrate how the crystalline architecture of polymer membranes can be tuned through their chemistry and thermal treatment in order to improve their conductivity and durability for commercial fuel cell performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗