Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel methods”

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 865 records · Page 48

Fully-Coupled Fluid-Structure Interaction Simulations of a Supersonic Parachute

A validated computational fluid-structure interaction method for simulating the complex interaction between the large deformation of very thin, highly deformable structures and compressible flows is extended to consider large-scale problems in supersonic flows using parallel computing. The coupled fluid-structure interaction system is solved in a partitioned, or weakly-coupled, manner. The foundations of the applied fluid-structure interaction method are a higher-order, block-structured Cartesian, sharp immersed boundary method for the compressible Navier-Stokes equations and a computational structural dynamics solver employing a geometrically nonlinear 3-node shell element based on the mixed interpolation of tensorial components formulation. The method is applied to large deformation fluid-structure interaction validation cases before being applied to the inflation of a supersonic parachute in the upper Martian atmosphere where the goal is to demonstrate the capabilities of the solver when considering large-scale problems in supersonic flows.

Boustani, Jonathan↗

Force-free magnetic fields - The magneto-frictional method

The problem under discussion is that of calculating magnetic field configurations in which the Lorentz force j x B is everywhere zero, subject to specified boundary conditions. We choose to represent the magnetic field in terms of Clebsch variables in the form B = grad alpha x grad beta. These variables are constant on any field line so that each field line is labeled by the corresponding values of alpha and beta. When the field is described in this way, the most appropriate choice of boundary conditions is to specify the values of alpha and beta on the bounding surface. We show that such field configurations may be calculated by a magneto-frictional method. We imagine that the field lines move through a stationary medium, and that each element of magnetic field is subject to a frictional force parallel to and opposing the velocity of the field line. This concept leads to an iteration procedure for modifying the variables alpha and beta, that tends asymptotically towards the force-free state. We apply the method first to a simple problem in two rectangular dimensions, and then to a problem of cylindrical symmetry that was previously discussed by Barnes and Sturrock (1972). In one important respect, our new results differ from the earlier results of Barnes and Sturrock, and we conclude that the earlier article was in error.

Yang, W. H.↗

Remote Sensing and Quantization of Analog Sensors

This method enables sensing and quantization of analog strain gauges. By manufacturing a piezoelectric sensor stack in parallel (physical) with a piezoelectric actuator stack, the capacitance of the sensor stack varies in exact proportion to the exertion applied by the actuator stack. This, in turn, varies the output frequency of the local sensor oscillator. The output, F(sub out), is fed to a phase detector, which is driven by a stable reference, F(sub ref). The output of the phase detector is a square waveform, D(sub out), whose duty cycle, t(sub W), varies in exact proportion according to whether F(sub out) is higher or lower than F(sub ref). In this design, should F(sub out) be precisely equal to F(sub ref), then the waveform has an exact 50/50 duty cycle. The waveform, D(sub out), is of generally very low frequency suitable for safe transmission over long distances without corruption. The active portion of the waveform, t(sub W), gates a remotely located counter, which is driven by a stable oscillator (source) of such frequency as to give sufficient digitization of t(sub W) to the resolution required by the application. The advantage to this scheme is that it negates the most-common, present method of sending either very low level signals (viz. direct output from the sensors) across great distances (anything over one-half meter) or the need to transmit widely varying higher frequencies over significant distances thereby eliminating interference [both in terms of beat frequency generation and in-situ EMI (electromagnetic interference)] caused by ineffective shielding. It also results in a significant reduction in shielding mass.

Strauss, Karl F.↗

Two-stream approximations to radiative transfer in planetary atmospheres - A unified description of existing methods and a new improvement

Existing two-stream approximations to radiative transfer theory for particulate media are shown to be represented by identical forms of coupled differential equations if the intensity is replaced by integrals of the intensity over hemispheres. One set of solutions thus suffices for all methods and provides convenient analytical comparisons. The equations also suggest modifications of the standard techniques so as to duplicate exact solutions for thin atmospheres and thus permit accurate determinations of the effects of typical aerosol layers. Numerical results for the plane albedos of plane-parallel atmospheres are given for conventional and modified Eddington approximations, conventional and modified two-point quadrature schemes, the hemispheric-constant method and the delta-function method, all for comparison with accurate discrete-ordinate solutions. A new two-stream approximation is introduced that reduces to the modified Eddington approximation in the limit of isotropic phase functions and to the exact solution in the limit of extreme anisotropic scattering. Comparisons of plane albedos and transmittances show the new method to be generally superior over a wide range of atmospheric conditions (including cloud and aerosol layers), especially in the case of nonconservative scattering.

Meador, W. E.↗

Geometrically Flexible and Efficient Flow Analysis of High Speed Vehicles Via Domain Decomposition, Part 1: Unstructured-Grid Solver for High Speed Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. This paper describes the key numerical methods and techniques that were found to be required to robustly obtain accurate solutions to hypersonic flows on non-hex-dominant unstructured grids. The methods and techniques described include: an augmented stencil, weighted linear least squares, cell-average gradient method, a robust multidimensional cell-average gradient-limiter process that is consistent with the augmented stencil of the cell-average gradient method and a cell-face gradient method that contains a cell skewness sensitive damping term derived using hyperbolic diffusion based concepts. A data-parallel matrix-based symmetric Gauss-Seidel point-implicit scheme, used to solve the governing equations, is described and shown to be more robust and efficient than a matrix-free alternative. In addition, a y+ adaptive turbulent wall boundary condition methodology is presented. This boundary condition methodology is deigned to automatically switch between a solve-to-the-wall and a wall-matching-function boundary condition based on the local y+ of the 1st cell center off the wall. The aforementioned methods and techniques are then applied to a series of hypersonic and supersonic turbulent flat plate unit tests to examine the efficiency, robustness and convergence behavior of the implicit scheme and to determine the ability of the solve-to-the-wall and y+ adaptive turbulent wall boundary conditions to reproduce the turbulent law-of-the-wall. Finally, the thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet flow-path is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the unstructured-grid solver for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗

Computation and parallel implementation for early vision

The problem of early vision is to transform one or more retinal illuminance images-pixel arrays-to image representations built out of such primitive visual features such as edges, regions, disparities, and clusters. These transformed representations form the input to later vision stages that perform higher level vision tasks including matching and recognition. Researchers developed algorithms for: (1) edge finding in the scale space formulation; (2) correlation methods for computing matches between pairs of images; and (3) clustering of data by neural networks. These algorithms are formulated for parallel implementation of SIMD machines, such as the Massively Parallel Processor, a 128 x 128 array processor with 1024 bits of local memory per processor. For some cases, researchers can show speedups of three orders of magnitude over serial implementations.

Gualtieri, J. Anthony↗

High precision computing with charge domain devices and a pseudo-spectral method therefor

The present invention enhances the bit resolution of a CCD/CID MVM processor by storing each bit of each matrix element as a separate CCD charge packet. The bits of each input vector are separately multiplied by each bit of each matrix element in massive parallelism and the resulting products are combined appropriately to synthesize the correct product. In another aspect of the invention, such arrays are employed in a pseudo-spectral method of the invention, in which partial differential equations are solved by expressing each derivative analytically as matrices, and the state function is updated at each computation cycle by multiplying it by the matrices. The matrices are treated as synaptic arrays of a neural network and the state function vector elements are treated as neurons. In a further aspect of the invention, moving target detection is performed by driving the soliton equation with a vector of detector outputs. The neural architecture consists of two synaptic arrays corresponding to the two differential terms of the soliton-equation and an adder connected to the output thereof and to the output of the detector array to drive the soliton equation.

Barhen, Jacob↗

Current distribution on a cylindrical antenna with parallel orientation in a lossy magnetoplasma

The current distribution and impedance of a thin cylindrical antenna with parallel orientation to the static magnetic field of a lossy magnetoplasma is calculated with the method of moments. The electric field produced by an infinitesimal current source is first derived. Results are presented for a wide range of plasma parameters. Reasonable answers are obtained for all cases except for the overdense hyperbolic case. A discussion of the numerical stability is included which not only applies to this problem but other applications of the method of moments.

Klein, C. A.↗

A comparison of five benchmarks

Five benchmark programs were obtained and run on the NASA Lewis CRAY X-MP/24. A comparison was made between the programs codes and between the methods for calculating performance figures. Several multitasking jobs were run to gain experience in how parallel performance is measured.

Huss, Janice E.↗

Digital and optical shape representation and pattern recognition; Proceedings of the Meeting, Orlando, FL, Apr. 4-6, 1988

The present conference discusses topics in pattern-recognition correlator architectures, digital stereo systems, geometric image transformations and their applications, topics in pattern recognition, filter algorithms, object detection and classification, shape representation techniques, and model-based object recognition methods. Attention is given to edge-enhancement preprocessing using liquid crystal TVs, massively-parallel optical data base management, three-dimensional sensing with polar exponential sensor arrays, the optical processing of imaging spectrometer data, hybrid associative memories and metric data models, the representation of shape primitives in neural networks, and the Monte Carlo estimation of moment invariants for pattern recognition.

Juday, Richard D.↗

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict 3-D ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region.

Potapczuk, M. G.↗

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict three-dimensional ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however, this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region.

Potapczuk, M. G.↗

Arc electrode interaction study

The project consisted of two parts: (1) the cathode interaction studies which were a continuation of previous work and had the objective of increasing our understanding of the microscopic phenomena controlling cathode erosion in arc jet thrusters, and (2) the studies of the anode attachment in arc jet thrusters. The cathode interaction studies consisted of (1) a continuation of some modeling work in which the previously derived model for the cathode heating was applied to some specific gases and electrode materials, and (2) experimental work in which various diagnostics was applied to the cathode. The specific diagnostics used were observation of the cathode tip during arcing using a Laser Strobe Video system in conjunction with a tele-microscope, a monochromator with an optical multichannel analyzer for the determination of the cathode temperature distribution, and various ex situ materials analysis methods. The emphasis of our effort was shifted to the cathode materials analysis because a parallel project was in place during the second half of 1993 with a visiting scientist pursuing arc electrode materials studies. As a consequence, the diagnostic investigations of the arc in front of the cathode had to be postponed to the first half of 1994, and we are presently preparing these measurements. The results of last year's study showed some unexpected effects influencing the cathode erosion behavior, such as increased erosion away from the cathode tip, and our understanding of these effects should improve our ability to control cathode erosion. The arc jet anode attachment studies concentrated on diagnostics of the instabilities in subsonic anode attachment arc jet thrusters, and were supplemental measurements to work which was performed by one of the authors who spent the summer as an intern at NASA Lewis Research Center. A summary of the results obtained during the internship are included because they formed an integral part of the study. Two tasks for 1994, the diagnostics of the anode closure phenomenon, and the use of arc jet thrusters for the deposition of c-BN, are being prepared.

Zhou, X.↗

Polarization from scattering in blobs

We describe a Monte Carlo method for calculating polarization due to scattering in various geometries. We apply the method in this paper to the case of scattering in a spherical cloud illuminated by parallel rays. The polarization is largest in small optical depth blobs, where single scattering dominates. The polarization decreases as optical depth increases until a limit is reached where scattering occurs on the surface, and increasing optical depth has little effect on the results. For a given optical depth, decreasing the albedo leads to an increase of the degree of polarization of the scattered light because single scattering increases relative to multiple scattering. As the optical depth of the blob increases, the scattered flux becomes increasingly backward throwing. In the high albedo blobs, the maximum polarization becomes skewed towards the illuminated hemisphere at high optical depths. We discuss the applications to polarization in supergiants and R CrB stars, and scattering in a clumpy interstellar medium.

Code, Arthur D.↗

New Approaches to HSCT Multidisciplinary Design and Optimization

New approaches to MDO have been developed and demonstrated during this project on a particularly challenging aeronautics problem- HSCT Aeroelastic Wing Design. To tackle this problem required the integration of resources and collaboration from three Georgia Tech laboratories: ASDL, SDL, and PPRL, along with close coordination and participation from industry. Its success can also be contributed to the close interaction and involvement of fellows from the NASA Multidisciplinary Analysis and Optimization (MAO) program, which was going on in parallel, and provided additional resources to work the very complex, multidisciplinary problem, along with the methods being developed. The development of the Integrated Design Engineering Simulator (IDES) and its initial demonstration is a necessary first step in transitioning the methods and tools developed to larger industrial sized problems of interest. It also provides a framework for the implementation and demonstration of the methodology. Attachment: Appendix A - List of publications. Appendix B - Year 1 report. Appendix C - Year 2 report. Appendix D - Year 3 report. Appendix E - accompanying CDROM.

Schrage, Daniel P.↗

Heart Fibrillation and Parallel Supercomputers

The Luo and Rudy 3 cardiac cell mathematical model is implemented on the parallel supercomputer CRAY - T3D. The splitting algorithm combined with variable time step and an explicit method of integration provide reasonable solution times and almost perfect scaling for rectilinear wave propagation. The computer simulation makes it possible to observe new phenomena: the break-up of spiral waves caused by intracellular calcium and dynamics and the non-uniformity of the calcium distribution in space during the onset of the spiral wave.

Kogan, B. Y.↗

Compute as Fast as the Engineers Can Think! ULTRAFAST COMPUTING TEAM FINAL REPORT

This report documents findings and recommendations by the Ultrafast Computing Team (UCT). In the period 10-12/98, UCT reviewed design case scenarios for a supersonic transport and a reusable launch vehicle to derive computing requirements necessary for support of a design process with efficiency so radically improved that human thought rather than the computer paces the process. Assessment of the present computing capability against the above requirements indicated a need for further improvement in computing speed by several orders of magnitude to reduce time to solution from tens of hours to seconds in major applications. Evaluation of the trends in computer technology revealed a potential to attain the postulated improvement by further increases of single processor performance combined with massively parallel processing in a heterogeneous environment. However, utilization of massively parallel processing to its full capability will require redevelopment of the engineering analysis and optimization methods, including invention of new paradigms. To that end UCT recommends initiation of a new activity at LaRC called Computational Engineering for development of new methods and tools geared to the new computer architectures in disciplines, their coordination, and validation and benefit demonstration through applications.

Biedron, R. T.↗

Genetically Engineered Microelectronic Infrared Filters

A genetic algorithm is used for design of infrared filters and in the understanding of the material structure of a resonant tunneling diode. These two components are examples of microdevices and nanodevices that can be numerically simulated using fundamental mathematical and physical models. Because the number of parameters that can be used in the design of one of these devices is large, and because experimental exploration of the design space is unfeasible, reliable software models integrated with global optimization methods are examined The genetic algorithm and engineering design codes have been implemented on massively parallel computers to exploit their high performance. Design results are presented for the infrared filter showing new and optimized device design. Results for nanodevices are presented in a companion paper at this workshop.

Cwik, Tom↗