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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Surface Roughness Effects from Additive Manufacturing in High Efficiency Gas Turbine Combustion Systems

In the past, clean combustion systems were characterized by high system development costs due to growing complexity and escalating manufacturing costs from conventional manufacturing processes. As a result, high efficiency concepts are difficult to design and more difficult to be cost-effectively manufactured. In recent years, with the introduction of additive manufacturing (AM) technologies, the rapid prototyping and mass production processes of clean combustion systems are promising to be significantly simplified with significant reduction in terms of engine manufacturing cost and engine energy cost. Compared with conventional manufacturing built parts, the AM process enabled a simpler design to be adopted for the nozzle, reducing the number of required braze and weld joints from twenty-five to just five. The resulting nozzle was 25% lighter and five times more durable and contributed to a 15% reduction in fuel burn in comparison with the previous model produced. Applying these improvements to a fleet of 5,000 turbofan engines at 150 kN thrust would result in fuel cost savings of $6B annually (at $5/gallon Jet-A fuel price, 700 gallon/hour fuel consumption rate, and 2,300 operational hours per year),and reduce CO2 emissions by more than 25 million tons per year. These advances support core Department of Energy (DOE)missions in energy efficiency, improving productivity, and environmental sustainability. Maximizing the benefit of these new capabilities will require high prediction capability of high speed turbulent flow with wall modeled Large Eddy Simulation (LES). GE Aviation maintains that advanced simulation technology and supercomputing is required in order to provide the appropriate boundary conditions to realize the maximum potential of AM and Ceramic Matrix Composite to reduce cooling flow, a first order penalty on the Brayton cycle. The algorithm developed using the ANSYS/Fluent software would provide the foundation for entirely new avenues of research and development with potential multi-billion dollar impact to the US economy, and significant reduction in carbon based emissions across the aerospace and power generation industries. The formulation of this iWLES (integral wall model for LES) is generic and allows to capture the changes in flow dynamics that have a significant impact on the wall bounded flow characteristics, such as swirler effective area, bulk swirl number, local pressure distribution, exit velocity profile, and turbulence kinetic energy profile. A periodic channel flow with rough flat plate is simulated using LES (Wall-Adapting Local Eddy-viscosity) model to verify the implementation of iWLES in the Fluent User Defined Function (UDF). As the flow fields are highly sensitive due to surface roughness of the wall bounded flows in the combustion systems, there is significant potential to conduct further LES studies focusing on the turbulence boundary layer interaction. Accurately capturing near wall physics will elucidate the impact of rough surfaces on combustor flow and aero-thermal interactions.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS1982

Synthesis of a modern/classical controller for large space structures

The paper is concerned with wide-sense stationary random vibrations of a representative large space structure subject to nonwhite excitations, and with synthesis of its controller amid modeling uncertainties due to truncated modes. Dynamics is analyzed by residue calculus in frequency domain. Steady state variance of a modal coordinate is shown to depend on, among other things, power spectral density of excitation at the modal frequency. Modal cost analysis is performed for an order reduction with a performance objective of attitude, attitude rate and energy control. Optimal and suboptimal output feedback controllers are compared. By treating deleted modes as an additive and a multiplicative perturbation in an idealized transfer function, the low frequency modes omitted according to modal cost analysis are found to cause instability. Finally, the size of the idealized model is shown to be governed by light damping and stringency in performance and stability specifications.

Hablani, H. B.↗
DOE OSTI2024

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Furthermore, adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Multiresolutional models of uncertainty generation and reduction

Kolmogorov's axiomatic principles of the probability theory, are reconsidered in the scope of their applicability to the processes of knowledge acquisition and interpretation. The model of uncertainty generation is modified in order to reflect the reality of engineering problems, particularly in the area of intelligent control. This model implies algorithms of learning which are organized in three groups which reflect the degree of conceptualization of the knowledge the system is dealing with. It is essential that these algorithms are motivated by and consistent with the multiresolutional model of knowledge representation which is reflected in the structure of models and the algorithms of learning.

Meystel, A.↗
DOE OSTI2021

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗
NASA NTRS1973

Dynamics of Flexible Spinning Satellites with Radial Wire Antennas

A dynamic analysis is presented for a spin stabilized spacecraft employing four radial wire antennas with tip masses, a configuration first employed in the IMP-J spacecraft. The use of wires in place of the usual booms represents the ultimate in weight reduction at the expanse of flexibility. The satellite is modelled as a 14 degree of freedom system, and the linearized equations of motion are found. The lowest order vibrational modes and natural frequencies of the gyroscopically coupled system are then determined. Because the satellite spin rate is decreased by antenna deployment, a spin-up maneuver is needed. The response of the time varying mode equations during spin-up is found, for the planar modes, in terms of Bessel functions and a Struve function of order -1/4. Because tables of the latter are not readily available, the particular solution is expressed in various forms including an infinite series of Bessel functions and a particularly useful asymptotic expansion.

Longman, R. W.↗
NASA NTRS1980

A method of order reduction for structural dynamics

This paper presents a new method of order reduction based on recent work of a similar nature applicable to system dynamics and control. This method provides a practical computational procedure for producing a condensed model which exactly preserves the slowest n1 modes (where n1 is the number of 'slow' eigenvalues) of the total n modes of the original model for almost any set of n1 degrees of freedom retained in the condensed model. The method can also be used to compute the eigensolutions corresponding to the n1 slowest modes of the original structural dynamics problem. The accuracy of the condensed model and the speed/accuracy performance of the eigensolver are compared with standard methods for a 90 DOF cantilevered plate.

Anderson, L. R.↗
NASA NTRS1992

Theoretical study of thermodynamic properties and reaction rates of importance in the high-speed research program

One of the primary goals of NASA's high-speed research program is to determine the feasibility of designing an environmentally safe commercial supersonic transport airplane. The largest environmental concern is focused on the amount of ozone destroying nitrogen oxides (NO(x)) that would be injected into the lower stratosphere during the cruise portion of the flight. The limitations placed on NO(x) emission require more than an order of magnitude reduction over current engine designs. To develop strategies to meet this goal requires first gaining a fundamental understanding of the combustion chemistry. To accurately model the combustor requires a computational fluid dynamics approach that includes both turbulence and chemistry. Since many of the important chemical processes in this regime involve highly reactive radicals, an experimental determination of the required thermodynamic data and rate constants is often very difficult. Unlike experimental approaches, theoretical methods are as applicable to highly reactive species as stable ones. Also our approximation of treating the dynamics classically becomes more accurate with increasing temperature. In this article we review recent progress in generating thermodynamic properties and rate constants that are required to understand NO(x) formation in the combustion process. We also describe our one-dimensional modeling efforts to validate an NH3 combustion reaction mechanism. We have been working in collaboration with researchers at LeRC, to ensure that our theoretical work is focused on the most important thermodynamic quantities and rate constants required in the chemical data base.

Langhoff, Stephen↗
NASA NTRS1973

Noise reduction of a tilt-rotor aircraft including effects on weight and performance

Various methods for far-field noise reduction of a tilt-rotor acoustic signature and the performance and weight tradeoffs which result from modification of the noise sources are considered in this report. In order to provide a realistic approach for the investigation, the Boeing tilt-rotor flight research aircraft (Model 222), was selected as the baseline. This aircraft has undergone considerable engineering development. Its rotor has been manufactured and tested in the Ames full-scale wind tunnel. The study reflects the current state-of-the-art of aircraft design for far-field acoustic signature reduction and is not based solely on an engineering feasibility aircraft. This report supplements a previous study investigating reduction of noise signature through the management of the terminal flight trajectory.

Gibs, J.↗
DOE OSTI2024

A Reduced-Order Model of a Nuclear Power Plant with Thermal Power Dispatch

This paper presents reduced-order modeling of thermal power dispatch (TPD) from a pressurized water reactor (PWR) for providing heat to nearby heat consuming industrial processes that seek to take advantage of nuclear heat to reduce carbon emissions. The reactor model includes the neutronics of the reactor core, thermal–hydraulics of the primary coolant cycle, and a three-lump model of the steam generator (SG). The secondary coolant cycle is represented with quasi-steady state mass and energy balance equations. The secondary cycle consists of a steam extraction system, high-pressure and low-pressure turbines, moisture separator and reheater, high-pressure and low-pressure feedwater heaters, deaerator, feedwater and condensate pumps, and a condenser. The steam produced by the SG is distributed between the turbines and the extraction steam line (XSL) that delivers steam to nearby industrial processes, such as production of clean hydrogen. The reduced-order simulator is verified by comparing predictions with results from separate validated steady-state and transient full-scope PWR simulators for TPD levels between 0% and 70% of the rated reactor power. All simulators indicate that the flow rate of steam in the main steam line and turbine systems decrease with increasing TPD, which causes a reduction in PWR electric power generation. The results are analyzed to assess the impact of TPD on system efficiency and feedwater flow control. Due to the simplicity of the proposed reduced-order model, it can be scaled to represent a PWR of any size with a few parametric changes. In the future, the proposed reduced-order model will be integrated into a power system model in a digital real-time simulator (DRTS) and physical hardware-in-the-loop simulations.

08 HYDROGEN↗
NASA NTRS1985

A method for high order linear system reduction and nonlinear system simplification

Least-squares-type algorithms for reducing the order of linear systems in the frequency domain and simplifying nonlinear systems in time domain are developed and demonstrated. The possible model structures are represented as nodes in a tree, and costs along the branches are assigned using the repeated-Gram-Schmidt orthogonalization procedure of Desrochers and Saridis (1980), permitting identification of the optimal n-term model by searching the tree to depth n, with no need for parameter identification. The efficiency and flexibility of the algorithms is shown in applications to the eighth-order linear system studied by Hsia (1972), a three-state eight-nonlinear-term aircraft-dynamics problem, and the related linear-controller problem (Garrard and Jordan, 1977).

Desrochers, A. A.↗
DOE OSTI2021

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗
DOE OSTI2020

Chromium isotope fractionation during reduction of Chromium(VI) by Iron(II/III)-bearing clay minerals

Chromium stable isotope ratios are used to trace the reduction of Cr(VI) to Cr(III) in both ancient and modern systems. However, quantitative interpretation of Cr isotopic signatures has been stymied by the large variability in isotopic fractionation factors for Cr(VI) reduction by different reductants. Here we determine Cr isotope fractionation factors during Cr(VI) reduction by Fe(II/III)-bearing clay minerals, which are abundant in subsurface environments. Several variables were tested: pH, total Fe content of the clay, and the fraction of reduced Fe within the clay (Fe(II)/Fe(total)). The latter controls the standard reduction potential of the clay. Our results demonstrate that neither pH nor total Fe content of the clay have major effects on isotopic fractionation. In contrast, as the effective standard reduction potential of the clay and thus the standard free energy of Cr(VI) reduction become more negative, Cr isotope fractionation factors decrease in magnitude from - 4.9 to -1.3‰ according to a linear free energy relationship. This linear free energy relationship can be predicted from Marcus electron transfer theory and allows first-order predictions of Cr isotope fractionation factors to be made from the standard reduction potential or Fe(II)/Fe(total) of a clay, potentially improving our ability to model Cr isotope signatures in geochemical systems. Chromium is the first isotope system to show such a linear free energy relationship over a diverse range of reductants, including both aqueous and solid-phase reductants, and may provide a model for determining other redox-driven kinetic isotope effects in environmentally important isotope systems.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Lattice Boltzmann model for conjugate heat transfer across thin walls

A lattice Boltzmann (LB) model with an efficient and accurate interface treatment for conjugate heat transfer across a thin wall between two different media is developed. The proposed interface treatment avoids fine meshing and computation within the thin layer; instead, the energy balance within the thin layer and the conjugate conditions on each interface are utilized to construct explicit updating schemes for the microscopic distribution functions of the LB model at the interior lattice nodes of the two media next to the thin layer. The proposed interface schemes reduce to the standard interface scheme for conjugate conditions in the literature in the limit of zero thickness of the thin layer, and thus it can be considered a more general interface treatment. A simplified version of the interface treatment is also proposed when the heat flux variation along the tangential direction of the thin layer is negligible. Three representative numerical tests are conducted to verify the applicability and accuracy of the proposed interface schemes. The results demonstrate that the intrinsic second-order accuracy of the LB model is preserved with the proposed interface schemes for thin layers with constant tangential fluxes, while for general situations with varying tangential fluxes, first-order accuracy is obtained. Finally, this interface treatment within the LB framework is attractive in conjugate heat transfer modeling involving thin layers for its simplicity, accuracy, and significant reduction in computational resources.

42 ENGINEERING↗
DOE OSTI2021

Space-Time Reduced-Order Modeling for Uncertainty Quantification

This work focuses on the space-time reduced-order modeling (ROM) method for solving large-scale uncertainty quantification (UQ) problems with multiple random coefficients. In contrast with the traditional space ROM approach, which performs dimension reduction in the spatial dimension, the space-time ROM approach performs dimension reduction on both the spatial and temporal domains, and thus enables accurate approximate solutions at a low cost. We incorporate the space-time ROM strategy with various classical stochastic UQ propagation methods such as stochastic Galerkin and Monte Carlo. Numerical results demonstrate that our methodology has significant computational advantages compared to state-of-the-art ROM approaches. By testing the approximation errors, we show that there is no obvious loss of simulation accuracy for space-time ROM given its high computational efficiency.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

An Eigensystem Realization Algorithm (ERA) for modal parameter identification and model reduction

A method, called the Eigensystem Realization Algorithm (ERA), is developed for modal parameter identification and model reduction of dynamic systems from test data. A new approach is introduced in conjunction with the singular value decomposition technique to derive the basic formulation of minimum order realization which is an extended version of the Ho-Kalman algorithm. The basic formulation is then transformed into modal space for modal parameter identification. Two accuracy indicators are developed to quantitatively identify the system modes and noise modes. For illustration of the algorithm, examples are shown using simulation data and experimental data for a rectangular grid structure.

Juang, J. N.↗
NASA NTRS2000

Effect of Negative Ions on the Conductivity of the Titan Atmosphere

In an earlier paper, Borucki et al (1987) calculated the electrical conductivity and electrical charge on aerosols in Titan's atmosphere due to the ionization by galactic cosmic rays and electron precipitation from Saturn's magnetosphere. The lower atmosphere was predicted to be substantially more conducting than the atmospheres of Earth and Venus because of the high concentration of free electrons. The prediction of a high conductivity is based on the lack of electrophillic species which form negative ions with low mobility and which reduce the number of free electrons. At that time, no molecular species capable of forming negative ions in concentrations sufficient to perturb the atmospheric conductivity were identified. Recently, E. Bakes and her colleagues have been investigating the formation of nitrogenous macromolecules using quantum mechanical methods. Their calculations indicate that the molecules will be highly electrophillic and are likely to be present in the atmosphere at mixing ratios of order 10(exp -7). This mixing ratio is sufficiently large that a substantial reduction in the conductivity is expected at altitudes below 100 km. Revision of the atmospheric model to accommodate the presence of negative ions and to increase the fidelity of the modeling will be described.

Borucki, W. J.↗
NASA NTRS2007

Model Checking a Byzantine-Fault-Tolerant Self-Stabilizing Protocol for Distributed Clock Synchronization Systems

This report presents the mechanical verification of a simplified model of a rapid Byzantine-fault-tolerant self-stabilizing protocol for distributed clock synchronization systems. This protocol does not rely on any assumptions about the initial state of the system. This protocol tolerates bursts of transient failures, and deterministically converges within a time bound that is a linear function of the self-stabilization period. A simplified model of the protocol is verified using the Symbolic Model Verifier (SMV) [SMV]. The system under study consists of 4 nodes, where at most one of the nodes is assumed to be Byzantine faulty. The model checking effort is focused on verifying correctness of the simplified model of the protocol in the presence of a permanent Byzantine fault as well as confirmation of claims of determinism and linear convergence with respect to the self-stabilization period. Although model checking results of the simplified model of the protocol confirm the theoretical predictions, these results do not necessarily confirm that the protocol solves the general case of this problem. Modeling challenges of the protocol and the system are addressed. A number of abstractions are utilized in order to reduce the state space. Also, additional innovative state space reduction techniques are introduced that can be used in future verification efforts applied to this and other protocols.

Malekpour, Mahyar R.↗