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 271 records · Page 15

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↗

Design of a Combined Proportional Integral Derivative Controller to Regulate the Temperature Inside a High-Temperature Tubular Solar Reactor

Solar fuels are proven to be promising candidates for thermochemical energy storage. However, the transient nature of solar radiation is an obstacle to maintaining a stable operational temperature inside a solar reactor. To overcome this challenge, the temperature of a solar reactor can be regulated by controlling the incoming solar radiation or the feedstock flowrate inside the reactor. In this work, a combined proportional integral derivative (PID) controller is implemented to regulate the temperature inside a high-temperature tubular solar reactor with counter-current flowing gas/particles. The control model incorporates two control systems to regulate incoming solar radiation and gas flow simultaneously. The design of the controller is based on a reduced-order numerical model of a high-temperature tubular solar reactor that is vertically oriented with an upward gas flow and downward particle flow. The reactor receives heat circumferentially through its wall over a finite segment of its length. Formulation of the heat transfer model is presented by applying the energy balance for the reactor tube and considering heat and mass transfer inside. A set of governing differential equations are solved numerically by using the finite volume method to obtain reactor wall, particles, and gas temperatures along the reactor length with various boundary conditions. Simulation results are used to tune the PID controller parameters by utilizing the Ziegler–Nichols tuning method. Both the simulation results and the controller performance are visualized on the LABVIEW platform. The controller is challenged to track different temperature setpoints with different scenarios of transient solar radiation. Here, the performance of the PID controller was compared to experimental results obtained from an industrial PID controller embedded in a 7 kW electric furnace. Results show that the combined PID controller is successful in maintaining a stable temperature inside the reactor by regulating the incoming solar radiation and the flowrate via small steady-state error and reasonable settling time and overshoot.

14 SOLAR ENERGY↗

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.↗

Scalability of high-performance PDE solvers

Performance tests and analyses are critical to effective high-performance computing software development and are central components in the design and implementation of computational algorithms for achieving faster simulations on existing and future computing architectures for large-scale application problems. In this article, we explore performance and space-time trade-offs for important compute-intensive kernels of large-scale numerical solvers for partial differential equations (PDEs) that govern a wide range of physical applications. We consider a sequence of PDE-motivated bake-off problems designed to establish best practices for efficient high-order simulations across a variety of codes and platforms. We measure peak performance (degrees of freedom per second) on a fixed number of nodes and identify effective code optimization strategies for each architecture. In addition to peak performance, we identify the minimum time to solution at 80% parallel efficiency. The performance analysis is based on spectral and p-type finite elements but is equally applicable to a broad spectrum of numerical PDE discretizations, including finite difference, finite volume, and h-type finite elements.

97 MATHEMATICS AND COMPUTING↗

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 PDE solutions in the least-squares norm are reported.

Wathen, A. J.↗

Analysis of the fluid flow and heat transfer in a thin liquid film in the presence and absence of gravity

The hydrodynamic and thermal behavior of a thin liquid film flowing over a solid horizontal surface is analyzed for both plane and radially spreading flows. The situations where the gravitational force is completely absent and where it is significant are analyzed separately and their practical relevance to a micro-gravity environment is discussed. In the presence of gravity, in addition to Reynolds number, the Froude number of the film is found to be an important parameter that determines the supercritical and subcritical flow regimes and any associated hydraulic jump. A closed-form solution is possible under some flow situations, whereas others require numerical integration of ordinary differential equations. The approximate analytical results are found to compare well with the available two-dimensional numerical solutions.

Rahman, M. M.↗

Period doubling and chaos in partial differential equations for thermosolutal convection

Numerical experiments on two-dimensional thermosolutal convection reveal a transition from periodic oscillations to chaos through a sequence of period-doubling bifurcations. Within the chaotic region there are narrow periodic windows. This is the first example of period-doubling in solutions of partial differential equations. A truncated model indicates that this behavior is associated with heteroclinic explosions.

Moore, D. R.↗