Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduced Order Models”

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,405 records · Page 78

Learning Hidden Structure in Multi-Fidelity Information Sources for Efficient Uncertainty Quantification (LDRD 218317)

This report summarizes the work done under the Laboratory Directed Research and Development (LDRD) project entitled "Learning Hidden Structure in Multi-Fidelity Information Sources for Efficient Uncertainty Quantification". In this project we investigated multi-fidelity strategies for fusing data from information sources of varying cost and accuracy. Most existing strategies exploit hierarchical relationships between models, for example that occur when different models are generated by refining a numerical discretization parameter. In this work we focused on encoding the relationships between information sources using directed acyclic graphs. The multi-fidelity networks can have general structure and represent a significantly greater variety of modeling relationships than recursive networks used in the current state literature. Numerical results show that a non-hierarchical multi-fidelity Monte Carlo strategy can reduce the cost of estimating uncertainty in predictions of a model of plasma expanding in a vacuum by almost two orders of magnitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A tentative unified sea model for scattering and emission

This work is an attempt to relate active and passive observations of the sea surface to each other, considering that more information can be obtained from a set of instruments when they are not just considered as a juxtaposition of separate and independent boxes, but as a single complex instrument. A careful review of several approaches to the solution of the scattering problem for a random rough surface has been performed. Inherent simplifications, and inconsistencies, are pointed out. To reduce the computation time to reasonable values, an approximate model is selected and some numerical results will be presented in order to illustrate the benefit that can be gained by a unified approach.

Guissard, A.↗

Ganymede and Callisto - Surface textural dichotomies and photometric analysis

Complete solar phase curves of the Ganymede and Callisto leading and trailing hemispheres, which have been obtained by reducing Voyager imaging observations and combining them with ground-based telescopic data, are presently fit to scattering models in order to derive hemispherical values of the single scattering albedo, the single particle phase function (SPPF), the compaction state (CS) of the optically active portion of the regolith, and the mean slope angle of macroscopic features. While Callisto's leading side is composed of particles that are more strongly backscattering than the trailing side, no hemispheric differences are found in the CS, surface roughness, or SPPF.

Buratti, Bonnie J.↗

The AXAF/XRS ADR - Engineering model

A spaceworthy Adiabatic Demagnetization Refrigerator is under development at Goddard Space Flight Center as part of the X-Ray Spectrometer (XRS), an instrument on the Advanced X-ray Astrophysics Facility (AXAF). XRS will employ an array of 32 microcalorimeters capable of detecting X-rays in the energy range of 0.3 - 10 keV. In order to achieve a desired resolution of 12 eV, these detectors must be operated at a temperature of 0.065 - 0.100 K. An ADR must be used to cool these detectors in space. A breadboard model was designed and built less than two years ago, and provided excellent results. We are presently at the development stage of the engineering model. Several changes have been made to the original design in order to improve the efficiency of the ADR, to reduce its weight, and to strengthen the salt pill suspension system so that the ADR can survive launch loads and have low sensitivity to microphonic inputs. We shall report on the results of these changes; what worked and what did not.

Serlemitsos, Aristides T.↗

Numerical simulation of the non-Newtonian mixing layer

This work is a continuing effort to advance our understanding of the effects of polymer additives on the structures of the mixing layer. In anticipation of full nonlinear simulations of the non-Newtonian mixing layer, we examined in a first stage the linear stability of the non-Newtonian mixing layer. The results of this study show that, for a fluid described by the Oldroyd-B model, viscoelasticity reduces the instability of the inviscid mixing layer in a special limit where the ratio (We/Re) is of order 1 where We is the Weissenberg number, a measure of the elasticity of the flow, and Re is the Reynolds number. In the present study, we pursue this project with numerical simulations of the non-Newtonian mixing layer. Our primary objective is to determine the effects of viscoelasticity on the roll-up structure. We also examine the origin of the numerical instabilities usually encountered in the simulations of non-Newtonian fluids.

Azaiez, Jalel↗

Large-Eddy/Lattice Boltzmann Simulations of Micro-blowing Strategies for Subsonic and Supersonic Drag Control

This report summarizes the progress made in the first 8 to 9 months of this research. The Lattice Boltzmann Equation (LBE) methodology for Large-eddy Simulations (LES) of microblowing has been validated using a jet-in-crossflow test configuration. In this study, the flow intake is also simulated to allow the interaction to occur naturally. The Lattice Boltzmann Equation Large-eddy Simulations (LBELES) approach is capable of capturing not only the flow features associated with the flow, such as hairpin vortices and recirculation behind the jet, but also is able to show better agreement with experiments when compared to previous RANS predictions. The LBELES is shown to be computationally very efficient and therefore, a viable method for simulating the injection process. Two strategies have been developed to simulate multi-hole injection process as in the experiment. In order to allow natural interaction between the injected fluid and the primary stream, the flow intakes for all the holes have to be simulated. The LBE method is computationally efficient but is still 3D in nature and therefore, there may be some computational penalty. In order to study a large number or holes, a new 1D subgrid model has been developed that will simulate a reduced form of the Navier-Stokes equation in these holes.

Menon, Suresh↗

Foreign Object Damage Identification in Turbine Engines

This report summarizes the collective work of a five-person team from different organizations examining the problem of detecting foreign object damage (FOD) events in turbofan engines from gas path thermodynamic and bearing accelerometer sensors, and determining the severity of damage to each component (diagnosis). Several detection and diagnostic approaches were investigated and a software tool (FODID) was developed to assist researchers detect/diagnose FOD events. These approaches include (1) fan efficiency deviation computed from upstream and downstream temperature/ pressure measurements, (2) gas path weighted least squares estimation of component health parameter deficiencies, (3) Kalman filter estimation of component health parameters, and (4) use of structural vibration signal processing to detect both large and small FOD events. The last three of these approaches require a significant amount of computation in conjunction with a physics-based analytic model of the underlying phenomenon the NPSS thermodynamic cycle code for approaches 1 to 3 and the DyRoBeS reduced-order rotor dynamics code for approach 4. A potential application of the FODID software tool, in addition to its detection/diagnosis role, is using its sensitivity results to help identify the best types of sensors and their optimum locations within the gas path, and similarly for bearing accelerometers.

Strack, William↗

Modeling Increased Complexity and the Reliance on Automation: FLightdeck Automation Problems (FLAP) Model

This paper highlights the development of a model that is focused on the safety issue of increasing complexity and reliance on automation systems in transport category aircraft. Recent statistics show an increase in mishaps related to manual handling and automation errors due to pilot complacency and over-reliance on automation, loss of situational awareness, automation system failures and/or pilot deficiencies. Consequently, the aircraft can enter a state outside the flight envelope and/or air traffic safety margins which potentially can lead to loss-of-control (LOC), controlled-flight-into-terrain (CFIT), or runway excursion/confusion accidents, etc. The goal of this modeling effort is to provide NASA's Aviation Safety Program (AvSP) with a platform capable of assessing the impacts of AvSP technologies and products towards reducing the relative risk of automation related accidents and incidents. In order to do so, a generic framework, capable of mapping both latent and active causal factors leading to automation errors, is developed. Next, the framework is converted into a Bayesian Belief Network model and populated with data gathered from Subject Matter Experts (SMEs). With the insertion of technologies and products, the model provides individual and collective risk reduction acquired by technologies and methodologies developed within AvSP.

Ancel, Ersin↗

An Optimal Control Modification to Model-Reference Adaptive Control for Fast Adaptation

This paper presents a method that can achieve fast adaptation for a class of model-reference adaptive control. It is well-known that standard model-reference adaptive control exhibits high-gain control behaviors when a large adaptive gain is used to achieve fast adaptation in order to reduce tracking error rapidly. High gain control creates high-frequency oscillations that can excite unmodeled dynamics and can lead to instability. The fast adaptation approach is based on the minimization of the squares of the tracking error, which is formulated as an optimal control problem. The necessary condition of optimality is used to derive an adaptive law using the gradient method. This adaptive law is shown to result in uniform boundedness of the tracking error by means of the Lyapunov s direct method. Furthermore, this adaptive law allows a large adaptive gain to be used without causing undesired high-gain control effects. The method is shown to be more robust than standard model-reference adaptive control. Simulations demonstrate the effectiveness of the proposed method.

Nguyen, Nhan T.↗

Why is it still too warm or cold in my house? Examining the relationships between energy efficient capital and household energy insecurity

Here, this paper examines the relationships between energy efficient (EE) capital technology and household energy insecurity in the United States. The theoretical model of these relationships employs household production theory to capture the demand for and production of household energy services, and a stochastic production frontier approach to describe how having access to and the usage intensity of EE capital technology could help alleviate inefficiency in the production of household energy services. A working hypothesis formulated from the theoretical model posits that having EE capital technology in the home will reduce the level of household energy insecurity experienced. The extent of energy insecurity experienced is inferred from an energy insecurity index value assigned to each household, generated via the application of a dichotomous Rasch model to questions contained in the 2015 Residential Energy Consumption Survey. Noting the potential simultaneous relationship that exists between a household having access to and the usage intensity of EE capital technology and the experience of being energy insecure, an instrumental variables approach was employed to estimate a series of ordered logit models. Results suggest access to EE capital technology in the form of Energy Star® appliances, Energy Star® windows, or a SMART thermostat does not reduce the probability of experiencing a greater level of energy insecurity. Nor does the usage intensity of EE capital. Thus, policy instruments designed to alleviate household energy insecurity may need to go beyond simply helping households obtain EE capital technology.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Spaceborne Doppler Radar Measurements of Rainfall: Correction of Errors Induced by Pointing Uncertainties

In this paper a sea surface radar echo spectral analysis technique to correct for the rainfall velocity error caused by radar-pointing uncertainty is presented. The correction procedure is quite straightforward when the radar is observing a homogeneous rainfall field. When nonuniform beam filling (NUBF) occurs and attenuating frequencies are used, however, additional steps are necessary in order to correctly estimate the antenna-pointing direction. This new technique relies on the application of the combined frequency-time (CFT) algorithm to correct for uneven attenuation effects on the observed sea surface Doppler spectrum. The performance of this correction technique was evaluated by a Monte Carlo simulation of the Doppler precipitation radar backscatter from high-resolution 3D rain fields (either generated by a cloud resolving numerical model or retrieved from airborne radar measurements). The results show that the antenna-pointing-induced error can, indeed, be reduced by the proposed technique in order to achieve 1 m s(exp -1) accuracy on rainfall vertical velocity estimates.

remote sensing↗

Maxwell's second- and third-order equations of transfer for non-Maxwellian gases

Condensed algebraic forms for Maxwell's second- and third-order equations of transfer are developed for the case of molecules described by either elastic hard spheres, inverse-power potentials, or by Bird's variable hard-sphere model. These hardly reduced, yet exact, equations provide a new point of origin, when using the moment method, in seeking approximate solutions in the kinetic theory of gases for molecular models that are physically more realistic than that provided by the Maxwell model. An important by-product of the analysis when using these second- and third-order relations is that a clear mathematical connection develops between Bird's variable hard-sphere model and that for the inverse-power potential.

Baganoff, D.↗

On the magnetic structure of the quiet transition region

Existing models of the quiet chromosphere-corona transition region predict a distribution of emission measure over temperature that agrees with observation for a T greater than about 100,000 K. These 'network' models assume that all magnetic field lines that emerge from the photosphere extend into and are in thermal contact with the corona. It is shown that the observed fine-scale structure of the photospheric magnetic network instead suggests a two-component picture in which magnetic funnels that open into the corona emerge from only a fraction of the network. The gas that makes up the hotter transition region is mostly contained within these funnels, as in standard models, but, because the funnels are more constricted in our picture, the heat flowing into the cooler transition region from the corona is reduced by up to an order of magnitude. The remainder of the network is occupied by a population of low-lying loops with lengths less than about 10,000 km. It is proposed that the cooler transition region is mainly located within such loops, which are magnetically insulated from the corona and must, therefore, be heated internally. The fine-scale structure of ultraviolet spectroheliograms is consistent with this proposal, and theoretical models of internally heated loops can explain the behavior of the emission measure below a T of aobut 100,000 K.

Rabin, D.↗

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulating Bone Loss in Microgravity Using Mathematical Formulations of Bone Remodeling

Most mathematical models of bone remodeling are used to simulate a specific bone disease, by disrupting the steady state or balance in the normal remodeling process, and to simulate a therapeutic strategy. In this work, the ability of a mathematical model of bone remodeling to simulate bone loss as a function of time under the conditions of microgravity is investigated. The model is formed by combining a previously developed set of biochemical, cellular dynamics, and mechanical stimulus equations in the literature with two newly proposed equations; one governing the rate of change of the area of cortical bone tissue in a cross section of a cylindrical section of bone and one governing the rate of change of calcium in the bone fluid. The mechanical stimulus comes from a simple model of stress due to a compressive force on a cylindrical section of bone which can be reduced to zero to mimic the effects of skeletal unloading in microgravity. The complete set of equations formed is a system of first order ordinary differential equations. The results of selected simulations are displayed and discussed. Limitations and deficiencies of the model are also discussed as well as suggestions for further research.

Pennline, James A.↗

A geometric approach to regulator and tracker design for an aerospace plane

The paper presents a nonlinear design approach drawing from singular perturbations, feedback linearization, and variable structure control, that leads to regulators with automatic gain scheduling which exhibit similar dynamic behavior over the entire flight envelope of the aerospace plane. Additionally, design approach provides for a systematic way to counter disturbance effects as well as modeling uncertainties. The unifying feature of the three nonlinear feedback control methodologies is that they all have a geometric interpretation. First, the translational dynamics are decomposed into reduced-order slow and fast dynamics by way of a formal singular perturbation analysis. After feedback linearization the fast dynamics are robustly stabilized via a variable structure control approach. The slow dynamics are stabilized using conventional proportional-integral compensation based on the nominal slow dynamics. A number of sample command and disturbance responses at opposite ends of the flight envelope are presented for a nonlinear aerospace plane model.

Van Buren, Mark A.↗

Towards Enhancing Coding Productivity for GPU Programming Using Static Graphs

The main contribution of this work is to increase the coding productivity of GPU programming by using the concept of Static Graphs. GPU capabilities have been increasing significantly in terms of performance and memory capacity. However, there are still some problems in terms of scalability and limitations to the amount of work that a GPU can perform at a time. To minimize the overhead associated with the launch of GPU kernels, as well as to maximize the use of GPU capacity, we have combined the new CUDA Graph API with the CUDA programming model (including CUDA math libraries) and the OpenACC programming model. We use as test cases two different, well-known and widely used problems in HPC and AI: the Conjugate Gradient method and the Particle Swarm Optimization. In the first test case (Conjugate Gradient) we focus on the integration of Static Graphs with CUDA. In this case, we are able to significantly outperform the NVIDIA reference code, reaching an acceleration of up to 11x thanks to a better implementation, which can benefit from the new CUDA Graph capabilities. In the second test case (Particle Swarm Optimization), we complement the OpenACC functionality with the use of CUDA Graph, achieving again accelerations of up to one order of magnitude, with average speedups ranging from 2x to 4x, and performance very close to a reference and optimized CUDA code. Our main target is to achieve a higher coding productivity model for GPU programming by using Static Graphs, which provides, in a very transparent way, a better exploitation of the GPU capacity. The combination of using Static Graphs with two of the current most important GPU programming models (CUDA and OpenACC) is able to reduce considerably the execution time w.r.t. the use of CUDA and OpenACC only, achieving accelerations of up to more than one order of magnitude. Finally, we propose an interface to incorporate the concept of Static Graphs into the OpenACC Specifications.

58 GEOSCIENCES↗

In-Situ Scanning Electron Microscope Experiments for Microscale Mechanical Testing and Validated Modeling of Fiber Reinforced Thermoplastics

A novel, in-situ, scanning electron microscope (SEM) mechanical testing capability for materials at the microscale which provides experimental validation to a machine learning (ML) toolset for full-field validation of physics-based micromechanics models is being developed by researchers at NASA Glenn Research Center. These are enabling technologies for the integration of multiscale digital twins for materials into system level models which will result in the improved performance, material discovery, reduced production cost and time, rapid characterization, and prognostic structural health monitoring (SHM) for materials and structures for extreme environments in support of NASA space exploration missions. In order to bridge the material structure-to-system gap for digital twins, physics-based models must be experimentally validated at multiple length scales. Seminal microscale experiments, conducted at the Air Force Research Laboratory (AFRL), were limited to transverse compression of single-layer, unidirectional thermoset polymer matrix composite (PMC) micropillar specimens [1]. The early phases of the current project followed those initial results and setup to reproduce the compression testing of PMC material on the custom-built piezoelectric actuated micromechanical testing rig built by MicroTesting Solutions LLC. In this work, samples of thermoplastic PMC material were first machined into 3 mm cubes, and then further machining and final milling was done using a Focused Ion Beam (FIB). The initial experiment was done on a pillar roughly 20 µm x 20 µm x 40 µm tall. Additional pillars were milled with final sizes ranging from 20 µm x 20 µm x 40 µm tall to 40 µm x 40 µm x 65 µm tall. A speckle pattern for in-situ full-field measurements using Digital Image Correlation (DIC) was applied with platinum, which was coated on the surface, and then the FIB was used to mill away some of the coating to produce an irregular pattern of Pt on the pillar surface. The samples were loaded into the custom testing rig and placed into the SEM and loaded under compression until failure. Images were collected in the SEM during testing. Post-processing of the images was conducted using DIC to obtain full-field displacement and strain measurements elucidating the role of the matrix as well as fiber-fiber interaction at the microscale within the composite subjected to compression loading well into the non-linear regime of the material. Moreover, the evolution of fiber-matrix debonding and matrix cracking is observed in-situ at the microscale. This data, along with images segmented with a newly developed ML toolset [2], was used to create and validate physics-based micromechanics models. An image of the failed micropillar is shown in Figure 1. The techniques developed in the initial compression experiment was tailored to the validation needs of the models and expanded to include different sized samples as well as possibly tension and fatigue.

Laura Wilson↗