SEARCH · Engineering Papers
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.
In this article, we make a few remarks on the role that attractors and inertial manifolds play in fluid mechanics problems. We then describe the role of incremental unknowns for approximating attractors and inertial manifolds when finite difference multigrid discretizations are used. The relation with direct numerical simulation and large eddy simulation is also mentioned.
Most reliability analysis techniques and tools assume that a system is used for a mission consisting of a single phase. However, multiple phases are natural in many missions. The failure rates of components, system configuration, and success criteria may vary from phase to phase. In addition, the duration of a phase may be deterministic or random. Recently, several researchers have addressed the problem of reliability analysis of such systems using a variety of methods. A new technique for phased-mission system reliability analysis based on Boolean algebraic methods is described. Our technique is computationally efficient and is applicable to a large class of systems for which the failure criterion in each phase can be expressed as a fault tree (or an equivalent representation). Our technique avoids state space explosion that commonly plague Markov chain-based analysis. A phase algebra to account for the effects of variable configurations and success criteria from phase to phase was developed. Our technique yields exact (as opposed to approximate) results. The use of our technique was demonstrated by means of an example and present numerical results to show the effects of mission phases on the system reliability.
We have monitored the evolution of the transient X-ray source GRO J0422+32 from approximately 2 weeks postdiscovery into its early decline phase at ultraviolet, optical, and radio wavelengths. Optical and ultraviolet spectra exhibit numerous, but relatively weak, high-excitation emission lines such as those arising from He II, N III, N V, and C IV superposed on an intrinsically blue continuum. High-resolution optical spectroscopy reveals line profiles which are double peaked, and in the case of the higher order Balmer lines, superposed on a broad absorption profile. The early outburst optical-ultraviolet continuum energy distribution is well represented by a two power-law fit with a break at approximately equal 4000 A. Radio observations with the Very Large Array (VLA) reveal a flat-spectrum source, slowly increasing in intensity at the earliest epochs observed, followed by an approximate power-law decay light curve with an index of -1. Light curves for each wavelength domain are presented and discussed. Notable are the multiple secondary outbursts seen in the optical more than 1 year postdiscovery, and spectral changes associated with secondary rises seen in the radio and UV. We find that the ultraviolet and optical characteristics of GRO J0422+32 as well as its radio evolution, are similar to other recent well-observed soft X-ray transients (also called X-ray novae) such as Cen X-4, A0620-00 (V616 Mon), and Nova Muscae 1991 (GS 1124-683), suggesting that GRO J0422+32 is also a member of that subclass of low-mass X-ray binaries. We present definitive astrometric determination of the source position, and place an upper limit of R approximately equals 20 from our analysis of the Palomar Observatory Sky Survey (POSS). Additionally, we derive distinct values for color excess from analysis of the optical (E(B-V) = 0.23) and ultraviolet (E(B-V) = 0.4) data, suggesting an intrinsic magnitude of 19-19.5 for the progenitor if it is mid-K dwarf. This leads to a likely range of 2.4-3.0 kpc for the source distance, which is consistent with our separate estimate of 2.4 +/- 0.4 kpc based on measurement of the NaD interstellar line profile. Adopting 2.4 kpc and E(B-V) = 0.23, the outburst absolute magnitude was M approximately equals 0.0, which is a typical value for this class of objects.
The problem of ion injection and acceleration at quasi perpendicular shocks has been the subject of some debate over the past two decades. It is widely known that these shocks efficiently accelerate particles that are well in the high-energy tail of the distribution. However, the issue of injection, or the acceleration of low-energy ions, has yet to reach a consensus. The fundamental issue is whether there is enough diffusion normal to the magnetic field for the particles to remain near the shock. Since transverse diffusion is a physical process that is not well understood in space plasmas, this is an important, and difficult issue to address. In this report, we will investigate the ion injection problem by performing test particle orbit integrations using synthesized turbulent fields. These fields are fully three-dimensional so that transverse diffusion is possible (cross-field diffusion is not possible in geometries where the electromagnetic fields are less than three dimensional). The synthesized fields are produced by superimposing a three-dimensional wave field on a background field. For completeness, we will compare the results from this model with the more well-established theories, such as the diffusive approximation and scatter-free shock drift acceleration. We will also compare these results with other numerical simulation techniques such as the well known hybrid simulation, and other test-particle calculations in which the shock fields are specified to have less than three dimensions. We will also discuss some recent relevant observations and how these compare with our results.
Experimental data on pressure distribution and heat transfer on a turbine airfoil were obtained over a range of Reynolds numbers from 0.75 to 7.0 x 10(exp 6) and a range of turbulence intensities from 1.8 to about 15%. The purpose of this study was to obtain fundamental heat transfer and pressure distribution data over a wide range of high Reynolds numbers and to extend the heat transfer data base to include the range or Reynolds numbers encountered in the Space Shuttle main engine turbopump turbines. The results of this study indicated that Reynolds number and turbulence intensity have a large effect on both the transition from laminar to turbulent flow and the resulting heat transfer. For a given turbulence intensity, heat transfer for all Reynolds numbers at the leading edge can generally be correlated with the Frossling number developed for lower Reynolds numbers. For a given turbulence intensity, heat transfer for the airfoil surfaces downstream of the leading edge can be approximately correlated with a dimensionless parameter. Comparisons of the experimental results were also made with a numerical solution from a two-dimensional Navier-Stokes code.
A research program is currently underway with the final objective of developing a fundamental understanding of the controlling mechanisms underlying the flammability diagrams of solid combustible materials and their derived fire properties. Given that there is a high possibility of an accidental fire occurring in a space-based facility, understanding the fire properties of materials that will be used in such facilities is of critical importance. With this purpose, the flammability diagrams of the materials, as those produced by the Lateral Ignition and Flame Spread Test (LIFT) apparatus and by a new forced flow device, the Forced Flow Ignition and Flame Spread Test (FIST) apparatus, will be obtained. The specific objective of the program is to apply the new flammability apparatus, which will more accurately reflect the potential ambient conditions of space-based environments, to the characterization of the materials for space applications. This paper presents a parametric study of oxidizer flow effects on the ignition curve of the flammability diagrams of PMMA. The dependence of the ignition delay time on the external radiant flux and either the sample width (LIFT) or the flow velocity (FIST) has been studied. Although preliminary, the results indicate that natural and forced convection flow changes, affect the characteristics of the ignition curves of the flammability diagrams. The major effect on the ignition time appears to be due to convective transfer variations at the fuel surface. At high radiant fluxes or high flow velocities, however, it appears that gas phase processes become increasingly important, affecting the overall ignition delay time. A numerical analysis of the solid fuel heating and pyrolysis has also been developed. The theoretical predictions approximate the experiments well for conditions in which the gas phase induction time is negligible.
Steady and unsteady viscous, three-dimensional flowfields are calculated using a thin layer approximation of Navier-Stokes equations in conjunction with Chimera overset grids. The finite-difference numerical scheme uses structured grids and a pentadiagonal flow solver called "OVERFLOW". The configuration of Boeing 747-200 has been chosen as one of configurations to be used as a platform for the SOFIA (Stratospheric Observatory For Infrared Astronomy). Initially, the steady flowfield of the full aircraft is calculated for the clean configuration (without a cavity to house telescope). This solution is then used to start the unsteady flowfield of a configuration containing cavity housing the observation telescope and its peripheral units. Analysis of unsteady flowfield in the cavity and its influence on the tail empennage, as well as the noise due to turbulence and optical quality of the flow are the main focus of this study. For the configuration considered here, the telescope housing cavity is located slightly downstream of the portwing. The entire flow-field is carefully constructed using 45 overset grids and consists of nearly 4 million grid points. All the computations axe done at one freestream flow condition of M(sub infinity) = 0.85, alpha = 2.5deg, and a Reynolds of Re = 1.85x10deg
The invention of oxygenic photosynthesis was likely accompanied by the introduction of large amounts of O2 and complementary reduced gases (chiefly CH4) into the atmosphere. To first approximation the venting of O2 and CH4 are stochiometrically linked. We therefore present a suite of numerical photochemical models that address the anoxic-oxic transition in an atmosphere driven by large linked inputs of biogenic 02 and CH4. We find in general that, in steady state, there are two solutions, one oxic and the other anoxic. The anoxic solution appears to be linearly stable. If volcanic SO2 fluxes are large, S disproportionates into oxidized (H2S04) and reduced (S8) exit channels. As elemental sulfur is insoluble it provides a means of preserving photochemical mass-independent fractionation (MIF). On the other hand, if the source of volcanic SO2 is smaller than today, all S can leave the atmosphere as S8. Under these conditions there would be no MIF signal. The oxic solution appears to be linearly unstable. In the oxic solutions S is invariably oxidized to sulfate, and the MIF signal would be absent. The transitional atmosphere is relatively unstable and is also the most photochemically active. Consequently it is the transitional atmosphere, not the oxic or anoxic atmospheres, that has the lowest CH4 levels and weakest greenhouse warming. As a practical matter we expect the transitional atmospheres to vary strongly in response to diurnal and seasonal biological forcing.
NASA climate reanalysis datasets from MERRA-2, distributed at the Goddard Earth Sciences Data and Information Services Center (GES DISC), have been used in broad research areas, such as climate variations, extreme weather, agriculture, renewable energy, and air quality, etc. The datasets contain numerous variables for atmosphere, land, and ocean, grouped into 95 products. The total archived volume is approximately 337 TB ( approximately 562K files) at the end of October 2017. Due to the large number of products and files, and large data volumes, it may be a challenge for a user to find and download the data of interest. The support team at GES DISC, working closely with the MERRA-2 science team, has created and is continuing to work on value added data services to best meet the needs of a broad user community. This presentation, using aerosol over Asia Monsoon as an example, provides an overview of the MERRA-2 data services at GES DISC, including: How to find the data? How many data access methods are provided? What are the best data access methods for me? How do download the subsetted (parameter, spatial, temporal) data and save in preferred spatial resolution and data format? How to visualize and explore the data online? In addition, we introduce a future online analytic tool designed for supporting application research, focusing on long-term hourly time-series data access and analysis.
A novel laboratory experimental design is described that will investigate the processing of dust grains in astrophysical shocks. Dust is a ubiquitous ingredient in the interstellar medium (ISM) of galaxies; however, its evolutionary cycle is still poorly understood. Especially shrouded in mystery is the efficiency of grain destruction by astrophysical shocks generated by expanding supernova remnants. While the evolution of these remnants is fairly well understood, the grain destruction efficiency in these shocks is largely unknown. The experiments described herein will fill this knowledge gap by studying the dust destruction efficiencies for shock velocities in the range of approximately 10-30 kilometers per second (microns per nanosecond), at which most of the grain destruction and processing in the ISM takes place. The experiments focus on the study of grain-grain collisions by accelerating small (1 millimeter) dust particles into a large (approximately 5-10 millimeter diameter) population; this simulates the astrophysical system well in that the more numerous, small grains impact and collide with the large population. Facilities that combine the versatility of high-power optical lasers with the diagnostic capabilities of X-ray free-electron lasers, e.g., the Matter in Extreme Conditions instrument at the SLAC (originally named Stanford Linear Accelerator Center) National Accelerator Laboratory, provide an ideal laboratory environment to create and diagnose dust destruction by astrophysically relevant shocks at the micron scale.
Pluto and Charon are most likely the remnants of a large number of objects that existed in the Uranus-Neptune region at early epochs of the solar system. Numerical integrations have shown that, in general, such objects were ejected from the planetary region on timescales of approximately 10(exp 7) years after Neptune and Uranus reached their current masses. It is thought that the Pluto-Charon system survived to current times without being dynamically removed in this way because it is trapped in a set of secular and mean motion resonances with Neptune. The best-known Pluto-Neptune orbit coupling is the 3:2 mean motion resonance discovered almost 30 years ago by C. Cohen and E. Hubbard. These workers showed that the resonance angle, delta is equivalent to 3(lambda(sub P)) - 2(lambda(sub N)) - omega-bar(sub P) where omega-bar(sub P) is the longitude of perihelion of the Pluto-Charon system, and lambda(sub N) and lambda(sub P) are the mean longitude of Neptune and Pluto-Charon respectively, librates about 180 deg with an amplitude, A(sub delta), of 76 deg. A numerical simulation project to map out the stability region of the 3:2 resonance is reported. The results of these simulations are important to understanding whether Pluto's long-term heliocentric stability requires only the 3:2 resonance, or whether it instead requires one or more of the other Pluto-Neptune resonances. Our study also has another important application. By investigating stability timescales as a function of orbital elements, we gain insight into the fraction of orbital phase space which the stable 3:2 resonance occupies. This fraction is directly related to the probability that the Pluto-Charon system (and possibly other small bodies) could have been captured into this resonance.
The lecture material will discuss the application of one-dimensional approximate Riemann solutions and high order accurate data reconstruction as building blocks for solving multi-dimensional hyperbolic equations. This building block procedure is well-documented in the nationally available literature. The relevant stability and convergence theory using positive operator analysis will also be presented. All participants in the minisymposium will be asked to solve one or more generic test problems so that a critical comparison of accuracy can be made among differing approaches.
The main drawback of Computational Fluid Dynamics (CFD) simulations has been the time and resource consuming nature which is not suitable for real-time applications. In this work, we first generated numerous CFD models of a given indoor space to obtain airspeed, temperature, and mean radiant temperature near an occupant as training data. Several artificial neural networks (ANN) models were trained using this CFD simulated data to approximate near real-time environmental conditions for a given occupant. This trained ANN model approach is a part of a real-time simulation of building operations using a combination of software and real hardware (HVAC equipment) approaches. The preliminary results suggest that the CFD- generated training data and the trained ANN model can accurately approximate such conditions in a real-time application, a method that has great potential in building simulation and building digital twin areas of research.
Static vibration and flutter analysis of a straight elastic axis blade was performed based on a finite element method solution. The total potential energy functional was formulated according to linear beam theory. The inertia and aerodynamic loads were formulated according to the blade absolute acceleration and absolute velocity vectors. In vibration analysis, the direction of motion of the blade during the first out-of-lane and first in-plane modes was examined; numerical results involve NASA/DOE Mod-0, McCauley propeller, north wind turbine and flat plate behavior. In flutter analysis, comparison cases were examined involving several references. Vibration analysis of a nonstraight elastic axis blade based on a finite element method solution was performed in a similar manner with the straight elastic axis blade, since it was recognized that a curved blade can be approximated by an assembly of a sufficient number of straight blade elements at different inclinations with respect to common system of axes. Numerical results involve comparison between the behavior of a straight and a curved cantilever beam during the lowest two in-plane and out-of-plane modes.
Density-reconstruction sharpens the baryon acoustic oscillations signal by undoing some of the smoothing incurred by nonlinear structure formation. Here, we present an analytical model for reconstruction based on the Zeldovich approximation, which for the first time includes a complete set of counterterms and bias terms up to quadratic order and can fit real and redshift-space data pre- and post-reconstruction data in both Fourier and configuration space over a wide range of scales. We compare our model to n-body data at $z = 0$ from the ${\tt DarkSky}$ simulation, finding sub-percent agreement in both real space and in the redshift-space power spectrum monopole out to $k = 0.4\ h$ Mpc$^{-1}$, and out to $k = 0.2\ h$ Mpc$^{-1}$ in the quadrupole, with comparable agreement in configuration space. We compare our model with several popular existing alternatives, updating existing theoretical findings for exponential damping in wiggle/no-wiggle splits of the BAO signal and discuss the usually-ignored effect of higher bias contributions on the reconstructed signal. In the appendices, we re-derive the former within our formalism, present exploratory results on higher-order corrections due to nonlinearities inherent to reconstruction, and present numerical techniques with which to calculate the redshift-space power spectrum of biased tracers within the Zeldovich approximation.
Structural Design and Analysis Toolset is a collection of approximately 26 Microsoft Excel spreadsheet programs, each of which performs calculations within a different subdiscipline of structural design and analysis. These programs present input and output data in user-friendly, menu-driven formats. Although these programs cannot solve complex cases like those treated by larger finite element codes, these programs do yield quick solutions to numerous common problems more rapidly than the finite element codes, thereby making it possible to quickly perform multiple preliminary analyses - e.g., to establish approximate limits prior to detailed analyses by the larger finite element codes. These programs perform different types of calculations, as follows: 1. determination of geometric properties for a variety of standard structural components; 2. analysis of static, vibrational, and thermal- gradient loads and deflections in certain structures (mostly beams and, in the case of thermal-gradients, mirrors); 3. kinetic energies of fans; 4. detailed analysis of stress and buckling in beams, plates, columns, and a variety of shell structures; and 5. temperature dependent properties of materials, including figures of merit that characterize strength, stiffness, and deformation response to thermal gradients
Differencing operators of arbitrarily high order can be constructed by interpolating a polynomial through a set of data followed by differentiation of this polynomial and finally evaluation of the polynomial at the point where a derivative approximation is desired. Furthermore, the interpolating polynomial can be constructed from algebraic, trigonometric, or, perhaps exponential polynomials. This paper begins with a comparison of such differencing operator construction. Next, the issue of proper grids for high order polynomials is addressed. Finally, an adaptive numerical method is introduced which adapts the numerical grid and the order of the differencing operator depending on the data. The numerical grid adaptation is performed on a Chebyshev grid. That is, at each level of refinement the grid is a Chebvshev grid and this grid is refined locally based on wavelet analysis.
The current drainage from a plasma through approximately 0.05 cm diameter holes in Kapton H polymide film, FEP type C, and quartz (7940 Corning glass) placed on a probe at voltages up to 2000 volts dc have been determined both analytically and experimentally. A Laplace field was used to numerically predict an upper limit for the drainage current. The measured current was less than the calculated current for the FEP and quartz, and the two were approximately in agreement for Kapton H. The measured currents were on the order of 5 microamperes with the calculated currents on the order of 10 microamperes.