Search NASA⌕ Search

SEARCH · Search NASA

Results for “computer experiments”

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

An implicit and stiffly stable finite element CFD algorithm for unsteady aerodynamics

A stable and accurate finite element CFD algorithm for hyperbolic/incompletely parabolic conservation law systems is described and verified. It combines a Taylor weak statement FEM, an optimal implicit Runge-Kutta time integration algorithm, and a matrix tensor product approximate factorization linear algebra procedure. The results of computational experiments show that the developed algorithm is robust.

Baker, A. J.↗

Minimizing distortion and internal forces in truss structures by simulated annealing

Inaccuracies in the length of members and the diameters of joints of large truss reflector backup structures may produce unacceptable levels of surface distortion and member forces. However, if the member lengths and joint diameters can be measured accurately it is possible to configure the members and joints so that root-mean-square (rms) surface error and/or rms member forces is minimized. Following Greene and Haftka (1989) it is assumed that the force vector f is linearly proportional to the member length errors e(sub M) of dimension NMEMB (the number of members) and joint errors e(sub J) of dimension NJOINT (the number of joints), and that the best-fit displacement vector d is a linear function of f. Let NNODES denote the number of positions on the surface of the truss where error influences are measured. The solution of the problem is discussed. To classify, this problem was compared to a similar combinatorial optimization problem. In particular, when only the member length errors are considered, minimizing d(sup 2)(sub rms) is equivalent to the quadratic assignment problem. The quadratic assignment problem is a well known NP-complete problem in operations research literature. Hence minimizing d(sup 2)(sub rms) is is also an NP-complete problem. The focus of the research is the development of a simulated annealing algorithm to reduce d(sup 2)(sub rms). The plausibility of this technique is its recent success on a variety of NP-complete combinatorial optimization problems including the quadratic assignment problem. A physical analogy for simulated annealing is the way liquids freeze and crystallize. All computational experiments were done on a MicroVAX. The two interchange heuristic is very fast but produces widely varying results. The two and three interchange heuristic provides less variability in the final objective function values but runs much more slowly. Simulated annealing produced the best objective function values for every starting configuration and was faster than the two and three interchange heuristic.

Kincaid, Rex K.↗

M51's spiral structure

The M51 system (NGC 5194/5195) provides an excellent problem both in spiral structure and in galaxy interactions. The authors present an analytic study of a computer experiment on the excitation mechanisms for M51's spiral arms and whether or not a halo is important for these mechanisms. This work extends previous numerical studies of the M51 system by including self-gravitation in a two component disk: gas and stars, and a dark halo. The analytic study provides two new observational constraints: the time (approx. 70 to 84 million years ago) and position angle of perigalacticon (300 degrees). By using these constraints and a simple conic approximation, the search for the companion's possible orbit is greatly simplified. This requires fewer N-body simulations than a fully self-gravitating orbit search.

Howard, S.↗

Structural factoring approach for analyzing stochastic networks

The problem of finding the distribution of the shortest path length through a stochastic network is investigated. A general algorithm for determining the exact distribution of the shortest path length is developed based on the concept of conditional factoring, in which a directed, stochastic network is decomposed into an equivalent set of smaller, generally less complex subnetworks. Several network constructs are identified and exploited to reduce significantly the computational effort required to solve a network problem relative to complete enumeration. This algorithm can be applied to two important classes of stochastic path problems: determining the critical path distribution for acyclic networks and the exact two-terminal reliability for probabilistic networks. Computational experience with the algorithm was encouraging and allowed the exact solution of networks that have been previously analyzed only by approximation techniques.

Hayhurst, Kelly J.↗

Spatial operator approach to flexible multibody system dynamics and control

The inverse and forward dynamics problems for flexible multibody systems were solved using the techniques of spatially recursive Kalman filtering and smoothing. These algorithms are easily developed using a set of identities associated with mass matrix factorization and inversion. These identities are easily derived using the spatial operator algebra developed by the author. Current work is aimed at computational experiments with the described algorithms and at modelling for control design of limber manipulator systems. It is also aimed at handling and manipulation of flexible objects.

Rodriguez, G.↗

Rain waves-wind waves interaction application to scatterometry

Modulation of a rain wave pattern by longer waves has been studied. An analytical model taking into account capillarity effects and obliquity of short waves has been developed. Modulation rates in wave number and amplitude have been computed. Experiments were carried out in a wave tank. First results agree with theoretical models, but higher values of modulation rates are measured. These results could be taken into account for understanding the radar response from the sea surface during rain.

Kharif, C.↗

Magma ocean: Mechanisms of formation

The thermal state of the Earth at the time relevant to formation of a magma ocean was dominated by the great impact that created the Moon. As shown in computer experiments, the iron in the impacting bodies quickly sank to the core of the proto-Earth, while a significant fraction of silicates was pushed far enough out beyond the geosynchronous limit to constitute the main material of the Moon. Most of any atmosphere would have been pushed aside, rather than being expelled in the impact. However, the energy remaining in the material not going to the core or expelled was still sufficient to raise its temperature some 1000's of degrees, enough to vaporize silicates and to generate a strong 'planetary wind': a hydrodynamic expansion carrying with it virtually all volatiles plus appreciable silicates. This expansion was violent and uneven in its most energetic stage, but probably the resulting magma ocean was global. The duration, until cooling, was sufficient for silicates to condense to melt and the duration was probably short. Comparison of the Earth and Venus indicates that the great impact was extraordinarily effective in removing volatiles from the proto-Earth; in particular, the enormous differences in primordial inert gases between the planets demand a catastrophic difference in origin circumstances. On the other hand, the comparison limits the amount of silicates lost by the Earth to a rather minor fraction; most of that expelled in the wind must have condensed soon enough for the silicate to fall back to Earth or be swept up by the proto-Moon. So the Earth was left with a magma ocean. The question is whether sufficient water was retained to constitute a steam atmosphere. Probably not, but unknowns affecting this question are the efficiencies of outgassing in great impacts and in subsequent convective churnings deep in the mantle. During the stage when mantle convection is turbulent, an appreciable fraction of volatiles were also retained at depth, perhaps in some mineral phases not yet well-defined. We still have primordial helium being outgassed.

Kaula, W. M.↗

Laser power beams obtained by the dynamic selection of emitting elements in an array

A method for partially correcting the spatial coherence of an array of laser beams is discussed. The method trades far-field intensity for simplistic implementation and array lifetime extension, resulting in a quasi-coherent array with features of adaptive optics. An algorithm for implementing the method is described and feasible options for improvement are presented. A computer experiment used the algorithm and fast Fourier transforms to calculate typical far-field intensity distributions of such arrays. The far-field distributions illustrate beam power concentration, beam steering, and (indirectly) extended lifetime of the array.

Williams, M. D.↗

High voltage spheres in an unmagnetized plasma - Fluid and PIC simulations

The basic physics involving the interaction of a high voltage sphere with a LEO-Plasma was investigated via computer experiments. Two approaches were used: the fluid description and the more general particle pushing method. Both techniques displayed qualitatively similar features. Thus, the initial time-dependent response of the plasma predicted by the fluid model, including the initial current surge and sheath formation, the formation of ion and electron density shells about the sphere with a double layer potential structure, and the subsequent propagation of the density shells away from the sphere. The PIC results also displayed an oscillatory character and slightly different time constants for the growth and decay mechanisms involved.

Thiemann, H.↗

Structure and dynamics of premixed flames in microgravity

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years with emphasis on the work performed since February 1992, the beginning of the current project. The focus of our research has been on investigating fundamental combustion questions concerning the propagation and extinction of gas-phase flames in microgravity and earth-gravity environments. Our approach to resolving these fundamental questions has been to use detailed time-dependent, multidimensional numerical models to perform carefully designed computational experiments. The basic questions we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein.

Kailasanath, K.↗

Preliminary Studies Leading Toward the Development of a LIDAR Bathymetry Mapping Instrument

The National Aeronautics and Space Administration (NASA) at Goddard Space Flight Center (GSFC) has developed a laser ranging device (LIDAR) which provides accurate and timely data of earth features. NASA/GSFC recently modified the sensor to include a scanning capability to produce LIDAR swaths. They have also integrated a Global Positioning System (GPS) and an Inertial Navigation System (INS) to accurately determine the absolute aircraft location and aircraft attitude (pitch, yaw, and roll), respectively. The sensor has been flown in research mode by NASA for many years. The LIDAR has been used in different configurations or modes to acquire such data as altimetry (topography), bathymetry (water depth), laser-induced fluorosensing (tracer dye movements, oil spills and oil thickness, chlorophyll and plant stress identification), forestry, and wetland discrimination studies. NASA and HARC are developing a commercial version of the instrument for topographic mapping applications. The next phase of the commercialization project will be to investigate other applications such as wetlands mapping and coastal bathymetry. In this paper we report on preliminary laboratory measurements to determine the feasibility of making accurate depth measurements in relatively shallow water (approximately 2 to 6 feet deep) using a LIDAR system. The LIDAR bathymetry measurements are relatively simple in theory. The water depth is determined by measuring the time interval between the water surface reflection and the bottom surface reflection signals. Depth is then calculated by dividing by the index of refraction of water. However, the measurements are somewhat complicated due to the convolution of the water surface return signal with the bottom surface return signal. Therefore in addition to the laboratory experiments, computer simulations of the data were made to show these convolution effects in the return pulse waveform due to: (1) water depth, and (2) changes in bottom surface reflectivity.

Hill, John M.↗

Optimization by decomposition in structural and multidisciplinary applications

An algorithm for a general, multilevel structural optimization by substructuring is derived, based on the linear decomposition concept that is rooted in the Bellman's Optimality Criterion enhanced with the optimum sensitivity derivatives used as a means to account for coupling among the subproblems, each of which is limited to optimization of a substructure. The algorithm applies also to those multidisciplinary problems whose subproblems form a hierarchy similar to that of substructures. In systems where the subproblems communicate with each other at the same level, the decomposition becomes non-hierarchic and the system may be optimized as a whole based on the derivatives of the system behavior with respect to the design variables computed by a method that bypasses finite differencing on the system analysis. When a multidisciplinary system includes a structure as its part, a hybrid, hierarchic/non-hierarchic decomposition applies. Numerical examples and references to computational experience accumulated to date illustrate the discussion.

Sobieszczanski-Sobieski, Jaroslaw↗

Analysis, approximation, and computation of a coupled solid/fluid temperature control problem

An optimization problem is formulated motivated by the desire to remove temperature peaks, i.e., 'hot spots', along the bounding surfaces of containers of fluid flows. The heat equation of the solid container is coupled to the energy equations for the fluid. Heat sources can be located in the solid body, the fluid, or both. Control is effected by adjustments to the temperature of the fluid at the inflow boundary. Both mathematical analyses and computational experiments are given.

Gunzburger, Max D.↗

[STS-41 Onboard 16mm Photography Quick Release]

This videotape features scenes of onboard activities. The videotape was shot by the crew. The scenes include the following: Ulysses' deployment, middeck experiments, computer workstations, and Earth payload bay views.

Source record↗

Introduction to the World Wide Web and Mosaic

This tutorial provides an introduction to some of the terminology related to the use of the World Wide Web and Mosaic. It is assumed that the user has some prior computer experience. References are included to other sources of additional information.

Youngblood, Jim↗

Numerical Conformal Mapping Using Cross-Ratios and Delaunay Triangulation

We propose a new algorithm for computing the Riemann mapping of the unit disk to a polygon, also known as the Schwarz-Christoffel transformation. The new algorithm, CRDT, is based on cross-ratios of the prevertices, and also on cross-ratios of quadrilaterals in a Delaunay triangulation of the polygon. The CRDT algorithm produces an accurate representation of the Riemann mapping even in the presence of arbitrary long, thin regions in the polygon, unlike any previous conformal mapping algorithm. We believe that CRDT can never fail to converge to the correct Riemann mapping, but the correctness and convergence proof depend on conjectures that we have so far not been able to prove. We demonstrate convergence with computational experiments. The Riemann mapping has applications to problems in two-dimensional potential theory and to finite-difference mesh generation. We use CRDT to produce a mapping and solve a boundary value problem on long, thin regions for which no other algorithm can solve these problems.

Driscoll, Tobin A.↗

Cloud identification using genetic algorithms and massively parallel computation

As a Guest Computational Investigator under the NASA administered component of the High Performance Computing and Communication Program, we implemented a massively parallel genetic algorithm on the MasPar SIMD computer. Experiments were conducted using Earth Science data in the domains of meteorology and oceanography. Results obtained in these domains are competitive with, and in most cases better than, similar problems solved using other methods. In the meteorological domain, we chose to identify clouds using AVHRR spectral data. Four cloud speciations were used although most researchers settle for three. Results were remarkedly consistent across all tests (91% accuracy). Refinements of this method may lead to more timely and complete information for Global Circulation Models (GCMS) that are prevalent in weather forecasting and global environment studies. In the oceanographic domain, we chose to identify ocean currents from a spectrometer having similar characteristics to AVHRR. Here the results were mixed (60% to 80% accuracy). Given that one is willing to run the experiment several times (say 10), then it is acceptable to claim the higher accuracy rating. This problem has never been successfully automated. Therefore, these results are encouraging even though less impressive than the cloud experiment. Successful conclusion of an automated ocean current detection system would impact coastal fishing, naval tactics, and the study of micro-climates. Finally we contributed to the basic knowledge of GA (genetic algorithm) behavior in parallel environments. We developed better knowledge of the use of subpopulations in the context of shared breeding pools and the migration of individuals. Rigorous experiments were conducted based on quantifiable performance criteria. While much of the work confirmed current wisdom, for the first time we were able to submit conclusive evidence. The software developed under this grant was placed in the public domain. An extensive user's manual was written and distributed nationwide to scientists whose work might benefit from its availability. Several papers, including two journal articles, were produced.

Buckles, Bill P.↗

Parallel Directionally Split Solver Based on Reformulation of Pipelined Thomas Algorithm

In this research an efficient parallel algorithm for 3-D directionally split problems is developed. The proposed algorithm is based on a reformulated version of the pipelined Thomas algorithm that starts the backward step computations immediately after the completion of the forward step computations for the first portion of lines This algorithm has data available for other computational tasks while processors are idle from the Thomas algorithm. The proposed 3-D directionally split solver is based on the static scheduling of processors where local and non-local, data-dependent and data-independent computations are scheduled while processors are idle. A theoretical model of parallelization efficiency is used to define optimal parameters of the algorithm, to show an asymptotic parallelization penalty and to obtain an optimal cover of a global domain with subdomains. It is shown by computational experiments and by the theoretical model that the proposed algorithm reduces the parallelization penalty about two times over the basic algorithm for the range of the number of processors (subdomains) considered and the number of grid nodes per subdomain.

Povitsky, A.↗