Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gridded data”

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 973 records · Page 54

Minimal upper mantle temperature variations consistent with observed heat flow and plate velocities

Heat-flow and plate-velocity measurements are used to model upper-mantle temperature variations; the results are intended to form the basis for a gravity-field-variation model and are also applicable to the interpretation of petrological and seismological data. A 5-deg grid is used, a depth of 280 km is chosen as the fully convecting level, and it is assumed that 85 percent of the global heat production (4.0 x 10 to the 13th W) comes from below this level. The velocity field is calculated by integrating the momentum equations down to 280 km (assuming laterally homogeneous viscosity and density) and then used to determine the temperature fields. The results are presented in graphs, tables, and maps. The largest lateral temperature differences are about 1500 C and occur in the top 20 km, where the largest (about 200 C/Myr) nonlinear terms of the energy equation are also found. Below 50 km, heat transfer becomes more convective than conductive, and the most significant temperature variation appears in the form of negative 'tongues' as cold as -825 C at about 100 km. Temperature variations of at least + or - 180 C are calculated at the fully convecting level.

Kaula, W. M.↗

Simulation of the earth's monthly average regional radiation balance derived from satellite measurements

Computer simulations of satellite-derived earth radiation parameters are used to determine the source and size of errors arising from measuring parameters over 1 month on a 2.5 deg x 2.5 deg longitude-latitude grid. November 1979 GOES data are used as a source of radiation parameter fields within each region. A mathematical model is given for the data-processing algorithms that are used to produce daily, monthly, and monthly hourly estimates of shortwave, longwave and net radiant exitance. Resulting retrieval errors are found to be due to insufficient knowledge of the temporal distribution of measurements. These errors, in combination with similar sampling errors produce biases in monthly net radiant exitance which are complex, regionally-dependent functions of the local time of the measurements. The 40 region standard errors of estimate for three satellites ranged from 2 to 3 W per sq. m., while single satellite errors ranged from 5 to 19 W per sq. m. A three-satellite system is therefore recommended to reduce biases.

Brooks, D. R.↗

An Eigensystem Realization Algorithm (ERA) for modal parameter identification and model reduction

A method, called the Eigensystem Realization Algorithm (ERA), is developed for modal parameter identification and model reduction of dynamic systems from test data. A new approach is introduced in conjunction with the singular value decomposition technique to derive the basic formulation of minimum order realization which is an extended version of the Ho-Kalman algorithm. The basic formulation is then transformed into modal space for modal parameter identification. Two accuracy indicators are developed to quantitatively identify the system modes and noise modes. For illustration of the algorithm, examples are shown using simulation data and experimental data for a rectangular grid structure.

Juang, J. N.↗

Supersonic flow computations using a rectangular-coordinate finite-volume method

A numerical procedure has been developed for the computation of supersonic flows over complex conical geometries. The full potential equation is solved using a finite-volume method with a non-body-fitted rectangular grid. The only mapping done is the transformation of the spherical cross-flow plane to a flat surface using a stereographic projection. A new procedure for very thin fins is described which does not require the resolution of the fin thickness. Applications for simple cones, conical wing-bodies, wave riders and finned geometries compare favorably with existing solutions with body-fitted grids and available experimental data.

Grossman, B.↗

Spectral energetics of the observed and simulated Northern Hemisphere general circulation during blocking episodes

The spectral energetics of the Northern Hemisphere circulation during blocking episodes of the FGGE year is investigated with gridded analyses of observational data and parallel simulation experiments. The purpose of this study is to describe the energetics distinctions of the observed and simulated blockings in the context of the general circulation and to assess the capability of the model to simulate blockings. In the observed circulation a pronounced winter blocking is developed and maintained by the nonlinear wave-wave interaction L(1) from the kinetic energy source for n = 3-10, where L(n) is the transfer of eddy kinetic energy from all other wavenumbers to wavenumber n. In the case of the double blocking in the winter, both L(1) and L(2) support the blocking. The kinetic energy source of n = 10 for upscale input at n = 1 and 2 is supported by the baroclinic conversion at n = 3-10. The simulated winter circulation shows strong baroclinic conversion at all wavenumbers, including ultralong waves. However, the simulation fails to produce pronounced blocking for the absence of L(1), and the converted energy cascades down to shorter waves. The wave-mean transfer of kinetic energy from the large-scale disturbances to the zonal mean component further prevents the accumulation of the kinetic energy at the ultralong waves. In contrast to the winter situation, the summer blocking seems to be directly supported by both L(4) and baroclinic conversion at other planetary-scale waves. Consequently, the summer circulation is better simulated than the winter circulation.

Kung, Ernest C.↗

Velocity-pressure integrated versus penalty finite element methods for high Reynolds number flows

Velocity-pressure integrated and consistent penalty finite element computations of high Reynolds number, laminar flows are presented. In both of the methods, the pressure has been interpolated using linear shape functions for a triangular element. The triangular element is contained inside the bi-quadratic isoparametric element. It has been reported previously that the pressure interpolation method, when used in the velocity-pressure integrated method, yielded accurate computational results for high Reynolds number flows. It is shown that use of the same pressure interpolation method in the consistent penalty finite element method yielded accurate velocity and pressure fields which were comparable to those obtained using the velocity-pressure integrated method. Accuracy of the two finite element methods has been demonstrated by comparing the computational results with available experimental data and/or fine-grid finite difference computational results. Advantages and disadvantages of the two methods are discussed on the basis of accuracy and convergence nature. Example problems considered include a lid-driven cavity flow for Reynolds number of 10,000, a laminar backward-facing step flow, a laminar flow through a nest of cylinders, and a channel flow with an internal blockage. A finite element computer program (NSFLOW/P) for the 2-D, incompressible Navier-Stokes equations is also presented.

Kim, Sang-Wook↗

Low frequency oscillations in total ozone measurements

Low frequency oscillations with periods of approximately one to two months are found in eight years of global grids of total ozone data from the Total Ozone Mapping Spectrometer (TOMS) satellite instrument. The low frequency oscillations corroborate earlier analyses based on four years of data. In addition, both annual and seasonal one-point correlation maps based on the 8-year TOMS data are presented. The results clearly show a standing dipole in ozone perturbations, oscillating with 35 to 50 day periods over the equatorial Indian Ocean-west Pacific region. This contrasts with the eastward moving dipole reported in other data sets. The standing ozone dipole appears to be a dynamical feature associated with vertical atmospheric motions. Consistent with prior analyses based on lower stratospheric temperature fields, large-scale standing patterns are also found in the extratropics of both hemispheres, correlated with ozone fluctuations over the equatorial west Pacific. In the Northern Hemisphere, a standing pattern is observed extending from the tropical Indian Ocean to the north Pacific, across North America, and down to the equatorial Atlantic Ocean region. This feature is most pronounced in the NH summer.

Gao, X. H.↗

Low-frequency oscillations in total ozone measurements

Low frequency oscillations with periods of approximately one to two months are found in eight years of global grids of total ozone data from the Total Ozone Mapping Spectrometer (TOMS) satellite instrument. The low frequency oscillations corroborate earlier analyses based on four years of data. In addition, both annual and seasonal one-point correlation maps based on the 8-year TOMS data are presented. The results clearly show a standing dipole in ozone perturbations, oscillating with 35 to 50 day periods over the equatorial Indian Ocean-west Pacific region. This contrasts with the eastward moving dipole reported in other data sets. The standing ozone dipole appears to be a dynamical feature associated with vertical atmospheric motions. Consistent with prior analyses based on lower stratospheric temperature fields, large-scale standing patterns are also found in the extratropics of both hemispheres, correlated with ozone fluctuations over the equatorial west Pacific. In the Northern Hemisphere, a standing pattern is observed extending from the tropical Indian Ocean to the north Pacific, across North America, and down to the equatorial Atlantic Ocean region. This feature is most pronounced in the NH summer.

Gao, X. H.↗

Navier-Stokes simulation of a close-coupled canard-wing-body configuration

The thin-layer Navier-Stokes equations are solved for the flow about a coplanar close-coupled canard-wing-body configuration at a transonic Mach number of 0.90 and at angles of attack ranging from 0 to 12 degrees. The influence of the canard on the wing flowfield, including canard-wing vortex interaction and wing vortex breakdown, is investigated. A study of canard downwash and canard leading-edge vortex effects, which are the primary mechanisms of the canard-wing interaction, is emphasized. Comparisons between the computations and experimental measurements of surface pressure coefficients, lift, drag and pitching moment data are favorable. A grid refinement study for configurations with and without canard shows that accurate results are obtained using a refined grid for angles of attack where vortex burst is present. At an angle of attack of approximately 12 deg, favorable canard-wing interaction which delays wing vortex breakdown is indicated by the computations and is in good agreement with experimental findings.

Tu, Eugene L.↗

Algorithms for contours depicting static electric fields during adverse weather conditions

A flexible and functional analytical tool is developed for the study of electric fields during adverse weather conditions. This tool is designed for use by members of the Atmospheric Science Group as part of their overall effort to appraise environmental conditions during these situations. It is also used to illustrate approaches open to those interested in the study of the physics of ambient electric field phenomena. Computer resources of KSC are coordinated with original software to produce contour interpretations of electric field data available from a grid of field mills spanning the region. Three model algorithms are presented and examples are given illustrating the system design, flexibility, and utility.

Rompala, John T.↗

Navier-Stokes prediction of large-amplitude delta-wing roll oscillations characterizing wing rock

High-incidence vortical flow about a 65-deg sweep delta wing undergoing static roll and large-amplitude, high-rate-of-roll oscillations is simulated numerically using the time-dependent, three-dimensional, Reynolds-averaged, Navier-Stokes equations. Turbulent computations are presented for static roll angles up through 42 degrees. The effects of roll angle on the vortex aerodynamics are discussed, and the solution accuracy is evaluated by comparison with experimental data. The effects of grid refinement and zonal boundary condition treatment on solution accuracy are assessed at zero roll angle. Numerical simulation of a forced periodic roll motion is also presented.

Chaderjian, Neal M.↗

Evaluation of a multigrid-based Navier-Stokes solver for aerothermodynamic computations

A multigrid acceleration technique developed for solving the three-dimensional Navier-Stokes equations is used for computing high Mach number flows over configurations of practical interest. An explicit multistage Runge-Kutta type of time-stepping scheme is used as the basic algorithm. Solutions are presented for a spherically blunted cone at Mach 10 and a modified shuttle orbiter at Mach 6. The computed surface heat-transfer distributions are shown to compare favorably with the experimental data. Effect of grid-refinement on computed heat-transfer distributions is also examined to assess the numerical accuracy of the computed solutions. The rapid convergence rate associated with multigrid schemes in previous applications at transonic speeds is observed at the higher Mach number flows investigated here.

Vatsa, Veer N.↗

Lower stratosphere waves during 1986-1989 southern springs

Global NMC height analyses and Nimbus 7 TOMS Grid-T Version 6 data are used to investigate the horizontal structure of ozone and 50-mb height waves for the 1986-1989 springs of the Southern Hemisphere. Empirical orthogonal functions (EOFs) are obtained for the ozone and height data independently to identify the meridional structure of observed long waves as well as their propagation characteristics. The leading EOFs for both the ozone and 50-mb height waves exhibit wavenumber 1 and 2 patterns. The correlation coefficients between ozone and 50-mb height principal components are high. They indicate high/low ozone values associated with high/low values of 50-mb heights. This distribution is consistent with the advection of ozone by planetary waves south of the ozone maximum, as shown by a simple linear wave model. The observed total ozone amplitudes for waves 1 and 2 for most years are consistent with the total ozone predictions by the wave model.

Nogues-Paegle, Julia↗

Implementation of the Baldwin-Barth turbulence model into the ZETA code and its diagnosis

The Baldwin-Barth turbulence model was implemented into Zeta, a time-accurate, zonal, integro-differential code for incompressible laminar and turbulent flows. The implementation procedure is patterned after the model subroutine in ARC2D. The results of ZETA with the Baldwin-Barth turbulence model were compared with experimental data, with ZETA using Baldwin-Lomax model, and with ARC2D using the Baldwin-Barth model. The Baldwin-Barth model subroutine was tested by inputting an ARC2D velocity solution of an NACA-0012 airfoil at R(sub e) = 3.9 x 10(exp 6) and alpha = 5 deg. The resultant turbulent viscosity and Reynolds stresses compared favorably with the original data. For the same grid having grid points inside the laminar sublayer, which is necessary due to the one-equation nature of the model, ZETA however predicts early separation. It was found that the current ZETA has problem with such a fine grid. Further work is in progress to solve this problem.

Low, Scott L.↗

GridTool: A surface modeling and grid generation tool

GridTool is designed around the concept that the surface grids are generated on a set of bi-linear patches. This type of grid generation is quite easy to implement, and it avoids the problems associated with complex CAD surface representations and associated surface parameterizations. However, the resulting surface grids are close to but not on the original CAD surfaces. This problem can be alleviated by projecting the resulting surface grids onto the original CAD surfaces. GridTool is designed primary for unstructured grid generation systems. Currently, GridTool supports VGRID and FELISA systems, and it can be easily extended to support other unstructured grid generation systems. The data in GridTool is stored parametrically so that once the problem is set up, one can modify the surfaces and the entire set of points, curves and patches will be updated automatically. This is very useful in a multidisciplinary design and optimization process. GridTool is written entirely in ANSI 'C', the interface is based on the FORMS library, and the graphics is based on the GL library. The code has been tested successfully on IRIS workstations running IRIX4.0 and above. The memory is allocated dynamically, therefore, memory size will depend on the complexity of geometry/grid. GridTool data structure is based on a link-list structure which allows the required memory to expand and contract dynamically according to the user's data size and action. Data structure contains several types of objects such as points, curves, patches, sources and surfaces. At any given time, there is always an active object which is drawn in magenta, or in their highlighted colors as defined by the resource file which will be discussed later.

Samareh-Abolhassani, Jamshid↗

Dynamics of the transition zone in coastal zone color scanner-sensed ocean color in the North Pacific during oceanographic spring

A transition zone in phytoplankton concentration running across the North Pacific basin at 30 deg to 40 deg north latitude corresponds to a basin-wide front in surface chlorophyll observed in a composite of coastal zone color scanner (CZCS) images for May, June, and July 1979-1986. This transition zone with low chlorophyll to the south and higher chlorophyll to the north can be simulated by a simple model of the concentration of phytoplankton, zooplankton, and dissolved nutrient (nitrate) in the surface mixed layer of the ocean applied to the North Pacific basin for the climatological conditions during oceanographic springtime (May, June, and July). The model is initialized with a 1 deg x 1 deg gridded estimate of wintertime (February, March, and April) mixed layer nitrate concentrations calculated from an extensive nutrient database and a similarly gridded mixed layer depth data set. Comparison of model predictions with CZCS data provides a means of evaluating the dynamics of the transition zone. We conclude that in the North Pacific, away from major boundary currents and coastal upwelling zones, wintertime vertical mixing determines the total nutrient available to the plankton ecosystem in the spring. The transition zone seen in basin-scale CZCS images is a reflection of the geographic variation in the wintertime mixed layer depth and the nitracline, leading to a latitudinal gradient in phytoplankton chlorophyll.

Glover, David M.↗

Effects of spatial order of accuracy on the computation of vortical flowfields

The objective of the present investigation is to assess the effect of the spatial order of accuracy used for the evaluation of the inviscid fluxes on the resolution of higher order quantitites, such as velocity gradients. The viscous terms are computed as second-order accurate with central difference formulas, even though for the explicit part of the algorithm higher order approximations may be used. A viscous/inviscid method is used, and the outer part of the flowfield is computed with the inviscid flow equations. The viscous boundary-layer type flow region close to the body surface is computed with an algebraic eddy viscosity model. Results obtained with the conservative and nonconservative formulations and the viscous/inviscid approach are compared with available experimental data. The effect of grid refinement on the accuracy of the solution is also presented.

Ekaterinaris, J. A.↗