Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model order reduction”

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 541 records · Page 30

Accelerating high-strain continuum-scale brittle fracture simulations with machine learning

Failure in brittle materials under dynamic loading conditions is a result of the propagation and coalescence of microcracks. Simulating this discrete crack evolution at the continuum level is computationally expensive or, in some cases, intractable, resulting in the need to make broad assumptions or neglect key physics. In this work, we have developed an approach using machine learning that overcomes the current inability to represent meso-scale physics at the macro-scale. Our approach leverages damage and stress data from a computationally expensive high-fidelity model that explicitly resolves microcrack behavior to build an inexpensive machine learning emulator. Once trained, the machine learning emulator is used to predict the evolution of crack length statistics, which then informs a continuum-scale constitutive model. This results in a significant speed-up of the workflow by four orders of magnitude. Both the machine learning emulator and the continuum-scale model are validated against the high-fidelity model and experimental data, respectively, showing excellent agreement. There are two key findings. The first is that we can reduce the dimensionality of the problem, establishing that the machine learning emulator only needs the length of the longest crack and one of the maximum stress components to capture the necessary physics. Another compelling finding is that the emulator can be trained in one experimental setting and transferred successfully to predict behavior in a different setting.

36 MATERIALS SCIENCE↗

Dual Phase Soft Magnetic Laminates for Low-cost, Non/Reduced-Rare-Earth Containing Electrical Machines

To accelerate the mass market adoption of electric drive vehicles, the key technology barriers in electric motors are (1) magnet cost and rare-earth element price volatility; (2) non-rare-earth electric motor performance; and (3) materials property optimization. The goal of this project was to address these barriers by advancing a unique and innovative dual phase soft magnetic material technology and demonstrating the material in a 30-kW synchronous reluctance motor without using any permanent magnet for electric vehicles. Dual phase magnetic materials offer the electric motor designer the ability to locally control the magnetic saturation level in a motor laminate, while at the same time enhancing the mechanical strength of the laminate material, resulting in an enhancement in motor performance and efficiency. Scalable dual phase soft magnetic laminates manufacturing technologies were developed in collaboration with multiple US manufacturers. 1000 lbs of alloy sheet with a thickness of 0.25mm and width of 280 mm was manufactured within the specifications. Batch sizes of up to 240 laminates per run were produced from the alloy sheet. Two prototype motors with dual phase soft magnetic laminates were designed, built, and tested. The major goal of building the subscale prototype as a pathway to develop scalable manufacturing technologies for the dual phase soft magnetic laminates was met. The additional goal of building and testing the subscale prototype in order to validate the calculated performance with the tested motor performance was also met. For the full-scale 30kW continuous power synchronous reluctance motor prototype, the tested performance met the targets in terms of continuous power at the operating speeds up to 8000 rpm. Post-test studies were conducted and the root causes for the discrepancy between the predicted and tested peak power, continuous power at high speed range, and efficiency were identified. Further modeling study showed that the dual phase rotor machine has a 27% higher torque to active weight ratio than an equivalent performance silicon steel rotor machine. Application space and multiple discussions with traction motor and electric vehicle manufacturers for commercialization of the dual phase soft magnetic material technology were identified and conducted. An initial cost model was established based on the developed manufacturing technologies with the US manufacturers. Future paths for further cost reduction were identified, including increasing the market volume by broadening the applications of the dual phase soft magnetic laminate technology for electric machines in other energy sections such as oil & gas, heating, ventilation, and air conditioning (HVAC), and power generation.

33 ADVANCED PROPULSION SYSTEMS↗

Mobility Energy Productivity Evaluation of Prediction-Based Vehicle Powertrain Control Combined with Optimal Traffic Management

Transportation vehicle and network system efficiency can be defined in two ways: 1) reduction of travel times across all the vehicles in the system, and 2) reduction in total energy consumed by all the vehicles in the system. The mechanisms to realize these efficiencies are treated as independent (i.e., vehicle and network domains) and, when combined, they have not been adequately studied to date. This research aims to integrate previously developed and published research on Predictive Optimal Energy Management Strategies (POEMS) and Intelligent Traffic Systems (ITS), to address the need for quantifying improvement in system efficiency resulting from simultaneous vehicle and network optimization. POEMS and ITS are partially independent methods which do not require each other to function but whose individual effectiveness may be affected by the presence of the other. In order to evaluate the system level efficiency improvements, the Mobility Energy Productivity (MEP) metric is used. MEP specifically measures the connectedness of a system while accounting for time and energy externalities of modes that provide mobility in a given location. A SUMO model is developed to reflect real traffic patterns in Fort Collins, Colorado and data is collected by a probe SUMO vehicle which is validated against data collected on a real vehicle driving the same routes through the city. Individual vehicle and system level efficiencies are calculated using SUMO outputs for scenarios which integrate POEMS and ITS independently as well as jointly. Results from application of POEMS and ITS show improvement in energy consumption and travel times respectively when compared to the respective baseline scenarios. Our conclusion is that there are promising synergistic benefits to travel time and energy efficiency when POEMS and ITS are combined.

ADVANCED PROPULSION SYSTEMS↗

Quasilinear modelling of collisional trapped electron modes

The simplification of collision operators is necessary for quasilinear turbulence modelling used with integrated modelling frameworks, such as the gyrokinetic code QuaLiKiz. The treatment of collisions greatly impacts the accuracy of trapped electron mode (TEM) modelling, which is necessary to predict the electron heat flux and the balance between inward and outward particle fluxes. In particular, accurate particle flux predictions are necessary to successfully model density peaking in the tokamak core. We explored two ways of improving collisional TEM model reduction for tokamak core plasmas. First, we carried out linear GENE simulations to study the complex interplay between collisions and trapped electrons. We then used these simulations to define an effective trapped fraction to characterize the collisional TEM based on two key parameters, the local inverse-aspect ratio $\epsilon$ and the collisionality $\nu ^\ast$. One aspect missing from analytical TEM research is that the collisional frequency and the bounce-transit frequency are both velocity dependent; this effective trapped fraction takes both into account. In doing so, we determined that two parameters are not enough to model the collisional TEM, as an additional third free parameter was necessary. We determined that this model, as currently formulated, is not suitable for integrated modelling purposes. Second, we directly improved QuaLiKiz's Krook operator, which relies on two free parameters. We determined that these parameters required adjustments against higher-fidelity collisional models. In order to improve density profile predictions when paired with integrated models, we refined the Krook operator by using GENE simulations as a higher-fidelity point of comparison. We then demonstrate strong improvement of density peaking predictions of QuaLiKiz within the integrated modelling framework JETTO.

fusion plasma↗

Many-Exciton Quantum Dynamics in a Ruddlesden–Popper Tin Iodide

We present a study on the many-body exciton interactions in a Ruddlesden–Popper tin halide, namely, (PEA) 2 SnI 4 (PEA = phenylethylammonium), using coherent two-dimensional electronic spectroscopy. The optical dephasing times of the third-order polarization observed in these systems are determined by exciton many-body interactions and lattice fluctuations. We investigate the excitation-induced dephasing (EID) and observe a significant reduction of the dephasing time with increasing excitation density as compared to its lead counterpart (PEA) 2 PbI 4 , which we have previously reported in a separate publication [J. Chem. Phys. 2020, 153, 164706]. Surprisingly, we find that the EID interaction parameter is four orders of magnitude higher in (PEA) 2 SnI 4 than that in (PEA) 2 PbI 4 . This increase in the EID rate may be due to exciton localization arising from a more statically disordered lattice in the tin derivative. This is supported by the observation of multiple closely spaced exciton states and the broadening of the linewidth with increasing population time (spectral diffusion), which suggests a static disordered structure relative to the highly dynamic lead-halide. Additionally, we find that the exciton nonlinear coherent lineshape shows evidence of a biexcitonic state with low binding energy (<10 meV) not observed in the lead system. We model the lineshapes based on a stochastic scattering theory that accounts for the interaction with a nonstationary population of dark background excitations. Our study provides evidence of differences in the exciton quantum dynamics between tin- and lead-based Ruddlesden–Popper metal halides (RPMHs) and links them to the exciton–exciton interaction strength and the static disorder aspect of the crystalline structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uncertainty Quantification and Sensitivity Analysis of Nuclear Construction Cost Reduction Pathways

High construction costs have plagued recent nuclear projects and they hamper the widespread deployment of nuclear technology. The Nuclear Cost Reduction Tool is a plant economic framework that quantifies the impact that various plant design and construction attributes have on construction costs and cost overruns and shows the positive effects of learning over a series of deployments. However, a downside of the current model is that all model outputs and capabilities are deterministic. To provide a more comprehensive view, this study evaluated the impact of model parameter uncertainty through a sensitivity analysis applied to 18 model parameters. This approach quantified the impact of model uncertainty on the output variables of Net Overnight Capital Cost (Net OCC), Construction Duration (CD), and Levelized Cost of Electricity (LCOE). Monte Carlo analysis revealed uncertainty distributions for these variables, showing that absolute uncertainty decreases over a series of deployments. A local sensitivity analysis showed that even small parameter perturbations (5%) can have a significant impact on project execution, highlighting areas that could reduce costs by billions across an order book of plants. The results of this study have improved the understanding of the model and identified the most impactful model parameters and construction attributes.

Levelized Cost of Electricity↗

tLaSDI: Thermodynamics-informed latent space dynamics identification

Here we propose a latent space dynamics identification method, namely tLaSDI, that embeds the first and second principles of thermodynamics. The latent variables are learned through an autoencoder as a nonlinear dimension reduction model. The latent dynamics are constructed by a neural network-based model that precisely preserves certain structures for the thermodynamic laws through the GENERIC formalism. An abstract error estimate is established, which provides a new loss formulation involving the Jacobian computation of autoencoder. The autoencoder and the latent dynamics are simultaneously trained to minimize the new loss. Computational examples demonstrate the effectiveness of tLaSDI, which exhibits robust generalization ability, even in extrapolation. In addition, an intriguing correlation is empirically observed between a quantity from tLaSDI in the latent space and the behaviors of the full-state solution.

97 MATHEMATICS AND COMPUTING↗

Enhancement of radiation detection performance with reduction of EH 6/7 deep levels in n-type 4H–SiC through thermal oxidation

We report the effect of EH 6/7 electron trap centers alone on the performance of high-resolution radiation detectors fabricated on n-type 4H–SiC epitaxial layers. A Schottky barrier detector (SBD) and a metal-oxide-semiconductor (MOS) capacitor detector fabricated using two sister samples derived from the same 50 μm 4H–SiC parent wafer exhibited widely different energy resolutions of 0.4% and 0.9% for 5486 keV alpha particles. An equivalent noise charge model analysis ruled out the effect of the detector capacitance and the leakage current on the resolution of the detectors. Deep level transient spectroscopic studies revealed the presence of two trapping centers in each detector within the temperature scan range 240–800 K. The Z 1/2 center, a potential electron trap, was detected in both the detectors in equal concentration, which suggested that the observed difference in the energy resolution is due to the presence of the other defect, the EH 6/7 center, in the SBD. The capture cross section of the EH 6/7 center was calculated to be three orders of magnitude higher than the second defect [a carbon antisite vacancy (CAV) center] observed in the MOS detector with an activation energy of 1.10 eV, which accounted for the enhanced electronic trapping in the SBD leading to its poor energy resolution. It has been proposed that the EH 6/7 centers in the SBD have likely been reconfigured to CAV pairs during the thermal growth of the silicon dioxide layer in the MOS detector. The proposed formation mechanism of CAV, a stable qubit state for quantum information processing, addresses the outstanding questions related to the role of defect dynamics in their formation.

61 RADIATION PROTECTION AND DOSIMETRY↗

Performance Assessment of OpenMP Compilers Targeting NVIDIA V100 GPUs

Heterogeneous systems are becoming increasingly prevalent. In order to exploit the rich compute resources of such systems, robust programming models are needed for application developers to seamlessly migrate legacy code from today’s systems to tomorrow’s. Over the past decade and more, directives have been established as one of the promising paths to tackle programmatic challenges on emerging systems. This work focuses on applying and demonstrating OpenMP offloading directives on five proxy applications. We observe that the performance varies widely from one compiler to the other; a crucial aspect of our work is reporting best practices to application developers who use OpenMP offloading compilers. While some issues can be worked around by the developer, there are other issues that must be reported to the compiler vendors. By restructuring OpenMP offloading directives, we gain an 18x speedup for the su3 proxy application on NERSC’s Cori system when using the Clang compiler, and a 15.7x speedup by switching max reductions to add reductions in the laplace mini-app when using the Cray-llvm compiler on Cori.

Davis, Josh↗

A New Track Trigger for Characterization of the Antiproton-Induced Background in the Mu2e Experiment

The Mu2e experiment at Fermilab will enable the search for the neutrinoless muon to electron conversion in the field of an Al nucleus, a charged lepton flavor violating process. If observed, there would be a clear indication of physics beyond the Standard Model. Mu2e aims to reach a single event sensitivity of $3 /times 10^{-17}$, improving from the previous limit by 4 orders of magnitude. This improvement relies on the development of trigger selection systems, designed to discard data from background-induced events by placing kinematic, topological cuts on a particle’s reconstructed track. One of the largest sources of background Mu2e faces is proton-antiproton annihilation. These annihilations produce a 2 GeV shower of particles, among which there could be an electron that mimics the conversion electron signal, with an expected number of 0.010 ± 0.010. The large uncertainty on this number is dominated by the systematic uncertainty associated with the theoretical production model. To better characterize this background, we have developed an antiproton trigger selection by taking advantage of the track multiplicity and topology of these events. We discuss the steps taken in this development and the first performance study of this trigger, evaluating the signal efficiency and background rate. This trigger is essential to enable a data-driven analysis targeting the reduction of the systematic uncertainty of the antiproton-induced background in the Mu2e experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multi-Electron Transfer at H-Terminated p-Si Electrolyte Interfaces: Large Photovoltages under Inversion Conditions

Here, photovoltages for hydrogen terminated p-Si(111) in an acetonitrile electrolyte were quantified with methyl viologen [1,1'-(CH 3 ) 2 -4,4'-bipyridinium](PF 6 ) 2 , abbreviated MV 2+ , and [Ru(bpy) 3 ](PF 6 ) 2 , where bpy is 2,2'-bipyridine, that respectively undergo two and three one-electron transfer reductions. The reduction potentials, E°, of the two MV 2+ reductions occurred at energies within the forbidden bandgap, while the three [Ru(bpy) 3 ] 2+ reductions occurred within the continuum of conduction band states. Bandgap illumination resulted in reduction that were more positive than that measured with a degenerately doped n + -Si demonstrative of a photovoltage, Vph, that increased in the order: MV 2+/+ (260 mV) < MV +/0 (400 mV) < Ru 2+/+ (530 mV) ~ Ru +/0 (540 mV) ~ Ru 0/– (550 mV). Pulsed 532 nm excitation generated electron-hole pairs whose dynamics were nearly constant under depletion conditions and increased markedly as the potential was raised or lowered. A long wavelength absorption feature assigned to conduction band electrons provided additional evidence for the presence of an inversion layer. Collectively, the data reveal that the most optimal photovoltage, as well at the longest electron-hole pair lifetime and the highest surface electron concentration, occur when E° lies energetically within the unfilled conduction band states where an inversion layer is present. The bell-shaped dependence for electron-hole pair recombination with the surface potential was predicted by the time-honored SRH model providing a clear indication that this interface provides access to all four bias conditions, i.e., accumulation, flat band, depletion, and inversion. The implications of these findings for photocatalysis applications and solar energy conversion are discussed.

14 SOLAR ENERGY↗

The strong coupling constant: state of the art and the decade ahead

Theoretical predictions for particle production cross sections and decays at colliders rely heavily on perturbative Quantum Chromodynamics (QCD) calculations, expressed as an expansion in powers of the strong coupling constant α$_{S}$. The current O(1%)uncertainty of the QCD coupling evaluated at the reference Z boson mass, αS(mZ2)=0.1179±0.0009, is one of the limiting factors to more precisely describe multiple processes at current and future colliders. A reduction of this uncertainty is thus a prerequisite to perform precision tests of the Standard Model as well as searches for new physics. This report provides a comprehensive summary of the state-of-the-art, challenges, and prospects in the experimental and theoretical study of the strong coupling. The current αS(mZ2)world average is derived from a combination of seven categories of observables: (i) lattice QCD, (ii) hadronic τ decays, (iii) deep-inelastic scattering and parton distribution functions fits, (iv) electroweak boson decays, hadronic final-states in (v) e$^{+}$e$^{−}$, (vi) e–p, and (vii) p–p collisions, and (viii) quarkonia decays and masses. We review the current status of each of these seven αS(mZ2)extraction methods, discuss novel α$_{S}$ determinations, and examine the averaging method used to obtain the world-average value. Each of the methods discussed provides a ‘wish list’ of experimental and theoretical developments required in order to achieve the goal of a per-mille precision on αS(mZ2)within the next decade.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Robustness of Deep Learning Classification to Adversarial Input on GPUs: Asynchronous Parallel Accumulation Is a Source of Vulnerability

The ability of machine learning (ML) classification models to resist small, targeted input perturbations—known as adversarial attacks—is a key measure of their safety and reliability. We show that floating-point non associativity (FPNA) coupled with asynchronous parallel programming on GPUs is sufficient to result in misclassification, without any perturbation to the input. Additionally, we show that this misclassification is particularly significant for inputs close to the decision boundary and that standard adversarial robustness results may be overestimated up to 4.6 when not considering machine-level details. We first study a linear classifier, before focusing on standard Graph Neural Network (GNN) architectures and datasets used in robustness assessments. We develop a novel black-box attack using Bayesian optimization to discover external workloads that can change the instruction scheduling which bias the output of reductions on GPUs and reliably lead to misclassification. Motivated by these results, we present a new learnable permutation (LP) gradient-based approach to learning floating-point operation orderings that lead to misclassifications. The LP approach provides a worst-case estimate in a computationally efficient manner, avoiding the need to run identical experiments tens of thousands of times over a potentially large set of possible GPU states or architectures. Finally, using instrumentation-based testing, we investigate parallel reduction ordering across different GPU architectures under external background workloads, when utilizing multi-GPU virtualization, and when applying power capping. Our results demonstrate that parallel reduction ordering varies significantly across architectures under the first two conditions, substantially increasing the search space required to fully test the effects of this parallel scheduler-based vulnerability. These results and the methods developed here can help to include machine-level considerations into adversarial robustness assessments, which can make a difference in safety and mission critical applications.

Shanmugavelu, Sanjif [Maxeler Technologies, a Groq↗

Circulating Current Reduction of Back-to-Back MMC with Advanced Grid-Support Functionalities for Medium-Voltage Applications

A Lyapunov energy function based circulating current reduction scheme for a back to back modular multilevel converter based system for a microgrid application has been proposed in this digest. In order to ensure no circulating or zero sequence currents within the modules, the dynamics from one leg of the converter to the other are evaluated. A Lyapunov energy function is then developed and utilizing the circulating current dynamics of the converter the overall control law is implemented. In this digest a decentralized hierarchical control architecture has been implemented for the main current control; where the outer loop or the central controller implements the higher level grid interconnection functions as per IEEE 1547-2018 standards and the inner loop based on synchronous reference frame implements the current control architecture. The inner loops' outputs are the commanded converter voltages either on either side which are augmented by the Lyapunov energy function based circulating current suppression architecture. To evaluate the efficacy of the overall system, the overall system is modeled around a 300kVA back to back SiC multilevel modular converter used as a direct-connected medium voltage (MV) grid intertie based on MATLAB/Simulink domain. Several important case studies are presented in this paper to prove the effectiveness of the proposed architecture. case study results are presented.

circulating current↗

Primordial black hole dark matter: A quantitative parameter sensitivity comparison across formation mechanisms and particle candidates

Primordial black holes (PBHs) in the asteroid-mass window ( 10 17 – 10 22 g ) can account for all of the dark matter without violating any observational constraint, yet are routinely dismissed as fine-tuned. I put that dismissal to the test by applying three complementary sensitivity measures uniformly across a broad landscape: three noninflationary PBH production mechanisms, six classes of inflationary PBH models, and seven particle dark matter benchmarks, all evaluated against the same observable target. Three distinct naturalness universality classes emerge, determined entirely by the analytic structure of the abundance map rather than by the nature of the dark matter candidate. Biased-domain-wall PBHs, in their least model-dependent (free- V b ) form, have the same low sensitivity, Δ = 4.5 , as off-resonance weakly interacting massive particles and freeze-in particles ( Δ = 2 ), a sensitivity that, because it is constant over the entire parameter space of the construction, also coincides trivially with its own Wilson-normalized average within that parameter space (Section Definition and conventions), an equivalence that concerns only the space over which Δ is computed and is not a naturalness statement about the construction as a whole; a further reduction to Δ = 2 is possible only under the additional, independently motivated but not required, assumption that the domain-wall bias is generated by Planck-suppressed operators; early matter-domination PBHs occupy an intermediate tier alongside coannihilating weakly interacting massive particles (WIMPs), unified by a structural identity in which the sensitivity measure equals the logarithm of the ratio of the formation scale to the matter–radiation equality scale; first-order phase transition PBHs, once the more accurate super-exponential collapse probability is used in place of the single-exponential approximation, instead belong to the same highly sensitive tier as resonant WIMP annihilation and single-field inflationary collapse, for a structurally distinct reason; single-field ultraslow-roll inflationary collapse is severely tuned for a distinct reason: a double exponential in which the power spectrum amplitude is itself exponentially sensitive to the inflaton potential coefficients, on top of the exponential collapse sensitivity of the abundance map. My main conclusion is that the claim that PBH dark matter is generically fine-tuned conflates the worst case with a landscape spanning every naturalness tier. The Barbieri-Giudice sensitivity computed here and the Wilson-normalized measure of Iovino and Riotto answer distinct and mutually consistent questions about the same construction, a distinction I clarify within the two-layer decomposition.

Profumo, Stefano [University of California, Santa ↗

Development of a high current density, high temperature superconducting cable for pulsed magnets

Abstract A low AC loss Rare Earth Barium Copper Oxide (REBCO) cable, based on the VIPER cable technology has been developed by Commonwealth Fusion Systems for use in high field, REBCO based tokamaks. The new cable is composed of partitioned and transposed copper ‘petals’ shaped to fit together in a circular pattern with each petal containing a REBCO tape stack and insulated from each other to reduce AC losses. A stainless steel jacket adds mechanical robustness—also serving as a vessel for solder impregnation—while a tube runs through the middle for cooling purposes. Additionally, fiber optic sensors are placed under the tape stacks for quench detection. To qualify this design, a series of experiments were conducted as part of the SPARC tokamak Central Solenoid Model Coil program—to retire the risks associated with full scale, fast ramping, high flux HTS Central Solenoid (CS) and Poloidal Field (PF) coils for tokamak fusion power plants and net energy demonstrators. These risk study and risk reduction experiments include (1) AC loss measurement and model validation in the range of ~5 T/s, (2) an IxB electromagnetic loading of over 850 kN/m at the cable level and up to 300 kN/m at the stack level, (3) a transverse compression resilience of over 350 MPa, (4) manufacturability at tokamak relevant speeds and scales, (5) cable to cable joint performance, (6) fiber optic based quench detection speed, accuracy, and feasibility, and (7) overall winding pack integration and magnet assembly. The result is a cable technology, now referred to as PIT VIPER, with AC losses that measure fifteen times lower (at ~5 T/s) than its predecessor technology; a 2% or lower degradation of critical current (Ic) at high IxB electromagnetic loads; no detectable Ic degradation up to 570 MPa of transverse compression on the cable unit cell; end to end magnet manufacturing, consistently producing Ic values within 7% of the model prediction; cable to cable joint resistances at 20 K on the order of ~15 nΩ; and fast, functional quench detection capabilities that do not involve voltage taps. This cable technology will be tested comprehensively in a Central Solenoid Model Coil to prove its readiness for compact, high field tokamak operation.

Sanabria, Charlie (ORCID:0000000150175309)↗

Solvent selection for biphasic extraction of 5-hydroxymethylfurfural via multiscale modeling and experiments

We introduce a comprehensive conceptual framework for selecting solvents for reactive extraction in biphasic organic-water systems and demonstrate it for the separation of HMF (5-hydroxylmethylfurfural), a platform chemical produced in the acid-catalyzed dehydration of hexoses. We first perform in silico screening of ~2500 solvents, from the ADFCRS-2018 database using the ADF COSMO-RS implementation, and classification, based on the solvent partition coefficient. We then determine experimentally the partition coefficients for HMF, fructose, and products of HMF rehydration (levulinic acid (LA), and formic acid (FA)), the mutual water-organic solvent solubilities, and the separation factors in >50 select solvents spanning multiple homologous series at room temperature and a typical reaction temperature with in situ sampling. We find that COSMO-RS is excellent for screening purposes (typical error in most cases within a factor of ~2). Increased temperatures lead to significant reduction in partitioning, and room temperature measurements are clearly inadequate for solvent selection. Upon down selecting classes of solvents based on separation performance, we perform experimental thermal stability and reaction compatibility studies of a small set of solvents at relevant reactive-extraction temperatures. We discover that many substituted phenols exhibit an order-of-magnitude increase in partitioning compared to conventional solvents due chiefly to hydrogen bond interactions and show the necessary stability but retain a significant fraction of water and LA, factors that need to be considered in technoeconomic analysis. In contrast, anilines, aldehydes, and acids are good to excellent regarding separation but incompatible with this specific reaction medium. Furthermore, this multifaceted framework can be extended to other biomass-derived products and processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cummins R-SOFC System Development

The overall purpose of this project was to reduce the Reversible-Solid Oxide Fuel Cell (R-SOFC) system cost by developing two technologies, an improved cell design and the incorporation of an ejector in the fuel recycle loop instead of a blower. A Simulink model of the baseline SOFC system was developed and calibrated with experimental test data. The R-SOFC system model was built by integrating GT Suite developed models of the steam generation components into the baseline Simulink SOFC system model. The ability to run the stack in SOEC operating mode was also added to the model. The system model was used to explore the ability of the R-SOFC system to meet operational constraints on Steam/Carbon ratio and H2 concentration on the fuel side electrode. A CFD ejector model was developed and used to explore a range of ejector design parameters, leading to the final ejector design that was prototyped for testing. A prototype steam ejector was first tested in a laboratory environment using room temperature air. The steam ejector was subsequently tested using the full hot recycle loop with all relevant heat exchangers and steam generation components. The test conditions utilized temperatures, pressures, and flow rates expected in an R-SOFC application. Throughout the experimental testing work, ejector performance test data was used to improve and then validate the CFD ejector model. A CFD cell model was developed and used to optimize thermal gradients, voltage, and cost of a new cell substrate design. A CFD comparison of co-flow and cross-flow cell designs informed the decision to use a co-flow design for the new cell substrate. Multiple rounds of CFD simulation were used to improve the cell design to minimize the variation in air and fuel distribution across cell channels and to minimize the variation in air and fuel distribution across different cells in the stack. A few prototypes of the new cell substrate design were produced and validated in a laboratory environment by thermally spraying and verifying that they met established manufacturing specifications. The cell manufacturing process was adjusted in order to bring these metrics within acceptable tolerances. Cummins’ internal calculations show that the new cell design reduces cost ~50% compared to the baseline cell, while the ejector + superheater/boiler concept reduces cost of the recycle loop by ~40%. The impact of these cost reductions on the cost of producing H2 will depend on the specific system where they are applied. Therefore, a Techno-Economic Analysis was completed using system cost as a variable, and showing how the NREL Current and Future system costs translate into H2 production cost.

Henrichsen, Lars↗