Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed solutions”

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 613 records · Page 34

Calculation of Flow Fields from Bow-Wave Profiles for the Downstream Region of Blunt-Nosed Circular Cylinders in Axial Hypersonic Flight

A method by which known bow-wave profiles may be analyzed to give the flow fields around blunt-nosed cylinders in axial hypersonic flow is presented. In the method, the assumption is made that the pressure distribution curve in a transverse plane is similar to that given by blast- wave theory. Numerical analysis based on the one-dimensional energy and continuity equations then leads to distributions of all the flow variables in the cross section, for either a perfect gas or a real gas. The entire flow field need not be solved. Attention can be confined to any desired station. The critical question is the validity of the above assumption. It is tested for the case of a hemisphere cylinder in flight at 20,000 ft/sec. The flow is analyzed for three stations along the cylindrical afterbody, and found to compare very closely with the results of an exact (inviscid) solution. The assumed form of the pressure distribution occurs at stations as close as 1.2 diameters to the body nose. However, it is suggested that the assumption may not apply this far forward in general, particularly when bodies of nonsmooth contour are considered.

Seiff, Alvin↗

Distributed interactive graphics applications in computational fluid dynamics

Implementation of two interactive, distributed graphics programs used in Computational Fluid Dynamics is discussed. Both programs run on a Cray 2 supercomputer and use a Silicon Graphics Iris workstation as the graphics front-end machine. The hardware and supporting software is from the Numerical Aerodynamic Simulation project. Using this configuration, the supercomputer does all of the numerically intensive work and the workstation allows the user to perform real-time interactive transformations on the displayed data. The first program was written originally as a distributed program which computes particle traces for fluid flow solutions existing on the supercomputer. The second is an older post-processing and plotting program which was modified to run in a distributed mode. Both programs have realized a large increase in capability as a distributed process. Some graphical results are presented.

Rogers, Stuart E.↗

Subsonic Gas Flow Past A Wing Profile

The use of the linearized equations of Chaplygin to calculate the subsonic flow of a gas permits solving the problem of the flow about a wing profile for absence and presence of circulation. The solution is obtained in a practical convenient form that permits finding all the required magnitudes for the gas flow (lift, lift moment velocity distribution over the profile, and critical Mach number). This solution is not expressed in simple closed form; for a certain simplifying assumption, however, the equations of Chaplygin can be reduced to equations with constant coefficients, and solutions are obtained by using only the mathematical apparatus of the theory of functions of a complex variable. The method for simplifying the equations was pointed out by Chaplygin himself. These applied similar equations to the solution of the flow problem and obtained a solution for the case of the absence of circulation.

Christianovich, S. A.↗

Characteristics of Coupled Nongray Radiating Gas Flows with Ablation Product Effects About Blunt Bodies During Planetary Entries

A computational method was developed for the fully-coupled solution of nongray, radiating gas flows with ablation product effects about blunt bodies during planetary entries. The treatment of radiation accounts for molecular band, continuum, and atomic line transitions with a detailed frequency dependence of the absorption coefficient. The ablation of the entry body was solved as part of the solution for a steady-state ablation process. The method was applied by results at typical conditions during entry to Venus. The radiative heating rates along the downstream region of the body can exceed the stagnation point value. The radiative heating to the body is attenuated in the boundary layer at the downstream region of the body and at the stagnation point of the body. A study of the radiation, inviscid flow about spherically capped, conical bodies during planetary entries shows that the nondimensional, radiative heating distributions are nonsimilar with entry conditions. Caution should be exercised in attempting to extrapolate results from known distributions to other entry conditions for which solutions have not yet been obtained.

Sutton, K.↗

Analysis of composite ablators using massively parallel computation

In this work, the feasibility of using massively parallel computation to study the response of ablative materials is investigated. Explicit and implicit finite difference methods are used on a massively parallel computer, the Thinking Machines CM-5. The governing equations are a set of nonlinear partial differential equations. The governing equations are developed for three sample problems: (1) transpiration cooling, (2) ablative composite plate, and (3) restrained thermal growth testing. The transpiration cooling problem is solved using a solution scheme based solely on the explicit finite difference method. The results are compared with available analytical steady-state through-thickness temperature and pressure distributions and good agreement between the numerical and analytical solutions is found. It is also found that a solution scheme based on the explicit finite difference method has the following advantages: incorporates complex physics easily, results in a simple algorithm, and is easily parallelizable. However, a solution scheme of this kind needs very small time steps to maintain stability. A solution scheme based on the implicit finite difference method has the advantage that it does not require very small times steps to maintain stability. However, this kind of solution scheme has the disadvantages that complex physics cannot be easily incorporated into the algorithm and that the solution scheme is difficult to parallelize. A hybrid solution scheme is then developed to combine the strengths of the explicit and implicit finite difference methods and minimize their weaknesses. This is achieved by identifying the critical time scale associated with the governing equations and applying the appropriate finite difference method according to this critical time scale. The hybrid solution scheme is then applied to the ablative composite plate and restrained thermal growth problems. The gas storage term is included in the explicit pressure calculation of both problems. Results from ablative composite plate problems are compared with previous numerical results which did not include the gas storage term. It is found that the through-thickness temperature distribution is not affected much by the gas storage term. However, the through-thickness pressure and stress distributions, and the extent of chemical reactions are different from the previous numerical results. Two types of chemical reaction models are used in the restrained thermal growth testing problem: (1) pressure-independent Arrhenius type rate equations and (2) pressure-dependent Arrhenius type rate equations. The numerical results are compared to experimental results and the pressure-dependent model is able to capture the trend better than the pressure-independent one. Finally, a performance study is done on the hybrid algorithm using the ablative composite plate problem. It is found that there is a good speedup of performance on the CM-5. For 32 CPU's, the speedup of performance is 20. The efficiency of the algorithm is found to be a function of the size and execution time of a given problem and the effective parallelization of the algorithm. It also seems that there is an optimum number of CPU's to use for a given problem.

Shia, David↗

Mathematical model for the Bridgman-Stockbarger crystal growing system

In a major technical breakthrough, a computer model for Bridgman-Stockbarger crystal growth was developed. The model includes melt convection, solute effects, thermal conduction in the ampule, melt, and crystal, and the determination of the curved moving crystal-melt interface. The key to the numerical method is the use of a nonuniform computational mesh which moves with the interface, so that the interface is a mesh surface. In addition, implicit methods are used for advection and diffusion of heat, concentration, and vorticity, for interface movement, and for internal gracity waves. This allows large time-steps without loss of stability or accuracy. Numerical results are presented for the interface shape, temperature distribution, and concentration distribution, in steady-state crystl growth. Solutions are presented for two test cases using water, with two different salts in solution. The two diffusivities differ by a factor of ten, and the concentrations differ by a factor of twenty.

Roberts, G. O.↗

Stability of distributed parameter systems.

A theorem is derived for the stability of solutions of general linear partial differential equations. A norm of the state space in the form of multiple integrals over the spatial domain is used for the Liapunov functional. Theorems and lemmas are also given for linear time-invariant constant coefficient distributed parameter systems, a class of nonlinear distributed parameter systems and for others with a Lure-type nonlinearity. The theorem conditions are similar to those known for corresponding ordinary differential equations but with operators replacing matrices.

Park, K. E.↗

Kinetic theory based new upwind methods for inviscid compressible flows

Two new upwind methods called the Kinetic Numerical Method (KNM) and the Kinetic Flux Vector Splitting (KFVS) method for the solution of the Euler equations have been presented. Both of these methods can be regarded as some suitable moments of an upwind scheme for the solution of the Boltzmann equation provided the distribution function is Maxwellian. This moment-method strategy leads to a unification of the Riemann approach and the pseudo-particle approach used earlier in the development of upwind methods for the Euler equations. A very important aspect of the moment-method strategy is that the new upwind methods satisfy the entropy condition because of the Boltzmann H-Theorem and suggest a possible way of extending the Total Variation Diminishing (TVD) principle within the framework of the H-Theorem. The ability of these methods in obtaining accurate wiggle-free solution is demonstrated by applying them to two test problems.

Deshpande, S. M.↗

Pickup protons and water ions at Comet Halley - Comparisons with Giotto observations

The cometary ion pickup process along the sun-comet line at Comet Halley is investigated using a quasi-linear diffusion model including both pitch angle and energy diffusion, adiabatic compression, and convective motion with the solar wind flow. The model results are compared with energetic ion distributions observed by instruments on board the Giotto spacecraft. The observed power spectrum index of magnetic turbulence (gamma) is 2-2.5. The present simulation shows that when gamma was 2, the calculated proton distributions were much more isotropic than the observed ones. The numerical solutions of the quasi-linear diffusion equations show that the isotropization of the pickup ion distribution, particularly at the pickup velocity, is not complete even close to the bow shock. Given the observed turbulence level, quasi-linear theory yields pickup ion energy distributions that agree with the observed ones quite well and easily produces energetic ions with energies up to hundreds of keV.

Ye, G.↗

A limiting analysis for edge effects in angle-ply laminates

A zeroth order solution for edge effects in angle ply composite laminates using perturbation techniques and a limiting free body approach was developed. The general method of solution for laminates is developed and then applied to the special case of a graphite/epoxy laminate. Interlaminar stress distributions are obtained as a function of the laminate thickness to width ratio h/b and compared to existing numerical results. The solution predicts stable, continuous stress distributions, determines finite maximum tensile interlaminar normal stress for two laminates, and provides mathematical evidence for singular interlaminar shear stresses.

Hsu, P. W.↗

On the equilibrium of a cylindrical plasma supported horizontally by magnetic fields in uniform gravity

The mechanical equilibrium of a cylinder of plasma suspended horizontally by magnetic fields in uniform gravity is examined. A set of exact solutions describing the equilibrium is derived assuming the plasma distribution to be cylindrically symmetric to obtain an exact force balance between plasma pressure, the Lorentz pressure, and gravity in space. The set of solutions considers a case of uniform temperature and cases where the temperature rises from zero at the center of the plasma cylinder to reach a constant asymptotic value outside the cylinder.

Lerche, I.↗

Visual Computing Environment

The Visual Computing Environment (VCE) is a NASA Lewis Research Center project to develop a framework for intercomponent and multidisciplinary computational simulations. Many current engineering analysis codes simulate various aspects of aircraft engine operation. For example, existing computational fluid dynamics (CFD) codes can model the airflow through individual engine components such as the inlet, compressor, combustor, turbine, or nozzle. Currently, these codes are run in isolation, making intercomponent and complete system simulations very difficult to perform. In addition, management and utilization of these engineering codes for coupled component simulations is a complex, laborious task, requiring substantial experience and effort. To facilitate multicomponent aircraft engine analysis, the CFD Research Corporation (CFDRC) is developing the VCE system. This system, which is part of NASA's Numerical Propulsion Simulation System (NPSS) program, can couple various engineering disciplines, such as CFD, structural analysis, and thermal analysis. The objectives of VCE are to (1) develop a visual computing environment for controlling the execution of individual simulation codes that are running in parallel and are distributed on heterogeneous host machines in a networked environment, (2) develop numerical coupling algorithms for interchanging boundary conditions between codes with arbitrary grid matching and different levels of dimensionality, (3) provide a graphical interface for simulation setup and control, and (4) provide tools for online visualization and plotting. VCE was designed to provide a distributed, object-oriented environment. Mechanisms are provided for creating and manipulating objects, such as grids, boundary conditions, and solution data. This environment includes parallel virtual machine (PVM) for distributed processing. Users can interactively select and couple any set of codes that have been modified to run in a parallel distributed fashion on a cluster of heterogeneous workstations. A scripting facility allows users to dictate the sequence of events that make up the particular simulation.

Lawrence, Charles↗

Quantum Key Distribution Applicability to Smart Grid Cybersecurity Systems

To meet the increasing demand for electricity and to have a more reliable and resilient electric grid against conventional and extreme events, grid modernization is more crucial now than ever before. This will require the development and deployment of devices that provide advanced communication capabilities. The overall efficiency, reliability, and resilience of the smart grid will be inextricably linked to the exchange of information between these devices. Unfortunately, the increased information flow will increase the potential attack surface and introduce new vulnerabilities. While a smarter grid will depend critically on information flow, these benefits will be accrued only if that information can be protected. Nowadays, information is secured in smart grids primarily through cryptography. However, with the increasing number of sophisticated attacks as well as the increasing computational power, the security of the “classical” cryptographic algorithms is threatened. Quantum information science offers solutions to this problem, specifically quantum key distribution (QKD), which provides a means for the generation and secure distribution of symmetric cryptographic keys. The security of QKD stems ultimately from the very nature of quantum physics. In this paper, we investigate the applicability of QKD to the various smart grid sectors and specific use cases. We have identified 18 smart grid use cases of interest for QKD suitability together with 7 QKD factors used for the assessment of the various use cases. For each use case, the impact to security of the loss of confidentiality, integrity, and/or availability is specified. In addition, the suitability of QKD is assessed for each use case with respect to multiple factors.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Precipitation growth in convective clouds

Analytical solutions to the equations of both the growth and motion of hailstones in updrafts and of cloud water contents which vary linearly with height were used to investigate hail growth in a model cloud. A strong correlation was found between the hail embyro starting position and its trajectory and final size. A simple model of the evolution of particle size distribution by coalescence and spontaneous and binary disintegrations was formulated. Solutions for the mean mass of the distribution and the equilibrium size distribution were obtained for the case of constant collection kernel and disintegration parameters. Azimuthal scans of Doppler velocity at a number of elevation angles were used to calculate high resolution vertical profiles of particle speed and horizontal divergence (the vertical air velocity) in a region of widespread precipitation trailing a mid-latitude squall line.

Srivastava, R. C.↗

Normal shock wave-turbulent boundary layer interactions in transonic flow near separation

The problem of a normal shock wave impinging on a flat plate, turbulent boundary layer is considered for the case where the external flow is transonic. Asymptotic methods are employed. It is shown that there are two outer regions, including the outer part of the boundary layer and the external flow, in which inviscid flow governing equations hold, and two regions near the wall, in which Reynolds and/or viscous stresses need be taken into account. The solutions in the outer regions lead to the pressure distribution on the wall, for which an analytical expression is presented, valid under those conditions when separation is imminent but has not yet occurred. The solutions valid in the inner regions lead to a separation criterion.

Adamson, T. C., Jr.↗

Distributed parallel processing applied to an implicit multigrid Euler/Navier-Stokes algorithm

An implicit multigrid algorithm for the solution of the Euler and Navier-Stokes equations has been implemented within the framework of multiple block-structured grids in which the physical domain is spatially decomposed into several blocks and the solution is advanced in parallel on each block. Utilities have been developed to implement such a scheme in a distributed computing environment. The multi-block algorithm is designed so that the explicit residual calculation is identical to that of single-block scheme, and therefore converged solutions for both schemes must be the same. To accelerate convergence, synchronous and asynchronous multigrid strategies are implemented. Significant speedups have been achieved in a multiple processor environment, while convergence rates similar to those of the single-block scheme are observed.

Tysinger, T. L.↗

Closed-form solution of decomposable stochastic models

Markov and semi-Markov processes are increasingly being used in the modeling of complex reconfigurable systems (fault tolerant computers). The estimation of the reliability (or some measure of performance) of the system reduces to solving the process for its state probabilities. Such a model may exhibit numerous states and complicated transition distributions, contributing to an expensive and numerically delicate solution procedure. Thus, when a system exhibits a decomposition property, either structurally (autonomous subsystems), or behaviorally (component failure versus reconfiguration), it is desirable to exploit this decomposition in the reliability calculation. In interesting cases there can be failure states which arise from non-failure states of the subsystems. Equations are presented which allow the computation of failure probabilities of the total (combined) model without requiring a complete solution of the combined model. This material is presented within the context of closed-form functional representation of probabilities as utilized in the Symbolic Hierarchical Automated Reliability and Performance Evaluator (SHARPE) tool. The techniques adopted enable one to compute such probability functions for a much wider class of systems at a reduced computational cost. Several examples show how the method is used, especially in enhancing the versatility of the SHARPE tool.

Sjogren, Jon A.↗

Investigations of the dynamical evolution of protoplanetary nebulae

The spectral energy distributions (SED's) of all the known Class I (protostellar) sources in the Taurus molecular cloud were modelled. The Tereby, Shu, & Cassen (1984, TSC) density distribution for a rotating, infalling envelope was adopted. The radiative equilibrium temperature distribution from the spherical average of the TSC density distribution was calculated. The resulting spherically-symmetric temperature distribution then provides the source function to obtain the emergent spectrum at a given inclination angle i from the formal solution of the transfer equation, using the exact density (opacity) distribution. Results showed that the SED's of the protostar candidates in Taurus can be reproduced with TSC models having infall rates close to the values predicted by the theory of isothermal cloud collapse. Flat Spectrum T Tauri Stars were studied. The mid- to far-infrared fluxes of 'flat spectrum' T Tauri stars can be explained by radiative equilibrium emission from infalling dusty envelopes. This explanation indicated that models employing 'active' disks, in which the temperature distribution is a parameterized power law, should be invoked with caution. Infall also naturally explains the scattered light nebulae detected around many flat-spectrum sources. Results showed that the SED's of the protostar candidates in Taurus can be reproduced with TSC models having infall rates close to the values predicted by the theory of isothermal cloud collapse. Flattened infalling envelope models are also being studied.

Hartmann, Lee W.↗