Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical differentiation”

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 217 records · Page 12

Goertler vortices in supersonic boundary layers

The problem of Goertler vortices in compressible boundary layers over concave walls is studied by solving parabolic partial differential equations. The numerical scheme used is based upon a fourth-order accurate compact difference scheme. The results show that both the velocity and temperature fluctuations may lead to a Goertler vortex. The vortex growth rates determined from the present method are found to differ somewhat from those given by normal mode solution. The effects of cooling, heating and pressure gradient are also studied. At a Mach number of 3.5, cooling has a small destabilizing effect. The effect of adverse pressure gradient is also found to be destabilizing on a Mach 8 boundary layer.

Spall, R. E.↗

GPS-aided gravimetry at 30 km altitude from a balloon-borne platform

A balloon-borne experiment, flown at 30 km altitude over New Mexico, was used to test dynamic differential Global Positioning System (GPS) tracking in support of gravimetry at high-altitudes. The experiment package contained a gravimeter (Vibrating String Accelerometer), a full complement of inertial instruments, a TI-4100 GPS receiver and a radar transponder. The flight was supported by two GPS receivers on the ground near the flight path. From the 8 hour flight, about a forty minute period was selected for analysis. Differential GPS phase measurements were used to estimate changes in position over the sample time interval, or average velocity. In addition to average velocity, differential positions and numerical averages of acceleration were obtained in three components. Gravitational acceleration was estimated by correcting for accelerations due to translational motion, ignoring all rotational effects.

Lazarewicz, Andrew R.↗

A new second-order integration algorithm for simulating mechanical dynamic systems

A new integration algorithm which has the simplicity of Euler integration but exhibits second-order accuracy is described. In fixed-step numerical integration of differential equations for mechanical dynamic systems the method represents displacement and acceleration variables at integer step times and velocity variables at half-integer step times. Asymptotic accuracy of the algorithm is twice that of trapezoidal integration and ten times that of second-order Adams-Bashforth integration. The algorithm is also compatible with real-time inputs when used for a real-time simulation. It can be used to produce simulation outputs at double the integration frame rate, i.e., at both half-integer and integer frame times, even though it requires only one evaluation of state-variable derivatives per integration step. The new algorithm is shown to be especially effective in the simulation of lightly-damped structural modes. Both time-domain and frequency-domain accuracy comparisons with traditional integration methods are presented. Stability of the new algorithm is also examined.

Howe, R. M.↗

Radiation transport around axisymmetric blunt body vehicles using a modified differential approximation

A moment method for computing 3D radiative transport in axisymmetric thermochemical nonequilibrium flows is developed. The method uses the P-1 approximation to reduce the governing system of integro-differential equations to a coupled set of partial differential equations. The numerical solution of these equations for realistic variations of the radiation properties is discussed. Representative results from the method are shown and compared to tangent slab calculations. The agreement between the transport methods is found to be about 10 percent in the stagnation region, with the difference increasing along the flank of the vehicle.

Hartung, Lin C.↗

IMP: A performance code

IMP (Integrated Mission Program) is a simulation language and code used to model present and future Earth, Moon, or Mars missions. The profile is user controlled through selection from a large menu of events and maneuvers. A Fehlberg 7/13 Runge-Kutta integrator with error and step size control is used to numerically integrate the differential equations of motion (DEQ) of three spacecraft, a main, a target, and an observer. Through selection, the DEQ's include guided thrust, oblate gravity, atmosphere drag, solar pressure, and Moon gravity effects. Guide parameters for thrust events and performance parameters of velocity changes (Delta-V) and propellant usage (maximum of five systems) are developed as needed. Print, plot, summary, and debug files are output.

Dauro, Vincent A., Sr.↗

Landstat 7 Views the World

The technical innovations of precise numerical radiometry, spectral differentiation and seasonally repetitive monitoring, introduced by the Landsat series of satellites, created a new measurement approach for terrestrial monitoring. The Landsat 7 mission, successfully initiated on April 15, 1999, reflects significant progress in realizing the scientific potential of this measurement strategy. Substantial improvements in calibration procedures, both prior to launch and during normal operations, have been accomplished to insure long-term stability in the acquired spectral radiometry. For the first time, Landsat 7 will be operated with a long-term data acquisition plan that will insure, barring national emergencies, that substantially cloud-free, seasonal coverage will be recorded and archived in the US for all land areas of the globe. The expected outcome of these efforts is a rapid improvement in our understanding of the Earth system, as well as conceptual knowledge that will underpin commercial application of this technology. The Cape Town 2000 Symposium will take place approximately nine months after the Landsat 7 satellite and instrument completed its engineering checkout phase and was placed in its final orbit in late June, 1999. Thus, this presentation will provide a timely status report, in a significant international setting, regarding the early on-orbit performance of the Landsat 7 spacecraft, the ETM+ instrument, and the long-term data acquisition strategy that has been implemented. The latest data display and visualization techniques will be utilized to present this significant new source of Earth observation data.

Williams, Darrel L.↗

Large Field Visualization with Demand-Driven Calculation

We present a system designed for the interactive definition and visualization of fields derived from large data sets: the Demand-Driven Visualizer (DDV). The system allows the user to write arbitrary expressions to define new fields, and then apply a variety of visualization techniques to the result. Expressions can include differential operators and numerous other built-in functions, ail of which are evaluated at specific field locations completely on demand. The payoff of following a demand-driven design philosophy throughout becomes particularly evident when working with large time-series data, where the costs of eager evaluation alternatives can be prohibitive.

Moran, Patrick J.↗

Meteorite Dunite Breccia MIL 03443: A Probable Crustal Cumulate Closely Related to Diogenites from the HED Parent Asteroid

There are numerous types of differentiated meteorites, but most represent either the crusts or cores of their parent asteroids. Ureilites, olivine-pyroxene-graphite rocks, are exceptions; they are mantle restites [1]. Dunite is expected to be a common mantle lithology in differentiated asteroids. In particular, models of the eucrite parent asteroid contain large volumes of dunite mantle [2-4]. Yet dunites are very rare among meteorites, and none are known associated with the howardite, eucrite, diogenite (HED) suite. Spectroscopic measurements of 4 Vesta, the probable HED parent asteroid, show one region with an olivine signature [5] although the surface is dominated by basaltic and orthopyroxenitic material equated with eucrites and diogenites [6]. One might expect that a small number of dunitic or olivine-rich meteorites might be delivered along with the HED suite. The 46 gram meteoritic dunite MIL 03443 (Fig. 1) was recovered from the Miller Range ice field of Antarctica. This meteorite was tentatively classified as a mesosiderite because large, dunitic clasts are found in this type of meteorite, but it was noted that MIL 03443 could represent a dunite sample of the HED suite [7]. Here I will present a preliminary petrologic study of two thin sections of this meteorite.

Mittlefehldt, David W.↗

A Survey of New Trends in Symbolic Execution for Software Testing and Analysis

Symbolic execution is a well-known program analysis technique which represents values of program inputs with symbolic values instead of concrete (initialized) data and executes the program by manipulating program expressions involving the symbolic values. Symbolic execution has been proposed over three decades ago but recently it has found renewed interest in the research community, due in part to the progress in decision procedures, availability of powerful computers and new algorithmic developments. We provide a survey of some of the new research trends in symbolic execution, with particular emphasis on applications to test generation and program analysis. We first describe an approach that handles complex programming constructs such as input data structures, arrays, as well as multi-threading. We follow with a discussion of abstraction techniques that can be used to limit the (possibly infinite) number of symbolic configurations that need to be analyzed for the symbolic execution of looping programs. Furthermore, we describe recent hybrid techniques that combine concrete and symbolic execution to overcome some of the inherent limitations of symbolic execution, such as handling native code or availability of decision procedures for the application domain. Finally, we give a short survey of interesting new applications, such as predictive testing, invariant inference, program repair, analysis of parallel numerical programs and differential symbolic execution.

Pasareanu, Corina S.↗

Smart Cameras for Remote Science Survey

Communication with remote exploration spacecraft is often intermittent and bandwidth is highly constrained. Future missions could use onboard science data understanding to prioritize downlink of critical features [1], draft summary maps of visited terrain [2], or identify targets of opportunity for followup measurements [3]. We describe a generic approach to classify geologic surfaces for autonomous science operations, suitable for parallelized implementations in FPGA hardware. We map these surfaces with texture channels - distinctive numerical signatures that differentiate properties such as roughness, pavement coatings, regolith characteristics, sedimentary fabrics and differential outcrop weathering. This work describes our basic image analysis approach and reports an initial performance evaluation using surface images from the Mars Exploration Rovers. Future work will incorporate these methods into camera hardware for real-time processing.

robotics↗

Experimental Simulations of Lunar Magma Ocean Crystallization: The Plot (But Not the Crust) Thickens

Numerical models of differentiation of a global-scale lunar magma ocean (LMO) have raised as many questions as they have answered. Recent orbital missions and sample studies have provided new context for a large range of lithologies, from the comparatively magnesian "purest anorthosite" reported by to Si-rich domes and spinel-rich clasts with widespread areal distributions. In addition, the GRAIL mission provided strong constraints on lunar crustal density and average thickness. Can this increasingly complex geology be accounted for via the formation and evolution of the LMO? We have in recent years been conducting extensive sets of petrologic experiments designed to fully simulate LMO crystallization, which had not been attempted previously. Here we review the key results from these experiments, which show that LMO differentiation is more complex than initial models suggested. Several important features expected from LMO crystallization models have yet to be reproduced experimentally; combined modelling and experimental work by our group is ongoing.

Draper, D. S.↗

Benchmark Performance Metrics of a Vane Propellant Management Device for a 0.15 m3 Liquid Hydrogen Tank

The use of cryogenic propellants has and will continue to play an integral role in manned-space exploration due the high specific impulses offered and its ubiquity through in-situ resource utilization. But guaranteeing vapor-free transfer of such low-surface tension liquids is difficult for traditional capillary-action propellant management devices (PMDs). Screen-channel liquid acquisition devices and compliant origami bladders are potential solutions, but to quantify the benefits these technologies offer, this paper present a cast study analyzing the performance metrics of an orthodox vane PMD as a benchmark for comparison. A 0.15 m3 liquid hydrogen tank at 20.3 K and 103 kPa was selected for study. Then assuming no body forces and no heat transfer (for simplicity), a steady-state 1-D differential equation was numerically solved in tandem with four possible wetted area configurations to yield an expulsion efficiency for a theoretically maximum inputted expulsion flow rate. Maximization of the flow rate was constrained by the onset of choked flow. Additional inputs, including vane height and vane number, were parametrically varied between 0.1 – 10 cm and 4 – 28 vanes, respectively.

Liquid Acquisition Device↗

Runge-Kutta Methods for Linear Ordinary Differential Equations

Three new Runge-Kutta methods are presented for numerical integration of systems of linear inhomogeneous ordinary differential equations (ODES) with constant coefficients. Such ODEs arise in the numerical solution of the partial differential equations governing linear wave phenomena. The restriction to linear ODEs with constant coefficients reduces the number of conditions which the coefficients of the Runge-Kutta method must satisfy. This freedom is used to develop methods which are more efficient than conventional Runge-Kutta methods. A fourth-order method is presented which uses only two memory locations per dependent variable, while the classical fourth-order Runge-Kutta method uses three. This method is an excellent choice for simulations of linear wave phenomena if memory is a primary concern. In addition, fifth- and sixth-order methods are presented which require five and six stages, respectively, one fewer than their conventional counterparts, and are therefore more efficient. These methods are an excellent option for use with high-order spatial discretizations.

Zingg, David W.↗

Optimal moving grids for time-dependent partial differential equations

Various adaptive moving grid techniques for the numerical solution of time-dependent partial differential equations were proposed. The precise criterion for grid motion varies, but most techniques will attempt to give grids on which the solution of the partial differential equation can be well represented. Moving grids are investigated on which the solutions of the linear heat conduction and viscous Burgers' equation in one space dimension are optimally approximated. Precisely, the results of numerical calculations of optimal moving grids for piecewise linear finite element approximation of partial differential equation solutions in the least squares norm.

Wathen, A. J.↗