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 253 records · Page 14

Diagnosing nuclear power plant pipe wall thinning due to flow accelerated corrosion using a passive, thermal non-destructive evaluation method: Feasibility assessment via numerical experiments

Flow accelerated corrosion (FAC) in nuclear power plant pipes is one of the leading causes of accidents, fatalities, damage and outages. Current FAC identification methods employ expensive sensing technology and are “active” methods, where the response of the piping system to an externally-generated thermal, mechanical or optical excitation must be measured. As a result, these techniques require a disruptive and time-consuming setup. Here we propose a method that utilizes pipe surface temperature measurements to passively monitor for FAC-induced pipe wall thinning without the need for expensive equipment or post-installation setup time. This diagnostic method utilizes a simulation data-driven diagnostic model to estimate the amount of thickness reduction in a pipe based on changes in measured steady-state pipe temperatures. In order to reduce the computational burden of generating large, simulation-based datasets, the behavior of the insulation of the pipe was modeled using a suitably calibrated heat transfer boundary parameter. Additionally, global sensitivity analysis was performed to determine system parameter(s), such as the temperature of water flowing inside the pipe, which significantly affect the steady state pipe wall temperature and could cause errors in diagnosis. Two diagnostic models, one using only the change in steady-state temperature as an indicator for FAC-induced pipe wall thinning and the other using water temperature as an additional diagnostic model input were evaluated for their ability to estimate thickness reductions in a pipe using simulated pipe wall temperature data. For the numerical experiments conducted in this work, both models estimated wall thickness with errors within 0.5 mm, indicating that the proposed technique can potentially be used as a low-cost, first-pass method for FAC monitoring.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Constrained Block Nonlinear Neural Dynamical Models

Neural network modules conditioned by known priors can be effectively trained and combined to represent systems with nonlinear dynamics. This work explores a novel formulation for data-efficient learning of deep control-oriented nonlinear dynamical models by embedding local model structure and constraints. The proposed method consists of neural network blocks that represent input, state, and output dynamics with constraints placed on the network weights and system variables. For handling partially observable dynamical systems, we utilize a state observer neural network to estimate the states of the system's latent dynamics. We evaluate the performance of the proposed architecture and training methods on system identification tasks for three nonlinear systems: a continuous stirred tank reactor, a two tank interacting system, and an aerodynamics body. Models optimized with a few thousand system state observations accurately represent system dynamics in open loop simulation over thousands of time steps from a single set of initial conditions. Experimental results demonstrate an order of magnitude reduction in open-loop simulation mean squared error for our constrained, block-structured neural models when compared to traditional unstructured and unconstrained neural network models.

Skomski, Elliott↗

Uncertainty propagation and sensitivity analysis for constrained optimization of nuclear waste vitrification

Abstract The vitrification of high‐level waste (HLW) by heating a mixture of glass‐forming chemicals (GFCs) with the waste can be improved using a constrained optimization problem. This study explores how different uncertainty propagation (UP) methods implemented with the optimization process can affect the glass formulation of nuclear waste glasses. UP is the effort of propagating uncertain inputs through a system to understand and quantify output distributions. Uncertainty intervals are crafted from output distributions to inform the optimization algorithm. UP is often implemented with Monte Carlo (MC) sampling for large nonlinear systems, which can be difficult to implement within a constrained optimization algorithm that requires derivative information. Other UP methods often used for optimization under uncertainty (OUU) can be designed to work within an established constrained optimization framework. Methods of UP are evaluated in this study including iterative sampling approaches, first‐order approximations, and surrogate modeling with machine learning (ML). A method of dimensional reduction based on global sensitivity analysis is introduced to support the UP methods for the large dimensionality of the problem. Analytical UP methods able to achieve similar optimums 10 times faster than the baseline MC approach, and produce 93.9% similar output distributions are reported.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Advanced System Thermal Fluids Solver Development for SAM

This work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key is the implementation of a high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. Leveraging the existing capabilities of the SAM code, significant code coverages were established in the finite volume method code. This in turn allows for a suite of test problems with different problem sizes and levels of complexity to be used to quantify the performance improvement of the finite volume method code. As evidently shown in this study, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE for the wide range of selected problems. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. In addition, for a complex reactor model, transient simulation was performed using the finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development. In this work, short-term priority development and testing items were identified, and long-term code adoption and integration plans were made for the eventual deployment of the finite volume method in the SAM code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interior noise control prediction study for high-speed propeller-driven aircraft

An analytical model was developed to predict the noise levels inside propeller-driven aircraft during cruise at M = 0.8. The model was applied to three study aircraft with fuselages of different size (wide body, narrow body and small diameter) in order to determine the noise reductions required to achieve the goal of an A-weighted sound level which does not exceed 80 dB. The model was then used to determine noise control methods which could achieve the required noise reductions. Two classes of noise control treatments were investigated: add-on treatments which can be added to existing structures, and advanced concepts which would require changes to the fuselage primary structure. Only one treatment, a double wall with limp panel, provided the required noise reductions. Weight penalties associated with the treatment were estimated for the three study aircraft.

Rennison, D. C.↗

State-space formulation of multi-shaker modal analysis

As modal testing techniques have improved in recent years, the demands placed on the accuracy of modal test results have increased. This paper describes a new multi-input, multi-output modal parameter estimation algorithm suitable for the identification of reduced-order models of structures. The method described is a frequency-domain method based on excitation and response spectra. The principal features of the algorithm are: (1) applicability to general linear, time-invariant systems, (2) direct use of multiple input and multiple response data, (3) identification of reduced-order models with minimum user interaction through the use of two automatic model-reduction techniques, (4) identification of a consistent set of modal parameters through the use of numerically stable least-squares techniques, and (5) complete user documentation.

Craig, Roy R., Jr.↗

Computational fluid dynamic control

A general technique is presented for modeling fluid, or gas, dynamic systems specifically for the development of control systems. The numerical methods which are generally used in computational fluid dynamics are borrowed to create either continuous-time or discrete-time models of the particular fluid system. The resulting equations can be either left in a nonlinear form, or easily linearized about an operating point. As there are typically very many states in these systems, the usual linear model reduction methods can be used on them to allow a low-order controller to be designed. A simple example is given which typifies many internal flow control problems. The resulting control is termed computational fluid dynamic control.

Hartley, Tom T.↗

On the Numerical Formulation of Parametric Linear Fractional Transformation (LFT) Uncertainty Models for Multivariate Matrix Polynomial Problems

Robust control system analysis and design is based on an uncertainty description, called a linear fractional transformation (LFT), which separates the uncertain (or varying) part of the system from the nominal system. These models are also useful in the design of gain-scheduled control systems based on Linear Parameter Varying (LPV) methods. Low-order LFT models are difficult to form for problems involving nonlinear parameter variations. This paper presents a numerical computational method for constructing and LFT model for a given LPV model. The method is developed for multivariate polynomial problems, and uses simple matrix computations to obtain an exact low-order LFT representation of the given LPV system without the use of model reduction. Although the method is developed for multivariate polynomial problems, multivariate rational problems can also be solved using this method by reformulating the rational problem into a polynomial form.

Belcastro, Christine M.↗

Wing Shape Sensing from Measured Strain

A new two step theory is investigated for predicting the deflection and slope of an entire structure using strain measurements at discrete locations. In the first step, a measured strain is fitted using a piecewise least squares curve fitting method together with the cubic spline technique. These fitted strains are integrated twice to obtain deflection data along the fibers. In the second step, computed deflection along the fibers are combined with a finite element model of the structure in order to extrapolate the deflection and slope of the entire structure through the use of System Equivalent Reduction and Expansion Process. The theory is first validated on a computational model, a cantilevered rectangular wing. It is then applied to test data from a cantilevered swept wing model.

system equivalent reduction and expansion process↗

Application of Quasi-Static Modal Analysis to an Orion Multi-Purpose Crew Vehicle Test

Bolted structural joints often exhibit load-dependent stiffness and energy dissipation that leads to nonlinear, amplitude dependent frequency and damping in the structure. As an alternative to direct integration of the nonlinear equations of motion, quasi-static modal analysis (QSMA) determines the dependence of frequency and damping on response amplitude using loading behavior from nonlinear static analyses. QSMA has previously been demonstrated to substantially reduce computational cost and maintain accuracy relative to full nonlinear dynamic simulation. This work explores the applicability of QSMA to a complex, large-scale aerospace structure. QSMA is employed to analyze a nonlinear model of test hardware developed to support the Orion Multi-Purpose Crew Vehicle program, which exhibited nonlinear behavior during dynamic testing at flight-like load levels. In addition to the extraction of amplitude-dependent frequency and damping curves, a Bouc-Wen hysteresis model was used in conjunction with the quasi-static results to develop nonlinear, uncoupled, time-domain modal equations of motion for the structure. Excellent agreement was observed between the reduced and full-order nonlinear models, encouraging future employment of QSMA to support accurate and efficient model reduction of structures with bolted joint nonlinearities.

QSMA↗

Enabling Mission Flexibility to Battery Driven Deep Space Endeavors With Generalized Battery-Health-Monitoring Using Physics-Based and Data-Driven Reduced-Order Models

The needs and requirements for an electrochemical energy storage for deep space exploration is well explored. It is often understood that different mission sites and environmental conditions require different battery chemistries or technologies. Additionally, various engineering solutions are deployed to overcome specific chemical challenges. One often overlooked need is the “health” monitoring of an electrochemical storage system. The term generalized health monitoring, as envisioned in this work, refers to the monitoring of various aspects such as electrode health, electrolyte health, reaction pathway health, cooling system health, sensor health, and BMS health [1]. Generalized health monitoring allows mission leads, engineers, and scientists to incorporate flexibility in mission designs, make on-the-fly mission changes, and extend the duration of science missions. Moreover, it enables automation and data-driven decision-making without compromising safety and performance. Recently, our group developed a hierarchy of thermal reduced-order models (TROM) by combining a physics-based modeling approach and data-driven model reduction techniques applied to flight data [2]. The resulting TROMs were found to be not only accurate but also identifiable from the flight data. Consequently, the coefficient of variance of the model parameters is small over the course of hundreds of flights, allowing for monitoring the parameter evolution trajectories as the battery ages and degrades. These parameters constitute the metrics of the generalized health of a battery. Monitoring their evolution allows such models to be used for anomaly detection and prognostics, improving early detection of abnormal behavior and thus enabling timely maintenance, longer battery life, and enhanced battery safety. For this presentation, the practicality of the thermal model will be validated on a pack of 14cells under various topology configurations such as 1S14P, 2P7S, 7S2P, and 1P14S. It is well known that manufacturing and non-uniform aging lead to variability in the performance of a cell, which is exacerbated by cell balancing during active load. Additionally, in extreme scenarios, the paramount objective is to complete the mission, regardless of the stresses on the battery. Topology-induced balancing issues further stress the battery. The goal of this study is to determine if the noise (identifiability) in the reduced-order thermal model parameters is sensitive to topology, cell spacing, cooling strategy, and manufacturing or age variability. The variability in cells is considered by assuming a multimodal distribution for microscopic parameters of a cell (such as porosity, tortuosity, reaction kinetics, volumetric thermal conductivity, and volumetric heat capacity). The compounded effect of manufacturing variability, topological selection, cooling strategies, and cell balancing ages each cell in a battery differently. The study aims to clarify whether the challenge in extracting maximum information depends on the minimum number of sensors or models used for data extraction.

Automation↗

A microfluidic system that replicates pharmacokinetic (PK) profiles in vitro improves prediction of in vivo efficacy in preclinical models

Test compounds used on in vitro model systems are conventionally delivered to cell culture wells as fixed concentration bolus doses; however, this poorly replicates the pharmacokinetic (PK) concentration changes seen in vivo and reduces the predictive value of the data. Herein, proof-of-concept experiments were performed using a novel microfluidic device, the Microformulator, which allows in vivo like PK profiles to be applied to cells cultured in microtiter plates and facilitates the investigation of the impact of PK on biological responses. We demonstrate the utility of the device in its ability to reproduce in vivo PK profiles of different oncology compounds over multiweek experiments, both as monotherapy and drug combinations, comparing the effects on tumour cell efficacy in vitro with efficacy seen in in vivo xenograft models. In the first example, an ERK1/2 inhibitor was tested using fixed bolus dosing and Microformulator-replicated PK profiles, in 2 cell lines with different in vivo sensitivities. The Microformulator-replicated PK profiles were able to discriminate between cell line sensitivities, unlike the conventional fixed bolus dosing. In a second study, murine in vivo PK profiles of multiple Poly(ADP-Ribose) Polymerase 1/2 (PARP) and DNA-dependent protein kinase (DNA-PK) inhibitor combinations were replicated in a FaDu cell line resulting in a reduction in cell growth in vitro with similar rank ordering to the in vivo xenograft model. Additional PK/efficacy insight into theoretical changes to drug exposure profiles was gained by using the Microformulator to expose FaDu cells to the DNA-PK inhibitor for different target coverage levels and periods of time. We demonstrate that the Microformulator enables incorporating PK exposures into cellular assays to improve in vitro–in vivo translation understanding for early therapeutic insight.

Singh, Dharaminder (ORCID:0000000346306403)↗

An experimental investigation of internal area ruling for transonic and supersonic channel flow

A simulated transonic rotor channel model was examined experimentally to verify the flow physics of internal area ruling. Pressure measurements were performed in the high speed wind tunnel at transonic speeds with Mach 1.5 and Mach 2 nozzle blocks to get an indication of the approximate shock losses. The results showed a reduction in losses due to internal area ruling with the Mach 1.5 nozzle blocks. The reduction in total loss coefficient was of the order of 17 percent for a high blockage model and 7 percent for a cut-down model.

Roberts, W. B.↗

Effect of bacterial cell addition on Fe(III) reduction and soil organic matter transformation in a farmland soil

The coupled biogeochemical processes of Fe(III) reduction and organic matter transformation profoundly impact terrestrial carbon cycling. However, little is known about the microbial role in soil organic matter (SOM) transformation during Fe(III) bio-reduction. Here we investigated the bio-reduction behavior of a black farmland soil and corresponding SOM transformation under circumneutral and anoxic conditions. A model dissimilatory Fe reducing bacterium, Geobacter sulfurreducens, was added in either live or dead form in order to enhance Fe(III) reduction in soil and to evaluate the accompanying transformation of SOM. The progress of Fe reduction was monitored and SOM transformation was characterized by various spectroscopy methods. Results showed that addition of either dead or live G. sulfurreducens cells increased the Fe(III) reduction rate and extent. Without cell addition, SOM transformation was insignificant within 13 days of incubation, only with some consumption of aliphatic/protein compounds, apparently due to their higher bio-degradability. Addition of dead or live cells resulted in more drastic SOM transformation, but through different mechanisms. With dead cell amendment, cell necromass and debris stimulated the activity of indigenous soil microbes by serving as extra carbon/energy sources. During Fe(III) reduction, the aliphatic/protein compounds were preferentially consumed by indigenous microbial communities, similar to the treatment without cell addition but with a greater extent of consumption. In comparison, addition of live G. sulfurreducens cells stimulated degradation of less bioavailable compounds, including more saturated and higher molecular weight SOM. It is possible that fast depletion of labile SOM by live G. sulfurreducens cells favored utilization of less bioavailable molecules (such as those with more aromatic structures) by the originally dormant species in native microbial community, suggesting an active role of live Fe(III)-reducing bacteria in affecting SOM transformation. In addition, fast assimilation of microbial related carbon, such as aromatic proteins and microbial byproducts, into SOM pools was also observed. Overall, our results suggest that addition of Fe(III)-reducing bacteria to a native soil can enhance Fe(III) reduction and accelerate the turnover of SOM through various mechanisms. The study provides new insights into coupled Fe(III) reduction and SOM transformation, as well as the “priming effect” of SOM using microbial cells as substrate.

54 ENVIRONMENTAL SCIENCES↗

Model based estimation of image depth and displacement

Passive depth and displacement map determinations have become an important part of computer vision processing. Applications that make use of this type of information include autonomous navigation, robotic assembly, image sequence compression, structure identification, and 3-D motion estimation. With the reliance of such systems on visual image characteristics, a need to overcome image degradations, such as random image-capture noise, motion, and quantization effects, is clearly necessary. Many depth and displacement estimation algorithms also introduce additional distortions due to the gradient operations performed on the noisy intensity images. These degradations can limit the accuracy and reliability of the displacement or depth information extracted from such sequences. Recognizing the previously stated conditions, a new method to model and estimate a restored depth or displacement field is presented. Once a model has been established, the field can be filtered using currently established multidimensional algorithms. In particular, the reduced order model Kalman filter (ROMKF), which has been shown to be an effective tool in the reduction of image intensity distortions, was applied to the computed displacement fields. Results of the application of this model show significant improvements on the restored field. Previous attempts at restoring the depth or displacement fields assumed homogeneous characteristics which resulted in the smoothing of discontinuities. In these situations, edges were lost. An adaptive model parameter selection method is provided that maintains sharp edge boundaries in the restored field. This has been successfully applied to images representative of robotic scenarios. In order to accommodate image sequences, the standard 2-D ROMKF model is extended into 3-D by the incorporation of a deterministic component based on previously restored fields. The inclusion of past depth and displacement fields allows a means of incorporating the temporal information into the restoration process. A summary on the conditions that indicate which type of filtering should be applied to a field is provided.

Damour, Kevin T.↗

Upper Arctic Ocean Properties and Relationships with Sea Ice in CMIP6 Historical Simulations

While current-generation CMIP and OMIP models have clear biases in their upper Arctic Ocean hydrography, it is less clear how these biases impact the models' ability to simulate the observed Arctic sea ice mean state and trends. In this study we seek to quantify cross-relationship between sea ice and ocean states in CMIP6 historical simulations and identify common model behaviors. Multi-model mean (MMM) simulations exhibit accelerated changes in the ice and ocean system since the late 20th century. Underlying the MMM is strong inter-model variation in the simulated ice and ocean mean states and their temporal variability including trends. Despite such inter-model differences, all models show a similar ratio between sea ice reduction and upper ocean warming such that models with higher ocean warming also have higher SIE reduction and vice versa. Our results also highlight the urgent needs of reliable Arctic Ocean observations or data products in order to better contextualize modeling results.

Upper Arctic Ocean↗

Hybrid learning techniques for scientific data reduction with performance guarantees

The research initiatives supported by the U.S. Department of Energy (DOE) Grant DE-SC0022265 are fundamentally aimed at pioneering advanced machine learning (ML) techniques for scientific data compression within high-performance computing (HPC) environments. This comprehensive body of work addresses the critical challenge posed by the exponential growth of data generated by scientific simulations in domains such as fusion energy, climate modeling, and computational fluid dynamics (CFD). A core objective is to develop compression algorithms that achieve substantial data reduction—often by orders of magnitude—while rigorously ensuring the fidelity of both the primary data (PD) and scientifically crucial derived quantities of interest (QoI). The methodologies deployed under this grant integrate sophisticated deep learning architectures, prominently featuring autoencoders, advanced generative models like conditional diffusion, and hybrid learning techniques. Key innovations include the development of Guaranteed Autoencoders (GAE) and the Guaranteed Conditional Diffusion with Tensor Correction (GCDTC) framework, which provide explicit, instance-level error bounds on reconstructed data. Furthermore, specialized strategies such as nonlinear constraint satisfaction are employed to preserve the integrity of QoI, a vital requirement for the trustworthiness of downstream scientific analyses. This research also focuses on the design and implementation of scalable, GPU-accelerated software pipelines that seamlessly integrate into existing HPC workflows, ensuring both computational efficiency and practical applicability. The CAESAR framework, for example, unifies foundation and generative models to create an adaptive and efficient compression solution for spatio-temporal scientific data. Collectively, these efforts represent a significant advancement in mitigating the scientific data deluge, enabling more effective data management, accelerated scientific discovery, and optimized utilization of HPC resources.

97 MATHEMATICS AND COMPUTING↗

Final report- UFL - RAPIDS2: A SciDAC Institute for Computer Science, Data, and Artificial Intelligence

The research initiatives supported by the U.S. Department of Energy (DOE) Grant DE-SC0022265 are fundamentally aimed at pioneering advanced machine learning (ML) techniques for scientific data compression within high-performance computing (HPC) environments. This comprehensive body of work addresses the critical challenge posed by the exponential growth of data generated by scientific simulations in domains such as fusion energy, climate modeling, and computational fluid dynamics (CFD). A core objective is to develop compression algorithms that achieve substantial data reduction—often by orders of magnitude—while rigorously ensuring the fidelity of both the primary data (PD) and scientifically crucial derived quantities of interest (QoI). The methodologies deployed under this grant integrate sophisticated deep learning architectures, prominently featuring autoencoders, advanced generative models like conditional diffusion, and hybrid learning techniques. Key innovations include the development of Guaranteed Autoencoders (GAE) and the Guaranteed Conditional Diffusion with Tensor Correction (GCDTC) framework, which provide explicit, instance-level error bounds on reconstructed data. Furthermore, specialized strategies such as nonlinear constraint satisfaction are employed to preserve the integrity of QoI, a vital requirement for the trustworthiness of downstream scientific analyses. This research also focuses on the design and implementation of scalable, GPU-accelerated software pipelines that seamlessly integrate into existing HPC workflows, ensuring both computational efficiency and practical applicability. The CAESAR framework, for example, unifies foundation and generative models to create an adaptive and efficient compression solution for spatio-temporal scientific data. Collectively, these efforts represent a significant advancement in mitigating the scientific data deluge, enabling more effective data management, accelerated scientific discovery, and optimized utilization of HPC resources.

97 MATHEMATICS AND COMPUTING↗