Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time scale separation”

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 469 records · Page 26

Control variates for lattice field theory

In most lattice field theories, correlators are plagued by a signal-to-noise problem of exponential difficulty in the time separation. We propose a method for improving the signal-to-noise ratio, in which control variates are systematically constructed from lattice Schwinger-Dyson relations. The method is demonstrated on various two-dimensional lattices in scalar field theory, and a strategy for scaling to larger systems is explored. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Synchronous and Concurrent Multidomain Computing Method for Cloud Computing Platforms

We present a numerical method for synchronous and concurrent solution of transient elastodynamics problem where the computational domain is divided into subdomains that may reside on separate computational platforms. Here, this work employs the variational multiscale discontinuous Galerkin (VMDG) method to develop interdomain transmission conditions for transient problems. The fine-scale modeling concept leads to variationally consistent coupling terms at the common interfaces. The method admits a large class of time discretization schemes, and decoupling of the solution for each subdomain is achieved by selecting any explicit algorithm. Numerical tests with a manufactured solution problem show optimal convergence rates. The energy history in a free vibration problem is in agreement with that of the solution from a monolithic computational domain.

97 MATHEMATICS AND COMPUTING↗

Global circuit variability from multiple stratospheric electrical measurements

The vertical component of the electric field, conductivity, and the derived current density from eight superpressure balloons are examined. Special emphasis was placed on the fair-weather, simultaneous measurements from widely spaced constant-altitude (26 km) balloons. The conductivity measurements were well organized by a simple ionization rate parameterization depending on the geomagnetic latitude. The variability of all the electrical parameters was found to be independent of the balloon separation. Much of the time the current density measurements were within 20 percent, for simultaneous flights, even though the balloons were up to 6000 km apart. Also, these simultaneous current density data show that the global current source was significantly variable on hourly and daily time scales. Finally, the use of the simultaneous current density data is discussed as a possible 'geoelectric index'. Advantages and limitations of such an index are discussed.

Norville, Kent↗

The effects of forcing on a single stream shear layer and its parent boundary layer

Forcing and its effect on fluid flows has become an accepted tool in the study and control of flow systems. It has been used both as a diagnostic tool, to explore the development and interaction of coherent structures, and as a method of controlling the behavior of the flow. A number of forcing methods have been used in order to provide a perturbation to the flow; among these are the use of an oscillating trailing edge, acoustically driven slots, external acoustic forcing, and mechanical piston methods. The effect of a planar mechanical piston forcing on a single stream shear layer is presented; it can be noted that this is one of the lesser studied free shear layers. The single stream shear layer can be characterized by its primary flow velocity scale and the thickness of the separating boundary layer. The velocity scale is constant over the length of the flow field; theta (x) can be used as a width scale to characterize the unforced shear layer. In the case of the forced shear layer the velocity field is a function of phase time and definition of a width measure becomes somewhat problematic.

Haw, Richard C.↗

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗

Homogeneity of spatial correlation statistics of tropical oceanic rainfall

The possibility of uniform horizontal correlation scales for tropical oceanic rainfall has been examined by a study of satellite-observed microwave data as a proxy measure of rain rates. From the brightness temperatures from the electrically scanning microwave radiometer on Nimbus 5 near nadir during the year 1974, the mean spatial autocorrelation function as a function of simultaneous pixel separation was calculated in each 5 by 5 deg grid box over the tropical Pacific and Atlantic for each season. The equal-time spatially lagged correlations were compared for geographical dependence to investigate the hypothesis of homogeneity. A simple model of the spatial statistics of the microwave brightness temperatures was used, consisting of a mixture of uncorrelated spatial white noise incoherently superimposed on a spatially correlated field (spatial red noise). The red noise signals are presumed to be generated by convective activity in the tropical atmosphere. The parameters of the red noise are consistent with this scheme over the tropical oceans, yielding a uniform spatial scale of about 50 km throughout.

Shin, Kyung-Sup↗

A Rapid Assessment of a Lunar Surface Laboratory Module

Johnson Space Center has piloted an innovation team known as the Forge to perform rapid turnaround, concurrent engineering studies. The Forge is intended to be able to pull together a team of strategically selected domain experts to innovatively solve a specified problem within a short period of time. The Forge’s inaugural project was the development of a Lunar Surface Laboratory Module concept. The lab module is intended as an option for the NASA Lunar Architecture Team and may allow for increased science productivity and return on investment than with only the Surface Habitat. The Forge team developed two concepts, a stationary laboratory and a mobile laboratory. Human factors evaluations have been conducted to determine the acceptability of the resulting design concepts. The human factors evaluations conducted include a card sort evaluation of the laboratories by science discipline and the System Usability Scale (SUS). The card sort focused subject matter expert expectations for the relative positioning of the laboratory functional areas, determining whether a given lab area should be adjacent to, or separated from, another. The SUS focused on layout efficiency. Both used human factors personnel and both life and physical scientists as subject matter experts. Each of these evaluations will focus primarily on objective data with limited subjective comments, reducing the time required for processing responses. The results of the card sort and SUS evaluations and their implications will be discussed. Based on the evaluation data, an overall acceptability of the Lunar Surface Laboratory Module produced by the Forge study can be determined. Lessons learned and recommendations for future use of human factors rapid assessment processes in early design trades will also be provided.

Forge↗

Conditions in the WR 140 Wind-Collision Region Revealed by the 1.083-mu m Hei Line Profile

We present spectroscopy of the P Cygni profile of the 1.083-μm He I line in the WC7+O5 colliding-wind binary (CWB) WR 140 (HD 193793), observed in 2008, before its periastron passage in 2009, and in 2016–2017, spanning the subsequent periastron passage. Both absorption and emission components showed strong variations. The variation of the absorption component as the O5 star was occulted by the wind-collision region (WCR) sets a tight constraint on its geometry. While the sightline to the O5 star traversed the WCR, the strength and breadth of the absorption component varied significantly on time-scales of days. An emission subpeak was observed on all our profiles. The variation of its radial velocity with orbital phase was shown to be consistent with formation in the WCR as it swung round the stars in their orbit. Modelling the profile gives a measure of the extent of the subpeak-forming region. In the phase range 0.93–0.99, the flux in the subpeak increased steadily, approximately inversely proportionally to the stellar separation, indicating that the shocked gas in the WCR where the line was formed was adiabatic. After periastron, the subpeak flux was anomalously strong and varied rapidly, suggesting formation in clumps downstream in the WCR. For most of the time, its flux exceeded the 2–10-keV X-ray emission, showing it to be a significant coolant of the shocked wind.

Peredur M Williams↗

Spatio–Temporal Machine Learning for Regional to Continental Scale Terrestrial Hydrology

Integrated hydrologic models can simulate coupled surface and subsurface processes but are computationally expensive to run at high resolutions over large domains. Here we develop a novel deep learning model to emulate subsurface flows simulated by the integrated ParFlow–CLM model across the contiguous US. We compare convolutional neural networks like ResNet and UNet run autoregressively against our novel architecture called the Forced SpatioTemporal RNN (FSTR). The FSTR model incorporates separate encoding of initial conditions, static parameters, and meteorological forcings, which are fused in a recurrent loop to produce spatiotemporal predictions of groundwater. We evaluate the model architectures on their ability to reproduce 4D pressure heads, water table depths, and surface soil moisture over the contiguous US at 1 km resolution and daily time steps over the course of a full water year. The FSTR model shows superior performance to the baseline models, producing stable simulations that capture both seasonal and event–scale dynamics across a wide array of hydroclimatic regimes. The emulators provide over 1,000× speedup compared to the original physical model, which will enable new capabilities like uncertainty quantification and data assimilation for integrated hydrologic modeling that were not previously possible. Our results demonstrate the promise of using specialized deep learning architectures like FSTR for emulating complex process–based models without sacrificing fidelity.

54 ENVIRONMENTAL SCIENCES↗

Numerical studies of major and minor stratospheric warmings caused by orographic forcing

A primitive equation spectral model using spherical harmonics is formulated to study dynamic interactions between the troposphere and stratosphere in association with sudden stratospheric warmings. Using sigma coordinates for five tropospheric layers and log-pressure coordinates for 26 stratospheric and mesospheric layers, separate model equations for each system are combined to form single matrix governing equations. The gradual introduction of large scale topography to balanced initial states representative of observed mean winter conditions in the Northern Hemisphere is used for the generation of planetary waves during 40-day time integrations. Results of these integrations indicate that stratospheric warmings can be simulated by this orographic forcing and that mean momentum flux divergence due to zonal mean motion appears to be an essential mechanism of these simulated sudden warmings. It was found that the strength of the polar night jet can be a determining factor whether a warming becomes 'major' or 'minor'.

Koermer, J. P.↗

Surface heating and patchiness in the coastal ocean off central California during a wind relaxation event

The difference between the temperature of the ocean at 4-cm and 2-m depth was continuously monitored during a cruise to the coastal transition zone off Point Arena, California, during June 1987. The two temperatures were coincident most of the time but diverged during one nearshore leg of the cruise where large temperature differences of up to 4.7 C were observed between the 4-cm and 2-m sensors, in areas which were separated by regions where the two temperatures were coincident as usual. The spatial scale of this 'patchy' thermal structure was about 5-10 km. A mixed layer model (Garwood, 1977) was used to simulate the near surface stratification when forced by the observed wind stress, surface heating, and optical clarity of the water. The model produced a thin strongly stratified surface layer at stations where exceptionally high turbidity was observed but did not produce such features otherwise. This simple model could not explain the horizontal patchiness in the thermal structure, which was likely due to patchiness in the near-surface chlorophyll distributions or to submesoscale variability of the surface wind stress.

Ramp, Steven R.↗

On the Computation of Sound by Large-Eddy Simulations

The effect of the small scales on the source term in Lighthill's acoustic analogy is investigated, with the objective of determining the accuracy of large-eddy simulations when applied to studies of flow-generated sound. The distribution of the turbulent quadrupole is predicted accurately, if models that take into account the trace of the SGS stresses are used. Its spatial distribution is also correct, indicating that the low-wave-number (or frequency) part of the sound spectrum can be predicted well by LES. Filtering, however, removes the small-scale fluctuations that contribute significantly to the higher derivatives in space and time of Lighthill's stress tensor T(sub ij). The rms fluctuations of the filtered derivatives are substantially lower than those of the unfiltered quantities. The small scales, however, are not strongly correlated, and are not expected to contribute significantly to the far-field sound; separate modeling of the subgrid-scale density fluctuations might, however, be required in some configurations.

Piomelli, Ugo↗

Better, Cheaper, Faster Molecular Dynamics

Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.

Pohorille, Andrew↗

The Sun's Impact on Climate

We provide an overview of the impact of the Sun on the Earth atmosphere and climate system, focused on heating of Earth's atmosphere and oceans. We emphasize the importance of the spectral measurements of SIM and SOLSTICE- that we must know how solar variations are distributed over ultraviolet, visible, and infrared wavelengths, since these have separate characteristic influences on Earth's ozone layer, clouds, and upper layers of the oceans. Emphasis is also given to understanding both direct and indirect influences of the Sun on the Earth, which involve feedbacks between Earth's stratosphere, troposphere, and oceans, each with unique time scales, dynamics, chemistry, and biology, interacting non-linearly. Especially crucial is the role of all three phases of water on Earth, water vapor being the primary greenhouse gas in the atmosphere, the importance of trace gases such as CO2 arising from their absorption in the "water vapor window" at 800 - 1250/cm (12.5 to 8 microns). Melting of polar ice is one major response to the post-industrial global warming, enhanced due to "ice-albedo" feedback. Finally, water in liquid form has a major influence due to cloud albedo feedback, and also due to the oceans' absorption of solar radiation, particularly at visible wavelengths, through the visible "liquid water window" that allows penetration of visible light deep into the mixed layer, while nearby ultraviolet and infrared wavelengths do not penetrate past the upper centimeter ocean surface skin layer. A large fraction of solar energy absorbed by the oceans goes into the latent heat of evaporation. Thus the solar heating of the atmosphere-ocean system is strongly coupled through the water cycle of evaporation, cloud formation, precipitation, surface runoff and ice formation, to Earth's energy budget and climate, each different climate component responding to variations in different solar spectral bands, at ultraviolet, visible and infrared wavelengths.

Cahalan, Robert↗

An Experimental Investigation of the Flow Over the Rear End of a Notchback Automobile Configuration

An experimental investigation of the flow over the rear end of a 0.16 scale notchback automobile configuration has been conducted in the NASA Langley Basic Aerodynamics Research Tunnel (BART). The objective of this work was to investigate the flow separation that occurs behind the backlight and obtain experimental data that can be used to understand the physics and time-averaged structure of the flow field. A three-component laser velocimeter was used to make non-intrusive, velocity measurements in the center plane and in a single cross-flow plane over the decklid. In addition to off-body measurements, flow conditions on the car surface were documented via surface flow visualization, boundary layer measurements, and surface pressures. The experimental data show several features previously identified by other researchers, but also reveal differences between the flow field associated with this particular configuration and the generally accepted models for the flow over a notchback rear end.

Jenkins, Luther N.↗

Centrifugal Sieve for Gravity-Level-Independent Size Segregation of Granular Materials

Conventional size segregation or screening in batch mode, using stacked vibrated screens, is often a time-consuming process. Utilization of centrifugal force instead of gravity as the primary body force can significantly shorten the time to segregate feedstock into a set of different-sized fractions. Likewise, under reduced gravity or microgravity, a centrifugal sieve system would function as well as it does terrestrially. When vibratory and mechanical blade sieving screens designed for terrestrial conditions were tested under lunar gravity conditions, they did not function well. The centrifugal sieving design of this technology overcomes the issues that prevented sieves designed for terrestrial conditions from functioning under reduced gravity. These sieves feature a rotating outer (cylindrical or conical) screen wall, rotating fast enough for the centrifugal forces near the wall to hold granular material against the rotating screen. Conventional centrifugal sieves have a stationary screen and rapidly rotating blades that shear the granular solid near the stationary screen, and effect the sieving process assisted by the airflow inside the unit. The centrifugal sieves of this new design may (or may not) have an inner blade or blades, moving relative to the rotating wall screen. Some continuous flow embodiments would have no inner auger or blades, but achieve axial motion through vibration. In all cases, the shearing action is gentler than conventional centrifugal sieves, which have very high velocity differences between the stationary outer screen and the rapidly rotating blades. The new design does not depend on airflow in the sieving unit, so it will function just as well in vacuum as in air. One advantage of the innovation for batch sieving is that a batch-mode centrifugal sieve may accomplish the same sieving operation in much less time than a conventional stacked set of vibrated screens (which utilize gravity as the primary driving force for size separation). In continuous mode, the centrifugal sieves can provide steady streams of fine and coarse material separated from a mixed feedstock flow stream. The centrifugal sieves can be scaled to any desired size and/or mass flow rate. Thus, they could be made in sizes suitable for small robotic exploratory missions, or for semi-permanent processing of regolith for extraction of volatiles of minerals. An advantage of the continuous-mode system is that it can be made with absolutely no gravity flow components for feeding material into, or for extracting the separated size streams from, the centrifugal sieve. Thus, the system is capable of functioning in a true microgravity environment. Another advantage of the continuous-mode system is that some embodiments of the innovation have no internal blades or vanes, and thus, can be designed to handle a very wide range of feedstock sizes, including occasional very large oversized pieces, without jamming or seizing up.

Walton, Otis R.↗

Quantitative Approach for the Accurate CFD Simulation of Hover in Turbulent Flow

Time-dependent Navier-Stokes simulations have been carried out for a V22 rotor in hover using an improved, low-diffusion, HLLE++ upwind algorithm in the OVERFLOW code. Emphasis is placed on lessons learned over the past decade regarding the effects of high-order spatial accuracy, grid resolution, and the use of detached eddy simulation in predicting the rotor performance, i.e., figure-of-merit (FM). A quick-start procedure is described together with a statistical measure of FM convergence that reduces hover computations by 5 times, similar in computational work for forward flight. Furthermore, Cartesian adaptive mesh refinement in the rotor wake revealed a complex turbulent flow with worm-like structures of various scales. These turbulent worms, found computationally more than a decade ago, have recently been observed in a separate German Aerospace Center (DLR) experiment. Moreover, adaptive mesh refinement has been used to resolve the tip-vortex to its correct physical size.

Neal M. Chaderjian↗

Performance of vortex generators in a Mach 2.5 low-bleed full scale 45-percent-internal-contraction axisymmetric inlet

Steady-state and dynamic flow characteristics associated with two sets of vortex generators having different mixing criteria were determined. The inlet performance with and without these vortex generators is presented. The vortex generators were successful in eliminating separation, increasing area-weighted total pressure recovery, and decreasing distortion. Transmission times obtained from cross-correlations of the wall static pressures and the diffuser exit total pressure showed no effect of the upstream flow characteristics on the diffuser exit pressures when generators were used. Without generators, separation occurred and the upstream pressure characteristics had immediate effects on the diffuser exit pressure characteristics.

Neumann, H. E.↗