Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel efficiency”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Transparent Conductive Ink for Additive Manufacturing

NASA analyzes, tests, packages, and fabricates electrical, electronic, and electromechanical (EEE) parts. Nanotechnology is listed in NASA's Technology Roadmap as a key area to invest for further development.1 This research project focused on using nanotechnology to improve electroluminescent lighting in terms of additive manufacturing and to increase energy efficiency. Specifically, this project's goal was to produce a conductive but transparent printable ink that can be sprayed on any surface for use as one of the electrodes in electroluminescent device design. This innovative work is known as thick film dielectric electroluminescent (TDEL) technology. TDEL devices are used for "backlighting, illumination, and identification due to their tunable color output, scalability, and efficiency" (I.K. Small, T.D. Rolin, and A.D. Shields, "3D Printed Electroluminescent Light Panels," NASA Fiscal Year 2017 Center Innovation Fund Proposal, unpublished data, 2017). These devices use a 'front-to-back' printing method, where the substrate is the transparent layer, and the dielectric and phosphor are layered on top. This project is a first step in the process of creating a 3D printable 'back-to-front' electroluminescent device. Back-to-front 3D-printed devices are beneficial because they can be printed onto different substrates and embedded in different surfaces, and the substrate is not required to be transparent, all because the light is emitted from the top surface through the transparent conductor. Advances in this area will help further development of printing TDEL devices on an array of different surfaces. Figure 1 demonstrates the layering of the two electrodes that are aligned in a parallel plate capacitor structure (I.K. Small, T.D. Rolin, and A.D. Shields, "3D Printed Electroluminescent Light Panels," NASA Fiscal Year 2017 Center Innovation Fund Proposal, unpublished data, 2017). Voltage is applied across the device, and the subsequent electron excitation results in light emission at the top layer.

Patlan, X. J.↗
NASA NTRS2023

Easy, Scalable Subsetting of GEDI Point Clouds

The GEDI Subsetter, a Python tool developed for NASA’s Multi-mission Algorithm and Analysis Platform (MAAP), optimizes the accessibility and visualization of GEDI point clouds by enabling users to efficiently subset data in a convenient, scalable manner. Complex science data often requires users to learn new software skills and handle many large files. Handling and cleaning large data sets is tedious and error-prone. These challenges significantly impede analysis. One of the goals of NASA's MAAP is to provide a platform that lowers the barrier to conducting research and analysis at scale. When a group of MAAP users wanted to conduct above-ground biomass estimation using GEDI data, we found that their existing workflow for leveraging GEDI data suffered from the barriers mentioned above. Furthermore, their workflow did not scale easily beyond a small number of granules. We found that existing tools related to GEDI data retrieval and subsetting were too limiting, so the GEDI Subsetter was written to support MAAP users’ needs. Being able to run many subsetting jobs simultaneously in the MAAP, and parallelizing the code itself, has led to significant speed improvements in obtaining relevant data, reducing subsetting time from hours to minutes. MAAP users can now more quickly and easily obtain only the data relevant to their research, by choosing which GEDI collection they want to work with (L1A, L2A, L2B, or L4A), and how they want to subset it, by specifying an area of interest, a temporal range, and relevant attributes. This has significantly reduced the feedback loop for users, allowing them to much more quickly subset GEDI data and begin their analysis. Although the GEDI Subsetter originally targeted users of the MAAP, it is generalized such that it can also be used outside of the MAAP and includes a command-line interface for convenience. Furthermore, with minor modifications, it should be possible to use it with non-GEDI data as the general pattern should be applicable to other sparse/track-based sensors.

Charles Daniels↗
NASA NTRS1993

Preliminary results: Root cause investigation of orbital anomalies and failures in NASA standard 50 ampere-hour nickel-cadmium batteries

Two lots of NASA standard 50 A.H. Ni-Cd battery cells, manufactured by Gates Aerospace Batteries and built into batteries by McDonnell Douglas, have experienced significant performance problems. The two lots were used on the Compton Gamma Ray Observatory and the Upper Atmosphere Research Satellite. Both of these satellites are Low Earth Orbital (LEO) satellites containing batteries on a parallel bus charged to NASA standard V/T curves using a NASA standard power regulator. The following preliminary conclusions were reached: (1) several plate and cell parameters have migrated within their spec limits over the years (in some cases, from one extreme to the other); (2) several parametric relationships, not generally monitored and therefore not under specification control, have also migrated over the years; (3) many of these changes appear to have taken place as a natural consequence of changes in GE/GAB materials and processes; (4) several of these factors may be 'conspiring' to aggravate known cell failure mechanisms (factors such as heavier plate, less teflon and/or less-uniform teflon, and less electrolyte) but all are still in spec (where specs exist); (5) the weight of the evidence collected to characterize the anomalies and to characterize the negative electrode itself, strongly suggests that alterations to the structure, composition, uniformity, and efficiency of the negative electrode are at the heart of the battery performance problems currently being experienced; and (6) further investigation at all levels (plate, cell, battery, and system) continues to be warranted.

Toft, Mark R.↗
NASA NTRS1993

Supersonic flow and shock formation in turbine tip gaps

Shock formation due to overexpansion of supersonic flow at the inlet to the tip clearance gap of a turbomachine has been studied. As the flow enters the tip gap, it accelerates around the blade pressure-side corner creating a region of minimum static pressure. The 'free streamline' separates from the wall at the corner; and, for Mach numbers greater than about 1.3, it curves back to intersect the blade tip. At this point, the freestream flow is abruptly turned parallel to the surface, giving rise to an oblique shock. The results are consistent with compressible sharp-edged orifice flow calculations found in the literature and with the theory of oblique shock wave formation in supersonic flow over a wedge. For freestream Mach numbers of 1.4 to 1.8, wave angles are 43 to 54 deg, and turning angles are 9 to 20 deg; as the Mach number increases, the angle of turn also increases. It appears that in a turbine, after separating from the inlet corner, the flow reattaches on the blade tip and an oblique shock is formed at 0.4-1.4 tip gap heights into the clearance gap. The resulting shock-boundary layer interaction may contribute to further enhancement of already high heat transfer to the blade tip in this region. This in turn could lead to higher blade temperatures and adversely affect blade life and turbine efficiency.

Moore, John↗
NASA NTRS1996

Development and Demonstration of an Ada Test Generation System

In this project we have built a prototype system that performs Feasible Path Analysis on Ada programs: given a description of a set of control flow paths through a procedure, and a predicate at a program point feasible path analysis determines if there is input data which causes execution to flow down some path in the collection reaching the point so that tile predicate is true. Feasible path analysis can be applied to program testing, program slicing, array bounds checking, and other forms of anomaly checking. FPA is central to most applications of program analysis. But, because this problem is formally unsolvable, syntactic-based approximations are used in its place. For example, in dead-code analysis the problem is to determine if there are any input values which cause execution to reach a specified program point. Instead an approximation to this problem is computed: determine whether there is a control flow path from the start of the program to the point. This syntactic approximation is efficiently computable and conservative: if there is no such path the program point is clearly unreachable, but if there is such a path, the analysis is inconclusive, and the code is assumed to be live. Such conservative analysis too often yields unsatisfactory results because the approximation is too weak. As another example, consider data flow analysis. A du-pair is a pair of program points such that the first point is a definition of a variable and the second point a use and for which there exists a definition-free path from the definition to the use. The sharper, semantic definition of a du-pair requires that there be a feasible definition-free path from the definition to the use. A compiler using du-pairs for detecting dead variables may miss optimizations by not considering feasibility. Similarly, a program analyzer computing program slices to merge parallel versions may report conflicts where none exist. In the context of software testing, feasibility analysis plays an important role in identifying testing requirements which are infeasible. This is especially true for data flow testing and modified condition/decision coverage. Our system uses in an essential way symbolic analysis and theorem proving technology, and we believe this work represents one of the few successful uses of a theorem prover working in a completely automatic fashion to solve a problem of practical interest. We believe this work anticipates an important trend away from purely syntactic-based methods for program analysis to semantic methods based on symbolic processing and inference technology. Other results demonstrating the practical use of automatic inference is being reported in hardware verification, although there are significant differences between the hardware work and ours. However, what is common and important is that general purpose theorem provers are being integrated with more special-purpose decision procedures to solve problems in analysis and verification. We are pursuina commercial opportunities for this work, and will use and extend the work in other projects we are engaged in. Ultimately we would like to rework the system to analyze C, C++, or Java as a key step toward commercialization.

Source record↗
NASA NTRS2011

Volumetric 3D Display System with Static Screen

Current display technology has relied on flat, 2D screens that cannot truly convey the third dimension of visual information: depth. In contrast to conventional visualization that is primarily based on 2D flat screens, the volumetric 3D display possesses a true 3D display volume, and places physically each 3D voxel in displayed 3D images at the true 3D (x,y,z) spatial position. Each voxel, analogous to a pixel in a 2D image, emits light from that position to form a real 3D image in the eyes of the viewers. Such true volumetric 3D display technology provides both physiological (accommodation, convergence, binocular disparity, and motion parallax) and psychological (image size, linear perspective, shading, brightness, etc.) depth cues to human visual systems to help in the perception of 3D objects. In a volumetric 3D display, viewers can watch the displayed 3D images from a completely 360 view without using any special eyewear. The volumetric 3D display techniques may lead to a quantum leap in information display technology and can dramatically change the ways humans interact with computers, which can lead to significant improvements in the efficiency of learning and knowledge management processes. Within a block of glass, a large amount of tiny dots of voxels are created by using a recently available machining technique called laser subsurface engraving (LSE). The LSE is able to produce tiny physical crack points (as small as 0.05 mm in diameter) at any (x,y,z) location within the cube of transparent material. The crack dots, when illuminated by a light source, scatter the light around and form visible voxels within the 3D volume. The locations of these tiny voxels are strategically determined such that each can be illuminated by a light ray from a high-resolution digital mirror device (DMD) light engine. The distribution of these voxels occupies the full display volume within the static 3D glass screen. This design eliminates any moving screen seen in previous approaches, so there is no image jitter, and has an inherent parallel mechanism for 3D voxel addressing. High spatial resolution is possible with a full color display being easy to implement. The system is low-cost and low-maintenance.

Geng, Jason↗
NASA NTRS1996

Multigrid with Overlapping Patches

Solving boundary value problems with optimal efficiency requires adaptivity and multilevel techniques. Previously, an implementation of the AFACx algorithm is presented that is based on rectangular Cartesian grids. This implementation does not allow for the over]ap of grids that lie on the same level of refinement. We investigate the case in which these grids overlap. A standard technique for overlapping grids is the Schwarz algorithm. Some ways of using the Schwarz algorithm in a standard multigrid scheme are presented. Also, a problem that arises in some situations with non-aligned, overlapping grids is described. This situation comes up in a natural way when the Schwarz algorithm is used as a relaxation scheme within a multilevel algorithm. We identify the reason for the bad convergence and show that by more sophisticated interpolation the difficulties can be overcome. Then we present a multiplicative Schwarz algorithm for a large number of grids that has a high potential for parallelization. Finally we give some numerical results for the FACx algorithm with overlapping grids on each refinement level. The implementation of the described codes uses C++ and the array class libraries A++ and P++. Using the A++/P++ programming environment, it was possible to move from a serial code to a parallel code within a few days.

Berndt, Markus↗
NASA NTRS1998

High-Density Plasma Reactors: Simulations for Design

The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.

Hash, David B.↗
NASA NTRS2003

Mechanical Properties of Triaxial Braided Carbon/Epoxy Composites

In an on-going effort to increase the safety and efficiency of turbine engines, the National Aeronautics and Space Administration is exploring lightweight alternatives to the metal containment structures that currently encase commercial jet engines. Epoxy reinforced with braided carbon fibers is a candidate structural material which may be suitable for an engine case. This paper reports flat-coupon mechanical-property experiments performed to compliment previously reported subcomponent impact testing and analytical simulation of containment structures. Triaxial-braid T700/5208 epoxy and triaxial-braid T700/M36 toughened epoxy composites were evaluated. Also, two triaxial-braid architectures (0 +/- 60 deg., 0 +/- 45 deg.) with the M36 resin were evaluated through tension, compression, and shear testing. Tensile behavior was compared between standard straight-sided specimens (ASTM D3039) and bowtie specimens. Both double-notch shear (ASTM D3846) and Iosepescu (ASTM D5379) tests were performed as well. The M36/0 +/- 45 deg. configuration yield the best response when measurements were made parallel to the axial tows. Conversely, the M36/0 +/- 60 deg. configuration was best when measurements were made perpendicular to the axial tows. The results were used to identify critical properties and to augment the analysis of impact experiments.

Bowman, C. L.↗
NASA NTRSDate not supplied

Radiative Heat Transfer Capability Implemented in OpenNCC for Conjugate Heat Transfer Applications

Thermal efficiency of gas turbine engine increases as the temperature and pressure at the combustor increases. Consequently, the materials used inside a combustor must survive an increasingly challenging environment. For this reason, accurate assessment of heat transfer is crucial for combustor design. While all three modes of heat transfer are present inside a combustor, the focus of this paper is the thermal radiation. Radiative heat transfer in a gas turbine combustors are particularly interesting from three reasons. Firstly, the radiative heat loss from the combustion region may affect the emission performance. Secondly, the cooling air will protect the liner from convection but not necessary from radiation. Finally, it is less frequently incorporated in CFD analysis than other forms of heat transfer. In this work, radiative heat transfer using discrete ordinate method has been incorporated in OpenNCC (a publicly releasable version of the National Combustion Code) developed at NASA Glenn Research Center. Aside from massively parallel computation capability using MPI and the ability to utilize unstructured mesh, the current implementation includes two types of spectral models, namely, the weighted some of gray gas model and the full spectrum correlated k-distribution model. After presenting the theory and the strategy of implementation, results of validation cases for gray gas and spectral models will be presented. While the implementation of the radiation solver is intended for gas turbine application, the radiation solver can run independently from the convection/combustion solver and the same theory can be applied to other application.

OpenNCC↗
NASA NTRSDate not supplied

Radiative Heat Transfer Capability Implemented in OpenNCC for Conjugate Heat Transfer Applications

Thermal efficiency of gas turbine engine increases as the temperature and pressure at the combustor increases. Consequently, the materials used inside a combustor must survive an increasingly challenging environment. For this reason, accurate assessment of heat transfer is crucial for combustor design. While all three modes of heat transfer are present inside a combustor, the focus of this paper is the thermal radiation. Radiative heat transfer in a gas turbine combustors are particularly interesting from three reasons. Firstly, the radiative heat loss from the combustion region may affect the emission performance. Secondly, the cooling air will protect the liner from convection but not necessary from radiation. Finally, it is less frequently incorporated in CFD analysis than other forms of heat transfer. In this work, radiative heat transfer using discrete ordinate method has been incorporated in OpenNCC (a publicly releasable version of the National Combustion Code) developed at NASA Glenn Research Center. Aside from massively parallel computation capability using MPI and the ability to utilize unstructured mesh, the current implementation includes two types of spectral models, namely, the weighted some of gray gas model and the full spectrum correlated k-distribution model. After presenting the theory and the strategy of implementation, results of validation cases for gray gas and spectral models will be presented. While the implementation of the radiation solver is intended for gas turbine application, the radiation solver can run independently from the convection/combustion solver and the same theory can be applied to other application.

OpenNCC↗
NASA NTRSDate not supplied

An Approach to Solving Enclosure Radiation Problems in A Multi-Physics Context

Thermal protection system analysis of complex features or damage sites can sometimes require modeling of high temperature enclosures. Implementing efficient and accurate view-factor algorithms required to model such problems is complex. The current work leverages the Non-equilibrium Radiation (NERO) software, which solves the radiation transport equation in a finite-volume scheme, to alleviating challenges often faced with view-factor calculations. By assuming heat transfer occurs only between grey bodies and that the medium is non-participating, computational cost of the method is significantly reduced. The enclosure physics are modeled through emitting and reflecting boundary conditions in NERO. The emitted radiative flux is dependent on the wall temperature which is a solution to the material response, obtained from Icarus, in this context. The Ares framework manages the time-advancement and exchange of the necessary data between the solvers. The surface energy balance is modified to account for the enclosure terms within the material response boundary condition. The methodology was verified against analytical solutions including radiating parallel plates, a hollow cylinder (shown in Fig. 1), and a hemisphere. Application of the methodology to inform the design of components of the Dragonfly system will be shown.

Ablation↗
NASA NTRS1990

On recursive least-squares filtering algorithms and implementations

In many real-time signal processing applications, fast and numerically stable algorithms for solving least-squares problems are necessary and important. In particular, under non-stationary conditions, these algorithms must be able to adapt themselves to reflect the changes in the system and take appropriate adjustments to achieve optimum performances. Among existing algorithms, the QR-decomposition (QRD)-based recursive least-squares (RLS) methods have been shown to be useful and effective for adaptive signal processing. In order to increase the speed of processing and achieve high throughput rate, many algorithms are being vectorized and/or pipelined to facilitate high degrees of parallelism. A time-recursive formulation of RLS filtering employing block QRD will be considered first. Several methods, including a new non-continuous windowing scheme based on selectively rejecting contaminated data, were investigated for adaptive processing. Based on systolic triarrays, many other forms of systolic arrays are shown to be capable of implementing different algorithms. Various updating and downdating systolic algorithms and architectures for RLS filtering are examined and compared in details, which include Householder reflector, Gram-Schmidt procedure, and Givens rotation. A unified approach encompassing existing square-root-free algorithms is also proposed. For the sinusoidal spectrum estimation problem, a judicious method of separating the noise from the signal is of great interest. Various truncated QR methods are proposed for this purpose and compared to the truncated SVD method. Computer simulations provided for detailed comparisons show the effectiveness of these methods. This thesis deals with fundamental issues of numerical stability, computational efficiency, adaptivity, and VLSI implementation for the RLS filtering problems. In all, various new and modified algorithms and architectures are proposed and analyzed; the significance of any of the new method depends crucially on specific application.

Hsieh, Shih-Fu↗
NASA NTRS1998

An Atmospheric General Circulation Model with Chemistry for the CRAY T3E: Design, Performance Optimization and Coupling to an Ocean Model

The design, implementation and performance optimization on the CRAY T3E of an atmospheric general circulation model (AGCM) which includes the transport of, and chemical reactions among, an arbitrary number of constituents is reviewed. The parallel implementation is based on a two-dimensional (longitude and latitude) data domain decomposition. Initial optimization efforts centered on minimizing the impact of substantial static and weakly-dynamic load imbalances among processors through load redistribution schemes. Recent optimization efforts have centered on single-node optimization. Strategies employed include loop unrolling, both manually and through the compiler, the use of an optimized assembler-code library for special function calls, and restructuring of parts of the code to improve data locality. Data exchanges and synchronizations involved in coupling different data-distributed models can account for a significant fraction of the running time. Therefore, the required scattering and gathering of data must be optimized. In systems such as the T3E, there is much more aggregate bandwidth in the total system than in any particular processor. This suggests a distributed design. The design and implementation of a such distributed 'Data Broker' as a means to efficiently couple the components of our climate system model is described.

Farrara, John D.↗
NASA NTRS2000

Web Exploration Tools for a Fast Federated Optical Survey Database

We implemented several new web-based tools to improve the efficiency and versatility of access to the APS Catalog of the POSS I (Palomar Observatory-National Geographic Sky Survey) and its associated image database. The most important addition was a federated database system to link the APS Catalog and image database into one Internet-accessible database. With the FDBS, the queries and transactions on the integrated database are performed as if it were a single database. We installed Myriad the FDBS developed by Professor Jaideep Srivastava and members of his group in the University of Minnesota Computer Science Department. It is the first system to provide schema integration, query processing and optimization, and transaction management capabilities in a single framework. The attached figure illustrates the Myriad architecture. The FDBS permits horizontal access to the data, not just vertical. For example, for the APS, queries can be made not only by sky position, but also by any parameter present in either of the databases. APS users will be able to produce an image of all the blue galaxies and stellar sources for comparison with x-ray source error ellipses from AXAF (X Ray Astrophysics Facility) (Chandra) for example. The FDBS is now available as a beta release with the appropriate query forms at our web site. While much of our time was occupied with adapting Myriad to the APS environment, we also made major changes in Star Base, our DBMS for the Catalog, at the web interface to improve its efficiency for issuing and processing queries. Star Base is now three times faster for large queries. Improvements were also made at the web end of the image database for faster access; although work still needs to be done to the image database itself for more efficient return with the FDBS. During the past few years, we made several improvements to the database pipeline that creates the individual plate databases queries by StarBase. The changes include improved positions especially for galaxies, using a new median centroider and integrated magnitudes for galaxies with an improved density-to-intensity calibration with a "sky" background subtraction. In the original version of StarBase the object classification fainter than 19.5-20.0 mag., was an extrapolation of the networks trained on brighter objects. We have used a new catalog of galaxies at the NGP to train a neural network on objects fainter than 20th mag. This improved classification is used in the new version of StarBase. We have also added a FITS table option for the returned data from queries on the object catalog. The APS image database includes images in both colors so we have added a tool for querying the image database in both colors simultaneously. The images can be displayed in parallel or blinked for comparison.

Humphreys, Roberta M.↗
NASA NTRS2003

Mechanical Properties of Triaxial Braided Carbon/Epoxy Composites

In an on-going effort to increase the safety and efficiency of turbine engines, the National Aeronautics and Space Administration is exploring lightweight alternatives to the metal containment structures that currently encase commercial jet engines. Epoxy reinforced with braided carbon fibers is a candidate structural material which may be suitable for an engine case. This paper reports flat-coupon mechanical-property experiments performed to compliment previously reported subcomponent impact testing and analytical simulation of containment structures. Triaxial-braid T700/5208 epoxy and triaxial-braid T700h436 toughened epoxy composites were evaluated. Also, two triaxial-braid architectures (0 degrees plus or minus 60 degrees, and 0 degrees plus or minus 45 degrees) with the M36 resin were evaluated through tension, compression, and shear testing. Tensile behavior was compared between standard straight-sided specimens (ASTM D3039) and bow-tie specimens. Both double-notch shear (ASTM D3846) and Iosepescu (ASTM D5379) tests were performed as well. The M36/O degrees plus or minus 45 degrees configuration yield the best response when measurements were made parallel to the axial tows. Conversely, the M36/0 degrees plus or minus 60 degrees configuration was best when measurements were made perpendicular to the axial tows. The results were used to identify critical properties and to augment the analysis of impact experiments.

Bowman, C. L.↗
NASA NTRS1994

Implementation and Assessment of Advanced Analog Vector-Matrix Processor

This paper discusses the design and implementation of an analog optical vecto-rmatrix coprocessor with a throughput of 128 Mops for a personal computer. Vector matrix calculations are inherently parallel, providing a promising domain for the use of optical calculators. However, to date, digital optical systems have proven too cumbersome to replace electronics, and analog processors have not demonstrated sufficient accuracy in large scale systems. The goal of the work described in this paper is to demonstrate a viable optical coprocessor for linear operations. The analog optical processor presented has been integrated with a personal computer to provide full functionality and is the first demonstration of an optical linear algebra processor with a throughput greater than 100 Mops. The optical vector matrix processor consists of a laser diode source, an acoustooptical modulator array to input the vector information, a liquid crystal spatial light modulator to input the matrix information, an avalanche photodiode array to read out the result vector of the vector matrix multiplication, as well as transport optics and the electronics necessary to drive the optical modulators and interface to the computer. The intent of this research is to provide a low cost, highly energy efficient coprocessor for linear operations. Measurements of the analog accuracy of the processor performing 128 Mops are presented along with an assessment of the implications for future systems. A range of noise sources, including cross-talk, source amplitude fluctuations, shot noise at the detector, and non-linearities of the optoelectronic components are measured and compared to determine the most significant source of error. The possibilities for reducing these sources of error are discussed. Also, the total error is compared with that expected from a statistical analysis of the individual components and their relation to the vector-matrix operation. The sufficiency of the measured accuracy of the processor is compared with that required for a range of typical problems. Calculations resolving alloy concentrations from spectral plume data of rocket engines are implemented on the optical processor, demonstrating its sufficiency for this problem. We also show how this technology can be easily extended to a 100 x 100 10 MHz (200 Cops) processor.

Gary, Charles K.↗
NASA NTRS2013

Time-Accurate Local Time Stepping and High-Order Time CESE Methods for Multi-Dimensional Flows Using Unstructured Meshes

With the wide availability of affordable multiple-core parallel supercomputers, next generation numerical simulations of flow physics are being focused on unsteady computations for problems involving multiple time scales and multiple physics. These simulations require higher solution accuracy than most algorithms and computational fluid dynamics codes currently available. This paper focuses on the developmental effort for high-fidelity multi-dimensional, unstructured-mesh flow solvers using the space-time conservation element, solution element (CESE) framework. Two approaches have been investigated in this research in order to provide high-accuracy, cross-cutting numerical simulations for a variety of flow regimes: 1) time-accurate local time stepping and 2) highorder CESE method. The first approach utilizes consistent numerical formulations in the space-time flux integration to preserve temporal conservation across the cells with different marching time steps. Such approach relieves the stringent time step constraint associated with the smallest time step in the computational domain while preserving temporal accuracy for all the cells. For flows involving multiple scales, both numerical accuracy and efficiency can be significantly enhanced. The second approach extends the current CESE solver to higher-order accuracy. Unlike other existing explicit high-order methods for unstructured meshes, the CESE framework maintains a CFL condition of one for arbitrarily high-order formulations while retaining the same compact stencil as its second-order counterpart. For large-scale unsteady computations, this feature substantially enhances numerical efficiency. Numerical formulations and validations using benchmark problems are discussed in this paper along with realistic examples.

Chang, Chau-Lyan↗