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 397 records · Page 22

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO reversible reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HNO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) value for HNO using state-of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +/- 0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗

Afferent innervation patterns of the saccule in pigeons

The innervation patterns of vestibular saccular afferents were quantitatively investigated in pigeons using biotinylated dextran amine as a neural tracer and three-dimensional computer reconstruction. Type I hair cells were found throughout a large portion of the macula, with the highest density observed in the striola. Type II hair cells were located throughout the macula, with the highest density in the extrastriola. Three classes of afferent innervation patterns were observed, including calyx, dimorph, and bouton units, with 137 afferents being anatomically reconstructed and used for quantitative comparisons. Calyx afferents were located primarily in the striola, innervated a number of type I hair cells, and had small innervation areas. Most calyx afferent terminal fields were oriented parallel to the anterior-posterior axis and the morphological polarization reversal line. Dimorph afferents were located throughout the macula, contained fewer type I hair cells in a calyceal terminal than calyx afferents and had medium sized innervation areas. Bouton afferents were restricted to the extrastriola, with multi-branching fibers and large innervation areas. Most of the dimorph and bouton afferents had innervation fields that were oriented dorso-ventrally but were parallel to the neighboring reversal line. The organizational morphology of the saccule was found to be distinctly different from that of the avian utricle or lagena otolith organs and appears to represent a receptor organ undergoing evolutionary adaptation toward sensing linear motion in terrestrial and aerial species.

NASA Discipline Neuroscience↗

Capillary Pressure Balance with Repeated Head-Up Tilt in Humans

We measured blood density (BD). plasma density (PD) and hematocrit (Ht) during repeated 70 deg head-up tilt (HUT) consisting of a pre-drink period with two supine (P and P3) and two HUT (P2 and P4) phases of 45 min each. At the end of P4 test subjects (N=8) drank 10 ml/kg body weight of isotonic (290 mOsm/kg) sodium chloride (ISO) or hypotonic (<10 mOsm/kg) unsweetened us (HYP), or nothing (control: CON). The following periods continued the supine (P5, P7) / upright (P6) sequence. Mass density of plasma and whole blood using a "mechanical oscillation" technique were measured from ear lobe blood. Tilting altered (p<0.001) BD, PD, and Ht in all consecutive Phases. BD and PD were different (p<0.05) between CON. HYP. and ISO Phases 6 and 7. The computed mass density of fluid (FD). which moved between intra. and extravascular compartments. averaged 1008 g/l and was the same in any conditions. In the non-drink treatment (P3, P5, P7), supine PV steadly decreased when compared to P1. Tilt-induced PV shifts ranged from 9.7 to 16.7% when compared to PV du. ring the respective previous Phases. After drinking, PV increased above CON values at the end of Phase 7 by 12.9% with ISO, and by 6.6% with HYP. Progressive hemoconcentration. occurred during non-drink supine; isotonic saline ingestion increased supine PV a control level but did not stop or reverse the decrease of upright hemoconcentration. Upright density values were not slowed with repeated upright tilting; decreased fluid loss occurred in consecutive upright periods. Upright capillary pressure balance was achieved with diminished plasma volume loss during repeated tilting.

Hinghofer-Szalkay, H. G.↗

Matlab GUI for a Fluid Mixer

The Test and Engineering Directorate at NASA John C. Stennis Space Center developed an interest to study the modeling, evaluation, and control of a liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem of a ground test facility. This facility carries out comprehensive ground-based testing and certification of liquid rocket engines including the Space Shuttle Main engine. A software simulation environment developed in MATLAB/SIMULINK (M/S) will allow NASA engineers to test rocket engine systems at relatively no cost. In the progress report submitted in February 2004, we described the development of two foundation programs, a reverse look-up application using various interpolation algorithms, a variety of search and return methods, and self-checking methods to reduce the error in returned search results to increase the functionality of the program. The results showed that these efforts were successful. To transfer this technology to engineers who are not familiar with the M/S environment, a four-module GUI was implemented allowing the user to evaluate the mixer model under open-loop and closed-loop conditions. The progress report was based on an udergraduate Honors Thesis by Ms. Jamie Granger Austin in the Department of Electrical Engineering and Computer Science at Tulane University, during January-May 2003, and her continued efforts during August-December 2003. In collaboration with Dr. Hanz Richter and Dr. Fernando Figueroa we published these results in a NASA Tech Brief due to appear this year. Although the original proposal in 2003 did not address other components of the test facility, we decided in the last few months to extend our research and consider a related pressurization tank component as well. This report summarizes the results obtained towards a Graphical User Interface (GUI) for the evaluation and control of the hydrogen mixer subsystem model and for the pressurization tank each taken individually. Further research would combine the two components - mixer and tank, for a more realistic simulation tool.

Barbieri, Enrique↗

Constraints and Opportunities in GCM Model Development

Over the past 30 years climate models have evolved from relatively simple representations of a few atmospheric processes to complex multi-disciplinary system models which incorporate physics from bottom of the ocean to the mesopause and are used for seasonal to multi-million year timescales. Computer infrastructure over that period has gone from punchcard mainframes to modern parallel clusters. Constraints of working within an ever evolving research code mean that most software changes must be incremental so as not to disrupt scientific throughput. Unfortunately, programming methodologies have generally not kept pace with these challenges, and existing implementations now present a heavy and growing burden on further model development as well as limiting flexibility and reliability. Opportunely, advances in software engineering from other disciplines (e.g. the commercial software industry) as well as new generations of powerful development tools can be incorporated by the model developers to incrementally and systematically improve underlying implementations and reverse the long term trend of increasing development overhead. However, these methodologies cannot be applied blindly, but rather must be carefully tailored to the unique characteristics of scientific software development. We will discuss the need for close integration of software engineers and climate scientists to find the optimal processes for climate modeling.

Schmidt, Gavin↗

Towards Accurate and Efficient Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from First Principles

Shape memory alloys (SMAs) can remember and recover their original shapes upon heating due to the existence of a reversible martensitic transition (MT) between the high-temperature austenite (A) and low-temperature martensite (M) phases. The martensitic transition temperature (MTT) is a crucial characteristic of an SMA. SMAs have a wide range of potential applications in aerospace, civil engineering, bioengineering, etc., but their operating temperatures are limited by the available SMAs. MTT can be tuned by alloying a binary with other metals, and the multicomponent NiTi-based SMAs have attracted tremendous research efforts recently. It is not efficient to employ the trial-and-error method alone due to the dramatically increased complexity and possibilities in compositions, and thus reliable theory and accurate computations play an indispensable role in creating SMAs with desirable properties.

Zhigang Wu↗

Multiobjective Multidisciplinary Optimization of Low-Boom Supersonic Transports Using Multifidelity Models

A multidisciplinary optimization (MDO) method has been developed to design a computational fluid dynamics (CFD) based low-boom configuration that can be obtained from a Pareto solution of a low-fidelity multiobjective MDO problem with mission constraints. This paper refines the developed MDO method by using multifidelity models for CFD-based multiobjective MDO and a better method for the system-level trade between the target low boom ground noise level and the overland range. The refined MDO method can generate a low-boom configuration that satisfies the mission requirements, has the lowest takeoff gross weight and the longest range for the low-boom mission, trims the low-boom cruise flight with fuel redistributions, and has a reversed equivalent area distribution closely matching a low-boom target with ground noise level below 70 PLdB. The validity of the refined MDO method is demonstrated by a design study of a low-boom supersonic transport that carries 40 passengers, flies a low-boom mission with cruise Mach of 1.7 and range of 3500 nm, and cruises overwater at Mach 1.8 with range of 3882 nm. Moreover, the refined MDO method eliminates the difference between the assumed cruise weight for CFD-based low-boom inverse design optimization and the estimated cruise weight of the optimal inverse design solution with respect to the mission requirements.

MDO↗

Isotropization and Evolution of Energy-Containing Eddies in Solar Wind Turbulence: Parker Solar Probe, Helios 1, ACE, WIND, and Voyager 1

We examine the radial evolution of correlation lengths perpendicular (l^ C) and parallel (  lC) to the magnetic-field direction, computed from solar wind magnetic-field data measured by Parker Solar Probe (PSP) during its first eight orbits, Helios 1, Advanced Composition Explorer (ACE), WIND, and Voyager 1 spacecraft. Correlation lengths are grouped by an interval’s alignment angle; the angle between the magnetic-field and solar wind velocity vectors (ΘBV). Parallel and perpendicular angular channels correspond to angles 0° < ΘBV < 40° and 50° < ΘBV < 90°, respectively. We observe an anisotropy in the inner heliosphere within 0.40 au, with  l l »^ 0.75 C C at 0.10 au. This anisotropy reduces with increasing heliocentric distance and the correlation lengths roughly isotropize within 1 au. Results from ACE and WIND support a reversal of the anisotropy, such that  l l »^ 1.29 C C at 1 au. The ratio does not appear to change significantly beyond 1 au, although the small number of parallel intervals in the Voyager data set precludes unambiguous conclusions from being drawn. This study provides insights regarding the radial evolution of the large, most energetic interacting turbulent fluctuations in the heliosphere. We also emphasize the importance of tracking the changes in sampling direction in PSP measurements as the spacecraft approaches the Sun, when using these data to study the radial evolution of turbulence. This can prove to be vital in understanding the more complex dynamics of the solar wind in the inner heliosphere and can assist in improving related simulations.

Manuel Enrique Cuesta↗

Large-scale structure of the heliospheric magnetic field - 1976-1991

A uniform set of computed large-scale coronal fields derived from photospheric field observations now exists for the interval 1976 to the present. Using this data base we can begin to compare Solar Cycle 21 with Solar Cycle 22. The character of the field and its evolution are the same, but many of the particulars are different. The polar field was about 25 percent stronger at solar minimum in 1986 than in 1976, consequently the heliospheric current sheet was confined to a narrower latitude range. In the rising phase of Cycle 21 a strong 4-sector pattern developed; the same pattern was not present in Cycle 22. Other differences include the timing of the polar field reversals and changes in latitudinal extent of the current sheet.

Hoeksema, J. T.↗

Symbolic generation of elastic rotor blade equations using a FORTRAN processor and numerical study on dynamic inflow effects on the stability of helicopter rotors

The process of performing an automated stability analysis for an elastic-bladed helicopter rotor is discussed. A symbolic manipulation program, written in FORTRAN, is used to aid in the derivation of the governing equations of motion for the rotor. The blades undergo coupled bending and torsional deformations. Two-dimensional quasi-steady aerodynamics below stall are used. Although reversed flow effects are neglected, unsteady effects, modeled as dynamic inflow are included. Using a Lagrangian approach, the governing equations are derived in generalized coordinates using the symbolic program. The program generates the steady and perturbed equations and writes into subroutines to be called by numerical routines. The symbolic program can operate on both expressions and matrices. For the case of hovering flight, the blade and dynamic inflow equations are converted to equations in a multiblade coordinate system by rearranging the coefficients of the equations. For the case of forward flight, the multiblade equations are obtained through the symbolic program. The final multiblade equations are capable of accommodating any number of elastic blade modes. The computer implementation of this procedure consists of three stages: (1) the symbolic derivation of equations; (2) the coding of the equations into subroutines; and (3) the numerical study after identifying mass, damping, and stiffness coefficients. Damping results are presented in hover and in forward flight with and without dynamic inflow effects for various rotor blade models, including rigid blade lag-flap, elastic flap-lag, flap-lag-torsion, and quasi-static torsion. Results from dynamic inflow effects which are obtained from a lift deficiency function for a quasi-static inflow model in hover are also presented.

Reddy, T. S. R.↗

Advances in Software Tools for Pre-processing and Post-processing of Overset Grid Computations

Recent developments in three pieces of software for performing pre-processing and post-processing work on numerical computations using overset grids are presented. The first is the OVERGRID graphical interface which provides a unified environment for the visualization, manipulation, generation and diagnostics of geometry and grids. Modules are also available for automatic boundary conditions detection, flow solver input preparation, multiple component dynamics input preparation and dynamics animation, simple solution viewing for moving components, and debris trajectory analysis input preparation. The second is a grid generation script library that enables rapid creation of grid generation scripts. A sample of recent applications will be described. The third is the OVERPLOT graphical interface for displaying and analyzing history files generated by the flow solver. Data displayed include residuals, component forces and moments, number of supersonic and reverse flow points, and various dynamics parameters.

Chan, William M.↗

Online model-based diagnosis to support autonomous operation of an advanced life support system

This article describes methods for online model-based diagnosis of subsystems of the advanced life support system (ALS). The diagnosis methodology is tailored to detect, isolate, and identify faults in components of the system quickly so that fault-adaptive control techniques can be applied to maintain system operation without interruption. We describe the components of our hybrid modeling scheme and the diagnosis methodology, and then demonstrate the effectiveness of this methodology by building a detailed model of the reverse osmosis (RO) system of the water recovery system (WRS) of the ALS. This model is validated with real data collected from an experimental testbed at NASA JSC. A number of diagnosis experiments run on simulated faulty data are presented and the results are discussed.

Non-NASA Center↗

Future Antarctic Bed Topography and Its Implications for Ice Sheet Dynamics

The Antarctic bedrock is evolving as the solid Earth responds to the past and ongoing evolution of the ice sheet. A recently improved ice loading history suggests that the Antarctic Ice Sheet (AIS) has generally been losing its mass since the Last Glacial Maximum. In a sustained warming climate, the AIS is predicted to retreat at a greater pace, primarily via melting beneath the ice shelves.We employ the glacial isostatic adjustment (GIA) capability of the Ice Sheet System Model (ISSM) to combine these past and future ice loadings and provide the new solid Earth computations for the AIS.We find that past loading is relatively less important than future loading for the evolution of the future bed topography. Our computations predict that the West Antarctic Ice Sheet (WAIS) may uplift by a few meters and a few tens of meters at years AD 2100 and 2500, respectively, and that the East Antarctic Ice Sheet is likely to remain unchanged or subside minimally except around the Amery Ice Shelf. The Amundsen Sea Sector in particular is predicted to rise at the greatest rate; one hundred years of ice evolution in this region, for example, predicts that the coastline of Pine Island Bay will approach roughly 45mmyr−1 in viscoelastic vertical motion. Of particular importance, we systematically demonstrate that the effect of a pervasive and large GIA uplift in the WAIS is generally associated with the flattening of reverse bed slope, reduction of local sea depth, and thus the extension of grounding line (GL) towards the continental shelf. Using the 3-D higher-order ice flow capability of ISSM, such a migration of GL is shown to inhibit the ice flow. This negative feedback between the ice sheet and the solid Earth may promote stability in marine portions of the ice sheet in the future.

ice dynamics↗

Comparisons of thermochemical nonequilibrium viscous flowfield predictions for AFE vehicle

Numerical solutions for the Aeroassist Flight Experiment vehicle were obtained from three methods at a trajectory point corresponding to the maximum aerodynamic heating. The flow regime and vehicle's speed require a viscous model of laminar flow and finite-rate chemical and thermal modeling. The computational domain covers both forebody and base such that the shock layer and near wake flowfield are included. Because of differences in computational grids, methods of solution, and models of rate equations, the results are generally in poor agreement. Temperature and species concentrations are strongly affected by the physical model equations and associated parameters. The chemistry model based on 11 species is found to yield lower translational temperature profile near stagnation than those from a seven-species model. The vibrational temperature varies according to the modeling details. All solutions indicate that strong neutral and/or molecular dissociation and weak ionization take place at the forebody and vibrational freezing is present in the afterbody expansion region where the vibrational temperature is higher than the translational temperature. Some forms of shear layer emanating from the aerobrake skirt coalesce in the region of reversed flow behind the vehicle.

Tam, Luen T.↗

Solution of mixed convection heat transfer from isothermal in-line fins

Transient and steady state combined natural and forced convective flows over two in-line finite thickness fins (louvers) in a vertical channel are numerically solved using two methods. The first method of solution is based on the 'Simple Arbitrary Lagrangian Eulerian' (SALE) technique which incorporates mainly two computational phases: (1) a Lagrangian phase in which the velocity field is updated by the effects of all forces, and (2) an Eulerian phase that executes all advective fluxes of mass, momentum and energy. The second method of solution uses the finite element code entitled FIDAP. In the first part, comparison of the results by FIDAP, SALE, and available experimental work were done and discussed for steady state forced convection over louvered fins. Good agreements were deduced between the three sets of results especially for the flow over a single fin. In the second part and in the absence of experimental literature, the numerical predictions were extended to the transient transports and to the opposing flow where pressure drop is reversed. Results are presented and discussed for heat transfer and pressure drop in assisting and opposing mixed convection flows.

Khalilollahi, Amir↗

Configuration and Intraseasonal Duration of Interannual Anomalies of the Great Plains Low-Level Jet

Despite the fact that the low-level jet of the southern Great Plains (the GPLLJ) of the U.S. is primarily a nocturnal phenomenon that virtually vanishes during the daylight hours, it is one of the most persistent and stable climatological features of the low-level continental flow during the warm-season months, May through August. We have used significant-level data to validate the skill of the GEOS-1 Data Assimilation System (DAS) in realistically detecting this jet and inferring its structure and evolution. We have then carried out a 15-year reanalysis with the GEOS-1 DAS to determine its climatology and mean diurnal cycle and to study its interannual variability. Interannual anomalies of the meridional flow associated with the GPLLJ are much smaller than the mean diurnal fluctuations, than random intraseasonal anomalies, and than the mean wind itself. There are three maxima of low-level meridional flow variance over the Great Plains and the Gulf of Mexico: a 1.2 m2 s-2 peak over the southeast Texas, to the east and south of the mean velocity peak, a 1.0 m2 s-2 peak over the western Gulf of Mexico, and a .8 m2 s-2 peak over the upper Great Plains (UGP), near the Nebraska/South Dakota border. Each of the three variance maxima corresponds to a spatially coherent, jet-like pattern of low-level flow interannual variability. There are also three dominant modes of interannual variability corresponding to the three variance maxima, but not in a simple one-to-one relationship. Cross-sectional profiles of mean southerly wind over Texas remain relatively stable and recognizable from year to year with only its eastward flank showing significant variability. This variability, however, exhibits a distinct, biennial oscillation during the first six to seven years of the reanalysis period and only then. This intermittent biennial oscillation (IBO, one of the three modes discussed in the previous paragraph) in the lowlevel flow is restricted to the region surrounding eastern Texas and is also evident in the NCEP/NCAR reanalysis data set from about 1978 to 1985 or 1986 and again from 1995 to 2000. It is evident as well in surface pressure in both the GEOS-1 and NCEP/NCAR sets. The interannual anomalies do not necessarily persist uniformly throughout an entire season, but can fluctuate from one part of the season to the next. To estimate the characteristic sub-seasonal time scales for coherence of these fluctuations, we have taken the weekly anomaly of low-level wind at each point of the domain from the climatological average for that given point and that given week of the season and computed the covariance of its fluctuations over all weeks and over all years with the weekly climatological anomaly of the meridional wind at each of the three reference points discussed above. The typical duration of a coherent interannual anomaly within a given warm season increases with decreasing latitude from 2 to 3 weeks over the UGP, to 6 to 7 weeks over eastern Texas. Coherence over the western Gulf of Mexico is intermediate between the two with a typical duration of 4 to 5 weeks. There appears to be evidence that the interannual anomalies over Texas the Gulf propagate to the UGP after a week and those over the Gulf propagate there after 2 to 3 weeks. There also appears to be some reverse propagation of interannual anomalies over the UGP to Texas and to the Gulf after a period of about one week. The interannual anomalies in southerly flow over eastern Texas seem to correlate well with interannual anomalies of surface temperature and (negative) ground wetness and over western Texas.

Helfand, H. M.↗

Analysis of plasma measurements for the Geotail mission

The CPI plasma measurements from the Geotail spacecraft are currently used by a number of scientists in support of varied research projects. Measurements from the CPI hot plasma analyzer have been processed to compute one-minute averages of plasma densities, temperatures, and velocities for the period 1 Oct. 1992 - 31 July 1993. These parameters are used in a preliminary survey of the magnetotail for distances from earth of 10 to 210 earth radii. The cold plasmas within the magnetotail drift towards the midplane and are thought to be a principal source for the hot plasma sheet. A remarkable result from Geotail is the observation of cold ion beams coexisting as distinct components in the presence of hot plasma-sheet plasmas. The development and evolution of plasmoids is a topic of considerable interest in studies of the magnetotail and magnetospheric substorms. A number of possible plasmoids have been identified in the Geotail data set based upon reversals of the Z component of the magnetic field as observed with the MGF instrumentation. An intensive study of the electron and ion velocity distributions and the plasma parameters for several of these events reveals unexpected features. Collaborative work with several Geotail instrument teams and other researchers is also ongoing.

Frank, Louis A.↗

Diffusive modification of primary zoning in olivine cores in Angrite LEW 87051

Angrite LEW 87051 consists of large olivine crystals set in a fine-grained groundmass that clearly represents a crystallized melt. A few olivines contain Ca-poor, Cr-rich cores that crystallized from a very different melt than the outer part of the crystals constituting the majority of olivine in LEW 87051. We evaluate a model in which the cores formed through fractional crystallization of one melt, then were incorporated into a different melt as xenocrysts, whereupon the original zoning patterns were modified by diffusion. Using a similar approach, we calculate zoning patterns for the cores that would result from perfect fractional crystallization, compare them with the observed zoning, and determine whether the differences could result from diffusive modification consistent with known diffusion rates for Ca, Mn, and Cr. Using distribution coefficients from the 1400 C, IW + 1 experiments, we computed CaO, Cr2O3, and MnO abundances in the hypothetical parent melt by inverting the olivine at the centers of the cores. We further assumed that the primary zoning profile for CaO is essentially unmodified, because the diffusion rate of Ca in olivine is slow. We carried out the fractional crystallization calculation until the calculated Ca content was that observed at break in zoning profiles at the outer edge of the core. We then normalized the distance of this calculated profile to the length of the observed profile in the olivine core and calculated profiles for MnO and Cr2O3. The CaO zoning profile agrees well with the observed profile. The observed MnO profile is slightly higher than the calculated profile near the edge. Diffusion calculations indicate that reversal of the general trend of primary zoning through diffusion would require that diffusion of Cr is 5-10x faster than that of Mn.

Mikouchi, T.↗