Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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 1,783 records · Page 99

Tropical Ocean Global Atmosphere (TOGA) Meteorological and Oceanographic Data Sets for 1985 and 1986

The Tropical Ocean Global Atmosphere (TOGA) Program is a component of the World Meteorological Organization (WMO)/International Council of Scientific Unions (ICSU) World Climate Research Program (WCRP). One of the objectives of TOGA, which began in 1985, is to determine the limits of predictability of monthly mean sea surface temperature variations in tropical regions. The TOGA program created a raison d'etre for an explosive growth of the tropical ocean observing system and a substantial improvement in numerical simulations from atmospheric and oceanic general circulation models. Institutions located throughout the world are involved in the TOGA-distributed active data archive system. The diverse TOGA data sets for 1985 and 1986, including results from general circulation models, are included on a CD-ROM. Variables on the CD-ROM are barometric pressure, surface air temperature, dewpoint temperature Cartesian components of surface wind, surface sensible and latent heat fluxes,Cartesian components of surface wind stress and of an index of surface wind stress, sea level, sea surface temperature, and depth profiles of temperature and current in the upper ocean. Some data sets are global in extent, some are regional and cover portions of an ocean basin. Data on the CD-ROM can be extracted with an Apple Macintosh or an IBM PC.

Halpern, D.↗

Efficient Agent-Based Cluster Ensembles

Numerous domains ranging from distributed data acquisition to knowledge reuse need to solve the cluster ensemble problem of combining multiple clusterings into a single unified clustering. Unfortunately current non-agent-based cluster combining methods do not work in a distributed environment, are not robust to corrupted clusterings and require centralized access to all original clusterings. Overcoming these issues will allow cluster ensembles to be used in fundamentally distributed and failure-prone domains such as data acquisition from satellite constellations, in addition to domains demanding confidentiality such as combining clusterings of user profiles. This paper proposes an efficient, distributed, agent-based clustering ensemble method that addresses these issues. In this approach each agent is assigned a small subset of the data and votes on which final cluster its data points should belong to. The final clustering is then evaluated by a global utility, computed in a distributed way. This clustering is also evaluated using an agent-specific utility that is shown to be easier for the agents to maximize. Results show that agents using the agent-specific utility can achieve better performance than traditional non-agent based methods and are effective even when up to 50% of the agents fail.

Agogino, Adrian↗

Three‐Dimensional Magnetic Reconnection Spreading in Current Sheets of Non‐Uniform Thickness

Abstract Magnetic reconnection in naturally occurring and laboratory settings often begins locally and elongates, or spreads, in the direction perpendicular to the plane of reconnection. Previous work has largely focused on current sheets with a uniform thickness, for which the predicted spreading speed for anti‐parallel reconnection is the local speed of the current carriers. We derive a scaling theory of three‐dimensional (3D) spreading of collisionless anti‐parallel reconnection in a current sheet with its thickness varying in the out‐of‐plane direction, both for spreading from a thinner to thicker region and a thicker to thinner region. We derive an expression for calculating the time it takes for spreading to occur for a current sheet with a given profile of its thickness. A key result is that when reconnection spreads from a thinner to a thicker region, the spreading speed in the thicker region is slower than both the Alfvén speed and the speed of the local current carriers by a factor of the ratio of thin to thick current sheet thicknesses. This is important because magnetospheric and solar observations have previously measured the spreading speed to be slower than previously predicted, so the present mechanism might explain this feature. We confirm the theory via a parametric study using 3D two‐fluid numerical simulations. We use the prediction to calculate the time scale for reconnection spreading in Earth's magnetotail during geomagnetic activity. The results are also potentially important for understanding reconnection spreading in solar flares and the dayside magnetopause of Earth and other planets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear dimension reduction for surrogate modeling using gradient information

We introduce a method for the nonlinear dimension reduction of a high-dimensional function $u:{\mathbb{R}}^d\rightarrow{\mathbb{R}}$, $d\gg 1$. Our objective is to identify a nonlinear feature map $g:{\mathbb{R}}^d\rightarrow{\mathbb{R}}^m$, with a prescribed intermediate dimension $m\ll d$, so that $u$ can be well approximated by $f\circ g$ for some profile function $f:{\mathbb{R}}^m\rightarrow{\mathbb{R}}$. We propose to build the feature map by aligning the Jacobian $\nabla g$ with the gradient $\nabla u$, and we theoretically analyze the properties of the resulting $g$. Once $g$ is built, we construct $f$ by solving a gradient-enhanced least squares problem. Our practical algorithm uses a sample $\{{\textbf{x}}^{(i)},u({\textbf{x}}^{(i)}),\nabla u({\textbf{x}}^{(i)})\}_{i=1}^N$ and builds both $g$ and $f$ on adaptive downward-closed polynomial spaces, using cross validation to avoid overfitting. We numerically evaluate the performance of our algorithm across different benchmarks, and explore the impact of the intermediate dimension $m$. We show that building a nonlinear feature map $g$ can permit more accurate approximation of $u$ than a linear $g$, for the same input data set.

97 MATHEMATICS AND COMPUTING↗

A Case For Intra-rack Resource Disaggregation in HPC

The expected halt of traditional technology scaling is motivating increased heterogeneity in high-performance computing (HPC) systems with the emergence of numerous specialized accelerators. As heterogeneity increases, so does the risk of underutilizing expensive hardware resources if we preserve today’s rigid node configuration and reservation strategies. This has sparked interest in resource disaggregation to enable finer-grain allocation of hardware resources to applications. However, there is currently no data-driven study of what range of disaggregation is appropriate in HPC. To that end, we perform a detailed analysis of key metrics sampled in NERSC’s Cori, a production HPC system that executes a diverse open-science HPC workload. In addition, we profile a variety of deep-learning applications to represent an emerging workload. We show that for a rack (cabinet) configuration and applications similar to Cori, a central processing unit with intra-rack disaggregation has a 99.5% probability to find all resources it requires inside its rack. In addition, ideal intra-rack resource disaggregation in Cori could reduce memory and NIC resources by 5.36% to 69.01% and still satisfy the worst-case average rack utilization.

97 MATHEMATICS AND COMPUTING↗

Transonic laminar boundary layers with surface curvature.

The effect of surface curvature (both longitudinal and transverse) and the associated pressure gradient across the flow is investigated analytically for a laminar boundary layer subjected to pressure gradients along the flow. Property variation which results from heat transfer and compressibility is taken into account. Numerical solutions of the boundary layer equations are obtained for locally similar sonic flow through the throat of a nozzle for a range of flow conditions and for various shaped nozzle surfaces with different amounts of wall cooling. A few solutions were also obtained for the analogous flow around the shoulder of a flat-faced body in a supersonic flow. The effect of various parameters that arise in the equations upon application of the Levy-Mangler transformation are investigated and discussed with respect to their influence on the velocity and total enthalpy profiles and the corresponding profile slopes at the surface to which the shear stress and heat transfer are related. An important finding is that at throat Reynolds numbers less than 100,000 the heat transfer parameter at a nozzle throat decreases as the throat radius of curvature decreases.

Back, L. H.↗

The influence of gravity wave breaking on the general circulation of the middle atmosphere

The zonal mean solstice circulation of the global middle atmosphere is simulated using a semi-spectral numerical model. Radiative heating and cooling is computed by the algorithm of Wehrbein and Leovy. Mechanical dissipation is represented by the gravity wave breaking parameterization of Lindzen. An inertial adjustment parameterization is used to prevent the development of inertially unstable meridional shears near the equator. It is shown that gravity wave drag and diffusion in the mesosphere can account for the observed large departure from radiative equilibrium in both summer and winter. Experiments incorporating a forced stationary wavenumber 1 disturbance indicate that planetary wave EP flux convergences, although they may modify the mean flow profile significantly, cannot provide the major source of mechanical dissipation in the winter mesosphere. A simulated sudden warming is accompanied by an equally strong mesospheric cooling. This cooling is caused primarily by the relaxation of the polar mesosphere toward radiative equilibrium when the easterly mean winds in the polar stratosphere induced by the sudden warming reduce the transmission of gravity waves into the mesosphere.

Holton, J. R.↗

Evidence for divergent plate-boundary characteristics and crustal spreading on Venus

Detailed examination of the topography and morphology of western Aphrodite Terra reveals numerous features that are similar to terrestrial divergent plate-boundary characteristics. Individual domains between fracture-zone-like discontinuities contain a variety of bilaterally symmetrical topographic elements that suggest that topographic features have been created at rise crests, rifted and separated, and moved laterally to their present symmetrical positions. The topographic and morphologic similarities, together with strikingly similar mirror-image map patterns on both sides of the rise axis, suggest that western Aphrodite Terra shares the characteristics of oceanic divergent plate boundaries, and is the site of crustal spreading on Venus. Topographic profiles are consistent with spreading rates of the order of several centimeters per year.

Head, James W., III↗

Pulse Detonation Engine Modeled

Pulse Detonation Engine Technology is currently being investigated at Glenn for both airbreathing and rocket propulsion applications. The potential for both mechanical simplicity and high efficiency due to the inherent near-constant-volume combustion process, may make Pulse Detonation Engines (PDE's) well suited for a number of mission profiles. Assessment of PDE cycles requires a simulation capability that is both fast and accurate. It should capture the essential physics of the system, yet run at speeds that allow parametric analysis. A quasi-one-dimensional, computational-fluid-dynamics-based simulation has been developed that may meet these requirements. The Euler equations of mass, momentum, and energy have been used along with a single reactive species transport equation, and submodels to account for dominant loss mechanisms (e.g., viscous losses, heat transfer, and valving) to successfully simulate PDE cycles. A high-resolution numerical integration scheme was chosen to capture the discontinuities associated with detonation, and robust boundary condition procedures were incorporated to accommodate flow reversals that may arise during a given cycle. The accompanying graphs compare experimentally measured and computed performance over a range of operating conditions for a particular PDE. Experimental data were supplied by Fred Schauer and Jeff Stutrud from the Air Force Research Laboratory at Wright-Patterson AFB and by Royce Bradley from Innovative Scientific Solutions, Inc. The left graph shows thrust and specific impulse, Isp, as functions of equivalence ratio for a PDE cycle in which the tube is completely filled with a detonable hydrogen/air mixture. The right graph shows thrust and specific impulse as functions of the fraction of the tube that is filled with a stoichiometric mixture of hydrogen and air. For both figures, the operating frequency was 16 Hz. The agreement between measured and computed values is quite good, both in terms of trend and magnitude. The error is under 10 percent everywhere except for the thrust value at an equivalence ratio of 0.8 in the left figure, where it is 14 percent. The simulation results shown were made using 200 numerical cells. Each cycle of the engine, approximately 0.06 sec, required 2.0 min of CPU time on a Sun Ultra2. The simulation is currently being used to analyze existing experiments, design new experiments, and predict performance in propulsion concepts where the PDE is a component (e.g., hybrid engines and combined cycles).

Paxson, Daniel E.↗

Numerical investigation of alpha particle confinement under the perturbation of neoclassical tearing modes and toroidal field ripple in CFETR

The confinement of alpha particles in burning plasma is a key issue in fusion reactor design, including particle interaction with instabilities. Here, we include two topics: the effect of neoclassical tearing modes (NTMs) and toroidal field ripple on alpha particle loss, and the assessment of particle redistribution under an NTM with a reduced model. We consider Chinese fusion engineering test reactor parameters, the alpha particle distribution given by TRANSP/NUBEAM and the NTM perturbation function given by the initial value code TM1. We show that the synergistic effect of the NTM and ripple is negligible; the particle loss fraction does not change with increasing NTM amplitude. Only passing particles are affected by the mode particle resonance, producing profile flattening but no increased loss because only trapped particles are influenced by ripple. To study alpha particle profile flattening, the work adopts an innovative method of phase vector rotation to determine regions of good and broken Kolmogorov–Arnold–Moser surfaces and equilibrates the particle density according to local stochasticity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reducing Uncertainty in Offshore Wind Energy Yield Estimates via a Metocean Reference Site

The offshore wind industry is burgeoning in the coastal waters of the United States, specifically along the Atlantic. For wind energy to be successful, reliable observations and model simulations are needed for resource assessment and forecasting. While many of these activities have already begun, there is currently an absence of observations at hub-height in these waters, with the closest available hub-height measurements usually taken onshore. Deployment of floating lidars has occurred through various federally funded projects, but only encapsulates time periods of a couple of years at best. Private industry is also beginning to leverage floating lidars, but this data is often proprietary, and not shared with the general public. In this work, we make the case for a metocean reference site for long-term offshore wind energy. Specifically, we quantify the impact of having a metocean reference site compared to other methods of determining hub-height winds and energy production. We use an offshore floating lidar to directly measure the wind resource, and compare these measurements to predictions derived from other widely-available surface meteorological variables. These prediction methods (vertical extrapolation, machine learning, and NWP output) produce a variety of vertical wind speed profiles, of which produce different energy yield estimates for a reference offshore turbine (Figure 1). While some methods perform reasonably well against the lidar, the uncertainty in these energy yield estimates has financial implications, further illustrating the need for long-term measurements in coastal waters.

machine learning↗

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗

Composition of the nighttime ionospheric F 1 region near the magnetic equator

The effects of vertical E x B transport on NO(+), O2(+) and O(+) densities in the nighttime equatorial ionospheric F 1 region are investigated. Ion densities are calculated as functions of altitude, latitude and local time by the numerical solution of coupled, time-dependent ion continuity equations, taking into account production, loss by charge exchange and dissociative recombination and transport by diffusion and E x B drift. The results of the calculations are compared with measurements of NO(+), O2(+) and O(+) ion densities obtained at low altitudes by a mass spectrometer on board the Atmospheric Explorer C satellite, and are found to be consistent with the observations, suggesting that in the equatorial region, vertical transport by E x B drift is primarily responsible for producing the observed NO(+), O2(+) and O(+) density profiles. In addition, the reaction of O2(+) with N(4S) is found to be an important sink for O2(+) and a source of NO(+) ions. Implications of the observed and calculated near constancy of electron and ion densities with altitude when NO(+) is the dominant ion on the growth of large-scale irregularities are also considered

Anderson, D. N.↗

Development and Integration of a Thermal Management Simulation for a Quadrotor Parallel Hybrid Propulsion System

This paper details the development of a propulsion system simulation for a six-passenger parallel hybrid quadrotor and utilizes the Numerical Propulsion System Simulation (NPSS) along with the NPSS Power System Library as the development environment. This simulation integrates an engine power plant with an electrical generation and distribution system and includes the required thermal management system. The thermal management system is comprised of liquid cooling loops that reject the heat load through air to coolant heat exchangers and utilizes a map-based performance estimation method. This method is developed within NPSS and detailed in this paper. The full system model is designed to predict system weight, range, and performance through a proposed mission profile. Results of the paper show an all-engine system maintains the best range, while a mostly electric system that utilizes an engine as a backup or a conditional power contributor offers range benefit.

Vertical lift and take off vehicle↗

Development and Integration of a Thermal Management Simulation for a Quadrotor Parallel Hybrid Propulsion System

This paper details the development of a propulsion system simulation for a six-passenger parallel hybrid quadrotor and utilizes the Numerical Propulsion System Simulation (NPSS) along with the NPSS Power System Library as the development environment. This simulation integrates an engine power plant with an electrical generation and distribution system and includes the required thermal management system. The thermal management system is comprised of liquid cooling loops that reject the heat load through air to coolant heat exchangers and utilizes a map-based performance estimation method. This method is developed within NPSS and detailed in this paper. The full system model is designed to predict system weight, range, and performance through a proposed mission profile. Results of the paper show an all-engine system maintains the best range, while a mostly electric system that utilizes an engine as a backup or a conditional power contributor offers range benefit.

Vertical lift and take off vehicle↗

Observation of Deep Convective Cloud-Top Height and Vertical Temperature Structure of Hurricane Using Hyperspectral Infrared Sounder and its Single-Field-View Retrieval Products

Hurricanes, severe tropical cyclones (TC), or typhoons are significant natural disasters that often result in substantial loss of life and property damage. Numerous studies have indicated that changes in TC intensity are closely linked to deep convective clouds (DCC), with stronger TCs typically exhibiting higher cloud top heights (CTH) compared to weaker TCs. The CTH can help determine if a tropical depression is at the onset of rapid intensification based on case studies. Therefore, accurate determination of TC CTH will be greatly helpful for monitoring TC development and studying TC dynamics. One traditional and most common method to derive CHT from satellite observations is using the thermal brightness temperature in atmospheric channels to match the sounding temperature profile. However, it was found that thermally derived CTH has a lower bias of approximately 1 km, and this bias tends to worsen for the tallest clouds. A new method using the hyperspectral infrared sounder will be presented. From the measurements of Cross-track Infrared Sounder (CrIS) on S-NPP and J-1, along with radiative transfer simulations, we identified the inverted-V spectral feature in the ozone (O3) band (near 9.6 μm) corresponding to high clouds. The depth of the inverted-V can be used to estimate the CTH. Since the depth is computed using the peak absorption O3 channel and the nearby most transparent O3 channel in this O3 band, the uncertainties associated with cloud emissivity and scattering by cloud particles in the traditional method can be ignored. From several hurricane case studies, we found that the CHT derived using this method can accurately capture the structure of the cloud tops in the eyewall, spiral rainbands, and surrounding regions. For example, Hurricane Dorian on September 2, 2019, showed a nicely outward-sloping and circular shape eye cloud in the early morning, but the circular shape of the eyewall cloud became distorted in the afternoon. For various hurricanes we examined, the distribution of CHT for the eyewall clouds differed significantly. To better study the thermodynamic structure of hurricane clouds, this research will analyze the vertical temperature profiles from a new single Field of View (SFOV) Sounder Atmospheric Products (SiFSAP), derived using CrIS and the Advanced Technology Microwave Sounder (ATMS) onboard SNPP and JPSS-1. SiFSAP has a spatial resolution of 15 km at nadir, which surpasses most global weather and climate models and other current operational sounding products. The combined use of ATMS and CrIS allows for retrievals near hurricane eyewalls and spiral rainbands. Wind fields from NASA’s Modern-Era Retrospective Analysis for Research and Applications Version-2 (MERRA-2) and ERA5 will be used to characterize transport, and comparisons between the model temperature and water vapor profiles with the corresponding SiFSAP products will also be provided.

SiFSAP↗

Impact of a Bluff Body on a Pressurized Oxy-fuel Flame

Designing an effective burner is vital for the development of pressurized oxy-fuel combustion technologies. In the present work, turbulent jet diffusion burners are adopted for a pressurized oxy-combustor, with a bluff-body employed to create a recirculation zone, thereby stabilizing the flame in such a combustor. The objective of this numerical study is to perform a systematic analysis of the characteristics of such a pressurized non-premixed flame. Specifically, a 15-bar pressurized oxy-fuel combustor of power 100 kWth is modeled by means of the Ansys FLUENT commercial platform, using the Reynolds-averaged Navier-Stokes (RANS) approach. The present work focuses on identifying the aerothermodynamic features of the pressurized oxy-fuel burner with a disk-shape bluff-body. It is shown that the fuel-to-oxidizer stream momentum ratio has a great impact on the temperature profile of the down-fired, co-axial, pressurized oxy-coal, diffusion flame. A parametric study of the blockage ratio of the burner identifies an optimal range for this pressurized burner.

Li, Lei↗

Numerical simulation of high Reynolds number transonic flows

A code has been developed for simulating high Reynolds number transonic flow fields of arbitrary configuration. An explicit finite-difference method with time splitting is used to solve the time-dependent equations for compressible turbulent flow. A nonorthogonal computational mesh of arbitrary configuration facilitates the description of the flow field. The code is applied to simulate the flow over an 18 percent thick circular-arc biconvex airfoil at zero angle of attack and free-stream Mach number of 0.775. A simple mixing-length model is used to describe the turbulence and chord Reynolds numbers of 1, 2, 4, and 10 million are considered. The solution describes in sufficient detail both the shock-induced and trailing-edge separation regions, and provides the profile and friction drag.

Deiwert, G. S.↗