Search NASA⌕ Search

SEARCH · Search NASA

Results for “space 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Finite difference grid generation by multivariate blending function interpolation

The General Interpolants Method (GIM) code which solves the multidimensional Navier-Stokes equations for arbitrary geometric domains is described. The geometry module in the GIM code generates two and three dimensional grids over specified flow regimes, establishes boundary condition information and computes finite difference analogs for use in the GIM code numerical solution module. The technique can be classified as an algebraic equation approach. The geometry package uses multivariate blending function interpolation of vector-values functions which define the shapes of the edges and surfaces bounding the flow domain. By employing blending functions which conform to the cardinality conditions the flow domain may be mapped onto a unit square (2-D) or unit cube (3-D), thus producing an intrinsic coordinate system for the region of interest. The intrinsic coordinate system facilitates grid spacing control to allow for optimum distribution of nodes in the flow domain.

Anderson, P. G.↗
NASA NTRS1981

Two-dimensional integrating matrices on rectangular grids

The use of integrating matrices in solving differential equations associated with rotating beam configurations is examined. In vibration problems, by expressing the equations of motion of the beam in matrix notation, utilizing the integrating matrix as an operator, and applying the boundary conditions, the spatial dependence is removed from the governing partial differential equations and the resulting ordinary differential equations can be cast into standard eigenvalue form. Integrating matrices are derived based on two dimensional rectangular grids with arbitrary grid spacings allowed in one direction. The derivation of higher dimensional integrating matrices is the initial step in the generalization of the integrating matrix methodology to vibration and stability problems involving plates and shells.

Lakin, W. D.↗
NASA NTRS1981

Implications of natural sources for the latitudinal gradients of NO/y/ in the unpolluted troposphere

The latitudinal distributions of tropospheric odd nitrogen (NO/y/) due to lightning discharges and stratospheric injection are studied by estimating latitude dependences of the sources and the rainout sink from data obtained in August-September, 1977, and May-June, 1978, during project GAMETAG. A one-dimensional (latitude dependent) model is developed for the annually, zonally, and vertically averaged concentration of trace species, NO(y) being treated as a single species, which decouples the NO(y) calculation from the rest of the model since NO(y) has no chemical sources or sinks in the troposphere. The equation governing the NO(y) distribution is solved using finite differences with 33 uniformly spaced grid points in the micron range of -1 to 1. The distribution due to the stratospheric source has a minimum in the tropics, and there are broad maxima centered in the 25-30 degree region in both hemispheres. The calculated NO(y) due to lightning has a maximum near 15 degrees S, and it decreases to nearly a fourth of the maximum value near the poles. The drop off from the maximum is sharper in the Southern Hemisphere in comparison with the Northern Hemisphere.

Hameed, S.↗
NASA NTRS1980

Doppler radar observation of the evolution of a thunderstorm

Triple-Doppler radar observations of the evolution of the three-dimensional structure of a thunderstorm on May 19, 1978 are analyzed. Continuous data were taken over a long period of the non-severe storm's lifetime as it passed through the radar and the portable automated mesonet network. A fairly low cloud top of 10 km and high reflectivities were observed, and horizontal rotations developed in the middle troposphere, which never reached lower levels. The cyclonic and anticyclonic circulations at mid-levels intensified after the maximum cell height was reached, and a high reflectivity maximum lasted 15-20 min aloft, with the core descending to the surface because of rapid fallout from the largest precipitation particles. Due to the fairly small scale of the updrafts, future use of smaller grid spacings is considered a necessity.

Heymsfield, G. M.↗
NASA NTRS1981

Comparison of the Nimbus-4 BUV ozone data with the Ames two-dimensional model

Predictions by the Ames two-dimensional model of altitude, latitude and seasonal variations in ozone distribution are compared with the first two years of Nimbus 4 backscattered ultraviolet (BUV) measurements in a preliminary attempt at model verification. The ozone observations consist of mixing ratios on the 1-, 2-, 5-, and 10-mbar pressure surfaces zonally and time averaged to obtain seasonal means for 1970 and 1971. The model is based on chemical reaction and photolysis rate constants recommended by the NASA Panel for Data Evaluation (1979), diurnally averaged for latitudes from 80 deg N to 80 deg S and altitudes from 0 to 60 km with 5 deg horizontal and 2.5 km vertical grid spacings. The large altitude, latitude and seasonal variations observed in the data are found to agree well with model predictions. Examination of the sensitivity of the model predictions to various assumed parameters indicates that improvements in agreement may be obtained by variations of the odd chlorine mixing ratio, odd-nitrogen level and transport parameters.

Borucki, W. J.↗
NASA NTRS1982

Numerical simulation of a conducting disk partially covered by an insulator interacting with a plasma

A numerical simulation has been made of the 'pinhole effect' which produces the characteristic 'S-shaped' current-voltage curve. A disk-shaped conducting probe immersed in a plasma is modeled using a particle-in-cell (PIC) code. A probe partially covered by a very thin insulating layer is considered, as well as a probe mounted on an insulating disk. The simulation uses a cylindrical particle mover and allows for a variable number of particles in the system. The simulation space grid uses three different mesh sizes, the coarsest being away from the probe and the finest near the probe, in order to accurately calculate the trajectories of the simulation particles contributing the current to the probe and to the surface charge density on the dielectric. The calculation of the electrostatic potential is done self-consistently using successive over-relaxation (SOR). Backscattering and secondary electron emission are included for the case of positive probe voltage.

Chaky, R. C.↗
NASA NTRS1983

Crustal geomagnetic field - Two-dimensional intermediate-wavelength spatial power spectra

Two-dimensional Fourier spatial power spectra of equivalent magnetization values are presented for a region that includes a large portion of the western United States. The magnetization values were determined by inversion of POGO satellite data, assuming a magnetic crust 40 km thick, and were located on an 11 x 10 array with 300 km grid spacing. The spectra appear to be in good agreement with values of the crustal geomagnetic field spatial power spectra given by McLeod and Coleman (1980) and with the crustal field model given by Serson and Hannaford (1957). The spectra show evidence of noise at low frequencies in the direction along the satellite orbital track (N-S). indicating that for this particular data set additional filtering would probably be desirable. These findings illustrate the value of two-dimensional spatial power spectra both for describing the geomagnetic field statistically and as a guide for diagnosing possible noise sources.

Mcleod, M. G.↗
NASA NTRS1983

Global mesoscale variability from collinear tracks of SEASAT altimeter data

Eight to nine sets of global collinear altimeter data with a cross-track grid spacing of approximately 900 km at the equator and 600 km at midlatitude have been obtained from the last 25 days of the 1978 SEASAT mission. Since the geoid is time-invariant, such observations can reveal sea surface height variability caused by dynamic ocean phenomena. Variations due to deep ocean mesoscale features are presently solved for by eliminating the longer wavelength deviations, yielding sea height profiles with an accuracy of a few cm, and a global variability map constructed from these data which reveal an exceptionally realistic perspective of mesoscale energetics. The dominance of exceedingly small variability over extensive regions of the oceans is noted. In both the Atlantic and the Pacific, the North Equatorial Current systems were clearly expressed as zonal bands of higher variability.

Cheney, R. E.↗
NASA NTRS1983

Three-dimensional viscous-flow computations using a directionally hybrid implicit-explicit procedure

A new, directionally dependent, hybrid numerical algorithm for solving the unsteady, three-dimensional Navier-Stokes equations has been developed and used to compute the viscous supersonic flow over complex configurations, which may generate local regions of embedded subsonic or streamwise separated flows or both. The new hybrid implicit-explicit algorithm is derived from the more general implicit Beam-Warming algorithm and is particularly suitable for viscous computations in which the grid spacing in the direction outward from the body is considerably smaller than the spacing in the other two directions. Numerical results obtained from both the hybrid and implicit schemes are presented and compared on the basis of numerical stability, convergence history, and computer and core memory requirements.

Rizk, Y. M.↗
NASA NTRS1983

The influence of convective activity on the vorticity budget

The influence of convective activity on the vorticity budget was determined during the AVE VII and AVE-SESAME I periods. This was accomplished by evaluating each term in the expanded vorticity equation with twisting and tilting and friction representing the residual of all other terms. Convective areas were delineated by use of radar summary charts. The influence of convective activity was established by analyzing contoured fields of each term as well as numerical values and profiles of the various terms in convective and nonconvective areas. Vertical motion was computed by the kinematic method, and all computations were performed over the central United States using a grid spacing of 158 km. The results show that, in convective areas in particular, the residual is of comparable magnitude to the horizontal advection and divergence terms, and therefore, cannot be neglected. In convective areas, the residual term represents a sink of vorticity below 500 mb and a strong source near 300 mb. In nonconvective areas, the residual was small in magnitude at all levels, but tended to be negative (vorticity sink) at 300 mb. The local change term, over convective areas, tended to be balanced by the residual term, and appeared to be a good indicator of development (vorticity becoming more cyclonic). Finally, the shape of the vertical profiles of the term in the budget equation agreed with those found by other investigators for easterly waves, but the terms were one order of magnitude larger than those for easterly waves.

Townsend, T. L.↗
NASA NTRS1984

Ocean Modeling of the North Atlantic

Present modeling of the North Atlantic is inadequate and can be improved in a number of ways. A number of important physical processes are listed in five categories from the viewpoints of how they are treated in isolation, how they are usually represented in present ocean basin models, and how they may be better represented in future models. In the first two categories of vertical boundary processes and internal vertical mixing, parameterizations exist which can easily be incorporated into models and which will have important effects on the simulated structure of the North Atlantic. For the third catagory (mesoscale eddy effects), adequate parameterizations do not exist; but the order of magnitude of the effects is known from observational and process-model studies. A horizontal grid spacing of 100 km or less in required to allow parameterizations with this order of magnitude, as well as to resolve the time-averaged ocean fields. In the fourth category of large scale transports improvements are suggested by way of increased vertical resolution and by the requirement that lateral mixing due to eddies takes place on isopycnal surfaces. Model incorporation of the latter phenomenta is underway. In the fifth category of miscellaneous high-latitude processes, formulations for the treatment of sea ice are available for use. However, the treatment of gravitational instability, which is crucial to deepwater formation in the Atlantic Ocean, will require additional refinements to account for the unresolved physics of chimney formations in the open ocean and buoyant plumes near ocean boundaries.

Albert J Semtner↗
NASA NTRS1984

Analysis of an infinite array of rectangular microstrip patches with idealized probe feeds

A solution is presented to the problem of an infinite array of microstrip patches fed by idealized current probes. The input reflection coefficient is calculated versus scan angle in an arbitrary scan plane, and the effects of substrate parameters and grid spacing are considered. It is pointed out that even when a Galerkin method is used the impedance matrix is not symmetric due to phasing through a unit cell, as required for scanning. The mechanism by which scan blindness can occur is discussed. Measurement results are presented for the reflection coefficient magnitude variation with angle for E-plane, H-plane, and D-plane scans, for various substrate parameters. Measured results from waveguide simulators are also presented, and the scan blindness phenomenon is observed and discussed in terms of forced surface waves and a modified grating lobe diagram.

Pozar, D. M.↗
NASA NTRS1985

A self-adaptive-grid method with application to airfoil flow

A self-adaptive-grid method is described that is suitable for multidimensional steady and unsteady computations. Based on variational principles, a spring analogy is used to redistribute grid points in an optimal sense to reduce the overall solution error. User-specified parameters, denoting both maximum and minimum permissible grid spacings, are used to define the all-important constants, thereby minimizing the empiricism and making the method self-adaptive. Operator splitting and one-sided controls for orthogonality and smoothness are used to make the method practical, robust, and efficient. Examples are included for both steady and unsteady viscous flow computations about airfoils in two dimensions, as well as for a steady inviscid flow computation and a one-dimensional case. These examples illustrate the precise control the user has with the self-adaptive method and demonstrate a significant improvement in accuracy and quality of the solutions.

Nakahashi, K.↗
NASA NTRS1986

A new approximate LU factorization scheme for the Reynolds-averaged Navier-Stokes equations

A new approximate LU factorization scheme is developed to solve the steady state Reynolds-averaged Navier-Stokes (NS) equations. Central differencing is used for both implicit and explicit operator and special care is taken to obtain well-conditioned factors on the implicit side. The scheme is then analyzed and optimized according to a simple linear analysis. It is unconditionally stable for the model hyperbolic equation in both two- and three-dimensions. However, the requirement for well-conditioned factors has essentially limited the effective time step the scheme can achieve. Supersonic and transonic three-dimensinal flows past a hemisphere cylinder are computed to demonstrate the convergence characteristics of the scheme. A good convergence rate is achieved for the inviscid case. Finally, an explicit eigenvector annihilation procedure is adopted successfully to remove the stiffness caused by the fine grid spacing for viscous flows.

Pan, D.↗
NASA NTRS1985

Design of an advanced flight planning system

The demand for both fuel conservation and four-dimensional traffic management require that the preflight planning process be designed to account for advances in airborne flight management and weather forecasting. The steps and issues in designing such an advanced flight planning system are presented. Focus is placed on the different optimization options for generating the three-dimensional reference path. For the cruise phase, one can use predefined jet routes, direct routes based on a network of evenly spaced grid points, or a network where the grid points are existing navaid locations. Each choice has its own problem in determining an optimum solution. Finding the reference path is further complicated by choice of cruise altitude levels, use of a time-varying weather field, and requiring a fixed time-of-arrival (four-dimensional problem).

Sorensen, J. A.↗
NASA NTRS1987

Leading-edge vortex solutions with large total pressure losses

Computations are presented for a Lambda = 75 deg delta wing in a supersonic freestream under two conditions which lead to leading-edge vortices. For one condition, analysis of the computed vortical flow reveals a closed streamline in the core. From varying computational parameters, it appears that this is due to truncation error of the convective derivatives. For the other condition, comparisons are made with wind-tunnel data, and good agreement is noted for pitot pressure distributions, flow angles on the symmetry plane, and the position of an embedded shock. Many of the aerodynamic parameters are shown to be insensitive to grid spacing.

Murman, Earll M.↗
NASA NTRS1989

An embedded mesh procedure for leading-edge vortex flows

A cell-vertex scheme is outlined for solving the flow about a delta wing with M (sub infinity) is greater than 1. Embedded regions of mesh refinement allow solutions to be obtained which have much higher resolution than those achieved to date. Effects of mesh refinement and artificial viscosity on the solutions are studied, to determine at what point leading-edge vortex solutions are grid-converged. A macroscale and a microscale for the size of the vortex are defined, and it is shown that the macroscale (which includes the wing surface properties) is converged on a moderately refined grid, while the microscale is very sensitive to grid spacing. The level of numerical diffusion in the core of the vortex is found to be substantial. Comparisons with the experiment are made for two cases which have transonic cross-flow velocities.

Powell, Kenneth G.↗
NASA NTRS1988

A 3-dimensional numerical simulation of the atmospheric injection of aerosols by a hypothetical basaltic fissure eruption

Researchers simulated the atmospheric response to a hypothetical basaltic fissure eruption using heating rates based on the Roza flow eruption. The simulation employs the Colorado State University Regional Atmospheric Model (RAMS) with scavenging effects. The numerical model is a three-dimensional non-hydrostatic time-split compressible cloud/mesoscale model. Explicit microphysics include prediction of cloud, rain, crystal, and hail precipitation types. Nucleation and phoretic scavenging are predicted assuming that the pollutant makes an effective cloud droplet nucleus. Smoke is carried as a passive tracer. Long and short wave radiation heating tendencies, including the effects of the smoke, are parameterized. The longwave emission by the lava surface is neglected in the parameterization and included as an explicit heating term instead. A regional scale domain of 100 x 100 km in the horizontal and 22 km high is used. The horizontal grid spacing is taken to be 2 km and the vertical spacing is taken to be 0.75 km. The initial atmospheric state is taken to be horizontally homogenous and based on the standard atmospheric sounding. The fissure is assumed to be 90 km long and oriented in a zig/zag pattern.

Tripoli, Gregory J.↗