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 127 records · Page 7

Physics-Based Modeling and Simulation of Self-Reacting Friction Stir Welding Using Computational Fluid Dynamics

A physics-based model was developed to simulate the behavior of material in a self-reacting friction stir welding (SR-FSW) process for the joining of metals. This steady-state model builds upon fundamental computational fluid dynamic (CFD) principles within Ansys Fluent to solve the discretized equations. The effective viscosity is calculated using a viscoplastic model using a Sheppard-Wright formulation of flow stress. Numerous advancements have been made in the incorporated physics including (1) temperature-dependent material properties; (2) locally adaptable flow and thermal boundary conditions; and (3) adapting material properties in nugget in response to microstructural changes. Simulation strategies to accelerate computation and improve numerical stability include adapting the mesh refinement and solver relaxation factors during simulation. The result is a highly robust and computationally efficient model capable of providing the material flow and temperature history across the domain. As material history determines the local microstructure and ultimately weld strength, an accurate and detailed physics-based model has the potential to accelerate SR-FSW process development. The model is highly adaptable to changes in process parameters, tool design, or alloy.

Process Modeling↗

Role of Dynamic Nucleation at Moving Boundaries in Phase and Microstructure Selection

Solidification microstructures that form under steady-state growth conditions (cells, dendrites, regular eutectics, etc.) are reasonably well understood in comparison to other, more complex microstructures, which form under intrinsically non-steady-state growth conditions due to the competition between the nucleation and growth of several phases. Some important practical examples in this latter class include microstructures forming in peritectic systems in highly undercooled droplets, and in strip cast stainless steels. Prediction of phase and microstructure selection in these systems has been traditionally based on (1) heterogeneous nucleation on a static interface, and (2) comparing the relative growth rate of different phase/microstructures under steady-state growth conditions. The formation of new phases, however, occurs via nucleation on, or ahead of, a moving boundary. In addition, the actual selection process is controlled by a complex interaction between the nucleation process and the growth competition between the nuclei and the pre-existing phase under non-steady-state conditions. As a result, it is often difficult to predict which microstructure will form and which phases will be selected under prescribed processing conditions. This research addresses this critical role of nucleation at moving boundaries in the selection of phases and solidification microstructures through quantitative experiments and numerical modeling in peritectic systems. In order to create a well characterized system in which to study this problem, we focus on the directional solidification of hypo- and hyper-peritectic alloys in the two-phase region, imposing a large enough ratio of temperature gradient/growth rate (G/V(sub p)) to suppress the morphological instability of both the parent (alpha) and peritectic (Beta) phases, i.e. each phase alone would grow as a planar front. Our combined experimental and theoretical results show that, already in this simplified case, the growth competition of these two phases leads to a rich variety of microstructures that depend sensitively upon the relative importance of nucleation, diffusion, and convection.

Karma, Alain↗

Advantages of multigrid methods for certifying the accuracy of PDE modeling

Numerical techniques for assessing and certifying the accuracy of the modeling of partial differential equations (PDE) to the user's specifications are analyzed. Examples of the certification process with conventional techniques are summarized for the three dimensional steady state full potential and the two dimensional steady Navier-Stokes equations using fixed grid methods (FG). The advantages of the Full Approximation Storage (FAS) scheme of the multigrid technique of A. Brandt compared with the conventional certification process of modeling PDE are illustrated in one dimension with the transformed potential equation. Inferences are drawn for how MG will improve the certification process of the numerical modeling of two and three dimensional PDE systems. Elements of the error assessment process that are common to FG and MG are analyzed.

Forester, C. K.↗

Magnetic Fluctuations in the Jovian Magnetosphere

The engine that drives the jovian magnetosphere is the mass added to the Io ion torus, accelerated to corotational velocities by field-aligned currents that couple the Io torus to the jovian ionosphere. The mass of the torus builds up to an amount that the magnetic forces cannot contain and the plasma, first slowly and then more rapidly, drifts outward. Numerous authors have treated this problem based first on the observations of the Pioneer 10 and 11 flybys; then on Voyager 1 and 2, and Ulysses; and finally most recently the Galileo orbiter. The initial observations revealed the now familiar magnetodisk, in which the field above and below the magnetic equator became quite radial in orientation and much less dipolar. The Galileo observations show this transformation to occur on average at 24 R(sub J) and to often be quite abrupt. These observations are consistent with outward transport of magnetized plasma that moves ever faster radially until about 50 R(sub J) on the nightside where the field lines stretch to the breaking point, reconnection occurs, and plasma and field islands are transported down the tail ultimately removing the mass from the magnetosphere that Io had deposited deep in the inner torus. The reconnection process creates empty flux tubes connected to Jupiter that are buoyant and thought to float inward and replace the flux carried out with the torus plasma. As described above, the jovian magnetosphere could very well be in a state of steady laminar circulation, but indeed it is not. The process is very unsteady and the wave levels can be very intense. The existence of these waves in turn can lead to processes that compete with the radial circulation pattern in removing plasma from the system. These waves can scatter particles so that they precipitate into the ionosphere. This process should be important in the Io torus where the atmospheric loss cone is relatively large and becomes less important as the loss cone decreases in size with radial distance. However, the Io torus is relatively quiet compared to the region outside the torus and it is not obvious without studying this scattering carefully whether the loss in the torus or out of the torus is greater and whether it can act rapidly enough to compete with the radial transport of ions to the tail in the life cycle of the mass added at Io. Closer to Io the ion cyclotron waves are most intense and possibly are associated with the losses in the Io flux tube. The waves are also diagnostic of both the Io atmospheric composition and the size and strength of the massloading process.

Russell, Christopher T.↗

Initial results from the NASA Lewis Bumpy Torus experiment

Initial results were obtained from low power operation of the NASA Lewis Bumpy Torus experiment, in which a steady-state ion heating method based on the modified Penning discharge is applied in a bumpy torus confinement geometry. The magnet facility consists of 12 superconducting coils, each 19 cm i.d. and capable of 3.0 T, equally spaced in a toroidal array 1.52 m in major diameter. A 18 cm i.d. anode ring is located at each of the 12 midplanes and is maintained at high positive potentials by a dc power supply. Initial observations indicate electron temperatures from 10 to 150 eV, and ion kinetic temperatures from 200 eV to 1200 eV. Two modes of operation were observed, which depend on background pressure, and have different radial density profiles. Steady state neutron production was observed. The ion heating process in the bumpy torus appears to parallel closely the mechanism observed when the modified Penning discharge was operated in a simple magnetic mirror field.

Roth, J. R.↗

Nonlinear combustion instability model in two- to three-dimensions

A model is developed on using quasi-steady state equations for the injection, atomization, and vaporization processes that consider the instantaneous local pressures, temperatures, density, and gas velocity vector associated with a 3-D wave. The coupled sets of equations are solved to determine frequency and growth rate of the oscillations as a function of wave amplitude. Calculated results are used to explain the observed instabilities in a 40 Klb LOX/Methane combustor.

Priem, Richard J.↗

Low power arcjet system spacecraft impacts

Application of electrothermal arcjets on communications satellites requires assessment of integration concerns identified by the user community. Perceived risks include plume contamination of spacecraft materials, induced arcing or electrostatic discharges between differentially charged spacecraft surfaces, and conducted and radiated electromagnetic interference (EMI) for both steady state and transient conditions. A Space Act agreement between Martin Marietta Astro Space, the Rocket Research Company, and NASA's Lewis Research Center was established to experimentally examine these issues. Spacecraft materials were exposed to an arcjet plume for 40 hours, representing 40 weeks of actual spacecraft life, and contamination was characterized by changes in surface properties. With the exception of the change in emittance of one sample, all measurable changes in surface properties resulted in acceptable end of life characteristics. Charged spacecraft samples were benignly and consistently reduced to ground potential during exposure to the powered arcjet plume, suggesting that the arcjet could act as a charge control device on spacecraft. Steady state EMI signatures obtained using two different power processing units were similar to emissions measured in a previous test. Emissions measured in UHF, S, C, Ku and Ka bands obtained a null result which verified previous work in the UHF, S, and C bands. Characteristics of conducted and radiated transient emissions appear within standard spacecraft susceptibility criteria.

Pencil, Eric J.↗

Structural dynamic measurement practices for turbomachinery at the NASA Lewis Research Center

Methods developed for measuring blade and rotor-shaft system response include optical systems, transient instruments, and special digital data processing equipment. Optical methods offer some distinct benefits for blade vibration measurement. Transient and steady state measurements of the response of rotor-shaft systems strongly affect analytical methods development. Digital computing systems allow processing of large volumes of high speed data from rotating blade sets. Also, digital systems develop useful vibration response signatures from randomly excited systems. Research facilities include the spin rig facility and the transient rotor response lab.

Kiraly, L. J.↗

Parameterization of planetary wave breaking in the middle atmosphere

A parameterization of planetary wave breaking in the middle atmosphere has been developed and tested in a numerical model which includes governing equations for a single wave and the zonal-mean state. The parameterization is based on the assumption that wave breaking represents a steady-state equilibrium between the flux of wave activity and its dissipation by nonlinear processes, and that the latter can be represented as linear damping of the primary wave. With this and the additional assumption that the effect of breaking is to prevent further amplitude growth, the required dissipation rate is readily obtained from the steady-state equation for wave activity; diffusivity coefficients then follow from the dissipation rate. The assumptions made in the derivation are equivalent to those commonly used in parameterizations for gravity wave breaking, but the formulation in terms of wave activity helps highlight the central role of the wave group velocity in determining the dissipation rate. Comparison of model results with nonlinear calculations of wave breaking and with diagnostic determinations of stratospheric diffusion coefficients reveals remarkably good agreement, and suggests that the parameterization could be useful for simulating inexpensively, but realistically, the effects of planetary wave transport.

Garcia, Rolando R.↗

Adaptive Optimization of Aircraft Engine Performance Using Neural Networks

Preliminary results are presented on the development of an adaptive neural network based control algorithm to enhance aircraft engine performance. This work builds upon a previous National Aeronautics and Space Administration (NASA) effort known as Performance Seeking Control (PSC). PSC is an adaptive control algorithm which contains a model of the aircraft's propulsion system which is updated on-line to match the operation of the aircraft's actual propulsion system. Information from the on-line model is used to adapt the control system during flight to allow optimal operation of the aircraft's propulsion system (inlet, engine, and nozzle) to improve aircraft engine performance without compromising reliability or operability. Performance Seeking Control has been shown to yield reductions in fuel flow, increases in thrust, and reductions in engine fan turbine inlet temperature. The neural network based adaptive control, like PSC, will contain a model of the propulsion system which will be used to calculate optimal control commands on-line. Hopes are that it will be able to provide some additional benefits above and beyond those of PSC. The PSC algorithm is computationally intensive, it is valid only at near steady-state flight conditions, and it has no way to adapt or learn on-line. These issues are being addressed in the development of the optimal neural controller. Specialized neural network processing hardware is being developed to run the software, the algorithm will be valid at steady-state and transient conditions, and will take advantage of the on-line learning capability of neural networks. Future plans include testing the neural network software and hardware prototype against an aircraft engine simulation. In this paper, the proposed neural network software and hardware is described and preliminary neural network training results are presented.

Simon, Donald L.↗

Subsonic flight test evaluation of a propulsion system parameter estimation process for the F100 engine

An adaptive-performance-seeking control system which optimizes the quasi-steady-state performance of the F-15 propulsion system is discussed. This paper presents flight- and ground-test evaluations of the propulsion system parameter-estimation process used by the performance seeking control system. The estimator consists of a compact propulsion system model and an extended Kalman filter. The extended Kalman filter estimates five engine component deviation parameters from measured inputs. The compact model uses measurements and Kalman-filter estimates as inputs to predict unmeasured propulsion parameters such as net propulsive force and fan stall margin. The ability to track trends and estimate absolute values of propulsion system parameters was demonstrated. For example, thrust stand results show a good correlation especially in trends between the performance seeking control estimated and measured thrust.

Orme, John S.↗

Cascaded coherent tracking systems with time-varying channels

A performance analysis is presented of a two-way coherent tracking system (of the type used for spacecraft tracking and navigation) in which the transmitted signals have pass through linear time-varying channels. The performance of the system is characterized by the steady state probability density function of the reduced phase error process in the second tracking loop, the system considered consisting of two first-order phase locked loops in cascade. While the log normal channel, which arises in communication through planetary atmospheres, was used as the channel model, the results can be extended to other channels, such as the Rice and Rayleigh channels.

Weber, W. J., III↗

On the maintenance of the Marine boundary layer-free tropospheric gradient of transient tracers

A simple model of the marine boundary layer is used to arrive at a steady-state budget for trace gases in this layer. Five processes maintain the tracer balance in the layer. The first is a large-scale horizontal convergence of the tracer; the second is the pumping by clouds of the tracer out of the boundary layer into the free troposphere; the third is penetration by the tracer from the free troposphere into the marine boundary from compensating subsidence in the cloud-free air surrounding cumulus clouds; the fourth is sea surface emission or deposition of the tracer; and the fifth is in situ photochemical production or loss of the tracer within the marine boundary layer. The budget equation is used in deriving an expression for the mixing ratio difference between the two regions. A scale analysis suggests that the mixing ratio difference between the marine boundary layer and the free troposphere can be an order of magnitude smaller in disturbed weather conditions than in suppressed conditions, showing that the difference in chemical concentrations between the free troposphere and the marine boundary layer cannot be rationalized solely on the basis of chemical measurements.

Gidel, L. T.↗

Excitation of nearly steady finite-amplitude barotropic waves

An exact nonlinear stationary solution is obtained for barotropic waves in a beta-plane channel and it is shown that it can be excited under a range of initial conditions. Results show that a finite-amplitude wave in a constant shear flow, given an initial phase tilt against the shear and a sufficient initial amplitude, interacts with the mean flow to produce a nearly steady state close to the exact stationary solution. This equilibration process involves nonlinear transients; in particular, as the flow equilibrates, the emergence of critical levels is accompanied by the neutralization of local mean vorticity gradients at these levels, thus allowing the solution to attain a nonsingular modal structure.

Hou, A. Y.↗

Advancements in silicon web technology

Low defect density silicon web crystals up to 7 cm wide are produced from systems whose thermal environments are designed for low stress conditions using computer techniques. During growth, the average silicon melt temperature, the lateral melt temperature distribution, and the melt level are each controlled by digital closed loop systems to maintain thermal steady state and to minimize the labor content of the process. Web solar cell efficiencies of 17.2 pct AM1 have been obtained in the laboratory while 15 pct efficiencies are common in pilot production.

Hopkins, R. H.↗

Distributed simulation of network protocols

Simulations of high speed network protocols are very CPU intensive operations requiring very long run times. Very high speed network protocols (Gigabit/sec rates) require longer simulation runs in order to reach a steady state, while at the same time requiring additional CPU processing for each unit of time because of the data rates for the traffic being simulated. As protocol development proceeds and simulations provide insights into any problems associated with the protocol, the simulation model often must be changed to generate additional or finer statistical performance information. Iterating on this process is very time consuming due to the required run times for the simulation models. The results of the efforts to distribute a high speed ring network protocol, Carrier Sensed Multiple Access/Ring Network (CSMA/RN), are presented.

Paterra, Frank↗

Interhemispheric asymmetry of the high-latitude ionospheric convection pattern

The assimilative mapping of ionospheric electrodynamics technique has been used to derive the large-scale high-latitude ionospheric convection patterns simultaneously in both northern and southern hemispheres during the period of January 27 to 29, 1992. When the interplanetary magnetic field (IMF) B(sub z) component is negative, the convection patterns in the southern hemisphere are basically the mirror images of those in the northern hemisphere. The total cross-polar cap potential drops in the two hemispheres are similar. When B(sub z) is positive and absolute value of B(sub y) greater than B(sub z), the convection configurations are mainly determined by B(sub y) and they may appear as normal 'two-cell' patterns in both hemispheres much as one would expect under southward IMF conditions. However, there is a significant difference in the cross-polar-cap potential drop between the two hemispheres, with the potential drop in the southern (summer) hemisphere over 50% larger than that in the northern (winter) hemisphere. As the ratio of absolute value of B(sub y)/B(sub z) decreases (less than one), the convection configuration in the two hemispheres may be significantly different, with reverse convection in the southern hemisphere and weak but disturbed convection in the northern hemisphere. By comparing the convection patterns with the corresponding spectrograms of precipitating particles, we interpret the convection patterns in terms of the concept of merging cells, lobe cells, and viscous cells. Estimates of the ' merging cell' potential drops, that is, the potential ascribed to the opening of the dayside field lines, are usually comparable between the two hemispheres, as they should be. The 'lobe cell' provides a potential between 8.5 and 26 kV and can differ greatly between hemispheres, as predicted. Lobe cells can be significant even for southward IMF, if absolute value of B(sub y) greater than the absolute value of B(sub z). To estimate the potential drop of the 'viscous cells,' we assume that the low-latitude boundary layer is on closed field lines. We find that this potential drop varies from case to case, with a typical value of 10 kV. If the source of these cells is truly a viscous interaction at the flank of the magnetopause, the process is likely spatially and temporally varying rather than steady state.

Lu, G.↗