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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Cart3D Simulations for the First AIAA Sonic Boom Prediction Workshop

Simulation results for the First AIAA Sonic Boom Prediction Workshop (LBW1) are presented using an inviscid, embedded-boundary Cartesian mesh method. The method employs adjoint-based error estimation and adaptive meshing to automatically determine resolution requirements of the computational domain. Results are presented for both mandatory and optional test cases. These include an axisymmetric body of revolution, a 69deg delta wing model and a complete model of the Lockheed N+2 supersonic tri-jet with V-tail and flow through nacelles. In addition to formal mesh refinement studies and examination of the adjoint-based error estimates, mesh convergence is assessed by presenting simulation results for meshes at several resolutions which are comparable in size to the unstructured grids distributed by the workshop organizers. Data provided includes both the pressure signals required by the workshop and information on code performance in both memory and processing time. Various enhanced techniques offering improved simulation efficiency will be demonstrated and discussed.

Predictions↗
NASA NTRS2012

A Newton-Raphson Method Approach to Adjusting Multi-Source Solar Simulators

NASA Glenn Research Center has been using an in house designed X25 based multi-source solar simulator since 2003. The simulator is set up for triple junction solar cells prior to measurements b y adjusting the three sources to produce the correct short circuit current, lsc, in each of three AM0 calibrated sub-cells. The past practice has been to adjust one source on one sub-cell at a time, iterating until all the sub-cells have the calibrated Isc. The new approach is to create a matrix of measured lsc for small source changes on each sub-cell. A matrix, A, is produced. This is normalized to unit changes in the sources so that Ax(delta)s = (delta)isc. This matrix can now be inverted and used with the known Isc differences from the AM0 calibrated values to indicate changes in the source settings, (delta)s = A ·'x.(delta)isc This approach is still an iterative one, but all sources are changed during each iteration step. It typically takes four to six steps to converge on the calibrated lsc values. Even though the source lamps may degrade over time, the initial matrix evaluation i s not performed each time, since measurement matrix needs to be only approximate. Because an iterative approach is used the method will still continue to be valid. This method may become more important as state-of-the-art solar cell junction responses overlap the sources of the simulator. Also, as the number of cell junctions and sources increase, this method should remain applicable.

Snyder, David B.↗
NASA NTRS1975

Abort-once-around entry dispersion corridor analysis

Abort-Once-Around (AOA) entry dispersion corridors were determined for Shuttle Mission 3A. These corridors are presented as plots of entry interface flight path angle versus range to target. The methods used to determine the corridors are discussed. Major dispersion sources are discussed and results presented. While specific trajectory inputs and constraints are subject to change, the dispersion corridors show the trends under consideration. The corridors presented show the delta V advantage of the two-burn over the one-burn AOA. The atmospheric dispersion study illustrates the need to target for the seasonal atmosphere. The 40 deg/30 deg angle of attack (alpha) (TPS design) profile did not provide adequate crossrange capability with worst case aerodynamic dispersions. This problem was alleviated with the change to a 38 deg/28 deg alpha profile. The backface temperatures calculated were generally higher than the present limits.

Kyle, H. C.↗
NASA NTRS1985

Some design characteristics of the AMPTE turn and orbit change maneuvers

The maneuvers carried out by the Active Magnetospheric Particle Tracer Explorers (AMPTE) including the Charge Composition Explorer (CCE) and the Ion Release Module (IRM) spacecraft are analyzed. Analytic and graphical methods are developed in order to carry out sensitivity analyses that helped design the nominal maneuvers, by taking into account errors in burn initiating time, motor performance, and spin axis pointing. A tradeoff analysis between errors in timing and Delta V magnitude is shown for the IRM orbit transfer, and a technique that allows for the determination of the attitude of spinner spacecraft by way of the observed Doppler shift resulting from an unbalanced turn is investigated.

Kechichian, J. A.↗
NASA NTRS1987

A fast radiation parameterization for atmospheric circulation models

Several currently available broadband parameterizations for longwave and shortwave radiation have been combined to produce a computationally fast radiation parameterization that is well suited for atmospheric circulation models. The main features of the parameterization are the ability to include overlapping partly cloudy layers in the longwave, the use of a delta-Eddington technique to treat clouds in the shortwave, and a computational structure that is amenable to vectorization on supercomputers. Selected results of off-line one-dimensional computations using the code have been compared with more rigorous methods as part of an international intercomparison program and found to be quite accurate.

HARSHVARDHAN↗
NASA NTRS1990

Vapor pressures of acetylene at low temperatures

The atmospheres of many of the outer planets and their satellites contain a large number of hydrocarbon species. In particular, acetylene (C2H2) has been identified at Jupiter, Saturn and its satellite Titan, Uranus and Neptune. In the lower atmospheres of these planets, where colder temperatures prevail, the condensation and/or freezing of acetylene is probable. In order to obtain accurate models of the acetylene in these atmospheres, it is necessary to have a complete understanding of its vapor pressures at low temperatures. Vapor pressures at low temperatures for acetylene are being determined. The vapor pressures are measured with two different techniques in order to cover a wide range of temperatures and pressures. In the first, the acetylene is placed in a sample tube which is immersed in a low temperature solvent/liquid nitrogen slush bath whose temperature is measured with a thermocouple. The vapor pressure is then measured directly with a capacitance manometer. For lower pressures, a second technique which was called the thin-film infrared method (TFIR) was developed. It involves measuring the disappearance rate of a thin film of acetylene at a particular temperature. The spectra are then analyzed using previously determined extinction coefficient values, to determine the disappearance rate R (where R = delta n/delta t, the number of molecules that disappear per unit time). This can be related to the vapor pressure directly. This technique facilitates measurement of the lower temperatures and pressures. Both techniques have been calibrated using CO2, and have shown good agreement with the existing literature data.

Masterson, C. M.↗
NASA NTRS1994

A rapid radiative transfer model for reflection of solar radiation

A rapid analytical radiative transfer model for reflection of solar radiation in plane-parallel atmospheres is developed based on the Sobolev approach and the delta function transformation technique. A distinct advantage of this model over alternative two-stream solutions is that in addition to yielding the irradiance components, which turn out to be mathematically equivalent to the delta-Eddington approximation, the radiance field can also be expanded in a mathematically consistent fashion. Tests with the model against a more precise multistream discrete ordinate model over a wide range of input parameters demonstrate that the new approximate method typically produces average radiance differences of less than 5%, with worst average differences of approximately 10%-15%. By the same token, the computational speed of the new model is some tens to thousands times faster than that of the more precise model when its stream resolution is set to generate precise calculations.

Xiang, X.↗
NASA NTRS2014

Design of the VISITOR Tool: A Versatile ImpulSive Interplanetary Trajectory OptimizeR

The design of trajectories for interplanetary missions represents one of the most complex and important problems to solve during conceptual space mission design. To facilitate conceptual mission sizing activities, it is essential to obtain sufficiently accurate trajectories in a fast and repeatable manner. To this end, the VISITOR tool was developed. This tool modularly augments a patched conic MGA-1DSM model with a mass model, launch window analysis, and the ability to simulate more realistic arrival and departure operations. This was implemented in MATLAB, exploiting the built-in optimization tools and vector analysis routines. The chosen optimization strategy uses a grid search and pattern search, an iterative variable grid method. A genetic algorithm can be selectively used to improve search space pruning, at the cost of losing the repeatability of the results and increased computation time. The tool was validated against seven flown missions: the average total mission (Delta)V offset from the nominal trajectory was 9.1%, which was reduced to 7.3% when using the genetic algorithm at the cost of an increase in computation time by a factor 5.7. It was found that VISITOR was well-suited for the conceptual design of interplanetary trajectories, while also facilitating future improvements due to its modular structure.

Corpaccioli, Luca↗
NASA NTRS1956

A method for calculation of free-space sound pressures near a propeller in flight including considerations of the chordwise blade loading

This report presents tabulated values of certain definite integral that are involved in the calculation of near-field propeller noise when the chordwise forces are assumed to be either uniform or of a Dirac delta type. The tabulations are over a wide range of operating conditions and are useful for estimating propeller noise when either the concept of an effective radius or radial distributions of forces are considered. Use of the tabulations is illustrated by several examples of calculated results for some specific propellers.

NOISE↗
NASA NTRS1972

Surface-flow, pressure, and heat-transfer studies on two conical delta wings at a Mach number of 6

An experimental investigation of the surface flow, pressures, and heat transfer on two conical delta wings having attached leading-edge shocks has been conducted at a Mach number of 6. The angle of attack was varied between 0 deg and 12 deg. The pressure data were compared with predictions obtained by the method-of-lines technique, and the heating data were compared with the heating levels predicted by the Spalding-Chi method.

Hefner, J. N.↗
NASA NTRS1984

Comparison of numerical techniques for integration of stiff ordinary differential equations arising in combustion chemistry

The efficiency and accuracy of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations are compared. The methods examined include two general-purpose codes, EPISODE and LSODE, and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an interactive solution of the algebraic energy conservation equation to compute the temperature does not result in significant errors. In addition, this method is more efficient than evaluating the temperature by integrating its time derivative. Significant reductions in computational work are realized by updating the rate constants (k = at(supra N) N exp(-E/RT) only when the temperature change exceeds an amount delta T that is problem dependent. An approximate expression for the automatic evaluation of delta T is derived and is shown to result in increased efficiency.

Radhakrishnan, K.↗
NASA NTRS1994

Sonic boom predictions using a solution-adaptive full-potential code

The TranAir full-potential code, which utilizes a solution-adaptive, unstructured Cartesian grid, was used to compute near-field pressure signatures of three generic geometries: (1) a cone-cylinder; (2) a low aspect ratio rectangular wing; and (3) a delta-wing/body. Pressure signatures computed by TranAir are extrapolated to near- and midfield distances using a quasilinear extrapolation technique. Results are compared with wind-tunnel data and, where possible, with extrapolated wind-tunnel data and Whitham-method results. In general, extrapolated TranAir data is in good agreement with experimental and Whitham-method pressure signatures, but peak pressures are slightly underpredicted. It is felt that data computed too close to the model caused some of the discrepancies. The flowfield grid structure was also found to be an important issue in obtaining accurate off-body data. Results were encouraging enough to proceed with the application of TranAir to conceptual configurations in support of the NASA high speed research program (HSRP).

Madson, Michael D.↗
NASA NTRS2003

Migration of Dust Particles from Comet 2P Encke

We investigated the migration of dust particles under the gravitational influence of all planets (except for Pluto), radiation pressure, Poynting-Robertson drag and solar wind drag for Beta equal to 0.002, 0.004, 0.01, 0.05, 0.1, 0.2, and 0.4. For silicate particles such values of Beta correspond to diameters equal to about 200, 100, 40, 9, 4, 2, and 1 microns, respectively. We used the Bulirsh-Stoer method of integration, and the relative error per integration step was taken to be less than lo-'. Initial orbits of the particles were close to the orbit of Comet 2P Encke. We considered initial particles near perihelion (runs denoted as Delta tsub o, = 0), near aphelion (Delta tsub o, = 0.5), and also studied their initial positions when the comet moved for Pa/4 after perihelion passage (such runs are denoted as Delta tsub o, =i 0.25), where Pa is the period of the comet. Variations in time T when perihelion was passed was varied with a step 0.1 day for series 'S' and with a step 1 day for series 'L'. For each Beta we considered N = 101 particles for "S" runs and 150 particles for "L" runs.

Ipatov, S. I.↗
NASA NTRS2006

Method of Estimating Continuous Cooling Transformation Curves of Glasses

A method is proposed for estimating the critical cooling rate and continuous cooling transformation (CCT) curve from isothermal TTT data of glasses. The critical cooling rates and CCT curves for a group of lithium disilicate glasses containing different amounts of Pt as nucleating agent estimated through this method are compared with the experimentally measured values. By analysis of the experimental and calculated data of the lithium disilicate glasses, a simple relationship between the crystallized amount in the glasses during continuous cooling, X, and the temperature of undercooling, (Delta)T, was found to be X = AR(sup-4)exp(B (Delta)T), where (Delta)T is the temperature difference between the theoretical melting point of the glass composition and the temperature in discussion, R is the cooling rate, and A and B are constants. The relation between the amount of crystallisation during continuous cooling and isothermal hold can be expressed as (X(sub cT)/X(sub iT) = (4/B)(sup 4) (Delta)T(sup -4), where X(sub cT) stands for the crystallised amount in a glass during continuous cooling for a time t when the temperature comes to T, and X(sub iT) is the crystallised amount during isothermal hold at temperature T for a time t.

Zhu, Dongmei↗
NASA NTRS2010

Guidance and Control System for a Satellite Constellation

A distributed guidance and control algorithm was developed for a constellation of satellites. The system repositions satellites as required, regulates satellites to desired orbits, and prevents collisions. 1. Optimal methods are used to compute nominal transfers from orbit to orbit. 2. Satellites are regulated to maintain the desired orbits once the transfers are complete. 3. A simulator is used to predict potential collisions or near-misses. 4. Each satellite computes perturbations to its controls so as to increase any unacceptable distances of nearest approach to other objects. a. The avoidance problem is recast in a distributed and locally-linear form to arrive at a tractable solution. b. Plant matrix values are approximated via simulation at each time step. c. The Linear Quadratic Gaussian (LQG) method is used to compute perturbations to the controls that will result in increased miss distances. 5. Once all danger is passed, the satellites return to their original orbits, all the while avoiding each other as above. 6. The delta-Vs are reasonable. The controller begins maneuvers as soon as practical to minimize delta-V. 7. Despite the inclusion of trajectory simulations within the control loop, the algorithm is sufficiently fast for available satellite computer hardware. 8. The required measurement accuracies are within the capabilities of modern inertial measurement devices and modern positioning devices.

Bryson, Jonathan Lamar↗
NASA NTRS1973

Unsteady transonic flow analysis for low aspect ratio, pointed wings.

Oswatitsch and Keune's parabolic method for steady transonic flow is applied and extended to thin slender wings oscillating in the sonic flow field. The parabolic constant for the wing was determined from the equivalent body of revolution. Laplace transform methods were used to derive the asymptotic equations for pressure coefficient, and the Adams-Sears iterative procedure was employed to solve the equations. A computer program was developed to find the pressure distributions, generalized force coefficients, and stability derivatives for delta, convex, and concave wing planforms.

Kimble, K. R.↗
NASA NTRS1986

Planform effects for low-fineness ratio multibody configurations at supersonic speeds

An experimental and theoretical investigation of planform effects on a low-fineness ratio multibody configuration has been conducted in NASA Langley Research Center's Unitary Plan Wind Tunnel at Mach numbers 1.6, 1.8, 2.0, and 2.16. Experimental and theoretical values of lift, drag, and pitching moment were obtained on three configurations which varied in outboard panel planform only. The three variations were at 65 deg delta, a 70/66 deg cranked arrow, and a 20 deg trapezoidal planform . The purpose of the study was to determine the effect of wing planform on the supersonic aerodynamics and to evaluate the ability of two existing linearized-theory aerodynamic methods to predict these effects. Experimental data showed that the planforms produced the lift, drag-due-to-lift, and pitching-moment characteristics typically found on single-body configurations. However, the data also showed that planform has a minimal influence on zero-lift drag, which is not typical of single-body configurations. Theoretical aerodynamic analysis indicated that codes based on linearized theory adequately predicted the effect of planform on the supersonic aerodynamics.

Mcmillin, S. N.↗
NASA NTRS1988

A new convergent point and distance modulus for the Hyades from radial velocities

The dynamics of the Hyades open star cluster is investigated analytically on the basis of the photoelectric radial velocities reported by Griffin et al. (1988). The methods used to determine the radial-velocity gradient, convergent point, and distance modulus are described, and the results are presented in tables and graphs and characterized in detail. Findings reported include convergent-point alpha = 98.2 + or - 1.1 deg and delta = 6.1 + or - 1.0 deg, space motion 48.0 + or - 0.27 km/s, distance to cluster center 45.4 + or - 2.1 pc, distance modulus 3.28 + or - 0.10, and cluster velocity dispersion 230 m/s.

Gunn, James E.↗