Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid codes”

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 505 records · Page 28

Plasma simulation using the massively parallel processor

Two dimensional electrostatic simulation codes using the particle-in-cell model are developed on the Massively Parallel Processor (MPP). The conventional plasma simulation procedure that computes electric fields at particle positions by means of a gridded system is found inefficient on the MPP. The MPP simulation code is thus based on the gridless system in which particles are assigned to processing elements and electric fields are computed directly via Discrete Fourier Transform. Currently, the gridless model on the MPP in two dimensions is about nine times slower that the gridded system on the CRAY X-MP without considering I/O time. However, the gridless system on the MPP can be improved by incorporating a faster I/O between the staging memory and Array Unit and a more efficient procedure for taking floating point sums over processing elements. The initial results suggest that the parallel processors have the potential for performing large scale plasma simulations.

Lin, C. S.↗

CFD Process Automation Using Overset Grids

This talk summarizes three applications of the overset grid method for CFD using some level of automated grid generation, flow solution and post-processing. These applications are 2D high-lift airfoil analysis (INS2D code), turbomachinery applications (ROTOR2/3 codes), and subsonic transport wing/body configurations (OVERFLOW code). These examples provide a forum for discussing the advantages and disadvantages of overset gridding for use in an automated CFD process. The goals and benefits of the automation incorporated in each application will be described, as well as the shortcomings of the approaches.

Buning, Pieter G.↗

Numerical studies of the fluid and optical fields associated with complex cavity flows

Numerical solutions for the flowfield about several cavity configurations have been computed using the Reynolds averaged Navier-Stokes equations. Comparisons between numerical and experimental results are made in two dimensions for free shear layers and a rectangular cavity, and in three dimensions for the transonic aero-window problem of the Stratospheric Observatory for Infrared Astronomy (SOFIA). Results show that dominant acoustic frequencies and magnitudes of the self excited resonant cavity flows compare well with the experiment. In addition, solution sensitivity to artificial dissipation and grid resolution levels are determined. Optical path distortion due to the flow field is modelled geometrically and is found to match the experiment. The fluid field was computed using a diagonalized scheme within an overset mesh framework. An existing code, OVERFLOW, was utilized with the additions of characteristic boundary condition and output routines required for reduction of the unsteady data. The newly developed code is directly applicable to a generalized three dimensional structured grid zone. Details are provided in a paper included in Appendix A.

Atwood, Christopher A.↗

An Improved Neutron Transport Algorithm for HZETRN

Long term human presence in space requires the inclusion of radiation constraints in mission planning and the design of shielding materials, structures, and vehicles. In this paper, the numerical error associated with energy discretization in HZETRN is addressed. An inadequate numerical integration scheme in the transport algorithm is shown to produce large errors in the low energy portion of the neutron and light ion fluence spectra. It is further shown that the errors result from the narrow energy domain of the neutron elastic cross section spectral distributions, and that an extremely fine energy grid is required to resolve the problem under the current formulation. Two numerical methods are developed to provide adequate resolution in the energy domain and more accurately resolve the neutron elastic interactions. Convergence testing is completed by running the code for various environments and shielding materials with various energy grids to ensure stability of the newly implemented method.

Slaba, Tony C.↗

Rotor Airloads Prediction Using Loose Aerodynamic Structural Coupling

This work couples a computational fluid dynamics (CFD) code and rotorcraft computational structural dynamics (CSD) code to calculate helicopter rotor airloads across a range of flight conditions. An iterative loose (weak) coupling methodology is used to couple the CFD and CSD codes on a per revolution, periodic basis. The CFD uses a high fidelity, Navier-Stokes, overset grid methodology with first principles-based wake capturing. Modifications are made to the CFD code for aeroelastic analysis. For a UH-60A Blackhawk helicopter, four challenging level flight conditions are computed: 1) low speed (u = 0.15) with blade-vortex interaction, 2) high speed (u = 0.37) with advancing blade negative lift, 3) high thrust with dynamic stall (u = 0.24), and 4) hover. Results are compared with UH-60A Airloads Program fight test data. Most importantly, for all cases the loose coupling methodology is shown to be stable, convergent, and robust with full coupling of normal force, pitching moment, and chord force. In comparison with flight test data, normal force and pitching moment magnitudes are in good agreement. For the high speed and dynamic stall cases a phase lag in comparison with the data is seen, nonetheless, the shapes of the curves are very good. Overall, the results are noteworthy improvement over lifting line aerodynamics used in rotorcraft comprehensive codes.

Potsdam, Mark↗

NASTRAN maintenance and enhancement experiences

The current capability is described which includes isoparametric elements, optimization of grid point sequencing, and eigenvalue routine. Overlay and coding errors were corrected for cyclic symmetry, transient response, and differential stiffness rigid formats. Error corrections and program enhancements are discussed along with developments scheduled for the current year and a brief description of analyses being performed using the program.

Schmitz, R. P.↗

Some prospects for the future of computational fluid dynamics

For the oncoming decade it is anticipated that new generations of high-capacity and high-speed computers will appear and that some of them will be dedicated to fluid dynamics applications. An attempt is made here to examine the potential for a new generation of numerical techniques that will accompany the computer advances. Possibilities for improvement in solution techniques, grid adaptations, turbulence approximations, language constructions, and general code robustness are considered.

Lomax, H.↗

Calculations Of Transonic Flow About A Wing

Report describes calculations of transonic airflows about wing in wind tunnel. Basic equations of flow used in study are Reynolds-averaged Navier-Stokes equations in strong conservation-law form. Equations of flow incorporated into finite-difference computer code called TNS (Transonic Navier-Stokes). Computational grid generated by solution of partial differential equations yielding smooth meshes conforming to surfaces of wing and wind tunnel.

Holst, Terry L.↗

Numerical and analytical study of fluid dynamic forces in seals and bearings

A numerical model based on a transformed, conservative form of the three dimensional Navier-Stokes equation and an analytical model based on lumped fluid parameters are presented and compared with studies of modeled rotor bearing seal systems. The rotor destabilizing factors are related to the rotative character of the flow field. It is shown that these destabilizing factors can be reduced through a decrease in the fluid average circumferential velocity. However, the rotative character of the flow field is a complex three dimensional system with bifurcated secondary flow patterns that significantly alter the fluid circumferential velocity. By transforming the Navier-Stokes equations to those for a rotating observer and using the numerical code PHOENICS-84 with a nonorthogonal body fitted grid, several numerical experiments were carried out to demonstrate the character of this complex flow field. In general, fluid injection and/or preswirl of the flow field opposing the shaft rotation significantly intensified these secondary recirculation zones and thus reduced the average circumferential velocity; injection or preswirl in the direction of rotation significantly weakened these zones.

Tam, L. T.↗

Fast Euler solver for transonic airfoils. I - Theory. II - Applications

Equations written in terms of generalized Riemann variables are presently integrated by inverting six bidiagonal matrices and two tridiagonal matrices, using an implicit Euler solver that is based on the lambda-formulation. The solution is found on a C-grid whose boundaries are very close to the airfoil. The fast solver is then applied to the computation of several flowfields on a NACA 0012 airfoil at various Mach number and alpha values, yielding results that are primarily concerned with transonic flows. The effects of grid fineness and boundary distances are analyzed; the code is found to be robust and accurate, as well as fast.

Dadone, Andrea↗

Computational Fluid Dynamics Program at NASA Ames Research Center

The Computational Fluid Dynamics (CFD) Program at NASA Ames Research Center is reviewed and discussed. The technical elements of the CFD Program are listed and briefly discussed. These elements include algorithm research, research and pilot code development, scientific visualization, advanced surface representation, volume grid generation, and numerical optimization. Next, the discipline of CFD is briefly discussed and related to other areas of research at NASA Ames including experimental fluid dynamics, computer science research, computational chemistry, and numerical aerodynamic simulation. These areas combine with CFD to form a larger area of research, which might collectively be called computational technology. The ultimate goal of computational technology research at NASA Ames is to increase the physical understanding of the world in which we live, solve problems of national importance, and increase the technical capabilities of the aerospace community. Next, the major programs at NASA Ames that either use CFD technology or perform research in CFD are listed and discussed. Briefly, this list includes turbulent/transition physics and modeling, high-speed real gas flows, interdisciplinary research, turbomachinery demonstration computations, complete aircraft aerodynamics, rotorcraft applications, powered lift flows, high alpha flows, multiple body aerodynamics, and incompressible flow applications. Some of the individual problems actively being worked in each of these areas is listed to help define the breadth or extent of CFD involvement in each of these major programs. State-of-the-art examples of various CFD applications are presented to highlight most of these areas. The main emphasis of this portion of the presentation is on examples which will not otherwise be treated at this conference by the individual presentations. Finally, a list of principal current limitations and expected future directions is given.

Holst, Terry L.↗

Recent advances in Runge-Kutta schemes for solving 3-D Navier-Stokes equations

A thin-layer Navier-Stokes has been developed for solving high Reynolds number, turbulent flows past aircraft components under transonic flow conditions. The computer code has been validated through data comparisons for flow past isolated wings, wing-body configurations, prolate spheroids and wings mounted inside wind-tunnels. The basic code employs an explicit Runge-Kutta time-stepping scheme to obtain steady state solution to the unsteady governing equations. Significant gain in the efficiency of the code has been obtained by implementing a multigrid acceleration technique to achieve steady-state solutions. The improved efficiency of the code has made it feasible to conduct grid-refinement and turbulence model studies in a reasonable amount of computer time. The non-equilibrium turbulence model of Johnson and King has been extended to three-dimensional flows and excellent agreement with pressure data has been obtained for transonic separated flow over a transport type of wing.

Vatsa, Veer N.↗

Studies of the charging of a thin dust layer in a plasma

The present study is designed to extend Wilson's model to include refinements in the actual operational design of the simulation code itself and to include additional physics in order to assess the resulting electric field structure and grain charges. Most of the initial phases of the work have focussed on the fact that the code, though in principle reliable, suffers from severe time constraints, requiring as must as 10-50 hours of CPU time for typical runs on a VAX computer. It is likely that by staging the code into steps in which early steps are carried out crudely and later steps successively fine tuned the total CPU time can be cut significantly. Preliminary results indicate the possibility to cut CPU time to as much as 1/3 or 1/5 the former time. The original Wilson code has been carefully examined, tested and documented. Shortcuts and staging systems have been introduced and tested which should make the code produce more accurate results with a finer grid system than would have been feasible heretofor. Continued testing and implementation of these refinements will be carried out at the author's home institution in the near future. In addition to the effort at reducing CPU time to a manageable amount, work is in progress to include a dust-size distribution instead of assuming all dust grains have the same size.

Peterson, Lennart R.↗

Mapping unstructured grid computations to massively parallel computers

Investigated here is this mapping problem: assign the tasks of a parallel program to the processors of a parallel computer such that the execution time is minimized. First, a taxonomy of objective functions and heuristics used to solve the mapping problem is presented. Next, we develop a highly parallel heuristic mapping algorithm, called Cyclic Pairwise Exchange (CPE), and discuss its place in the taxonomy. CPE uses local pairwise exchanges of processor assignments to iteratively improve an initial mapping. A variety of initial mapping schemes are tested and recursive spectral bipartitioning (RSB) followed by CPE is shown to result in the best mappings. For the test cases studied here, problems arising in computational fluid dynamics and structural mechanics on unstructured triangular and tetrahedral meshes, RSB and CPE outperform methods based on simulated annealing. Much less time is required to do the mapping and the results obtained are better. Compared with random and naive mappings, RSB and CPE reduce the communication time two fold for the test problems used. Finally, we use CPE in two applications on a CM-2. The first application is a data parallel mesh-vertex upwind finite volume scheme for solving the Euler equations on 2-D triangular unstructured meshes. CPE is used to map grid points to processors. The performance of this code is compared with a similar code on a Cray-YMP and an Intel iPSC/860. The second application is parallel sparse matrix-vector multiplication used in the iterative solution of large sparse linear systems of equations. We map rows of the matrix to processors and use an inner-product based matrix-vector multiplication. We demonstrate that this method is an order of magnitude faster than methods based on scan operations for our test cases.

Hammond, Steven Warren↗

Multidisciplinary modeling and simulation of a generic hypersonic vehicle

The modeling and simulation of hypersonic vehicles which involve the interaction of a number of major technical disciplines including aerodynamics, structures, heat transfer, propulsion, and control engineering are described. Finite element numerical formulations of individual disciplines including their appropriate integration for multidisciplinary simulation are presented. An efficient unstructured grid generation strategy, evolved by modifying an existing code, is developed.

Gupta, K. K.↗

Computational/experimental analysis of three low sonic boom configurations with design modifications

The Euler code, designated AIRPLANE, which uses an unstructured tetrahedral mesh was used to compute near-field sonic boom pressure signatures on three modern low sonic boom configurations: the Mach 2, Mach 3, and Haglund models. The TEAM code which uses a multi-zoned structured grid was used to calculate pressure signatures for the Mach 2 model. The computational pressure signatures for the Mach 2 and Mach 3 models are compared with recent experimental data. The computed pressure signatures were extracted at distances less than one body length below the configuration and extrapolated to the experimental distance. The Mach 2 model was found to have larger overpressures off-ground-track than on-ground-track in both computational and experimental results. The correlations with the experiment were acceptable where the signatures were not contaminated by instrumentation and model-support hardware. AIRPLANE was used to study selected modifications to improve the overpressures of the Mach 2 model.

Cliff, Susan E.↗

Status of the CNRS-LCSR program on high pressure droplet vaporization and burning

Depending on the surrounding flow and thermodynamic conditions, a single droplet may experience several gasification regimes, ranging from the envelope flame regime to pure vaporization. In practical situations, such as rocket propulsion or diesel combustion, the size distribution of droplets is, at best, bimodal, so that the possibility exists for the simultaneous presence of various regimes. For example, very small droplets are transported by the gas phase with zero relative velocity. This picture validates then the spherical symmetry hypothesis applied to the droplet and to the diffusion flame enveloping it. On the other hand, for larger droplets, a relative velocity exists due to drag forces. The most important influence of forced convection on droplet burning is the possibility to extinguish globally the envelope flame, or to establish a flame stabilized in the wake region. The burning rates of these regimes differ strongly. The characteristic time of droplet gasification is also influenced by the surrounding pressure and temperature. A parametric investigation of single droplet burning regimes is then helpful in providing the necessary physical ideas for sub-grid models used in spray combustion numerical prediction codes. The CNRS-LCSR experimental program on droplet vaporization and burning deals with these various regimes: stagnant and convective monocomponent droplet burning convective mono and bicomponent droplet vaporization; high temperature convective mono and biocomponent droplet vaporization; burning regimes of hydrazine and hydroxyl-ammonium-nitrate based monopropellant droplets and the vaporization regimes of liquid oxygen droplets. Studies on interacting droplets and on liquid aluminum droplets will start in the near future. The influence of high pressure is a common feature of all these studies. This paper summarizes the status of the CNRS-LCSR program on the effects of high pressure on monocomponent single droplet burning and vaporization, and some recent results obtained under normal and reduced gravity conditions with suspended droplets are presented. In the work described here, parabolic flights of an aircraft is used to create a reduced gravity environment of the order of 10(exp -2) g.

Chauveau, Christian↗

High pressure droplet burning experiments in reduced gravity

A parametric investigation of single droplet gasification regimes is helpful in providing the necessary physical ideas for sub-grid models used in spray combustion numerical prediction codes. A research program has been initiated at the LCSR to explore the vaporization regimes of single and interacting hydrocarbon and liquid oxygen droplets under high pressure conditions. This paper summarizes the status of the LCSR program on the high pressure burning of single fuel droplets; recent results obtained under normal and reduced gravity conditions with suspended droplets are presented. In the work described here, parabolic flights of the CNES Caravelle is used to create a reduced gravity environment of the order of 10(exp -2) g(sub O). For all the droplet burning experiments reported here, the suspended droplet initial diameters are scattered around 1.5 mm; and the ambient air temperature is 300 K. The ambient pressure is varied between 0.1 MPa and 12 MPa. Four fuels are investigated: methanol (Pc = 7.9 MPa), n-heptane (Pc = 2.74 MPa), n-hexane (Pc = 3.01 MPa) and n-octane (Pc = 2.48 MPa).

Chauveau, Christian↗