SEARCH · Search NASA

Results for “convection simulation parallel computation”

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 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Large-eddy simulations of compressible convection on massively parallel computers

We report preliminary implementation of the large-eddy simulation (LES) technique in 2D simulations of compressible convection carried out on the CM-2 massively parallel computer. The convective flow fields in our simulations possess structures similar to those found in a number of direct simulations, with roll-like flows coherent across the entire depth of the layer that spans several density scale heights. Our detailed assessment of the effects of various subgrid scale (SGS) terms reveals that they may affect the gross character of convection. Yet, somewhat surprisingly, we find that our LES solutions, and another in which the SGS terms are turned off, only show modest differences. The resulting 2D flows realized here are rather laminar in character, and achieving substantial turbulence may require stronger forcing and less dissipation.

Xie, Xin↗

NASA NTRS1995

Parallel Computation for Natural Convection in Cavities

Parallel computation for thermal convective flows in cavities with adiabatic horizontal boundaries and driven by differential heating of the two vertical end walls, is investigated using supercomputers. A parallel computation code has been implemented using a finite-difference method with a multigrid elliptic solver and a Dufort-Frankel scheme. The domain decomposition techniques are discussed in detail. The parallel code is numerically stable, computationally efficient, and portable to various parallel architectures which support either PVM or NX libraries for communications. Finally, numerical results for various Rayleigh numbers and Prandtl numbers are presented.

convection simulation parallel computation↗

NASA NTRS1994

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗

NASA NTRS1996

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

This research program dealt with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in January 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled three-component problem were developed during 1994 and 1995. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor tor parallel versions of ENG10 was developed. During 1995 and 1996 we developed the capability tor the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames. Benchmark results were presented at the 1196 Computational Aeroscience meeting.

Felippa, C. A.↗

NASA NTRS1981

Numerical simulation of natural convection in a spherical container due to cooling at the center (idealization of the Lal/Kroes experiment)

Natural convection in a spherical container with cooling at the center was numerically simulated using a numerical fluid dynamics computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes crystal growing experiment to be performed on Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗

NASA NTRS1981

Manufacturing in space: Fluid dynamics numerical analysis

Natural convection in a spherical container with cooling at the center was numerically simulated using the Lockheed-developed General Interpolants Method (GIM) numerical fluid dynamic computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes growing experiments to be performed on board Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. For a temperature difference of 10 C from the cooling sting at the center to the container surface, and a gravitional loading of 0.000001 g a computed maximum fluid velocity of about 2.4 x 0.00001 cm/sec was reached after about 250 sec. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗

NASA NTRS1997

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

Applications are described of high-performance computing methods to the numerical simulation of complete jet engines. The methodology focuses on the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field elements. New partitioned analysis procedures to treat this coupled three-component problem were developed. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. The NASA-sponsored ENG10 program was used for the global steady state analysis of the whole engine. This program uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 was developed as well as the capability for the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames.

Felippa, C. A.↗

NASA NTRS2004

Hurricane Forecasting with the High-resolution NASA Finite-volume General Circulation Model

A high-resolution finite-volume General Circulation Model (fvGCM), resulting from a development effort of more than ten years, is now being run operationally at the NASA Goddard Space Flight Center and Ames Research Center. The model is based on a finite-volume dynamical core with terrain-following Lagrangian control-volume discretization and performs efficiently on massive parallel architectures. The computational efficiency allows simulations at a resolution of a quarter of a degree, which is double the resolution currently adopted by most global models in operational weather centers. Such fine global resolution brings us closer to overcoming a fundamental barrier in global atmospheric modeling for both weather and climate, because tropical cyclones and even tropical convective clusters can be more realistically represented. In this work, preliminary results of the fvGCM are shown. Fifteen simulations of four Atlantic tropical cyclones in 2002 and 2004 are chosen because of strong and varied difficulties presented to numerical weather forecasting. It is shown that the fvGCM, run at the resolution of a quarter of a degree, can produce very good forecasts of these tropical systems, adequately resolving problems like erratic track, abrupt recurvature, intense extratropical transition, multiple landfall and reintensification, and interaction among vortices.

Atlas, R.↗

NASA NTRS1993

Modeling and new equipment definition for the vibration isolation box equipment system

Our MSAD-funded research project is to provide numerical modeling support for the VIBES (Vibration Isolation Box Experiment System) which is an IML2 flight experiment being built by the Japanese research team of Dr. H. Azuma of the Japanese National Aerospace Laboratory. During this reporting period, the following have been accomplished: A semi-consistent mass finite element projection algorithm for 2D and 3D Boussinesq flows has been implemented on Sun, HP And Cray Platforms. The algorithm has better phase speed accuracy than similar finite difference or lumped mass finite element algorithms, an attribute which is essential for addressing realistic g-jitter effects as well as convectively-dominated transient systems. The projection algorithm has been benchmarked against solutions generated via the commercial code FIDAP. The algorithm appears to be accurate as well as computationally efficient. Optimization and potential parallelization studies are underway. Our implementation to date has focused on execution of the basic algorithm with at most a concern for vectorization. The initial time-varying gravity Boussinesq flow simulation is being set up. The mesh is being designed and the input file is being generated. Some preliminary 'small mesh' cases will be attempted on our HP9000/735 while our request to MSAD for supercomputing resources is being addressed. The Japanese research team for VIBES was visited, the current set up and status of the physical experiment was obtained and ongoing E-Mail communication link was established.

Sani, Robert L.↗

NASA NTRS1991

Computational study of three dimensional viscous flow through a turbine cascade using a multi-domain spectral technique

The three dimensional viscous flow through a planar turbine cascade is numerically simulated by direct solution of the incompressible Navier-Stokes equations. Flow dependence in the spanwise direction is represented by direct expansion in Chebyshev polynomials, while the discretization on planes parallel to the endwalls is accomplished using the spectral element method. Elemental mapping from the physical to the computational space uses an algebraic mapping technique. A fractional time stepping method that consists of an explicit nonlinear convective step, an implicit pressure correction step, and an implicit viscous step is used to advance the Navier-Stokes equations forward in time. Results computed at moderate Reynolds numbers show a three dimensional endwall flow separation, a midspan separation of the blade suction surface boundary layer, and other three-dimensional features such as the presence of a saddle point flow in the endwall region. In addition, the computed skin friction lines are shown to be orthogonal to the surface vorticity lines, demonstrating the accuracy achievable in the present method.

Renaud, Earl W.↗

NASA NTRS1994

Head-on parallel blade-vortex interaction

An experimental and computational study was carried out to investigate the parallel head-on blade-vortex interaction (BVI) and its noise generation mechanism. A shock tube, with an enlarged test section, was used to generate a compressible starting vortex which interacted with a target airfoil. The dual-pulsed holographic interferometry (DPHI) technique and airfoil surface pressure measurements were employed to obtain quantitative flow data during the BVI. A thin-layer Navier-Stokes code (BV12D), with a high-order upwind-biased scheme and a multizonal grid, was also used to simulate numerically the phenomena occurring in the head-on BVI. The detailed structure of a convecting vortex was studied through independent measurements of density and pressure distributions across the vortex center. Results indicate that, in a strong head-on BVI, the opposite pressure peaks are generated on both sides of the leading edge as the vortex approaches. Then, as soon as the vortex passes by the leading edge, the high-pressure peak suddenly moves toward the low-peak-reducing in magnitude as it moves--simultaneously giving rise to the initial sound wave. In both experiment and computation, it is shown that the viscous effect plays a significant role in head-on BVIs.

Lee, Soogab↗

NASA NTRS1977

Parallel computation of unsteady, three-dimensional, chemically reacting, nonequilibrium flow using a time-split finite-volume method on the Illiac IV

A description is presented of the split finite-volume method which is a viable numerical procedure for performing with the aid of a modern special purpose vector computer numerical simulation studies of complicated flow fields, including chemical reactions, about geometrically complex bodies. Such numerical studies are needed for the development of atmospheric entry vehicles such as the space shuttle. The equations which are approximated are quite general and can be used in studies of combustion, pollution, and other chemically reacting flow phenomena, where convective transport effects dominate the influence of radiative, viscous, and other transport mechanisms. The shock perturbed flow about a shuttle orbiter flying at a large angle of attack during atmospheric entry is illustrated. The method uses a time splitting of the convection differencing operator to achieve efficient data management.

Reinhardt, W. A.↗

NASA NTRS1999

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work is limited to very simple models. In this project, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. The proposed work is based on previously funded NASA work that allowed us to develop a two-dimensional numerical method for boiling flows and to demonstrate the ability of the method to simulate film boiling. While numerical simulations of multi-fluid flows have been advanced in a major way during the last five years, or so, similar capability for flows with phase change are still in their infancy. Although the feasibility of the proposed approach has been demonstrated, it has yet to be extended and applied to fully three-dimensional simulations. Here, a fully three-dimensional, parallel, grid adaptive code will be developed. The numerical method will be used to study nucleate boiling in microgravity, with particular emphasis on two aspects of the problem: 1) Examination of the growth of bubbles at a wall nucleation site and the instabilities of rapidly growing bubbles. Particular emphasis will be put on accurately capturing the thin wall layer left behind as a bubble expands along a wall, on computing instabilities on bubble surfaces as bubbles grow, and on quantifying the effects of both these phenomena on heat transfer; and 2) Examination of the effect of shear flow on bubble growth and heat transfer.

Tryggvason, Gretar↗

NASA NTRS2001

Computation of Coupled Thermal-Fluid Problems in Distributed Memory Environment

The thermal-fluid coupling problems are very important to aerospace and engineering applications. Instead of analyzing heat transfer and fluid flow separately, this study merged two well-accepted engineering solution methods, SINDA for thermal analysis and FDNS for fluid flow simulation, into a unified multi-disciplinary thermal fluid prediction method. A fully conservative patched grid interface algorithm for arbitrary two-dimensional and three-dimensional geometry has been developed. The state-of-the-art parallel computing concept was used to couple SINDA and FDNS for the communication of boundary conditions through PVM (Parallel Virtual Machine) libraries. Therefore, the thermal analysis performed by SINDA and the fluid flow calculated by FDNS are fully coupled to obtain steady state or transient solutions. The natural convection between two thick-walled eccentric tubes was calculated and the predicted results match the experiment data perfectly. A 3-D rocket engine model and a real 3-D SSME geometry were used to test the current model, and the reasonable temperature field was obtained.

Wei, H.↗

NASA NTRS1984

Numerical simulation of the interaction of a vortex with stationary airfoil in transonic flow

A perturbation form of an implicit conservative, noniterative numerical algorithm for the two-dimensional thin layer Navier-Stokes and Euler equations is used to compute the interaction flow-field of a vortex with stationary airfoil. A Lamb-like analytical vortex having a finite core is chosen to interact with a thick (NACA 0012) and a thin (NACA 64A006) airfoil independently in transonic flow. Two different configurations of vortex interaction are studied, viz., (1) when the vortex is fixed at one location in the flowfield, and (2) when the vortex is convecting past the airfoil at freestream velocity. Parallel computations of this interacting flowfield are also done using a version of the Transonic Small Disturbance Code (ATRAN2). A special treatment of the leading edge region for thin airfoils is included in this code. With this, the three methods gave qualitatively similar results for the weaker interactions considered in this study. However, the strongest interactions considered proved to be beyond the capabilities of the small disturbance code. The results also show a far greater influence of the vortex on the airfoil flowfield when the vortex is stationary than when it is convecting with the flow.

Srinivasan, G. R.↗

NASA NTRS1984

Three-dimensional Navier-Stokes calculations of multiple interacting vortex rings

Results from a finite-difference Navier-Stokes code for three-dimensional, unsteady, vortical flows in unbounded domains are presented and analyzed in this paper. The vortical flows presented are representative of vortex rings and other closed vortical tubes or structures in fluid mechanics. Such structures are important elements in fluid flows such as jets, atmospheric turbulence, and the far-field wakes of aircraft, and studies of their interaction may aid in an understanding of complex fluid flows. The paper demonstrates that computational methods can be used as a viable alternative or supplement to experimental techniques for studying the physics of vortex flows. The separate visualization of vortex stretching, convection, and diffusion is presented in this paper for a single elliptical vortex ring.The calculations employ a truncated series expansion technique to simulate the unbounded nature of the fluid flow with a finite computational domain, which is a more accurate technique than the conventional freestream boundary specification. The numerical divergence of the three-dimensional vorticity field is considered as a useful estimate of truncation error, and the use of a kinetic energy decay law as a calculation check is demonstrated. Results from the Navier-Stokes code are presented for the unsteady motion of two and four vortex rings along parallel axes, and the results agree qualitatively with experimental flow visualization.

Chamberlain, J. P.↗

NASA NTRS1988

Computer simulation of the Farley-Buneman instability and anomalous electron heating in the auroral ionosphere

The nonlinear saturation of the Farley-Buneman instability in a collisional plasma is studied by a 2 1/2 dimensional electrostatic particle simulation which includes inelastic and elastic collisions of electrons and elastic collision of ions with neutrals. In the simulation, a uniform convection electric field is applied externally so that the relative velocity between the electrons and ions is greater than the ion sound speed and destabilizes the instability. A nonlinear frequency shift from higher to lower frequencies and diffusion of the wave spectrum in two dimensional wave number space are found. It is found that the dominant mechanism for electron heating is due to an enhanced effective electron collision frequency and hence enhanced resistive heating as suggested by Primdahl (1986) and not due to the heating of electrons by the electric field of the waves parallel to the magnetic field. For the ionospheric conditions discussed by Schlegel and St.-Maurice (1981)an anomalous heating rate of about 4 x 10 to the 7th W/cu cm is found.

Machida, S.↗

NASA NTRS1994

Steepened structures in equatorial spread F. 2: Theory

This paper investigates the properties of a one-dimensional fluid model of plasma convection in the equatorial F region ionosphere. The model equations are similar in form to Burgers equation except for additional higher-order spatial derivatives. Like Burgers equation, solution to the model have the form of propagating, shocklike structures. Numerical simulations of the model closely resemble the steepened structures observed by sounding rocket plasma density probes within equatorial spread F. Simulated denstiy power spectra, like the spectra computed from in situ data, seem to possess power law forms with a break at wavelengths of about 100 m. The precise wavenumber of the spectral break is determined by the ambipolar diffusion coefficient. The model predicts that electric field fluctuations perpendicular to the direction of plasma steepening should be proportional to the plasma density fluctuations. Electric field fluctuations parallel to the steepening will be due primarily to the ambipolar field and have a Boltzmann relationship with density (square of the absolute value of delta E) approximately equal to (K(exp 2))(square of the absolute value of (delta n/n)). At wavelengths less than about 300 m, the ambipolar field should be the dominant component of the total field intensity.

Hysell, D. L↗