Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer Graphics”

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 811 records · Page 45

An MCNP6 Version of the Visual Editor - Paper 2

The Visual Editor is a graphical user interface for working with the MCNP computer code developed at Los Alamos National Laboratory. It consists of a graphical user interface written in C++ that is linked directory to the MCNP Fortran code. Work on the original Visual Editor started in 1992 and was first published in 1993. The original version of the Visual Editor was created on Unix/Linux operating systems and was linked to the MCNP computer code. The Visual Editor was adapted to the Windows operating system in 1999. In 2008 the Visual Editor was adapted to work with the MCNPX computer code. This paper discusses the latest upgrade of the Visual Editor to be compatible with the MCNP6 computer code. Two graphical user interfaces have been created, the Visual MCNP6 Plotter and the Visual MCNP6 Editor. The Visual MCNP6 Plotter includes all of the standard plotting features available in the default MCNP6 plotting package including geometry plotting, tally plotting, cross section plotting, mesh tally plotting and weight window mesh plotting. The Visual MNCP6 Editor has additional capabilities to allow for the creation of MCNP input files with graphical tools. The Visual MCNP6 Editor allows user to create and modify the input file within the Editor. Modifications are not allowed in the Visual MCNP6 Plotter. Both graphical user interfaces allow for advanced plotting features including 3D geometry plots and particle track plotting. Complete details on these tools along with their current status can be found at the Visual Editor website (www.mcnpvised.com). (author)

61 RADIATION PROTECTION AND DOSIMETRY↗

Towards the clinical implementation of iterative low-dose cone-beam CT reconstruction in image-guided radiation therapy: Cone/ring artifact correction and multiple GPU implementation

Purpose: Compressed sensing (CS)-based iterative reconstruction (IR) techniques are able to reconstruct cone-beam CT (CBCT) images from undersampled noisy data, allowing for imaging dose reduction. However, there are a few practical concerns preventing the clinical implementation of these techniques. On the image quality side, data truncation along the superior–inferior direction under the cone-beam geometry produces severe cone artifacts in the reconstructed images. Ring artifacts are also seen in the half-fan scan mode. On the reconstruction efficiency side, the long computation time hinders clinical use in image-guided radiation therapy (IGRT). Methods: Image quality improvement methods are proposed to mitigate the cone and ring image artifacts in IR. The basic idea is to use weighting factors in the IR data fidelity term to improve projection data consistency with the reconstructed volume. In order to improve the computational efficiency, a multiple graphics processing units (GPUs)-based CS-IR system was developed. The parallelization scheme, detailed analyses of computation time at each step, their relationship with image resolution, and the acceleration factors were studied. The whole system was evaluated in various phantom and patient cases. Results: Ring artifacts can be mitigated by properly designing a weighting factor as a function of the spatial location on the detector. As for the cone artifact, without applying a correction method, it contaminated 13 out of 80 slices in a head-neck case (full-fan). Contamination was even more severe in a pelvis case under half-fan mode, where 36 out of 80 slices were affected, leading to poorer soft tissue delineation and reduced superior–inferior coverage. The proposed method effectively corrects those contaminated slices with mean intensity differences compared to FDK results decreasing from ∼497 and ∼293 HU to ∼39 and ∼27 HU for the full-fan and half-fan cases, respectively. In terms of efficiency boost, an overall 3.1 × speedup factor has been achieved with four GPU cards compared to a single GPU-based reconstruction. The total computation time is ∼30 s for typical clinical cases. Conclusions: The authors have developed a low-dose CBCT IR system for IGRT. By incorporating data consistency-based weighting factors in the IR model, cone/ring artifacts can be mitigated. A boost in computational efficiency is achieved by multi-GPU implementation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗

Alloy Design Workbench-Surface Modeling Package Developed

NASA Glenn Research Center's Computational Materials Group has integrated a graphical user interface with in-house-developed surface modeling capabilities, with the goal of using computationally efficient atomistic simulations to aid the development of advanced aerospace materials, through the modeling of alloy surfaces, surface alloys, and segregation. The software is also ideal for modeling nanomaterials, since surface and interfacial effects can dominate material behavior and properties at this level. Through the combination of an accurate atomistic surface modeling methodology and an efficient computational engine, it is now possible to directly model these types of surface phenomenon and metallic nanostructures without a supercomputer. Fulfilling a High Operating Temperature Propulsion Components (HOTPC) project level-I milestone, a graphical user interface was created for a suite of quantum approximate atomistic materials modeling Fortran programs developed at Glenn. The resulting "Alloy Design Workbench-Surface Modeling Package" (ADW-SMP) is the combination of proven quantum approximate Bozzolo-Ferrante-Smith (BFS) algorithms (refs. 1 and 2) with a productivity-enhancing graphical front end. Written in the portable, platform independent Java programming language, the graphical user interface calls on extensively tested Fortran programs running in the background for the detailed computational tasks. Designed to run on desktop computers, the package has been deployed on PC, Mac, and SGI computer systems. The graphical user interface integrates two modes of computational materials exploration. One mode uses Monte Carlo simulations to determine lowest energy equilibrium configurations. The second approach is an interactive "what if" comparison of atomic configuration energies, designed to provide real-time insight into the underlying drivers of alloying processes.

Abel, Phillip B.↗

Computational Diagnostic Techniques for Electromagnetic Scattering: Analytical Imaging, Near Fields, and Surface Currents

This paper presents three techniques and the graphics implementations which can be used as diagnostic aides in the design and understanding of scattering structures: Imaging, near fields, and surface current displays. The imaging analysis is a new bistatic k space approach which has potential for much greater information than standard experimental approaches. The near field and current analysis are implementations of standard theory while the diagnostic graphics displays are implementations exploiting recent computer engineering work station graphics libraries.

Hom, Kam W.↗

IMP: Interactive mass properties program. Volume 1: Program description

A method of computing a weights and center of gravity analysis of a flight vehicle using interactive graphical capabilities of the Adage 340 computer is described. The equations used to calculate area, volume, and mass properties are based on elemental surface characteristics. The input/output methods employ the graphic support of the Adage computer. Several interactive program options are available for analyzing the mass properties of a vehicle. These options are explained.

Stewart, W. A.↗

Accelerating population balance-Monte Carlo simulation for coagulation dynamics from the Markov jump model, stochastic algorithm and GPU parallel computing

This paper proposes a comprehensive framework for accelerating population balance-Monte Carlo (PBMC) simulation of particle coagulation dynamics. By combining Markov jump model, weighted majorant kernel and GPU (graphics processing unit) parallel computing, a significant gain in computational efficiency is achieved. The Markov jump model constructs a coagulation-rule matrix of differentially-weighted simulation particles, so as to capture the time evolution of particle size distribution with low statistical noise over the full size range and as far as possible to reduce the number of time loopings. Here three coagulation rules are highlighted and it is found that constructing appropriate coagulation rule provides a route to attain the compromise between accuracy and cost of PBMC methods. Further, in order to avoid double looping over all simulation particles when considering the two-particle events (typically, particle coagulation), the weighted majorant kernel is introduced to estimate the maximum coagulation rates being used for acceptance–rejection processes by single-looping over all particles, and meanwhile the mean time-step of coagulation event is estimated by summing the coagulation kernels of rejected and accepted particle pairs. The computational load of these fast differentially-weighted PBMC simulations (based on the Markov jump model) is reduced greatly to be proportional to the number of simulation particles in a zero-dimensional system (single cell). Finally, for a spatially inhomogeneous multi-dimensional (multi-cell) simulation, the proposed fast PBMC is performed in each cell, and multiple cells are parallel processed by multi-cores on a GPU that can implement the massively threaded data-parallel tasks to obtain remarkable speedup ratio (comparing with CPU computation, the speedup ratio of GPU parallel computing is as high as 200 in a case of 100 cells with 10 000 simulation particles per cell). These accelerating approaches of PBMC are demonstrated in a physically realistic Brownian coagulation case. The computational accuracy is validated with benchmark solution of discrete-sectional method. The simulation results show that the comprehensive approach can attain very favorable improvement in cost without sacrificing computational accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The implementation of the graphics of program EAGLE: A numerical grid generation code on NASA Langley SNS computer system

Program EAGLE (Eglin Arbitrary Geometry Implicit Euler) Numerical Grid Generation System is a composite (multi-block) algebraic or elliptic grid generation system designed to discretize the domain in and/or around any arbitrarily shaped three dimensional regions. This system combines a boundary conforming surface generation scheme and includes plotting routines designed to take full advantage of the DISSPLA Graphics Package (Version 9.0). Program EAGLE is written to compile and execute efficiently on any Cray machine with or without solid state disk (SSD) devices. Also, the code uses namelist inputs which are supported by all Cray machines using the FORTRAN compiler CFT77. The namelist inputs makes it easier for the user to understand the inputs and operation of Program EAGLE. EAGLE's numerical grid generator is constructed in the following form: main program, EGG (executive routine); subroutine SURFAC (surface generation routine); subroutine GRID (grid generation routine); and subroutine GRDPLOT (grid plotting routines). The EAGLE code was modified to use on the NASA-LaRC SNS computer (Cray 2S) system. During the modification a conversion program was developed for the output data of EAGLE's subroutine GRID to permit the data to be graphically displayed by IRIS workstations, using Plot3D. The code of program EAGLE was modified to make operational subroutine GRDPLOT (using DI-3000 Graphics Software Packages) on the NASA-LaRC SNS Computer System. How to implement graphically, the output data of subroutine GRID was determined on any NASA-LaRC graphics terminal that has access to the SNS Computer System DI-300 Graphics Software Packages. A Quick Reference User Guide was developed for the use of program EAGLE on the NASA-LaRC SNS Computer System. One or more application program(s) was illustrated using program EAGLE on the NASA LaRC SNS Computer System, with emphasis on graphics illustrations.

Houston, Johnny L.↗

Large Eddy Simulation of a Supercritical Fuel Jet in Cross Flow using GPU-Acceleration

Fuel is routinely injected into combustion chambers at elevated pressures resulting in formation of supercritical fluid flow. The dispersion of fuel and the combustion processes are significantly different in high pressure conditions as the raditional liquid atomization and evaporation are replaced by diffusion-dominated mixing. Modeling supercritical flow is challenging due to the complex thermodynamics and transport properties of the fluids and in most cases high fidelity numerical simulations are prohibitively expensive. However, numerical investigations must incorporate the flow properties appropriate for the regime for accurate description of the physical processes involved. In this paper, we consider real-fluid properties and simulate the turbulent mixing of fuel jet in cross flow at high ambient pressure. In order to make efficient use of the computer resources and reduce the time to solution, without compromising the resolution requirements, the real-fluid property routines are accelerated via Graphics Processing Unit (GPU) computation. The numerical simulations performed at different resolutions suggest that the high fidelity computations with refined mesh are essential to capture the non-linear dynamics of vortex evolution downstream of the fuel jet and signifies the necessity for an accelerated code. The code is accelerated using the Kokkos C++ library and is demonstrated to be about 2 times faster than the non-accelerated version. Considerable computational cost associated with the thermodynamics routines is reduced in the accelerated version. The results also suggest the resource requirements to compute turbulent closures and the need for acceleration of the routines associated with the sub-grid scale closures.

Gottiparthi, Kalyana C.↗

A Study of Complex Deep Learning Networks on High Performance, Neuromorphic, and Quantum Computers

Current Deep Learning models use highly optimized convolutional neural networks (CNN) trained on large graphical processing units (GPU)-based computers with a fairly simple layered network topology, i.e., highly connected layers, without intra-layer connections. Complex topologies have been proposed, but are intractable to train on current systems. Building the topologies of the deep learning network requires hand tuning, and implementing the network in hardware is expensive in both cost and power.In this paper, we evaluate deep learning models using three different computing architectures to address these problems: quantum computing to train complex topologies, high performance computing (HPC) to automatically determine network topology, and neuromorphic computing for a low-power hardware implementation. Due to input size limitations of current quantum computers we use the MNIST dataset for our evaluation.The results show the possibility of using the three architectures in tandem to explore complex deep learning networks that are untrainable using a von Neumann architecture. We show that a quantum computer can find high quality values of intra-layer connections and weights, while yielding a tractable time result as the complexity of the network increases; a high performance computer can find optimal layer-based topologies; and a neuromorphic computer can represent the complex topology and weights derived from the other architectures in low power memristive hardware.This represents a new capability that is not feasible with current von Neumann architecture. It potentially enables the ability to solve very complicated problems unsolvable with current computing technologies.

97 MATHEMATICS AND COMPUTING↗

A Study of Complex Deep Learning Networks on High Performance, Neuromorphic, and Quantum Computers

Current Deep Learning models use highly optimized convolutional neural networks (CNN) trained on large graphical processing units (GPU)-based computers with a fairly simple layered network topology, i.e., highly connected layers, without intra-layer connections. Complex topologies have been proposed, but are intractable to train on current systems. Building the topologies of the deep learning network requires hand tuning, and implementing the network in hardware is expensive in both cost and power. In this paper, we evaluate deep learning models using three different computing architectures to address these problems: quantum computing to train complex topologies, high performance computing (HPC) to automatically determine network topology, and neuromorphic computing for a low-power hardware implementation. Due to input size limitations of current quantum computers we use the MNIST dataset for our evaluation. The results show the possibility of using the three architectures in tandem to explore complex deep learning networks that are untrainable using a von Neumann architecture. We show that a quantum computer can find high quality values of intra-layer connections and weights, while yielding a tractable time result as the complexity of the network increases; a high performance computer can find optimal layer-based topologies; and a neuromorphic computer can represent the complex topology and weights derived from the other architectures in low power memristive hardware. This represents a new capability that is not feasible with current von Neumann architecture. It potentially enables the ability to solve very complicated problems unsolvable with current computing technologies.

Potok, Thomas E↗

Mesoscopic Modeling and Rapid Simulation of Incremental Changes in Epidemic Scenarios on GPUs

In simulation-based studies and analyses of epidemics, a major challenge lies in resolving the conflict between fidelity of models and the speed of their simulation. Another related challenge arises in dealing with the large number of what–if scenarios that need to be explored. Here, we describe new computational methods that together provide an approach to dealing with both challenges. A mesoscopic modeling approach is described that strikes a middle ground between macroscopic models based on coupled differential equations and microscopic models built on fine-grained behaviors at the individual entity level. The mesoscopic approach offers the ability to incorporate complex compositions of multiple layers of dynamics even while retaining the potential for aggregate behaviors at varying levels. It also is an excellent match to the accelerator-based architectures of modern computing platforms in which graphical processing units (GPUs) can be exploited for fast simulation via the parallel execution mode of single instruction multiple thread (SIMT). The challenge of simulating a large number of scenarios is addressed via a method of sharing model state and computation across a tree of what–if scenarios that are localized, incremental changes to a large base simulation. A combination of the mesoscopic modeling approach and the incremental what–if scenario tree evaluation has been implemented in the software on modern GPUs. Synthetic simulation scenarios are presented to demonstrate the computational characteristics of our approach. Results from the experiments with large population data, including USA, UK, and India, illustrate the modeling methodology and computational performance on thousands of synthetically generated what–if scenarios. Execution of our implementation scaled to 8192 GPUs of supercomputing platforms demonstrates the ability to rapidly evaluate what–if scenarios several orders of magnitude faster than the conventional methods.

97 MATHEMATICS AND COMPUTING↗

Robot Display and Control program (ROBDAC)

The ROBOT display and control program (ROBDAC) to be written for the Adage 340 computer is designed to communicate with the Univac 1110 computer with the aid of the graphic support communications programs. ROBDAC and the Adage computer act as an input/output device for the program ROBOT, which is run on the Univac 1110. ROBDAC controls the input to ROBOT by displaying images of the inputs, providing an editing technique for these inputs, shipping the input data to the Univac, and then commanding ROBOT to start execution. The output is done in a similar manner after receiving the output data from the Univac. This data will be displayed numerically as well as by plotted graphs. All of the input/output will be done interactively with practically 'finger tip' control. A general program flow diagram of ROBDAC is included.

Stewart, W. A.↗

GPU Accelerated Prognostics

Prognostic methods enable operators and maintainers to predict the future performance for critical systems. However, these methods can be computationally expensive and may need to be performed each time new information about the system becomes available. In light of these computational requirements, we have investigated the application of graphics processing units (GPUs) as a computational platform for real-time prognostics. Recent advances in GPU technology have reduced cost and increased the computational capability of these highly parallel processing units, making them more attractive for the deployment of prognostic software. We present a survey of model-based prognostic algorithms with considerations for leveraging the parallel architecture of the GPU and a case study of GPU-accelerated battery prognostics with computational performance results.

Prognostics↗

GPU-based simulation of the two-dimensional unstable structure of gaseous oblique detonations

In this paper, the two-dimensional structure of unstable oblique detonations induced by the wedge from a supersonic combustible gas flow is simulated using the reactive Euler equations with a one-step Arrhenius chemistry model. A wide range of activation energy of the combustible mixture is considered. Computations are performed on the Graphical Processing Unit (GPU) to reduce the simulation runtimes. A large computational domain covered by a uniform mesh with high grid resolution is used to properly capture the development of instabilities and the formation of different transverse wave structures. After the initiation point, where the oblique shock transits into a detonation, an instability begins to manifest and in all cases, the left-running transverse waves first appear, followed by the subsequent emergence of right-running transverse waves forming the dual-head triple point structure. This study shows that for low activation energies, a long computational length must be carefully considered to reveal the unstable surface due to the slow growth rate of the instability. For high activation energies, the flow behind the unstable oblique detonation features the formation of unburnt gas pockets and strong vortex-pressure wave interaction resulting in a chaotic-like vortical structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

STARS: A general-purpose finite element computer program for analysis of engineering structures

STARS (Structural Analysis Routines) is primarily an interactive, graphics-oriented, finite-element computer program for analyzing the static, stability, free vibration, and dynamic responses of damped and undamped structures, including rotating systems. The element library consists of one-dimensional (1-D) line elements, two-dimensional (2-D) triangular and quadrilateral shell elements, and three-dimensional (3-D) tetrahedral and hexahedral solid elements. These elements enable the solution of structural problems that include truss, beam, space frame, plane, plate, shell, and solid structures, or any combination thereof. Zero, finite, and interdependent deflection boundary conditions can be implemented by the program. The associated dynamic response analysis capability provides for initial deformation and velocity inputs, whereas the transient excitation may be either forces or accelerations. An effective in-core or out-of-core solution strategy is automatically employed by the program, depending on the size of the problem. Data input may be at random within a data set, and the program offers certain automatic data-generation features. Input data are formatted as an optimal combination of free and fixed formats. Interactive graphics capabilities enable convenient display of nodal deformations, mode shapes, and element stresses.

Gupta, K. K.↗

EVA/manned systems

Viewgraphs on extravehicular activity/manned systems are presented. Objectives of crewstation design include (1) development of human-computer interface technology and graphical presentations, including multi-dimensional visual and aural displays; (2) provide a technology base for autonomous vision and other perceptual systems, virtual workstation technology, and computational vision systems; and (3) develop databases and models of human strength, motion, and body positions in microgravity environments. Extravehicular technology is addressed.

Jenkins, James P.↗

Predicted vibration responses of Apollo structure and effects of pressure correlation lengths on response

Progressive wove and reverberant acoustic fields ore often used in the laboratory for qualification testing of aerospace structures that are subjected to unsteady aerodynamic environments which occur during flight through the atmosphere. Acoustic simulation of these in-flight environments i s achieved i f the structural vibration response to acoustic excitation is equivalent, over the frequency range of interest, to the response caused by aerodynamic turbulence. To determine the degree of response simulation that could be expected for such qualification testing of the Apollo Spacecraft, responses of a segment of the Apollo structure were estimated theoretically for both acoustic and flight environments. The vibration analysis presented in this report was performed for the Spacecraft Lunar Module Adaptor; and the three environments treated are turbulence at Mach 1 and Mach 2, a reverberant acoustic field, and a modified acoustic progressive wave field. The progressive wave field is assumed to be generated within a specially designed shroud of sixteen axially oriented ducts which is constructed around and contoured to the external skin of the vehicle. Each duct is open on the side adjacent to the skin and is driven acoustically at one end by an independent noise source, thus allowing for different pressure correlation patterns around the circumference of the vehicle. The purpose of the analysis was to determine the effects of surface pressure correlation lengths on the response characteristics of representative Apollo structure, and to select on this basis, the optimum acoustic environment for structural vibration qualification tests of the Apollo in the Spacecraft Acoustic Laboratory at MSC. The method of analysis employed consists of determining, independently, the space average, mean-square acceleration spectral density for each of the significant modes of vibration of an equivalent cylindrical shell, and then expressing the total structural response as a linear summation of these modal responses. Approximately 570 classical modes of a pinned-end cylinder were used in order to find the spectral response characteristics of the structure throughout the frequency range of 10 - 1,000 Hz. The effective forcing functions or joint acceptances of the various modes of vibration were computed and are presented graphically for each of the pressure excitation environments. A brief development of the response equations used is presented in the report. Space average acceleration response spectra were computed for twelve different cases including aerodynamic turbulence at Mach 1 and Mach 2, reverberant acoustic field, five duct correlation patterns, two cases of an axially damped progressive wave field, and two cases in which the structural damping constant was varied from the expected value. The response spectra were computed and plotted by a high speed digital computer far a flat excitation pressure spectrum; and the response spectra were numerically integrated to give one-third octave bond average responses. A complete set of digital computer programs has been developed for both flat plate and cylindrical shell structures for the above environments, and each requires run times of three minutes or less. Thus, the work presented herein, along with the results, constitutes a practical advancement to the existing state-of-the-art of vibration prediction.

R W White↗