Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid following”

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 73 records · Page 4

Tracer water transport and subgrid precipitation variation within atmospheric general circulation models

A capability is developed for monitoring tracer water movement in the three-dimensional Goddard Institute for Space Science Atmospheric General Circulation Model (GCM). A typical experiment with the tracer water model follows water evaporating from selected grid squares and determines where this water first returns to the Earth's surface as precipitation or condensate, thereby providing information on the lateral scales of hydrological transport in the GCM. Through a comparison of model results with observations in nature, inferences can be drawn concerning real world water transport. Tests of the tracer water model include a comparison of simulated and observed vertically-integrated vapor flux fields and simulations of atomic tritium transport from the stratosphere to the oceans. The inter-annual variability of the tracer water model results is also examined.

Koster, Randal D.↗

Isoparametric finite element methods to estimate discontinuous coefficients in two-dimensional elliptic equations

A summary of findings is presented on the problem of estimating discontinuous coefficients, including the locations of interfaces at which discontinuities occur, for elliptic equations on a two-dimensional domain. Finite-element approximations are constructed that are based on an underlying grid that adaptively moves to follow locations of interfaces as they are iteratively updated in the numerical estimation process.

Lamm, Patricia K.↗

Does the resistive tearing instability nonlinearly evolve to a fast reconnection mode

A fundamental problem in applying linear tearing instability theory to the rapid processes (particle acceleration, heating) in flares was the characteristically slow rate of reconnection. This problem can be at least partially overcome if the tearin mode nonlinearly evolves to a regime in which the reconnection rate is substantially enhanced, such as that for the Petschek configuration. This possibility was often suggested, and some numerical simulations appear to provide support for such a view. Numerical simulation are used to study the nonlinear evolution of the tearing stability and show that a fast Petschek-like regime may not be achieved. This conclusion follows when there are sufficient grid points within the diffusion region to completely resolve the nonlinear dynamic interactions in the diffusion layer. When the numerical resolution is not adequate, the solution does appear to approach a Petschek configuration. The resolved solution contains reverse flow vortices and current sheets, terminated with a current reversal, similar to those obtained by Syrovatsky (JEPT, 33, 933, 1971).

Steinolfson, R. S.↗

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

Turbulent combustion can not be simulated adequately by conventional moment closure turbulent models. The probability density function (PDF) method offers an attractive alternative: in a PDF model, the chemical source terms are closed and do not require additional models. Because the number of computational operations grows only linearly in the Monte Carlo scheme, it is chosen over finite differencing schemes. A grid dependent Monte Carlo scheme following J.Y. Chen and W. Kollmann has been studied in the present work. It was found that in order to conserve the mass fractions absolutely, one needs to add further restrictions to the scheme, namely alpha(sub j) + gamma(sub j) = alpha(sub j - 1) + gamma(sub j + 1). A new algorithm was devised that satisfied this restriction in the case of pure diffusion or uniform flow problems. Using examples, it is shown that absolute conservation can be achieved. Although for non-uniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Hsu, Andrew T.↗

TPSAS-NF1676L-33601-DND

The Clouds and the Earth's Radiant Energy System (CERES) project now has over 18 years of an accurately observed top-of-the-atmosphere (TOA) flux record for climate monitoring and diagnostic studies. CERES provides the climate community the following parameters: coincident instantaneous 1? gridded CERES observed TOA fluxes, computed profile and surface fluxes, as well as MODIS cloud and aerosol retrievals. Clouds and radiation interaction one key factor that dominates climate feedbacks and is also the most difficult problem with large uncertainty. To further advance our understanding of the cloud-radiation interaction, the climate community need data with accurate fluxes and their associated cloud properties for both observational and modelling study. The new CERES FluxByCldTyp data product is produced for this purpose. The flux product combines for the first time CERES measured TOA fluxes along with the associated MODIS cloud properties, which links radiative flux directly to a specific cloud type. This was achieved by computing sub-footprint fluxes for the clear-sky and cloudy portions of the CERES footprint. The spatially distributed cloud properties within the CERES footprint were retrieved from 2-km MODIS pixels, that were stratified by cloud type. The sub-footprint fluxes were estimated from MODIS radiances based on empirically derived narrowband to broadband coefficients. The combined sub-footprint fluxes can then be validated with the observed footprint flux. This presentation will focus on the algorithm development of the flux-by-cloud-type product and its validation.

Moguo Sun↗

TPSAS-NF1676L-33753-DND

The Clouds and the Earth's Radiant Energy System (CERES) project now has over 18 years of an accurately observed top-of-the-atmosphere (TOA) flux record for climate monitoring and diagnostic studies. CERES provides the climate community the following parameters: coincident instantaneous 1 gridded CERES observed TOA fluxes, computed profile and surface fluxes, as well as MODIS cloud and aerosol retrievals. Clouds and radiation interaction one key factor that dominates climate feedbacks and is also the most difficult problem with large uncertainty. To further advance our understanding of the cloud-radiation interaction, the climate community need data with accurate fluxes and their associated cloud properties for both observational and modelling study. The new CERES FluxByCldTyp data product is produced for this purpose. The flux product combines for the first time CERES measured TOA fluxes along with the associated MODIS cloud properties, which links radiative flux directly to a specific cloud type. This was achieved by computing sub-footprint fluxes for the clear-sky and cloudy portions of the CERES footprint. The spatially distributed cloud properties within the CERES footprint were retrieved from 2-km MODIS pixels, that were stratified by cloud type. The sub-footprint fluxes were estimated from MODIS radiances based on empirically derived narrowband to broadband coefficients. The combined sub-footprint fluxes can then be validated with the observed footprint flux. This presentation will focus on the algorithm development of the flux-by-cloud-type product and its validation. The paper will also discuss CERES FlxbyCldTyp simulator and its application.

Moguo Sun↗

Radar analyses of mesoscale meteorological phenomena during the AVE/VAS correlation field experiment

Radar data were collected during the selected critical times of the field experiment and subsequently analyzed for specific items. The analyses of the radar data for the AVE/VAS experiment provided statistically significant values for each of the 10 Km by 10 Km grids within radar range. The resulting information for correlation with satellite data included the following derived items that were averaged for each grid area: (1) rainfall rate in mm/hr; (2) dBZ (reflectivity) values; (3) accumulated rainfall values per hour; (4) accumulated rainfall values for a 6-hour period; (5) vertically integrated liquid water content per square meter; and (6) vertical height of the radar axis at the midpoint of each grid. Additional products derived from radar data are being investigated. An example of one such product is the derivation of the errors in integrated rainfall with different sampling periods. This is of significance for correlation with satellite data in that normally a step-function type of rainfall rate is used to derive the total rainfall over a period.

Huebner, G. L.↗

Some Numerical Experiments on Detonation Wave Propagation

In this paper we present the results of a series of numerical experiments done on the propagation and initiation of a detonation wave. The calculations are performed in one-dimension, with considerable grid resolution. Of particular interest are the following questions: (1) the nature of periodic and chaotic instabilities generated by the wave; (2) the influence of the grid resolution on these instabilities; (3) the influence of the 'quality' of the numerical scheme; and (4) the influence of 'noise'. In the calculations, we use a second-order Total Variation Diminishing (TVD) scheme as the basic numerical method, with grid spacings as low as a fraction of a micron. Detonations waves are generated at the closed end of a tube, and allowed to propagate for approximately 20 cm. The required energy for successful initiation of the detonation will be measured for different cases of grid resolution and numerical schemes. A modified version of the TVD scheme has also been devised, which allows for much lower numerical diffusion of the radical species in the exponentially growing region behind the shock. The effect of this modification will be demonstrated. Oscillations in peak pressure and induction length are seen to develop in some cases: the oscillations can go through a sequence of modes, from a regular, high frequency mode to a low frequency mode with period doubling. A chaotic regime can also be obtained. General conclusions on the quality of algorithms will be presented. We will also discuss the performance of a version of the code developed on the IBM SP2 parallel computer.

Cambier, Jean-Luc↗

A Data Parallel Multizone Navier-Stokes Code

We have developed a data parallel multizone compressible Navier-Stokes code on the Connection Machine CM-5. The code is set up for implicit time-stepping on single or multiple structured grids. For multiple grids and geometrically complex problems, we follow the "chimera" approach, where flow data on one zone is interpolated onto another in the region of overlap. We will describe our design philosophy and give some timing results for the current code. The design choices can be summarized as: 1. finite differences on structured grids; 2. implicit time-stepping with either distributed solves or data motion and local solves; 3. sequential stepping through multiple zones with interzone data transfer via a distributed data structure. We have implemented these ideas on the CM-5 using CMF (Connection Machine Fortran), a data parallel language which combines elements of Fortran 90 and certain extensions, and which bears a strong similarity to High Performance Fortran (HPF). One interesting feature is the issue of turbulence modeling, where the architecture of a parallel machine makes the use of an algebraic turbulence model awkward, whereas models based on transport equations are more natural. We will present some performance figures for the code on the CM-5, and consider the issues involved in transitioning the code to HPF for portability to other parallel platforms.

Jespersen, Dennis C.↗

Performance of Titanium Optics on a NASA 30 Cm Ion Thruster

The results of performance tests with two titanium optics sets are presented and compared to those of molybdenum optics. All tests were conducted on a 30 cm ion thruster that was nearly identical to the NASA Solar Electric Propulsion Technology Application Readiness (NSTAR) thruster design. Optics performance tests were conducted over a thruster input power range of 0.5 to 4.6 kW. Optics performance including impingement-limited total voltages, electron backstreaming limits, screen grid ion transparencies, near-field beam current density profiles, beam divergence angles, and beam divergence thrust correction factors were determined throughout this power range. The impingement-limited total voltages for titanium optics were within 10 to 55 V of those for molybdenum optics. Electron backstreaming limit magnitude as a function of peak beam current density for both molybdenum and titanium optics were within a few volts of each other, indicating similar hot grid gaps for these two grid materials during steady-state operation. Beam divergence half-angles at 90 percent of the total beam current and thrust correction factors for both titanium optics sets were within one degree and one percent, respectively, of those for molybdenum optics. When thruster power was increased to 2.3 kW immediately following discharge ignition, the titanium screen grid came into contact with the accelerator grid within five minutes of ignition. Relative to molybdenum, titanium's larger thermal expansion and smaller thermal conductivity likely caused the screen grid to thermally expand more relative to the accelerator grid during startup.

Soulas, George C.↗

A Navier-Strokes Chimera Code on the Connection Machine CM-5: Design and Performance

We have implemented a three-dimensional compressible Navier-Stokes code on the Connection Machine CM-5. The code is set up for implicit time-stepping on single or multiple structured grids. For multiple grids and geometrically complex problems, we follow the 'chimera' approach, where flow data on one zone is interpolated onto another in the region of overlap. We will describe our design philosophy and give some timing results for the current code. A parallel machine like the CM-5 is well-suited for finite-difference methods on structured grids. The regular pattern of connections of a structured mesh maps well onto the architecture of the machine. So the first design choice, finite differences on a structured mesh, is natural. We use centered differences in space, with added artificial dissipation terms. When numerically solving the Navier-Stokes equations, there are liable to be some mesh cells near a solid body that are small in at least one direction. This mesh cell geometry can impose a very severe CFL (Courant-Friedrichs-Lewy) condition on the time step for explicit time-stepping methods. Thus, though explicit time-stepping is well-suited to the architecture of the machine, we have adopted implicit time-stepping. We have further taken the approximate factorization approach. This creates the need to solve large banded linear systems and creates the first possible barrier to an efficient algorithm. To overcome this first possible barrier we have considered two options. The first is just to solve the banded linear systems with data spread over the whole machine, using whatever fast method is available. This option is adequate for solving scalar tridiagonal systems, but for scalar pentadiagonal or block tridiagonal systems it is somewhat slower than desired. The second option is to 'transpose' the flow and geometry variables as part of the time-stepping process: Start with x-lines of data in-processor. Form explicit terms in x, then transpose so y-lines of data are in-processor. Form explicit terms in y, then transpose so z-lines are in processor. Form explicit terms in z, then solve linear systems in the z-direction. Transpose to the y-direction, then solve linear systems in the y-direction. Finally transpose to the x direction and solve linear systems in the x-direction. This strategy avoids inter-processor communication when differencing and solving linear systems, but requires a large amount of communication when doing the transposes. The transpose method is more efficient than the non-transpose strategy when dealing with scalar pentadiagonal or block tridiagonal systems. For handling geometrically complex problems the chimera strategy was adopted. For multiple zone cases we compute on each zone sequentially (using the whole parallel machine), then send the chimera interpolation data to a distributed data structure (array) laid out over the whole machine. This information transfer implies an irregular communication pattern, and is the second possible barrier to an efficient algorithm. We have implemented these ideas on the CM-5 using CMF (Connection Machine Fortran), a data parallel language which combines elements of Fortran 90 and certain extensions, and which bears a strong similarity to High Performance Fortran. We make use of the Connection Machine Scientific Software Library (CMSSL) for the linear solver and array transpose operations.

Jespersen, Dennis C.↗

A conservative approach for flow field calculations on multiple grids

In the computation of flow fields about complex configurations, it is very difficult to construct body-fitted coordinate systems. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the 'multiple grids' or 'zonal grids' approach and its applications are investigated in this study. The method follows the conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-stage Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Some steady state solutions of the Euler equations are presented and discussed.

Kathong, M.↗

A conservative approach for flow field calculations on multiple grids

In the computation of flow fields about complex configurations, it is very difficult to construct body-fitted coordinate systems. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the multiple grids or zonal grids approach and its applications are investigated in this study. The method follows the conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-state Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Some steady state solutions of the Euler equations are presented and discussed.

Kathong, Monchai↗

Multiblock grid generation for jet engine configurations

The goal was to create methods for generating grids with minimal human intervention that are applicable to a wide range of problems and compatible with existing numerical methods and with existing and proposed computers. The following topics that are related to multiblock grid generation are briefly covered in viewgraph form: finding a domain decomposition, dimensioning grids, grid smoothing, manipulating grids and decompositions, and some specializations for jet engine configurations.

Stewart, Mark E. M.↗

Computer graphics for management: An abstract of capabilities and applications of the EIS system

The Executive Information Services (EIS) system, developed as a computer-based, time-sharing tool for making and implementing management decisions, and including computer graphics capabilities, was described. The following resources are available through the EIS languages: centralized corporate/gov't data base, customized and working data bases, report writing, general computational capability, specialized routines, modeling/programming capability, and graphics. Nearly all EIS graphs can be created by a single, on-line instruction. A large number of options are available, such as selection of graphic form, line control, shading, placement on the page, multiple images on a page, control of scaling and labeling, plotting of cum data sets, optical grid lines, and stack charts. The following are examples of areas in which the EIS system may be used: research, estimating services, planning, budgeting, and performance measurement, national computer hook-up negotiations.

Solem, B. J.↗

Recent Developments in Grid Generation and Force Integration Technology for Overset Grids

Recent developments in algorithms and software tools for generating overset grids for complex configurations are described. These include the overset surface grid generation code SURGRD and version 2.0 of the hyperbolic volume grid generation code HYPGEN. The SURGRD code is in beta test mode where the new features include the capability to march over a collection of panel networks, a variety of ways to control the side boundaries and the marching step sizes and distance, a more robust projection scheme and an interpolation option. New features in version 2.0 of HYPGEN include a wider range of boundary condition types. The code also allows the user to specify different marching step sizes and distance for each point on the surface grid. A scheme that takes into account of the overlapped zones on the body surface for the purpose of forces and moments computation is also briefly described, The process involves the following two software modules: MIXSUR - a composite grid generation module to produce a collection of quadrilaterals and triangles on which pressure and viscous stresses are to be integrated, and OVERINT - a forces and moments integration module.

Chan, William M.↗

Progress in Grid Generation: From Chimera to DRAGON Grids

Hybrid grids, composed of structured and unstructured grids, combines the best features of both. The chimera method is a major stepstone toward a hybrid grid from which the present approach is evolved. The chimera grid composes a set of overlapped structured grids which are independently generated and body-fitted, yielding a high quality grid readily accessible for efficient solution schemes. The chimera method has been shown to be efficient to generate a grid about complex geometries and has been demonstrated to deliver accurate aerodynamic prediction of complex flows. While its geometrical flexibility is attractive, interpolation of data in the overlapped regions - which in today's practice in 3D is done in a nonconservative fashion, is not. In the present paper we propose a hybrid grid scheme that maximizes the advantages of the chimera scheme and adapts the strengths of the unstructured grid while at the same time keeps its weaknesses minimal. Like the chimera method, we first divide up the physical domain by a set of structured body-fitted grids which are separately generated and overlaid throughout a complex configuration. To eliminate any pure data manipulation which does not necessarily follow governing equations, we use non-structured grids only to directly replace the region of the arbitrarily overlapped grids. This new adaptation to the chimera thinking is coined the DRAGON grid. The nonstructured grid region sandwiched between the structured grids is limited in size, resulting in only a small increase in memory and computational effort. The DRAGON method has three important advantages: (1) preserving strengths of the chimera grid; (2) eliminating difficulties sometimes encountered in the chimera scheme, such as the orphan points and bad quality of interpolation stencils; and (3) making grid communication in a fully conservative and consistent manner insofar as the governing equations are concerned. To demonstrate its use, the governing equations are discretized using the newly proposed flux scheme, AUSM+, which will be briefly described herein. Numerical tests on representative 2D inviscid flows are given for demonstration. Finally, extension to 3D is underway, only paced by the availability of the 3D unstructured grid generator.

Liou, Meng-Sing↗

Specificity and Transfer in Learning How to Follow Navigation Instructions

We report a series of experiments that use a navigation task in which instructions for navigating in a space displayed as grids on a computer screen are given to subjects who then attempt to follow them by mouse clicking on the grids. The navigation task was broken down into component dimensions (e.g., presentation mode of the instructions, length of the instructions, characteristics of the display, size of the grids, response type). For each task dimension, one condition was used at training and the same or another condition was used at test. Each task dimension was examined in terms of two measures. One measure provided an index of transfer (i.e., better performance at test than at training when test and training involved different conditions), and the other provided an index of specificity (i.e., better performance at test when training and test conditions were the same than when training and test conditions were different). By and large, these two indices were complementary, so there was evidence of either transfer or specificity but not both. For one dimension transfer but no specificity was evident, and for another dimension specificity but no transfer was evident. For the remaining dimensions, however, there was asymmetrical transfer, with transfer evident for some conditions and specificity evident for others. The findings are interpreted within the procedural reinstatement framework. They have practical implications concerning how to optimize training and how much fidelity to the testing situation is necessary when training.

Healy, Alice F.↗