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 1,045 records · Page 58

On the stabilization of three-dimensional boundary layers by suction and cooling

A significant reduction in the drag of transonic aircraft can be achieved by using an active method of boundary-layer control to maintain laminar flow on the aerodynamic surfaces. The stabilizing influence of suction for both favorable and adverse pressure gradients is demonstrated by means of the self-similar incompressible three-dimensional boundary layers on yawed wedges. Profile instability is measured by the maximum amplification rate of fixed-frequency disturbances computed according to linearized, locally parallel, spatial stability theory. Suction is found to be more effective in controlling Tollmien-Schlichting instability than stationary cross-flow disturbances. The effectiveness of surface cooling as a method of stabilization is compared with suction for the boundary layers on a transonic 35 deg swept wing of infinite span. It is found from compressible stability theory that when the surface is cooled to a uniform temperature such that the maximum cross-flow velocity is reduced by the same amount as with a uniform suction distribution, the stabilizing effect of cooling on cross-flow disturbances is less than with suction.

Mack, L. M.↗

Trajectory optimization for real-time guidance. I - Time-varying LQR on a parallel processor

A key algorithmic element of a real-time trajectory optimization hardware/software implementation, the quadratic program (QP) solver element, is presented. The purpose of the effort is to make nonlinear trajectory optimization fast enough to provide real-time commands during guidance of a vehicle such as an aeromaneuvering orbiter. Many methods of nonlinear programming require the solution of a QP at each iteration. In the trajectory optimization case the QP has a special dynamic programming structure, a LQR-like structure. QP algorithm speed is increased by taking advantage of this special structure and by parallel implementation.

Psiaki, Mark L.↗

An Overview of the NCC Spray/Monte-Carlo-PDF Computations

This paper advances the state-of-the-art in spray computations with some of our recent contributions involving scalar Monte Carlo PDF (Probability Density Function), unstructured grids and parallel computing. It provides a complete overview of the scalar Monte Carlo PDF and Lagrangian spray computer codes developed for application with unstructured grids and parallel computing. Detailed comparisons for the case of a reacting non-swirling spray clearly highlight the important role that chemistry/turbulence interactions play in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that some of the combustion occurs in a predominantly premixed-flame environment and the rest in a predominantly diffusion-flame environment. However, the non-PDF solution predicts wrongly for the combustion to occur in a vaporization-controlled regime. Near the premixed flame, the Monte Carlo particle temperature distribution shows two distinct peaks: one centered around the flame temperature and the other around the surrounding-gas temperature. Near the diffusion flame, the Monte Carlo particle temperature distribution shows a single peak. In both cases, the computed PDF's shape and strength are found to vary substantially depending upon the proximity to the flame surface. The results bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations. Finally, we end the paper by demonstrating the computational viability of the present solution procedure for its use in 3D combustor calculations by summarizing the results of a 3D test case with periodic boundary conditions. For the 3D case, the parallel performance of all the three solvers (CFD, PDF, and spray) has been found to be good when the computations were performed on a 24-processor SGI Origin work-station.

Raju, M. S.↗

The Dynamics of Some Iterative Implicit Schemes

The global asymptotic nonlinear behavior of some standard iterative procedures in solving nonlinear systems of algebraic equations arising from four implicit linear multistep methods (LMMs) in discretizing 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations is analyzed using the theory of dynamical systems. With the aid of parallel Connection Machines (CM-2 and CM-5), the associated bifurcation diagrams as a function of the time step, and the complex behavior of the associated 'numerical basins of attraction' of these iterative implicit schemes are revealed and compared. Studies showed that all of the four implicit LMMs exhibit a drastic distortion and segmentation but less shrinkage of the basin of attraction of the true solution than standard explicit methods. The numerical basins of attraction of a noniterative implicit procedure mimic more closely the basins of attraction of the differential equations than the iterative implicit procedures for the four implicit LMMs.

Yee, H. C.↗

Radiative interactions in multi-dimensional chemically reacting flows using Monte Carlo simulations

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The amount and transfer of the emitted radiative energy in a finite volume element within a medium are considered in an exact manner. The spectral correlation between transmittances of two different segments of the same path in a medium makes the statistical relationship different from the conventional relationship, which only provides the non-correlated results for nongray methods is discussed. Validation of the Monte Carlo formulations is conducted by comparing results of this method of other solutions. In order to further establish the validity of the MCM, a relatively simple problem of radiative interactions in laminar parallel plate flows is considered. One-dimensional correlated Monte Carlo formulations are applied to investigate radiative heat transfer. The nongray Monte Carlo solutions are also obtained for the same problem and they also essentially match the available analytical solutions. the exact correlated and non-correlated Monte Carlo formulations are very complicated for multi-dimensional systems. However, by introducing the assumption of an infinitesimal volume element, the approximate correlated and non-correlated formulations are obtained which are much simpler than the exact formulations. Consideration of different problems and comparison of different solutions reveal that the approximate and exact correlated solutions agree very well, and so do the approximate and exact non-correlated solutions. However, the two non-correlated solutions have no physical meaning because they significantly differ from the correlated solutions. An accurate prediction of radiative heat transfer in any nongray and multi-dimensional system is possible by using the approximate correlated formulations. Radiative interactions are investigated in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The governing equations are based on the fully elliptic Navier-Stokes equations. Chemical reaction mechanisms were described by a finite rate chemistry model. The correlated Monte Carlo method developed earlier was employed to simulate multi-dimensional radiative heat transfer. Results obtained demonstrate that radiative effects on the flowfield are minimal but radiative effects on the wall heat transfer are significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and nozzle size on the radiative and conductive wall fluxes.

Liu, Jiwen↗

UCLA parallel PIC framework

The UCLA Parallel PIC Framework (UPIC) has been developed to provide trusted components for the rapid construction of new, parallel Particle-in-Cell (PIC) codes. The Framework uses object-based ideas in Fortran95, and is designed to provide support for various kinds of PIC codes on various kinds of hardware. The focus is on student programmers. The Framework supports multiple numerical methods, different physics approximations, different numerical optimizations and implementations for different hardware. It is designed with "defensive" programming in mind, meaning that it contains many error checks and debugging helps. Above all, it is designed to hide the complexity of parallel processing. It is currently being used in a number of new Parallel PIC codes.

Norton, Charles D.↗

Monte Carlo investigation of thrust imbalance of solid rocket motor pairs

The Monte Carlo method of statistical analysis is used to investigate the theoretical thrust imbalance of pairs of solid rocket motors (SRMs) firing in parallel. Sets of the significant variables are selected using a random sampling technique and the imbalance calculated for a large number of motor pairs using a simplified, but comprehensive, model of the internal ballistics. The treatment of burning surface geometry allows for the variations in the ovality and alignment of the motor case and mandrel as well as those arising from differences in the basic size dimensions and propellant properties. The analysis is used to predict the thrust-time characteristics of 130 randomly selected pairs of Titan IIIC SRMs. A statistical comparison of the results with test data for 20 pairs shows the theory underpredicts the standard deviation in maximum thrust imbalance by 20% with variability in burning times matched within 2%. The range in thrust imbalance of Space Shuttle type SRM pairs is also estimated using applicable tolerances and variabilities and a correction factor based on the Titan IIIC analysis.

Sforzini, R. H.↗

Unsteady supersonic aerodynamic theory by the method of potential gradient

A generalized solution of the hyperbolic wave equation has been derived. The method used has been further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals have been evaluated in closed form, while numerical integration methods are suggested for treating more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient because line integrations along the other two sides need only be considered. Furthermore, all the integrands vanish on the Mach cone and the need for determining the hyperbolic curves of intersection of the cone with the lifting surface is avoided. Generalized aerodynamic coefficients for three AGARD planforms have been calculated and compared with the available results.

Jones, W. P.↗

Measurement of transient strain and surface temperature on simulated turbine blades using noncontacting techniques

Noncontacting techniques were used to measure strain and temperature in thermally cycled simulated turbine blades. An electro-optical extensometer was used to measure the displacement between parallel targets mounted on the leading edge of the blades throughout a complete heating and cooling cycle. An infrared photographic pyrometry method was used to measure blade steady state surface temperature. The blade was cyclically heated and cooled by moving it into and out of a Mach 1 hot-gas stream. Transient leading edge strain and steady state surface temperature distributions are presented for blades of three different configurations.

Calfo, F. D.↗

A decentralized square root information filter/smoother

A number of developments has recently led to a considerable interest in the decentralization of linear least squares estimators. The developments are partly related to the impending emergence of VLSI technology, the realization of parallel processing, and the need for algorithmic ways to speed the solution of dynamically decoupled, high dimensional estimation problems. A new method is presented for combining Square Root Information Filters (SRIF) estimates obtained from independent data sets. The new method involves an orthogonal transformation, and an information matrix filter 'homework' problem discussed by Schweppe (1973) is generalized. The employed SRIF orthogonal transformation methodology has been described by Bierman (1977).

Bierman, G. J.↗

Capacitive Displacement Sensor With Frequency Readout

Simple displacement-measuring circuit senses capacitance between two parallel conducting plates and produces output signal, with frequency proportional to distance between plates. Principle of circuit provides advantages over other methods because of frequency-encoded output and high linearity. Used to measure displacements.

Fritsch, Klaus↗

Monte-Carlo modeling of polar wind photoelectron distributions with anomalous heat flux

In situ measurements above the polar cap by the DE satellites show asymmetric field-aligned electron velocity distributions in the photoelectron energy range 5-60 eV. The formation of these anisotropic distributions is described by a collisional kinetic equation that includes the effects of Coulomb collisions (via a Fokker-Planck operator), those of a parallel electric field, and the magnetic mirror effects of the geomagnetic field. This equation is solved using a Monte Carlo simulation method. The simulation illustrates the formation of the portion of the observed electron distributions originating in the ionosphere and makes it possible to evaluate the significance of the suprathermal electrons at higher altitude. Both observations and calculations indicate that these electrons possess an energy flux which, by its effect on the polar-wind electric field, may influence the dynamics of the polar-wind outflow.

Yasseen, F.↗

Tomographic Measurement Of Laser-Bored Holes

Nondestructive technique detects internal variations in arrays of small holes. Inspection method checks laser-bored holes for accuracy. Combines computed tomography and digital laminography. Both types of views made at many parallel planes within plate. System prints out tables of measured and standard deviation of diameter at all planes for each hole.

Willenberg, James D.↗

Apparatus Splits Glass Tubes Longitudinally

Tubes split into half cylinders by hot-wire/thermal-shock method. Tube to be cut placed on notched jig in apparatus. Nichrome wire stretched between arms of pivoted carriage and oriented parallel to notch. Wire heated by electrical current while resting on tube. After heating for about 1 minute for each millimeter of thickness of glass, tube quenched in water and split by resulting thermal shock. Apparatus used to split tubes in sizes ranging from 3/8 in. in diameter by 1 in. long to 1 1/2 in. in diameter by 4 in. long.

Shaw, Ernest↗

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The primary objective of this research is to extend current capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first two years of this research have been concentrated on a priori investigations of single-point Probability Density Function (PDF) methods for providing subgrid closures in reacting turbulent flows. In the efforts initiated in the third year, our primary focus has been on performing actual LES by means of PDF methods. The approach is based on assumed PDF methods and we have performed extensive analysis of turbulent reacting flows by means of LES. This includes simulations of both three-dimensional (3D) isotropic compressible flows and two-dimensional reacting planar mixing layers. In addition to these LES analyses, some work is in progress to assess the extent of validity of our assumed PDF methods. This assessment is done by making detailed companions with recent laboratory data in predicting the rate of reactant conversion in parallel reacting shear flows. This report provides a summary of our achievements for the first six months of the third year of this program.

Givi, P.↗

Aerodynamic optimization studies on advanced architecture computers

The approach to carrying out multi-discipline aerospace design studies in the future, especially in massively parallel computing environments, comprises of choosing (1) suitable solvers to compute solutions to equations characterizing a discipline, and (2) efficient optimization methods. In addition, for aerodynamic optimization problems, (3) smart methodologies must be selected to modify the surface shape. In this research effort, a 'direct' optimization method is implemented on the Cray C-90 to improve aerodynamic design. It is coupled with an existing implicit Navier-Stokes solver, OVERFLOW, to compute flow solutions. The optimization method is chosen such that it can accomodate multi-discipline optimization in future computations. In the work , however, only single discipline aerodynamic optimization will be included.

Chawla, Kalpana↗

Some Experiences with Nonoverlapping Schur Complement Parallel Preconditioning for CFD Calculations

In this work we consider solving matrices which arise from the discretization of advection-diffusion field equations on arbitrary triangulated domains using stabilized numerical methods. The talk will discuss several candidate matrix preconditioning algorithms based on the 2 x 2 block factorization induced by an apriori partitioning of the triangulated domain. Application of the 2 x 2 block preconditioner requires the formation and inversion of the Schur complement submatrix. We consider several strategies for simplifying this task: incomplete Schur complement factorizations, drop tolerance element filling, Schur complement probing, and localized Schur complement inversion. Numerical results will be shown comparing performance and efficiency of these approximations. The matrix preconditioner has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The remainder of the talk will show numerous examples in CFD to demonstrate the efficiency and robustness of the techniques.

Barth, Timothy J.↗

Noninvasive Medical Diagnostics & Treatment Using Ultrasonics

In parallel to the industrial application of NDE to flaw detection and material property determination, the medical community has succesfully adapted such methods to the noninvasaive diagnostics and treatment of many conditions and disorders of the human body.

NDE noninvasive diagnostics destruction of blood c↗