Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Unstructured Grid”

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 667 records · Page 37

Annual Research Briefs - 1996, Center for Turbulence Research

Partial contents include: (1) Simulation of jet impingement heat transfer with the k-e-V2 model; (2) A new methodology for turbulence modelers using Direct Numerical Simulation (DNS) database analysis; (3) A particle representation model for the deformation of homogeneous turbulence; (4) Modeling of inhomogeneous compressible turbulence using a two-scale statistical theory; (5) Lewis number effects on partially premixed flames; (6) Conditional moment closure of mixing and reaction in turbulent non-premixed combustion; (7) Dynamic models for LES of turbulent front propagation with a spectral method; (8) Jet noise using large eddy simulation; (9) Towards numerical simulations of trailing-edge aeroacoustics; (10) A model for the onset of vortex breakdown; (11) Compressing an elliptic vortex: transition to turbulence by tumble breakdown; (12) A new method for the adaptive control of vortex-wall interactions; (13) Symmetries in turbulent boundary layer flows; (14) Near-wall models in large eddy simulations of flow behind a backward-facing step; (15) Assessment of the subgrid-scale models at low and high Reynolds numbers; (16) Large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids; (17) Progress on LES of flow past a circular cylinder; (18) Experimental investigation of flow through an asymmetric plane diffuser; (19) Progress in the large-eddy simulation of an asymmetric plane diffuser; (20) A numerical study of self-similarity in a turbulent plane wake using large-eddy simulation; (21) On the generation of turbulent inflow conditions for boundary layer simulations, and (22) A new class of finite difference schemes.

COMPUTATIONAL FLUID DYNAMICS↗

A Least-Squares Finite Element Method for Electromagnetic Scattering Problems

The least-squares finite element method (LSFEM) is applied to electromagnetic scattering and radar cross section (RCS) calculations. In contrast to most existing numerical approaches, in which divergence-free constraints are omitted, the LSFF-M directly incorporates two divergence equations in the discretization process. The importance of including the divergence equations is demonstrated by showing that otherwise spurious solutions with large divergence occur near the scatterers. The LSFEM is based on unstructured grids and possesses full flexibility in handling complex geometry and local refinement Moreover, the LSFEM does not require any special handling, such as upwinding, staggered grids, artificial dissipation, flux-differencing, etc. Implicit time discretization is used and the scheme is unconditionally stable. By using a matrix-free iterative method, the computational cost and memory requirement for the present scheme is competitive with other approaches. The accuracy of the LSFEM is verified by several benchmark test problems.

Wu, Jie↗

Research in Computational Aeroscience Applications Implemented on Advanced Parallel Computing Systems

Improving the numerical linear algebra routines for use in new Navier-Stokes codes, specifically Tim Barth's unstructured grid code, with spin-offs to TRANAIR is reported. A fast distance calculation routine for Navier-Stokes codes using the new one-equation turbulence models is written. The primary focus of this work was devoted to improving matrix-iterative methods. New algorithms have been developed which activate the full potential of classical Cray-class computers as well as distributed-memory parallel computers.

Wigton, Larry↗

Assessment of an Euler-Interacting Boundary Layer Method Using High Reynolds Number Transonic Flight Data

Flight-measured high Reynolds number turbulent-flow pressure distributions on a transport wing in transonic flow are compared to unstructured-grid calculations to assess the predictive ability of a three-dimensional Euler code (USM3D) coupled to an interacting boundary layer module. The two experimental pressure distributions selected for comparative analysis with the calculations are complex and turbulent but typical of an advanced technology laminar flow wing. An advancing front method (VGRID) was used to generate several tetrahedral grids for each test case. Initial calculations left considerable room for improvement in accuracy. Studies were then made of experimental errors, transition location, viscous effects, nacelle flow modeling, number and placement of spanwise boundary layer stations, and grid resolution. The most significant improvements in the accuracy of the calculations were gained by improvement of the nacelle flow model and by refinement of the computational grid. Final calculations yield results in close agreement with the experiment. Indications are that further grid refinement would produce additional improvement but would require more computer memory than is available. The appendix data compare the experimental attachment line location with calculations for different grid sizes. Good agreement is obtained between the experimental and calculated attachment line locations.

Bonhaus, Daryl L.↗

Tetrahedral Finite-Volume Solutions to the Navier-Stokes Equations on Complex Configurations

A review of the algorithmic features and capabilities of the unstructured-grid flow solver USM3Dns is presented. This code, along with the tetrahedral grid generator, VGRIDns, is being extensively used throughout the U.S. for solving the Euler and Navier-Stokes equations on complex aerodynamic problems. Spatial discretization is accomplished by a tetrahedral cell-centered finite-volume formulation using Roe's upwind flux difference splitting. The fluxes are limited by either a Superbee or MinMod limiter. Solution reconstruction within the tetrahedral cells is accomplished with a simple, but novel, multidimensional analytical formula. Time is advanced by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the near-wall region of the boundary layer. The issues of accuracy and robustness of USM3Dns Navier-Stokes capabilities are addressed for a flat-plate boundary layer, and a full F-16 aircraft with external stores at transonic speed.

Frink, Neal T.↗

On Convergence Acceleration Techniques for Unstructured Meshes

A discussion of convergence acceleration techniques as they relate to computational fluid dynamics problems on unstructured meshes is given. Rather than providing a detailed description of particular methods, the various different building blocks of current solution techniques are discussed and examples of solution strategies using one or several of these ideas are given. Issues relating to unstructured grid CFD problems are given additional consideration, including suitability of algorithms to current hardware trends, memory and cpu tradeoffs, treatment of non-linearities, and the development of efficient strategies for handling anisotropy-induced stiffness. The outlook for future potential improvements is also discussed.

Mavriplis, Dimitri J.↗

A New Axisymmetric MHD Model of the Interaction of the Solar Wind with Venus

A new two-dimensional axisymmetric MHD model is used to study the interaction of the solar wind with Venus under conditions where the interplanetary field is approximately aligned with the solar wind velocity. This numerical model solves the MHD transport equations for density, velocity, pressure, and magnetic field on an adaptively refined, unstructured grid system. This use of an adaptive grid allows high spatial resolution in regions of large density/velocity gradients and yet can be run on a workstation. The actual grid sizes vary from about 0.06 R(sub v) near the bowshock to 2 R(sub v) in the unperturbed solar wind. The results of the calculations are compared with observed magnetic field values obtained from the magnetometer on the Pioneer Venus Orbiter, at a time when the angle between the solar wind velocity vector and the interplanetary magnetic field (IMF) was only 7.6 deg. Good qualitative agreement between the observed and calculated field behavior is found. The overall results suggest that the induced magnetotail disappears when the IMF is radial for an extended time period and implies that it weakens when the field rotated through a near-radial orientation.

DeZeeuw, Darren L.↗

Forced Mixer Nozzle Optimization

Computational fluid dynamic (CFD) and computational acoustic analyses (CAA) were performed for a TFE731-40 compound nozzle, a TFE731-60 mixer nozzle and an Energy Efficient Engine (E(sup 3)) mixer nozzle for comparison with available data. The CFD analyses were performed with a three dimensional, Navier-Stokes solution of the flowfield on an unstructured grid using the RAMPANT program. The CAA analyses were performed with the NASA Glenn MGB program using a structured grid. A successful aerodynamic solution for the TFE731-40 compound nozzle operating statically was obtained, simulating an engine operating on a test stand. Analysis of the CFD results of the TFE731-40 with the MGB program produced predicted sound power levels that agree quite well with the measured data front full-scale static engine tests. Comparison of the predicted sound pressure with the data show good agreement near the jet axis, but the noise levels are overpredicted at angles closer to the inlet. The predicted sound power level for the TFE731-60 did not agree as well with measured static engine data as the TFE731-40. Although a reduction in the predicted noise level due to the mixed flow was observed, the reduction was not as significant as the measured data. The analysis of the V2 mixer from the E(sup 3) study showed that peak temperatures predicted in the mixer exit flowfield were within 5 percent of the values measured by the exit probes. The noise predictions of the V2 mixer nozzle tended to be 3-5 dB higher in peak noise level than the measurements. In addition, the maximum frequency of the noise was also overpredicted. An analysis of the 3 candidate mixer nozzle configurations demonstrated the feasibility of using centerbody lobes and porosity to improve mixing efficiency. A final configuration was designed with a predicted thermal mixing efficiency that was 5 percent higher than the 3 candidate mixers. The results of the MGB noise calculations show that the final design will exceed the design goal of a 3 dB reduction in noise as compared to the baseline TFE731-40.

Sheoran, Yogi↗

Parallel Implementation of the Discontinuous Galerkin Method

This paper describes a parallel implementation of the discontinuous Galerkin method. Discontinuous Galerkin is a spatially compact method that retains its accuracy and robustness on non-smooth unstructured grids and is well suited for time dependent simulations. Several parallelization approaches are studied and evaluated. The most natural and symmetric of the approaches has been implemented in all object-oriented code used to simulate aeroacoustic scattering. The parallel implementation is MPI-based and has been tested on various parallel platforms such as the SGI Origin, IBM SP2, and clusters of SGI and Sun workstations. The scalability results presented for the SGI Origin show slightly superlinear speedup on a fixed-size problem due to cache effects.

Baggag, Abdalkader↗

Parallel Visualization of Large-Scale Aerodynamics Calculations: A Case Study on the Cray T3E

This paper reports the performance of a parallel volume rendering algorithm for visualizing a large-scale, unstructured-grid dataset produced by a three-dimensional aerodynamics simulation. This dataset, containing over 18 million tetrahedra, allows us to extend our performance results to a problem which is more than 30 times larger than the one we examined previously. This high resolution dataset also allows us to see fine, three-dimensional features in the flow field. All our tests were performed on the Silicon Graphics Inc. (SGI)/Cray T3E operated by NASA's Goddard Space Flight Center. Using 511 processors, a rendering rate of almost 9 million tetrahedra/second was achieved with a parallel overhead of 26%.

Ma, Kwan-Liu↗

CFD Predictions of Sonic-Boom Characteristics for Unmodified and Modified SR-71 Configurations

Shaped sonic-boom signatures refer to signatures that look something other than the typical N-waves. Shaped sonic-boom signatures such as "flat-top," "ramp-type," or "hybrid-type" waveforms have been shown to reduce the subjective loudness without requiring reductions in overpressure peaks. The shaping of sonic-boom signatures requires increasing the shock rise time and changes in frequency spectra. So far, a flat-top waveform was shown to be achievable in wind tunnels; however, the influence of long propagation distance and real atmosphere on shaped signatures should be addressed using flight tests. Two different approaches have been proposed for sonic-boom minimization flight tests. The first approach, proposed by Eagle Aerospace, is for a flight test using a modified BQM-34 "FIREBEE" remotely piloted vehicle. The 30-foot long FIREBEE has a steady state flight condition at the Mach number and altitude of interest, and it can be recovered by helicopter from the water. As an alternative approach, a modified SR-71 vehicle has been proposed by the McDonnell Douglas Corporation. Benefits of the SR-71 include its variable geometry supersonic inlets, small cockpit bulge, higher Mach number capabilities, slender design, and longer length (105 foot). The present investigation addresses the sonic-boom analysis for the second vehicle.The objective of the current investigation is to assess the feasibility of a modified SR-71 configuration, with McDonnell Douglas-designed fuselage modifications, intended to produce shaped sonic-boom signatures on the ground. The present study describes the use of a higher-order computational fluid dynamics (CFD) method to predict the sonic-boom characteristics for both unmodified and modified SR-71 configurations. An Euler unstructured grid methodology is used to predict the near-field, three-dimensional pressure patterns generated by both SR-71 models. The computed near-field pressure signatures are extrapolated to specified distances below the aircraft down to impingement on the ground using the code MDBOOM. Comparisons of the near-field pressure signatures with available flight-test data are presented in the current paper.

Fouladi, Kamran↗

Unstructured Navier-Stokes Analysis of Full TCA Configuration

This paper presents an Unstructured Navier-Stokes Analysis of Full TCA (Technology Concept Airplane) Configuration. The topics include: 1) Motivation; 2) Milestone and approach; 3) Overview of the unstructured-grid system; 4) Results on full TCA W/B/N/D/E configuration; 5) Concluding remarks; and 6) Future directions.

Frink, Neal T.↗

Simulated Inlet Unstart and Nacelle/Diverter Effects for the Boeing Reference H Configuration

The Boeing Reference H configuration was tested in the NASA Ames 9x7 Supersonic Wind Tunnel. A simulated unstarted inlet was evaluated as well as the aerodynamic performance of the configuration with and without nacelle and diverter components. These experimental results were compared with computational results from the unstructured grid Euler flow solver AIRPLANE. The comparisons between computational and experimental results were good, and demonstrated that the Euler code is capable of efficiently and accurately predicting the changes in the aerodynamic coefficients associated with inlet unstart and the effects of the nacelle and diverter components.

Cliff, Susan E.↗

Advanced Multigrid Solvers for Fluid Dynamics

The main objective of this project has been to support the development of multigrid techniques in computational fluid dynamics that can achieve "textbook multigrid efficiency" (TME), which is several orders of magnitude faster than current industrial CFD solvers. Toward that goal we have assembled a detailed table which lists every foreseen kind of computational difficulty for achieving it, together with the possible ways for resolving the difficulty, their current state of development, and references. We have developed several codes to test and demonstrate, in the framework of simple model problems, several approaches for overcoming the most important of the listed difficulties that had not been resolved before. In particular, TME has been demonstrated for incompressible flows on one hand, and for near-sonic flows on the other hand. General approaches were advanced for the relaxation of stagnation points and boundary conditions under various situations. Also, new algebraic multigrid techniques were formed for treating unstructured grid formulations. More details on all these are given below.

Brandt, Achi↗

National Combustion Code: Parallel Implementation and Performance

The National Combustion Code (NCC) is being developed by an industry-government team for the design and analysis of combustion systems. CORSAIR-CCD is the current baseline reacting flow solver for NCC. This is a parallel, unstructured grid code which uses a distributed memory, message passing model for its parallel implementation. The focus of the present effort has been to improve the performance of the NCC flow solver to meet combustor designer requirements for model accuracy and analysis turnaround time. Improving the performance of this code contributes significantly to the overall reduction in time and cost of the combustor design cycle. This paper describes the parallel implementation of the NCC flow solver and summarizes its current parallel performance on an SGI Origin 2000. Earlier parallel performance results on an IBM SP-2 are also included. The performance improvements which have enabled a turnaround of less than 15 hours for a 1.3 million element fully reacting combustion simulation are described.

Quealy, A.↗

Performance Modeling of an Experimental Laser Propelled Lightcraft

A computational plasma aerodynamics model is developed to study the performance of an experimental laser propelled lightcraft. The computational methodology is based on a time-accurate, three-dimensional, finite-difference, chemically reacting, unstructured grid, pressure- based formulation. The underlying physics are added and tested systematically using a building-block approach. The physics modeled include non-equilibn'um thermodynamics, non-equilibrium air-plasma finite-rate kinetics, specular ray tracing, laser beam energy absorption and equi refraction by plasma, non-equilibrium plasma radiation, and plasma resonance. A series of transient computations are performed at several laser pulse energy levels and the simulated physics are discussed and compared with those of tests and literature. The predicted coupling coefficients for the lightcraft compared reasonably well with those of tests conducted on a pendulum apparatus.

Wang, Ten-See↗

Parallel Tetrahedral Mesh Adaptation with Dynamic Load Balancing

The ability to dynamically adapt an unstructured grid is a powerful tool for efficiently solving computational problems with evolving physical features. In this paper, we report on our experience parallelizing an edge-based adaptation scheme, called 3D_TAG. using message passing. Results show excellent speedup when a realistic helicopter rotor mesh is randomly refined. However. performance deteriorates when the mesh is refined using a solution-based error indicator since mesh adaptation for practical problems occurs in a localized region., creating a severe load imbalance. To address this problem, we have developed PLUM, a global dynamic load balancing framework for adaptive numerical computations. Even though PLUM primarily balances processor workloads for the solution phase, it reduces the load imbalance problem within mesh adaptation by repartitioning the mesh after targeting edges for refinement but before the actual subdivision. This dramatically improves the performance of parallel 3D_TAG since refinement occurs in a more load balanced fashion. We also present optimal and heuristic algorithms that, when applied to the default mapping of a parallel repartitioner, significantly reduce the data redistribution overhead. Finally, portability is examined by comparing performance on three state-of-the-art parallel machines.

Oliker, Leonid↗

Ordering Unstructured Meshes for Sparse Matrix Computations on Leading Parallel Systems

The ability of computers to solve hitherto intractable problems and simulate complex processes using mathematical models makes them an indispensable part of modern science and engineering. Computer simulations of large-scale realistic applications usually require solving a set of non-linear partial differential equations (PDES) over a finite region. For example, one thrust area in the DOE Grand Challenge projects is to design future accelerators such as the SpaHation Neutron Source (SNS). Our colleagues at SLAC need to model complex RFQ cavities with large aspect ratios. Unstructured grids are currently used to resolve the small features in a large computational domain; dynamic mesh adaptation will be added in the future for additional efficiency. The PDEs for electromagnetics are discretized by the FEM method, which leads to a generalized eigenvalue problem Kx = AMx, where K and M are the stiffness and mass matrices, and are very sparse. In a typical cavity model, the number of degrees of freedom is about one million. For such large eigenproblems, direct solution techniques quickly reach the memory limits. Instead, the most widely-used methods are Krylov subspace methods, such as Lanczos or Jacobi-Davidson. In all the Krylov-based algorithms, sparse matrix-vector multiplication (SPMV) must be performed repeatedly. Therefore, the efficiency of SPMV usually determines the eigensolver speed. SPMV is also one of the most heavily used kernels in large-scale numerical simulations.

Oliker, Leonid↗