Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large-scale 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 235 records · Page 13

Turbulent magnetic reconnection

The results of two-dimensional spectral-method MHD simulations of the effect of turbulence on magnetic-reconnection processes are presented in extensive graphs and tables and characterized in detail. In a periodic sheet-pinch configuration, early-onset, persistent turbulence is found to occur when a low level of broadband fluctuations is introduced into the initial data. Features observed include small-scale unsteady coherent electric-current and vorticity structures, enhanced viscous and resistive dissipation, unsteady and spatially asymmetric fluid flow, large-scale and small-scale magnetic islands, multiple X points, large-amplitude electric-field fluctuations, and spectral characteristics indicating that the turbulence is broadband in both spatial directions and similar to homogeneous MHD turbulence. The importance of reconnection for studies of solar, magnetospheric, and laboratory plasmas in inhomogeneous magnetic fields is indicated.

Matthaeus, W. H.↗

A comparative study of rapidly and slowly rotating dynamical regimes in a terrestrial general circulation model

As a preliminary step in the development of a general circulation model for general planetary use, a simplified vesion of thef GISS Model I GCM has been run at various rotation periods to investigate differences between the dynamical regimes of rapidly and slowly rotating planets. To isolate the dynamical processes, the hydrologic cycle is suppressed and the atmosphere is forced with perpetual annual mean solar heating. All other parameters except the rotation period remain fixed at their terrestrial values. Experiments were conducted for rotation periods of 2/3, 1, 2, 4, 8, 16, 64 and 256 days. The results are in qualitative agreement with similar experiments carried out previously with other GCMs and with certain aspects of one Venus GCM simulation. As rotation rate decreases, the energetics shifts from baroclinc to quasi-barotropic when the Rossby radius of deformation reaches planetary scale. The Hadley cell expands poleward and replaces eddies as the primary mode of large-scale heat transport. Associated with this is a poleward shift of the baroclinic zone and jet stream and a reduction of the equator-pole temperature contrast. Midlatitude jet strength peaks at 8 days period, as does the weak positive equatorial zonal wind which occurs at upper levels at all rotation periods. Eddy momentum transport switches from poleward to equatorward at the same period. Tropospheric mean static stability generally increases in the tropics and decreases in midlatitudes as rotation rate decreases, but the global mean static stability is independent of rotation rate. The peak in the eddy kinetic energy spectrum shifts toward lower wavenumbers, reaching wavenumber 1 at a period of 8 days. Implications of these results for the dynamics of Venus and Titan are discussed. Specifically, it is suggested that the extent of low-level convection determines whether the Gierasch mechanism contributes significantly to equatorial superrotation on these planets.

Del Genio, Anthony D.↗

A quantitative approach to the topology of large-scale structure

A quantitative measure of the topology of large-scale structure: the genus of density contours in a smoothed density distribution, is described and applied. For random phase (Gaussian) density fields, the mean genus per unit volume exhibits a universal dependence on threshold density, with a normalizing factor that can be calculated from the power spectrum. If large-scale structure formed from the gravitational instability of small-amplitude density fluctuations, the topology observed today on suitable scales should follow the topology in the initial conditions. The technique is illustrated by applying it to simulations of galaxy clustering in a flat universe dominated by cold dark matter. The technique is also applied to a volume-limited sample of the CfA redshift survey and to a model in which galaxies reside on the surfaces of polyhedral 'bubbles'. The topology of the evolved mass distribution and 'biased' galaxy distribution in the cold dark matter models closely matches the topology of the density fluctuations in the initial conditions. The topology of the observational sample is consistent with the random phase, cold dark matter model.

Gott, J. Richard, III↗

The topology of large-scale structure. I - Topology and the random phase hypothesis

Many models for the formation of galaxies and large-scale structure assume a spectrum of random phase (Gaussian), small-amplitude density fluctuations as initial conditions. In such scenarios, the topology of the galaxy distribution on large scales relates directly to the topology of the initial density fluctuations. Here a quantitative measure of topology - the genus of contours in a smoothed density distribution - is described and applied to numerical simulations of galaxy clustering, to a variety of three-dimensional toy models, and to a volume-limited sample of the CfA redshift survey. For random phase distributions the genus of density contours exhibits a universal dependence on threshold density. The clustering simulations show that a smoothing length of 2-3 times the mass correlation length is sufficient to recover the topology of the initial fluctuations from the evolved galaxy distribution. Cold dark matter and white noise models retain a random phase topology at shorter smoothing lengths, but massive neutrino models develop a cellular topology.

Weinberg, David H.↗

Double inflation - A possible resolution of the large-scale structure problem

A model is presented for the large-scale structure of the universe in which two successive inflationary phases resulted in large small-scale and small large-scale density fluctuations. This bimodal density fluctuation spectrum in an Omega = 1 universe dominated by hot dark matter leads to large-scale structure of the galaxy distribution that is consistent with recent observational results. In particular, large, nearly empty voids and significant large-scale peculiar velocity fields are produced over scales of about 100 Mpc, while the small-scale structure over less than about 10 Mpc resembles that in a low-density universe, as observed. Detailed analytical calculations and numerical simulations are given of the spatial and velocity correlations.

Turner, Michael S.↗

Simultaneous density and electric field fluctuation spectra associated with velocity shears in the auroral oval

A detailed study is presented of simultaneous density and electric field fluctuation spectra over a large-scale length range seen in association with large structured convective plasma flows, field-aligned currents, and particle precipitation at high latitudes. The data were obtained for two Dynamics Explorer 2 orbits at two different altitudes within the F region and the topside ionosphere. The observations are compared with results of nonlinear simulations of shear flow-driven instabilities and predictions based on two-dimensional turbulence arguments, with particular reference to the Kelvin-Helmholtz process.

Basu, Sunanda↗

Simulation of ocean SAR images via phase history generation

A method for simulating an ocean synthetic-aperture radar (SAR) image is illustrated for a simple internal wave current pattern. The method calculates both the amplitude image and the radar signal history. The simulation model used consists of six stages. In the first stage, the full wave spectrum in two spatial coordinates is calculated from the wind speed and direction using the action spectral density equation. In the second stage, the pixel size is selected and used as the basis for dividing the spectrum into large and small scale motions. Realizations for the large-scale ocean surface height and velocity are then calculated. In the third stage, the sensor wavelength and geometry are used to calculate the small scale statistics (radar cross section, coherence time, root-mean square (RMS), radial velocity and RMS slope). In the fourth stage, an autoregressive method is used to generate a realization of the surface reflectivity history that is consistent with the radial velocity and the radial velocity variance. In the fifth stage, the signal history is generated by summing the reflectivities at the proper times with the appropriate antenna weighting function for the SAR. As a consequence of this process, speckle is automatically included in the signal. In the final stage, the SAR image can be created by using any of the traditional ways to process the signal history, including variable focusing and multilook processing.

Bennett, John R.↗

Large-scale experimental technology with remote sensing in land surface hydrology and meteorology

Two field experiments to study atmospheric and land surface processes and their interactions are summarized. The Hydrologic-Atmospheric Pilot Experiment, which tested techniques for measuring evaporation, soil moisture storage, and runoff at scales of about 100 km, was conducted over a 100 X 100 km area in France from mid-1985 to early 1987. The first International Satellite Land Surface Climatology Program field experiment was conducted in 1987 to develop and use relationships between current satellite measurements and hydrologic, climatic, and biophysical variables at the earth's surface and to validate these relationships with ground truth. This experiment also validated surface parameterization methods for simulation models that describe surface processes from the scale of vegetation leaves up to scales appropriate to satellite remote sensing.

Brutsaert, Wilfried↗

The formation of topological defects in phase transitions

It was argued, and fought through numerical work that the results of non-dynamical Monte Carlo computer simulations cannot be applied to describe the formation of topological defects when the correlation length at the Ginzburg temperature is significantly smaller than the horizon size. To test the current hypothesis that infinite strings at formation are essentially described by Brownian walks of size the correlation length at the Ginzburg temperature, fields at the Ginzburg temperature were equilibrated. Infinite structure do not exist in equilibrium for reasonable definitions of the Ginzburg temperature, and horizons must be included in a proper treatment. A phase transition, from small-scale to large-scale string or domain wall structure, is found to occur very close to the Ginzburg temperature, in agreement with recent work. The formation process of domain walls and global strings were investigated through the breaking of initially ordered states. To mimic conditions in the early Universe, cooling times are chosen so that horizons exist in the sample volume when topological structure formation occurs. The classical fields are evolved in real-time by the numerical solution of Langevin equations of motion on a three dimensional spatial lattice. The results indicate that it is possible for most of the string energy to be in small loops, rather than in long strings, at formation.

Hodges, Hardy M.↗

A numerical study of a class of TVD schemes for compressible mixing layers

At high Mach numbers the two-dimensional time-developing mixing layer develops shock waves, positioned around large-scale vortical structures. A suitable numerical method has to be able to capture the inherent instability of the flow, leading to the roll-up of vortices, and also must be able to capture shock waves when they develop. Standard schemes for low speed turbulent flows, for example spectral methods, rely on resolution of all flow-features and cannot handle shock waves, which become too thin at any realistic Reynolds number. The performance of a class of second-order explicit total variation diminishing (TVD) schemes on a compressible mixing layer problem was studied. The basic idea is to capture the physics of the flow correctly, by resolving down to the smallest turbulent length scales, without resorting to turbulence or sub-grid scale modeling, and at the same time capture shock waves without spurious oscillations. The present study indicates that TVD schemes can capture the shocks accurately when they form, but (without resorting to a finer grid) have poor accuracy in computing the vortex growth. The solution accuracy depends on the choice of limiter. However a larger number of grid points are in general required to resolve the correct vortex growth. The low accuracy in computing time-dependent problems containing shock waves as well as vortical structures is partly due to the inherent shock-capturing property of all TVD schemes. In order to capture shock waves without spurious oscillations these schemes reduce to first-order near extrema and indirectly produce clipping phenomena, leading to inaccuracy in the computation of vortex growth. Accurate simulation of unsteady turbulent fluid flows with shock waves will require further development of efficient, uniformly higher than second-order accurate, shock-capturing methods.

Sandham, N. D.↗

Breakdown of vertically propagating two-dimensional gravity waves forced by orography

The propagation of orographic gravity waves into an atmosphere with exponentially decreasing density is simulated with a two-dimensional, nonlinear, time-dependent numerical model. After the stationary wave is established over the mountain, the model predicts that wave breaking causes a large reduction of the vertical momentum flux in the flow, not only at levels where wave breaking is present, but also far below the lowest occurrence of overturning. More than half of the decrease in momentum flux is explained by the presence of large amplitude, downward propagating waves, which are generated in regions of wave breaking. The downward propagating waves appear almost simultaneously with overturning, and have nonzero phase speeds, suggesting a strongly nonlinear generation mechanism that depends on local wave properties. The generation of these downward propagating waves is a robust process, insensitive to mountain height, mountain width, or density scale height. These results have important implications for observational studies of orographically generated waves as well as for schemes that seek to parameterize the effects of orogoraphy in large-scale models.

Bacmeister, Julio T.↗

The effect of latent heat release on synoptic-to-planetary wave interactions and its implication for satellite observations: Theoretical modeling

Simple models are being developed to simulate interaction of planetary and synoptic-scale waves incorporating the effects of large-scale topography; eddy heat and momentum fluxes (or nonlinear dynamics); radiative heating/cooling; and latent heat release (precipitation) in synoptic-scale waves. The importance of latent heat release is determined in oceanic storm tracks for temporal variability and time-mean behavior of planetary waves. The model results were compared with available observations of planetary and synoptic-scale wave variability and time-mean circulation. The usefulness of monitoring precipitation in oceanic storm tracks by satellite observing systems was ascertained. The modeling effort includes two different low-order quasi-geostrophic models-time-dependent version and climatological mean version. The modeling also includes a low-order primitive equation model. A time-dependent, multi-level version will be used to validate the two-level Q-G models and examine effects of spherical geometry.

Branscome, Lee E.↗

Large-scale advanced propeller blade pressure distributions: Prediction and data

Two Euler analyses techniques, finite difference and finite volume, were employed to predict the blade surface pressure distributions of a large scale advanced propeller. The predicted pressure distributions were compared with wind tunnel data. Both techniques produced blade pressure distributions that are in fairly good agreement with the data over the range of test Mach numbers of 0.2 to 0.78. However, the numerical simulations fail to predict correctly the measured pressure distributions for the low Mach number, high power case. The data indicate the presence of a leading edge vortex for this case. A discussion of the compressibility effects is also presented.

Nallasamy, M.↗

The stability of a collisionless cosmological shell

The P3 M technique is used here to simulate the evolution of collisionless shells in an Omega = 1 universe. Starting from the spherical similarity solution, a bootstrap technique is used to follow the evolution over very large expansion factors. It is found that the overall structure follows the similarity solution for a long period during which bound clumps grow within the shell. At late times the growth of structure depends on induced velocity perturbations in material outside the shell. If such perturbations are suppressed, structure on the shell becomes self-similar. When induced motions in the background medium are included, the evolution at late times is dominated by large-scale modes as predicted by linear stability analysis. The stable final state appears to consist of one or two massive clumps on the edge of a spherical void. The possible application of these results to the origin of galaxies and large-scale structure is discussed.

White, Simon D. M.↗

Multilayer sea ice in the Arctic - Model- and satellite-derived

The behavior of sea ice is a strong function of ice strength which, in the central Arctic, depends on the relative fractions of multiyear and first-year ice coverage. A dynamic-thermodynamic sea ice model has been configured to include both multiyear and first-year ice. The model has been run for a 7-year period, 1979-1985, for which a satellite microwave record of multiyear sea ice coverage is available. The simulated concentrations of multiyear ice are slightly higher than the satellite-derived concentrations, although the total coverage of simulated multiyear ice agrees well with the satellite-derived coverage. Both the simulated and the satellite-derived coverages in the early 1980s show large interannual variations that represent potentially large perturbations in the mass balance of Arctic pack ice. These fluctuations are generally consistent in the model and satellite data, and they are generally consistent with year-to-year differences in the large-scale wind forcing. In at least some cases the wintertime distribution of multiyear ice is influenced by the pattern of ice drift as far back as the previous summer. The late December coverage of multiyear ice in the Alaskan region appears to foreshadow summer ice conditions in the same region. The short-term response to an intense synoptic system is considerably larger in the satellite data than in the model results.

Walsh, John E.↗

Compiler-directed cache management in multiprocessors

The necessity of finding alternatives to hardware-based cache coherence strategies for large-scale multiprocessor systems is discussed. Three different software-based strategies sharing the same goals and general approach are presented. They consist of a simple invalidation approach, a fast selective invalidation scheme, and a version control scheme. The strategies are suitable for shared-memory multiprocessor systems with interconnection networks and a large number of processors. Results of trace-driven simulations conducted on numerical benchmark routines to compare the performance of the three schemes are presented.

Cheong, Hoichi↗

On the large-scale structures formed by wakes of open cosmic strings

Large-scale structures of the universe have been variously described as sheetlike, filamentary, cellular, bubbles or spongelike. Recently cosmic strings became one of viable candidates for a galaxy formation scenario, and some of the large-scale structures seem to be simply explained by the open cosmic strings. According to this scenario, sheets are wakes which are traces of moving open cosmic strings where dark matter and baryonic matter have accumulated. Filaments are intersections of such wakes and high density regions are places where three wakes intersect almost orthogonally. The wakes formed at t sub eq become the largest surface density among all wakes, where t sub eq is the epoch when matter density equals to radiation density. If we assume that there is one open cosmic string per each horizon, then it can be explained that the typical distances among wakes, filaments and clusters are also approx. 10(exp 2) Mpc. This model does not exclude a much more large scale structure. Open cosmic string may move even now and accumulate cold dark matter after its traces. However, the surface density is much smaller than the ones formed at t sub eq. From this model, it is expected that the typical high density region will have extended features such as six filaments and three sheets and be surrounded by eight empty regions (voids). Here, the authors are mainly concerned with such structures and have made numerical simulations for the formation of such large scale structures.

Hara, Tetsuya↗

Transition to hard turbulence in thermal convection at infinite Prandtl number

Direct numerical simulations of two-dimensional high Rayleigh (Ra) number, base-heated thermal convection in large aspect-ratio boxes are presented for infinite Prandtl number fluids, as applied to the earth's mantle. A transition is characaterized in the flow structures in the neighborhood of Ra between 10 to the 7th and 10 to the 8th. These high Ra flows consist of large-scale cells with strong intermittent, boundary-layer instabilities. For Ra exceeding 10 to the 7th it is found that the heat-transfer mechanism changes from one characterized by mushroom-like plumes to one consisting of disconnected ascending instabilities, which do not carry with them all the thermal anomaly from the bottom boundary layer. Plume-plume collisions become much more prominent in high Ra situations and have a tendency of generating a pulse-like behavior in the fixed plume. This type of instability represents a distinct mode of heat transfer in the hard turbulent regime. Predictions of this model can be used to address certain issues concerning the mode of time-dependent convection in the earth's mantle.

Hansen, Ulrich↗