SEARCH · Search NASA

Results for “numerical algorithm”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1993

A generic efficient adaptive grid scheme for rocket propulsion modeling

The objective of this research is to develop an efficient, time-accurate numerical algorithm to discretize the Navier-Stokes equations for the predictions of internal one-, two-dimensional and axisymmetric flows. A generic, efficient, elliptic adaptive grid generator is implicitly coupled with the Lower-Upper factorization scheme in the development of ALUNS computer code. The calculations of one-dimensional shock tube wave propagation and two-dimensional shock wave capture, wave-wave interactions, shock wave-boundary interactions show that the developed scheme is stable, accurate and extremely robust. The adaptive grid generator produced a very favorable grid network by a grid speed technique. This generic adaptive grid generator is also applied in the PARC and FDNS codes and the computational results for solid rocket nozzle flowfield and crystal growth modeling by those codes will be presented in the conference, too. This research work is being supported by NASA/MSFC.

Mo, J. D.↗

NASA NTRS1993

Modeling of Non-Spherical Droplet Dynamics

A two-dimensional time-dependent computer code based on the modified Arbitrary Lagrangian Eulerian (ALE) technique, has been developed to simulate non-spherical droplet dynamics and evaporation under convective flows at real rocket combustion chamber conditions. The equations of mass, momentum, energy and species are simultaneously solved for both liquid and gas phases with an accurate dynamic interface tracking. The jump boundary conditions across the deforming droplet surface are obtained by applying the integral forms of conservation of mass, momentum, and energy. At each time step, the interface geometry and flow properties at the droplet surface are implicitly solved by satisfying the interface boundary conditions. A Lagrangian technique was developed to track the arbitrarily moving interface between the liquid droplet and the external gas. An elliptic grid generator is adopted to dynamically reconstruct grids both inside and outside the droplet surface. This code has been used to study droplet oscillation, droplet deformation/breakup, nonspherical droplet evaporation in both low and high pressure convective flows. This presentation briefly describes the numerical algorithm for modeling of the nonspherical droplet dynamics and demonstrates the representative simulation results of nonspherical droplet evaporation at low and high pressure convective flows. Potential applications of this code to rocket combustor design and performance predictions are discussed.

Deng, Zheng-Tao↗

NASA NTRS1993

Code Validation Study for Base Flows

New and old rocket launch concepts recommend the clustering of motors for improved lift capability. The flowfield of the base region of the rocket is very complex and can contain high temperature plume gases. These hot gases can cause catastrophic problems if not adequately designed for. To assess the base region characteristics, advanced computational fluid dynamics (CFD) is being used. As a precursor to these calculations the CFD code requires validation on base flows. The primary objective of this code validation study was to establish a high level of confidence in predicting base flows with the USA CFD code. USA has been extensively validated for fundamental flows and other applications. However, base heating flows have a number of unique characteristics so it was necessary to extend the existing validation for this class of problems. In preparation for the planned NLS 1.5 Stage base heating analysis, six case sets were studied to extend the USA code validation data base. This presentation gives a cursive review of three of these cases. The cases presented include a 2D axi-symmetric study, a 3D real nozzle study, and a 3D multi-species study. The results of all the studies show good general agreement with data with no adjustments to the base numerical algorithms or physical models in the code. The study proved the capability of the USA code for modeling base flows within the accuracy of available data.

Ascoli, Edward P.↗

NASA NTRS1995

Time accurate application of the MacCormack 2-4 scheme on massively parallel computers

Many recent computational efforts in turbulence and acoustics research have used higher order numerical algorithms. One popular method has been the explicit MacCormack 2-4 scheme. The MacCormack 2-4 scheme is second order accurate in time and fourth order accurate in space, and is stable for CFL's below 2/3. Current research has shown that the method can give accurate results but does exhibit significant Gibbs phenomena at sharp discontinuities. The impact of adding Jameson type second, third, and fourth order artificial viscosity was examined here. Category 2 problems, the nonlinear traveling wave and the Riemann problem, were computed using a CFL number of 0.25. This research has found that dispersion errors can be significantly reduced or nearly eliminated by using a combination of second and third order terms in the damping. Use of second and fourth order terms reduced the magnitude of dispersion errors but not as effectively as the second and third order combination. The program was coded using Thinking Machine's CM Fortran, a variant of Fortran 90/High Performance Fortran, and was executed on a 2K CM-200. Simple extrapolation boundary conditions were used for both problems.

Hudson, Dale A.↗

NASA NTRS1995

FPCAS3D User's guide: A three dimensional full potential aeroelastic program, version 1

The FPCAS3D computer code has been developed for aeroelastic stability analysis of bladed disks such as those in fans, compressors, turbines, propellers, or propfans. The aerodynamic analysis used in this code is based on the unsteady three-dimensional full potential equation which is solved for a blade row. The structural analysis is based on a finite-element model for each blade. Detailed explanations of the aerodynamic analysis, the numerical algorithms, and the aeroelastic analysis are not given in this report. This guide can be used to assist in the preparation of the input data required by the FPCAS3D code. A complete description of the input data is provided in this report. In addition, six examples, including inputs and outputs, are provided.

Bakhle, Milind A.↗

NASA NTRS1995

Film temperatures in the presence of cavitation

Numerical algorithms are developed and implemented for the treatment of laminar lubricating-film temperatures associated with cavitated regions. The reformation front, with its film-content discontinuity and flow reversal, is given special attention. Computational economy is achieved through the use of Lobatto-point locations for flow-property determinations.

Elrod, Harold G.↗

NASA NTRS1991

Recent advances in Euler and Navier-Stokes methods for calculating helicopter rotor aerodynamics and acoustics

This paper outlines some recent advances in the application of the Euler and Navier-Stokes computational fluid dynamics methods to analyze nonlinear problems of helicopter aerodynamics and acoustics. A complete flowfield simulation of helicopters is currently not feasible with these methods. However, the use of the state-of-the-art numerical algorithms in conjunction with powerful supercomputers, like the Cray-2, have enabled notable progress to be made in modeling several individual components of this complex flow in hover and forward flight.

Srinivasan, G. R.↗

NASA NTRS1993

Use of networked workstations for parallel nonlinear structural dynamic simulations of rotating bladed-disk assemblies

The principal objective of this research is to investigate, develop and demonstrate coarse-grained, parallel-processing strategies for nonlinear dynamic simulations for rotating bladed-disk assemblies. The parallel -processing strategies addressed include numerical algorithms for parallel nonlinear solutions and techniques to effect load balancing among processors. The parallel environment employed is a distributed-memory, coarse-grained one consisting of networked workstations. A parallel explicit time integration method has been implemented for transient nonlinear solutions of rotationg bladed-disk assemblies. Automatic domain partitioning techniques have been investigated for load balancing among processors. Advanced computing environments, data structures and interactive computer graphics all contribute to an integrated parallel finite element analysis system to facilitate more efficient and powerful dynamic simulations.

Hsieh, Shang-Hsien↗

NASA NTRS1995

Time-domain approach for the transient responses in stratified viscoelastic Earth models

We have developed the numerical algorithm for the computation of transient viscoelastic responses in the time domain for a radially stratified Earth model. Stratifications in both the elastic parameters and the viscosity profile have been considered. The particular viscosity profile employed has a viscosity maximum with a constrast of O(100) in the mid lower mantle. The distribution of relaxation times reveals the presence of a continuous spectrum situated between O(100) and O(exp 4) years. The principal mode is embedded within this continuous spectrum. From this initial-value approach we have found that for the low degree harmonics the non-modal contributions are comparable to the modal contributions. For this viscosity model the differences between the time-domain and normal-mode results are found to decrease strongly with increasing angular order. These calculations also show that a time-dependent effective relaxation time can be defined, which can be bounded by the relaxation times of the principal modes.

Hanyk, L.↗

NASA NTRS1995

Electromagnetic scattering analysis of a three-dimensional-cavity-backed aperture in an infinite ground plane using a combined finite element method/method of moments approach

A combined finite element method/method of moments (FEM/MoM) approach is used to analyze the electromagnetic scattering properties of a three-dimensional-cavity-backed aperture in an infinite ground plane. The FEM is used to formulate the fields inside the cavity, and the MoM (with subdomain bases) in both spectral and spatial domains is used to formulate the fields above the ground plane. Fields in the aperture and the cavity are solved using a system of equations resulting from the combination of the FEM and the MoM. By virtue of the FEM, this combined approach is applicable to all arbitrarily shaped cavities with inhomogeneous material fillings, and because of the subdomain bases used in the MoM, the apertures can be of any arbitrary shape. This approach leads to a partly sparse and partly full symmetric matrix, which is efficiently solved using a biconjugate gradient algorithm. Numerical results are presented to validate the analysis.

Reddy, C. J.↗

NASA NTRS1995

Study of nonlinear MHD equations governing the wave propagation in twisted coronal loops

The solar corona, modelled by a low beta, resistive plasma slab, sustains MHD wave propagations due to shearing footpoint motions in the photosphere. By using a numerical algorithm the excitation and nonlinear development of MHD waves in twisted coronal loops are studied. The plasma responds to the footpoint motion by sausage waves if there is no twist. The twist in the magnetic field of the loop destroys initially developed sausage-like wave modes and they become kinks. The transition from sausage to kink modes is analyzed. The twist brings about mode degradation producing high harmonics and this generates more complex fine structures. This can be attributed to several local extrema in the perturbed velocity profiles. The Alfven wave produces remnants of the ideal 1/x singularity both for zero and non-zero twist and this pseudo-singularity becomes less pronounced for larger twist. The effect of nonlinearity is clearly observed by changing the amplitude of the driver by one order of magnitude. The magnetosonic waves also exhibit smoothed remnants of ideal logarithmic singularities when the frequency of the driver is correctly chosen. This pseudo-singularity for fast waves is absent when the coronal loop does not undergo any twist but becomes pronounced when twist is included. On the contrary, it is observed for slow waves even if there is no twist. Increasing the twist leads to a higher heating rate of the loop. The larger twist shifts somewhat uniformly distributed heating to layers inside the slab corresponding to peaks in the magnetic field strength.

Parhi, S.↗

NASA NTRS1997

Parallel Computation of Unsteady Flows on a Network of Workstations

Parallel computation of unsteady flows requires significant computational resources. The utilization of a network of workstations seems an efficient solution to the problem where large problems can be treated at a reasonable cost. This approach requires the solution of several problems: 1) the partitioning and distribution of the problem over a network of workstation, 2) efficient communication tools, 3) managing the system efficiently for a given problem. Of course, there is the question of the efficiency of any given numerical algorithm to such a computing system. NPARC code was chosen as a sample for the application. For the explicit version of the NPARC code both two- and three-dimensional problems were studied. Again both steady and unsteady problems were investigated. The issues studied as a part of the research program were: 1) how to distribute the data between the workstations, 2) how to compute and how to communicate at each node efficiently, 3) how to balance the load distribution. In the following, a summary of these activities is presented. Details of the work have been presented and published as referenced.

Source record↗

NASA NTRS1999

Parallelizing OVERFLOW: Experiences, Lessons, Results

The computer code OVERFLOW is widely used in the aerodynamic community for the numerical solution of the Navier-Stokes equations. Current trends in computer systems and architectures are toward multiple processors and parallelism, including distributed memory. This report describes work that has been carried out by the author and others at Ames Research Center with the goal of parallelizing OVERFLOW using a variety of parallel architectures and parallelization strategies. This paper begins with a brief description of the OVERFLOW code. This description includes the basic numerical algorithm and some software engineering considerations. Next comes a description of a parallel version of OVERFLOW, OVERFLOW/PVM, using PVM (Parallel Virtual Machine). This parallel version of OVERFLOW uses the manager/worker style and is part of the standard OVERFLOW distribution. Then comes a description of a parallel version of OVERFLOW, OVERFLOW/MPI, using MPI (Message Passing Interface). This parallel version of OVERFLOW uses the SPMD (Single Program Multiple Data) style. Finally comes a discussion of alternatives to explicit message-passing in the context of parallelizing OVERFLOW.

Jespersen, Dennis C.↗

NASA NTRS2000

Heating and Large Scale Dynamics of the Solar Corona

The effort was concentrated in the areas: coronal heating mechanism, unstructured adaptive grid algorithms, numerical modeling of magnetic reconnection in the MRX experiment: effect of toroidal magnetic field and finite pressure, effect of OHMIC heating and vertical magnetic field, effect of dynamic MESH adaption.

Schnack, Dalton D.↗

NASA NTRS2000

VCSEL Applications and Simulation

This viewgraph presentation gives an overview of Vertical Cavity Surface Emitting Laser (VCSEL) simulation and its applications. Details are given on the optical interconnection in information technology of VCSEL, the formulation of the simulation, its numeric algorithm, and the computational results.

Cheung, Samson↗

NASA NTRS1994

Navier-Stokes Solutions of Rotor and Rotor-Body Flows

This paper will review the advances made recently in the Navier-Stokes CFD methods to simulate aerodynamics and aeroacoustics of helicopter rotors and rotor-body flows. Although a complete flowfield simulation of full helicopter is currently not feasible with these methods, impressive gains have been made in analyzing individual components of this complex problem in a very detailed manner. The use of the state-of-the-art numerical algorithms in solution methods, in conjunction with powerful supercomputers, like the Cray-2, have enabled noticeable progress to be made in modeling viscous-inviscid interactions, blade-vortex interactions, tip-vortex: simulation and wake effects, as well as high speed impulsive noise in hover and forward flight for isolated rotor blades. This paper will critically evaluate the presently available Euler and Navier-Stokes methods, both finite-difference and finite volume methods using structured and unstructured grids for helicopter applications for accuracy, suitability, and computational efficiency. The review will also include the recent progress made using overset grids to model rotor-body flows. All the material for this review will be drawn from the published material shown below.

Srinivasan, G. R.↗

NASA NTRS1995

High Performance Parallel Computational Nanotechnology

At a recent press conference, NASA Administrator Dan Goldin encouraged NASA Ames Research Center to take a lead role in promoting research and development of advanced, high-performance computer technology, including nanotechnology. Manufacturers of leading-edge microprocessors currently perform large-scale simulations in the design and verification of semiconductor devices and microprocessors. Recently, the need for this intensive simulation and modeling analysis has greatly increased, due in part to the ever-increasing complexity of these devices, as well as the lessons of experiences such as the Pentium fiasco. Simulation, modeling, testing, and validation will be even more important for designing molecular computers because of the complex specification of millions of atoms, thousands of assembly steps, as well as the simulation and modeling needed to ensure reliable, robust and efficient fabrication of the molecular devices. The software for this capacity does not exist today, but it can be extrapolated from the software currently used in molecular modeling for other applications: semi-empirical methods, ab initio methods, self-consistent field methods, Hartree-Fock methods, molecular mechanics; and simulation methods for diamondoid structures. In as much as it seems clear that the application of such methods in nanotechnology will require powerful, highly powerful systems, this talk will discuss techniques and issues for performing these types of computations on parallel systems. We will describe system design issues (memory, I/O, mass storage, operating system requirements, special user interface issues, interconnects, bandwidths, and programming languages) involved in parallel methods for scalable classical, semiclassical, quantum, molecular mechanics, and continuum models; molecular nanotechnology computer-aided designs (NanoCAD) techniques; visualization using virtual reality techniques of structural models and assembly sequences; software required to control mini robotic manipulators for positional control; scalable numerical algorithms for reliability, verifications and testability. There appears no fundamental obstacle to simulating molecular compilers and molecular computers on high performance parallel computers, just as the Boeing 777 was simulated on a computer before manufacturing it.

Saini, Subhash↗

NASA NTRS2001

Vorticity Transport and Wave Emission in the Protoplanetary Nebula

Higher order numerical algorithms (4th order in time, 3rd order in space) are applied to the Euler/Energy equations and are used to examine vorticity transport and wave motion in a non-self gravitating, initially isentropic Keplerian disk. In this talk we will examine the response of the nebula to an isolated vortex with a circulation about equal to the rotation rate of Jupiter. The vortex is located on the 4 AU circle and the nebula is simulated from 1 to 24 AU. We show that the vortex emits pressure-supported density and Rossby-type wave packets before it decays within a few orbits. The acoustic density waves evolve into weak (non entropy preserving) shock waves that propagate over the entire disk. The Rossby waves remain in the vicinity of the initial vortex disturbance, but are rapidly damped. Temporal frequencies and spatial wavenumbers are derived using the simulation data and compared with analytical dispersion relations from the linearized Euler/Energy equations.

Davis, S. S.↗