Search NASA⌕ Search

SEARCH · Search NASA

Results for “Higher order calculations”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Ion abundances and implications for photochemistry in Comets Halley (1986 III) and Bradfield (1987 XXIX)

The Ohio State University Image Dissector Scanner on the Perkins 1.8-m telescope at the Lowell Observatory was used to record spectra of the plasma tails of Comets P/Halley (1986 III) and P/Bradfield (1987 XXIX, also 1987s). The ionic species CO(+), N2(+), CH(+), and H2O(+) were identified in these spectra, and column densities for them were calculated from measured fluxes. The observed N2(+)/H2O(+) ratios are at least an order of magnitude lower and the observed CH(+)/H2O(+) ratios are a factor of 100 higher than theoretical results. The abundance ratio N2/CO was derived in the plasma tail of Bradfield from N2(+) and CO(+) data, and found to be an order of magnitude higher than a value measured in Halley. The relative ion abundances of CH(+), N2(+), CO(+), and H2O(+) in Halley are consistent with in situ measurements obtained from the Giotto spacecraft.

Lutz, Barry L.↗
NASA NTRS2019

Photoionization Calculations of the Radiation Force Due To Spectral Lines in AGNs

One of the main mechanisms that could drive mass outflows in active galactic nuclei (AGNs) is radiation pressure due to spectral lines. Although straightforward to understand, the actual magnitude of the radiation force is challenging to compute because the force depends on the physical conditions in the gas, as well as the strength, spectral energy distribution (SED), and geometry of the radiation field. We present results from our photoionization and radiation transfer calculations of the force multiplier, M(ξ, t), using the same radiation field to compute the gas photoionization and thermal balance. We assume low gas density (n = 104 per cu.cm) and column density (N ≤ 1017 per sq.cm), a Boltzmann distribution for the level populations, and the Sobolev approximation. Here, we describe results for two SEDs corresponding to an unobscured and obscured AGN in NGC 5548. Our main results are the following: (1) although M(ξ, t) starts to decrease with ξ for ξ gsim 1 as shown by others, this decrease in our calculations is relatively gradual and could be nonmonotonic as M(ξ, t) can increase by a factor of a few for ξ ≈ 10–1000; (2) at these same ξ for which the multiplier is higher than in previous calculations, the gas is thermally unstable by the isobaric criterion; (3) non-local thermodynamic equilibrium effects reduce M(t, ξ) by over two orders of magnitude for ξ gsim 100. The dynamical consequence of result (1) is that line driving can be important for ξ as high as 1000 when the LTE approximation holds, while result (2) provides a natural cloud formation mechanism that may account for the existence of narrow line regions. Result (3) suggests that line driving may not be important for ξ gsim 100 in tenuous plasma.

Dannen, Randall C.↗
NASA NTRS1996

Ring Current Development During Storm Main Phase

The development of the ring current ions in the inner magnetosphere during the main phase of a magnetic storm is studied. The temporal and spatial evolution of the ion phase space densities in a dipole field are calculated using a three dimensional ring current model, considering charge exchange and Coulomb losses along drift paths. The simulation starts with a quiet time distribution. The model is tested by comparing calculated ion fluxes with Active Magnetospheric Particle Tracer Explorers/CCE measurement during the storm main phase on May 2, 1986. Most of the calculated omnidirectional fluxes are in good agreement with the data except on the dayside inner edge (L less than 2.5) of the ring current, where the ion fluxes are underestimated. The model also reproduces the measured pitch angle distributions of ions with energies below 10 keV. At higher energy, an additional diffusion in pitch angle is necessary in order to fit the data. The role of the induced electric field on the ring current dynamics is also examined by simulating a series of substorm activities represented by stretching and collapsing the magnetic field lines. In response to the impulsively changing fields, the calculated ion energy content fluctuates about a mean value that grows steadily with the enhanced quiescent field.

Fok, Mei-Ching↗
NASA NTRS1977

Interactive calculation procedure for supersonic flows

An interactive procedure was developed for supersonic viscous flows that can be used for either two-dimensional or axisymmetric configurations. The procedure is directed to supersonic internal flows as well as those supersonic external flows that require consideration of mutual interaction between the outer flow and the boundary layer flow. The flow field is divided into two regions: an inner region which is highly viscous and mostly subsonic and an outer region where the flow is supersonic and in which viscous effects are small but not negligible. For the outer region a numerical solution is obtained by applying the method of characteristics to a system of equations which includes viscous and conduction transport terms only normal to the streamlines. The inner region is treated by a system of equations of the boundary layer type that includes higher order effects such as longitudinal and transverse curvature and normal pressure gradients. These equations are coupled and solved simultaneously in the physical coordinates by using an implicit finite difference scheme. This system can also be used to calculate laminar and turbulent boundary layers using a scalar eddy viscosity concept.

Tassa, Y.↗
NASA NTRS1976

New potentials for conventional aircraft when powered by hydrogen-enriched gasoline

Hydrogen enrichment for aircraft piston engines is under study in a new NASA program. The objective of the program is to determine the feasibility of inflight injection of hydrogen in general aviation aircraft engines to reduce fuel consumption and to lower emission levels. A catalytic hydrogen generator will be incorporated as part of the air induction system of a Lycoming turbocharged engine and will generate hydrogen by breaking down small amounts of the aviation gasoline used in the normal propulsion system. This hydrogen will then be mixed with gasoline and compressed air from the turbocharger before entering the engine combustion chamber. The special properties of the hydrogen-enriched gasoline allow the engine to operate at ultralean fuel/air ratios, resulting in higher efficiencies and hence less fuel consumption. This paper summarizes the results of a systems analysis study. Calculations assuming a Beech Duke aircraft indicate that fuel savings on the order of 20% are possible. An estimate of the potential for the utilization of hydrogen enrichment to control exhaust emissions indicates that it may be possible to meet the 1979 Federal emission standards.

Menard, W. A.↗
NASA NTRS1984

An experimental study of a supercritical trailing-edge flow

An experimental study has been conducted of a transonic, turbulent, high-Reynolds-number blunt trailing-edge flow. The model shape and the surface pressure distribution are characteristics of a modern supercritical airfoil under shock-free conditions. Reynolds number and pressure gradient scaling of the boundary layer are relevant to airfoil applications. The data set is exceptionally accurate and consistent, with the momentum balance accounting for the flux of momentum to within 1 percent, except in the immediate vicinity of the blunt trailing edge. The experimental flow exhibits strong viscous-inviscid interaction and higher-order boundary-layer effects including strong adverse streamwise pressure gradient, significant normal pressure gradients associated with surface and streamline curvature, and significant wake curvature. Navier-Stokes calculations with a two-equation K-epsilon turbulence model predict the correct pressure distribution which demonstrates the utility of these engineering tools. The experiment approaches separation at the strailing edge. However, in comparison to the experiment, the calculations predict too high skin friction and insufficient displacement thickness growth. An analysis of the turbulent and mean flow fields reveals the turbulence model defects are likely in modeling the dissipation source and sink terms, and in the eddy viscosity relation.

Brown, J. L.↗
NASA NTRS1988

The photodissociation of CO in circumstellar envelopes

The CO photodissociation rate for the unshielded ISM is calculated using recent laboratory results which confirm that photodissociation occurs by way of line absorption. A value of 2.0 x 10 to the -10th/s, an order of magnitude higher than the rate used in the past, is obtained. The new rate and a treatment of the radiative transfer and shielding are used to develop a theory for the CO abundance in the circumstellar envelopes of cool, evolved stars, and results are presented on the spatial variation of CO, C, and C(+). It is shown that these distributions play important roles in determining the observational properties of circumstellar envelopes.

Mamon, G. A.↗
NASA NTRS2013

Studies of Inviscid Flux Schemes for Acoustics and Turbulence Problems

Five different central difference schemes, based on a conservative differencing form of the Kennedy and Gruber skew-symmetric scheme, were compared with six different upwind schemes based on primitive variable reconstruction and the Roe flux. These eleven schemes were tested on a one-dimensional acoustic standing wave problem, the Taylor-Green vortex problem and a turbulent channel flow problem. The central schemes were generally very accurate and stable, provided the grid stretching rate was kept below 10%. As near-DNS grid resolutions, the results were comparable to reference DNS calculations. At coarser grid resolutions, the need for an LES SGS model became apparent. There was a noticeable improvement moving from CD-2 to CD-4, and higher-order schemes appear to yield clear benefits on coarser grids. The UB-7 and CU-5 upwind schemes also performed very well at near-DNS grid resolutions. The UB-5 upwind scheme does not do as well, but does appear to be suitable for well-resolved DNS. The UF-2 and UB-3 upwind schemes, which have significant dissipation over a wide spectral range, appear to be poorly suited for DNS or LES.

Morris, Chris↗
NASA NTRS2018

High Fidelity Collision Probabilities Estimated Using Brute Force Monte Carlo Simulations

The NASA Conjunction Assessment Risk Analysis team has implemented new software to estimate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm employs a brute force Monte Carlo (BFMC) method which differs from most other methods because it uses orbital states and covariances propagated from their orbit determination epoch times using the full set of Astrodynamics Support Workstation higher order theory models, including the High Accuracy Satellite Drag Model. This paper de-scribes the BFMC algorithm, presents comparisons of BFMC P (sub c) estimates to those calculated using other methods, and discusses the implications for conjunction risk assessment.

Probabilities↗
NASA NTRS1961

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗
NASA NTRS1999

High Order Approximations for Compressible Fluid Dynamics on Unstructured and Cartesian Meshes

The development of high-order accurate numerical discretization techniques for irregular domains and meshes is often cited as one of the remaining challenges facing the field of computational fluid dynamics. In structural mechanics, the advantages of high-order finite element approximation are widely recognized. This is especially true when high-order element approximation is combined with element refinement (h-p refinement). In computational fluid dynamics, high-order discretization methods are infrequently used in the computation of compressible fluid flow. The hyperbolic nature of the governing equations and the presence of solution discontinuities makes high-order accuracy difficult to achieve. Consequently, second-order accurate methods are still predominately used in industrial applications even though evidence suggests that high-order methods may offer a way to significantly improve the resolution and accuracy for these calculations. To address this important topic, a special course was jointly organized by the Applied Vehicle Technology Panel of NATO's Research and Technology Organization (RTO), the von Karman Institute for Fluid Dynamics, and the Numerical Aerospace Simulation Division at the NASA Ames Research Center. The NATO RTO sponsored course entitled "Higher Order Discretization Methods in Computational Fluid Dynamics" was held September 14-18, 1998 at the von Karman Institute for Fluid Dynamics in Belgium and September 21-25, 1998 at the NASA Ames Research Center in the United States. During this special course, lecturers from Europe and the United States gave a series of comprehensive lectures on advanced topics related to the high-order numerical discretization of partial differential equations with primary emphasis given to computational fluid dynamics (CFD). Additional consideration was given to topics in computational physics such as the high-order discretization of the Hamilton-Jacobi, Helmholtz, and elasticity equations. This volume consists of five articles prepared by the special course lecturers. These articles should be of particular relevance to those readers with an interest in numerical discretization techniques which generalize to very high-order accuracy. The articles of Professors Abgrall and Shu consider the mathematical formulation of high-order accurate finite volume schemes utilizing essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) reconstruction together with upwind flux evaluation. These formulations are particularly effective in computing numerical solutions of conservation laws containing solution discontinuities. Careful attention is given by the authors to implementational issues and techniques for improving the overall efficiency of these methods. The article of Professor Cockburn discusses the discontinuous Galerkin finite element method. This method naturally extends to high-order accuracy and has an interpretation as a finite volume method. Cockburn addresses two important issues associated with the discontinuous Galerkin method: controlling spurious extrema near solution discontinuities via "limiting" and the extension to second order advective-diffusive equations (joint work with Shu). The articles of Dr. Henderson and Professor Schwab consider the mathematical formulation and implementation of the h-p finite element methods using hierarchical basis functions and adaptive mesh refinement. These methods are particularly useful in computing high-order accurate solutions containing perturbative layers and corner singularities. Additional flexibility is obtained using a mortar FEM technique whereby nonconforming elements are interfaced together. Numerous examples are given by Henderson applying the h-p FEM method to the simulation of turbulence and turbulence transition.

Barth, Timothy↗
NASA NTRS2010

Uncertainties Associated with Theoretically Calculated N2-Broadened Half-Widths of H2O Lines

With different choices of the cut-offs used in theoretical calculations, we have carried out extensive numerical calculations of the N2-broadend Lorentzian half-widths of the H2O lines using the modified Robert-Bonamy formalism. Based on these results, we are able to thoroughly check for convergence. We find that, with the low-order cut-offs commonly used in the literature, one is able to obtain converged values only for lines with large half-widths. Conversely, for lines with small half-widths, much higher cut-offs are necessary to guarantee convergence. We also analyse the uncertainties associated with calculated half-widths, and these are correlated as above. In general, the smaller the half-widths, the poorer the convergence and the larger the uncertainty associated with them. For convenience, one can divide all H2O lines into three categories, large, intermediate, and small, according to their half-width values. One can use this division to judge whether the calculated half-widths are converged or not, based on the cut-offs used, and also to estimate how large their uncertainties are. We conclude that with the current Robert- Bonamy formalism, for lines in category lone can achieve the accuracy requirement set by HITRAN, whereas for lines in category 3, it 'is impossible to meet this goal.

Ma, Q.↗
NASA NTRS1976

Power conservation for reflector antennas with truncated feed patterns

Reflector-antenna calculations using idealised truncated feed patterns can lead to incorrect values of total secondary radiated power. It is demonstrated that this discrepancy is due to the presence of higher-order spherical modes incident on the reflector. When a proper spherical-wave expansion of the incident field is used, in conjunction with the physical-optics technique, to determine the scattered field, the total power of the scattered field will equal the power radiated by the feed.

Rush, W. V. T.↗
NASA NTRS1993

Noise levels and their effects on Shuttle crewmembers' performance: Operational concerns

When excessive, noise can result in sleep interference, fatigue, interference with verbal communication, and hearing damage. Shuttle crewmembers are exposed to noise throughout their mission. The contribution of noise to decrements in crew performance over these extended exposure durations was the focus of this study. On the STS-40/SLS-1, mission noise levels were evaluated through the use of a sound level meter and a crew questionnaire. Crewmembers noted that sleep, concentration, and relaxation were negatively impacted by high noise levels. Speech Interference Levels (SIL's), calculated from the sound level measurements, suggested that crewmembers were required to raise their voice in order to be heard. No difficulty detecting caution and warning alarms was noted. The higher than desirable noise levels in Spacelab were attributed to flight specific payloads for which acoustic waivers were granted. It is recommended that current noise levels be reduced in Spacelab and the Orbiter Middeck especially as longer missions are planned for the buildup of Space Station Freedom. Levels of NC 50 are recommended in areas where speech communication is required and NC 40 in sleep areas. These levels are in accordance with the NASA Man-Systems Integration Standards. Measurements proposed for subsequent orbiter missions are discussed.

Koros, Anton S.↗
NASA NTRS1995

Latitudinal distribution of black carbon soot in the upper troposphere and lower stratosphere

Black carbon soot from the upper troposphere and lower stratosphere has been systematically collected at latitudes from 90 deg N to 45 deg S. The measured latitudinal distribution of this soot at 10 to 11 km altitude is found to covary with commercial air traffic fuel use, suggesting that aircraft fuel combustion at altitude is the principal source. In addition, at latitudes where the commercial air traffic is high, measured black carbon soot values are high even at 20 km altitude, suggesting that aircraft-generated soot injected just above the tropopause may be transported to higher altitudes. During the volcanically influenced period in which these samples were collected, the number abundances, total mass, and calculated total surface area of black carbon soot are 2-3 orders of magnitude lower than similar measures of sulfuric acid aerosol. During volcanically quiescent periods, the calculated total surface area of black carbon soot aerosol is of the same order of magnitude as that of the background sulfuric acid aerosol. It appears from this comparison that black carbon soot is only capable of influencing lower stratosphere or upper troposphere chemistry during periods when the aerosol budget is not dominated by volcanic activity. It remains to determine the extent to which black carbon soot particles act as nuclei for sulfuric acid aerosol formation. However, mass balance calculations suggest that aircraft soot injected at altitude does not represent a significant source of condensation nuclei for sulfuric acid aerosols.

Blake, David F.↗
NASA NTRS2018

Analysis of the Honeywell Uncertified Research Engine (HURE) with Ice Crystal Cloud Ingestion at Simulated Altitudes: Public Version

The Honeywell Uncertified Research Engine (HURE), a research version of a turbofan engine that never entered production, was tested in the NASA Propulsion System Laboratory (PSL), an altitude test facility at the NASA Glenn Research Center. The PSL is a facility that is equipped with water spray bars capable of producing an ice cloud consisting of ice particles, having a controlled particle diameter and concentration in the air flow. In preparation for testing of the HURE, numerical analysis of flow and ice particle thermodynamics was performed on the compression system of the turbofan engine to predict operating conditions that could potentially result in a risk of ice accretion due to ice crystal ingestion. The results of those analyses formed the basis of the test matrix. The goal of the test matrix was to have ice accrete in two regions of the compression system: region one, which consists of the fan-stator through the inlet guide vane (IGV), and region two which is the first stator within the high pressure compressor. The predictive analyses were performed with the mean line compressor flow modeling code (COMDES-MELT) which includes an ice particle model. Together these comprise a one-dimensional icing tool. The HURE engine was tested in PSL with the ice cloud over the range of operating conditions of altitude, ambient temperature, simulated flight Mach number, and fan speed with guidance from the analytical predictions. The engine was fitted with video cameras at strategic locations within the engine compression system flow path where ice was predicted to accrete, in order to visually confirm ice accretion when it occurred. In addition, traditional compressor instrumentation such as total pressure and temperature probes, static pressure taps, and metal temperature thermocouples were installed in targeted areas where the risk of ice accretion was expected. The current research focuses on the analysis of the data that was obtained after testing the HURE engine in PSL with ice crystal ingestion. The computational method was enhanced by computing key parameters through the fan-stator at multiple spanwise locations, in order to increase the fidelity with the current mean-line method. In addition, other sources of heat (non-adiabatic walls) were suspected to be the cause of accretion near the splitter-lip and shroud. Since there were no thermocouples near the splitter, a simple order of magnitude heat transfer model was implemented to estimate the wall temperature. Future analyses will require a higher fidelity thermal analysis of the compression system metal walls to accurately calculate the total heat flux to the ice particle. For many data points analyzed, there were differences between the thermodynamic system model and the measured test data that may partially be responsible for uncertainties with the results of the current analyses.

Turbomachinery↗
NASA NTRS1990

Experimental and analytical evaluation of efficiency of helicopter planetary stage

The efficiency of a helicopter transmission planetary stage was studied both experimentally and analytically. Experiments were done by using a back-to-back, test-and-slave arrangement. The experiments were a parametric study of the effects of operating conditions on efficiency. In order to enhance the analysis, a model was developed that calculates the power required for the meshing gears to displace oil trapped between the gear teeth. In general, the analysis predicted higher efficiencies than were measured. The results of this study were compared with those of other studies.

Krantz, Timothy L.↗
NASA NTRS1991

The propagation characteristics of the plate modes of acoustic emission waves in thin aluminum plates and thin graphite/epoxy composite plates and tubes

Acoustic emission was interpreted as modes of vibration in plates. Classical plate theory was used to predict dispersion curves for the two fundamental modes and to calculate the shapes of flexural waveforms produced by vertical step function loading. There was good agreement between theoretical and experimental results for aluminum. Composite materials required the use of a higher order plate theory (Reissner-Mindlin) to get good agreement with the measured velocities. Four composite plates with different laminate stacking sequences were studied. The dispersion curves were determined from phase spectra of the time dependent waveforms. Plate modes were shown to be useful for determining the direction of source motion. Aluminum plates were loaded by breaking a pencil lead against their surface. By machining slots at angles to the plane of a plate, the direction in which the force acted was varied. Changing the source motion direction produced regular variations in the waveforms. To demonstrate applicability beyond simple plates, waveforms produced by lead breaks on a thin walled composite tube were also shown to be interpretable as plate modes. The tube design was based on the type of struts proposed for Space Station Freedom's trussed structures.

Prosser, William H.↗