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,477 records · Page 82

A study of pilot modeling in multi-controller tasks

A modeling approach, which utilizes a matrix of transfer functions to describe the human pilot in multiple input, multiple output control situations, is studied. The approach used was to extend a well established scalar Wiener-Hopf minimization technique to the matrix case and then study, via a series of experiments, the data requirements when only finite record lengths are available. One of these experiments was a two-controller roll tracking experiment designed to force the pilot to use rudder in order to coordinate and reduce the effects of aileron yaw. One model was computed for the case where the signals used to generate the spectral matrix are error and bank angle while another model was computed for the case where error and yaw angle are the inputs. Several anomalies were observed to be present in the experimental data. These are defined by the descriptive terms roll up, break up, and roll down. Due to these algorithm induced anomalies, the frequency band over which reliable estimates of power spectra can be achieved is considerably less than predicted by the sampling theorem.

Whitbeck, R. F.↗

A point-particle-based hydride shell-shedding model for ejecta particle transport in reactive environments

A shock wave passing over a rough or perturbed metal surface will induce a limiting case of Richtmyer–Meshkov instability and will cause small particles to eject from the surface and transport into the surrounding medium. These particles are known as ejecta and can be either solid or liquid in nature. Recent experiments have shown that liquid cerium ejecta clouds exhibit unexpected non-monotonic acceleration behaviors as well as temperature plateaus after a brief temperature rise if they are transporting in a chemically reactive, hydrogen-based medium while they act as expected in an inert medium. This work details a point-particle model developed for reactive cerium ejecta transport, which attempts to account for these new physics through the behavior of a developing solid hydride shell, which is believed to form as a product of the reaction. The overall model incorporates the effects of the reaction on the particle properties as well as the effects of potential shedding of the shell into sub-micrometer scale flakes and potential phase change of the hydride if the ejecta particles reach the melt point of the hydride layer. The model is tested by performing simulations of the original motivating experiments and comparing quantities, such as ejected mass, velocimetry, and temperature profiles, against the experimental data. While the model is able to capture many general features of the observed anomalies, some inaccuracies still exist. These point to both missing physics in the model (such as a deuterium adsorption mechanism on the hydride layer) as well as a lack of knowledge of certain material properties (such as the strength of cerium hydride to determine dynamic fracture thicknesses) needed to fully reduce the uncertainties in the model by up to an order of magnitude and perform a true attempt at model validation.

97 MATHEMATICS AND COMPUTING↗

LDRD-DR: Accelerating Nuclear Fuel Qualification through Integrated Multiscale and Multiphysics Models [Slides]

Nuclear fuel qualification implies ensuring that a new fuel can safely operate in a reactor under both normal and accident conditions. Historically, extensive irradiation campaigns probe fuel performance relationships and interconnections in an integral sense, which is very time-consuming and costly. The presentation focus is on introducing advanced modeling and simulation into the qualification process in order to reduce the qualification timeline by targeting fewer integral tests

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

ACCELERATED FITTING OF STELLAR SPECTRA

Stellar spectra are often modeled and fitted by interpolating within a rectilinear grid of synthetic spectra to derive the stars’ labels: stellar parameters and elemental abundances. However, the number of synthetic spectra needed for a rectilinear grid grows exponentially with the label space dimensions, precluding the simultaneous and self-consistent fitting of more than a few elemental abundances. Shortcuts such as fitting subsets of labels separately can introduce unknown systematics and do not produce correct error covariances in the derived labels. In this paper we present a new approach—Convex Hull Adaptive Tessellation (chat)—which includes several new ideas for inexpensively generating a sufficient stellar synthetic library, using linear algebra and the concept of an adaptive, data-driven grid. A convex hull approximates the region where the data lie in the label space. A variety of tests with mock data sets demonstrate that chat can reduce the number of required synthetic model calculations by three orders of magnitude in an eight-dimensional label space. The reduction will be even larger for higher dimensional label spaces. In chat the computational effort increases only linearly with the number of labels that are fit simultaneously. Around each of these grid points in the label space an approximate synthetic spectrum can be generated through linear expansion using a set of “gradient spectra” that represent flux derivatives at every wavelength point with respect to all labels. These techniques provide new opportunities to fit the full stellar spectra from large surveys with 15–30 labels simultaneously.

79 ASTRONOMY AND ASTROPHYSICS↗

Foreground bias from parametric models of far-IR dust emission

We use simple toy models of far-IR dust emission to estimate the accuracy to which the polarization of the cosmic microwave background can be recovered using multi-frequency fits, if the parametric form chosen for the fitted dust model differs from the actual dust emission. Commonly used approximations to the far-IR dust spectrum yield CMB residuals comparable to or larger than the sensitivities expected for the next generation of CMB missions, despite fitting the combined CMB + foreground emission to precision 0.1% or better. The Rayleigh–Jeans approximation to the dust spectrum biases the fitted dust spectral index by Δβ{sub d}=0.2 and the inflationary B-mode amplitude by Δr=0.03. Fitting the dust to a modified blackbody at a single temperature biases the best-fit CMB by Δr>0.003 if the true dust spectrum contains multiple temperature components. A 13-parameter model fitting two temperature components reduces this bias by an order of magnitude if the true dust spectrum is in fact a simple superposition of emission at different temperatures, but fails at the level Δr=0.006 for dust whose spectral index varies with frequency. Restricting the observing frequencies to a narrow region near the foreground minimum reduces these biases for some dust spectra but can increase the bias for others. Data at THz frequencies surrounding the peak of the dust emission can mitigate these biases while providing a direct determination of the dust temperature profile.

79 ASTRONOMY AND ASTROPHYSICS↗

Performance improvement of robots using a learning control scheme

Many applications of robots require that the same task be repeated a number of times. In such applications, the errors associated with one cycle are also repeated every cycle of the operation. An off-line learning control scheme is used here to modify the command function which would result in smaller errors in the next operation. The learning scheme is based on a knowledge of the errors and error rates associated with each cycle. Necessary conditions for the iterative scheme to converge to zero errors are derived analytically considering a second order servosystem model. Computer simulations show that the errors are reduced at a faster rate if the error rate is included in the iteration scheme. The results also indicate that the scheme may increase the magnitude of errors if the rate information is not included in the iteration scheme. Modification of the command input using a phase and gain adjustment is also proposed to reduce the errors with one attempt. The scheme is then applied to a computer model of a robot system similar to PUMA 560. Improved performance of the robot is shown by considering various cases of trajectory tracing. The scheme can be successfully used to improve the performance of actual robots within the limitations of the repeatability and noise characteristics of the robot.

Krishna, Ramuhalli↗

Computational modelling of dump combustors flowfield

A computational model aimed at predicting the flowfield of dump combustors is presented. The turbulent combustion model is based on the conserved scalar approach and on a convenient specification of its probability density function, which reduces the computation of the mean density to a closed form. Turbulence is modeled by means of the k-epsilon model. The averaged conservation equations are solved by a technique based on a staggered grid and on the SIMPLE solver. The computational model is applied to a simple dump combustor to assess the computer time requirements and accuracy. The turbulent combustion model is shown to reduce the computer time by an order of magnitude when compared to evaluating the mean density by numerical quadrature.

Lentini, D.↗

The three-point function in an ensemble of numerical simulations

The three-point function in Fourier space is evaluated for an ensemble of 2D numerical simulations with nine initial power spectra distinguished by spectral index and cutoff. To remove the main dependences on scale and time, results are presented as the reduced amplitude Q in the hierarchical model. To lowest nonvanishing order in perturbation theory, Q is a constant, independent of length scale, of time, and of initial spectrum. The results show that in the nonlinear regime of evolution, normalizing to Q does remove the main variations, but systematic dependences on spectral index and cutoff remain at late times that depend on properties of the initial spectrum. Models with initial power cut off on a scale that has been evolved far into the nonlinear regime are indistinguishable from models with no cutoff.

Fry, J. N.↗

Foreground Bias from Parametric Models of Far-IR Dust Emission

We use simple toy models of far-IR dust emission to estimate the accuracy to which the polarization of the cosmic microwave background can be recovered using multi-frequency fits, if the parametric form chosen for the fitted dust model differs from the actual dust emission. Commonly used approximations to the far-IR dust spectrum yield CMB residuals comparable to or larger than the sensitivities expected for the next generation of CMB missions, despite fitting the combined CMB plus foreground emission to precision 0.1 percent or better. The Rayleigh-Jeans approximation to the dust spectrum biases the fitted dust spectral index by (Delta)(Beta)(sub d) = 0.2 and the inflationary B-mode amplitude by (Delta)(r) = 0.03. Fitting the dust to a modified blackbody at a single temperature biases the best-fit CMB by (Delta)(r) greater than 0.003 if the true dust spectrum contains multiple temperature components. A 13-parameter model fitting two temperature components reduces this bias by an order of magnitude if the true dust spectrum is in fact a simple superposition of emission at different temperatures, but fails at the level (Delta)(r) = 0.006 for dust whose spectral index varies with frequency. Restricting the observing frequencies to a narrow region near the foreground minimum reduces these biases for some dust spectra but can increase the bias for others. Data at THz frequencies surrounding the peak of the dust emission can mitigate these biases while providing a direct determination of the dust temperature profile.

cosmic microwave background↗

Predictive performance modeling framework for a novel enclosed particle receiver configuration and application for thermochemical energy storage

Concentrating solar power (CSP) plants can provide dispatchable power with a thermal energy storage capability for increased renewable-energy grid penetration. Particle-based CSP systems permit higher temperatures, and thus, potentially higher solar-to-electric efficiency than state-of-the-art molten-salt heat-transfer systems. This paper describes a detailed numerical analysis framework for estimating the performance of a novel, geometrically complex, enclosed particle receiver design. The receiver configuration uses arrays of small tubular absorbers to collect and subsequently transfer solar energy to a flowing particulate medium. The enclosed nature of the receiver design renders it amenable to either an inert heat-transfer medium, or a reactive heat-transfer medium that requires a controllable ambient environment. The numerical analysis framework described in this study is demonstrated for the case of thermal reduction of CaCr 0.1 Mn 0.9 O 3-$\delta$ for thermochemical energy storage. The modeling strategy consists of Monte Carlo ray tracing for absorbed solar-energy distributions from a surround heliostat field, computational fluid dynamics modeling of small-scale local tubular arrays, surrogate response surfaces that approximately capture simulated tubular array performance, a quasi-two-dimensional reduced-order description of counter-flow reactive solids and purge gas, and a radiative exchange model applied to embedded-cavity structures at the size scale of the full receiver. In this work we apply the numerical analysis strategy to a single receiver configuration, but the framework can be generically applicable to alternative enclosed designs. In conclusion, we assess sensitivity of receiver performance to surface optical properties, heat-transfer coefficients, solids outlet temperature, and purge-gas feed rates, and discuss the significance of model assumptions and results for future receiver development.

14 SOLAR ENERGY↗

The ANMLite Language and Logic for Specifying Planning Problems

We present the basic concepts of the ANMLite planning language. We discuss various aspects of specifying a plan in terms of constraints and checking the existence of a solution with the help of a model checker. The constructs of the ANMLite language have been kept as simple as possible in order to reduce complexity and simplify the verification problem. We illustrate the language with a specification of the space shuttle crew activity model that was constructed under the Spacecraft Autonomy for Vehicles and Habitats (SAVH) project. The main purpose of this study was to explore the implications of choosing a robust logic behind the specification of constraints, rather than simply proposing a new planning language.

Butler, Ricky W.↗

Cycle-averaged dynamics of a periodically driven, closed-loop circulation model

Time-varying elastance models have been used extensively in the past to simulate the pulsatile nature of cardiovascular waveforms. Frequently, however, one is interested in dynamics that occur over longer time scales, in which case a detailed simulation of each cardiac contraction becomes computationally burdensome. In this paper, we apply circuit-averaging techniques to a periodically driven, closed-loop, three-compartment recirculation model. The resultant cycle-averaged model is linear and time invariant, and greatly reduces the computational burden. It is also amenable to systematic order reduction methods that lead to further efficiencies. Despite its simplicity, the averaged model captures the dynamics relevant to the representation of a range of cardiovascular reflex mechanisms. c2004 Elsevier Ltd. All rights reserved.

Non-NASA Center↗

Turbulence modeling of shock separated boundary-layer flows

Computations of transonic and hypersonic shock-separated boundary-layer flows using zero-equation (algebraic), one-equation (kinetic energy), and two-equation (kinetic energy plus length scale) turbulence eddy viscosity models are described and compared with measurements. The computations make use of a new Navier-Stokes computer algorithm that has reduced computing times by one to two orders of magnitude. The algorithm, and how the turbulence models are incorporated into it, are described. Results for the transonic flow show that the unmodified one-equation model is superior to the zero-equation model in skin-friction predictions. For the hypersonic flow, a highly modified one-equation model that accurately predicts surface pressure and heat transfer is described. Preliminary two-equation model results are also presented.

Coakley, T. J.↗

Model Reduction for Fast Assessment of Grid Impact of High Penetration PV

A key challenge in properly analyzing the impact of solar photovoltaics (PV) during planning, interconnection and operation in the distribution grid is computational complexity. Such studies require the use of Quasi-Steady-Time-Series (QSTS) power flow to capture the complex interaction of PV with voltage regulating equipment and protection devices and a large number of scenario studies are needed to capture location specific constraints and uncertainties. In addition, a typical distribution feeder can have thousands of nodes and there can be thousands of such feeders in a given utility system. To solve this problem, a new and systematic method for distribution network reduction is proposed. The objective is to reduce the size of distribution systems by more than 80% while preserving the accuracy of power flow results. To this end, an extended Ward-equivalent method has been proposed to create an equivalent circuit. Moreover, load and PV aggregation based on the backward sweep method has been developed to further simplify the feeder. The proposed method is tested with a modified IEEE 123-node distribution test system. The simulation results demonstrate the high efficiency and accuracy of the proposed approach for distribution network reduction.

reduced order systems, load flow, load modeling, v↗

Thermal design, testing, and analysis of the Cosmic Background Explorer's external calibrator

The thermal design of a calibrator used to correlate the collected data with an absolute reference onboard of the Cosmic Background Explorer is considered. Since flight conditions could not be simulated in ground testing, the design was qualified by analysis and testing. The requirements imposed on the calibrator are outlined, along with spacecraft hardware modifications aimed at reducing the thermal load on the calibrator. A geometric mathematical model created in order to obtain view factors from the cryogenic instruments to the spacecraft and dewar is discussed. It is shown that the calibrator met or exceeded all of its performance requirements, and that the modifications to such energy sources as the thermal shield and differential microwave radiometers proved effective in reducing the radiative energy the calibrator received from the spacecraft to acceptable levels.

Mosier, Carol L.↗

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗

Gaussian approximation potential for amorphous Si : H

Hydrogenation of amorphous silicon (a–Si : H) is critical for reducing defect densities, passivating midgap states and surfaces, and improving photoconductivity in silicon-based electro-optical devices. Modeling the atomic-scale structure of this material is critical to understanding these processes, which in turn is needed to describe c–Si/a–Si : H heterojunctions that are at the heart of modern solar cells with world-record efficiency. Density functional theory (DFT) studies achieve the required high accuracy but are limited to moderate system sizes of 100 atoms or so by their high computational cost. Simulations of amorphous materials have been hindered by this high cost because large structural models are required to capture the medium-range order that is characteristic of such materials. Empirical potential models are much faster, but their accuracy is not sufficient to correctly describe the frustrated local structure. Data-driven, machine-learned interatomic potentials have broken this impasse and have been highly successful in describing a variety of amorphous materials in their elemental phase. Here, we extend the Gaussian approximation potential (GAP) for silicon by incorporating the interaction with hydrogen, thereby significantly improving the degree of realism with which amorphous silicon can be modeled. We show that our Si : H GAP enables the simulation of hydrogenated silicon with an accuracy very close to DFT but with computational expense and run times reduced by several orders of magnitude for large structures. Here, we demonstrate the capabilities of the Si : H GAP by creating models of hydrogenated liquid and amorphous silicon and showing that their energies, forces, and stresses are in excellent agreement with DFT results, and their structure as captured by bond and angle distributions are in agreement with both DFT and experiments.

36 MATERIALS SCIENCE↗

Variational quantum simulation of the critical Ising model with symmetry averaging

Here we investigate the use of deep multiscale entanglement renormalization ansatz (DMERA) circuits as a variational ansatz. We use the exactly solvable one-dimensional critical transverse-field Ising model as a test bed. Numerically exact simulation of the quantum circuit ansatz can in this case be carried out to hundreds of qubits by exploiting efficient classical algorithms for simulating matchgate circuits. We find that, for this system, the DMERA strongly outperforms a standard quantum approximate optimization algorithm (QAOA)–style ansatz, and that a major source of systematic error in correlation functions approximated using the DMERA is the breaking of the translational and Kramers-Wannier symmetries of the transverse-field Ising model. We are able to reduce this error by up to four orders of magnitude by symmetry averaging, without incurring additional cost in qubits or circuit depth. Here, we propose that this technique for mitigating systematic error could be applied to noisy intermediate-scale quantum (NISQ) simulations of physical systems with other symmetries.

1-dimensional spin chains↗