Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel simulation”

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 1,063 records · Page 59

The Production of Jets From Magnetic Accretion Disks: Simulation of the Blandford-Payne Mechanism

We have performed magnetohydrodynamic (MRD) simulations of the production of jets from magnetized accretion disks with a factor of 5 greater extent in space and time, and with more models, than any study published so far. We find that jets are produced by such disks in a broad range of parameter space, and by at least two different mechanisms. We also are able to follow the propagation of the jet well beyond the accretion disk into the region of hydrodynamic collimation. The code used is our MHD simulation code FLOW (K. Lind, D. Payne, D. Meier, and R. Blandford, 1989), converted to run on Caltech's massively parallel Intel Touchstone Delta supercomputer. Some of these models may be directly applicable to observed radio sources.

Blandford-Payne↗

Large Amplitude Whistler Waves and Electron Acceleration in the Earth's Radiation Belts: A Review of STEREO and Wind Observations

One of the critical problems for understanding the dynamics of Earth's radiation belts is determining the physical processes that energize and scatter relativistic electrons. We review measurements from the Wind/Waves and STEREO S/Waves waveform capture instruments of large amplitude whistler-mode waves. These observations have provided strong evidence that large amplitude (100s mV/m) whistler-mode waves are common during magnetically active periods. The large amplitude whistlers have characteristics that are different from typical chorus. They are usually nondispersive and obliquely propagating, with a large longitudinal electric field and significant parallel electric field. We will also review comparisons of STEREO and Wind wave observations with SAMPEX observations of electron microbursts. Simulations show that the waves can result in energization by many MeV and/or scattering by large angles during a single wave packet encounter due to coherent, nonlinear processes including trapping. The experimental observations combined with simulations suggest that quasilinear theoretical models of electron energization and scattering via small-amplitude waves, with timescales of hours to days, may be inadequate for understanding radiation belt dynamics.

Cattell, Cynthia↗

The Safe Removal of Frozen Air from the Annulus of an LH2 Storage Tank

Large Liquid Hydrogen (LH2) storage tanks are vital infrastructure for NASA. Eventually, air may leak into the evacuated and perlite filled annular region of these tanks. Although the vacuum level is monitored in this region, the extremely cold temperature causes all but the helium and neon constituents of air to freeze. A small, often unnoticeable pressure rise is the result. As the leak persists, the quantity of frozen air increases, as does the thermal conductivity of the insulation system. Consequently, a notable increase in commodity boil-off is often the first indicator of an air leak. Severe damage can result from normal draining of the tank. The warming air will sublimate which will cause a pressure rise in the annulus. When the pressure increases above the triple point, the frozen air will begin to melt and migrate downward. Collection of liquid air on the carbon steel outer shell may chill it below its ductility range, resulting in fracture. In order to avoid a structural failure, as described above, a method for the safe removal of frozen air is needed. A thermal model of the storage tank has been created using SINDA/FLUINT modeling software. Experimental work is progressing in an attempt to characterize the thermal conductivity of a perlite/frozen nitrogen mixture. A statistical mechanics model is being developed in parallel for comparison to experimental work. The thermal model will be updated using the experimental/statistical mechanical data, and used to simulate potential removal scenarios. This paper will address methodologies and analysis techniques for evaluation of two proposed air removal methods.

cryogenic insulation↗

Development of a Grid-Independent Geos-Chem Chemical Transport Model (v9-02) as an Atmospheric Chemistry Module for Earth System Models

The GEOS-Chem global chemical transport model (CTM), used by a large atmospheric chemistry research community, has been re-engineered to also serve as an atmospheric chemistry module for Earth system models (ESMs). This was done using an Earth System Modeling Framework (ESMF) interface that operates independently of the GEOSChem scientific code, permitting the exact same GEOSChem code to be used as an ESM module or as a standalone CTM. In this manner, the continual stream of updates contributed by the CTM user community is automatically passed on to the ESM module, which remains state of science and referenced to the latest version of the standard GEOS-Chem CTM. A major step in this re-engineering was to make GEOS-Chem grid independent, i.e., capable of using any geophysical grid specified at run time. GEOS-Chem data sockets were also created for communication between modules and with external ESM code. The grid-independent, ESMF-compatible GEOS-Chem is now the standard version of the GEOS-Chem CTM. It has been implemented as an atmospheric chemistry module into the NASA GEOS- 5 ESM. The coupled GEOS-5-GEOS-Chem system was tested for scalability and performance with a tropospheric oxidant-aerosol simulation (120 coupled species, 66 transported tracers) using 48-240 cores and message-passing interface (MPI) distributed-memory parallelization. Numerical experiments demonstrate that the GEOS-Chem chemistry module scales efficiently for the number of cores tested, with no degradation as the number of cores increases. Although inclusion of atmospheric chemistry in ESMs is computationally expensive, the excellent scalability of the chemistry module means that the relative cost goes down with increasing number of cores in a massively parallel environment.

CTM↗

Applying Orbital Multi-Angle Photopolarimetric Observations to Study Properties of Aerosols in the Earth's Atmosphere: Implications of Measurements in the 1.378 µm Spectral Channel to Retrieve Microphysical Characteristics and Composition of Stratospheric Aerosols

We analyze the possibilities of orbital photopolarimetric measurements to study properties of aerosols in the Earth's atmosphere. As an example, we consider the case when such measurements are performed within a narrow spectral channel centered at 1.378 µm that allows to retrieve microphysical characteristics of stratospheric aerosols separately from those of tropospheric aerosols. We consider the case of stratospheric aerosols caused by volcanic eruption, and adopt the model of the stratosphere in the form of a homogeneous plane-parallel layer composed of polydisperse spherical particles. We use numerically exact solutions of the vector radiative transfer equation to theoretically simulate measurements carried out at various numbers of scattering angles, including: (i) radiance measurements alone; (ii) polarization measurements alone; and (iii) radiance and polarization measurements together. The results of computations show that the simultaneous use of radiance and polarization measurements at a sufficiently large number of scattering angles enables one to retrieve the optical thickness, effective radius, and refractive index of aerosols with adequate accuracy. We demonstrate how the accuracy of the derived values of the optical parameters of aerosols depends on the accuracy of measurements of the intensity and polarization of the reflected light, optical thickness of aerosol layer itself, effective radius of aerosols, width of the particle size distribution, and number of viewing angles.

Aerosols↗

NASA Langley FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

This work summarizes the results of FUN3D analyses conducted for the 1st AIAA Stability and Control Workshop on behalf of participants from the NASA Langley Research Center. The workshop was created to establish best practices for the prediction of stability and control derivatives using computational fluid dynamics and assess the limitations of these methods when those best practices are applied. The inaugural workshop considered the ONERA version of the NASA/Boeing Common Research Model, which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel data at small sideslip angles remain unpublished and served as ‘blind’ data for computational comparisons. The present research generated workshop test case data using the NASA FUN3D code, which is a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Steady- state numerical simulations were conducted for workshop test cases investigating the following: grid convergence, Mach number effect on static stability, wind tunnel sting increments, static stability-derivative calculations, and a sideslip angle sweep. Results were generated for two series of unstructured, mixed-element grids, one set provided by the workshop and another set created using the HeldenMesh grid generation software. The results provided include total- and component-level breakdowns of the force and moment coefficients, in addition to sectional pressure distributions for the wing and tail components for comparisons to wind tunnel data.

CFD↗

Spacecraft contamination investigation by direct simulation Monte Carlo - Contamination on UARS/HALOE

A three-dimensional version of the direct simulation Monte Carlo method is adapted to study the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite, with the ultimate goal of estimating the cumulative level of molecular deposits on critical optical surfaces of the Halogen Occultation Experiment. Emphasis is placed on describing the development of preprocessing codes that simplify the incorporation of geometry and outgassing data into the main run code. A novel use of parallel processing software allows use of multiple processors to solve different regions of the overall flowfield. Initial results relating to contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the Halogen Occultation telescope's aperture plane are presented, along with data related to code performance.

Rault, Didier F. G.↗

Oblique Instability of Quasi-Parallel Whistler Waves in the Presence of Cold and Warm Electron Populations

Whistler waves propagating nearly parallel to the ambient magnetic field experience a nonlinear instability due to transverse currents when the background plasma has a population of sufficiently low energy electrons. Intriguingly, this nonlinear process may generate oblique electrostatic waves, including whistlers near the resonance cone with properties resembling oblique chorus waves in the Earth’s magnetosphere. Focusing on the generation of oblique whistlers, earlier analysis of the instability is extended here to the case where low-energy background plasma consists of both a “cold” population with energy of a few eV and a “warm” electron component with energy of the order of 100 eV. This is motivated by spacecraft observations in the Earth’s magnetosphere where oblique chorus waves were shown to interact resonantly with the warm electrons. The main new results are: 1) the instability producing oblique electrostatic waves is sensitive to the shape of the electron distribution at low energies. In the whistler range of frequencies, two distinct peaks in the growth rate are typically present for the model considered: a peak associated with the warm electron population at relatively low wavenumbers and a peak associated with the cold electron population at relatively high wavenumbers; 2) overall, the instability producing oblique whistler waves near the resonance cone persists (with a reduced growth rate) even in the cases where the temperature of the cold population is relatively high, including cases where cold population is absent and only the warm population is included; 3) particle-in-cell simulations show that the instability leads to heating of the background plasma and formation of characteristic plateau and beam features in the parallel electron distribution function in the range of energies resonant with the instability. The plateau/beam features have been previously detected in spacecraft observations of oblique chorus waves. However, they have been attributed to external sources and have been proposed to be the mechanism generating oblique chorus. In the present scenario, the causality link is reversed and the instability generating oblique whistler waves is shown to be a possible mechanism for formation of the plateau and beam features.

Vadim Roytershteyn↗

Molecular Dynamics Simulation of Carbon Nanotube Based Gears

We used molecular dynamics to investigate the properties and design space of molecular gears fashioned from carbon nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. One gear was powered by forcing the atoms near the end of the buckytube to rotate, and a second gear was allowed.to rotate by keeping the atoms near the end of its buckytube on a cylinder. The meshing aromatic gear teeth transfer angular momentum from the powered gear to the driven gear. A number of gear and gear/shaft configurations were simulated. Cases in vacuum and with an inert atmosphere were examined. In an extension to molecular dynamics technology, some simulations used a thermostat on the atmosphere while the hydrocarbon gear's temperature was allowed to fluctuate. This models cooling the gears with an atmosphere. Results suggest that these gears can operate at up to 50-100 gigahertz in a vacuum or inert atmosphere at room temperature. The failure mode involves tooth slip, not bond breaking, so failed gears can be returned to operation by lowering temperature and/or rotation rate. Videos and atomic trajectory files in xyz format are presented.

Han, Jie↗

Structural, Physical, and Compositional Analysis of Lunar Simulants and Regolith

Relative to the prior manned Apollo and unmanned robotic missions, planned Lunar initiatives are comparatively complex and longer in duration. Individual crew rotations are envisioned to span several months, and various surface systems must function in the Lunar environment for periods of years. As a consequence, an increased understanding of the surface environment is required to engineer and test the associated materials, components, and systems necessary to sustain human habitation and surface operations. The effort described here concerns the analysis of existing simulant materials, with application to Lunar return samples. The interplay between these analyses fulfills the objective of ascertaining the critical properties of regolith itself, and the parallel objective of developing suitable stimulant materials for a variety of engineering applications. Presented here are measurements of the basic physical attributes, i.e. particle size distributions and general shape factors. Also discussed are structural and chemical properties, as determined through a variety of techniques, such as optical microscopy, SEM and TEM microscopy, Mossbauer Spectroscopy, X-ray diffraction, Raman microspectroscopy, inductively coupled argon plasma emission spectroscopy and energy dispersive X-ray fluorescence mapping. A comparative description of currently available stimulant materials is discussed, with implications for more detailed analyses, as well as the requirements for continued refinement of methods for simulant production.

Greenberg, Paul↗

Better, Cheaper, Faster Molecular Dynamics

Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.

Pohorille, Andrew↗

Numerical simulations of double layers and auroral electric fields

Recent one-dimensional and two-dimensional numerical simulations of double layers (DLs) in the electric fields of the auroral plasma are reviewed, with reference to observational data. It is found that two-dimensional DLs driven by current sheets of finite thickness have different characteristics, depending on whether the layer thickness is less than or much greater than the ion gyroradius: When thickness is less than ion gyroradius, V-shaped DLs form with nearly equal parallel and perpendicular potential drops; when layer thickness is much greater than ion gyroradius the major parallel potential drop occurs outside the current sheet and the perpendicular electric fields are localized at the edges of the current sheet. It is shown that some features of the simulated fields, such as the amplitudes and scale lengths, are qualitatively similar to those observed in space.

Singh, N.↗

Predictions of lithium interactions with earth's bow shock in the presence of wave activity

The results of a test-particle simulation studying the movement of a lithium tracer ion injected upstream of the bow shock are reported. Wave activity consists of parallel and antiparallel propagating Alfven waves characterized by a frequency power spectrum within a frequency or range of amplitudes defined separately in the upstream and downstream regions. The results show that even a moderate level of wave activity can substantially change the results obtained in the absence of waves. Among the effects observed are: (1) increased ion transmission; (2) both the average energy gain and spread about the average are increased for transmitted and reflected particles; (3) the average final pitch angle for transmitted particles tends to 90 deg, and the spread of reflected particles is reduced; and (4) the spatial dispersion of the ions on the bow shock after a single encounter is increased.

Decker, R. B.↗

User's manual for the two-dimensional transputer graphics toolkit

The user manual for the 2-D graphics toolkit for a transputer based parallel processor is presented. The toolkit consists of a package of 2-D display routines that can be used for the simulation visualizations. It supports multiple windows, double buffered screens for animations, and simple graphics transformations such as translation, rotation, and scaling. The display routines are written in occam to take advantage of the multiprocessing features available on transputers. The package is designed to run on a transputer separate from the graphics board.

Ellis, Graham K.↗

General-Purpose Data-Formatting Input/Output System

Multiplexable input/output (I/O) system developed as interface between host computer and real-time, fixed-base simulator cockpit. Unit designed to be general-purpose interface and operated through virtually any 8- or 16-bit, transistor/transistor-logic-level, parallel-I/O port with or without handshaking. Accepts byte-coded data in form of data blocks from host computer and routes data to various system I/O modules like discrete outputs, lamp drivers, and seven-segment-display drivers. Formats and routes data from input modules back to host computer. Other modules designed to drive alphanumeric displays and communications interfaces and to serve as analog-to-digital and digital-to-analog converter modules.

Busquets, Anthony M.↗

A further study of cumulus interactions and mergers - Three-dimensional simulations with trajectory analysis

Cumulus processes involved in the interaction and merging of clouds under the influence of different imposed conditions (including large-scale lifting forcing, environmental wind shear, and cloud microphysical processes) were studied using simulations with a three-dimensional model. The design of the study was to generate several convective clouds randomly inside the model domain, and then to observe and analyze the interactions and merging between the simulated clouds. Ten merged clouds were identified. Seven of these, each involving two previously separated clouds, generally lie along a line parallel to the initial environmental wind shear vector, while one (also a two-cloud system) lies along a line perpendicular to the wind shear vector prior to merging. The remaining two merging systems involve three parent clouds each; they are a combination of parallel and perpendicular cells. The merging mechanisms associated with three-cloud merging cases are studied by examining the temperature, pressure, and wind fields prior to, during, and following the merging.

Tao, Wei-Kuo↗

Space applications of artificial intelligence; Proceedings of the Annual Goddard Conference, Greenbelt, MD, May 16, 17, 1989

Theoretical and implementation aspects of AI systems for space applications are discussed in reviews and reports. Sections are devoted to planning and scheduling, fault isolation and diagnosis, data management, modeling and simulation, and development tools and methods. Particular attention is given to a situated reasoning architecture for space repair and replace tasks, parallel plan execution with self-processing networks, the electrical diagnostics expert system for Spacelab life-sciences experiments, diagnostic tolerance for missing sensor data, the integration of perception and reasoning in fast neural modules, a connectionist model for dynamic control, and applications of fuzzy sets to the development of rule-based expert systems.

Rash, James L.↗

A survey of parallel programming tools

This survey examines 39 parallel programming tools. Focus is placed on those tool capabilites needed for parallel scientific programming rather than for general computer science. The tools are classified with current and future needs of Numerical Aerodynamic Simulator (NAS) in mind: existing and anticipated NAS supercomputers and workstations; operating systems; programming languages; and applications. They are divided into four categories: suggested acquisitions, tools already brought in; tools worth tracking; and tools eliminated from further consideration at this time.

Cheng, Doreen Y.↗