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 793 records · Page 44

Ultrafast Electron Transfer in Au-Cyanobacteria Hybrid for Solar to Chemical Production

The rise of inorganic-biohybrid organisms for solar-to-chemical production has spurred mechanistic investigations into the dynamics of the biotic-abiotic interface to drive the development of next-generation hybrid systems. The model system, cyanobacteria-gold nanoparticle hybrids, combines a light harvester with a photosynthetic bacterium to drive the reduction of CO2 to glycerol with improved efficiency and increased glycerol production by 14.6%, in comparison to cyanobacteria only. In this work, we report insights into this unique photochemical behavior and propose a charge-transfer pathway from Au nanoparticle to cyanobacteria. Transient absorption (TA) spectroscopy revealed that photoexcited electron transfer rates are on the order of a few ps to the potential electron acceptor in photosystem II. This work represents a promising platform to utilize a conventional spectroscopic methodology to extract insights from more complex biotic-abiotic hybrid systems.

charge transfer↗

Anolyte Enhances Catalyst Utilization and Ion Transport Inside a CO 2 Electrolyzer Cathode

Electrochemical CO 2 reduction is a promising technology to capture and convert CO 2 to valuable chemicals. High Faradaic efficiencies of CO 2 reduction products are achieved with zero-gap alkaline CO 2 electrolyzers with a supporting electrolyte at the anode (anolyte). Herein, we investigate the effect of anolyte on the electrode properties such as catalyst utilization, ionic accessibility etc. of a CO 2 reduction cathode using electrochemical techniques and cell configurations that avoid the complexities related to co-electrolysis. Using 1M KOH as the anolyte and a Cu gas-diffusion-electrode with low Nafion content as the model CO 2 reduction electrode, we find that electrode capacitance (proxy for electrochemically active surface area) and ionic conductivity inside the cathode increase approximately 4 and 447 times, respectively, in presence of KOH. Liquid anolyte wets the electrode’s pore structure more efficiently than capillary condensation of feed water vapor. The ionomer coverage is very low, and its distribution inside the electrode is highly fragmented. Surface ion conduction mechanisms inside the electrode are orders of magnitude lower than the bulk ion conduction in presence of anolyte. This study shows that when an anolyte (e.g., KOH) is used, catalyst utilization and ionic accessibility inside the electrode increase significantly.

30 DIRECT ENERGY CONVERSION↗

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗

Recent Results from CHAMP Tracking and Accelerometer Data Analysis

The CHAMP mission's unique combination of sensors and orbit configuration will enable unprecedented improvements in modeling and understanding the Earth's static gravity field and its temporal variations. CHAMP is the first of two missions (GRACE to be launched in the early part of 02') that combine a new generation of Global Positioning System (GPS) receivers, a high precision three-axis accelerometer, and star cameras for the precision attitude determination. In order to isolate the gravity signal for science investigations, it is necessary to perform a detailed reduction and analysis of the GPS and Satellite Laser Ranging (SLR) tracking data in conjunction with the accelerometer and attitude data. Precision orbit determination based on the GPS and SLR tracking data will isolate the orbit perturbations, while the accelerometer data will be used to distinguish the non-gravitational forces from those due to the geopotential (static, and time varying). In preparation for the CHAMP and GRACE missions, extensive modifications have been made to NASA/GSFC's GEODYN orbit determination software to enable the simultaneous reduction of spacecraft tracking (e.g. GPS and SLR), three-axis accelerometer and precise attitude data. Several weeks of CHAMP tracking and accelerometer data have been analyzed and the results will be presented. Precision orbit determination analysis based on tracking data alone in addition to results based on the simultaneous reduction of tracking and accelerometer data will be discussed. Results from a calibration of the accelerometer will be presented along with the results from various orbit determination strategies.

Luthcke, S. B.↗

Anthropogenic extremely low volatility organics (ELVOCs) Govern the Growth of Molecular Clusters over the Southern Great Plains during the Springtime

New particle formation (NPF) and growth govern cloud condensation nuclei (CCN) concentrations in many regions. The mechanisms governing the nucleation of molecular clusters vary substantially in different regions of the atmosphere. Additionally, the growth of these clusters from ~2 to 20 nm sizes is often governed by the availability of extremely low volatility organic vapours (ELVOCs). While the pathways to ELVOC formation from the oxidation of biogenic monoterpenes with ozone is better understood, the chemical and mechanistic pathways for ELVOC formation from oxidation of anthropogenic organics are not well understood. We integrate measurements and three-dimensional regional model simulations with the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to understand the processes governing new particle formation and growth and secondary organic aerosol (SOA) formation during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems (HI-SCALE) field campaign at the Southern Great Plains (SGP) observatory in Oklahoma, and contrast it with a site within the Bankhead National Forest (BNF), Alabama in Southeast USA, where 5-year long measurements will begin in 2024. Simulations show that nucleation rates are at least an order of magnitude higher at SGP compared to BNF during the springtime days (April 28 and May 14, 2016), largely due to lower H2SO4 concentrations at BNF, which are needed for nucleation. In addition, the larger CS at BNF (compared to SGP) increase the loss of molecular clusters by coagulation to pre-existing particles. Among the 8 different nucleation mechanisms in WRF-Chem, we find that the amine+H2SO4 nucleation mechanism dominates at the SGP site, while the pure organic ion induced nucleation mechanism dominates over BNF. Through various WRF-Chem sensitivity simulations, we find that anthropogenic ELVOCs are critical for explaining the growth of newly formed particles and the resulting number size distribution observed near the surface at the SGP site during the daytime. In addition, we show that treating organic particles as semisolid, with strong diffusion-limited uptake of organic vapours, brings model predictions into closer agreement with the observed evolution of particle size distribution. Simulations also predict that anthropogenic SOA, formed by the oxidation of aromatic volatile organic compounds (VOCs), is the dominant organic aerosol component at SGP, while biogenic SOA dominates particle composition at the BNF site in Southeast USA on these days.

Shrivastava, ManishKumar B.↗

Flight and Ground Operations in Support of Airframe Noise Reduction Tests

The National Aeronautics and Space Administration (NASA) Acoustic Research Measurements (ARM) project was established to evaluate via flight tests the noise reduction benefits of the Adaptive Compliant Trailing Edge (ACTE) technology along with various main landing gear noise reduction concepts. The ACTE replaces the original Fowler flaps on the NASA SubsoniC Research Aircraft Testbed (SCRAT), thus creating a seamless trailing edge that provides significant noise abatement benefits. The various main landing gear noise reduction concepts are grouped under the LAnding Gear noisE Reduction (LAGER) task and consist of fairings placed on the main landing gear along with two separate treatments applied to the main landing gear wheel well cavities. This paper discusses the tasks necessary to prepare each of these technologies for the ARM flights. The LAGER hardware was taken from model-scale concepts tested in wind tunnels to flight hardware, which had to be cleared as airworthy for the ARM flights. The ACTE flaps were initially intended to be removed from the SCRAT prior to the start of the ARM project. Retaining the ACTE flaps on the aircraft for a longer period of time to support the ARM flights resulted in additional inspections and considerations since the ACTE flaps were flown longer and at certain flight conditions for longer periods of time than initially analyzed. The flight and ground operations required for the ARM tests required extensive coordination among multiple groups and organizations in order to be successful. This paper provides an overview of the hardware development, ground operations, and flight operations which went into acquiring the desired acoustic measurements. In general, the flights were successful and demonstrated the noise reduction benefits of the ACTE flaps, the LAGER gear fairings, and the LAGER gear cavity treatments.

ground operations↗

Online adaptive basis refinement and compression for reduced-order models via vector-space sieving

In many applications, projection-based reduced-order models (ROMs) have demonstrated the ability to provide rapid approximate solutions to high-fidelity full-order models (FOMs). However, there is no a priori assurance that these approximate solutions are accurate; their accuracy depends on the ability of the low-dimensional trial basis to represent the FOM solution. As a result, ROMs can generate inaccurate approximate solutions, e.g., when the FOM solution at the online prediction point is not well represented by training data used to construct the trial basis. To address this fundamental deficiency of standard model-reduction approaches, this work proposes a novel online-adaptive mechanism for efficiently enriching the trial basis in a manner that ensures convergence of the ROM to the FOM, yet does not incur any FOM solves. The mechanism is based on the previously proposed adaptive $h$-refinement method for ROMs, but improves upon this work in two crucial ways. First, the proposed method enables basis refinement with respect to any orthogonal basis (not just the Kronecker basis), thereby generalizing the refinement mechanism and enabling it to be tailored to the physics characterizing the problem at hand. Second, the proposed method provides a fast online algorithm for periodically compressing the enriched basis via an efficient proper orthogonal decomposition (POD) method, which does not incur any operations that scale with the FOM dimension. We report these two features allow the proposed method to serve as (1) a failsafe mechanism for ROMs, as the method enables the ROM to satisfy any prescribed error tolerance online (even in the case of inadequate training), and (2) an efficient online basis-adaptation mechanism, as the combination of basis enrichment and compression enables the basis to adapt online while controlling its dimension.

42 ENGINEERING↗

Comparing three generations of D-Wave quantum annealers for minor embedded combinatorial optimization problems

Abstract Quantum annealing (QA) is a novel type of analog computation that aims to use quantum mechanical fluctuations to search for optimal solutions of Ising problems. QA in the transverse Ising model, implemented on D-Wave quantum processing units, are available as cloud computing resources. In this study we report concise benchmarks across three generations of D-Wave quantum annealers, consisting of four different devices, for the NP-hard discrete combinatorial optimization problems unweighted maximum clique and unweighted maximum cut on random graphs. The Ising, or equivalently quadratic unconstrained binary optimization, formulation of these problems do not require auxiliary variables for order reduction, and their overall structure and weights are not highly variable, which makes these problems simple test cases to understand the sampling capability of current D-Wave quantum annealers. All-to-all minor embeddings of size 52, with relatively uniform chain lengths, are used for a direct comparison across the Chimera, Pegasus, and Zephyr device topologies. A grid-search over annealing times and the minor embedding chain strengths is performed in order to determine the level of reasonable performance for each device and problem type. Experiment metrics that are reported are approximation ratios for non-broken chain samples, chain break proportions, and time-to-solution for the maximum clique problem instances. How fairly the quantum annealers sample optimal maximum cliques, for instances which contain multiple maximum cliques, is quantified using entropy of the measured ground state distributions. The newest generation of quantum annealing hardware, which has a Zephyr hardware connectivity, performed the best overall with respect to approximation ratios and chain break frequencies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sunspot models with Alfven wave emission

Sunspot models have been computed on the assumption that the missing flux is transported by undissipated Alfven waves. In order to estimate the flux of these waves, we propose an extension of Opik's cellular model of convection to include the effects of a vertical magnetic field on horizontal gas flow. Horizontal motions, not vertical motions, are impeded more or less severely depending on the electrical conductivity, and this reduces the convective flux. These motions, however, shake the field lines; this is assumed to be a source of Alfven waves, compensating for the reduction in the convective flux. The free parameter D/H (ratio of cell diameter to cell depth) is adjusted such that the total sum of radiative, convective, and Alfven wave fluxes remains constant at all depths and equals the undisturbed solar flux entering the spot from below. Our treatment is speculative to some extent, but it leads to a unique spot model for a given field strength.

Mullan, D. J.↗

Robust highly durable solid oxide fuel cell cathodes – Improved materials compatibility & self-regulating surface chemistry

Solid oxide fuel cells (SOFCs) are electrochemical conversion devices that directly transform hydrogen or hydrocarbon fuels to electricity, with energy efficiencies as high as 90%, coupled with reduced emissions. Several factors, however, remain to be addressed when considering scale-up of SOFC technology, including the need to overcome decreased performance due to sluggish rates of the oxygen reduction reaction (ORR) at the cathode under reduced temperatures and susceptibility to degradation in performance from surface poisoning e.g. from chromia, while limiting the use of critical raw materials (lanthanides and transition metals) present in high performing mixed ionic electronic conducting electrodes like (La,Sr)CoO 3 (LSC). In this project we explored the key descriptors for determining ORR activity in SOFC electrodes and tried to recover performance degradation by applying them to SOFC electrodes. In order to do this, we first selected a model mixed ionic electronic conducting (MIEC) oxide, Pr-doped CeO 2 (Pr 0.1 Ce 0.9 O 2-δ , PCO), which is a chemically stable fluorite and free of inherent poison sources (e.g. Sr segregation in LSC) that potentially react with external impurities such as Cr-species vaporized from the interconnect. The three approaches originally planned in this project are as follows: 1) evaluation of scavenger exsolution characteristics, 2) study of scavengers gettering efficacy following Cr and Si poisoning and 3) integration of new compositions into porous electrodes. Among them, exceptional progress has been made in 2) and 3), especially understanding the role of surface infiltrants in impacting electrode performance and degradation of PCO materials. We found that the Smith acidity scale for binary oxides serves as a powerful descriptor for tuning and predicting the oxygen exchange kinetics on MIEC PCO surfaces. As a result, with infiltration with binary oxides, ranging from strongly basic (Li 2 O) to strongly acidic (SiO 2 ) onto the surface of porous PCO, it was possible to systematically vary the chemical surface exchange coefficient (k chem ) by 6 orders of magnitude! L i2 O increased k chem by nearly 1,000 times over that of pristine PCO, while SiO 2 decreased k chem by nearly the same factor. Strikingly, although the pre-exponential of k chem scales linearly with the acidity of the infiltrated binary oxide, there is nearly no change in the activation energy. With this insight, we attributed the origin of these dramatic changes in k chem values to the systematic increase and decrease in the surface electron density induced by infiltrated binary oxides. More interestingly, although both Cr 2 O 3 and SiO 2 were determined to be acidic by Smith, suggesting that this feature could likely be the primary reason that these compounds serve to poison the ORR on SOFC cathodes, the effect of poisoning could be subsequently tuned by adding multiple infiltrants and controlling their relative surface acidities. We also systematically examined the effect of serial infiltration of basic and acidic oxides. It turned out that serial infiltration of Li not only recovers approximately 20-fold degraded k chem of PCO by acidic Cr 2 O 3 but its k chem is enhanced even beyond that of the non-infiltrated PCO by more than three orders of magnitude. This was further verified with a screen-printing PCO symmetric cell in terms of the electrode performance (area-specific resistance, ASR) related to approach 3). These observations point to acidity as a key descriptor not only in tuning and predicting the ORR activity of SOFC cathodes that largely determines the overall performance of SOFC, but in mitigating and reactivating poisoned electrode performance. This work provides novel guidelines for making the electrode performance much more active and robust in SOFCs, which can further be applied to all applications requiring oxygen exchange reaction, such as electrolyzers, permeation membranes and gas sensors.

30 DIRECT ENERGY CONVERSION↗

Toward predictive permeabilities: Experimental measurements and multiscale simulation of methanol transport in Nafion

A polymer membrane's permeability to solutes determines its suitability for various applications: a permeability value is essential for predicting performance in diverse contexts. Using aqueous methanol permeation through Nafion as an example, we describe a methodology for determining membrane permeability that accounts for boundary layer effects and the possibility of swelling. For the materials and apparatus used herein, analysis of a permeance measurement and computational fluid dynamics simulations show that the mass transfer boundary layer is on the order of ones to tens of microns. Additionally, the data are used to develop and validate a multiscale model describing solute permeation through a hydrated membrane as a series of physical mechanistic steps: reversible adsorption from solution at the membrane interface, diffusion driven by a concentration gradient within the membrane, and reversible desorption into solution at the opposite membrane interface. The validated model is used to predict methanol transport across a solar-driven CO 2 reduction device and to assess the impact of polymer changes on the measured value. The approach of combining experimental data, computational fluid dynamics, and the mechanistic multiscale model is expected to provide more accurate analysis of membrane permeation data in cases with polymer swelling or unusual device geometries, among others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energy conserving quadrature based dimensionality reduction for nonlinear hydrodynamics problems

Projection based dimensionality reduction is used to approximate the full discretiza tion of a PDE problem in order to reduce the cost of its numerical simulation while staying faithful to the full dynamics. For nonlinear problems, hyperreduction methods are additionally needed to remove the dependence of the nonlinear terms on the full problem’s size. In this project, we develop an energy conserving hyperreduction method and apply it to Eulerian-Lagrangian hydrodynamics problems, demonstrating that the reduced model successfully accelerates the numerical simulations while conserving the problem’s total energy with high accuracy.

97 MATHEMATICS AND COMPUTING↗

Ocean Drifter Project

Four drifting buoys were fabricated from mostly local hardware plus animal locating transmitters. The first two working models were dropped from a Wallops aircraft, to test the mechanical integrity and to test the ability of the whole system to withstand the drop and provide usable data. There are two temperature sensors on board, one to measure sea surface temperature and the other to measure air temperature. An onboard computer controls the transmissions, sending three water temperature readings for each air temperature reading in order to differentiate one from the other after data reduction. Guidelines for the development and test of optical drifting buoys include: air deployable and expendable, reasonable cost for multiple deployment, 1 year lifetime, use of existing technology and a minimal validation parameter set.

Mason, C. D.↗

Perspective on the impact of weightlessness on calcium and bone metabolism

As humans venture into space to colonize the moon and travel to distant planets in the 21st century, they will be confronted with a bone disease that could potentially limit their space exploration activities or put them at risk for fracture when they return to earth. It is now recognized that an unloading of the skeleton, either due to strict bed rest or in zero gravity, leads on average to a 1%-2% reduction in bone mineral density at selected skeletal sites each month. The mechanism by which unloading of the skeleton results in rapid mobilization of calcium stores from the skeleton is not fully understood, but it is thought to be related to down regulation in PTH and 1,25-dihydroxyvitamin D3 production. Bone modeling and mineralization in chick embryos is not affected by microgravity, suggesting that bone cells adapt and ultimately become addicted to gravity in order to maintain a structurally sound skeleton. Strategies need to be developed to decrease microgravity-induced bone resorption by either mimicking gravity's effect on bone metabolism, or enhancing physically or pharmacologically bone formation in order to preserve astronauts' bone health.

short duration↗

The role of topography in geodetic gravity field modelling

Masses associated with the topography, bathymetry, and its isostatic compensation are a dominant source of gravity field variations, especially at shorter wavelengths. On global scales the topographic/isostatic effects are also significant, except for the lowest harmonics. In practice, though, global effects need not be taken into account as such effects are included in the coefficients of the geopotential reference fields. On local scales, the short-wavelength gravity variations due to the topography may, in rugged terrain, be an order of magnitude larger than other effects. In such cases, explicit or implicit terrain reduction procedures are mandatory in order to obtain good prediction results. Such effects may be computed by space-domain integration or by fast Fourier transformation (FFT) methods. Numerical examples are given for areas of the Canadian Rockies. In principle, good knowledge of the topographic densities is required to produce the smoothest residual field. Densities may be determined from sample measurements or by gravimetric means, but both are somewhat troublesome methods in practice. The use of a standard density, e.g., 2.67 g/cu cm, may often yield satisfactory results and may be put within a consistent theoretical framework. The independence of density assumptions is the key point of the classical Molodensky approach to the geodetic boundary value problem. The Molodensky solutions take into account that land gravity field observations are done on a non-level surface. Molodensky's problem may be solved by integral expansions or more effective FFT methods, but the solution should not be intermixed with the use of terrain reductions. The methods are actually complimentary and may both be required in order to obtain the smoothest possible signal, least prone to aliasing and other effects coming from sparse data coverage, typical of rugged topography.

Forsberg, R.↗

NASA Disasters Program: A Partner for Coastal Disaster Response, Resiliency, and Risk Reduction

Summary: The NASA Applied Sciences: Disasters Program has kicked off a pilot project as “Mid-Atlantic Resiliency Demonstration Study, Communities at Intensive Risk” in FY17 and FY18; Efforts focus on a holistic, interdisciplinary approach to integrate a full suite of capabilities from numerical weather prediction, tide and surge modeling at high spatial resolution and urban/neighborhood scales, and remote sensing capabilities for water detection and depth estimation; Future efforts - Partner with Federal Agencies, academia professional societies to complement and augment ongoing work using Earth observations unique to NASA - Translate hazards to risk in order to understand local and regional impacts and minimize vulnerability.

Bell, Jordan↗

Iterative Stability Enforcement in Adaptive Antoulas–Anderson Algorithms for \({\boldsymbol{\mathcal{H}_2}}\) Model Reduction

This paper presents an extension of the Adaptive-Antoulas-Anderson (AAA) algorithm for rational modelling. Specifically, our new stable multi-input multi-output AAA (smiAAA) algorithm builds rational approximations of multi-input signals with a common set of stable poles. A new methodology is presented for iteratively enforcing stability constraints on the poles. We demonstrate the strengths of this approach compared to the stability enforcement in the FastAAA algorithm. Results using the smiAAA algorithm are compared with the commonly used Vector Fitting algorithm and the more recently published RKFIT algorithm. Vector Fitting and RKFIT both require the user to input the number of poles to use in the approximations. If the final approximation is not accurate enough, the user must re-start Vector Fitting or RKFIT with a larger number of poles and/or a new starting location for the poles. In contrast, the smiAAA algorithm is designed to allow the user to simply input the desired accuracy of the approximations, and the necessary number of poles is detected automatically. This permits users to produce approximations of a desired accuracy with no knowledge about the underlying order of the system being approximated, preventing the algorithm from ever needing to be rerun. An additional feature for preventing extraneous poles from being returned by AAA is also discussed. The cause of these extraneous poles is efficiently detected and removed by our presented methodology. In conclusion, the examples presented demonstrate that smiAAA can efficiently produce approximations of similar or better accuracy than Vector Fitting and RKFIT while requiring less input from the user.

97 MATHEMATICS AND COMPUTING↗