Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “delta method”

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 469 records · Page 26

Stable isotope analysis at the molecular level: A new approach for determining the origins of amino acids in the Murchison meteorite

A combined gas chromatography/isotope ratio mass spectrometry (GC/IRMS) method has been developed that permits the direct stable carbon isotope analysis of N(O)-trifluoroacetyl-isopropyl esters of individual amino acids and their respective enantiomers at nanomole abundances. Calculation of the original delta C-13 values of the amino acids is accomplished via a correction for the carbon introduced during the derivatization process. Previous GC/IRMS analyses of individual amino acids in the non-hydrolyzed water extract of an interior sample of a Murchison meteorite stone revealed an enrichment in C-13 relative to terrestrial organic matter, in agreement with previous findings for bulk extracts. The range of amino acid delta C-13 values (+5 to +30 per mill, PDB) suggests possible kinetic effects during synthesis. In this study, an apparent kinetic isotope effect was also observed for the amino acid products of a spark discharge experiment. These preliminary resutls are supportive of a similar mechanism for the abiotic synthesis of amino acids in the Murchison meteorite.

Engel, M. H.↗

Numerical Modeling of Liquid-Vapor Phase Change

We implemented a two- and three-dimensional finite difference/front tracking technique to solve liquid-vapor phase change problems. The mathematical and the numerical features of the method were explained in great detail in our previous reports, Briefly, we used a single formula representation which incorporated jump conditions into the governing equations. The interfacial terms were distributed as singular terms using delta functions so that the governing equations would be the same as conventional conservation equations away from the interface and in the vicinity of the interface they would provide correct jump conditions. We used a fixed staggered grid to discretize these equations and an unstructured grid to explicitly track the front. While in two dimensions the front was simply a connection of small line segments, in three dimensions it was represented by a connection of small triangular elements. The equations were written in conservative forms and during the course of computations we used regriding to control the size of the elements of the unstructured grid. Moreover, we implemented a coalescence in two dimensions which allowed the merging of different fronts or two segments of the same front when they were sufficiently close. We used our code to study thermocapillary migration of bubbles, burst of bubbles at a free surface, buoyancy-driven interactions of bubbles, evaporation of drops, rapid evaporation of an interface, planar solidification of an undercooled melt, dendritic solidification, and a host of other problems cited in the reference.

Esmaeeli, Asghar↗

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'prop-burn-prop' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons are made with three methods using modified linearized covariance propagation. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. The second method does not model the maneuver but includes a process noise matrix to account for the uncertainty in its magnitude. The third method, which is essentially a hybrid of the first two, includes the nominal portion of the maneuver via the state transition matrix and uses a process noise matrix to account for the magnitude uncertainty. The first method is unable to produce the final orbit covariance except in the case of zero maneuver uncertainty. The second method yields good accuracy for the final covariance matrix but fails to model the final orbital state accurately. Agreement between the simulated covariance data produced by this method and the Monte Carlo truth data fell within 0.5-2.5 percent over a range of maneuver sizes that span two orders of magnitude (0.1-20 m/s). The third method, which yields a combination of good accuracy in the computation of the final covariance matrix and correct accounting for the presence of the maneuver in the nominal orbit, is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, PC. However, applications for the two other methods exist and are briefly discussed. Although the process model ("prop-burn-prop") that was studied is very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗

Method and apparatus for contour mapping using synthetic aperture radar

By using two SAR antennas spaced a known distance, B, and oriented at substantially the same look angle to illuminate the same target area, pixel data from the two antennas may be compared in phase to determine a difference delta phi from which a slant angle theta is determined for each pixel point from an equation Delta phi = (2 pi B/lambda)sin(theta - alpha), where lambda is the radar wavelength and alpha is the roll angle of the aircraft. The height, h, of each pixel point from the aircraft is determined from the equation h = R cos theta, and from the known altitude, a, of the aircraft above sea level, the altitude (elevation), a', of each point is determined from the difference a - h. This elevation data may be displayed with the SAR image by, for example, quantizing the elevation at increments of 100 feet starting at sea level, and color coding pixels of the same quantized elevation. The distance, d, of each pixel from the ground track of the aircraft used for the display may be determined more accurately from the equation d = R sin theta.

Goldstein, R. M.↗

Discrete-continuous variable structural synthesis using dual methods

Approximation concepts and dual methods are extended to solve structural synthesis problems involving a mix of discrete and continuous sizing type of design variables. Pure discrete and pure continuous variable problems can be handled as special cases. The basic mathematical programming statement of the structural synthesis problem is converted into a sequence of explicit approximate primal problems of separable form. These problems are solved by constructing continuous explicit dual functions, which are maximized subject to simple nonnegativity constraints on the dual variables. A newly devised gradient projection type of algorithm called DUAL 1, which includes special features for handling dual function gradient discontinuities that arise from the discrete primal variables, is used to find the solution of each dual problem. Computational implementation is accomplished by incorporating the DUAL 1 algorithm into the ACCESS 3 program as a new optimizer option. The power of the method set forth is demonstrated by presenting numerical results for several example problems, including a pure discrete variable treatment of a metallic swept wing and a mixed discrete-continuous variable solution for a thin delta wing with fiber composite skins.

Schmit, L. A.↗

The lateral-directional characteristics of a 74-degree Delta wing employing gothic planform vortex flaps

The low speed lateral/directional characteristics of a generic 74 degree delta wing body configuration employing the latest generation, gothic planform vortex flaps was determined. Longitudinal effects are also presented. The data are compared with theoretical estimates from VORSTAB, an extension of the Quasi vortex lattice Method of Lan which empirically accounts for vortex breakdown effects in the calculation of longitudinal and lateral/directional aerodynamic characteristics. It is indicated that leading edge deflections of 30 and 40 degrees reduce the magnitude of the wing effective dihedral relative to the baseline for a specified angle of attack or lift coefficient. For angles of attack greater than 15 degrees, these flap deflections reduce the configuration directional stability despite improved vertical tail effectiveness. It is shown that asymmetric leading edge deflections are inferior to conventional ailerons in generating rolling moments. VORSTAB calculations provide coarse lateral/directional estimates at low to moderate angles of attack. The theory does not account for vortex flow induced, vertical tail effects.

Grantz, A. C.↗

Viscous vortical flow calculations over delta wings

Two approaches to calculate turbulent vortical flows over delta wing configurations are illustrated. The first is for a simple delta wing at low speeds using the boundary layer approximation to treat the effects of the secondary separation. The second is for the supersonic case of a generic fighter using the NASA Ames parabolized Navier/Stokes method. Test/theory comparisons are given in both cases.

Blom, G.↗

Spectroscopic Diagnostics and an Arc Jet Heated Air Plasma

Spectral radiation measurements were made in the range of 200 to 900 nm across a section of the plenum of an arc jet wind tunnel using a series of optical fibers. The spectra contained line radiation from Oxygen and Nitrogen atoms and molecular radiation from N2(+), N2, and NO. Abel inversion technique is used to obtain radial distribution of the spectra. The analysis yielded radial profiles of the electronic excitation, vibrational and rotational temperatures of the flow field. Spectral fitting yielded branching ratios for different vibrational and rotational bands. Relatively mild flow conditions, i.e. enthalpy and mass flow rate, were used for prolonged measurements of up to and over two hours to establish the best experimental methods of temperature determinations. Signal to noise was improved by at least an order of magnitude enabling the molecular vibrational band heads of N2(+) (first negative system), N2 (second positive system), and NO (beta, gamma, delta, and epsilon systems) to be resolved in the lower ultraviolet wavelength regions. The increased signal to noise ratio also enabled partial resolution of the rotational lines of N2(+) and N2 in certain regions of minimal overlap. Comparison of the spectra with theoretical models such as the NEQAIR2 code are presented and show potential for fitting the spectra when reliable calibration is performed for the complete wavelength range.

Mack, Larry Howard, Jr.↗

A First Look at Carbon and Oxygen Stable Isotope Measurements of Martian Atmospheric C02 by the Phoenix Lander

Precise stable isotope measurements of the CO2 in the martian atmosphere have the potential to provide important constraints for our understanding of the history of volatiles, the carbon cycle, current atmospheric processes, and the degree of water/rock interaction on Mars. The isotopic composition of the martian atmosphere has been measured using a number of different methods (Table 1), however a precise value (<1%) has yet to be achieved. Given the elevated 13C values measured in carbonates in martian meteorites it has been supposed that the martian atmosphere was enriched in delta(sup 13)C. This was supported by measurements of trapped CO2 gas in EETA 79001[2] which showed elevated delta(sup 13)C values (Table 1). More recently, Earth-based spectroscopic measurements of the martian atmosphere have measured the martian CO2 to be depleted in delta(sup 13)C relative to CO2 in the terrestrial atmosphere. The spectroscopic measurements performed by Krasnopolsky et al. were reported with approx.2% uncertainties which are much smaller than the Viking measurements, but still remain very large in comparison to the magnitude of carbon and oxygen isotope fractionations under martian surface conditions. The Thermal Evolved Gas Analyzer (TEGA) instrument on the Mars Phoenix Lander included a magnetic sector mass spectrometer (EGA) which had the goal of measuring the isotopic composition of martian atmospheric CO2 to within 0.5%. The mass spectrometer is a miniature magnetic sector instrument intended to measure both the martian atmosphere as well as gases evolved from heating martian soils. Ions produced in the ion source are drawn out by a high voltage and focused by a magnetic field onto a set of collector slits. Four specific trajectories are selected to cover the mass ranges, 0.7 - 4, 7 - 35, 14 - 70, and 28 - 140 Da. Using four channels reduces the magnitude of the mass scan and provides simultaneous coverage of the mass ranges. Channel electron multiplier (CEM) detectors that operate in the pulse counting mode detect the ion beams.

Niles, P.B.↗

Trajectory specification for high capacity air traffic control

Method and system for analyzing and processing information on one or more aircraft flight paths, using a four-dimensional coordinate system including three Cartesian or equivalent coordinates (x, y, z) and a fourth coordinate .delta. that corresponds to a distance estimated along a reference flight path to a nearest reference path location corresponding to a present location of the aircraft. Use of the coordinate .delta., rather than elapsed time t, avoids coupling of along-track error into aircraft altitude and reduces effects of errors on an aircraft landing site. Along-track, cross-track and/or altitude errors are estimated and compared with a permitted error bounding space surrounding the reference flight path.

Paielli, Russell A.↗

Requirements for Large Eddy Simulation Computations of Variable-Speed Power Turbine Flows

Variable-speed power turbines (VSPTs) operate at low Reynolds numbers and with a wide range of incidence angles. Transition, separation, and the relevant physics leading to them are important to VSPT flow. Higher fidelity tools such as large eddy simulation (LES) may be needed to resolve the flow features necessary for accurate predictive capability and design of such turbines. A survey conducted for this report explores the requirements for such computations. The survey is limited to the simulation of two-dimensional flow cases and endwalls are not included. It suggests that a grid resolution necessary for this type of simulation to accurately represent the physics may be of the order of Delta(x)+=45, Delta(x)+ =2 and Delta(z)+=17. Various subgrid-scale (SGS) models have been used and except for the Smagorinsky model, all seem to perform well and in some instances the simulations worked well without SGS modeling. A method of specifying the inlet conditions such as synthetic eddy modeling (SEM) is necessary to correctly represent the inlet conditions.

Grid↗

Optimal Fragmentation and Dispersion of Hazardous Near-Earth Objects

The complex problem of protecting the Earth from the possibility of a catastrophic impact by a hazardous near-Earth object (NEO) has been recently reassessed in [1]. In a letter on NEOs from the White House Office of Science and Technology Policy (OSTP) to the U.S. Senate and Congress in 2010, the White House OSTP strongly recommended that NASA take the lead in conducting research activities for NEO detection, characterization, and deflection technologies. Furthermore, President Obama's new National Space Policy specifically directs NASA to "pursue capabilities, in cooperation with other departments, agencies, and commercial partners, to detect, track, catalog, and characterize NEOs to reduce the risk of harm to humans from an unexpected impact on our planet." The Planetary Defense Task Force of the NASA Advisory Council also recommended that the NASA Office of the Chief Technologist (OCT) begin efforts to investigate asteroid deflection techniques. With national interest growing in the United States, the NEO threat detection and mitigation problem was recently identified as one of NASA's Space Technology Grand Challenges. An innovative solution to NASA's NEO Impact Threat Mitigation Grand Challenge problem was developed through a NIAC Phase I study (9/16/11 - 9/15/12), and it will be further investigated for a NIAC Phase II study (9/10/12 - 9/9/14). Various NEO deflection technologies, including nuclear explosions, kinetic impactors, and slow-pull gravity tractors, have been proposed and examined during the past two decades. Still, there is no consensus on how to reliably deflect or disrupt hazardous NEOs in a timely manner. It is expected that the most probable mission scenarios will have a mission lead time much shorter than 10 years, so the use of nuclear explosives becomes the most feasible method for planetary defense. Direct intercept missions with a short warning time will result in arrival closing velocities of 10-30 kilometers per second with respect to the target asteroid. Given such a large arrival delta V requirement, a rendezvous mission to the target asteroid is infeasible with existing launch vehicles. Furthermore, state-of-the-art penetrating subsurface nuclear explosion technology limits the penetrator's impact velocity to less than approximately 300 meters per second because higher impact velocities prematurely destroy the nuclear fuzing mechanisms. Therefore, significant advances in hypervelocity nuclear interceptor/ penetrator technology must be achieved to enable a last-minute nuclear disruption mission with intercept velocities as high as 30 kilometers per second. Consequently, a HAIV (Hypervelocity Asteroid Intercept Vehicle) mission architecture (Figure 1.1), which blends a hypervelocity kinetic impactor with a subsurface nuclear explosion for optimal fragmentation and dispersion of hazardous NEOs, has been developed through a Phase I study, and it will be further developed and validated through a Phase II study.

Asteroid↗

Flutter analysis and testing of pairs of aerodynamically interfering delta wings

To examine the effect on flutter of the aerodynamic interference between pairs of closely spaced delta wings, several structurally uncoupled 1/80th-scale models were studied by experiment and analysis. Flutter test boundaries run in a 26-in transonic blowdown wind tunnel were compared with subsonic analytical results generated using the doublet lattice method. Trends for several combinations of vertical and longitudinal wing separation showed that flutter speeds can be significantly lowered in closely spaced configurations. For some configurations, a new flutter mechanism, characterized by coupling of the flexible modes from both surfaces at a distinctive flutter frequency, was predicted and observed.

Chipman, R. R.↗

MAGSAT science investigations

A preliminary magnetization model for the U.S. was produced from MAGSAT data. The double grid processing method resulted in a 100 km source spacing. While some spurious features are present in the map, much of the detail appears to be real. The determination of the optional source spacing for inversion of MAGSAT delta B data is in progress. A fine attitude vector data set was obtained for comparison with the theoretical vector field associated with the magnetization model.

Source record↗

Creep crack-growth: A new path-independent integral (T sub c), and computational studies

The development of valid creep fracture criteria is considered. Two path-independent integral parameters which show some degree of promise are the C* and (Delta T)sub c integrals. The mathematical aspects of these parameters are reviewed by deriving generalized vector forms of the parameters using conservation laws which are valid for arbitrary, three dimensional, cracked bodies with crack surface tractions (or applied displacements), body forces, inertial effects, and large deformations. Two principal conclusions are that (Delta T)sub c has an energy rate interpretation whereas C* does not. The development and application of fracture criteria often involves the solution of boundary/initial value problems associated with deformation and stresses. The finite element method is used for this purpose. An efficient, small displacement, infinitesimal strain, displacement based finite element model is specialized to two dimensional plane stress and plane strain and to power law creep constitutive relations. A mesh shifting/remeshing procedure is used for simulating crack growth. The model is implemented with the quartz-point node technique and also with specially developed, conforming, crack-tip singularity elements which provide for the r to the n-(1+n) power strain singularity associated with the HRR crack-tip field. Comparisons are made with a variety of analytical solutions and alternate numerical solutions for a number of problems.

Stonesifer, R. B.↗

Delta function excitation of waves in the earth's ionosphere

Excitation of the earth's ionosphere by delta function current sheets is considered, and the temporal and spatial evolution of wave packets is analyzed for a two-component collisional F2 layer. Approximations of an inverse Fourier-Laplace transform via saddle point methods provide plots of typical wave packets. These illustrate cold plasma wave theory and may be used as a diagnostic tool since it is possible to relate specific features, e.g., the frequency of a modulation envelope, to plasma parameters such as the electron cyclotron frequency. It is also possible to deduce the propagation path length and orientation of a remote radio beacon.

Vidmar, R. J.↗

An evaluation of Global Positioning System data for Landsat-4 orbit determination

The Navstar Global Positioning System (GPS) navigation data obtained from an experimental GPS receiver/processor package (GPSPAC) onboard the Landsat-4 spacecraft are evaluated for their accuracy and reliability for use in Landsat-4 orbit determination. Different orbit determination scenarios and methods are considered. The parameters estimated include a subset of three Landsat-4 clock parameters and an atmospheric drag coefficient, in addition to the orbital elements. It is found that Landsat-4 orbit solutions based on delta pseudorange data generally agree with the definitive solutions to the 50-meter level.

Fang, B. T.↗