Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model order reduction”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows

In the present work, a general formulation is proposed to implement the contact angle boundary conditions for the second-order Phase-Field models, which is applicable to N-phase (N ≥ 2) moving contact line problems. To remedy the issue of mass change due to the contact angle boundary condition, a source term or Lagrange multiplier is added to the original second-order Phase-Field models, which is determined by the consistent and conservative volume distribution algorithm so that the summation of the order parameters and the consistency of reduction are not influenced. To physically couple the proposed formulation to the hydrodynamics, especially for large-density-ratio problems, the consistent formulation is employed. The reduction-consistent conservative Allen-Cahn models are chosen as examples to illustrate the application of the proposed formulation. The numerical scheme that preserves the consistency and conservation of the proposed formulation is employed to demonstrate its effectiveness. Results produced by the proposed formulation are in good agreement with the exact and/or asymptotic solutions. The proposed method captures complex dynamics of moving contact line problems having large density ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Bayes goes fast: Uncertainty quantification for a covariant energy density functional emulated by the reduced basis method

A covariant energy density functional is calibrated using a principled Bayesian statistical framework informed by experimental binding energies and charge radii of several magic and semi-magic nuclei. The Bayesian sampling required for the calibration is enabled by the emulation of the high-fidelity model through the implementation of a reduced basis method (RBM)—a set of dimensionality reduction techniques that can speed up demanding calculations involving partial differential equations by several orders of magnitude. The RBM emulator we build—using only 100 evaluations of the high-fidelity model—is able to accurately reproduce the model calculations in tens of milliseconds on a personal computer, an increase in speed of nearly a factor of 3,300 when compared to the original solver. Besides the analysis of the posterior distribution of parameters, we present model calculations for masses and radii with properly estimated uncertainties. We also analyze the model correlation between the slope of the symmetry energy L and the neutron skin of 48 Ca and 208 Pb. The straightforward implementation and outstanding performance of the RBM makes it an ideal tool for assisting the nuclear theory community in providing reliable estimates with properly quantified uncertainties of physical observables. Such uncertainty quantification tools will become essential given the expected abundance of data from the recently inaugurated and future experimental and observational facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1994

Local dynamic subgrid-scale models in channel flow

The dynamic subgrid-scale (SGS) model has given good results in the large-eddy simulation (LES) of homogeneous isotropic or shear flow, and in the LES of channel flow, using averaging in two or three homogeneous directions (the DA model). In order to simulate flows in general, complex geometries (with few or no homogeneous directions), the dynamic SGS model needs to be applied at a local level in a numerically stable way. Channel flow, which is inhomogeneous and wall-bounded flow in only one direction, provides a good initial test for local SGS models. Tests of the dynamic localization model were performed previously in channel flow using a pseudospectral code and good results were obtained. Numerical instability due to persistently negative eddy viscosity was avoided by either constraining the eddy viscosity to be positive or by limiting the time that eddy viscosities could remain negative by co-evolving the SGS kinetic energy (the DLk model). The DLk model, however, was too expensive to run in the pseudospectral code due to a large near-wall term in the auxiliary SGS kinetic energy (k) equation. One objective was then to implement the DLk model in a second-order central finite difference channel code, in which the auxiliary k equation could be integrated implicitly in time at great reduction in cost, and to assess its performance in comparison with the plane-averaged dynamic model or with no model at all, and with direct numerical simulation (DNS) and/or experimental data. Other local dynamic SGS models have been proposed recently, e.g., constrained dynamic models with random backscatter, and with eddy viscosity terms that are averaged in time over material path lines rather than in space. Another objective was to incorporate and test these models in channel flow.

Cabot, William H.↗
NASA NTRS1979

Acoustic duct liner optimization using finite elements

Determination of a mathematical model which may be used for obtaining optimum acoustic liners in aeroengine ducts was the goal of this study. Starting with a baseline finite element model which consumed too much computer time to permit optimum liner evaluation for realistic ducts and source frequencies, a survey was conducted of candidate methods for reducing computational effort. A highly effective boundary exclusive decomposition method was developed which enables L-U decomposition of a major partition of the global matrix to be performed once only. Timing experiments using the technique in duct liner optimization show a reduction of between one and two orders of magnitude in the central processor unit (CPU) time of the baseline model.

Abrahamson, A. L.↗
NASA NTRS2006

Utilizing a Common Model Architecture for the Simulation of On-Orbit Human Spaceflight Operations

Advances in computing and simulation technology during the 1980 s promoted the development of multiple simulations for on-orbit human spaceflight operations at the National Aeronautics and Space Administration (NASA) Johnson Space Center (JSC). By the late 1980 s, it became increasingly clear that redundancy risked leading to inefficiencies in dissemination and incorporation of models across simulations, institutional organizations, and spaceflight programs. Software development tools, such as the Trick Simulation Environment, began to emerge providing a generalized framework for the development, integration, and operation of simulations. Throughout the 1990s, engineering, operations, and training simulations began to migrate to Trick in order to take advantage of a common simulation environment. Due to forecasted reductions in funding to the Shuttle and International Space Station (ISS) programs in the early 2000s, the common framework that arose from this simulation migration was formalized as the Common Model Architecture (CMA). The CMA was established in order to further reduce redundancy, while promoting model sharing and enhancing model integrity. A minimal set of standards for data flows, functional interactions, and model organization based on Trick were defined for those features found to be most important in facilitating model interaction. Model-providing organizations within NASA JSC took ownership of models for their specific areas of expertise, and agreements were reached on common model sharing and usage. This paper describes the evolution of commonality enabled by the Trick Simulation Environment at JSC leading to the CMA, the standards defined by the CMA, a recent CMA simulation example, and potential application of the CMA to NASA s new Exploration program.

Leslie J Quiocho↗
NASA NTRS2001

Parametric Study of a YAV-8B Harrier in Ground Effect Using Time-Dependent Navier-Stokes Computations

Flow simulations using the time-dependent Navier-Stokes equations remain a challenge for several reasons. Principal among them are the difficulty to accurately model complex flows, and the time needed to perform the computations. A parametric study of such complex problems is not considered practical due to the large cost associated with computing many time-dependent solutions. The computation time for each solution must be reduced in order to make a parametric study possible. With successful reduction of computation time, the issue of accuracy, and appropriateness of turbulence models will become more tractable.

Shishir, Pandya↗
NASA NTRS1994

Flutter suppression digital control law design and testing for the AFW wind tunnel model

The design of a control law for simultaneously suppressing the symmetric and antisymmetric flutter modes of a sting mounted fixed-in-roll aeroelastic wind-tunnel model is described. The flutter suppression control law was designed using linear quadratic Gaussian theory, and it also involved control law order reduction, a gain root-locus study, and use of previous experimental results. A 23 percent increase in the open-loop flutter dynamic pressure was demonstrated during the wind-tunnel test. Rapid roll maneuvers at 11 percent above the symmetric flutter boundary were also performed when the model was in a free-to-roll configuration.

Mukhopadhyay, Vivek↗
NASA NTRS1978

Characterization of the effects of metallic impurities on silicon solar cell performance

The effects of controlled concentrations of secondary impurities (Al, C, Ca, Cr, Cu, Fe, Mg, Mn, Mo, Ta, Ti, V, Zn, and Zr), incorporated alone or in combinations into Czochralski and float zone crystals, on the performance of silicon solar cells were investigated. Impurity concentrations ranged from 10 to the 11th to 10 to the 17th/cu cm. Solar cells were fabricated by a conventional diffusion process and characterized by computer reduction of current-voltage data. The data suggest that performance loss primarily results from reduction of the base diffusion length. On the basis of this assumption, a first order analytic model which predicts cell performance as a function of impurity concentrations is developed. Calculated performance parameters are in good agreement with observation, except for some impurities, such as Fe, Cu and Ni, which degrade cells via recombination and defects in the junction space-charge region. N-base devices are less affected by impurities, although degradation mechanisms appear to be the same as in p-devices. There seems to be very limited interaction between impurities in multiply-doped cells.

Davis, J. R.↗
NASA NTRS1981

Computational methods of robust controller design for aerodynamic flutter suppression

The development of Riccati iteration, a tool for the design and analysis of linear control systems is examined. First, Riccati iteration is applied to the problem of pole placement and order reduction in two-time scale control systems. Order reduction, yielding a good approximation to the original system, is demonstrated using a 16th order linear model of a turbofan engine. Next, a numerical method for solving the Riccati equation is presented and demonstrated for a set of eighth order random examples. A literature review of robust controller design methods follows which includes a number of methods for reducing the trajectory and performance index sensitivity in linear regulators. Lastly, robust controller design for large parameter variations is discussed.

Anderson, L. R.↗
NASA NTRS1989

Computational design of low aspect ratio wing-winglet configurations for transonic wind-tunnel tests

Computational designs were performed for three different low aspect ratio wing planforms fitted with nonplanar winglets; one of the three configurations was selected to be constructed as a wind tunnel model for testing in the NASA LaRC 8-foot transonic pressure tunnel. A design point of M = 0.8, C(sub L) is approximate or = to 0.3 was selected, for wings of aspect ratio equal to 2.2, and leading edge sweep angles of 45 deg and 50 deg. Winglet length is 15 percent of the wing semispan, with a cant angle of 15 deg, and a leading edge sweep of 50 deg. Winglet total area equals 2.25 percent of the wing reference area. The design process and the predicted transonic performance are summarized for each configuration. In addition, a companion low-speed design study was conducted, using one of the transonic design wing-winglet planforms but with different camber and thickness distributions. A low-speed wind tunnel model was constructed to match this low-speed design geometry, and force coefficient data were obtained for the model at speeds of 100 to 150 ft/sec. Measured drag coefficient reductions were of the same order of magnitude as those predicted by numerical subsonic performance predictions.

Kuhlman, John M.↗
NASA NTRS1992

Model reduction using new optimal Routh approximant technique

An optimal Routh approximant of a single-input single-output dynamic system is a reduced-order transfer function of which the denominator is obtained by the Routh approximation method while the numerator is determined by minimizing a time-response integral-squared-error (ISE) criterion. In this paper, a new elegant approach is presented for obtaining the optimal Routh approximants for linear time-invariant continuous-time systems. The approach is based on the Routh canonical expansion, which is a finite-term orthogonal series of rational basis functions, and minimization of the ISE criterion. A procedure for combining the above approach with the bilinear transformation is also presented in order to obtain the optimal bilinear Routh approximants of linear time-invariant discrete-time systems. The proposed technique is simple in formulation and is amenable to practical implementation.

Hwang, Chyi↗
NASA NTRS1995

Synthesis of low color, atomic oxygen resistant polyimides

The purpose of this project was to develop low color, atomic oxygen resistant polyimides for potential applications on spacecraft in low earth orbit. The material is needed in order to protect satellites and spacecraft from the gases and radiation found at those altitudes. Phosphorous containing polyimides have been shown to be especially resistant to corrosion and weight loss under oxygen plasma. Unfortunately the color of these phosphorous containing polyimides is still too high for them to be good heat insulators. While they are not as effective as teflon, the current material of choice. polyimides are much less dense than teflon and would be especially valuable if they could be made with low color. The approach taken was to synthesize a monomer which would contain the elements needed for giving the final polyimide its desired properties. In particular the monomer should incorporate phosphine or phosphine oxides and have bulky side groups to block any color forming charge transfer structures. The target molecule, 3,5-di-(trifluoromethylphenyl)-bis(3-aminophenyl) phosphine oxide, (containing both a phosphine oxide group and a bulky ditrifluoromethylphenyl group) was synthesized via three reactions in overall yield of 21 percent. In addition, a model compound, bis(3-phenylamine) phenyl phosphine oxide, was synthesized two different ways in order to establish the conditions for the nitration of phosphine oxides and their reduction to the amine. Finally, a trisubstituted phosphine oxide was synthesized. In all, seven phosphorus containing organic compounds were synthesized, purified and characterized. The model compound was reacted with oxydiphthalic anhydride to form a polyamic acid with inherent viscosity of 0.34. This material was cast into a film and heated, forming a normally colored fairly strong polyimide with a Tg of 240 C. The target compound was reacted with 6-fluorodiphthalic anhydride to give a polyamic acid with inherent viscosity of 0.19 and cast to give a heavily cracked colored film with a Tg of 230 C.

MacInnes, Dave↗
DOE OSTI2020

Closure Learning for Nonlinear Model Reduction Using Deep Residual Neural Network

Developing accurate, efficient, and robust closure models is essential in the construction of reduced order models (ROMs) for realistic nonlinear systems, which generally require drastic ROM mode truncations. We propose a deep residual neural network (ResNet) closure learning framework for ROMs of nonlinear systems. The novel ResNet-ROM framework consists of two steps: (i) In the first step, we use ROM projection to filter the given nonlinear system and construct a spatially filtered ROM. This filtered ROM is low-dimensional, but is not closed. (ii) In the second step, we use ResNet to close the filtered ROM, i.e., to model the interaction between the resolved and unresolved ROM modes. We emphasize that in the new ResNet-ROM framework, data is used only to complement classical physical modeling (i.e., only in the closure modeling component), not to completely replace it. We also note that the new ResNet-ROM is built on general ideas of spatial filtering and deep learning and is independent of (restrictive) phenomenological arguments, e.g., of eddy viscosity type. The numerical experiments for the 1D Burgers equation show that the ResNet-ROM is significantly more accurate than the standard projection ROM. The new ResNet-ROM is also more accurate and significantly more efficient than other modern ROM closure models.

Xie, Xuping↗
NASA NTRS1992

Design, test, and evaluation of three active flutter suppression controllers

Three control law design techniques for flutter suppression are presented. Each technique uses multiple control surfaces and/or sensors. The first method uses traditional tools (such as pole/zero loci and Nyquist diagrams) for producing a controller that has minimal complexity and which is sufficiently robust to handle plant uncertainty. The second procedure uses linear combinations of several accelerometer signals and dynamic compensation to synthesize the model rate of the critical mode for feedback to the distributed control surfaces. The third technique starts with a minimum-energy linear quadratic Gaussian controller, iteratively modifies intensity matrices corresponding to input and output noise, and applies controller order reduction to achieve a low-order, robust controller. The resulting designs were implemented digitally and tested subsonically on the active flexible wing wind-tunnel model in the Langley Transonic Dynamics Tunnel. Only the traditional pole/zero loci design was sufficiently robust to errors in the nominal plant to successfully suppress flutter during the test. The traditional pole/zero loci design provided simultaneous suppression of symmetric and antisymmetric flutter with a 24-percent increase in attainable dynamic pressure. Posttest analyses are shown which illustrate the problems encountered with the other laws.

Adams, William M., Jr.↗
NASA NTRS1979

Modeling of switching regulator power stages with and without zero-inductor-current dwell time

State-space techniques are employed to derive accurate models for the three basic switching converter power stages: buck, boost, and buck/boost operating with and without zero-inductor-current dwell time. A generalized procedure is developed which treats the continuous-inductor-current mode without dwell time as a special case of the discontinuous-current mode when the dwell time vanishes. Abrupt changes of system behavior, including a reduction of the system order when the dwell time appears, are shown both analytically and experimentally. Merits resulting from the present modeling technique in comparison with existing modeling techniques are illustrated.

Lee, F. C. Y.↗
DOE OSTI2022

Energetically consistent model reduction for metriplectic systems

The metriplectic formalism is useful for describing complete dynamical systems which conserve energy and produce entropy. This creates challenges for model reduction, as the elimination of high-frequency information will generally not preserve the metriplectic structure which governs long-term stability of the system. Based on proper orthogonal decomposition, a provably convergent metriplectic reduced-order model is formulated which is guaranteed to maintain the algebraic structure necessary for energy conservation and entropy formation. Further, numerical results on benchmark problems show that the proposed method is remarkably stable, leading to improved accuracy over long time scales at a moderate increase in cost over naive methods.

42 ENGINEERING↗
NASA NTRS1995

Gas-phase chemistry in dense interstellar clouds including grain surface molecular depletion and desorption

We present time-dependent models of the chemical evolution of molecular clouds which include depletion of atoms and molecules onto grain surfaces and desorption, as well as gas-phase interactions. We have included three mechanisms to remove species from the grain mantles: thermal evaporation, cosmic-ray-induced heating, and photodesorption. A wide range of parameter space has been explored to examine the abundance of species present both on the grain mantles and in the gas phase as a function of both position in the cloud (visual extinction) and of evolutionary state (time). The dominant mechanism that removes molecules from the grain mantles is cosmic-ray desorption. At times greater than the depletion timescale, the abundances of some simple species agree with abundances observed in the cold dark cloud TMC-1. Even though cosmic-ray desorption preserves the gas-phase chemistry at late times, molecules do show significant depletions from the gas phase. Examination of the dependence of depletion as a function of density shows that when the density increases from 10(exp 3)/cc to 10(exp 5)/cc several species including HCO(+), HCN, and CN show gas-phase abundance reductions of over an order of magnitude. The CO: H2O ratio in the grain mantles for our standard model is on the order of 10:1, in reasonable agreement with observations of nonpolar CO ice features in rho Ophiuchus and Serpens. We have also examined the interdependence of CO depletion with the space density of molecular hydrogen and binding energy to the grain surface. We find that the observed depletion of CO in Taurus in inconsistent with CO bonding in an H2O rich mantle, in agreement with observations. We suggest that if interstellar grains consist of an outer layer of CO ice, then the binding energies for many species to the grain mantle may be lower than commonly used, and a significant portion of molecular material may be maintained in the gas phase.

Bergin, E. A.↗
NASA NTRS1987

Flexible system model reduction and control system design based upon actuator and sensor influence functions

A model reduction technique based on aggregation with respect to sensor and actuator influence functions rather than modes is presented for large systems of coupled second-order differential equations. Perturbation expressions which can predict the effects of spillover on both the reduced-order plant model and the neglected plant model are derived. For the special case of collocated actuators and sensors, these expressions lead to the derivation of constraints on the controller gains that are, given the validity of the perturbation technique, sufficient to guarantee the stability of the closed-loop system. A case study demonstrates the derivation of stabilizing controllers based on the present technique. The use of control and observation synthesis in modifying the dimension of the reduced-order plant model is also discussed. A numerical example is provided for illustration.

Yam, Yeung↗