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 631 records · Page 35

CFD code evaluation for internal flow modeling

Research on the computational fluid dynamics (CFD) code evaluation with emphasis on supercomputing in reacting flows is discussed. Advantages of unstructured grids, multigrids, adaptive methods, improved flow solvers, vector processing, parallel processing, and reduction of memory requirements are discussed. As examples, researchers include applications of supercomputing to reacting flow Navier-Stokes equations including shock waves and turbulence and combustion instability problems associated with solid and liquid propellants. Evaluation of codes developed by other organizations are not included. Instead, the basic criteria for accuracy and efficiency have been established, and some applications on rocket combustion have been made. Research toward an ultimate goal, the most accurate and efficient CFD code, is in progress and will continue for years to come.

Chung, T. J.↗

NAS Grid Benchmarks: A Tool for Grid Space Exploration

We present an approach for benchmarking services provided by computational Grids. It is based on the NAS Parallel Benchmarks (NPB) and is called NAS Grid Benchmark (NGB) in this paper. We present NGB as a data flow graph encapsulating an instance of an NPB code in each graph node, which communicates with other nodes by sending/receiving initialization data. These nodes may be mapped to the same or different Grid machines. Like NPB, NGB will specify several different classes (problem sizes). NGB also specifies the generic Grid services sufficient for running the bench-mark. The implementor has the freedom to choose any specific Grid environment. However, we describe a reference implementation in Java, and present some scenarios for using NGB.

Frumkin, Michael↗

Statistical Analysis of CFD Solutions From the Fifth AIAA Drag Prediction Workshop

A graphical framework is used for statistical analysis of the results from an extensive N-version test of a collection of Reynolds-averaged Navier-Stokes computational fluid dynamics codes. The solutions were obtained by code developers and users from North America, Europe, Asia, and South America using a common grid sequence and multiple turbulence models for the June 2012 fifth Drag Prediction Workshop sponsored by the AIAA Applied Aerodynamics Technical Committee. The aerodynamic configuration for this workshop was the Common Research Model subsonic transport wing-body previously used for the 4th Drag Prediction Workshop. This work continues the statistical analysis begun in the earlier workshops and compares the results from the grid convergence study of the most recent workshop with previous workshops.

Statistical analysis↗

Development and application of computational aerothermodynamics flowfield computer codes

Research was performed in the area of computational modeling and application of hypersonic, high-enthalpy, thermo-chemical nonequilibrium flow (Aerothermodynamics) problems. A number of computational fluid dynamic (CFD) codes were developed and applied to simulate high altitude rocket-plume, the Aeroassist Flight Experiment (AFE), hypersonic base flow for planetary probes, the single expansion ramp model (SERN) connected with the National Aerospace Plane, hypersonic drag devices, hypersonic ramp flows, ballistic range models, shock tunnel facility nozzles, transient and steady flows in the shock tunnel facility, arc-jet flows, thermochemical nonequilibrium flows around simple and complex bodies, axisymmetric ionized flows of interest to re-entry, unsteady shock induced combustion phenomena, high enthalpy pulsed facility simulations, and unsteady shock boundary layer interactions in shock tunnels. Computational modeling involved developing appropriate numerical schemes for the flows on interest and developing, applying, and validating appropriate thermochemical processes. As part of improving the accuracy of the numerical predictions, adaptive grid algorithms were explored, and a user-friendly, self-adaptive code (SAGE) was developed. Aerothermodynamic flows of interest included energy transfer due to strong radiation, and a significant level of effort was spent in developing computational codes for calculating radiation and radiation modeling. In addition, computational tools were developed and applied to predict the radiative heat flux and spectra that reach the model surface.

Venkatapathy, Ethiraj↗

A survey of grid generation techniques in computational fluid dynamics

The use of numerically-generated body-conforming coordinate systems has become well-established in computational fluid dynamics, allowing general codes to be written which can treat boundaries of arbitrary shape. There is now a rapidly increasing interest in the dynamic coupling of the grid with the physical solution so that the grid lines continually move to achieve concentration in regions of strong variation of the physical solution and alignment with shocks, flame fronts, etc. Such adaptive coordinate systems have been developed using a variety of techniques for sensing and reacting to the developing regions of concentration and/or alignment. This paper reviews the techniques and applications of these dynamically adaptive systems and other systems.

Thompson, J. F.↗

A method for solving the transonic full-potential equation for general configurations

A method is developed for solving the full-potential equation for two-dimensional and axisymmetric flow which retains the grid and boundary condition simplicity of the transonic small-disturbance codes. The method is based on a finite-volume formulation of the mass conservation equation in a Cartesian coordinate system, and is an extension of the method of Purvis and Burkhalter (1979). This finite-volume approach, combined with the simple boundary treatment, is shown to result in a highly robust method applicable to a wide range of geometries and flow conditions. The accuracy of the method is demonstrated for general geometries in two-dimensional and axisymmetric flows. The use of this method results in significant gains in convergence rate over the vertical-line over-relaxation scheme by incorporating an AF2-type algorithm (Ballhaus et al., 1978). It is suggested that the simplicity of this method shold allow a relatively easy extension to complex geometries in three-dimensional flows, and complex two-dimensional configurations such as multielement airfoils should be amenable to this method.

Wedan, B.↗

Evaluation of aerothermal modeling computer programs

Various computer programs based upon the SIMPLE or SIMPLER algorithm were studied and compared for numerical accuracy, efficiency, and grid dependency. Four two-dimensional and one three-dimensional code originally developed by a number of research groups were considered. In general, the accuracy and computational efficieny of these TEACH type programs were improved by modifying the differencing schemes and their solvers. A brief description of each program is given. Error reduction, spline flux and second upwind differencing programs are covered.

Hsieh, K. C.↗

A fast upwind solver for the Euler equations on three-dimensional unstructured meshes

An upwind scheme is presented for solving the three-dimensional Euler equations on unstructured tetrahedral meshes. Spatial discretization is accomplished by a cell-centered finite-volume formulation using flux-difference splitting. Higher-order differences are formed by a novel cell reconstruction process which results in computational times per cell comparable to those of structured codes. The approach yields highly resolved solutions in regions of smooth flow while avoiding oscillations across shocks without explicit limiting. Solutions are advanced in time by a 3-stage Runge-Kutta time-stepping scheme with convergence accelerated to steady state by local time stepping and implicit residual smoothing. Solutions are presented for a range of configurations in the transonic speed regime to demonstrate code accuracy, speed, and robustness. The results include an assessment of grid sensitivity and convergence acceleration by mesh sequencing.

Frink, Neal T.↗

Results from computational analysis of a mixed compression supersonic inlet

A numerical study was performed to simulate the critical flow through a supersonic inlet. This flowfield has many phenomena such as shock waves, strong viscous effects, turbulent boundary layer development, boundary layer separations and mass flow suction through the walls (bleed). The computational tools used in this study were two full Navier-Stokes (FNS) codes. The supersonic inlet that was analyzed in this study is the Variable Diameter Centerbody (VDC), inlet. This inlet is a candidate concept for the next generation supersonic transport. Application of the code to the inlet geometry involved effort in generating an efficient grid geometry and specifying boundary conditions, particularly in the bleed region and at the outflow boundary. Results for a critical inlet operation compare favorably to Method of Characteristics predictions and experimental data.

Saunders, J. D.↗

Numerical simulations of the convective flame in white dwarfs

A first step toward better understanding of the mechanism driving convective flames in exploding white dwarfs is presented. The propagation of the convective flame is examined using a two-dimensional implicit hydrodynamical code. The large scales of the instability are captured by the grid while the scales that are smaller than the grid resolution are approximated by a mixing-length approximation. It is found that largescale perturbations (of order of the pressure scale height) do grow significantly during the expansion, leading to a very nonspherical burning front. The combustion rate is strongly enhanced (compared to the unperturbed case) during the first second, but later the expansion of the star suppresses the flame speed, leading to only partial incineration of the nuclear fuel. Our results imply that large-scale perturbations by themselves are not enough to explain the mechanism by which convective flames are driven, and a study of the whole spectrum of relevant perturbations is needed. The implications of these preliminary results on future simulations, in the context of current models for Type Ia supernovae, are discussed.

Livne, Eli↗

Simulated rarefied aerodynamics of the Magellan spacecraft during aerobraking

Aerodynamic loads upon the Magellan spacecraft during aerobraking through the atmosphere of Venus are computed at off-design attitudes with a direct simulation Monte Carlo (DSMC) particle method. Simulated rarefied flows at nominal altitudes near 140 km and an entry speed of 8.6 km/s were compared to simulated and analytic free molecular results. Aerodynamic moments, forces, and heating for rarefied entry at all attitudes were 7-10 percent below free molecular results. All moments acted to restore the vehicle to its nominal zero-pitch, zero-yaw attitude. Suggested canting of the solar panels is an innovative configuration to assess gas-surface interaction during aerobraking. The resulting roll torques about the central body-axis as predicted in rarefied flow simulations were nearly twice that predicted for free molecular flow, although differences became less distinct for thermal accommodation coefficients well below unity. Roll torques increased dramatically with reduced accommodation coefficients employed in the simulation. In the DSMC code, periodic free-molecule boundary conditions and a coarse computational grid and body resolution served to minimize the simulation size and cost while retaining solution validity.

Haas, Brian L.↗

Aircraft system aft telescope cavity configuration study for Stratospheric Observatory for Infrared Astronomy (SOFIA), phase 2

The SOFIA Aircraft System (AS) Phase 1 Study was a broad-based study which addressed itself to satisfying technical and programmatic requirements by drawing from existing technology and applying cost-efficient commercial approaches to the aircraft modification. In this SOFIA AS Phase 2 Study, five critical areas of the aircraft were selected for more detailed investigation: forward pressure bulkhead, aft bulkhead, 'free' shell to bulkhead interface, shell cut-out to bulkhead interface, and flooring. The in-depth investigation of these areas upon a finite element model (FEM), with a fine grid model in areas of particular interest, is discussed. The FEM code used is called 'STRAP' and was developed by the engineering firm, Rasmussen and Associates. STRAP is NASTRAN compatible to within 1%. The loads applied to the model were approximated from known 747 envelope conditions. The areas of investigation, and a section through the fuselage is shown. The thrust of this investigation was to develop the design concepts conceived under phase 1 to the point where detailed design could be undertaken with a high level of confidence.

Source record↗

Numerical Simulations of Wing-Body Junction Flows

The goal of the research project is to contribute to the optimized design of fan bypass systems in advanced turbofan engines such as the Advanced Ducted Propulsors (ADP). The immediate objective is to perform numerical simulation of duct-strut interactions to elucidate the loss mechanisms associated with this configuration that is characteristic of ADP. These numerical simulations would complement an experimental study being undertaken at Purdue University. As the first step in the process, a numerical study of wing-body junction flow is being undertaken as it shares a number of characteristics with the duct-strut interaction flow. The presence of the characteristic horseshoe vortex and the associated secondary flow are the salient features that contribute to making this flow a challenge to predict numerically. The simulations will be performed with the NPARC code on the CRAY Y-MP platform at LeRC. The grids for the simulation have been generated using an algebraic mapping technique with a multisurface algorithm.

Krishnamurthy, R.↗

Comparison of Implicit Schemes for the Incompressible Navier-Stokes Equations

For a computational flow simulation tool to be useful in a design environment, it must be very robust and efficient. To develop such a tool for incompressible flow applications, a number of different implicit schemes are compared for several two-dimensional flow problems in the current study. The schemes include Point-Jacobi relaxation, Gauss-Seidel line relaxation, incomplete lower-upper decomposition, and the generalized minimum residual method preconditioned with each of the three other schemes. The efficiency of the schemes is measured in terms of the computing time required to obtain a steady-state solution for the laminar flow over a backward-facing step, the flow over a NACA 4412 airfoil, and the flow over a three-element airfoil using overset grids. The flow solver used in the study is the INS2D code that solves the incompressible Navier-Stokes equations using the method of artificial compressibility and upwind differencing of the convective terms. The results show that the generalized minimum residual method preconditioned with the incomplete lower-upper factorization outperforms all other methods by at least a factor of 2.

Rogers, Stuart E.↗

Research Efforts in Development of NPARC 2D/3D CFD Codes

The objective of the research was to develop a capability in the NPARC computational fluid dynamics (CFD) code to efficiently solve for unsteady airflows with moving geometry and grids. The application of interest was the unsteady flow in a high-speed aircraft inlet operating at the supercritical condition in which a terminal shock resides within the diffuser.

Slater, John W.↗

Application of Chimera Navier-Stokes Code for High Speed Flows

The primary task for this year was performed in support of the "Trailblazer" project. The purpose of the task was to perform an extensive CFD study of the shock boundary-layer interaction between the engine-diverters and the primary body surfaces of the Trailblazer vehicle. Information gathered from this study would be used to determine the effectiveness of the diverters in preventing the boundary-layer coming off of the vehicle forebody from entering the main engines. The PEGSUS code was used to define the "holes" and "boundaries" for each grid. Two sets of CFD calculations were performed.Extensive post-processing of the results was performed.

Ajmani, Kumud↗

Information Environments

The objective of GRC CNIS/IE work is to build a plug-n-play infrastructure that provides the Grand Challenge Applications with a suite of tools for coupling codes together, numerical zooming between fidelity of codes and gaining deployment of these simulations onto the Information Power Grid. The GRC CNIS/IE work will streamline and improve this process by providing tighter integration of various tools through the use of object oriented design of component models and data objects and through the use of CORBA (Common Object Request Broker Architecture).

Follen, Gregory J.↗

Numerical Models of Human Circulatory System under Altered Gravity: Brain Circulation

A computational fluid dynamics (CFD) approach is presented to model the blood flow through the human circulatory system under altered gravity conditions. Models required for CFD simulation relevant to major hemodynamic issues are introduced such as non-Newtonian flow models governed by red blood cells, a model for arterial wall motion due to fluid-wall interactions, a vascular bed model for outflow boundary conditions, and a model for auto-regulation mechanism. The three-dimensional unsteady incompressible Navier-Stokes equations coupled with these models are solved iteratively using the pseudocompressibility method and dual time stepping. Moving wall boundary conditions from the first-order fluid-wall interaction model are used to study the influence of arterial wall distensibility on flow patterns and wall shear stresses during the heart pulse. A vascular bed modeling utilizing the analogy with electric circuits is coupled with an auto-regulation algorithm for multiple outflow boundaries. For the treatment of complex geometry, a chimera overset grid technique is adopted to obtain connectivity between arterial branches. For code validation, computed results are compared with experimental data for steady and unsteady non-Newtonian flows. Good agreement is obtained for both cases. In sin-type Gravity Benchmark Problems, gravity source terms are added to the Navier-Stokes equations to study the effect of gravitational variation on the human circulatory system. This computational approach is then applied to localized blood flows through a realistic carotid bifurcation and two Circle of Willis models, one using an idealized geometry and the other model using an anatomical data set. A three- dimensional anatomical Circle of Willis configuration is reconstructed from human-specific magnetic resonance images using an image segmentation method. The blood flow through these Circle of Willis models is simulated to provide means for studying gravitational effects on the brain circulation under auto-regulation.

Kim, Chang Sung↗