Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state process”

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 487 records · Page 27

Constructing an Efficient Self-Tuning Aircraft Engine Model for Control and Health Management Applications

Self-tuning aircraft engine models can be applied for control and health management applications. The self-tuning feature of these models minimizes the mismatch between any given engine and the underlying engineering model describing an engine family. This paper provides details of the construction of a self-tuning engine model centered on a piecewise linear Kalman filter design. Starting from a nonlinear transient aerothermal model, a piecewise linear representation is first extracted. The linearization procedure creates a database of trim vectors and state-space matrices that are subsequently scheduled for interpolation based on engine operating point. A series of steady-state Kalman gains can next be constructed from a reduced-order form of the piecewise linear model. Reduction of the piecewise linear model to an observable dimension with respect to available sensed engine measurements can be achieved using either a subset or an optimal linear combination of "health" parameters, which describe engine performance. The resulting piecewise linear Kalman filter is then implemented for faster-than-real-time processing of sensed engine measurements, generating outputs appropriate for trending engine performance, estimating both measured and unmeasured parameters for control purposes, and performing on-board gas-path fault diagnostics. Computational efficiency is achieved by designing multidimensional interpolation algorithms that exploit the shared scheduling of multiple trim vectors and system matrices. An example application illustrates the accuracy of a self-tuning piecewise linear Kalman filter model when applied to a nonlinear turbofan engine simulation. Additional discussions focus on the issue of transient response accuracy and the advantages of a piecewise linear Kalman filter in the context of validation and verification. The techniques described provide a framework for constructing efficient self-tuning aircraft engine models from complex nonlinear simulations.Self-tuning aircraft engine models can be applied for control and health management applications. The self-tuning feature of these models minimizes the mismatch between any given engine and the underlying engineering model describing an engine family. This paper provides details of the construction of a self-tuning engine model centered on a piecewise linear Kalman filter design. Starting from a nonlinear transient aerothermal model, a piecewise linear representation is first extracted. The linearization procedure creates a database of trim vectors and state-space matrices that are subsequently scheduled for interpolation based on engine operating point. A series of steady-state Kalman gains can next be constructed from a reduced-order form of the piecewise linear model. Reduction of the piecewise linear model to an observable dimension with respect to available sensed engine measurements can be achieved using either a subset or an optimal linear combination of "health" parameters, which describe engine performance. The resulting piecewise linear Kalman filter is then implemented for faster-than-real-time processing of sensed engine measurements, generating outputs appropriate for trending engine performance, estimating both measured and unmeasured parameters for control purposes, and performing on-board gas-path fault diagnostics. Computational efficiency is achieved by designing multidimensional interpolation algorithms that exploit the shared scheduling of multiple trim vectors and system matrices. An example application illustrates the accuracy of a self-tuning piecewise linear Kalman filter model when applied to a nonlinear turbofan engine simulation. Additional discussions focus on the issue of transient response accuracy and the advantages of a piecewise linear Kalman filter in the context of validation and verification. The techniques described provide a framework for constructing efficient self-tuning aircraft engine models from complex nonlinear simulatns.

Armstrong, Jeffrey B.↗

ALTKAL: An optimum linear filter for GEOS-3 altimeter data

ALTKAL is a computer program designed to smooth sea surface height data obtained from the GEOS 3 altimeter, and to produce minimum variance estimates of sea surface height and sea surface slopes, along with their standard derivations. The program operates by processing the data through a Kalman filter in both the forward and backward directions, and optimally combining the results. The sea surface height signal is considered to have a geoid signal, modeled by a third order Gauss-Markov process, corrupted by additive white noise. The governing parameters for the signal and noise processes are the signal correlation length and the signal-to-noise ratio. Mathematical derivations of the filtering and smoothing algorithms are presented. The smoother characteristics are illustrated by giving the frequency response, the data weighting sequence and the transfer function of a realistic steady-state smoother example. Based on nominal estimates for geoidal undulation amplitude and correlation length, standard deviations for the estimated sea surface height and slope are 12 cm and 3 arc seconds, respectively.

Fang, B. T.↗

A theory of stationary long waves. III - Quasi-normal modes in a singular waveguide

The present paper deals with the fundamental issue of whether one can treat waves as normal modes when critical surface, where the phase speed of the wave matches the zonal wind speed, are present. In particular the question of whether a Rossby critical level (such as the zero-wind line for stationary waves) is absorbing or reflecting is raised and subsequently addressed. It is found that the critical level is never totally absorbing; Rossby waves are partially reflected even if the critical layer is dominated by dissipative processes. The relevance of nonlinearity in planetary-scale Rossby wave critical layers is also discussed and it is found to be the dominant mechanism. With the relative magnitudes of nonlinearity versus viscosity relevant to the earth's atmosphere it is found that the steady-state critical level should be almost perfectly reflecting to incident Rossby waves. Consequently, normal-mode solutions can be found; the quantization condition for these waves is also derived.

Tung, K. K.↗

Mephisto: Interfacial Destabilization in Metal Alloys

The destabilizing mechanisms at a solidification interface were studied to obtain information on the kinetics and morphologies in the transient and steady state, and to separate the influences of liquid phase instabilities from interfacial instabilities. A differential seebeck voltage measurements technique was developed to provide a continuous record of the solid-liquid interface temperature as the solidification rate is varied to determine the kinetic coefficients. Signal processing and noise suppression techniques allow nonovolt precision which corresponds to mK accuracy for the interfacial temperature.

Favier, J. J.↗

Dispersal of Disks Around Young Stars: Constraints on Kuiper Belt Formation

We review the evidence pertaining to the lifetimes of planet-forming disks and discuss possible disk dispersal mechanisms: 1) viscous accretion of material onto the central source; 2) close stellar encounters; 3) stellar winds; and 4) photoevaporation by ultraviolet radiation. We focus on 3) and 4) and describe the quasi-steady state appearance and the overall evolution of disks under the influence of winds and radiation from the central star and of radiation from external OB stars. Viscous accretion likely dominates disk dispersal in the inner disk (r approx. <= 10 AU), while photoevaporation is the principal process of disk dispersal outside of r approx. >= 10 AU for low mass stars. Disk dispersal timescales are compared and discussed in relation to theoretical estimates for planet formation timescales. Photoevaporation may explain the large differences in the hydrogen content of the giant planets in the solar system. The commonly held belief that our early sun's stellar wind dispersed the solar nebula is called into question. Finally, we study the constraints that the evaporation of the outer disk has on the formation of Kuiper belts in extrasolar planetary systems.

Hollenbach, David↗

Evolution of the Electron Yield Curves of Insulators as a Function of Impinging Electron Fluence and Energy

Electron emission and concomitant change accumulation near the surface of insulators in central to understanding spacecraft charging. We present a study of changes in electron emission yields as a result of internal charge build up due to electron dose. Evolution of total, backscattered and secondary yield results over a broad range of incident energies are presented for two representative insulators, Kapton(trademark) and Al2O3. Reliable yield curves for uncharged insulators are measured and quantifiable changes in yield are observed due to < 100 fC/sq mm fluences. We find excellent agreement with a phenomenological argument based on insulator charging predicted by the yield curve; this includes a decrease in the rate of change of the yield as incident energies approach the crossover energies and as accumulated internal charge reduces the landing energy to asymptotically approach a steady state surface charge and unity yield. We also find that the exponential decay of yield curves with fluence exhibit an energy dependent decay constant, alpha(E). Finally, we discuss physics based models for this energy dependence. To understand fluence and energy dependence of these charging processes requires knowledge of how charge is deposited within the insulator, the mechanisms for charge trapping and transport within the insulator, and how the profile of trapped charge affects the transport and emission of charges from insulators.

Sim, Alec↗

Promoted Metals Combustion at Ambient and Elevated Temperatures

Promoted combustion testing of materials, Test 17 of NASA STD-6001, has been used to assess metal propensity to burn in oxygen rich environments. An igniter is used at the bottom end of a rod to promote ignition, and if combustion is sustained, the burning progresses from the bottom to the top of the rod. The physical mechanisms are very similar to the upward flammability test, Test 1 of NASA STD-6001. The differences are in the normal environmental range of pressures, oxygen content, and sample geometry. Upward flammability testing of organic materials can exhibit a significant transitional region between no burning to complete quasi-state burning. In this transitional region, the burn process exhibits a probabilistic nature. This transitional region has been identified for metals using the promoted combustion testing method at ambient initial temperatures. The work given here is focused on examining the transitional region and the quasi-steady burning region both at conventional ambient testing conditions and at elevated temperatures. A new heated promoted combustion facility and equipment at Marshall Space Flight Center have just been completed to provide the basic data regarding the metals operating temperature limits in contact with oxygen rich atmospheres at high pressures. Initial data have been obtained for Stainless Steel 304L, Stainless Steel 321, Haynes 214, and Inconel 718 at elevated temperatures in 100-percent oxygen atmospheres. These data along with an extended data set at ambient initial temperature test conditions are examined. The pressure boundaries of acceptable, non-burning usage is found to be lowered at elevated temperature.

Engel, Carl D.↗

Bubble Dynamics on a Heated Surface

In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Kassemi, M.↗

Control-augmented structural synthesis

A control-augmented structural synthesis methodology is presented in which the structural member sizes and active control system feedback gains are treated simultaneously as independent design variables. Multiple static and harmonic dynamic loading conditions are considered. Constraints are imposed on static displacements, natural frequencies, and the magnitudes of the steady-state dynamic displacements and actuator forces to ensure acceptable system behavior. Side constraints imposed on the design variables protect against the generation of unrealizable designs. Example problems are presented that demonstrate the method and underscore the importance of integrating the structural and active control system design process.

Lust, R. V.↗

Diamond and diamondlike carbon as wear-resistant, self-lubricating coatings for silicon nitride

Recent work on the friction and wear properties of as-deposited fine-grain diamond, polished coarse-grain diamond, and as-deposited diamondlike carbon (DLC) films in humid air at a relative humidity of approximately 40 percent and in dry nitrogen is reviewed. Two types of chemical vapor deposition (CVD) processes are used to deposit diamond films on silicon nitride (Si3N4) substrates: microwave-plasma and hot-filament. Ion beams are used to deposit DLC films of Si3N4 substrates. The diamond and DLC films in sliding contact with hemispherical bare Si3N4 pins have low steady-state coefficients of friction (less than 0.2) and low wear rates (less than 10(exp -7) mm(exp 2)/N-m), and thus, can be used effectively as wear-resistant, self-lubricating coatings for Si3N4 in the aforementioned two environments.

Miyoshi, Kazuhisa↗

Phase Transition in Dusty Plasmas: A Microphysical Description

Dust grains immersed in plasma discharges acquire a large negative charge and settle into a dust cloud at the edge of the sheath. In this region, the plasma ions stream toward the electrode at a velocity u approx. cs=(T(sub e)/m(sub i))(exp 1/2). Experimentally at sufficiently high gas pressure P, the random kinetic energy of the grains is damped by gas friction, and the grains are strongly coupled and self-organize into a crystalline configuration. For lower pressures despite the dissipation of grain kinetic energy to gas friction, the dust grains reach a steady-state kinetic temperature T(sub d) which is much larger than the temperature of any other component in the plasma. T(sub d) is so large that the dust acts like a fluid. We have used the dynamically shielded dust (DSD) model to simulate these physical processes. We find that the known experimental features are nicely reproduced in the simulations, and that additional features are revealed. In the figure we plot the variation of T(sub d) as P is continuously varied in a DSD code run. A marked difference is evident between the critical pressure P(sub m) for the melting transition as P is decreased, and the critical pressure P(sub c) for the condensation transition as P is increased. For P(sub m) is less than P is less than P(sub c), mixed phase states are seen. This hysteresis occurs because the instability which triggers melting is different from the instability that heats the dust in the fluid phase and inhibits freezing. At low pressure, the dust is subject to a two-stream instability with the ions. This instability is responsible for the high temperature of the dust at low pressure. The basic physics underlying the melting transition has been elucidated in a series of papers. We are developing a first-principles analytic approach to the melting transition, which embodies the same physics that is present in the DSD code.

Joyce, Glenn↗

Data and Modeling of Dendrites Subject to A Step Change in Pressure (TDSE)

There is considerable interest in dendritic solidification because of the influence dendrites have in the determination of microstructure, and thereby in the physical properties of cast metals and alloys. Current theories and models of dendritic growth generally couple diffusion effects in the melt with the physics of the interface. Data and subsequent analysis prior of the tip growth speed and radii of thermal succinonitrile dendrites in the near-convection free, on-orbit, free-fall environment demonstrate that these theories yield predictions that are reasonably in agreement with the results of experiment. However, data and analysis for assessing the interfacial physics component of theory are not sufficiently detailed or definitive. To study fundamental aspects of dendritic interface stability, we are measuring and modeling the kinetics and morphology of dendrites as they evolve from one well-defined steady state at a pre-set supercooling, through a transient stage, to a different well-defined steady state. More specifically, we subject succinonitrile dendrites, growing under steady-state conditions, to a rapid change in pressure. This leads to a rapid change in thermal driving force from the corresponding change in both the equilibrium melting temperature due to the Clapeyron effect, and a change in the far-field temperature due to adiabatic temperature changes in the bulk liquid and solid. Subsequently, we observe transformations from a well-characterized initial state into a new steady-state. Initial data reveal that the dendrite tip velocity changes almost as fast as the pressure charges, while the tip radius changes occur more slowly, taking from 10 60 seconds depending on the size of the step change and the final supercooling. Computer modeling of this process shows both agreements and disagreements with the experimental data. In making these observations and measurements, we are gaining new understandings of interfacial dynamics and state-selection physics.

Koss, Matthew B.↗

An Overview of the Patch Integral Method (PIM), a New Heat Transfer Analysis Tool for Hypersonic Wind Tunnel Facilities at NASA Langley

NASA Langley’s hypersonic wind tunnels are heavily leveraged for planetary missions. The data collection method in these tunnels is thermography, and surface temperature measurements of the model surface are collected and reduced to produce surface heating data, as seen in Fig 1. However, during model injection, no temperature data are collected, and thus conventional, integral heat transfer methods cannot be used to solve for surface heating. A method was developed in the 1990’s to reduce this data despite the data gap, known as the step approximation method. The method assumes that the film coefficient behaves as a step function, the model is semi-infinite, and thermal properties are constant. With these simplifying assumptions, a Laplace transform can be performed to result in an equation that takes the initial temperature of the model and a temperature at some point in time to back out the film coefficient at that time. This is the method that is used in the current thermographic data reduction software, IHEAT. While computationally light-weight, the step approximation has several issues associated with it. The time-history of temperature is not accounted for, which is vital as heat transfer is an integral process. Additionally, the required semi-infinite assumption is unnecessary and might be violated during runtime. Thermal variation of material properties can have a sizeable impact on heating results and are not modeled by the method. This method also takes multiple seconds to “collapse” to a steady state, which is undesirable from both a facility and data reduction standpoint. The method is also very sensitive to the “effective time” approximation, an approximation of when heating instantaneously starts (which is a nonphysical simplification), and a small variation in this value can result in an error in heating results.

J. S. Cheatwood↗

Quantifying the Energy Budget in the Solar Wind from 13.3 to 100 Solar Radii

A variety of energy sources, ranging from dynamic processes, such as magnetic reconnection and waves, to quasi steady terms, such as plasma pressure, may contribute to the acceleration of the solar wind. We utilize a combination of charged particle and magnetic field observations from the Parker Solar Probe (PSP) to attempt to quantify the steady-state contribution of the proton pressure, the electric potential, and the wave energy to the solar wind proton acceleration observed by PSP between 13.3 and ∼100 solar radii (R☉). The proton pressure provides a natural kinematic driver of the outflow. The ambipolar electric potential acts to couple the electron pressure to the protons, providing another definite proton acceleration term. Fluctuations and waves, while inherently dynamic, can act as an additional effective steady-state pressure term. To analyze the contributions of these terms, we utilize radial binning of single-point PSP measurements, as well as repeated crossings of the same stream at different distances on individual PSP orbits (i.e., fast radial scans). In agreement with previous work, we find that the electric potential contains sufficient energy to fully explain the acceleration of the slower wind streams. On the other hand, we find that the wave pressure plays an increasingly important role in the faster wind streams. The combination of these terms can explain the continuing acceleration of both slow and fast wind streams beyond 13.3 R☉.

J S Halekas↗

G-189A analytical simulation of the integrated waste management-water system using radioisotopes for thermal energy

An analytical simulation of the RITE-Integrated Waste Management and Water Recovery System using radioisotopes for thermal energy was prepared for the NASA-Manned Space Flight Center (MSFC). The RITE system is the most advanced concept water-waste management system currently under development and has undergone extended duration testing. It has the capability of disposing of nearly all spacecraft wastes including feces and trash and of recovering water from usual waste water sources: urine, condensate, wash water, etc. All of the process heat normally used in the system is produced from low penalty radioisotope heat sources. The analytical simulation was developed with the G189A computer program. The objective of the simulation was to obtain an analytical simulation which can be used to (1) evaluate the current RITE system steady state and transient performance during normal operating conditions, and also during off normal operating conditions including failure modes; and (2) evaluate the effects of variations in component design parameters and vehicle interface parameters on system performance.

Coggi, J. V.↗

Impulsive phase of solar flares. I - Characteristics of high energy electrons

The models describing the operative physical processes occurring during the impulsive phase of a solar flare can be broadly classified as thermal or nonthermal models. The current investigation is part of a series of studies which analyze the general characteristics of the 'nonthermal' models. The results are presented of an analysis regarding the steady state distribution of energetic electrons within the flare plasma. The time independent Fokker-Planck equation is used to determine the steady state distribution of the electrons in the ambient plasma. The model investigated encompasses thick-target, thin-target, and trap aspects of nonthermal hard X-ray models. High energy electrons with some initial energy spectrum and pitch angle distribution are injected at a height above the photosphere either into a closed magnetic loop structure or into an open field structure. Some approximate analytic solutions are discussed and the results of a numerical analysis are presented.

Leach, J.↗

Dynamics of magnetosphere-ionosphere coupling including turbulent transport

The dynamics of magnetosphere-ionosphere coupling has been investigated by means of a two-dimensional two-fluid MHD model including anomalous resistivity. When field-aligned current is generated on auroral field lines, the disturbance propagates toward the ionosphere in the form of a kinetic Alfven wave. When the current exceeds a critical value, microscopic turbulence is produced, which modifies the propagation of the Alfven wave. This process is modeled by a nonlinear collision frequency, which increases with the excess of the drift velocity over the critical value. The system evolves toward an electrostatic structure, with the perpendicular electric field having a shorter scale than the field-aligned current. The approach to a steady state is strongly dependent on the presence or absence of the turbulence and on the boundary conditions imposed in the generator. As current is increased or scale size is decreased, the turbulent region reflects and absorbs most of the Alfven wave energy, decoupling the generator from the ionosphere.

Lysak, R. L.↗

Overview of Sensitivity Analysis and Shape Optimization for Complex Aerodynamic Configurations

This paper presents a brief overview of some of the more recent advances in steady aerodynamic shape-design sensitivity analysis and optimization, based on advanced computational fluid dynamics. The focus here is on those methods particularly well- suited to the study of geometrically complex configurations and their potentially complex associated flow physics. When nonlinear state equations are considered in the optimization process, difficulties are found in the application of sensitivity analysis. Some techniques for circumventing such difficulties are currently being explored and are included here. Attention is directed to methods that utilize automatic differentiation to obtain aerodynamic sensitivity derivatives for both complex configurations and complex flow physics. Various examples of shape-design sensitivity analysis for unstructured-grid computational fluid dynamics algorithms are demonstrated for different formulations of the sensitivity equations. Finally, the use of advanced, unstructured-grid computational fluid dynamics in multidisciplinary analyses and multidisciplinary sensitivity analyses within future optimization processes is recommended and encouraged.

Newman, Perry A.↗