SEARCH · Search NASA
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.
Explore the source record for details and available documents.
While projection-based reduced order models can reduce the dimension of full order solutions, the resulting reduced models may still contain terms that scale with the full order dimension. Hyper-reduction techniques are sampling-based methods that further reduce this computational complexity by approximating such terms with a much smaller dimension. The goal of this work is to introduce the points selection algorithm developed by Shin and Xiu as a hyper-reduction method. The selection algorithm was originally proposed as a stochastic collocation method for uncertainty quantification. Since the algorithm aims at maximizing a quantity $\mathcal{S}$ that measures both the column orthogonality and the determinant, we refer to the algorithm as S-OPT. Numerical examples are provided to demonstrate the performance of S-OPT and to compare its performance with a gappy proper orthogonal decomposition (POD) algorithm. Here, we found that using the S-OPT algorithm is shown to predict the full order solutions with higher accuracy than gappy POD especially when the number of sampling points is small, although we note that S-OPT shows slow asymptotic convergence with respect to the number of samples for some applications, e.g., Lagrangian hydrodynamics.
Experimental dynamic substructuring is a means whereby a mathematical model for a substructure can be obtained experimentally and then coupled to a model for the rest of the assembly to predict the response. Recently, various methods have been proposed that use a transmission simulator to overcome sensitivity to measurement errors and to exercise the interface between the substructures; including the Craig-Bampton, Dual Craig-Bampton, and Craig-Mayes methods. This work compares the advantages and disadvantages of these reduced order modeling strategies for two dynamic substructuring problems. The methods are first used on an analytical beam model to validate the methodologies. Then they are used to obtain an experimental model for structure consisting of a cylinder with several components inside connected to the outside case by foam with uncertain properties. This represents an exceedingly difficult structure to model and so experimental substructuring could be an attractive way to obtain a model of the system.
Reduced-order models (ROMs) are a widely used and powerful approach to reducing the complexity of predictive physics-based numerical simulations for a wide range of applications, including electronics and fluid mechanics such as geologic CO 2 sequestration (GCS). ROMs are critical for optimization, sensitivity analysis, model calibration and uncertainty quantification where full-order models cannot be feasibly executed many times. Traditional approaches generate a single ROM for each simulated response (e.g., CO 2 injection rates, pH changes) based on a set of training simulations. Here, we demonstrate that a single ROM can display excellent overall predictive statistics, but have predictions that dramatically and unacceptably deviate from simulator responses especially when the response variable has a large range (i.e., vary over multiple orders of magnitude). For example, we show that a traditional statistically-high-performing GCS ROM (coefficient of determination R 2 of 0.99) can have average absolute relative errors of over 200%. To address this, we propose a new and novel approach where a set of sub-ROMs are generated to overcome the potential pitfalls in traditional single ROM development. The effectiveness of the proposed approach—the ROMster framework—is demonstrated using a case study of predicted CO 2 injection rates for GCS. We find our approach is a robust and general framework for ROM development, reducing the average “error” from 200% to only 4% in our case study.
Applications of reduced order modeling (ROM) to support analysis of complex reactor behavior using high fidelity simulations have developed rapidly in recent years. Reduction implies any computational approach aiming to reduce the cost of the simulation, especially for situations involving repeated executions such as probabilistic risk assessment and uncertainty quantification applications. This article presents a novel non-intrusive methodology to render reduction for multi-physics models by taking advantage of the combined reduction introduced by each sub-physics in the simulation. Next, a surrogate model is constructed in terms of the reduced dimensions. A key component of the proposed methodology is to upper-bound the errors resulting from the reduction to ensure its reliability for subsequent engineering applications. To implement and demonstrate the proposed ROM algorithm, the INL’s MAMMOTH environment is employed to analyze the level of reduction in the coupled radiation-thermal transport modeling of a 2D quarter fuel pin in a light water reactor spectrum. MAMMOTH couples the neutronics model of Rattlesnake module and the fuel performance model of BISON module. Here, the results show that the reduction obtained with coupled physics is more significant than that with individual sub-physics models.
Biomass pyrolysis is a promising approach for producing economic and environmentally-friendly fuels and bioproducts. Biomass pyrolysis experiments show that feedstock properties have a significant impact on product yields and composition. Scientists are developing detailed chemical reaction mechanisms to capture the relationships between biomass composition and pyrolysis yields. These mechanisms can be computationally intensive. In this study, we investigate the use of a machine learning reduced order model (ROM) for assessing the costs and emissions of a pyrolysis biorefinery. Here, we developed a Kriging-based ROM to predict pyrolysis yields of 314 feedstock samples based on the results of a detailed chemical kinetic pyrolysis mechanism. The ROM is integrated into a chemical process model for calculating mass and energy yields in a commercial-scale (2000 tonne/day) biorefinery. The ROM estimated biofuel yields of 65 to 130 gallons per ton of dry biomass. This results in biofuel minimum fuel-selling prices of $2.62 to $5.43 per gallon and emissions of -13.62 to 145 kg of CO 2 per MJ. The ROM achieved an average mean square error of 1.8e-20 and a mean absolute error of 0.53%. These results suggest that ROMs can facilitate rapid feedstock screening for biorefinery systems.
Several reduced-order modeling strategies have been developed to create low-order models of geometrically nonlinear structures from detailed finite element models, allowing one to compute the dynamic response of the structure at a dramatically reduced cost. But, the parameters of these reduced-order models are estimated by applying a series of static loads to the finite element model, and the quality of the reduced-order model can be highly sensitive to the amplitudes of the static load cases used and to the type/number of modes used in the basis. Our paper proposes to combine reduced-order modeling and numerical continuation to estimate the nonlinear normal modes of geometrically nonlinear finite element models. Not only does this make it possible to compute the nonlinear normal modes far more quickly than existing approaches, but the nonlinear normal modes are also shown to be an excellent metric by which the quality of the reduced-order model can be assessed. Hence, the second contribution of this work is to demonstrate how nonlinear normal modes can be used as a metric by which nonlinear reduced-order models can be compared. Moreover, various reduced-order models with hardening nonlinearities are compared for two different structures to demonstrate these concepts: a clamped–clamped beam model, and a more complicated finite element model of an exhaust panel cover.
Microstructure evolution in the electrodes of solid oxide fuel cell is an important degradation mechanism which reduces active sites for redox reaction and the electric conductivity. Phase field models for microstructure evolution simulation are usually expensive for large scale simulations. In this paper, a reduced-order coarsening model is developed using dynamic discrepancy reduced modeling, which reduces the model order by inserting Gaussian process stochastic functions into the dynamic equations of Ostwald ripening. The reduced order model has been calibrated on a dataset generated by a phase field model that has been well validated to experiments. A validating dataset has also been generated with which the model prediction show good agreement. This model is further applied to predict long term microstructure evolution in different SOFC electrodes. Finally, this work is the first attempt of building a degradation model of SOFC using data science techniques.
Not provided.
Abstract not provided.
Novel 2D materials, e.g., graphene oxide (GO), are attractive building blocks in the design of advanced materials due to their reactive chemistry, which can enhance interfacial interactions while providing good in-plane mechanical properties. Recent studies have hypothesized that the randomly distributed two-phase microstructure of GO, which arises due to its oxidized chemistry, leads to differences in nano- vs meso–scale mechanical responses. However, this effect has not been carefully studied using molecular dynamics due to computational limitations. Herein, a continuum mechanics model, formulated based on density functional based tight binding (DFTB) constitutive results for GO nano-flakes, is establish for capturing the effect of oxidation patterns on the material mechanical properties. GO is idealized as a continuum heterogeneous two-phase material, where the mechanical response of each phase, graphitic and oxidized, is informed from DFTB simulations. A finite element implementation of the model is validated via MD simulations and then used to investigate the existence of GO representative volume elements (RVE). We find that for the studied GO, an RVE behavior arises for monolayer sizes in excess to 40 nm. Moreover, we reveal that the response of monolayers with two main different functional chemistries, epoxide-rich and hydroxyl–rich, present distinct differences in mechanical behavior. In addition, we explored the role of defect density in GO, and validate the applicability of the model to larger length scales by predicting membrane deflection behavior, in close agreement with previous experimental and theoretical observations. As such the work presents a reduced order modeling framework applicable in the study of mechanical properties and deformation mechanisms in 2D multiphase materials.
Traditional linear subspace reduced order models (LS-ROMs) are able to accelerate physical simulations, in which the intrinsic solution space falls into a subspace with a small dimension, i.e., the solution space has a small Kolmogorov n-width. However, for physical phenomena not of this type, e.g., any advection-dominated flow phenomena such as in traffic flow, atmospheric flows, and air flow over vehicles, a low-dimensional linear subspace poorly approximates the solution. To address cases such as these, we have developed a fast and accurate physicsinformed neural network ROM, namely nonlinear manifold ROM (NM-ROM), which can better approximate high-fidelity model solutions with a smaller latent space dimension than the LSROMs. Our method takes advantage of the existing numerical methods that are used to solve the corresponding full order models. The efficiency is achieved by developing a hyper-reduction technique in the context of the NM-ROM. Numerical results show that neural networks can learn a more efficient latent space representation on advection-dominated data from 1D and 2D Burgers’ equations. A speedup of up to 2.6 for 1D Burgers’ and a speedup of 11.7 for 2D Burgers’ equations are achieved with an appropriate treatment of the nonlinear terms through a hyper-reduction technique. Finally, a posteriori error bounds for the NM-ROMs are derived that take account of the hyper-reduced operators.
Traditional linear subspace reduced order models (LS-ROMs) are able to accelerate physical simulations in which the intrinsic solution space falls into a subspace with a small dimension, i.e., the solution space has a small Kolmogorov n-width. However, for physical phenomena not of this type, e.g., any advection-dominated flow phenomena such as in traffic flow, atmospheric flows, and air flow over vehicles, a lowdimensional linear subspace poorly approximates the solution. To address cases such as these, we have developed a fast and accurate physics-informed neural network ROM, namely nonlinear manifold ROM (NM-ROM), which can better approximate high-fidelity model solutions with a smaller latent space dimension than the LS-ROMs. Our software takes advantage of the existing numerical methods that are used to solve the corresponding full order models. The efficiency is achieved by developing a hyper-reduction technique in the context of the NM-ROM. Numerical results show that neural networks can learn a more efficient latent space representation on advection-dominated data from 1D and 2D Burgers' equations. A speedup of up to 2.6 for 1D Burgers' and a speedup of 11.7 for 2D Burgers' equations are achieved with an appropriate treatment of the nonlinear terms through a hyper-reduction technique.
Optimization-based coupling (OBC) is an attractive alternative to traditional Lagrange multiplier approaches in multiple modeling and simulation contexts. However, application of OBC to time-dependent problems has been hindered by the computational cost of finding the stationary points of the associated Lagrangian, which requires primal and adjoint solves. This issue can be mitigated by using OBC in conjunction with computationally efficient reduced order models (ROMs). To demonstrate the potential of this combination, in this paper, we develop an optimization-based ROM-ROM coupling for a transient advection-diffusion transmission problem. We pursue the “optimize-then-reduce” path toward solving the minimization problem at each time step and solve reduced space adjoint system of equations, where the main challenge in this formulation is the generation of adjoint snapshots and reduced bases for the adjoint systems required by the optimizer. One of the main contributions of the paper is a new technique for an efficient adjoint snapshot collection for gradient-based optimizers in the context of optimization-based ROM-ROM couplings. In conclusion, we present numerical studies demonstrating the accuracy of the approach along with comparison between various approaches for selecting a reduced order basis for the adjoint systems, including decay of snapshot energy, average iteration counts, and timings.
The 2016 Parameterized Reduced Order Modeling (PROM) Workshop was held in June, 2016, in Albuquerque, NM. This workshop included 30 researchers who took part in a two day discussion regarding the state of the art for PROMs, complimentary reduced order modeling (ROM) theories, and discussion of the future directions of PROM research. The goals of the workshop were three-fold: to assess the relative accuracy, efficiency, and merits of the different PROM methods; to discuss the state of the art for ROMs and how PROMs can benefit from these advances; and to define the pressing challenges for PROMs and a path for future research collaborations.