Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid size”

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 37 records · Page 2

Net diffusivity in ocean general circulation models with nonuniform grids

The numerical vertical diffusivity K(num), embedded in a numerical ocean general circulation model with nonuniform vertical grid, is estimated. It is shown that in a downwelling region, K(num) is negative for a grid with grid size increasing with depth. When the grid size increment, or the downward vertical velocity, is large, K(num) may exceed the vertical diffusivity specified and may result in a negative effective vertical diffusivity. Therefore care needs to be taken to specify the vertical diffusivity in a numerical model with nonuniform grid, and a lower bound is generally imposed in order to avoid an unphysical negative value. Some possible effects of the negative effective diffusivity are discussed.

Yin, F. L.↗

Accuracy and speed in computing the Chebyshev collocation derivative

We studied several algorithms for computing the Chebyshev spectral derivative and compare their roundoff error. For a large number of collocation points, the elements of the Chebyshev differentiation matrix, if constructed in the usual way, are not computed accurately. A subtle cause is is found to account for the poor accuracy when computing the derivative by the matrix-vector multiplication method. Methods for accurately computing the elements of the matrix are presented, and we find that if the entities of the matrix are computed accurately, the roundoff error of the matrix-vector multiplication is as small as that of the transform-recursion algorithm. Results of CPU time usage are shown for several different algorithms for computing the derivative by the Chebyshev collocation method for a wide variety of two-dimensional grid sizes on both an IBM and a Cray 2 computer. We found that which algorithm is fastest on a particular machine depends not only on the grid size, but also on small details of the computer hardware as well. For most practical grid sizes used in computation, the even-odd decomposition algorithm is found to be faster than the transform-recursion method.

Don, Wai-Sun↗

Characterization of silicon-gate CMOS/SOS integrated circuits processed with ion implantation

The procedure used to generate MEBES masks and produce test wafers from the 10X Mann 1600 Pattern Generator Tape using existing CAD utility programs and the MEBES machine in the RCA Solid State Technology Center are described. The test vehicle used is the MSFC-designed SC102 Solar House Timing Circuit. When transforming the Mann 1600 tapes into MEBES tapes, extreme care is required in order to obtain accurate minimum linewidths when working with two different coding systems because the minimum grid sizes may be different for the two systems. The minimum grid sizes are 0.025 mil for MSFC Mann 1600 and 0.02 mil for MEBES. Some snapping to the next grid is therefore inevitable, and the results of this snapping effect are significant when submicron lines are present. However, no problem was noticed in the SC102 circuit because its minimum linewidth is 0.3 mil (7.6 microns). MEBES masks were fabricated and wafers were processed using the silicon-gate CMOS/SOS and aluminum-gate COS/MOS processing.

Woo, D. S.↗

Problem size, parallel architecture and optimal speedup

The communication and synchronization overhead inherent in parallel processing can lead to situations where adding processors to the solution method actually increases execution time. Problem type, problem size, and architecture type all affect the optimal number of processors to employ. The numerical solution of an elliptic partial differential equation is examined in order to study the relationship between problem size and architecture. The equation's domain is discretized into n sup 2 grid points which are divided into partitions and mapped onto the individual processor memories. The relationships between grid size, stencil type, partitioning strategy, processor execution time, and communication network type are analytically quantified. In so doing, the optimal number of processors was determined to assign to the solution, and identified (1) the smallest grid size which fully benefits from using all available processors, (2) the leverage on performance given by increasing processor speed or communication network speed, and (3) the suitability of various architectures for large numerical problems.

Nicol, David M.↗

Problem size, parallel architecture, and optimal speedup

The communication and synchronization overhead inherent in parallel processing can lead to situations where adding processors to the solution method actually increases execution time. Problem type, problem size, and architecture type all affect the optimal number of processors to employ. The numerical solution of an elliptic partial differential equation is examined in order to study the relationship between problem size and architecture. The equation's domain is discretized into n sup 2 grid points which are divided into partitions and mapped onto the individual processor memories. The relationships between grid size, stencil type, partitioning strategy, processor execution time, and communication network type are analytically quantified. In so doing, the optimal number of processors was determined to assign to the solution, and identified (1) the smallest grid size which fully benefits from using all available processors, (2) the leverage on performance given by increasing processor speed or communication network speed, and (3) the suitability of various architectures for large numerical problems.

Nicol, David M.↗

Turbulent compressible convection in a deep atmosphere. III - Tests on the validity and limitation of the numerical approach

The effects of various approximations on the results of three-dimensional computation of turbulent compressible convection are examined and related to the phenomenon of bubble breaking in stellar convection. The problem of viscosity and the numerical issue of efficiency versus accuracy in such computations are discussed, and a model is formulated and tested with regard to turbulent viscosity and accuracy. Results on the convergence to a statistically stationary state, the effects of time accuracy on the statistics, the effects of the adjustable constant parameter c(mu), and the effects of aspect ratio, horizontal grid size, and vertical grid distribution are discussed. A preliminary look at the three-dimensional flow structure is given.

Chan, K. L.↗

Interaction of decaying trailing vortices in ground shear

The drift of trailing vortices in a cross wind near to the ground is modeled by an unsteady, two-dimensional, rotational flowfield with a concentration of large vorticity in spots having finite total strength. The problem is analyzed by a combination of the method of matched asymptotic analyses for the decay of the vortical spots and the Euler solution for the unsteady rotational flow. A primary feature of the numerical scheme is that the grid size and time step depend only the length and velocity scales of the background flow and is independent of the effective core size of a vortical spot which can be much smaller than the grid size. Numerical results are presented to demonstrated the strong interaction between the trajectjories of the vortical spots and the redistribution of vorticity in the background flowfield.

Liu, C. H.↗

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

Direct numerical simulations (DNS) and large-eddy simulations (LES) simulations of flow through a periodic channel with a constriction are performed using the dynamic Smagorinsky model at two Reynolds numbers of 2800 and 10595. The LES equations are solved using higher order compact schemes. DNS are performed for the lower Reynolds number case using a fine grid and the data are used to validate the LES results obtained with a coarse and a medium size grid. LES simulations are also performed for the higher Reynolds number case using a coarse and a medium size grid. The results are compared with an existing reference data set. The DNS and LES results agreed well with the reference data. Reynolds stresses, sub-grid eddy viscosity, and the budgets for the turbulent kinetic energy are also presented. It is found that the turbulent fluctuations in the normal and spanwise directions have the same magnitude. The turbulent kinetic energy budget shows that the production peaks near the separation point region and the production to dissipation ratio is very high on the order of five in this region. It is also observed that the production is balanced by the advection, diffusion, and dissipation in the shear layer region. The dominant term is the turbulent diffusion that is about two times the molecular dissipation.

Balakumar, P.↗

Motion and interaction of decaying trailing vortices in spanwise shear wind

A simulation is presented of the drift of trailing vortices in a cross-wind near the ground by an unsteady, two-dimensional, rotational flow field with a concentration of large vorticity in vortical spots (having a finite but small effective size and finite total strength). The problem is analyzed by a combination of the method of matched asymptotic analyses for the decay of the vortical spots and the Euler solution for the unsteady rotational flow. Using the method of averaging, a special numerical method is developed in which the grid size and time step depend only on the length and velocity scales of the background flow and are independent of the effective core size of a vortical spot. The core size can be much smaller than the grid size, whereas the peak velocity in the core is inversely propertional to the spot size. Numerical results are presented to demonstrate the strong interaction between the trajectories of the vortical spots and the change of the vorticity distribution in the background flow field.

Liu, C. H.↗

CUDA Simulation of Homogeneous, Incompressible Turbulence

We discuss very fast Compute Unified Device Architecture (CUDA) simulations of ideal homogeneous incompressible turbulence based on Fourier models. These models have associated statistical theories that predict that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. Prior numerical simulations have shown that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We review the theoretical basis of this "broken ergodicity" as applied to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence. Our new simulations examine the phenomenon of broken ergodicity through very long time and large grid size runs performed on a state-of-the-art CUDA platform. Results comparing various CUDA hardware configurations and grid sizes are discussed. NS and MHD results are compared.

Morin, Lee↗

Polar microwave brightness temperatures from Nimbus-7 SMMR: Time series of daily and monthly maps from 1978 to 1987

A time series of daily brightness temperature gridded maps (October 25, 1978 through August 15, 1987) were generated from all ten channels of the Nimbus-7 Scanning Multichannel Microwave Radiometer orbital data. This unique data set can be utilized in a wide range of applications including heat flux, ocean circulation, ice edge productivity, and climate studies. Two sets of data in polar stereographic format are created for the Arctic region: one with a grid size of about 30 km on a 293 by 293 array similar to that previously utilized for the Nimbus-5 Electrically Scanning Microwave Radiometer, while the other has a grid size of about 25 km on a 448 by 304 array identical to what is now being used for the DMSP Scanning Multichannel Microwave Imager. Data generated for the Antaractic region are mapped using the 293 by 293 grid only. The general technique for mapping, and a quality assessment of the data set are presented. Monthly and yearly averages are also generated from the daily data and sample geophysical ice images and products derived from the data are given. Contour plots of monthly ice concentrations derived from the data for October 1978 through August 1987 are presented to demonstrate spatial and temporal detail which this data set can offer, and to show potential research applications.

Comiso, Josefino C.↗

Nonlinear vibrations of rectangular plates.

A finite-difference method is developed to determine the large amplitude dynamic responses of thin elastic plates subjected to uniform pressure pulse-type loads. Four different sets of boundary conditions are considered. Some specific problems are solved. The results are compared with approximate solutions obtained by Yamaki (1961). The numerical method presented provides an accurate and efficient approximate solution to the problem, and should be useful as a check on other approximate methods. The grid-size and the time-step necessary for obtaining numerical stability depend on the particular problem. For many cases the method converges rapidly and a rather large grid-size and time-step is adequate.

Bayles, D. J.↗

Influence of numerical dissipation in computing supersonic vortex-dominated flows

Steady supersonic vortex-dominated flows are solved using the unsteady Euler equations for conical and three-dimensional flows around sharp- and round-edged delta wings. The computational method is a finite-volume scheme which uses a four-stage Runge-Kutta time stepping with explicit second- and fourth-order dissipation terms. The grid is generated by a modified Joukowski transformation. The steady flow solution is obtained through time-stepping with initial conditions corresponding to the freestream conditions, and the bow shock is captured as a part of the solution. The scheme is applied to flat-plate and elliptic-section wings with a leading edge sweep of 70 deg at an angle of attack of 10 deg and a freestream Mach number of 2.0. Three grid sizes of 29 x 39, 65 x 65 and 100 x 100 have been used. The results for sharp-edged wings show that they are consistent with all grid sizes and variation of the artificial viscosity coefficients. The results for round-edged wings show that separated and attached flow solutions can be obtained by varying the artificial viscosity coefficients. They also show that the solutions are independent of the way time stepping is done. Local time-stepping and global minimum time-steeping produce same solutions.

Kandil, O. A.↗

Numerical simulation of precipitation induced downbursts

Using the Terminal Area Simulation System (TASS), numerical simulations of downburst structure and sensitivity, based on vertical profiles of environmental temperature, humidity and wind velocity observed during the June and August 1982 JAWS project, are presented. Two-dimensional axisymmetric simulations examining downburst evolution, structure and sensitivity, assume a 40-m constant grid size on a 10-km diameter 5-km deep cylindrical domain. The three-dimensional experiment, examining the effects of vertical wind shear and other asymmetrical aspects of the downpour, assumed a 500-m constant horizontal grid size and a 35 x 35 x 18.5-km area. In the downburst primary structure depicted, outflow speeds were found to be sensitive to environmental temperature and humidity, as well as to precipitation radius and intensity. A vortex ring was found to propagate downwards, and maximum outflow winds occur when the vortex ring first reaches the surface.

Proctor, F. H.↗

Convective Systems Over the South China Sea: Cloud-Resolving Model Simulations

The South China Sea Monsoon Experiment (SCSMEX) was conducted in May-June 1998. One of its major objectives is to better understand the key physical processes for the onset and evolution of the summer monsoon over Southeast Asia and southern China. Multiple observation platforms (e.g., upper-air soundings, Doppler radar, ships, wind profilers, radiometers, etc.) during SCSMEX provided a first attempt at investigating the detailed characteristics of convective storms and air pattern changes associated with monsoons over the South China Sea region. SCSMEX also provided rainfall estimates which allows for comparisons with those obtained from the Tropical Rainfall Measuring Mission (TRMM), a low earth orbit satellite designed to measure rainfall from space. The Goddard Cumulus Ensemble (GCE) model (with 1-km grid size) is used to understand and quantify the precipitation processes associated with the summer monsoon over the South China Sea. This is the first (loud-resolving model used to simulate precipitation processes in this particular region. The GCE-model results captured many of the observed precipitation characteristics because it used a fine grid size. For example, the temporal variation of the simulated rainfall compares quite well to the sounding-estimated rainfall variation. The time and domain-averaged temperature (heating/cooling) and water vapor (drying/ moistening) budgets are in good agreement with observations. The GCE-model-simulated rainfall amount also agrees well with TRMM rainfall data. The results show there is more evaporation from the ocean surface prior to the onset of the monsoon than after the on-~et of monsoon when rainfall increases. Forcing due to net radiation (solar heating minus longwave cooling) is responsible for about 25% of the precipitation in SCSMEX The transfer of heat from the ocean into the atmosphere does not contribute significantly to the rainfall in SCSMEX. Model sensitivity tests indicated that total rain production is reduced 17-18% in runs neglecting the ice phase. The SCSMEX results are compared to other GCE-model-simulated weather systems that developed during other field campaigns (i.e., west Pacific warm pool region, eastern Atlantic region and central USA). Large-scale forcing vie temperature and water vapor tendency, is the major energy source for net condensation in the tropical cases. The effects of large-scale cooling exceed that of large-scale moistening in the west pacific warm pool region and eastern Atlantic region. For SCSMEX, however, the effects of large-scale moistening predominate. Net radiation and sensible and latent hc,it fluxes play a much more important role in the central USA.

Tao, Wei-Kuo↗

Session on validation of coupled models

The session on validation of coupled models is reviewed. The current use of a mesoscale model with a grid size of 20-km during STORM-FEST in 1992 has proven to be extremely valuable. The availability of forecast products at a much higher temporal and spatial resolution was very helpful for mesoscale forecasting, mission planning, and the guidance of research aircraft. Recent numerical simulation of ocean cyclones and mesoscsle convective systems using nonhydrostatic cloud/mesoscale models with a grid size as small as 2-km have demonstrated the potential of these models for predicting mesoscale convective systems, squall lines, hurricane rainbands, mesoscale gravity waves, and mesoscale frontal structures embedded within an extratropical cyclone. Although mesoscale/cloud scale models have demonstrated strong potential for use in operational forecasting, very limited quantitative evaluation (and verification) of these models were performed. As a result, the accuracy, the systematic biases, and the useful forecasts limits were not properly defined for these models. Also, no serious attempts were made to use these models for operational prediction of mesoscale convective systems.

Kuo, Bill↗

Computation of saddle point of attachment

Low-speed flows over a cylinder mounted on a flat plate are studied numerically in order to confirm the existence of a saddle point of attachment in the flow before an obstacle, to analyze the flow characteristics near the saddle point theoretically, and to address the significance of the saddle point of attachment to the construction of external flow structures, the interpretation of experimental surface oil-flow patterns, and the theoretical definition of three-dimensional flow separation. Two numerical codes, one for an incompressible flow and another for a compressible flow, are used for various Mach numbers, Reynolds numbers, grid sizes, and numbers of grid points. It is pointed out that the potential presence of a saddle point of attachment means that a line of 'oil accumulation' from both sides of a skin-friction line emanating outward from a saddle point can be either a line of separation or a line of attachment.

Hung, Ching-Mao↗