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 559 records · Page 31

Computation of space shuttle flow fields using noncentered finite-difference schemes.

Second- and third-order, noncentered finite-difference schemes are described for the numerical solution of the hyperbolic equations of fluid dynamics. The advantages of noncentered methods over the more conventional centered schemes are: simpler programming logic, nonhomogeneous terms are easily included, and generalization to multidimensional problems is direct. Second- and third-order methods are compared with regard to dissipative and dispersive errors and shock-capturing ability. These schemes are then used in a shock-capturing technique to determine the inviscid, supersonic flow field surrounding space shuttle vehicles (SSV). Resulting flow fields about typical pointed and blunted, delta-winged SSVs at angle of attack are presented and compared with experiment.

Kutler, P.↗

A digital audio/video interleaving system

A method of interleaving an audio signal with its associated video signal for simultaneous transmission or recording, and the subsequent separation of the two signals, is described. Comparisons are made between the new audio signal interleaving system and the Skylab Pam audio/video interleaving system, pointing out improvements gained by using the digital audio/video interleaving system. It was found that the digital technique is the simplest, most effective and most reliable method for interleaving audio and/or other types of data into the video signal for the Shuttle Orbiter application. Details of the design of a multiplexer capable of accommodating two basic data channels, each consisting of a single 31.5-kb/s digital bit stream are given. An adaptive slope delta modulation system is introduced to digitize audio signals, producing a high immunity of work intelligibility to channel errors, primarily due to the robust nature of the delta-modulation algorithm.

Richards, R. W.↗

A strawman targeting technique for spinner spacecraft in the presence of pointing uncertainty

A graphical method that allows the selection of the attitude orientation of a spinner spacecraft prior to the application of a single fixed magnitude velocity change to effect a desired orbit transfer is presented. The selection of the final pointing vector along or opposite to which the desired Delta-V is applied, takes into account pointing errors due to initial attitude determination uncertainty and execution errors introduced during the rhumb line precession of the spin axis. The characterization of the target space in terms of any achieved orbit parameters of interest as well as the final pointing uncertainty is sufficient to determine the right final attitude that insures satisfaction of mission requirements.

Kechichian, J. A.↗

A user's manual for the Electromagnetic Surface Patch code: ESP version 3

This report serves as a user's manual for Version III of the Electromagnetic Surface Patch Code or ESP code. ESP is user-oriented, based on the method of moments (MM) for treating geometries consisting of an interconnection of thin wires and perfectly conducting polygonal plates. Wire/plate junctions must be about 0.1 lambda or more from any plate edge. Several plates may intersect along a common edge. Excitation may be by either a delta-gap voltage generator or by a plane wave. The thin wires may have finite conductivity and also may contain lumped loads. The code computes most of the usual quantities of interest such as current distribution, input impedance, radiation efficiency, mutual coupling, far zone gain patterns (both polarizations) and radar-cross-section (both/cross polarizations).

Newman, E. H.↗

An overview of selected NASP aeroelastic studies at the NASA Langley Research Center

Following an initial discussion of the NASP flight environment, the results of recent aeroelastic testing of NASP-type highly swept delta-wing models in Langley's Transonic Dynamics Tunnel (TDT) are summarized. Subsonic and transonic flutter characteristics of a variety of these models are described, and several analytical codes used to predict flutter of these models are evaluated. These codes generally provide good, but conservative predictions of subsonic and transonic flutter. Also, test results are presented on a nonlinear transonic phenomena known as aileron buzz which occurred in the wind tunnel on highly swept delta wings with full-span ailerons. An analytical procedure which assesses the effects of hypersonic heating on aeroelastic instabilities (aerothermoelasticity) is also described. This procedure accurately predicted flutter of a heated aluminum wing on which experimental data exists. Results are presented on the application of this method to calculate the flutter characteristics of a fine-element model of a generic NASP configuration. Finally, it is demonstrated analytically that active controls can be employed to improve the aeroelastic stability and ride quality of a generic NASP vehicle flying at hypersonic speeds.

Spain, Charles V.↗

Divergence Free High Order Filter Methods for Multiscale Non-ideal MHD Flows

Low-dissipative high order filter finite difference methods for long time wave propagation of shock/turbulence/combustion compressible viscous MHD flows has been constructed. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (Delta . B) numerical error in the sense that no standard divergence cleaning is required. For certain 2-D MHD test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Yee, H. C.↗

Kepler Observations of Rapid Optical Variability in the BL Lac Object W2r192+42

We present the first Kepler monitoring of a strongly variable BL Lac, W2R1926+42. The light curve covers 181 days with approx. 0.2% errors, 30 minute sampling and >90% duty cycle, showing numerous delta-I/I > 25% flares over timescales as short as a day. The flux distribution is highly skewed and non-Gaussian. The variability shows a strong rms-flux correlation with the clearest evidence to date for non-linearity in this relation. We introduce a method to measure periodograms from the discrete autocorrelation function, an approach that may be well-suited to a wide range of Kepler data. The periodogram is not consistent with a simple power-law, but shows a flattening at frequencies below 7x10(exp -5) Hz. Simple models of the power spectrum, such as a broken power law, do not produce acceptable fits, indicating that the Kepler blazar light curve requires more sophisticated mathematical and physical descriptions than currently in use.

galaxies↗

Arcminute fluctuations in the microwave background from clusters of galaxies

A method for computing arcmin microwave fluctuations produced by Compton scattering of the cosmic background photons by hot electrons in clusters of galaxies is described. Microwave images of the sky for a range of Omega and primordial fluctuation spectral index n are generated which are then 'observed' to determine Delta T/T in precisely the same manner as actual observations to determine if the cluster-induced fluctuations are consistent with the measured upper limit. The geometry used by Uson and Wilkinson (1984) in the NRAO experiment and Readhead et al. (1989) in the OVRO experiment are applied to the simulated images. The 95 percent confidence lower limit for Omega is found to be about 1/10 for n = -1 (which approximates the CDM mass spectrum for clusters), while for n = 0 it is 1/7; for n = +1 the limit is 1/5 if the gas density profile extends to five core radii.

Markevitch, M.↗

A Parallel Incompressible Navier-Stokes Solver with Multigrid Iterations

We developed a parallel, numerically accurate and stable, and computationally efficient finate-difference incompressible Navier-Stokes (N-S) fluid flow solver. The solver runs on both sequential and massively parallel computers. The numerical method used here is a second-order projection method on a staggered grid. The code is highly modular and it can be used either as a stand-alone flow solver and or a template code which can be adapted or expanded to a specific application. Numerical results and parallel performances of our code on Intel Delta and Paragon are reported.

Navier-Stokes solver↗

Estimating Mangrove Aboveground Biomass from Airborne LiDAR Data: a Case Study from the Zambezi River Delta

Mangroves are ecologically and economically important forested wetlands with the highest carbon (C) density of all terrestrial ecosystems. Because of their exceptionally large C stocks and importance as a coastal buffer, their protection and restoration has been proposed as an effective mitigation strategy for climate change. The inclusion of mangroves in mitigation strategies requires the quantification of C stocks (both above and belowground) and changes to accurately calculate emissions and sequestration. A growing number of countries are becoming interested in using mitigation initiatives, such as REDD+ (reducing emissions from deforestation and forest degradation), in these unique coastal forests. However, it is not yet clear how methods to measure C traditionally used for other ecosystems can be modified to estimate biomass in mangroves with the precision and accuracy needed for these initiatives. Airborne Lidar (ALS) data has often been proposed as the most accurate way for larger scale assessments but the application of ALS for coastal wetlands is scarce, primarily due to a lack of contemporaneous ALS and field measurements. Here, we evaluated the variability in field and Lidar-based estimates of aboveground biomass (AGB) through the combination of different local and regional allometric models and standardized height metrics that are comparable across spatial resolutions and sensor types, the end result being a simplified approach for accurately estimating mangrove AGB at large scales and determining the uncertainty by combining multiple allometric models.We then quantified wall-to-wall AGB stocks of a tall mangrove forest in the Zambezi Delta, Mozambique. Our results indicate that the Lidar H100 height metric correlates well with AGB estimates, with R(exp 2) between 0.80 and 0.88 and RMSE of 33% or less. When comparing Lidar H100 AGB derived from three allometric models, mean AGB values range from 192 Mg ha(exp −1) up to 252 Mg ha(exp −1).We suggest the best model to predict AGB was based on the East Africa specific allometry and a power-based regression that used Lidar H100 as the height input with an R(exp 2) of 0.85 and an RMSE of 122 Mg ha(exp −1) or 33%. The total AGB of the Lidar inventoried mangrove area (6654 ha) was 1 350 902 Mg with a mean AGB of 203 Mg ha(exp −1) +/-166 Mg ha(exp −1). Because the allometry suggested here was developed using standardized height metrics, it is recommended that the models can generate AGB estimates using other remote sensing instruments that are more readily accessible over other mangrove ecosystems on a large scale, and as part of future carbon monitoring efforts in mangroves.

Fatoyinbo, Temilola↗

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗

Path planning for assembly of strut-based structures

A path planning method with collision avoidance for a general single chain nonredundant or redundant robot is proposed. Joint range boundary overruns are also avoided. The result is a sequence of joint vectors which are passed to a trajectory planner. A potential field algorithm in joint space computes incremental joint vectors delta-q = delta-q(sub a) + delta-q(sub c) + delta-q(sub r). Adding delta-q to the robot's current joint vector leads to the next step in the path. Delta-q(sub a) is obtained by computing the minimum norm solution of the underdetermined linear system J delta-q(sub a) = x(sub a) where x(sub a) is a translational and rotational force vector that attracts the robot to its goal position and orientation. J is the manipulator Jacobian. Delta-q(sub c) is a collision avoidance term encompassing collisions between the robot (links and payload) and obstacles in the environment as well as collisions among links and payload of the robot themselves. It is obtained in joint space directly. Delta-q(sub r) is a function of the current joint vector and avoids joint range overruns. A higher level discrete search over candidate safe positions is used to provide alternatives in case the potential field algorithm encounters a local minimum and thus fails to reach the goal. The best first search algorithm A* is used for graph search. Symmetry properties of the payload and equivalent rotations are exploited to further enlarge the number of alternatives passed to the potential field algorithm.

Muenger, Rolf↗

The Effect of Infrastructure Sharing in Estimating Operations Cost of Future Space Transportation Systems

NASA and the aerospace industry are extremely serious about reducing the cost and improving the performance of launch vehicles both manned or unmanned. In the aerospace industry, sharing infrastructure for manufacturing more than one type spacecraft is becoming a trend to achieve economy of scale. An example is the Boeing Decatur facility where both Delta II and Delta IV launch vehicles are made. The author is not sure how Boeing estimates the costs of each spacecraft made in the same facility. Regardless of how a contractor estimates the cost, NASA in its popular cost estimating tool, NASA Air force Cost Modeling (NAFCOM) has to have a method built in to account for the effect of infrastructure sharing. Since there is no provision in the most recent version of NAFCOM2002 to take care of this, it has been found by the Engineering Cost Community at MSFC that the tool overestimates the manufacturing cost by as much as 30%. Therefore, the objective of this study is to develop a methodology to assess the impact of infrastructure sharing so that better operations cost estimates may be made.

Sundaram, Meenakshi↗

Small Launch Vehicle Trade Space Definition: Development of a Zero Level Mass Estimation Tool with Trajectory Validation

Recent high level interest in the capability of small launch vehicles has placed significant demand on determining the trade space these vehicles occupy. This has led to the development of a zero level analysis tool that can quickly determine the minimum expected vehicle gross liftoff weight (GLOW) in terms of vehicle stage specific impulse (Isp) and propellant mass fraction (pmf) for any given payload value. Utilizing an extensive background in Earth to orbit trajectory experience a total necessary delta v the vehicle must achieve can be estimated including relevant loss terms. This foresight into expected losses allows for more specific assumptions relating to the initial estimates of thrust to weight values for each stage. This tool was further validated against a trajectory model, in this case the Program to Optimize Simulated Trajectories (POST), to determine if the initial sizing delta v was adequate to meet payload expectations. Presented here is a description of how the tool is setup and the approach the analyst must take when using the tool. Also, expected outputs which are dependent on the type of small launch vehicle being sized will be displayed. The method of validation will be discussed as well as where the sizing tool fits into the vehicle design process.

Waters, Eric D.↗

Four-body trajectory optimization

A collection of typical three-body trajectories from the L1 libration point on the sun-earth line to the earth is presented. These trajectories in the sun-earth system are grouped into four distinct families which differ in transfer time and delta V requirements. Curves showing the variations of delta V with respect to transfer time, and typical two and three-impulse primer vector histories, are included. The development of a four-body trajectory optimization program to compute fuel optimal trajectories between the earth and a point in the sun-earth-moon system are also discussed. Methods for generating fuel optimal two-impulse trajectories which originate at the earth or a point in space, and fuel optimal three-impulse trajectories between two points in space, are presented. A brief qualitative comparison of these methods is given. An example of a four-body two-impulse transfer from the Li libration point to the earth is included.

Pu, C. L.↗

Utility of correlation techniques in gravity and magnetic interpretation

Internal correspondence uses Poisson's Theorem in a moving-window linear regression analysis between the anomalous first vertical derivative of gravity and total magnetic field reduced to the pole. The regression parameters provide critical information on source characteristics. The correlation coefficient indicates the strength of the relation between magnetics and gravity. Slope value gives delta j/delta sigma estimates of the anomalous source. The intercept furnishes information on anomaly interference. Cluster analysis consists of the classification of subsets of data into groups of similarity based on correlation of selected characteristics of the anomalies. Model studies are used to illustrate implementation and interpretation procedures of these methods, particularly internal correspondence. Analysis of the results of applying these methods to data from the midcontinent and a transcontinental profile shows they can be useful in identifying crustal provinces, providing information on horizontal and vertical variations of physical properties over province size zones, validating long wavelength anomalies, and isolating geomagnetic field removal problems.

Chandler, V. W.↗

Optimum two-impulse transfers for preliminary interplanetary trajectory design

A method for numerically determining the optimum two-impulse transfer between two positions in two different heliocentric orbits is described. The formulas are based on earlier work by McCue and Bender (1965) for optimum two-impulse transfers between two orbits about the earth. The contours of minimum total Delta-V(Delta-V to launch plus Delta-V to achieve a given final orbit at the target body) are plotted on axes of true anomaly of the launch body at launch and true anomaly of the target body at arrival. The optimum transfers obtained are 'time-open' and the Delta-V from the initial orbit to the final orbit is found using patched conics. These 'prime rib' plots are useful for selecting initial conditions in the early phases of mission design.

Hulkower, N. D.↗