SEARCH · Engineering Papers

Results for “numerical relativity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1991

A new polarized atmospheric radiative transfer model

A plane-parallel polarized radiative transfer model is described. The model is used to compute the radiance exiting a vertically inhomogeneous atmosphere containing randomly-oriented particles. Both solar and thermal sources of radiation are considered. A direct method of incorporating the polarized scattering information is combined with the doubling and adding method to produce a relatively simple formulation. Several numerical results are presented for verification and comparison.

Evans, K. F.↗

NASA NTRS1991

Some characteristics of small-scale turbulence in a turbulent duct flow

The fine-scale structure of turbulence in a fully developed turbulent duct flow is examined by considering the 3D velocity derivative field obtained from direct numerical simulations at two relatively small Reynolds numbers. The magnitudes of all mean-square derivatives (normalized by wall variables) increase with the Reynolds number, the increase being largest at the wall. These magnitudes are not consistent with the assumption of local isotropy except perhaps near the duct center-line. When the assumption of local isotropy is relaxed to one of local axisymmetry, or invariance with respect to rotation about a coordinate axis (here chosen in the streamwise direction), satisfactory agreement is indicated by the data outside the wall region. Support for axisymmetry is demonstrated by anisotropy invariant maps of the dissipation and vorticity tensors.

Antonia, R. A.↗

NASA NTRS1993

NiCd cell reliability in the mission environment

This paper summarizes an effort by Gates Aerospace Batteries (GAB) and the Reliability Analysis Center (RAC) to analyze survivability data for both General Electric and GAB NiCd cells utilized in various spacecraft. For simplicity sake, all mission environments are described as either low Earth orbital (LEO) or geosynchronous Earth orbit (GEO). 'Extreme value statistical methods' are applied to this database because of the longevity of the numerous missions while encountering relatively few failures. Every attempt was made to include all known instances of cell-induced-failures of the battery and to exclude battery-induced-failures of the cell. While this distinction may be somewhat limited due to availability of in-flight data, we have accepted the learned opinion of the specific customer contacts to ensure integrity of the common databases. This paper advances the preliminary analysis reported upon at the 1991 NASA Battery Workshop. That prior analysis was concerned with an estimated 278 million cell-hours of operation encompassing 183 satellites. The paper also cited 'no reported failures to date.' This analysis reports on 428 million cell hours of operation emcompassing 212 satellites. This analysis also reports on seven 'cell-induced-failures.'

Denson, William K.↗

NASA NTRS1992

Structural relations for time-dependent intermediate shocks

A quantitative description of time-dependent (MHD) intermediate shocks is formulated. In non-coplanar Riemann problems, time-dependent 2 yields 3 intermediate shocks evolve in time as a localized self-similar structure whose strength decreases as 1/sq rt t, and whose width expands as sq rt t. We derive structural relations, similar to Rankine-Hugoniot relations, between the plasma properties and the magnitude of the transverse magnetic field, which we hope can help identify time-dependent 2 yields 3 intermediate shocks at the magnetopause, in the solar wind, and elsewhere in space. The analytic forms of the structural relations agree with our numerical results.

Wu, C. C.↗

NASA NTRS1994

CFD analyses for advanced pump design

As one of the activities of the NASA/MSFC Pump Stage Technology Team, the present effort was focused on using CFD in the design and analysis of high performance rocket engine pumps. Under this effort, a three-dimensional Navier-Stokes code was used for various inducer and impeller flow field calculations. An existing algebraic grid generation procedure was-extended to allow for nonzero blade thickness, splitter blades, and hub/shroud cavities upstream or downstream of the (main) blades. This resulted in a fast, robust inducer/impeller geometry/grid generation package. Problems associated with running a compressible flow code to simulate an incompressible flow were resolved; related aspects of the numerical algorithm (viz., the matrix preconditioning, the artificial dissipation, and the treatment of low Mach number flows) were addressed. As shown by the calculations performed under the present effort, the resulting code, in conjunction with the grid generation package, is an effective tool for the rapid solution of three-dimensional viscous inducer and impeller flows.

Dejong, F. J.↗

NASA NTRS1999

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

NASA NTRS2003

The Filtered Abel Transform and Its Application in Combustion Diagnostics

Many non-intrusive combustion diagnosis methods generate line-of-sight projections of a flame field. To reconstruct the spatial field of the measured properties, these projections need to be deconvoluted. When the spatial field is axisymmetric, commonly used deconvolution method include the Abel transforms, the onion peeling method and the two-dimensional Fourier transform method and its derivatives such as the filtered back projection methods. This paper proposes a new approach for performing the Abel transform method is developed, which possesses the exactness of the Abel transform and the flexibility of incorporating various filters in the reconstruction process. The Abel transform is an exact method and the simplest among these commonly used methods. It is evinced in this paper that all the exact reconstruction methods for axisymmetric distributions must be equivalent to the Abel transform because of its uniqueness and exactness. Detailed proof is presented to show that the two dimensional Fourier methods when applied to axisymmetric cases is identical to the Abel transform. Discrepancies among various reconstruction method stem from the different approximations made to perform numerical calculations. An equation relating the spectrum of a set of projection date to that of the corresponding spatial distribution is obtained, which shows that the spectrum of the projection is equal to the Abel transform of the spectrum of the corresponding spatial distribution. From the equation, if either the projection or the distribution is bandwidth limited, the other is also bandwidth limited, and both have the same bandwidth. If the two are not bandwidth limited, the Abel transform has a bias against low wave number components in most practical cases. This explains why the Abel transform and all exact deconvolution methods are sensitive to high wave number noises. The filtered Abel transform is based on the fact that the Abel transform of filtered projection data is equal to an integral transform of the original projection data with the kernel function being the Abel transform of the filtering function. The kernel function is independent of the projection data and can be obtained separately when the filtering function is selected. Users can select the best filtering function for a particular set of experimental data. When the kernal function is obtained, it can be used repeatedly to a number of projection data sets (rovs) from the same experiment. When an entire flame image that contains a large number of projection lines needs to be processed, the new approach significantly reduces computational effort in comparison with the conventional approach in which each projection data set is deconvoluted separately. Computer codes have been developed to perform the filter Abel transform for an entire flame field. Measured soot volume fraction data of a jet diffusion flame are processed as an example.

Simons, Stephen N.↗

NASA NTRS2004

Understanding the Global Structure and Evolution of Coronal Mass Ejections in the Solar Wind

This report summarizes the technical progress made during the first six months of the second year of the NASA Living with a Star program contract Understanding the global structure and evolution of coronal mass ejections in the solar wind, between NASA and Science Applications International Corporation, and covers the period November 18, 2003 - May 17,2004. Under this contract SAIC has conducted numerical and data analysis related to fundamental issues concerning the origin, intrinsic properties, global structure, and evolution of coronal mass ejections in the solar wind. During this working period we have focused on a quantitative assessment of 5 flux rope fitting techniques. In the following sections we summarize the main aspects of this work and our proposed investigation plan for the next reporting period. Thus far, our investigation has resulted in 6 refereed scientific publications and we have presented the results at a number of scientific meetings and workshops.

Riley, Pete↗

NASA NTRS1992

An isotopic biogeochemical study of the Green River oil shale

Thirty-five different samples from three different sulfur cycles were examined in this stratigraphically oriented study of the Shell 22x-l well (U.S.G.S. C177 core) in the Piceance Basin, Colorado. Carbon isotopic compositions of constituents of Green River bitumens indicate mixing of three main components: products of primary photoautotrophs and their immediate consumers (delta approximately -30% vs PDB), products of methanotrophic bacteria (delta approximately -85%), and products of unknown bacteria (delta approximately -40%). For individual compounds synthesized by primary producers, delta-values ranged from -28 to -32%. 13C contents of individual primary products (beta-carotane, steranes, acyclic isoprenoids, tricyclic triterpenoids) were not closely correlated, suggesting diverse origins for these materials. 13C contents of numerous hopanoids were inversely related to sulfur abundance, indicating that they derived both from methanotrophs and from other bacteria, with abundances of methanotrophs depressed when sulfur was plentiful in the paleoenvironment. gamma-Cerane coeluted with 3 beta(CH3),17 alpha(H),21 beta(H)-hopane, but delta-values could be determined after deconvolution. gamma-Cerane (delta approximately -25%) probably derives from a eukaryotic heterotroph grazing on primary materials, the latter compound (delta approximately -90%) must derive from methanotrophic organisms. 13C contents of n-alkanes in bitumen differed markedly from those of paraffins generated pyrolytically. Isotopic and quantitative relationships suggest that alkanes released by pyrolysis derived from a resistant biopolymer of eukaryotic origin and that this was a dominant constituent of total organic carbon.

Non-NASA Center↗

NASA NTRS1998

Multi-Dimensional Astrophysical Structural and Dynamical Analysis: I. Development of a Nonlinear Finite Element Approach

A new field of numerical astrophysics is introduced which addresses the solution of large, multidimensional structural or slowly-evolving problems (rotating stars, interacting binaries, thick advective accretion disks, four dimensional spacetimes, etc.), as well as rapidly evolving systems.

astrophysics↗

NASA NTRS2010

Three-Body Abrasion Testing Using Lunar Dust Simulants to Evaluate Surface System Materials

Numerous unexpected operational issues relating to the abrasive nature of lunar dust, such as scratched visors and spacesuit pressure seal leaks, were encountered during the Apollo missions. To avoid reoccurrence of these unexpected detrimental equipment problems on future missions to the Moon, a series of two- and three-body abrasion tests were developed and conducted in order to begin rigorously characterizing the effect of lunar dust abrasiveness on candidate surface system materials. Two-body scratch tests were initially performed to examine fundamental interactions of a single particle on a flat surface. These simple and robust tests were used to establish standardized measurement techniques for quantifying controlled volumetric wear. Subsequent efforts described in the paper involved three-body abrasion testing designed to be more representative of actual lunar interactions. For these tests, a new tribotester was developed to expose samples to a variety of industrial abrasives and lunar simulants. The work discussed in this paper describes the three-body hardware setup consisting of a rotating rubber wheel that applies a load on a specimen as a loose abrasive is fed into the system. The test methodology is based on ASTM International (ASTM) B611, except it does not mix water with the abrasive. All tests were run under identical conditions. Abraded material specimens included poly(methyl methacrylate) (PMMA), hardened 1045 steel, 6061-T6 aluminum (Al) and 1018 steel. Abrasives included lunar mare simulant JSC- 1A-F (nominal size distribution), sieved JSC-1A-F (<25 m particle diameter), lunar highland simulant NU-LHT-2M, alumina (average diameter of 50 m used per ASTM G76), and silica (50/70 mesh used per ASTM G65). The measured mass loss from each specimen was converted using standard densities to determine total wear volume in cm3. Abrasion was dominated by the alumina and the simulants were only similar to the silica (i.e., sand) on the softer materials of aluminum and PMMA. The nominal JSC- 1A-F consistently showed more abrasion wear than the sieved version of the simulant. The lunar dust displayed abrasivity to all of the test materials, which are likely to be used in lunar landing equipment. Based on this test experience and pilot results obtained, recommendations are made for systematic abrasion testing of candidate materials intended for use in lunar exploration systems and in other environments with similar dust challenges.

Kobrick, Ryan L.↗

NASA NTRS2013

Enhanced Schapery Theory Software Development for Modeling Failure of Fiber-Reinforced Laminates

Progressive damage and failure analysis (PDFA) tools are needed to predict the nonlinear response of advanced fiber-reinforced composite structures. Predictive tools should incorporate the underlying physics of the damage and failure mechanisms observed in the composite, and should utilize as few input parameters as possible. The purpose of the Enhanced Schapery Theory (EST) was to create a PDFA tool that operates in conjunction with a commercially available finite element (FE) code (Abaqus). The tool captures the physics of the damage and failure mechanisms that result in the nonlinear behavior of the material, and the failure methodology employed yields numerical results that are relatively insensitive to changes in the FE mesh. The EST code is written in Fortran and compiled into a static library that is linked to Abaqus. A Fortran Abaqus UMAT material subroutine is used to facilitate the communication between Abaqus and EST. A clear distinction between damage and failure is imposed. Damage mechanisms result in pre-peak nonlinearity in the stress strain curve. Four internal state variables (ISVs) are utilized to control the damage and failure degradation. All damage is said to result from matrix microdamage, and a single ISV marks the micro-damage evolution as it is used to degrade the transverse and shear moduli of the lamina using a set of experimentally obtainable matrix microdamage functions. Three separate failure ISVs are used to incorporate failure due to fiber breakage, mode I matrix cracking, and mode II matrix cracking. Failure initiation is determined using a failure criterion, and the evolution of these ISVs is controlled by a set of traction-separation laws. The traction separation laws are postulated such that the area under the curves is equal to the fracture toughness of the material associated with the corresponding failure mechanism. A characteristic finite element length is used to transform the traction-separation laws into stress-strain laws. The ISV evolution equations are derived in a thermodynamically consistent manner by invoking the stationary principle on the total work of the system with respect to each ISV. A novel feature is the inclusion of both pre-peak damage and appropriately scaled, post-peak strain softening failure. Also, the characteristic elements used in the failure degradation scheme are calculated using the element nodal coordinates, rather than simply the square root of the area of the element.

Pineda, Evan J.↗

NASA NTRS2015

A Triboelectric Sensor Array for Electrostatic Studies on the Lunar Surface

The moons electrostatic environment requires careful consideration in the development of future lunar landers. Electrostatically charged dust was well documented during the Apollo missions to cause thermal control, mechanical, and visibility issues. The fine dust particles that make up the surface are electrostatically charged as a result of numerous charging mechanisms. The relatively dry conditions on the moon creates a prime tribocharging environment during surface operations. The photoelectric effect is dominant for lunar day static charging, while plasma electrons are the main contributor for lunar night electrostatic effects. Electrostatic charging is also dependent on solar intensity, Earth-moon relative positions, and cosmic ray flux. This leads to a very complex and dynamic electrostatic environment that must be studied for the success of long term lunar missions.In order to better understand the electrostatic environment of planetary bodies, Kennedy Space Center, in previous collaboration with the Jet Propulsion Laboratory, has developed an electrostatic sensor suite. One of the instruments included in this package is the triboelectric sensor array. It is comprised of strategically selected materials that span the triboelectric series and that also have previous spaceflight history. In this presentation, we discuss detailed testing with the triboelectric sensor array performed at Kennedy Space Center. We will discuss potential benefits and use cases of this low mass, low cost sensor package, both for science and for mission success.

Charging↗

NASA NTRS2018

Dynamics and Control of a Tethered Enhanced Gravity Tractor Performing Asteroid Deflection

The dynamics and control of an Enhanced Gravity Tractor (EGT) augmented with a tether for deflecting an asteroid are studied. A conventional EGT consists of collected asteroidal mass collocated with the spacecraft. Because of the presence of a tether, the collected mass is placed where the EGT would have been without a tether, and the spacecraft is placed farther away from the asteroid. Doing so improves the fuel efficiency and safety margin of the EGT operation without significantly sacrificing the gravitational attraction between the asteroid and the EGT. The tether is modeled as a series of particles connected by spring-dashpot systems. Physical properties of the tether are selected to be similar to those of the SPECTRA-1000, Kevlar-29, and Kevlar-49 fibers. It is assumed that control is applied only to the spacecraft, and there is no active control associated with the collected mass. A Proportional-Derivative (PD) controller is employed to maintain the spacecraft and the collected mass at desired positions relative to the asteroid. Numerical simulations of tethered EGT operations at 2008 EV5, Itokawa, Apophis, and a fictitious ellipsoidal asteroid are performed. It is demonstrated that a PD controller is capable of accomplishing the control objectives. The gravity gradient and the control force keep the tether stretched throughout a normal tethered EGT operation, and the load on the tether is well within the design limit of the tether material. While including multiple particles in the tether model is essential in capturing details of tether vibration, the number of particles does not significantly affect the motions of the collected mass and the spacecraft. In addition, the distance from the asteroid mass center to the collected mass should be chosen judiciously in the case of a rotating slender asteroid; some distance ranges should be avoided as excessive lateral oscillations can be excited by resonance between the asteroid rotation and tether pendular motion.

Shen, Haijun↗

NASA NTRS2020

Properties and Astrophysical Implications of the 150 M Binary Black Hole Merger GW190521

The gravitational-wave signal GW190521 is consistent with a binary black hole (BBH) merger source at redshift 0.8 with unusually high component masses, 85 (+21,-14) Mꙩ and 66 (+17,-18) Mꙩ, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theory, in the approximate range 65–120 Mꙩ. The probability that at least one of the black holes in GW190521 is in that range is 99.0%. The final mass of the merger (142 (+28,-16) Mꙩ) classifies it as an intermediate-mass black hole. Under the assumption of a quasi-circular BBH coalescence, we detail the physical properties of GW190521’s source binary and its post-merger remnant, including component masses and spin vectors. Three different waveform models, as well as direct comparison to numerical solutions of general relativity, yield consistent estimates of these properties. Tests of strong-field general relativity targeting the merger-ringdown stages of the coalescence indicate consistency of the observed signal with theoretical predictions. We estimate the merger rate of similar systems to be 0.13 (+0.30,-0.11)Gp/cu.c yr. We discuss the astrophysical implications of GW190521 for stellar collapse and for the possible formation of black holes in the pair-instability mass gap through various channels: via (multiple) stellar coalescences, or via hierarchical mergers of lower-mass black holes in star clusters or in active galactic nuclei. We find it to be unlikely that GW190521 is a strongly lensed signal of a lower-mass black hole binary merger. We also discuss more exotic possible sources for GW190521, including a highly eccentric black hole binary, or a primordial black hole binary.

R. Abbott↗

NASA NTRS1995

[Kelvin-Helmholtz instability in protostellar jets]

NASA grant NAG 5 2866, funded by the Astrophysics Theory Program, enabled the study the Kelvin-Helmholtz instability in protostellar jets. In collaboration with co-investigator Philip Hardee, the PI derived the analytic dispersion relation for the instability in including a cooling term in the energy equation which was modeled as one of two different power laws. Numerical solutions to this dispersion relation over a wide range of perturbation frequencies, and for a variety of parameter values characterizing the jet (such as Mach number, and density ratio) were found It was found that the growth rates and wavelengths associated with unstable roots of the dispersion relation in cooling jets are significantly different than those associated with adiabatic jets, which have been studied previously. In collaboration with graduate student Jianjun Xu (funded as a research associate under this grant), hydrodynamical simulations were used to follow the growth of the instability into the nonlinear regime. It was found that asymmetric surface waves lead to large amplitude, sinusoidal distortions of the jet, and ultimately to disruption Asymmetric body waves, on the other hand, result in the formation of shocks in the jet beam in the nonlinear regime. In cooling jets, these shocks lead to the formation of dense knots and filaments of gas within the jet. For sufficiently high perturbation frequencies, however, the jet cannot respond and it remains symmetric. Applying these results to observed systems, such as the Herbig-Haro jets HH34, HH111 and HH47 which have been observed with the Hubble Space Telescope, we predicted that some of the asymmetric structures observed in these systems could be attributed to the K-H modes, but that perturbations on timescales associated with the inner disk (about 1 year) would be too rapid to cause disruption. Moreover, it was discovered that weak shock 'spurs' in the ambient gas produced by ripples in the jet surface due to nonlinear, modes of surface and/or body waves could accelerate the ambient gas to low velocity. This latter effect represents a new mechanism by which supersonic jets can accelerate low velocity outflows.

Stone, James↗

NASA NTRS1995

Wave propagation and scattering in computational aeroacoustics

This paper describes predictions of model problems in computational aeroacoustics. Two problem classes are considered: The first (Workshop Category 3, Problem 1 and 2) considers two-dimensional wave propagation and non-reflecting boundary conditions in the presence of a mean flow. The second (Workshop Category 4, Problem 1) examines wall boundary conditions. For the last problem we introduce the Impedance Mismatch Method (IMM) to treat the solid wall boundaries. In this method the solid wall is simulated using a wall region in which the characteristic impedance is set to a different value to that in the fluid region. This method has advantages over traditional solid wall boundary conditions including simplicity of coding, speed of computations, and the ability to treat curved boundaries efficiently. Several numerical examples are given in addition to the Workshop Problems. The discretization of the Euler equations is performed in all cases with a Dispersion-Relation-Preserving (DRP) algorithm. The numerical results are compared with either analytical solutions or solutions obtained using traditional solid wall boundary conditions.

Chung, Cathy↗