Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reduced Order Modeling”

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 1,351 records · Page 75

Improved single-particle phase-space distributions for viscous fluid dynamic models of relativistic heavy ion collisions

Comparison of hydrodynamic calculations with experimental data inevitably requires a model for converting the fluid to particles. For shear viscous fluids, several different conversion models have been proposed, such as the quadratic Grad corrections, the Strickland-Romatschke (SR) ansatz, self-consistent shear corrections from linearized kinetic theory, or corrections from the relaxation time approach. We show, via comparison to nonlinear 2 → 2 kinetic theory evolution, that even for low specific shear viscosities 0.03 ≲ $\textit{η/s}$ ≲ 0.2 , these four models have varying accuracy and still significant errors in reproducing particle phase space densities from hydrodynamic fields. Moreover, we demonstrate that the overall reconstruction error of additive shear viscous $f = f_{eq} + δf$ models dramatically reduces by up to one order of magnitude, if one ensures through exponentiation that $\textit{f}$ is always positive. Further, we also illustrate how even more accurate viscous $\textit{δf}$ models can be constructed by incorporating the first time derivatives of hydrodynamic fields to account for the past evolution of the system. Such time derivatives are readily available in hydrodynamic simulations, though usually not included in the output. Though our comparisons are limited to a one-component massless system undergoing a (0 + 1)-dimensional (0 + 1D) longitudinal boost-invariant expansion, we expect that the improved models will be useful in 2 + 1D and 3 + 1D hydrodynamic studies as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Selecting Charged Current Muon Neutrino Interactions on Argon with the Short-Baseline Near Detector

The uncertainties associated with neutrino-nucleus cross section modelling need to be reduced in order to measure neutrino oscillation parameters, the mass hierarchy and CP violation to a high degree of precision. SBND, a 112 ton liquid argon time projection chamber, will be the near detector of the short-baseline neutrino program at Fermilab. Once data taking begins in 2022, it will provide flux constraints for sterile neutrino searches and produce world leading neutrino-argon cross sections with up to seven million neutrino events in three years. This thesis describes the simulation, calibration and reconstruction tools developed to measure charged current muon neutrino interactions.Backgrounds from cosmic ray muons, neutrino interactions outside of the detector and different neutrino interaction channels were identified and removed as part of this process.The signal events were selected with purity of 90% and an efficiency of 67%, providing a strong starting point for meas urements of exclusive channels. The expected statistical and systematic uncertainties on rate predictions and cross section measurements at the full 6.6×10$^{20}$ proton on target exposure were evaluated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Two time scale stabilization of systems with output feedback

The problem of constant gain output feedback regulator design for linear systems with ill-conditioned dynamics is considered in the context of singular perturbation theory. A design approach is developed in which gains can be separately calculated to stabilize reduced-order slow and fast subsystem models. By employing the notion of combined control and observation spillover suppression, conditions are derived assuring that these gains will stabilize the full-order system, assuming sufficient frequency separation between the slow and fast subsystems. An LQ design procedure is described in which the spillover suppression conditions are satisfied by adjoining penalty functions to the subsystem performance indices. The theory is demonstrated in a controller design for a flexible space structure.

Moerder, D. D.↗

Robust Assignment Of Eigensystems For Flexible Structures

Improved method for placement of eigenvalues and eigenvectors of closed-loop control system by use of either state or output feedback. Applied to reduced-order finite-element mathematical model of NASA's MAST truss beam structure. Model represents deployer/retractor assembly, inertial properties of Space Shuttle, and rigid platforms for allocation of sensors and actuators. Algorithm formulated in real arithmetic for efficient implementation. Choice of open-loop eigenvector matrix and its closest unitary matrix believed suitable for generating well-conditioned eigensystem with small control gains. Implication of this approach is that element of iterative search for "optimal" unitary matrix appears unnecessary in practice for many test problems.

Juang, Jer-Nan↗

High-order state space simulation models of helicopter flight mechanics

This paper describes the formulation and validation of a high-order linearized mathematical model of helicopter flight mechanics, which includes rotor flap and lag degrees of freedom as well as inflow dynamics. The model is extracted numerically from an existing nonlinear, blade element, time simulation model. Extensive modifications in the formulation and solution process of the nonlinear model, required for a theoretically rigorous linearization, are described in detail. The validation results show that the linearized model successfully captures the coupled rotor-fuselage dynamics in the frequency band most critical for the design of advanced flight control systems. Additional results quantify the extent to which the order of the model can be reduced without loss of fidelity.

Kim, Frederick D.↗

Application of unsymmetric block Lanczos vectors in system identification

This paper demonstrates a new system identification approach of using Lanczos coordinates in place of modal coordinates. Identified experimental Lanczos vectors can be directly used in many structural dynamics analysis applications. A multi-input, multi-output frequency-domain technique was used to extract system matrices and an unsymmetric block Lanczos algorithm was used to reduce the order of the experimental model. A cantilever beam example showed promising results, indicating that a new system identification approach using Lanczos coordinates is worthy of further study.

Kim, H. M., Jr.↗

High-order state space simulation models of helicopter flight mechanics

This paper describes the formulation and validation of a high-order linearized mathematical model of helicopter flight mechanics, which includes rotor flap and lag degrees of freedom as well as inflow dynamics. The model is extracted numerically from an existing nonlinear, blade element, real-time simulation model. Extensive modifications in the formulation and solution process of the nonlinear model, required for a theoetically rigorous linearization, are described in detail. The validation results show that the linearized model successfully captures the coupled rotor-fuselage dynamics in the frequency band most critical for the design of advanced flight control systems. Additional results quantify the extent to which the order of the model can be reduced without loss of fidelity.

Kim, Frederick D.↗

Modeling to Improve the Risk Reduction Process for Command File Errors

The Jet Propulsion Laboratory has learned that even innocuous errors in the spacecraft command process can have significantly detrimental effects on a space mission. Consequently, such Command File Errors (CFE), regardless of their effect on the spacecraft, are treated as significant events for which a root cause is identified and corrected. A CFE during space mission operations is often the symptom of imbalance or inadequacy within the system that encompasses the hardware and software used for command generation as well as the human experts and processes involved in this endeavor. As we move into an era of increased collaboration with other NASA centers and commercial partners, these systems become more and more complex. Consequently, the ability to thoroughly model and analyze CFEs formally in order to reduce the risk they pose is increasingly important. In this paper, we summarize the results of applying modeling techniques previously developed to the DAWN flight project. The original models were built with the input of subject matter experts from several flight projects. We have now customized these models to address specific questions for the DAWN flight project and formulating use cases to address their unique mission needs. The goal of this effort is to enhance the project's ability to meet commanding reliability requirements for operations and to assist them in managing their Command File Errors.

spacecraft↗

Confinement, reduced entanglement, and spin-glass order in a random quantum spin-ice model

In this work, we study an effective spin model derived perturbatively from a random transverse-field Ising model on the pyrochlore lattice. The model consists of spin configurations on the pyrochlore lattice, restricted to the spin-ice subspace, with spins interacting with random Ising exchange couplings as well as ring exchanges along the hexagons of the lattice. This model is studied by exact diagonalization up to N=64 site systems. Furthermore, we calculate spin-glass correlation functions and local entanglement entropy S T between spins in a single tetrahedron and the rest of the system. We find that the model undergoes two phase transitions. At weak randomness the model is in a quantum spin-ice phase where S T =ln 6. Increasing randomness, at low transverse fields, first leads to a frozen phase, with long-range spin-glass order and S T =ln 2 corresponding to the Cat states associated with Ising order. Further increase in randomness leads to a random resonating-hexagon phase with a frozen backbone of spins and a broad distribution of entanglement entropies. The implications of these studies for non-Kramers rare-earth pyrochlores are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Learning compact physics‐aware delayed photocurrent models using dynamic mode decomposition

Abstract Radiation‐induced photocurrent in semiconductor devices can be simulated using complex physics‐based models, which are accurate, but computationally expensive. This presents a challenge for implementing device characteristics in high‐level circuit simulations where it is computationally infeasible to evaluate detailed models for multiple individual circuit elements. In this work we demonstrate a procedure for learning compact delayed photocurrent models that are efficient enough to implement in large‐scale circuit simulations, but remain faithful to the underlying physics. Our approach utilizes dynamic mode decomposition (DMD), a system identification technique for learning reduced‐order discrete‐time dynamical systems from time series data based on singular value decomposition. To obtain physics‐aware device models, we simulate the excess carrier density induced by radiation pulses by solving numerically the ambipolar diffusion equation, then use the simulated internal state as training data for the DMD algorithm. Our results show that the significantly reduced‐order delayed photocurrent models obtained via this method accurately approximate the dynamics of the internal excess carrier density—which can be used to calculate the induced current at the device boundaries—while remaining compact enough to incorporate into larger circuit simulations.

Hanson, Joshua↗

Modular design attitude control system

The problem of applying modular attitude control to a rigid body - flexible suspension model of a flexible spacecraft with some state variables inaccessible was addressed by developing a sequence of single axis models and generating a series of reduced state linear observers of minimum order to reconstruct those scalar variables that were inaccessible. The specific single axis models treated consisted of two, three and four rigid bodies, respectively, interconnected by a flexible shaft passing through the mass centers of the bodies. Reduced state linear observers of all orders up to one less than the total number of scalar state variables were generated for each of the three single axis models cited.

Chichester, F. D.↗

Modular design attitude control system

The problem of applying modular attitude control to a rigid body-flexible suspension model of a flexible spacecraft with some state variables inaccessible was addressed by developing a sequence of single axis models and generating a series of reduced state linear observers of minimum order to reconstruct those scalar state variables that were inaccessible. The specific single axis models treated consisted of two, three and four rigid bodies, respectively, interconnected by a flexible shaft passing through the mass centers of the bodies. Reduced state linear observers of all orders up to one less than the total number of scalar state variables were generated for each of the three single axis models cited.

Chichester, F. D.↗

A new stepwise carbon cycle data assimilation system using multiple data streams to constrain the simulated land surface carbon cycle

Here, large uncertainties in land surface models (LSMs) simulations still arise from inaccurate forcing, poor description of land surface heterogeneity (soil and vegetation properties), incorrect model parameter values and incomplete representation of biogeochemical processes. The recent increase in the number and type of carbon cycle-related observations, including both in situ and remote sensing measurements, has opened a new road to optimize model parameters via robust statistical model–data integration techniques, in order to reduce the uncertainties of simulated carbon fluxes and stocks. In this study we present a carbon cycle data assimilation system that assimilates three major data streams, namely the Moderate Resolution Imaging Spectroradiometer (MODIS)-Normalized Difference Vegetation Index (NDVI) observations of vegetation activity, net ecosystem exchange (NEE) and latent heat (LE) flux measurements at more than 70 sites (FLUXNET), as well as atmospheric CO 2 concentrations at 53 surface stations, in order to optimize the main parameters (around 180 parameters in total) of the Organizing Carbon and Hydrology in Dynamics Ecosystems (ORCHIDEE) LSM (version 1.9.5 used for the Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations). The system relies on a stepwise approach that assimilates each data stream in turn, propagating the information gained on the parameters from one step to the next. Overall, the ORCHIDEE model is able to achieve a consistent fit to all three data streams, which suggests that current LSMs have reached the level of development to assimilate these observations. The assimilation of MODIS-NDVI (step 1) reduced the growing season length in ORCHIDEE for temperate and boreal ecosystems, thus decreasing the global mean annual gross primary production (GPP). Using FLUXNET data (step 2) led to large improvements in the seasonal cycle of the NEE and LE fluxes for all ecosystems (i.e., increased amplitude for temperate ecosystems). The assimilation of atmospheric CO 2 , using the general circulation model (GCM) of the Laboratoire de Météorologie Dynamique (LMDz; step 3), provides an overall constraint (i.e., constraint on large-scale net CO 2 fluxes), resulting in an improvement of the fit to the observed atmospheric CO 2 growth rate. Thus, the optimized model predicts a land C (carbon) sink of around 2.2 PgC yr -1 (for the 2000–2009 period), which is more compatible with current estimates from the Global Carbon Project (GCP) than the prior value. The consistency of the stepwise approach is evaluated with back-compatibility checks. The final optimized model (after step 3) does not significantly degrade the fit to MODIS-NDVI and FLUXNET data that were assimilated in the first two steps, suggesting that a stepwise approach can be used instead of the more “challenging” implementation of a simultaneous optimization in which all data streams are assimilated together. Most parameters, including the scalar of the initial soil carbon pool size, changed during the optimization with a large error reduction. This work opens new perspectives for better predictions of the land carbon budgets.

58 GEOSCIENCES↗

Burst Pressure Solutions of Thin and Thick-Walled Cylindrical Vessels

Pressure vessels (PVs) are widely used in the energy industry. Accurate burst pressure is critical to structural design and safe operation for both thin and thick-walled PVs. The traditional strength theories utilized a single-parameter material property, such as the yield stress or the ultimate tensile stress (UTS) to develop failure models for determining the yield or ultimate pressure carrying capacity in the PV design. The UTS-based Barlow formula is a typical burst pressure model developed from the Tresca strength theory that provides the basis for developing regulation rules and failure models for different industry design codes, such as ASME BPVC, ASME B31.3, and ASME B31G, among others. In order to reduce the conservatism of the Tresca strength model, ASME BPVC recently adapted failure models developed from the von Mises strength theory for the PV design and analysis. It has been commonly accepted that the burst pressure of pipelines depends on the UTS and strain hardening exponent, n, of the pipeline steel. An average shear stress yield theory was thus developed, and the Zhu-Leis solution of burst pressure was obtained as a function of UTS and n for thin-walled line pipes. Experiments showed that the Zhu-Leis solution provides an accurate, reliable prediction of burst pressure for defect-free thin-walled pipes. In order to extend the Zhu–Leis solution to thick-walled cylindrical PVs, this paper defined three new flow stresses, modified the traditional strength theories, and obtained three new burst pressure solutions that are valid for both thin and thick-walled cylindrical vessels. The proposed flow stresses are able to describe the tensile strength and the plastic flow response of PVs for a strain hardening steel. The associated strength theories were then developed in terms of the Tresca, von Mises, and Zhu-Leis yield criteria. From these new strength theories, three burst pressure solutions were obtained for thick-walled cylinders, where the von Mises solution is an upper bound prediction, the Tresca solution is a lower bound prediction, and the Zhu-Leis solution is an intermediate prediction of burst pressure for thick-walled cylinders. Lastly, the proposed burst pressure solutions were evaluated and validated by two large datasets of full-scale burst tests for thick-walled tubes and for thin-walled pipes.

42 ENGINEERING↗

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗

OPTICON: Pro-Matlab software for large order controlled structure design

A software package for large order controlled structure design is described and demonstrated. The primary program, called OPTICAN, uses both Pro-Matlab M-file routines and selected compiled FORTRAN routines linked into the Pro-Matlab structure. The program accepts structural model information in the form of state-space matrices and performs three basic design functions on the model: (1) open loop analyses; (2) closed loop reduced order controller synthesis; and (3) closed loop stability and performance assessment. The current controller synthesis methods which were implemented in this software are based on the Generalized Linear Quadratic Gaussian theory of Bernstein. In particular, a reduced order Optimal Projection synthesis algorithm based on a homotopy solution method was successfully applied to an experimental truss structure using a 58-state dynamic model. These results are presented and discussed. Current plans to expand the practical size of the design model to several hundred states and the intention to interface Pro-Matlab to a supercomputing environment are discussed.

Peterson, Lee D.↗

Computation of Turbulent Wake Flows in Variable Pressure Gradient

Transport aircraft performance is strongly influenced by the effectiveness of high-lift systems. Developing wakes generated by the airfoil elements are subjected to strong pressure gradients and can thicken very rapidly, limiting maximum lift. This paper focuses on the effects of various pressure gradients on developing symmetric wakes and on the ability of a linear eddy viscosity model and a non-linear explicit algebraic stress model to accurately predict their downstream evolution. In order to reduce the uncertainties arising from numerical issues when assessing the performance of turbulence models, three different numerical codes with the same turbulence models are used. Results are compared to available experimental data to assess the accuracy of the computational results.

Duquesne, N.↗

Using Co-Optimized Machine Learned Manifolds for Modeling Chemically Reacting Flows

Chemically reacting flows play a key role in a wide range of engineered systems, from chemical and polymer processing to combustion-based energy conversion technologies. Simulations of these flows involve solving a coupled set of partial differential equations for mass, momentum, energy, and all relevant chemical species in the system. Chemical reaction pathways may be extremely complex and involve hundreds or more intermediate species, with reactions that occur over timescales varying by several orders of magnitude - presenting a significant numerical stiffness challenge. The combination of these factors makes simulation of chemically reacting flows vastly more expensive than nonreactive simulations, and often makes direct solution of the governing equations intractable. It is necessary to apply lower-fidelity models in place of the detailed governing equations in order to reduce computational cost to enable reacting flow simulation tools to be used in the engineering design process. Many of the models employed for this purpose are based on reducing the dimension of the thermochemical state, motivated by the observation that the observed thermochemical states in a system lie on a low-dimensional manifold in thermochemical state space. This behavior occurs due to the fast equilibration of certain reactive and transport processes, and physics-based manifold models rely on idealized assumptions about the balance of timescales and the way in which chemistry and transport are coupled. In this work, we apply a novel method for data-driven manifold-based modeling that can leverage data from high-fidelity reacting flow simulations to improve model accuracy in cases where the physics-based modeling assumptions break down. The approach is designed to be broadly applicable across chemically reacting flow systems but is applied here to turbulent combustion modeling.

machine learning↗