Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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,315 records · Page 73

NH 4 OH Looping with Membrane CO 2 Absorber and Distributed Stripper for Enhanced Algae Growth

The University of Kentucky Center for Applied Energy (UK CAER) has devised a unique, integrated CO2 capture and utilization technology. CO2 from coal-fired power generation flue gas is first captured at half the operating cost of a typical aqueous CO2 capture system (CCS), distributed in an aqueous stream and then fixed by algae in bioreactors where the algae production is increased by 50% over that with a typical intermittent nutrient feeding system. Lower CCS operating cost is achieved by eliminating the flue gas pretreatment step for cooling and SO2 removal, eliminating steam extraction from the power generation steam cycle for solvent regeneration, and eliminating CO2 compression. Higher algae production is achieved by continuous, just-in-time nutrient feed to the bioreactors directly from a distributed solvent regenerator, which maintains the bioreactor pH for optimum growth. The process starts with a uniquely configured membrane absorber, where the flue gas is indirectly contacted with an ammonium hydroxide (NH4OH) solvent. Dissolved NH3 is attractive for both CO2 capture and as an algae nutrient. For CO2 capture it is inexpensive, has a low regeneration energy, is thermally- and oxidatively-stable and has a viscosity near that of water, which makes is easy to transport. Numerous studies have shown that the scrubbing capacity of NH3 is approximately 0.9-1.2 kg of CO2/kg of NH3, with a CO2 removal efficiency of ~99% and half the solvent regeneration energy than that of 30 wt% MEA[1, 2, 3]. NH3 is attractive as an algae nutrient due to its low cost. The rich NH4OH solvent is pumped to a set of distributed regenerators which are co-located with the algae bioreactors. Solvent pumping, transport and distribution reduces the balance of plant (BOP) cost compared to a typical aqueous CCS related to the flue gas duct and boost fan required to transport the flue gas. The energy required for the distributed solvent regeneration is supplied by solar-thermal panels eliminating the need for steam extraction from the power generation steam cycle. After solvent regeneration, the product stream contains both the CO2 captured from the flue gas and volatized NH3 from the solvent. This product stream is fed directly to the bioreactors, eliminating the need for compression of the CO2 stream. The relative amounts of CO2 and NH3 in the product stream are adjusted and controlled by a controlling the regeneration conditions (pressure and temperature). The continuous feed of the right ratio of nutrients overcomes the typical inhibition of algae growth resulting from frequent pH swings in the bioreactor due to unbalanced (intermittent) feeding systems for CO2 and N. Also, because the regenerators will operate at pressure and be located in close proximity to the bioreactors, there is no worry about pressure drop when sparging the gas into the algae. Sparging produces small bubbles which is beneficial for mass transfer efficiency. One known challenge when using an NH4OH solvent is high NH3 emission. Hydrophobic membranes are used for CO2 capture using an aqueous NH3 solution[4, 5] without the direct contact between flue gas and aqueous solution. Additionally, UK CAER CO2 capture and utilization process manages NH3 slip in three extra measures. First, NH3 slip is minimized by working with minimal species partial pressure, which is proportional to the concentration in the liquid. Hence, lowering the capture solvent concentration will lower the NH3 partial pressure. Second, UK CAER’s previous work has demonstrated that the addition of Zn2+ into NH3 solutions to chelate the NH3 can reduce NH3 volatility. Third, the configuration of the membrane CO2 absorber utilizes condensed water from the flue gas to continually wash the gas-side of the membrane to reduce fouling and recapture NH3 slip. Additional details about the UK CAER unique, integrated CO2 capture and utilization technology will be presented along with technology development plans. Diao, N., Q. Li, and Z. Fang. 2004. Heat transfer in ground heat exchangers with groundwater advection. International Journal of Thermal Sciences. 43: 1203-1211, He, Q., M. Chen, L. Meng, K. Liu, and W. Pan. 2004. Study on Carbon Dioxide Removal from Flue Gas by Absorption of Aqueous Ammonia. Western Kentucky University. Yeh, A.C., and H. Bai. 1999. Comparison of ammonia and monoethanolamine solvents to reduce CO2 greenhouse gas emissions. The Science of the Total Environment. 228: 121-133, Villeneuve, K., D. Roizard, J.C. Remigy, M. Iacono, and S. Rode. 2018. CO2 capture by aqueous ammonia with hollow fiber membrane contactors: Gas phase reactions and performance stability. Separation and Purification Technology, 199: 189-197, Toro Molina, C., and C. Bouallou. 2016. Carbon dioxide absorption by ammonia intensified with membrane contactors. Clean Techn Environ Policy 18, 2133–2146 (2016)

20 FOSSIL-FUELED POWER PLANTS↗

Recent Developments in Security-Constrained AC Optimal Power Flow: Overview of Challenge 1 in the ARPA-E Grid Optimization Competition

In “Recent Developments in Security-Constrained AC Optimal Power Flow: Overview of Challenge 1 in the ARPA-E Grid Optimization Competition,” we review the state of the art in practical algorithms for scheduling power-systems operations in the short term and the results of the recent competition organized by the U.S. Advanced Research Projects Agency–Energy. We explain the mixed-integer nonlinear formulation used in the competition for nonspecialists in electrical engineering, the context and organization of the competition, and the performance of competitors. We find that the collective approaches and results of competitors provide support for efforts to move nonlinear optimization techniques into industrial applications, as they have proven to be a robust and efficient alternative to current linear approximation techniques.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Josephson Traveling Wave Parametric Amplifiers to Enable Future Neutrino Mass Measurements

Technologies for “quantum sensing" have the promise of dramatically improved sensitivity in some experiments that explore nuclear physics at the boundaries of our current knowledge. In particular, the Josephson Traveling Wave Parametric Amplifier (JTWPA) could enable future experiments such as Project 8 that use Cyclotron Radiation Emission Spectroscopy (CRES) to measure the absolute neutrino mass scale in tritium beta decays. CRES consists of detecting a signal of approximately one femtowatt of power radiated by electrons trapped in a 1-T magnetic field. The frequency of that radiation, around 27 GHz for a tritium endpoint electron, is inversely proportional to the electron energy. One of the significant challenges of a CRES experiment is the detection of minimal signal power levels. State-of-the-art transistor-based technologies will not likely achieve sufficient signal-to-noise ratio (SNR) in a competitive neutrino mass experiment. However, emerging quantum technologies have been demonstrated to operate at the fundamental lower limits of noise permitted by the laws of quantum mechanics. Of those technologies, JTWPAs appear to satisfy future performance needs with speed gain, bandwidth, and very low noise. We propose to bring together a group of neutrino physicists, experts in the needs of future neutrino mass measurements, with leaders in the field of quantum information science and JTWPA technology to develop JTWPAs appropriate for neutrino mass measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessing Options for Particulate Measurement Tools to Advance Wood Heaters (CRADA Final Report)

For approximately 40 years, particulate matter (PM) from residential wood heaters has been certified in North America using an EPA dilution tunnel with a filter-based gravimetric measurement. While reasonably repeatable, this approach reports only a single value for the entire burn sequence and does not provide information on PM at different points within the sequence, which would support improved wood heater designs. PM monitors that report with high time resolution offer a potential solution, but their performance for various types of heaters and across heater operating conditions has not been extensively evaluated. In this 18-month collaborative project, the Hearth, Patio & Barbecue Association (HPBA) and Lawrence Berkeley National Laboratory (LBNL) tested five continuous PM instruments spanning three measurement principles—light-scattering, inertial mass (TEOM), and triboelectric—against gravimetric and size-resolved references on three residential heater types (cordwood, catalytic cordwood, and pellet) across startup, high-burn, and low-burn operating phases.

42 ENGINEERING↗

(Oxy)hydroxides Formed on Aluminum Fuel Materials After Irradiation and Long-Term Wet Storage - 20354

The aluminum cladding of research-reactor fuel experiences general corrosion when in contact with water during in-reactor service and post-discharge wet storage, resulting in the formation of adherent aluminum (oxy)hydroxide films. These (oxy)hydroxides contain chemically-bound water that poses challenges for extended dry storage due to the risk of thermal or radiolytic decomposition releasing free water and/or hydrogen and oxygen gases. This study describes characterization of the (oxy)hydroxides present on several aluminum materials used in reactor operation and subsequently stored wet in the L-Basin storage facility at the Savannah River Site (SRS) for an extended period. Characterization data providing insight into the loading, composition, and morphology of (oxy)hydroxides to be expected on service-exposed aluminum cladding provides valuable benchmarks for designing adequate drying and dry-storage approaches. This work is part of a broader investigation to address knowledge gaps and technical data needs for dry storage of aluminum-clad spent nuclear fuel (ASNF), which included in-lab growth of (oxy)hydroxide films on aluminum alloy substrates to investigate formation behavior, investigation of drying methods to remove existing (oxy)hydroxides from ASNF cladding, and measurement of radiolytic yield of hydrogen from (oxy)hydroxide powders and films. In this study, (oxy)hydroxide films were characterized for three aluminum-alloy materials used in reactors and subsequently stored wet for up to approximately 40 years in L-Basin at SRS: one cropping from a Missouri University Research Reactor (MURR) fuel element (Al-6061 alloy), one cropping from a Universal Sleeve Housing (USH) (Al-6063 alloy), and one Mark-16B fuel assembly (either Al-6061 or Al- 6063). The USH and Mark-16B were used in SRS production reactors. Characterization of the as-received (oxy)hydroxides included scanning electron microscopy (SEM) in both plan-view and cross-section to characterize the (oxy)hydroxide layer's morphology, thickness, and structure. X-ray diffraction (XRD) was used to identify the chemical composition and distinguish between the various aluminum (oxy)hydroxides known to form under reactor and storage conditions. XRD analysis revealed both bayerite (Al(OH)3) and boehmite (AlOOH) on the surface of the MURR and USH samples, as well as bayerite, boehmite, and gibbsite (another Al(OH)3 polymorph) on the surface of the Mark-16B sample. The aluminum trihydroxides, bayerite and gibbsite, are typically associated with corrosion in low-temperature (<80 deg. C) water, while boehmite is expected to form at higher water temperature (>80 deg. C). The presence of bayerite on the USH, which is believed to have operated close to 90 deg. C, suggests that boehmite formed during in-reactor exposure was not protective against further hydroxide growth in low-temperature wet storage. Cross-section scanning electron microscopy (SEM) showed total (oxy)hydroxide layer thicknesses of ∼5- 10 μm for the MURR and ∼5-15 μm for the Mark-16B. The thickness of the USH's (oxy)hydroxide layer was indiscernible by the current mounting and imaging method, despite plan-view SEM and XRD confirming the presence of an (oxy)hydroxide layer. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Converse to Lieb–Robinson Bounds in One Dimension Using Index Theory

Abstract Unitary dynamics with a strict causal cone (or “light cone”) have been studied extensively, under the name of quantum cellular automata (QCAs). In particular, QCAs in one dimension have been completely classified by an index theory. Physical systems often exhibit only approximate causal cones; Hamiltonian evolutions on the lattice satisfy Lieb–Robinson bounds rather than strict locality. This motivates us to study approximately locality preserving unitaries (ALPUs). We show that the index theory is robust and completely extends to one-dimensional ALPUs. As a consequence, we achieve a converse to the Lieb–Robinson bounds: any ALPU of index zero can be exactly generated by some time-dependent, quasi-local Hamiltonian in constant time. For the special case of finite chains with open boundaries, any unitary satisfying the Lieb–Robinson bound may be generated by such a Hamiltonian. We also discuss some results on the stability of operator algebras which may be of independent interest.

Physics↗

Statistical fracture behavior of doped UO 2 using a ball-on-ring equibiaxial flexure test method

Metal oxide dopants, such as titanium and chromium oxides, have garnered considerable attention for their potential to increase grain size (≥ 30 µm) in UO 2 fuel, purportedly enhancing fission gas retention during reactor operation. Fuel performance is significantly impacted by fuel fracture behavior, so it is important to understand the effects of enhanced grain size and dopant content on UO 2 fuel fracture. UO 2 pellets were doped with 0.1 wt% TiO 2 and 0.3 wt% Cr 2 O 3 to alter density and grain size. Inductively coupled plasma mass spectroscopy measured dopant levels pre- and post-sintering. X-ray diffraction revealed lattice changes and microstrain via Rietveld refinement. Field emission scanning electron microscopy determined grain sizes of approximately 30 µm for TiO 2 doping and 7 µm for Cr 2 O 3 doping. Transverse rupture strength tests were performed on over 30 samples per dataset to obtain characteristic strength and Weibull modulus. Results indicate no statistical difference in fracture strength between 0.1 wt% TiO 2 doped UO 2 and undoped UO 2 , while 0.3 wt% Cr 2 O 3 doped UO 2 exhibited a 20% decrease in fracture strength. Doped UO 2 samples also showed reduced Weibull modulus compared to undoped UO 2 , suggesting increased scatter in fracture strength. This study's findings suggest that titanium and chromium oxide doping in UO 2 , regardless of grain size, induce residual stresses, decreasing fracture strength and increasing variability in fracture behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Automated co-adding and energy calibration of large array microcalorimeter data with zero sample knowledge

State-of-the-art microcalorimeter spectrometers now contain large detector arrays with hundreds of individual pixels. Each individual pixel outputs a unique and non-linear response with respect to deposited energy. This work describes a pattern-recognition algorithm to combine these responses into a single energy-calibrated histogram, referred to as co-adding pixels. Photo-peaks from different pixels are matched together based upon how well the match aligns the centroids and heights of neighboring peaks. This usually results in around 100 co-adding calibration points from 30 to 300 keV for a several day acquisition of plutonium items with masses between 0.5 and 10 grams. An additional algorithm energy-calibrates this co-added spectrum using the fluoresced K x-ray emissions from a tantalum absorber and inherent x-ray escape peaks from the tin absorbers. Both algorithms operate without knowledge of the source and are fully automated. This work presents results from the acquisitions of high and low burnup plutonium, 10% enriched uranium, a 153 Gd calibration source, and a 57 Co+ 166m Ho calibration source. In all measurements, resolution defined as the full-width at half-maximum (FWHM) of photo-peaks is preserved between the individual pixel and co-added spectra at around 65 eV for incident photon energies between 60 and 208 keV. The energy calibration algorithm is approximate and yields a calibration curve off by an average of around 200 eV for incident photon energies between 60 and 208 keV.

47 OTHER INSTRUMENTATION↗

Characterisation of divertor detachment onset in JET-ILW hydrogen, deuterium, tritium and deuterium–tritium low-confinement mode plasmas

Measurements of the ion currents to and plasma conditions at the low-field side (LFS) divertor target plate in low-confinement mode plasmas in the JET ITER-like wall materials configuration show that the core plasma density required to detach the LFS divertor plasma is independent of the hydrogenic species protium, deuterium and tritium, and a 40 %/60 % deuterium–tritium mixture. This observation applies to a divertor plasma configuration with the LFS strike line connected to the horizontal part of the LFS divertor chosen because of its superior diagnostic coverage. The finding is independent of the operational status of the JET cryogenic pump. The electron temperature (T e ) at the LFS strike line was markedly reduced from 25 eV to 5 eV over a narrow range of increasing core plasma density, and observed to be between 2 eV and 3 eV at the onset of detachment. The electron density (n e ) peaks across the LFS plasma when T e at the target plate is 1 eV, and spatially moves to the X-point for higher core densities. The density limit was found approximately 20 % higher in protium than in tritium and deuterium–tritium plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CFL Optimized Forward–Backward Runge–Kutta Schemes for the Shallow-Water Equations

Abstract We present the formulation and optimization of a Runge–Kutta-type time-stepping scheme for solving the shallow-water equations, aimed at substantially increasing the effective allowable time step over that of comparable methods. This scheme, called FB-RK(3,2), uses weighted forward–backward averaging of thickness data to advance the momentum equation. The weights for this averaging are chosen with an optimization process that employs a von Neumann–type analysis, ensuring that the weights maximize the admittable Courant number. Through a simplified local truncation error analysis and numerical experiments, we show that the method is at least second-order in time for any choice of weights and exhibits low dispersion and dissipation errors for well-resolved waves. Further, we show that an optimized FB-RK(3,2) can take time steps up to 2.8 times as large as a popular three-stage, third-order strong stability-preserving Runge–Kutta method in a quasi-linear test case. In fully nonlinear shallow-water test cases relevant to oceanic and atmospheric flows, FB-RK(3,2) outperforms SSPRK3 in admittable time step by factors roughly between 1.6 and 2.2, making the scheme approximately twice as computationally efficient with little to no effect on solution quality. Significance Statement The purpose of this work is to develop and optimize time-stepping schemes for models relevant to oceanic and atmospheric flows. Specifically, for the shallow-water equations we optimize for schemes that can take time steps as large as possible while retaining solution quality. We find that our optimized schemes can take time steps between 1.6 and 2.2 times larger than schemes that cost the same number of floating point operations, translating directly to a corresponding speedup. Our ultimate goal is to use these schemes in climate-scale simulations.

54 ENVIRONMENTAL SCIENCES↗

Atmospheric Drivers of Wind Turbine Blade Leading Edge Erosion: Review and Recommendations for Future Research

Leading edge erosion (LEE) of wind turbine blades causes decreased aerodynamic performance leading to lower power production and revenue and increased operations and maintenance costs. LEE is caused primarily by materials stresses when hydrometeors (rain and hail) impact on rotating blades. The kinetic energy transferred by these impacts is a function of the precipitation intensity, droplet size distributions (DSD), hydrometeor phase and the wind turbine rotational speed which in turn depends on the wind speed at hub-height. Hence, there is a need to better understand the hydrometeor properties and the joint probability distributions of precipitation and wind speeds at prospective and operating wind farms in order to quantify the potential for LEE and the financial efficacy of LEE mitigation measures. However, there are relatively few observational datasets of hydrometeor DSD available for such locations. Here, we analyze six observational datasets from spatially dispersed locations and compare them with existing literature and assumed DSD used in laboratory experiments of material fatigue. We show that the so-called Best DSD being recommended for use in whirling arm experiments does not represent the observational data. Neither does the Marshall Palmer approximation. We also use these data to derive and compare joint probability distributions of drivers of LEE; precipitation intensity (and phase) and wind speed. We further review and summarize observational metrologies for hydrometeor DSD, provide information regarding measurement uncertainty in the parameters of critical importance to kinetic energy transfer and closure of data sets from different instruments. A series of recommendations are made about research needed to evolve towards the required fidelity for a priori estimates of LEE potential.

17 WIND ENERGY↗

Core Physics Characteristics of Extended Enrichment and High Burnup Boiling Water Reactor Fuel

This paper presents the highlights of boiling water reactor (BWR) core physics studies performed at Oak Ridge National Laboratory as part of a series of studies conducted to compare low-enriched uranium (LEU) with LEU+ fuel. The studies analyzed isotopic fuel content, lattice parameters (Phase 1), and core physics (Phase 2) to identify challenges in operation, storage, and transportation for BWRs and pressurized water reactors (PWRs). Because of a lack of publicly available lattice and core designs for modern BWR fuel assemblies and reactor cores, several optimized lattice designs were generated, and different core loading strategies were investigated. Twelve optimized lattice designs with 235 U enrichments ranging from 1.6% to 9% and gadolinia loadings ranging from 3 to 8 wt% were used to model axial enrichment and geometry variations in fuel assemblies for core designs. Each core shares a common set of approximations in design and analysis to allow for consistent comparisons between LEU and LEU+ fuel. The objective is to highlight anticipated changes in core behavior with respect to the reference LEU core. The results of this study show that the differences in LEU and LEU+ core reactor physics characteristics are less significant than the differences in lattice physics characteristics reported in the Phase 1 studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Non-intrusive reduced-order modeling for dynamical systems with spatially localized features

This work presents a non-intrusive reduced-order modeling framework for dynamical systems with spatially localized features characterized by slow singular value decay. The proposed approach builds upon two existing methodologies for reduced and full-order non-intrusive modeling, namely Operator Inference (OpInf) and sparse Full-Order Model (sFOM) inference. We decompose the domain into two complementary subdomains that exhibit fast and slow singular value decay. The dynamics of the subdomain exhibiting slow singular value decay are learned with sFOM while the dynamics with intrinsically low dimensionality on the complementary subdomain are learned with OpInf. The resulting, coupled OpInf-sFOM formulation leverages the computational efficiency of OpInf and the high resolution of sFOM, and thus enables fast non-intrusive predictions for conditions beyond those sampled in the training data set. A novel regularization technique with a closed-form solution based on the Gershgorin disk theorem is introduced to promote stable sFOM and OpInf models. We also provide a data-driven indicator for subdomain selection and ensure solution smoothness over the interface via a post-processing interpolation step. We evaluate the efficiency of the approach in terms of offline and online speedup through a quantitative, parametric computational cost analysis. We demonstrate the coupled OpInf-sFOM formulation for two test cases: a one-dimensional Burgers’ model for which accurate predictions beyond the span of the training snapshots are presented, and a two-dimensional parametric model for the Pine Island Glacier ice thickness dynamics, for which the OpInf-sFOM model achieves an average prediction error on the order of 1% with an online speedup factor of approximately 8$\times$ compared to the numerical simulation.

42 ENGINEERING↗

Additive manufacturing process parameter determination for a new Fe-C-Cu alloy

As a potential replacement for stainless steel alloys commonly used to print parts with laser powder bed fusion, a new Cu precipitation strengthened ferrous alloy, with composition Fe-0.2C-6Cu (wt%), was recently developed. This material is Co- and Ni-free, printable, has mechanical properties comparable to that of high strength stainless steels (approximately 1300 MPa UTS), and is cost-advantaged relative to existing low-alloy steels for additive manufacturing of parts with complex shapes. To gain traction for broader application, optimal laser powder bed fusion (LPBF) processing conditions to produce nearly defect-free parts without compromising mechanical strength are needed. Here, an optimal processing window based upon laser speeds that achieve minimal porosity was quantified via comparisons of measured melt pool penetration, scan speeds, printed material density and laser beam energy density on a commercial LPBF unit operating at 350 W, 80 μm hatch spacing, and 50 μm layer thickness based upon optimal parameters for 17-4PH. Furthermore, the window is 300 to 500 mm/s to achieve >99.8 % dense Fe-0.2C-6Cu (wt%) specimens. Other defects such as burning, spatter, and lack of fusion were also avoided in this window. The window was validated with in-situ synchrotron X-ray imaging that enabled visualization of the vapor cavity, melting, and solidification during a single laser track scan on a miniature powder bed sample. In addition, in-situ infrared imaging provided temperature fields, cooling rate and solidification range, and confirmed minimal printing defects and spatter within the processing window.

36 MATERIALS SCIENCE↗

Sensitivity Analysis and Effective Parametrization of PEM Fuel Cell Models

The cost of proton-exchange-membrane fuel cells (PEMFCs) remains a major hurdle in large-scale commercialization of this technology. To improve their performance and reduce cost, novel materials and electrode designs are continuously envisioned, e.g., non-PGM catalyst layers, ultra-thin Pt/Pt-Ni based catalyst layers, structured ionomer arrays or NSTF catalyst layers.1 Understanding the impact of these improvement strategies can be extremely time and cost intensive due to complex physical phenomena and large design space. We have previously developed a PEMFC modeling framework2 which has been a time and cost effective tool for understanding and optimizing the complex multi-physics phenomena within PEMFCs; however, several of the cell parameters used in the modeling have large spread in measured data.3 Furthermore, several transport parameters such as water adsorption kinetics have not been accurately measured and the approximations are spread over several orders of magnitude. These uncertainties cause problems in ascertaining accuracy of the modeling approach and they reduce the predictive power of the numerical models. The aim of this work is to identify the sensitivity of PEFC numerical model outputs to various input parameters. The previously in-house developed MEA modeling framework2 is used for PEMFC modeling. The sensitivity of the model outputs with respect to inputs parameters is obtained by analyzing the condition numbers for different output-input pairs at varying operating conditions. An example of the sensitivity analysis is shown in Figure 1. The condition numbers are obtained for the entire possible range of input parameters at varying operating conditions to identify the most crucial parameters of the PEMFC model. Based on our preliminary analysis, parameters related to kinetics (exchange current density and ECSA) and heat/water management in electrodes and ionomer (ionomer fraction, thermal conductivity) are most crucial. One of the major goals of this work is to identify the most crucial set of parameters towards which the model shows maximum sensitivity. This will guide future experimentalists to measure these properties with higher accuracy. Furthermore, the sensitivity analysis will also enable us to optimize the PEMFC performance by selectively targeting the most sensitive parameters and thereby making the largest impact. Acknowledgements The work is funded under the Fuel Cell Performance and Durability Consortium (FC-PAD), by the Fuel Cell Technologies Office (FCTO), Office of Energy Efficiency and Renewable Energy (EERE), of the U.S. Department of Energy under contract number DE-AC02-05CH11231. The authors would like to thank Nathan Craig at Robert Bosch LLC for his valuable input in designing the sensitivity analysis. The authors would also like to thank Giovanna Bucci and Matthias Hanauer at Robert Bosch for their valuable inputs and discussion. References P. K. Sinha, W. Gu, A. Kongkanand and E. Thompson, J. Electrochem. Soc., 158, B831 (2011). L. M. Pant, M. R. Gerhardt, N. Macauley, R. Mukundan, R. L. Borup and A. Z. Weber, Electrochim. Acta, 326, 134963 (2019). R. Vetter and J. O. Schumacher, ArXiv181110091 Phys. (2018). Figure 1

Pant, Lalit↗

Machine-Learning-Based Rotating Detonation Engine Diagnostics: Evaluation for Application in Experimental Facilities

Real-time monitoring of combustion behavior is a crucial step toward actively controlled rotating detonation engine (RDE) operation in laboratory and industrial environments. Various machine learning methods have been developed to advance diagnostic efficiencies from conventional postprocessing efforts to real-time methods. Here this work evaluates and compares conventional techniques alongside convolutional neural network (CNN) architectures trained in previous studies, including image classification, object detection, and time series classification, according to metrics affecting diagnostic feasibility, external applicability, and performance. Real-time, capable diagnostics are deployed and evaluated using an altered experimental setup. Image-based CNNs are applied to externally provided images to approximate dataset restrictions. Image classification using high-speed chemiluminescence images and time series classification using high-speed flame ionization and pressure measurements achieve classification speeds enabling real-time diagnostic capabilities, averaging laboratory-deployed diagnostic feedback rates of 4–5 Hz. Object detection achieves the most refined resolution of 20 μs in postprocessing. Image and time series classification require the additional correlation of sensor data, extending their time-step resolutions to 80 ms. Comparisons show that no single diagnostic approach outperforms its competitors across all metrics. This finding justifies the need for a machine learning portfolio containing a host of networks to address specific needs throughout the RDE research community.

33 ADVANCED PROPULSION SYSTEMS↗