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

Applying Dynamical Systems Theory to Optimize Libration Point Orbit Stationkeeping Maneuvers for WIND

NASA's WIND mission has been operating in a large amplitude Lissajous orbit in the vicinity of the interior libration point of the Sun-Earth/Moon system since 2004. Regular stationkeeping maneuvers are required to maintain the orbit due to the instability around the collinear libration points. Historically these stationkeeping maneuvers have been performed by applying an incremental change in velocity, or (delta)v along the spacecraft-Sun vector as projected into the ecliptic plane. Previous studies have shown that the magnitude of libration point stationkeeping maneuvers can be minimized by applying the (delta)v in the direction of the local stable manifold found using dynamical systems theory. This paper presents the analysis of this new maneuver strategy which shows that the magnitude of stationkeeping maneuvers can be decreased by 5 to 25 percent, depending on the location in the orbit where the maneuver is performed. The implementation of the optimized maneuver method into operations is discussed and results are presented for the first two optimized stationkeeping maneuvers executed by WIND.

WIND↗
NASA NTRS2013

Tuning and Robustness Analysis for the Orion Absolute Navigation System

The Orion Multi-Purpose Crew Vehicle (MPCV) is currently under development as NASA's next-generation spacecraft for exploration missions beyond Low Earth Orbit. The MPCV is set to perform an orbital test flight, termed Exploration Flight Test 1 (EFT-1), some time in late 2014. The navigation system for the Orion spacecraft is being designed in a Multi-Organizational Design Environment (MODE) team including contractor and NASA personnel. The system uses an Extended Kalman Filter to process measurements and determine the state. The design of the navigation system has undergone several iterations and modifications since its inception, and continues as a work-in-progress. This paper seeks to show the efforts made to-date in tuning the filter for the EFT-1 mission and instilling appropriate robustness into the system to meet the requirements of manned space ight. Filter performance is affected by many factors: data rates, sensor measurement errors, tuning, and others. This paper focuses mainly on the error characterization and tuning portion. Traditional efforts at tuning a navigation filter have centered around the observation/measurement noise and Gaussian process noise of the Extended Kalman Filter. While the Orion MODE team must certainly address those factors, the team is also looking at residual edit thresholds and measurement underweighting as tuning tools. Tuning analysis is presented with open loop Monte-Carlo simulation results showing statistical errors bounded by the 3-sigma filter uncertainty covariance. The Orion filter design uses 24 Exponentially Correlated Random Variable (ECRV) parameters to estimate the accel/gyro misalignment and nonorthogonality. By design, the time constant and noise terms of these ECRV parameters were set to manufacturer specifications and not used as tuning parameters. They are included in the filter as a more analytically correct method of modeling uncertainties than ad-hoc tuning of the process noise. Tuning is explored for the powered-flight ascent phase, where measurements are scarce and unmodelled vehicle accelerations dominate. On orbit, there are important trade-off cases between process and measurement noise. On entry, there are considerations about trading performance accuracy for robustness. Process Noise is divided into powered flight and coasting ight and can be adjusted for each phase and mode of the Orion EFT-1 mission. Measurement noise is used for the integrated velocity measurements during pad alignment. It is also used for Global Positioning System (GPS) pseudorange and delta- range measurements during the rest of the flight. The robustness effort has been focused on maintaining filter convergence and performance in the presence of unmodeled error sources. These include unmodeled forces on the vehicle and uncorrected errors on the sensor measurements. Orion uses a single-frequency, non-keyed GPS receiver, so the effects due to signal distortion in Earth's ionosphere and troposphere are present in the raw measurements. Results are presented showing the efforts to compensate for these errors as well as characterize the residual effect for measurement noise tuning. Another robustness tool in use is tuning the residual edit thresholds. The trade-off between noise tuning and edit thresholds is explored in the context of robustness to errors in dynamics models and sensor measurements. Measurement underweighting is also presented as a method of additional robustness when processing highly accurate measurements in the presence of large filter uncertainties.

Holt, Greg N.↗
NASA NTRS1989

Determination of longitudinal aerodynamic derivatives using flight data from an icing research aircraft

A flight test was performed with the NASA Lewis Research Center's DH-6 icing research aircraft. The purpose was to employ a flight test procedure and data analysis method, to determine the accuracy with which the effects of ice on aircraft stability and control could be measured. For simplicity, flight testing was restricted to the short period longitudinal mode. Two flights were flown in a clean (baseline) configuration, and two flights were flown with simulated horizontal tail ice. Forty-five repeat doublet maneuvers were performed in each of four test configurations, at a given trim speed, to determine the ensemble variation of the estimated stability and control derivatives. Additional maneuvers were also performed in each configuration, to determine the variation in the longitudinal derivative estimates over a wide range of trim speeds. Stability and control derivatives were estimated by a Modified Stepwise Regression (MSR) technique. A measure of the confidence in the derivative estimates was obtained by comparing the standard error for the ensemble of repeat maneuvers, to the average of the estimated standard errors predicted by the MSR program. A multiplicative relationship was determined between the ensemble standard error, and the averaged program standard errors. In addition, a 95 percent confidence interval analysis was performed for the elevator effectiveness estimates, C sub m sub delta e. This analysis identified the speed range where changes in C sub m sub delta e could be attributed to icing effects. The magnitude of icing effects on the derivative estimates were strongly dependent on flight speed and aircraft wing flap configuration. With wing flaps up, the estimated derivatives were degraded most at lower speeds corresponding to that configuration. With wing flaps extended to 10 degrees, the estimated derivatives were degraded most at the higher corresponding speeds. The effects of icing on the changes in longitudinal stability and control derivatives were adequately determined by the flight test procedure and the MSR analysis method discussed herein.

Ranaudo, R. J.↗
NASA NTRS1978

A computer program for calculating aerodynamic characteristics of low aspect-ratio wings with partial leading-edge separation

The necessary information for using a computer program to predict distributed and total aerodynamic characteristics for low aspect ratio wings with partial leading-edge separation is presented. The flow is assumed to be steady and inviscid. The wing boundary condition is formulated by the Quasi-Vortex-Lattice method. The leading edge separated vortices are represented by discrete free vortex elements which are aligned with the local velocity vector at midpoints to satisfy the force free condition. The wake behind the trailing edge is also force free. The flow tangency boundary condition is satisfied on the wing, including the leading and trailing edges. The program is restricted to delta wings with zero thickness and no camber. It is written in FORTRAN language and runs on CDC 6600 computer.

Mehrotra, S. C.↗
NASA NTRS1981

Computation of simple three-dimensional wing-vortex interaction in transonic flow

An implicit delta form finite-difference algorithm for Euler equations in conservation law form has been used in preliminary calculations of three-dimensional wing-vortex interactions. Both steady and unsteady transonic flow wing-vortex interactions are computed. The computations themselves are meant to guide upcoming wind tunnel experiments of the same flow field. Various modifications to the numerical method that are intended to improve computational efficiency are also described and tested in both two- and three-dimensions.

Srinivasan, G. R.↗
NASA NTRS1994

Method for ambiguity resolution in range-Doppler measurements

A method for resolving range and Doppler target ambiguities when the target has substantial range or has a high relative velocity in which a first signal is generated and a second signal is also generated which is coherent with the first signal but at a slightly different frequency such that there exists a difference in frequency between these two signals of Delta f(sub t). The first and second signals are converted into a dual-frequency pulsed signal, amplified, and the dual-frequency pulsed signal is transmitted towards a target. A reflected dual-frequency signal is received from the target, amplified, and changed to an intermediate dual-frequency signal. The intermediate dual-frequency signal is amplified, with extracting of a shifted difference frequency Delta f(sub r) from the amplified intermediate dual-frequency signal done by a nonlinear detector. The final step is generating two quadrature signals from the difference frequency Delta f(sub t) and the shifted difference frequency Delta f(sub r) and processing the two quadrature signals to determine range and Doppler information of the target.

Heymsfield, Gerald M.↗
NASA NTRSDate not supplied

Exploring Unstructured Mesh Adaptation for Hybrid Reynolds-Averaged Navier–Stokes/Large Eddy Simulation

Mesh adaptation methods for the Reynolds-averaged Navier–Stokes (RANS) equations are rapidly maturing and beginning to impact the design of aerospace vehicles. RANS turbulence modeling improvements have slowed and may stagnate. Wall-modeled large eddy simulation (LES) and hybrid RANS/LES (HRLES) may provide an improved modeling capability but require specialized expertise to construct appropriate meshes and are considered be too ex-pensive for routine practical use. The realization of the CFD Vision 2030 Study includes improving geometry linkage, mesh generation/adaptation, and turbulence modeling/resolving methods for automated management of errors and uncertainties of physics-based, predictive modeling that can set the stage for ensuring a vehicle is in compliance with a regulation or specification (i.e., certification or qualification by analysis). An exploration of mesh adaptation for HRLES is performed to document synergies and challenges between mesh adaptation and HRLES. Vortex breakdown over a delta wing is examined to show the improvement of HRLES over RANS turbulence modeling approaches. A high lift configuration is shown to demonstrate complex geometry capability. Research and development opportunities are identified to advocate for continuing investments that may allow HRLES to enter routine practical use as a tool for aerospace vehicle analysis and design.

Michael A Park↗
NASA NTRS1995

A magnetospheric magnetic field model with flexible current systems driven by independent physical parameters

A tilt-dependent magnetic field model of the Earth's magnetosphere with variable magnetopause standoff distance is presented. Flexible analytic representations for the ring and cross-tail currents, each composed of the elements derived from the Tsyganenko and Usmanov (1982) model, are combined with the fully shielded vacuum dipole configurations of Voigt (1981). Although the current sheet does not warp in the y-z plane, changes in the shape and position of the neutral sheet with dipole tilt are consistent with both MHD equilibrium theory and observations. In addition, there is good agreement with observed Delta B profiles and the average equatorial contours of magnetic field magnitude. While the dipole field is rigorously shielded within the defined magnetopause, the ring and cross-tails currents are not similarly confined, consequently, the model's region of validity is limited to the inner magnetosphere. The model depends on four independent external parameters. We present a simple but limited method of simulating several substorm related magnetic field changes associated with the disrupion of the near-Earth cross-tail current sheet and collapse of the midnight magnetotail field region. This feature further facilitates the generation of magnetic field configuration time sequences useful in plasma convection simulations of real magnetospheric events.

Hilmer, Robert V.↗
NASA NTRS1998

Design and Calibration of the X-33 Flush Airdata Sensing (FADS) System

This paper presents the design of the X-33 Flush Airdata Sensing (FADS) system. The X-33 FADS uses a matrix of pressure orifices on the vehicle nose to estimate airdata parameters. The system is designed with dual-redundant measurement hardware, which produces two independent measurement paths. Airdata parameters that correspond to the measurement path with the minimum fit error are selected as the output values. This method enables a single sensor failure to occur with minimal degrading of the system performance. The paper shows the X-33 FADS architecture, derives the estimating algorithms, and demonstrates a mathematical analysis of the FADS system stability. Preliminary aerodynamic calibrations are also presented here. The calibration parameters, the position error coefficient (epsilon), and flow correction terms for the angle of attack (delta alpha), and angle of sideslip (delta beta) are derived from wind tunnel data. Statistical accuracy of' the calibration is evaluated by comparing the wind tunnel reference conditions to the airdata parameters estimated. This comparison is accomplished by applying the calibrated FADS algorithm to the sensed wind tunnel pressures. When the resulting accuracy estimates are compared to accuracy requirements for the X-33 airdata, the FADS system meets these requirements.

Whitmore, Stephen A.↗
NASA NTRS1984

The numerical solution of the Helmholtz equation for wave propagation problems in underwater acoustics

The Helmholtz Equation (-delta-K(2)n(2))u=0 with a variable index of refraction, n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. A numerical algorithm was developed and a computer code implemented that can effectively solve this equation in the intermediate frequency range. The equation is discretized using the finite element method, thus allowing for the modeling of complicated geometrices (including interfaces) and complicated boundary conditions. A global radiation boundary condition is imposed at the far field boundary that is exact for an arbitrary number of propagating modes. The resulting large, non-selfadjoint system of linear equations with indefinite symmetric part is solved using the preconditioned conjugate gradient method applied to the normal equations. A new preconditioner is developed based on the multigrid method. This preconditioner is vectorizable and is extremely effective over a wide range of frequencies provided the number of grid levels is reduced for large frequencies. A heuristic argument is given that indicates the superior convergence properties of this preconditioner.

Bayliss, A.↗
NASA NTRS1985

Numerical solutions of Navier-Stokes equations for a Butler wing

The flow field is simulated on the surface of a given delta wing (Butler wing) at zero incident in a uniform stream. The simulation is done by integrating a set of flow field equations. This set of equations governs the unsteady, viscous, compressible, heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. These equations are solved by the finite difference method, and results obtained for high-speed freestream conditions are compared with theoretical and experimental results. In this study, the Navier-Stokes equations are solved numerically. These equations are unsteady, compressible, viscous, and three-dimensional without neglecting any terms. The time dependency of the governing equations allows the solution to progress naturally for an arbitrary initial initial guess to an asymptotic steady state, if one exists. The equations are transformed from physical coordinates to the computational coordinates, allowing the solution of the governing equations in a rectangular parallel-piped domain. The equations are solved by the MacCormack time-split technique which is vectorized and programmed to run on the CDC VPS 32 computer.

Abolhassani, J. S.↗
NASA NTRS1986

Extensions of the concept of suction analogy to prediction of vortex lift effect

Flow field data for a double delta wing at low speed were used to determine the location of a vortex action point. The result was found to be consistent with what was determined for a delta wing. In supersonic flow, the action point location was determined empirically. For a wing with rounded leading edges, an assumption for initial vortex separation was shown to be equivalent to initial leading edge bubble separation for airfoils. A theoretical formulation by the section analogy to determine the delayed vortex separation on a cambered wing with rounded leading edges was presented. The method of suction analogy was further shown to be applicable to predicting the body vortex lift.

Lan, C. E.↗
NASA NTRS1994

Light stable isotope analysis of meteorites by ion microprobe

The main goal was to develop the necessary secondary ion mass spectrometer (SIMS) techniques to use a Cameca ims-4f ion microprobe to measure light stable isotope ratios (H, C, O and S) in situ and in non-conducting mineral phases. The intended application of these techniques was the analysis of meteorite samples, although the techniques that have been developed are equally applicable to the investigation of terrestrial samples. The first year established techniques for the analysis of O isotope ratios (delta O-18 and delta O-17) in conducting mineral phases and the measurement of S isotope ratios (delta S-34) in a variety of sulphide phases. In addition, a technique was developed to measure delta S-34 values in sulphates, which are insulators. Other research undertaken in the first year resulted in SIMS techniques for the measurement of wide variety of trace elements in carbonate minerals, with the aim of understanding the nature of alteration fluids in carbonaceous chondrites. In the second year we developed techniques for analyzing O isotope ratios in nonconducting mineral phases. These methods are potentially applicable to the measurement of other light stable isotopes such as H, C and S in insulators. Also, we have further explored the analytical techniques used for the analysis of S isotopes in sulphides by analyzing troilite in a number of L and H ordinary chondrites. This was done to see if there was any systematic differences with petrological type.

Mcsween, Harry Y., Jr.↗
NASA NTRS1995

Reduction, analysis, and properties of electric current systems in solar active regions

The specific attraction and, in large part, the significance of solar magnetograms lie in the fact that they give the most important data on the electric currents and the nonpotentiality of active regions. Using the vector magnetograms from the Marshall Space Flight Center (MSFC), we employ a unique technique in the area of data analysis for resolving the 180 deg ambiguity in order to calculate the spatial structure of the vertical electric current density. The 180 deg ambiguity is resolved by applying concepts from the nonlinear multivariable optimization theory. The technique is shown to be of particular importance in very nonpotential active regions. The characterization of the vertical electric current density for a set of vector magnetograms using this method then gives the spatial scale, locations, and magnitude of these current systems. The method, which employs an intermediate parametric function which covers the magnetogram and which defines the local `preferred' direction, minimizes a specific functional of the observed transverse magnetic field. The specific functional that is successful is the integral of the square of the vertical current density. We find that the vertical electric current densities have common characteristics for the extended bipolar (beta) (gamma) (delta)-regions studied. The largest current systems have j(sub z)'s which maximizes around 30 mA/sq m and have a linear decreasing distribution to a diameter of 30 Mn.

Gary, G. Allen↗
NASA NTRS1996

The Contradiction Between the Measurement Theory of Quantum Mechanics and the Theory that the Velocity of Any Particle Can Not be Larger than the Velocity of Light

By the measurement theory of quantum mechanics and the method of Fourier transform,we proved that the wave function psi(x,y,z,t)= (8/((2(pi)(2L(exp (1/2)))(exp 3))(Phi(L,t,x)Phi(L,t,y)Phi(L,t,z)). According to the theory that the velocity of any particle can not be larger than the velocity of light and the Born interpretation, when absolute value of delta greater than (ct+ L),Phi(L,t,delta) = 0. But according to the calculation, we proved that for some delta, even if absolute value of delta is greater than (ct+L), Phi(L,t,delta) is not equal to 0.

Shen, Y.↗
NASA NTRS2002

Progress Toward Generation of a Navier-Stokes Database for a Harrier in Ground Effect

The Harrier YAV-8B aircraft is capable of vertical and short-field take-off and landing (V/STOL) by directing its four exhaust nozzles toward the ground, or conventional flight by rotating its nozzles into a horizontal position. The British Royal Air Force and the United States Marine Corps have used this aircraft for more than 30 years to provide a quick reaction time for troop support, and reduce the need for long runways. The success of this powered-lift (PL) vehicle has also prompted the more recent design of the Joint Strike Fighter (JSF). However there are significant safety issues that must be addressed when operating a PL vehicle in close proximity to the ground. Hot Gas Ingestion (HGI) by the inlets can result in a rapid loss of powered lift; and high-speed jet flows along the ground plane can induce low pressures underneath the vehicle, causing a 'suck-down' effect. Under these conditions, departure from controlled flight may occur. Moreover, unsteady ground vortices and jet fountains can affect the aircraft,s controllability and its proximity to ground troops. The viscous, time-dependent flow fields of PL vehicles are difficult to accurately and efficiently predict using Computational Fluid Dynamics (CFD). A number of researchers have used the time-dependent Reynolds-averaged Navier-Stokes (RANS) equations to compute flows for single and multiple jets in a cross-flow. A few have added some geometric complexity to the problem by computing flows for jet-augmented delta wings near a ground plane. Smith et.al. computed for the first time a single RANS solution about a simplified Harrier. This geometry included a fuselage, wing, leading edge root extension (LERX), inlets, and exhaust nozzles. All of these investigations cite two practical problems with computing these flows: 1) the need for improved solution accuracy; and, 2) the need for faster solution methods. We view the need for faster solution methods as key to improving the solution accuracy and making this class of computation more routine. One can hardly refine grids, explore the use of advanced turbulence models, and generate databases when it takes weeks of dedicated computer time for a single solution. Chaderjian, Ahmad, Pandya, and Murman have focused on reducing the time-to-solution for this very difficult and complex problem through process automation and exploitation of parallel computing. They began with the Harrier geometry reported, and added a deflected wing flap and empennage for greater realism. To date more than 80 solutions have been carried out. This paper will describe this process and progress made in reducing the time required to generate a simple longitudinal force and moment database for a Harrier in ground effect. It shows a typical snap-shot from an unsteady streakline animation, where fluid particles are colored by temperature. The ground vortex and a jet-fountain vortex are highlighted. It also shows a similar streakline image, where HGI occurs due to the vehicle in close proximity to the ground. It is show the mean lift coefficient as a function of angle of attack and height. The angle of attack range was 4 deg less than or = alpha less than or = 10 deg with an increment of 1 degree, and the height range was 10 ft less than or = h less than or = 30ft with an increment of 5 feet. This 35 solution database was extended to over 2500 cases using a monotone cubic-spline interpolation procedure. The suck-down effect (reduction of lift near the ground) is highlighted in the figure. The "cushion effect," the conventional reduction of lift as the vehicle moves out of ground effect, is also indicated. All 35 RANS solutions were obtained using 952 Silicon Graphics Origin 2000 and 3000 processors in dedicated mode for one week. Typically, 112 processors were assigned to each case. Some other cases used fewer processors to utilize all available CPUS. The final paper will report on the automation of the solution process, including: grid generation, job monitoring, solution completion criteria, and post processing. Moreover, improvements in parallel efficiency for a dual time-step algorithm for the RANS equations will also be presented. Results will be discussed in detail using unsteady streakline flow visualization to correlate unsteady flow structures with dominant aerodynamic frequencies. The stability derivatives, CL, and CL, will also be presented.

Chaderjian, Neal M.↗
NASA NTRS1990

Vortex Flow Aerodynamics

The AGARD Fluid Dynamics Panel sponsored this symposium to provide an update review of the aerodynamic effects of separation-induced vortex flows on the design and off-design performance of fighter aircraft. A better understanding is needed to predict and control these vortex flows throughout the flight envelope at subsonic, transonic, and supersonic speeds, and especially during high lift operations for take-offs, landings, and sustained and instantaneous maneuvers. Topics of interest included vortex development and burst, modeling and validation of the full range of analytic methods, slender-body vortex flows at high angles of attack, vortex control and management, and unsteady vortex flow effects. In addition to the symposium papers, the results of the round table discussion are presented.

delta wings↗
NASA NTRS1992

Three-dimensional time-marching aeroelastic analyses using an unstructured-grid Euler method

Modifications to a three dimensional, implicit, upwind, unstructured-grid Euler code for aeroelastic analysis of complete aircraft configurations are described. The modifications involve the addition of the structural equations of motion for their simultaneous time integration with the governing flow equations. The paper presents a detailed description of the time marching aeroelastic procedure and presents comparisons with experimental data to provide an assessment of the capability. Flutter results are shown for an isolated 45 degree swept-back wing and a supersonic transport configuration with a fuselage, clipped delta wing, and two identical rearward-mounted nacelles. Comparisons between computed and experimental flutter characteristics show good agreement, giving confidence in the accuracy of the aeroelastic capability that was developed.

Rausch, Russ D.↗