Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid simulation”

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 271 records · Page 15

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗

Simulating Transonic Buffet Aerodynamics for the Boeing Transonic Truss-Braced Wing Aircraft

The flow past the Boeing Transonic Truss-Braced Wing (TTBW) aircraft is simulated with a Hybrid Reynolds-Averaged Navier Stokes Large-Eddy Simulations (HRLES) turbulence modeling approach to study transonic buffet onset and evaluate the predictive capability of the numerical approach. All simulations used structured overset curvilinear grids for the free-air configuration, with all simulations performed using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework. Aerodynamic loads and surface pressure obtained from the HRLES solutions are compared to results from unsteady RANS (URANS) and experiments conducted at NASA Ames 11-by 11-Foot Transonic Wind Tunnel. Overall good agreement is obtained in the predicted loads and surface pressure. Unsteady pressure data from HRLES and URANS simulations is utilized to compute power spectral density (PSD) to predict transonic buffet onset and compared to the Kulite data obtained from the experiment. Comparisons of the PSD spectrum also show reasonable agreement with the wind tunnel data, with HRLES and URANS predicting buffet onset at a slightly earlier angle-of-attack than the experiment.

AATT↗

Simulation Of Flow In A Turbine Cascade

Report discusses application of lower-upper (LU) implicit numerical-integration scheme, together with interactive grid-generation method, to simulation of inviscid, compressible flow in two-dimensional model of turbine cascade. Study part of continuing effort to develop techniques for computation of complicated three-dimensional flows in turbomachinery.

Choo, Yung K.↗

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante↗

Analysis of subgrid-scale eddy viscosity with use of results from direct numerical simulations

Without resort to any modeling, subgrid-scale eddy viscosity is computed from the results of high-resolution (64 cubed and 128 cubed grid points) direct numerical simulations of three-dimensional homogeneous isotropic decaying turbulence. In these simulations the eddy viscosity peaks sharply at the cutoff wavenumber, in rough agreement with the results of Kraichnan (1976). In addition, in the low-wavenumber range the eddy viscosity may be negative, contrary to the generally accepted concept of subgrid-scale eddy viscosity. Some possible explanations of this behavior are discussed.

Domaradzki, J. Andrzej↗

High-Fidelity Simulation of Turbulent Flow Past a Gaussian Bump

A spanwise-periodic computation of a turbulent flow past a Gaussian bump is performed in the form of a hybrid direct numerical simulation and wall-resolved large-eddy simulation. A fourth-order spatially-accurate flow solver is employed to perform the simulation, using 10.2 billion grid points for a Reynolds number of 170000 based on the bump height. The key findings from the simulation are reported in the acceleration and deceleration flow regions associated with the bump shape. Significant anisotropy in the normal Reynolds stresses, along both the wall-normal and streamwise directions, is observed within the acceleration region. The ratio between the Reynolds shear stress and turbulent kinetic energy in that region also experiences significant deviations from the norms of a zero pressure gradient turbulent boundary layer. The chosen Reynolds number generates strong flow separation in the adverse pressure gradient region, which is in contrast with a previous simulation at half the Reynolds number that only indicated incipient separation. An internal layer generated in the acceleration region evolves into a free shear layer that develops in the deceleration region and separates. Proper modeling of this inner layer appears crucial to predict the flow separation. Surface curvature effects on the attached flow development are also discussed.

Turbulent Boundary Layer↗

Similar Estimates of Temperature Impacts on Global Wheat Yield by Three Independent Methods

The potential impact of global temperature change on global crop yield has recently been assessed with different methods. Here we show that grid-based and point-based simulations and statistical regressions (from historic records), without deliberate adaptation or CO2 fertilization effects, produce similar estimates of temperature impact on wheat yields at global and national scales. With a 1 C global temperature increase, global wheat yield is projected to decline between 4.1% and 6.4%. Projected relative temperature impacts from different methods were similar for major wheat-producing countries China, India, USA and France, but less so for Russia. Point-based and grid-based simulations, and to some extent the statistical regressions, were consistent in projecting that warmer regions are likely to suffer more yield loss with increasing temperature than cooler regions. By forming a multi-method ensemble, it was possible to quantify 'method uncertainty' in addition to model uncertainty. This significantly improves confidence in estimates of climate impacts on global food security.

Climate impacts↗

Computation Of Flow About A Helicopter Rotor

Improved method developed to simulate numerically flow about helicopter rotor with multiple blades. Computes vortical wake beneath rotor without using simplified ad hoc mathematical models of effects in wake. Patched grids used in numerical simulation of flow about rotor blades. Cylindrical configuration simplifies interpolation at interface between rotating and stationary grids.

Chen, C. L.↗

Grid Erosion Modeling of the NEXT Ion Thruster Optics

Results from several different computational studies of the NEXT ion thruster optics are presented. A study of the effect of beam voltage on accelerator grid aperture wall erosion shows a non-monotonic, complex behavior. Comparison to experimental performance data indicates improvements in simulation of the accelerator grid current, as well as very good agreement with other quantities. Also examined is the effect of ion optics choice on the thruster life, showing that TAG optics provide better margin against electron backstreaming than NSTAR optics. The model is used to predict the change in performance with increasing accelerator grid voltage, showing that although the current collected on the accel grid downstream face increases, the erosion rate decreases. A study is presented for varying doubly-ionized Xenon current fraction. The results show that performance data is not extremely sensitive to the current fraction.

Ernhoff, Jerold W.↗

Simulation of an Isolated Tiltrotor in Hover with an Unstructured Overset-Grid RANS Solver

An unstructured overset-grid Reynolds Averaged Navier-Stokes (RANS) solver, FUN3D, is used to simulate an isolated tiltrotor in hover. An overview of the computational method is presented as well as the details of the overset-grid systems. Steady-state computations within a noninertial reference frame define the performance trends of the rotor across a range of the experimental collective settings. Results are presented to show the effects of off-body grid refinement and blade grid refinement. The computed performance and blade loading trends show good agreement with experimental results and previously published structured overset-grid computations. Off-body flow features indicate a significant improvement in the resolution of the first perpendicular blade vortex interaction with background grid refinement across the collective range. Considering experimental data uncertainty and effects of transition, the prediction of figure of merit on the baseline and refined grid is reasonable at the higher collective range- within 3 percent of the measured values. At the lower collective settings, the computed figure of merit is approximately 6 percent lower than the experimental data. A comparison of steady and unsteady results show that with temporal refinement, the dynamic results closely match the steady-state noninertial results which gives confidence in the accuracy of the dynamic overset-grid approach.

Lee-Rausch, Elizabeth M.↗

Navier-Stokes Simulation Of Turbine Rotor/Stator Interaction

Patched-grid method used to effect three-dimensional Navier-Stokes numerical simulation of flow interacting with both rotor and stator in turbine. Involves use of two or more grids or sets of grids to solve finite-difference versions of Navier-Stokes or other equations of fluid flow. Ability to perform such calculations contributes to understanding of complicated turbomachinery flows and to consequent improvements in designs of turbomachinery.

Rai, Man Mohan↗

Effects of Satellite Sampling Configurations on Derived Gridded Fields

Various configurations of a scanning satellite instrument are simulated by sampling realistic fields of nitrous oxide. Synoptic grids are computed from the resulting simulated orbital data and compared to the original sampled data fields. Results are compared with those obtained by flying a simulated satellite over low-resolution fields and fields that are static in time. Although increasing the number of instrument scan positions does provide more information along an orbital swath, using more than three to five scan positions does not significantly increase the accuracy of global synoptic grids using the gridding techniques described here.

Lait, Leslie R.↗

Cause and Cure - Deterioration in Accuracy of CFD Simulations With Use of High-Aspect-Ratio Triangular Tetrahedral Grids

Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD re-searchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions and also cause numerical instability. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where triangular/tetrahedral elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identities the reason behind the difficulties in use of such high-aspect ratio triangular/tetrahedral elements is presented here. As will be shown, it turns out that the degree of accuracy deterioration of gradient computation involving a triangular element is hinged on the value of its shape factor Gamma def = sq sin Alpha1 + sq sin Alpha2 + sq sin Alpha3, where Alpha1; Alpha2 and Alpha3 are the internal angles of the element. In fact, it is shown that the degree of accuracy deterioration increases monotonically as the value of Gamma decreases monotonically from its maximal value 9/4 (attained by an equilateral triangle only) to a value much less than 1 (associated with a highly obtuse triangle). By taking advantage of the fact that a high-aspect ratio triangle is not necessarily highly obtuse, and in fact it can have a shape factor whose value is close to the maximal value 9/4, a potential solution to avoid accuracy deterioration of gradient computation associated with a high-aspect ratio triangular grid is given. Also a brief discussion on the extension of the current mathematical framework to the tetrahedral-grid case along with some of the practical results of this extension is also provided. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, the effectiveness of the gradient evaluation procedures within the CESE framework (that have their basis on the analysis presented here) to produce accurate and stable results on such high-aspect ratio meshes is also showcased.

CESE Method↗

Regional Climate Simulation of the Anomalous Events of 1998 using a Stretched-Grid GCM with Multiple Areas of Interest

The GEOS (Goddard Earth Observing System) stretched-grid (SG) GCM developed and thoroughly tested over the last few years, is used for simulating the major anomalous regional climate events of 1998. The anomalous regional climate events are simulated simultaneously during the 13 months long (November-1997 - December-1998) SG-GCM simulation due to using the new SG-design with multiple (four) areas of interest. The following areas/regions of interest (one at each global quadrant) are implemented: U.S./Northern Mexico, the El-Nino/Brazil area, India-China, and Eastern Indian Ocean/Australia.

Fox-Rabinovitz, M. S.↗

Managing Complexity in Multidisciplinary Visualization

As high performance computing technology progresses, computational simulations are becoming more advanced in their capabilities. In the computational aerosciences domain, single discipline steady-state simulations computed on a single grid are far from the state-of-the-art. In their place are complex, time-dependent multidisciplinary simulations that attempt to model a given geometry more realistically. The product of these multidisciplinary simulations is a massive amount of data stored in different formats, grid topologies, units of measure, etc., as a result of the differences in the simulated physical domains. In addition to the challenges posed by setting up and performing the simulation, additional challenges exist in analyzing computational results. Visualization plays an important role in the advancement of multidisciplinary simulations. To date, visualization has been used to aid in the interpretation of large amounts of simulation data. Because the human visual system is effective in digesting a large amount of information presented graphically, visualization has helped simulation scientists to understand complex simulation results. As these simulations become even more complex, integrating several different physical domains, visualization will be critical to digest the massive amount of information. Another important role for visualization is to provide a common communication medium from which the domain scientists can use to develop, debug, and analyze their work. Multidisciplinary analyses are the next step in simulation technology, not only in computational aerosciences, but in many other areas such as global climate modeling. Visualization researchers must understand and work towards the challenges posed by multidisciplinary simulation scenarios. This paper addresses some of these challenges, describing technologies that must be investigated to create a useful visualization analysis tool for domain scientists.

Miceli, Kristina D.↗

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.↗

Multi-Fidelity Computational Analysis of a Quiet Single-Main Rotor Helicopter for Air Taxi Operations

A multi-fidelity computational study has been conducted on NASA’s quiet single-main rotor helicopter concept for urban air mobility using computational fluid dynamics and rotorcraft comprehensive analysis tools. Two approaches have been followed in order to simulate the rotor flows: (1) a rotor disk model, which models the rotors as source terms, and (2) a high-fidelity approach using high-order accurate schemes and dual-time stepping, which simulates the rotor with its individual rotating blade grids. The delayed detached-eddy simulation model has been employed. A loose-coupling approach between the flow solver and a rotorcraft comprehensive code is utilized to include vehicle trim and motion variables. The quiet single-main rotor helicopter is one of the conceptual designs intended to focus and guide NASA’s research activities in support of aircraft development for vertical take-off and landing air taxi operations.

RVLT↗

High-Fidelity Computational Analysis Of Nasa’s Air Taxi Concepts

High-fidelity Computational Fluid Dynamics (CFD) simulations have been carried out for NASA’s air taxi concepts. The high-fidelity approach uses high-order accurate schemes and dual-time stepping, which simulates the rotor with its individual rotating blade grids. The delayed detached-eddy simulation model has been employed. A loose-coupling approach between the flow solver and a rotorcraft comprehensive code is utilized to include vehicle trim and motion variables.NASA’s air taxi conceptual designs are intended to focus and guide NASA’s research activities in support of aircraft development for vertical take-off and landing air taxi operations

ARMD↗