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 775 records · Page 43

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↗

Large-scale aeroacoustic research feasibility and conceptual design of test-section inserts for the Ames 80- by 120-foot wind tunnel

An engineering feasibility study was made of aeroacoustic inserts designed for large-scale acoustic research on aircraft models in the 80 by 120 foot Wind Tunnel at NASA Ames Research Center. The advantages and disadvantages of likely designs were analyzed. Results indicate that the required maximum airspeed leads to the design of a particular insert. Using goals of 200, 150, and 100 knots airspeed, the analysis indicated a 30 x 60 ft open-jet test section, a 40 x 80 ft open jet test section, and a 70 x 100 ft closed test section with enhanced wall lining, respectively. The open-jet inserts would be composed of a nozzle, collector, diffuser, and acoutic wedges incorporated in the existing 80 x 120 test section. The closed test section would be composed of approximately 5 ft acoustic wedges covered by a porous plate attached to the test section walls of the existing 80 x 120. All designs would require a double row of acoustic vanes between the test section and fan drive to attenuate fan noise and, in the case of the open-jet designs, to control flow separation at the diffuser downstream end. The inserts would allow virtually anechoic acoustic studies of large helicopter models, jets, and V/STOL aircraft models in simulated flight. Model scale studies would be necessary to optimize the aerodynamic and acoustic performance of any of the designs. In all designs studied, the existing structure would have to be reinforced. Successful development of acoustically transparent walls, though not strictly necessary to the project, would lead to a porous-wall test section that could be substituted for any of the open-jet designs, and thereby eliminate many aerodynamic and acoustic problems characteristic of open-jet shear layers. The larger size of the facility would make installation and removal of the insert components difficult. Consequently, scheduling of the existing 80 x 120 aerodynamic test section and scheduling of the open-jet test section would likely be made on an annual or longer basis. The enhanced wall-lining insert would likely be permanent. Although the modifications are technically feasible, the economic practicality of the project was not evaluated.

Soderman, Paul T.↗

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↗

A physically based fractional cloudiness parameterization

A bulk PBL model with a simple internal vertical structure and a simple second-order closure has been developed for use as a PBL parameterization in a large-scale model. A method to determine the value of sigma, the fractional area covered by rising motion, is given. Large-eddy simulations of Moeng (1984, 1986) are used to validate various aspects of the parameterization.

Randall, David A.↗

Photophoretic Propulsion Enabling Mesosphere Exploration NIAC Phase I Final Report

This Phase I report presents a comprehensive study on photophoretic flyers—innovative, ultralight, solar-powered vehicles that harness photophoretic forces generated via Knudsen pumping to achieve sustained flight in the mesosphere (50–80 km altitude). By integrating advanced materials such as nanocardboard— characterized by its extremely low areal density (~1 g/m²) and high bending stiffness—with ultrathin light-absorbing coatings, the project converts incident solar radiation directly into a directed thrust. Extensive experimental investigations, coupled with high-fidelity computational fluid dynamics (CFD) simulations using ANSYS Fluent, validate the concept across various three-dimensional geometries, including spherical, conical, and rocket-shaped configurations. These simulations bridge the gap between free-molecular and continuum flow regimes, demonstrating that optimized designs can generate lift forces sufficient to support kilogram-scale payloads even in low-pressure environments. At the heart of this innovation is the use of Knudsen pumping, where temperature gradients across porous surfaces induce directional gas flow, creating a modest overpressure that provides lift. The report introduces an analytical framework that interpolates between the well-known low-Reynolds number drag regime and the high-Reynolds number momentum theory. This model accurately predicts lift based on design parameters such as microchannel dimensions, porous wall geometry, areal density, and nozzle exit area. For instance, simulations indicate that 10-meter-scale structures with carefully engineered porous walls can achieve the necessary pressure differential to support scientifically significant payloads (~1 kg). The study also explores a hybrid propulsion approach that combines solar buoyancy with photophoretic lift. Initially, solar heating creates a buoyant force that elevates the flyer to mesospheric altitudes. Once in the optimal pressure range, the photophoretic mechanism—powered by Knudsen pumping—takes over as the primary source of lift, ensuring stable, long-duration flight. This dual-mode operation not only facilitates the deployment of photophoretic flyers but also broadens the potential applications for mesospheric exploration. In addition to propulsion, the report investigates the integration of photophoretic thrusters for trajectory control of existing research balloons in the upper stratosphere. By dynamically adjusting the nozzle orientation and controlling the flow-through velocity, these thrusters provide precise maneuverability, enabling the flyers to counteract atmospheric disturbances and adjust their flight paths in real time. For example, a photophoretic thruster approximately 7.5 by 7.5 meters in size could be unfolded below a payload gondola of a 60 million-cubic-foot zero-pressure balloon. Such a thruster can provide horizontal speed control of approximately 1 m/s using only sunlight and no moving parts (except those needed to track the Sun and control the jet direction). Importantly, photophoretic thrusters operate more efficiently at higher altitudes, which is complementary to known trajectory control techniques, such as propellers and tethered wings, which are more effective at lower altitudes. Finally, the report identifies three scientific research thrusts where mesospheric aircraft technology can have a profound impact: atmospheric tides, characterization of gravity waves, and investigation of mesospheric instabilities. Overall, the findings of this Phase I project represent a significant advancement in photophoretic propulsion technology. By demonstrating that large-scale, ultralight structures can be powered solely by solar radiation—via carefully engineered Knudsen pumping—this work lays a robust foundation for scalable, near-space flight architectures. Future refinements in material fabrication, structural optimization, and integrated trajectory control are expected to further enhance performance, paving the way for operational demonstrations that could revolutionize atmospheric science, remote sensing, and communication networks.

Knudsen Pump↗

CoLoRe-2LPT: Lyman-$α$ mock catalogues for the validation of DESI cosmological analyses

The Lyman-$α$ (Ly$α$) forest has become a crucial probe for studying the large-scale structure of the universe at high redshift ($z > 2$), providing powerful constraints on Baryon Acoustic Oscillations (BAO) and the full-shape (FS) clustering of matter. As a key ingredient for upcoming BAO and FS analyses, we present a new generation of fast cosmological Ly$α$ mocks based on second-order Lagrangian perturbation theory (2LPT). These new mocks significantly improve upon previous log-normal approaches, both at accurately capturing small scale clustering and at recovering the non-linear broadening of the BAO peak. They are able to reproduce Ly$α$ statistics within $10\%$ of the latest DESI measurement; including the Ly$α$ bias and the redshift-space distortion $β$ parameter, mean transmitted flux, and 1D power spectrum. The corresponding quasar (QSO) clustering is also improved with respect to previous approaches, calibrated against high-resolution Abacus simulations, recovering the observational QSO linear bias to less than $5\%$ and improving redshift-space distortions via 2LPT velocities and the addition of Fingers-of-God effects. Furthermore, these mocks incorporate high column density systems and metal lines, allowing us to explore the effects and systematics induced by these astrophysical contaminants. This new set of mocks has been key for enhancing the modeling and validation of the DESI DR2 Ly$α$ full shape cosmological analysis. This work provides a physically motivated and computationally efficient tool for simulating current and next-generation Ly$α$ surveys and validating FS and BAO analysis.

Bernal, M.F. Ruiz-Herrera [Madrid, CIEMAT] (ORCID:↗