Empirical and Analytical Eigenfunctions as Bases for Reduced Order Methods in Heat Transfer
Explore the source record for details and available documents.
Bates, Jakob Graham [Sandia National Laboratories ↗
SEARCH · Engineering Papers
Results for “reduced order methods”
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.
Reduced-Order Modeling Method for Phase-Change Thermal Energy Storage Heat Exchangers
Thermal energy storage can facilitate the effective utilization of renewable energy. To speed up the design process of thermal energy storage devices, it is critical to develop fast and accurate modeling methods for phase change material embedded heat exchangers (PCM HXs). This study developed and compared two approximation-assisted reduced-order PCM HX models for the simulation of thermal storage components and systems, which were verified against a validated finite-volume model. They are a pure black-box model, and a grey-box model based on the Number of Transfer Units (or effectiveness-NTU) approach. We used the reduced-order models to predict the performance of standalone PCM HX for 1000 cases with different fluid inlet conditions and HX designs. We also integrated the models with a vapor compression system to predict the compressor energy consumption and total charging time under various conditions. The results show that overall, the black-box model gave more accurate results than the grey-box model. On average, the mean absolute deviation in the PCM HX fluid outlet temperature was 0.05 K and 0.1 K for the black-box model and grey-box model, respectively. The grey-box model gave larger temperature deviation toward the end of the phase change process, due to the simplified two-node PCM representation. On the system level, the mean absolute deviation in compressor energy consumption was 0.2% and 0.3% for the black-box model and grey-box model, respectively. Their mean absolute deviation in total charging time was 1.1% and 2.6%, respectively. In terms of computation efficiency, the system simulation speedup ratio gained by using the reduced-order models was 11 to 57. System simulation time was decreased from an average of 1465 s to an average of 59 s. This shows that the proposed reduced-order modeling methods can be used to predict system performance of PCM HXs with less than 3% accuracy penalty, and 25 times less computational time than finite-volume models, enabling faster design and evaluation of PCM thermal storage devices.
25 ENERGY STORAGE↗
E3SM Atmosphere surrogate construction using machine learning and reduced order modeling methods.
Abstract not provided.
Chowdhary, Kamaljit↗
A reduced-order shifted boundary method for parametrized incompressible Navier–Stokes equations
Here, we investigate a projection-based reduced order model of the steady incompressible Navier–Stokes equations for moderate Reynolds numbers. In particular, we construct an “embedded” reduced basis space, by applying proper orthogonal decomposition to the Shifted Boundary Method, a high-fidelity embedded method recently developed. We focus on the geometrical parametrization through level-set geometries, using a fixed Cartesian background geometry and the associated mesh. This approach avoids both remeshing and the development of a reference domain formulation, as typically done in fitted mesh finite element formulations. Two-dimensional computational examples for one and three parameter dimensions are presented to validate the convergence and the efficacy of the proposed approach.
"embedded" basis↗
Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations
In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.
97 MATHEMATICS AND COMPUTING↗
Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications
This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.
97 MATHEMATICS AND COMPUTING↗
Variance-Reduced Accelerated First-Order Methods: Central Limit Theorems and Confidence Statements
In this paper, we consider a strongly convex stochastic optimization problem and propose three classes of variable sample-size stochastic first-order methods: (i) the standard stochastic gradient descent method, (ii) its accelerated variant, and (iii) the stochastic heavy-ball method. In each scheme, the exact gradients are approximated by averaging across an increasing batch size of sampled gradients. We prove that when the sample size increases at a geometric rate, the generated estimates converge in mean to the optimal solution at an analogous geometric rate for schemes (i)–(iii). Based on this result, we provide central limit statements, whereby it is shown that the rescaled estimation errors converge in distribution to a normal distribution with the associated covariance matrix dependent on the Hessian matrix, the covariance of the gradient noise, and the step length. If the sample size increases at a polynomial rate, we show that the estimation errors decay at a corresponding polynomial rate and establish the associated central limit theorems (CLTs). Under certain conditions, we discuss how both the algorithms and the associated limit theorems may be extended to constrained and nonsmooth regimes. As a result, we provide an avenue to construct confidence regions for the optimal solution based on the established CLTs and test the theoretical findings on a stochastic parameter estimation problem.
Lei, Jinlong↗
Learning physics-based reduced-order models from data using nonlinear manifolds
Here we present a novel method for learning reduced-order models of dynamical systems using nonlinear manifolds. First, we learn the manifold by identifying nonlinear structure in the data through a general representation learning problem. The proposed approach is driven by embeddings of low-order polynomial form. A projection onto the nonlinear manifold reveals the algebraic structure of the reduced-space system that governs the problem of interest. The matrix operators of the reduced-order model are then inferred from the data using operator inference. Numerical experiments on a number of nonlinear problems demonstrate the generalizability of the methodology and the increase in accuracy that can be obtained over reduced-order modeling methods that employ a linear subspace approximation.
97 MATHEMATICS AND COMPUTING↗
DT-HYDRO
The software solves the time dependent, one dimensional (1D) coupled mass and momentum balance equations governing the elastic flow of water through the penstock, turbine and draft tube in a hydroelectric facility using a high order finite volume based method. The numerical method is based on the Kurganov-Tadmor central method paired with the Monotonic Upstream-centered Scheme for Conservation Laws (MUSCL). This solution method accurately resolves the fast transient behavior of the flow, including water hammer. Additionally, the software estimates the full 3D flow field within the turbine chamber in real time, a feat that is made possible by leveraging pre-computed CFD results by utilizing a reduced order modeling method based on an efficient randomized singular value decomposition (SVD) driven proper orthogonal decomposition (POD) with POD-mode weight regression. The reduced order model of the 3D flow is directly coupled to the 1D elastic flow model so the entire flow field through the penstock and turbine system is resolved quickly and with high fidelity.
Gurecky, William [Oak Ridge National Laboratory (O↗
Bayesian operator inference for data-driven reduced-order modeling
This work proposes a Bayesian inference method for the reduced-order modeling of time-dependent systems. Informed by the structure of the governing equations, the task of learning a reduced-order model from data is posed as a Bayesian inverse problem with Gaussian prior and likelihood. The resulting posterior distribution characterizes the operators defining the reduced-order model, hence the predictions subsequently issued by the reduced-order model are endowed with uncertainty. The statistical moments of these predictions are estimated via a Monte Carlo sampling of the posterior distribution. Since the reduced models are fast to solve, this sampling is computationally efficient. Furthermore, the proposed Bayesian framework provides a statistical interpretation of the regularization term that is present in the deterministic operator inference problem, and the empirical Bayes approach of maximum marginal likelihood suggests a selection algorithm for the regularization hyperparameters. The proposed method is demonstrated on two examples: the compressible Euler equations with noise-corrupted observations, and a single-injector combustion process.
97 MATHEMATICS AND COMPUTING↗
Non-autoregressive time-series methods for stable parametric reduced-order models
Advection-dominated dynamical systems, characterized by partial differential equations, are found in applications ranging from weather forecasting to engineering design where accuracy and robustness are crucial. There has been significant interest in the use of techniques borrowed from machine learning to reduce the computational expense and/or improve the accuracy of predictions for these systems. These rely on the identification of a basis that reduces the dimensionality of the problem and the subsequent use of time series and sequential learning methods to forecast the evolution of the reduced state. Often, however, machine-learned predictions after reduced-basis projection are plagued by issues of stability stemming from incomplete capture of multiscale processes as well as due to error growth for long forecast durations. To address these issues, we have developed a non-autoregressive time series approach for predicting linear reduced-basis time histories of forward models. In particular, we demonstrate that non-autoregressive counterparts of sequential learning methods such as long short-term memory (LSTM) considerably improve the stability of machine-learned reduced-order models. Further, we evaluate our approach on the inviscid shallow water equations and show that a non-autoregressive variant of the standard LSTM approach that is bidirectional in the principal component directions obtains the best accuracy for recreating the nonlinear dynamics of partial observations. Moreover-and critical for many applications of these surrogates-inference times are reduced by three orders of magnitude using our approach, compared with both the equation-based Galerkin projection method and the standard LSTM approach.
97 MATHEMATICS AND COMPUTING↗
Control oriented model order reduction for the hydraulic system of an autonomous wheel loader
Control design for an autonomous wheel loader requires a high-fidelity and low-order dynamics model. The objective of order reduction is to decrease the number of states in the model while maintaining performance comparable to the original model. Of all the power components in a wheel loader, the model order reduction for steering and working hydraulic systems is most demanding due to its nonlinear and stiff dynamical nature. This paper describes a physics-inspired model order reduction method that can reduce the model order by close to thirty percent. It is achieved by discarding orders with fast dynamics and consolidating multiple states into less new states. Furthermore, this method can also be extended to other types of off-road vehicles such as excavators, bulldozers, etc. Simulation results demonstrate that the reduced-order model exhibits very similar performance to the full-order model under typical operating conditions of wheel loaders, with an output error of less than six percent.
42 ENGINEERING↗
Offline Maximizing Minimally Invasive Proper Orthogonal Decomposition for Reduced-Order Modeling of S n Radiation Transport
Deterministic solutions to the Sn radiation transport equation can be computationally expensive to calculate. Reduced-order modeling enables efficient approximation of the full-order model (FOM) solution. We propose a novel method for constructing reduced-order models (ROMs) of the S n radiation transport equation, offline maximizing minimally invasive (OMMI) proper orthogonal decomposition (POD). POD uses the method of snapshots to create a reduced-order basis for constructing an ROM. Minimally invasive POD leverages the sweep infrastructure existing in deterministic transport codes to create a POD-based ROM, even when infeasible by traditional methods. Offline maximizing minimally invasive proper orthogonal decomposition (OMMI-POD) extends minimally invasive POD by performing sweeps offline, therefore maximizing the potential speedup. OMMI-POD does so by creating a library of reduced systems from a training set. This library of reduced systems is then interpolated to provide a rapid approximate solution of the S n radiation transport equation. The model is evaluated on a set of test problems, achieving a low error with a 466 times speedup over the FOM. Also presented is a study of the effect of sampling method on the performance of OMMI-POD, specifically comparing naive uniform sampling to the more accurate and computationally expensive greedy sampling.
97 MATHEMATICS AND COMPUTING↗
Non-linear Wire Rope Isolator Model to Enhance Transportation Simulation of Fragile Equipment
Abstract Introduction Wire rope isolators are often used as means of vibration isolation for the transportation of fragile machine equipment. When deployed as suspension systems in a transportation scenario, a simple spring–damper model is often used to predict the dynamic loads that act on the structure being transported. However, this simple model has proved to be too simplistic to be used in the development of an accurate numerical model that can predict the vibration levels experienced by the inner components of fragile equipment during its transportation. This paper describes the experimental tests conducted on a wire rope isolator used for the transportation of the prototype SSR1 cryomodule. Materials and methods A hysteretic Bouc–Wen model has been used to analytically describe the force–deformation relationship of the wire rope isolators. The developed model of the isolator has been implemented in a larger model to simulate the actual transportation of the prototype SSR1 cryomodule, a section of the new PIP-II linear accelerator under construction at Fermilab. A series of multibody dynamic simulations with rigid and flexible components was used to numerically determine the acceleration of some critical components. Results and conclusions An actual experimental transportation was simulated using two numerical models: the developed Bouc–Wen model and a conventional spring–damper model. It is shown how the Bouc–Wen formulation of the isolator characteristics drastically improves the correspondence between experimental and simulated results if compared to a spring–damper model, especially in the range of 0–30 Hz which is the most critical for transportation problems.
42 ENGINEERING↗
Bridging Scales in Bioenergy and Catalysis: A Review of Mesoscale Modeling Applications, Methods, and Future Directions
Between the molecular and reactor scales, which are familiar to the chemical engineering community, lies an intermediate regime, here termed the “mesoscale,” where transport phenomena and reaction kinetics compete on similar time scales. Bioenergy and catalytic processes offer particularly important examples of mesoscale phenomena owing to their multiphase nature and the complex, highly variable porosity characteristic of biomass and many structured catalysts. In this review, we overview applications and methods central to mesoscale modeling as they apply to reaction engineering of biomass conversion and catalytic processing. A brief historical perspective is offered to put recent advances in context. Applications of mesoscale modeling are described, and several specific examples from biomass pyrolysis and catalytic upgrading of bioderived intermediates are highlighted. Methods including reduced order modeling, finite element and finite volume approaches, geometry construction and import, and visualization of simulation results are described; in each category, recent advances, current limitations, and areas for future development are presented. Owing to improved access to high-performance computational resources, advances in algorithm development, and sustained interest in reaction engineering to sustainably meet societal needs, we conclude that a significant upsurge in mesoscale modeling capabilities is on the horizon that will accelerate design, deployment, and optimization of new bioenergy and catalytic technologies.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Vibration of a cryomodule’s coupler during road transportation: Multibody analysis and experimental validation
The road transportation of critical components for cryogenic cryomodules is a challenging task. Indeed, the shipping and the handling of fundamental sub-assemblies can expose these components to harmful dynamic loads that can jeopardize their structural integrity. This article has the purpose of developing a finite element multibody model capable of assessing the acceleration on a test coupler for a cavity of the Single Spoke Resonator 1 subjected to a typical road transportation trip. A multibody model consisting of rigid bodies and finite element–imported flexible bodies has been created allowing a comparison of the simulated and experimental acceleration. Afterwards, an experimental on-road test, whose set-up is described in this article, has been carried out at Fermilab with instrumented components. Being able to effectively simulate and validate the dynamic effects on such sub-assemblies makes it possible, for the future, to simulate more complex structures subject to dynamic loadings, such as the entire cryomodules during shipping. Furthermore, the results from the simulations can be used to guide the design of new suspension systems suitable for the reduction mitigation of structural vibration during the road transportation of cryomodules.
42 ENGINEERING↗
Investigation of finite-volume methods to capture shocks and turbulence spectra in compressible flows
The aim of the present paper is to provide a comparison between several finite-volume methods of different numerical accuracy: the second-order Godunov method with PPM interpolation and the high-order finite-volume WENO method. In this work, the results show that while on a smooth problem the high-order method performs better than the second-order one, when the solution contains a shock all the methods collapse to first-order accuracy. In the context of the decay of compressible homogeneous isotropic turbulence with shocklets, the actual overall order of accuracy of the methods reduces to second-order, despite the use of fifth-order reconstruction schemes at cell interfaces. Most important, results in terms of turbulent spectra are similar regardless of the numerical methods employed, except that the PPM method fails to provide an accurate representation in the high-frequency range of the spectra. It is found that this specific issue comes from the slope-limiting procedure and a novel hybrid PPM/WENO method is developed that has the ability to capture the turbulent spectra with the accuracy of a high-order method, but at the cost of the second-order Godunov method. Overall, it is shown that virtually the same physical solution can be obtained much faster by refining a simulation with the second-order method and carefully chosen numerical procedures, rather than running a coarse high-order simulation. Our results demonstrate the importance of evaluating the accuracy of a numerical method in terms of its actual spectral dissipation and dispersion properties on mixed smooth/shock cases, rather than by the theoretical formal order of convergence rate.
97 MATHEMATICS AND COMPUTING↗