SEARCH · Engineering Papers
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.
Since its' launch on August 10, 1992, the TOPEX/Poseidon satellite hs successfully observed the earth's ocean circulation using a combination of precision orbit determination (POD) and dual-frequency radar altimetry.
The advances in Precise Orbit Determination (POD) over the past three decades have been driven in large measure by the increasing demands of satellite altimetry missions. Since the launch of Seasat in 1978, both tracking-system technologies and orbit modeling capabilities have evolved considerably. The latest in a series of precise (TOPEX-class) altimeter missions is the Ocean Surface Topography Mission (OSTM, also Jason-2). GPS-based orbit solutions for this mission are accurate to 1-cm (radial RMS) within 3-5 hrs of real time. These GPS-based orbit products provide the basis for a near-real time sea-surface height product that supports increasingly diverse applications of operational oceanography and climate forecasting.
A cube satellite (CubeSat) is an evolving and emerging technology that gives a novice or advanced researcher relatively affordable access to space research experiments and applications. The initial CubeSat standard was created in 1999 by California Polytechnic State University, San Luis Obispo and Stanford University’s Space Systems Development Lab to facilitate direct access to space for university students. This initial CubeSat standard has now been adopted by hundreds of organizations worldwide and includes not only universities, educational institutions, but private firms and government organizations. Dozens of CubeSats have been launched since 2003 and have come from more than 29 states in the United States. The CubeSat standard facilitates frequent and affordable access to space with launch opportunities available on most launch vehicles. CubeSats are a class of research spacecraft called nanosatellites and are built to standard CubeSat Units or U dimensions of 10 by 10 by 10 cm and are formally classified as 1U, 2U, 3U, or 6U in size. Most CubeSats are deployed from a Poly-Picosatellite Orbital Deployer called a P-POD. Partnerships among NASA, U.S. industry, and educational institutions are being formed to build upon existing successful CubeSat initiatives with a goal to expand and include launching 50 small satellites from 50 states within the next several years. An extensive and detailed literature review that includes over 830 citations has been conducted to provide a comprehensive resource on both NASA and non-NASA CubeSat experiments and applications that can serve as a guide for background information on CubeSats as well as a valuable resource of lessons learned from CubeSats that have been launched in the past. CubeSats are currently being launched from all over the world on different launch vehicle platforms. Some organizations providing launch opportunities are California Polytechnic State University (http://www.cubesat.org/contactus), ISISPACE Group (https://www.isispace.nl/), Nanoracks (http://nanoracks.com/), Spaceflight Industries, Inc. (http://spaceflight.com/), TriSept Corporation (https://trisept.com/), and Tyvak Nano-Satellite Systems, Inc. (https://www.tyvak.com/).
After a decade of preparation, the GPS precise orbit determinatin (POD) experiment on TOPEX/Poseidon is now yielding definite results.
In this work, we will discuss observations from images acquired from a time-resolved, high-speed self-aligned focusing Schlieren campaign that was performed at the CW-22 facility at NASA Glenn Research Center to understand the dynamic behavior of thick trailing-edge high-pressure turbine blades simulating a ceramic matrix composite (CMC) construction at high inlet turbulence conditions. For the CMC-9 blade, which has a trailing edge thickness of 9% of the axial chord, we identified a regime where an excessive total pressure loss (loss anomaly) is observed only for a narrow range of Reynolds numbers at a fixed exit Mach number of 0.74. The loss anomaly is qualitatively explained by our images, which were taken at 10 distinct blade Reynolds numbers spanning a factor of 6. The images show a significant increase in energy related to the oscillations due to transonic vortex shedding at the Reynolds numbers related to the high loss conditions. This increased energy leads to the formation of strong acoustic waves that turn into shock waves at the highest loss conditions. From our observations stemming from Spectral POD analysis of the high-speed images, we see the acoustic waves produced by the trailing edge vortex shedding exist in all conditions tested; but the shedding frequency has a very slight trend upwards as the Reynolds number is increased, varying about 6% in the range tested. Considering this variation of shedding frequency as a function of Reynolds number, which is well-established for other bluff-body flows, we stipulate there may be a potential feedback mechanism involving an acoustic information transfer path across neighboring blades that may explain why only a narrow range of Reynolds numbers displays strong, shock-forming vortex shedding. We consider a few feedback paths and examine the timing based on the mean flow field from a high-resolution LES simulation. It appears that all feedback mechanisms are viable, presenting an integer number of delay cycles with respect to disturbances generated at the trailing edge. Most noteworthy, however, is the acoustic/shock-boundary layer interaction mechanism at the blade geometric throat. Based on our analysis, this mechanism likely plays an important role in realistic, curved turbine blade passages.
The understanding, modeling and control of the cycle-to-cycle variation (CCV) in the modern internal combustion engine (ICE) is a key scientific challenge to achieve stable engine operation. High CCV in the engine combustion chamber may contribute to partial burn, misfire and knock, which adversely affects the engine performance and may potentially damage the engine. The objective of the current study is to leverage high-fidelity numerical simulations to improve the understanding of the causes of CCV. Using the massively parallel code, Nek5000, multi-cycle, wall-resolved large-eddy simulations (LES) were performed for the General Motors (GM), Transparent Combustion Chamber (TCC-III) optical engine under motored operating conditions. Further, the large-scale structures of the in-cylinder flow were investigated using a triple proper orthogonal decomposition (POD) technique to explore the characteristics of different parts of the flow and their contributions to CCV. The kinetic energy of the subset of flow structures were determined and correlated between the intake and compression strokes. The insights from the analysis of the large-scale flow structures will be used to assist the development of improved engine designs with reduced CCV and enhance the engine performance.
Natural convection in porous media is a highly nonlinear multiphysical problem relevant to many engineering applications (e.g., the process of CO 2 sequestration). Here, we extend and present a non-intrusive reduced order model of natural convection in porous media employing deep convolutional autoencoders for the compression and reconstruction and either radial basis function (RBF) interpolation or artificial neural networks (ANNs) for mapping parameters of partial differential equations (PDEs) on the corresponding nonlinear manifolds. To benchmark our approach, we also describe linear compression and reconstruction processes relying on proper orthogonal decomposition (POD) and ANNs. Further, we present comprehensive comparisons among different models through three benchmark problems. The reduced order models, linear and nonlinear approaches, are much faster than the finite element model, obtaining a maximum speed-up of 7 × 10 6 because our framework is not bound by the Courant–Friedrichs–Lewy condition; hence, it could deliver quantities of interest at any given time contrary to the finite element model. Our model’s accuracy still lies within a relative error of 7% in the worst-case scenario. We illustrate that, in specific settings, the nonlinear approach outperforms its linear counterpart and vice versa. We hypothesize that a visual comparison between principal component analysis (PCA) and t-Distributed Stochastic Neighbor Embedding (t-SNE) could indicate which method will perform better prior to employing any specific compression strategy.
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.
Hypofractionated radiation therapy can be used to treat patients with muscle-invasive bladder cancer unable to have radical therapy. Toxicity is a key concern, but adaptive plan-of the day (POD) image-guided radiation therapy delivery could improve outcomes by minimizing the volume of normal tissue irradiated. The HYBRID trial assessed the multicenter implementation, safety, and efficacy of this strategy.
We present an efficient data-driven regression approach for constructing reduced-order models (ROMs) of reaction-diffusion systems exhibiting pattern formation. The ROMs are learned non-intrusively from available training data of physically accurate numerical simulations. The method can be applied to general nonlinear systems through the use of polynomial model form, while not requiring knowledge of the underlying physical model, governing equations, or numerical solvers. The process of learning ROMs is posed as a low-cost least-squares problem in a reduced-order subspace identified via Proper Orthogonal Decomposition (POD). Numerical experiments on classical pattern-forming systems–including the Schnakenberg and Mimura–Tsujikawa models–demonstrate that higher-order surrogate models significantly improve prediction accuracy while maintaining low computational cost. The proposed method provides a flexible, non-intrusive model reduction framework, well suited for the analysis of complex spatio-temporal pattern formation phenomena.
Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.
This paper compares pressure reconstruction strategies in Model Order Reduction for incompressible flows with convective heat transfer. The Navier-Stokes equation are reduced along with the passive scalar transport equation for the temperature using the POD-Galerkin technique. Six different pressure reconstruction methods are evaluated, two of which are novel to the best of the authors’ knowledge. Accurate pressure reconstruction is key to avoid error buildup in when solving for the conservation of linear momentum at the reduced level. The six approaches are compared using Direct Numerical Simulations of convective heat exchange processes in a 3D Backward Facing Step with a heated cylinder. Additionally, when comparing time-averaged metrics, we observe that the reconstruction methods that approximate the reduced pressure field using techniques borrowed from full order models (mechanical analogy, pressure Poisson, and velocity supremizers) yield higher errors than the methods that seek to stabilize the reduced systems (reduced residual stabilization, artificial divergence, and Uzawa operator).
Fluctuating surface pressures on a bluff body exposed to a boundary layer flow generally are characterized as a spatiotemporally varying random field. In this paper, a dynamic mode decomposition (DMD) was applied to extract dominant features embedded in these random pressure fields. Utilizing an unsupervised machine learning algorithm, spatial modes and their temporal variations were grouped into different clusters at scales, e.g., macro, meso, and micro. A proper orthogonal decomposition (POD) of the experimental data was carried out to observe commonalities and distinctive perspectives each decomposition offers. Here, a comprehensive examination of the DMD/POD for their convergence criteria, data sufficiency, and modal components analysis was conducted. The physical interpretation of the spatiotemporal pressure field based on these decomposition schemes was discussed. At different scales, the DMD modes can capture the evolution of aerodynamic features, e.g., convection of vortices (or vortex tubes) and other structures. The distribution of energy among these three broad scales also reflects an energy cascade in pressure fluctuations akin to turbulence.
Partitioned methods allow one to build a simulation capability for coupled problems by reusing existing single-component codes. In so doing, partitioned methods can shorten code development and validation times for multiphysics and multiscale applications. In this work, we consider a scenario in which one or more of the “codes” being coupled are projection-based reduced order models (ROMs), introduced to lower the computational cost associated with a particular component. We simulate this scenario by considering a model interface problem that is discretized independently on two non-overlapping subdomains. Here we then formulate a partitioned scheme for this problem that allows the coupling between a ROM “code” for one of the subdomains with a finite element model (FEM) or ROM “code” for the other subdomain. The ROM “codes” are constructed by performing proper orthogonal decomposition (POD) on a snapshot ensemble to obtain a low-dimensional reduced order basis, followed by a Galerkin projection onto this basis. The ROM and/or FEM “codes” on each subdomain are then coupled using a Lagrange multiplier representing the interface flux. To partition the resulting monolithic problem, we first eliminate the flux through a dual Schur complement. Application of an explicit time integration scheme to the transformed monolithic problem decouples the subdomain equations, allowing their independent solution for the next time step. We show numerical results that demonstrate the proposed method’s efficacy in achieving both ROM-FEM and ROM-ROM coupling.
A novel plasma actuator, the Linear Counter-flow using a Point Embedded Electrode (LCPEE), is developed for the prevention of dynamic stall for a sinusoidal pitching movement between α = 4° and α = 18°. The LCPEE is implemented on a NACA0012 airfoil and tested at a Reynolds Number of Re c = 2 × 105 and reduced frequency of k = π/16. Prior investigations using a standard linear actuator showed that the exposed electrode introduced perturbations passively which delayed dynamic stall when the actuator was off. For the LCPEE actuator, when turned off, there is least passive delay. When the LCPEE is turned on at St f = 50, the dynamic stall is prevented for the sinusoidal pitching motion of the airfoil. The LCPEE actuator is also tested for the same sinusoidal motion between α = 6° and α = 20°. Four cases are considered for the higher α range of motion: actuator off, actuator on at St f = 50 with a sinusoidal input waveform, actuator on at St f = 50 with a triangular input waveform, and actuator on at Stf = 100 with a sinusoidal input waveform. Specifically for the last case, experiment shows no flow reversal demonstrating the efficacy of LCPEE in controlling the dynamic stall. The effects of LCPEE on the flow energy distribution have also been studied by using proper orthogonal decomposition (POD) method.
Modern wind turbines experience uneven inflow conditions across the rotor, due to the ambient flow’s shear and thermal stratification. Such conditions alter the shape and length of turbine wakes and thus impact the loads and power generation of downstream turbines. To this end, understanding the spatial evolution of the individual wakes under different atmospheric conditions is key to controlling and optimising turbine arrays. With this numerical study we aim to obtain a better understanding of the fundamental physics governing the near-wake dynamics of wind turbines under shear and thermal stability, by examining their tip-vortex breakup mechanisms. Our approach considers scale-resolving simulations of a single turbine wake under a linear shear profile as well as the application of harmonic tip perturbations to trigger flow instabilities. For the subsequent analysis we use the proper orthogonal decomposition (POD) method to extract coherent structures from the flow, and we also calculate mean kinetic energy fluxes to quantify each coherent structure’s contribution to wake recovery. The wake’s helical spiral is found to hinder wake recovery for all studied ambient flow conditions, whereas the mutual inductance instability has positive MKE flux leading to an enhanced wake recovery. Finally, the ambient shear has the largest impact on the local MKE flux with respect to downstream location by changing the shape of the curve and location of extrema, whereas thermal stratification has only a minimal impact on the magnitude of the near-wake local MKE flux distribution.
Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. Here, in this paper, we develop, implement, and test two model-driven emulators based on the (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples (“snapshots”) for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires many fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.
Dynamic load models add significant complexity to bulk power system time-domain simulations. The complexity is due to the large number of ordinary differential equations (ODEs) introduced by the dynamic load components such as induction motors. It is challenging to derive reduced-order models (ROMs) for dynamic loads due to the nonlinear functions in their governing equations. This paper applies the discrete empirical interpolation method enhanced proper orthogonal decomposition (DEIM-POD) to approximate the full dynamic load model with the ROM that minimizes the projection error of the nonlinear functions in dynamic load ODEs onto their dominant modes. This approach only requires evaluation of nonlinear functions at selected observation points. The observation points selected by DEIM also provide information for screening critical load buses where dynamic load model parameters contribute the most to the accuracy of ROM across multiple contingencies. The proposed approach is validated on IEEE 9-bus, WECC 179-bus and 2384-bus Polish systems.