Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computing”

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 325 records · Page 18

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

Recent work in the aerothermodynamics community has focused on the development of reduced-order models for thermo-chemical non-equilibrium which avoid the restrictive assumptions of multi-temperature models and the prohibitive cost associated with State-to-State (StS) models. In the present work, this is accomplished by lumping energy states together and assuming groups of states are roughly in equilibrium. As a result, the non-equilibrium behavior of a gas can be captured at a reduced computational cost from a full StS simulation. In this work, we present a hybrid grouping model for studying energy transfer and dissociation in a mixture of nitrogen molecules due to N2-N2 reactions. This is accomplished by making use of a grouping strategy informed by data from the N2-N StS kinetic data. However, due to the massive computational cost associated with constructing StS data for the N2-N2 system, the kinetic data for the hybrid grouping model are calculated using the quasi-classical trajectory (QCT) method by sampling states for trajectory within the groups. This general framework is called the Maximum-Entropy Quasi-Classical Trajectory (ME-QCT) method. The primary challenge associated with this method is that rates for reverse grouped reactions cannot be obtained through detailed balance at a group level, due to the variation of group internal temperatures. To construct the full model for N2-N2 grouped kinetics using the ME-QCT method, detailed balance is invoked at the microscopic level, allowing for the calculation of the full kinetic data from QCT. Results will be presented using the full ME-QCT model for the N2-N2 system in an isothermal and isochoric reactor simulation. In addition, simple CFD test cases for a one-dimensional standing shock and for a quasi-one-dimensional nozzle will be used for demonstration of the ME-QCT method. This method allows for the calculation of non-equilibrium behavior for the N2-N2 system without the prohibitive cost of a full StS simulation. Moreover, it enables the construction of a unified model for the dissociating and recombining non-equilibrium flows.

nonequilibrium modeling↗

Fast Query-Optimized Kernel-Machine Classification

A recently developed algorithm performs kernel-machine classification via incremental approximate nearest support vectors. The algorithm implements support-vector machines (SVMs) at speeds 10 to 100 times those attainable by use of conventional SVM algorithms. The algorithm offers potential benefits for classification of images, recognition of speech, recognition of handwriting, and diverse other applications in which there are requirements to discern patterns in large sets of data. SVMs constitute a subset of kernel machines (KMs), which have become popular as models for machine learning and, more specifically, for automated classification of input data on the basis of labeled training data. While similar in many ways to k-nearest-neighbors (k-NN) models and artificial neural networks (ANNs), SVMs tend to be more accurate. Using representations that scale only linearly in the numbers of training examples, while exploring nonlinear (kernelized) feature spaces that are exponentially larger than the original input dimensionality, KMs elegantly and practically overcome the classic curse of dimensionality. However, the price that one must pay for the power of KMs is that query-time complexity scales linearly with the number of training examples, making KMs often orders of magnitude more computationally expensive than are ANNs, decision trees, and other popular machine learning alternatives. The present algorithm treats an SVM classifier as a special form of a k-NN. The algorithm is based partly on an empirical observation that one can often achieve the same classification as that of an exact KM by using only small fraction of the nearest support vectors (SVs) of a query. The exact KM output is a weighted sum over the kernel values between the query and the SVs. In this algorithm, the KM output is approximated with a k-NN classifier, the output of which is a weighted sum only over the kernel values involving k selected SVs. Before query time, there are gathered statistics about how misleading the output of the k-NN model can be, relative to the outputs of the exact KM for a representative set of examples, for each possible k from 1 to the total number of SVs. From these statistics, there are derived upper and lower thresholds for each step k. These thresholds identify output levels for which the particular variant of the k-NN model already leans so strongly positively or negatively that a reversal in sign is unlikely, given the weaker SV neighbors still remaining. At query time, the partial output of each query is incrementally updated, stopping as soon as it exceeds the predetermined statistical thresholds of the current step. For an easy query, stopping can occur as early as step k = 1. For more difficult queries, stopping might not occur until nearly all SVs are touched. A key empirical observation is that this approach can tolerate very approximate nearest-neighbor orderings. In experiments, SVs and queries were projected to a subspace comprising the top few principal- component dimensions and neighbor orderings were computed in that subspace. This approach ensured that the overhead of the nearest-neighbor computations was insignificant, relative to that of the exact KM computation.

Mazzoni, Dominic↗

Fully Associative, Nonisothermal, Potential-Based Unified Viscoplastic Model for Titanium-Based Matrices

A number of titanium matrix composite (TMC) systems are currently being investigated for high-temperature air frame and propulsion system applications. As a result, numerous computational methodologies for predicting both deformation and life for this class of materials are under development. An integral part of these methodologies is an accurate and computationally efficient constitutive model for the metallic matrix constituent. Furthermore, because these systems are designed to operate at elevated temperatures, the required constitutive models must account for both time-dependent and time-independent deformations. To accomplish this, the NASA Lewis Research Center is employing a recently developed, complete, potential-based framework. This framework, which utilizes internal state variables, was put forth for the derivation of reversible and irreversible constitutive equations. The framework, and consequently the resulting constitutive model, is termed complete because the existence of the total (integrated) form of the Gibbs complementary free energy and complementary dissipation potentials are assumed a priori. The specific forms selected here for both the Gibbs and complementary dissipation potentials result in a fully associative, multiaxial, nonisothermal, unified viscoplastic model with nonlinear kinematic hardening. This model constitutes one of many models in the Generalized Viscoplasticity with Potential Structure (GVIPS) class of inelastic constitutive equations.

Source record↗

Solid rocket motor aft field joint flow field analysis

An efficient Navier-Stokes analysis was successfully applied to simulate the complex flow field in the vicinity of a slot in a solid rocket motor with segment joints. The capability of the computer code to resolve the flow near solid surfaces without using a wall function assumption was demonstrated. In view of the complex nature of the flow field in the vicinity of the slot, this approach is considered essential. The results obtained from these calculations provide valuable design information, which would otherwise be extremely difficult to obtain. The results of the axisymmetric calculations indicate the presence of a region of reversed axial flow at the aft-edge of the slot and show the over-pressure in the slot to be only about 10 psi. The results of the asymmetric calculations indicate that a pressure asymmetry more than two diameters downstream of the slot has no noticeable effect on the flow field in the slot. They also indicate that the circumferential pressure differential caused in the slot due to failure of a 15 deg section of the castable inhibitor will be approximately 1 psi.

Sabnis, Jayant S.↗

Ablative Rocket Deflector Testing and Computational Modeling

A deflector risk mitigation program was recently conducted at the NASA Stennis Space Center. The primary objective was to develop a database that characterizes the behavior of industry-grade refractory materials subjected to rocket plume impingement conditions commonly experienced on static test stands. The program consisted of short and long duration engine tests where the supersonic exhaust flow from the engine impinged on an ablative panel. Quasi time-dependent erosion depths and patterns generated by the plume impingement were recorded for a variety of different ablative materials. The erosion behavior was found to be highly dependent on the material s composition and corresponding thermal properties. For example, in the case of the HP CAST 93Z ablative material, the erosion rate actually decreased under continued thermal heating conditions due to the formation of a low thermal conductivity "crystallization" layer. The "crystallization" layer produced near the surface of the material provided an effective insulation from the hot rocket exhaust plume. To gain further insight into the complex interaction of the plume with the ablative deflector, computational fluid dynamic modeling was performed in parallel to the ablative panel testing. The results from the current study demonstrated that locally high heating occurred due to shock reflections. These localized regions of shock-induced heat flux resulted in non-uniform erosion of the ablative panels. In turn, it was observed that the non-uniform erosion exacerbated the localized shock heating causing eventual plume separation and reversed flow for long duration tests under certain conditions. Overall, the flow simulations compared very well with the available experimental data obtained during this project.

Allgood, Daniel C.↗

The numerical calculation of laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for several separated flows and the accuracy of the method is verified. A computer program listing and complete solution case are included.

Klineberg, J. M.↗

On laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for a number of separated flows and the accuracy of the method is verified.

Klineberg, J. M.↗

Thickening oscillation of a delta wing using Navier-Stokes and Navier-displacement equations

The problem of unsteady, supersonic, locally conical, vortical flow around a delta wing undergoing thickening oscillation is solved using the unsteady, thin-layer Navier-Stokes equations and the unsteady linearized Navier-displacement equations. The unsteady, thin-layer Navier-Stokes equations are solved using the implicit approximate-factorization finite-volume scheme to compute the conservative components of the flow vector field. With the conservative components known at any time step, the linearized, Navier-displacement equations are solved using the alternating, direction-implicit scheme to obtain the grid points displacements due to known displacement boundary conditions. A grid-displacements limiter, in the form of a low mesh Reynolds number, is used to limit grid-folding in regions of highly reversed flow.

Chuang, Hsin-Kung A.↗

Wake Velocimetry of a Sphere-Cone Model in a Mach 10 Air Freestream using FLEET

Velocity measurements in the wake behind a sphere-cone-shaped vehicle wind tunnel model in a Mach 10 hypersonic air flow using femtosecond laser electronic excitation tagging (FLEET) are reported in this work. The FLEET measurements were performed at 1 kHz using a femtosecond laser centered at 800 nm and an ultrahigh-speed intensified camera system. The results reported here are for a single freestream condition corresponding to approximately Re ∞ /L = 3.6∙10 6 m -1 . The performance of the velocimetry instrument was evaluated in the quiescent test section at conditions relevant to the wake measurements. For velocity measurements in the wake, the FLEET measurement volume was translated to several radial locations from the model centerline in the vertical direction and at a streamwise location corresponding to approximately 12 mm from the payload. Reverse flow with mean velocities ranging from -12 m/s to -48 m/s were observed in the recirculating region of the wake, whereas hypersonic velocities in the range of 1079 m/s to 1183 m/s were observed in the free shear layer. A velocity uncertainty methodology is also outlined and applied for the reported measurements. The velocity data reported in this work is expected to play a significant role in the validation of computational codes modeling the hard-to-predict separated hypersonic wake.

FLEET↗

Wake Velocimetry of a Sphere-Cone Model in a Mach 10 Air Freestream using FLEET

Velocity measurements in the wake behind a sphere-cone-shaped vehicle wind tunnel model in a Mach 10 hypersonic air flow using femtosecond laser electronic excitation tagging (FLEET) are reported in this work. The FLEET measurements were performed at 1 kHz using a femtosecond laser centered at 800 nm and an ultrahigh-speed intensified camera system. The results reported here are for a single freestream condition corresponding to approximately Re ∞ /L = 3.6∙10 6 m -1 . The performance of the velocimetry instrument was evaluated in the quiescent test section at conditions relevant to the wake measurements. For velocity measurements in the wake, the FLEET measurement volume was translated to several radial locations from the model centerline in the vertical direction and at a streamwise location corresponding to approximately 12 mm from the payload. Reverse flow with mean velocities ranging from -12 m/s to -48 m/s were observed in the recirculating region of the wake, whereas hypersonic velocities in the range of 1079 m/s to 1183 m/s were observed in the free shear layer. A velocity uncertainty methodology is also outlined and applied for the reported measurements. The velocity data reported in this work is expected to play a significant role in the validation of computational codes modeling the hard-to-predict separated hypersonic wake.

FLEET↗

Temporal behavior of density perturbations in the polar wind

The time-dependent continuity and momentum equations for H(+) are solved in order to study the temporal evolution of density perturbations in the supersonic, collisionless polar wind. The results indicate that the sophisticated computer models of small-scale plasma expansions are indeed applicable to large-scale flow problems. The important prediction of the production of superthermal ions from small-scale simulations by the process of plasma expansion is borne out by the calculations in this paper. Energetic ions are produced for each type of density perturbation through the plasma expansion process. In the case of an extended density depletion, the forward-reverse shock pair forms when the depletion is very strong, but otherwise only the forward shock forms. In the case of localized density depletions, the density cavity evolves into two rarefaction waves joined by a forward-reverse shock pair. In the case of localized density enhancements, the shock pair evolves from compressive waves.

Singh, N.↗

Application of a reversible chemical reaction system to solar thermal power plants

Three distributed dish solar thermal power systems using various applications of SO2/SO3 chemical energy storage and transport technology were comparatively assessed. Each system features various roles for the chemical system: (1) energy storage only, (2) energy transport, or (3) energy transport and storage. These three systems were also compared with the dish-Stirling, using electrical transport and battery storage, and the central receiver Rankine system, with thermal storage, to determine the relative merit of plants employing a thermochemical system. As an assessment criterion, the busbar energy costs were compared. Separate but comparable solar energy cost computer codes were used for distributed receiver and central receiver systems. Calculations were performed for capacity factors ranging from 0.4 to 0.8. The results indicate that SO2/SO3 technology has the potential to be more cost effective in transporting the collected energy than in storing the energy for the storage capacity range studied (2-15 hours)

Hanseth, E. J.↗

Cross-B convection of artificially created, negative-ion clouds and plasma depressions - Low-speed flow

A negative-ion, positive-ion plasma produced by the release of an electron attachment chemical into the F region becomes electrically polarized by collisions with neutrals moving across magnetic field lines. The resulting electric field causes E x B drift of the two ion species and the residual electrons. The cross-field flow of the modified ionosphere is computed using a two-dimensional numerical simulation which includes electron attachment and mutual neutralization chemistry, self-consistent electric fields, and three-species plasma transport. The velocity of the plasma is initially in the direction of the neutral wind because the negative-ion cloud is a Pedersen conductivity enhancement. As the positive and negative ions react, the Pedersen conductivity becomes depressed below the ambient value and the velocity of the plasma reverses direction. A plasma hole remains after the positive and negative ions have mutually neutralized. The E x B gradient drift instability produces irregularities on the upwind edge of the hole.

Bernhardt, Paul A.↗

Numerical, experimental, and theoretical study of convective instability of flows over pointed bodies at incidence

A study is conducted to investigate whether the behavior of the asymmetric mean flow observed on pointed bodies of revolution at incidence remains consistent with the presence of a convective instability of an original symmetric flow, even in the incidence range where virtually bistable behavior of the asymmetric flow is observed. By means of a retractable wire located near the tip, it is determined experimentally that for all angles of attack tested (30 to 60 degs), changing the size or location of the controlled disturbance results in a finite change in the asymmetric flow field, even to the extent of reversing the sign of the side force or becoming almost symmetric. The process is reversible; returning the wire to an original position likewise restores the corresponding flow field and mean side force. The effect of wire location (roll angle and distance from the tip) as well as angle of attack and flow conditions are evaluated experimentally by means of flow visualization and side-force measurements for a generic ogive-cylinder model in a low-speed wind tunnel. Evaluation of the results in the light of computational observations and theoretical considerations yields an affirmative answer to the question posed.

Degani, David↗

Operator versus computer control of adaptive automation

Adaptive automation refers to real-time allocation of functions between the human operator and automated subsystems. The article reports the results of a series of experiments whose aim is to examine the effects of adaptive automation on operator performance during multi-task flight simulation, and to provide an empirical basis for evaluations of different forms of adaptive logic. The combined results of these studies suggest several things. First, it appears that either excessively long, or excessively short, adaptation cycles can limit the effectiveness of adaptive automation in enhancing operator performance of both primary flight and monitoring tasks. Second, occasional brief reversions to manual control can counter some of the monitoring inefficiency typically associated with long cycle automation, and further, that benefits of such reversions can be sustained for some time after return to automated control. Third, no evidence was found that the benefits of such reversions depend on the adaptive logic by which long-cycle adaptive switches are triggered.

Hilburn, Brian↗

Spectral synthesis in the ultraviolet. III - The spectral morphology of normal stars in the mid-ultraviolet

The morphology of 218 mid-UV spectra of stars ranging from O through K in spectral type is examined. Several new line and continuum indices are defined and their usefulness as temperature, luminosity, and metallicity discriminants is discussed. Mid-UV stellar continua are found to be markedly affected by abundance. A UV excess, delta(2600-V), is computed which is more sensitive by a factor of 10 to (Fe/H) than is delta(U-B). The relative strength of spectral lines in the mid-UV is not as strongly affected by abundance. Mid-UV spectra are much more sensitive to the temperature of the stellar population than to either metallicity or the dwarf/giant ratio. Mg II 2800 shows unexpected behavior, displaying no sensitivity to abundance for cool stars and a reversed sensitivity in FG dwarfs such that metal-poor stars have stronger Mg II strengths at the same temperature than more metal-rich stars.

Fanelli, Michael N.↗

A spectral algorithm for envelope reduction of sparse matrices

A new algorithm for reducing the envelope of a sparse matrix is presented. This algorithm is based on the computation of eigenvectors of the Laplacian matrix associated with the graph of the sparse matrix. A reordering of the sparse matrix is determined based on the numerical values of the entries of an eigenvector of the Laplacian matrix. Numerical results show that the new reordering algorithm can in some cases reduce the envelope by more than a factor of two over the current standard algorithms such as Gibbs-Poole-Stockmeyer (GPS) or SPARSPAK's reverse Cuthill-McKee (RCM).

Barnard, Stephen T.↗

Software Safety Risk in Legacy Safety-Critical Computer Systems

Safety-critical computer systems must be engineered to meet system and software safety requirements. For legacy safety-critical computer systems, software safety requirements may not have been formally specified during development. When process-oriented software safety requirements are levied on a legacy system after the fact, where software development artifacts don't exist or are incomplete, the question becomes 'how can this be done?' The risks associated with only meeting certain software safety requirements in a legacy safety-critical computer system must be addressed should such systems be selected as candidates for reuse. This paper proposes a method for ascertaining formally, a software safety risk assessment, that provides measurements for software safety for legacy systems which may or may not have a suite of software engineering documentation that is now normally required. It relies upon the NASA Software Safety Standard, risk assessment methods based upon the Taxonomy-Based Questionnaire, and the application of reverse engineering CASE tools to produce original design documents for legacy systems.

Hill, Janice↗