Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “graph analysis”

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 703 records · Page 39

Development of a calibrated software reliability model for flight and supporting ground software for avionic systems

The object of this project was to develop and calibrate quantitative models for predicting the quality of software. Reliable flight and supporting ground software is a highly important factor in the successful operation of the space shuttle program. The models used in the present study consisted of SMERFS (Statistical Modeling and Estimation of Reliability Functions for Software). There are ten models in SMERFS. For a first run, the results obtained in modeling the cumulative number of failures versus execution time showed fairly good results for our data. Plots of cumulative software failures versus calendar weeks were made and the model results were compared with the historical data on the same graph. If the model agrees with actual historical behavior for a set of data then there is confidence in future predictions for this data. Considering the quality of the data, the models have given some significant results, even at this early stage. With better care in data collection, data analysis, recording of the fixing of failures and CPU execution times, the models should prove extremely helpful in making predictions regarding the future pattern of failures, including an estimate of the number of errors remaining in the software and the additional testing time required for the software quality to reach acceptable levels. It appears that there is no one 'best' model for all cases. It is for this reason that the aim of this project was to test several models. One of the recommendations resulting from this study is that great care must be taken in the collection of data. When using a model, the data should satisfy the model assumptions.

Lawrence, Stella↗

Effects of Heat Transfer Coefficient Variation on Nuclear Thermal Propulsion Engine Performance

A physics-based Nuclear Thermal Propulsion (NTP) Testing Reference Design (TRD) power balance model was coded in Simulink to investigate engine performance for various design and parameter modifications. Since the primary mode of heat transfer in NTP engines is convective, the convective heat transfer coefficient (HTC) is a key parameter that requires accurate representation. The industry standard Westinghouse correlation has an uncertainty of ±20% which was investigated in this study. The results showed that a 20% decrease in the HTC led to a 4.14% increase in maximum fuel temperature while a 20% decrease in the HTC led to a 1.81% decrease in maximum fuel temperature suggesting that narrowing the uncertainty of this correlation through experimental work would be a critical step in the development of NTP engines. Furthermore, a maximum fuel temperature relationship with specific impulse was developed for the TRD engine which showed potential engine operation between specific impulse values of 715 and 900 seconds with minimal changes to the engine design. This graph could be useful for high level vehicle performance estimations for fuel types with different maximum operating temperatures.

Heat Transfer Coefficient↗

Historical Pattern Analysis of Global Geothermal Power Capacity Development

Between 1913 and 1958, Italy was the only country with an operational geothermal power plant until New Zealand installed its first plant in 1958. At present, 24 countries are involved in the geothermal power market, and they have a combined installed capacity of 16,127 GW. This study analyzes the historical patterns of geothermal power capacity in the world and in individual countries to investigate the ideal global geothermal development pattern by examining the annual cumulative capacity (ACC) and the annual capacity addition (ACA) graphs of the historical development of geothermal power capacity in 24 countries. First, the global patterns are analyzed using these graphs in five periods (1945-1957, 1958-1976, 1977-1991, 1992-2002, and 2003-2020) that are marked by a series of characteristics of ACA peaks separated by two major troughs. Then, five characteristic patterns are developed in five periods globally. These patterns correspond to the early-stage linear, the first acceleration, the first steady-state linear, the second acceleration, and the second steady-state linear developments. A positive relationship exists between global patterns and the 5-year shifted oil-price curve: two major factors influenced global development: 1) increasing oil prices and increasing awareness of global climate change, and 2) global development of geothermal power. Last, we investigate these patterns in each country. The top ten countries, which comprise 93.3% of the world's total installed capacity are separated into five groups based on the availability and characteristics of patterns globally developed in five periods. Group-1 (the United States) has an installed capacity of 3,794 MW, Group-2 (Mexico and Philippines) 963-1935 MW, Group-3 (New Zealand, Italy, Iceland, and Japan) 601-1,037 MW, and Group-4 (Indonesia, Kenya, and Turkiye) 944-2,356 MW. The remaining 14 countries (6.7%), which are called Group 5, are still in an immature stage and have installed capacities of 7-262 MW and are not involved in pattern analysis. Overall, geothermal power in the world is in its third stage of development, which had its peak development after 1977. A fourth development peak may be expected to occur after this through business-as-usual cases. The biggest barrier to the development of the global geothermal power market is the risk associated with exploration and drilling. If risk mitigation systems and funds are employed, the growth of geothermal power production projects could accelerate.

GEOTHERMAL ENERGY↗

An analysis of satellite observations of the O I EUV dayglow

Observations of the latitudinal variation of the 1356-A, 1304-A, and 989-A O I earth dayglow, obtained with the USAF STP78-1 satellite at solar maximum on March 21, 1979, are reported and compared with the predictions of thermospheric models and published laboratory measurements of electron-impact excitation cross sections. The results are presented in extensive tables and graphs and characterized in detail, with particular attention to the 1304-A solar flux and the 989-A branching ratio. The STP78-1 data (as well as simultaneous AE-E satellite measurements ) are shown to indicate exospheric temperatures that were 220-240-K higher, but varied significantly less with latitude, than those predicted by the MSIS-83 model (Hedin, 1983) for the day of the observations.

Link, R.↗

Evaluation of Methods for Causal Discovery in Hydrometeorological Systems

Understanding causal relations is of utmost importance in hydrology and climate research for systems identification, prediction, and understanding systems behavior in a changing climate. Traditionally, researchers in hydrometeorology attempted to study causal questions by conducting controlled experiments using numerical models. This approach, however, in most cases of interest provides uncertain results because the models are approximate representation of the natural system. An alternative approach that has recently drawn significant attention in several fields is to infer causal relations from purely observational data. It possesses several traits to its utility particularly in hydrometeorology due to the rapid accumulation of in situ and remotely sensed data records. The first objective of this study is to present a brief description of four causal discovery methods (Granger causality, Transfer Entropy, graph-based algorithms, and Convergent Cross Mapping) with special emphasis on the assumptions on which they are built. Second, using synthetic data generated from a hydrological model, we assess their performance in retrieving causal information taking into account sensitivity to sample size and presence of noise. Last, we use causal analysis to examine and formulate hypotheses on causal drivers of evapotranspiration in a shrubland region during summer and winter seasons. An interpretation of the hypotheses based on canopy seasonal dynamics and evapotranspiration processes is presented. It is hoped that the results presented here can be useful in guiding researchers studying hydrometeorological systems as to which causal method is most appropriate to the characteristics of the system under study.

54 ENVIRONMENTAL SCIENCES↗

Impedance-Aware Graph Convolutional Networks for Voltage Estimation in Active Distribution Networks

Voltage estimation plays a key role in ensuring the effective control and reliability of distribution networks. However, traditional machine learning methods often fail to capture the details of the distribution network’s topology. To overcome this challenge, graph convolutional networks (GCN) have emerged as an alternative. Graph convolutional networks inherently capture the topology of the grid, utilizing correlations to achieve precise voltage estimation. Other machine learning models and conventional GCNs fail to account for the distribution line characteristics found in the real world, limiting their effectiveness. This paper proposes an advanced variant of GCN called the Impedance-Aware Graph Convolutional Network (IA-GCN). The IA-GCN layer incorporates the magnitude of the impedance into the graph convolution mechanism, allowing it to capture topological nuances and provide valuable insights into node interrelationships by considering impedance as an intrinsic dimension. The performance of the IA-GCN layer is then compared with that of GCN and GraphSAGE layers through a surrogate model for voltage estimation. The performance analysis demonstrates that IA-GCN outperforms GCN by reducing the MAE by 87.55% and improving the R-squared value by 98%.

Ravi, Abhijith↗

Gas-core reactor power transient analysis.

The nuclear fuel in the gas-core reactor concept is a ball of uranium plasma radiating thermal photons. The photons are met by an inflowing hydrogen stream, which is seeded with submicron size, depleted uranium particles. A 'wall-burnout' condition exists if the thermal photons can reach the cavity liner because of insufficient absorption by the hydrogen. An analysis was conducted in order to determine the time for which the maximum steady state reactor power could be exceeded without damage to the cavity liner due to burnout. Wall-burnout time as a function of the power increase above the initial steady state condition is shown in a graph.

Kascak, A. F.↗

Research investigation of the physical interactions and phenomena associated with hypervelocity sub-micron particles

Signals from impact ionization plasmas were studied as a means of performing microparticle composition analysis. Impact ionization signal response was measured in a time-of-flight (TOF) system for lanthanum hexaboride, carbonyl iron, and aluminum microparticle impacts on a tantalum target, primarily in the 1 - 8 km/s velocity range. Oscilloscope photographs of representative ion TOF signal response are given for each material studied. Graphs and histograms are presented of the total charge collected as well as the charge collected in each observed ion mass group. Data show that ion signals consist primarily of the lower ionization potential elements over the 1 - 8 km/s range.

Roy, N. L.↗

Evaluation of flow topology from numerical data

Results obtained from numerical calculations and modern (optical) diagnostics are often too complicated for manual inspection, manipulation and display. A simpler but still accurate description of these results is needed to facilitate data understanding. The paper discusses preliminary investigations into methods for the decomposition of 2-D and 3-D fluid flow data bases into elementary structures for purposes of description, analysis and comparison. An approach which involves the development of scene-like representation of the flow topology is presented. Using features such as critical points and dividing streamlines as a basis, a representation of the global topology of the flow is generated. The topology is then represented by a graph with the various structures represented by the nodes and their relationships in the flow by the connecting lines of the graph. Once the flow field has been placed in this form, it can be studied and compared with other data sets using techniques of syntactic pattern recognition or displayed using 3-D graphics techniques.

Hesselink, Lambertus↗

The computation of non-equilibrium chemically-reacting flows

Numerical techniques for the analysis of flows with nonequilibrium thermodynamics and chemistry are developed and demonstrated. Particular attention is given to the formulation of the thermodynamic and chemical gas models, the derivation of the governing equations, and the solution procedures. Results for typical problems involving shock tubes and supersonic diffusers are presented in extensive graphs and briefly characterized.

Grossman, B.↗

Effects of Heat Flux, Oxygen Concentration and Glass Fiber Volume Fraction on Pyrolysate Mass Flux from Composite Solids

Experimental work on the effects of heat flux, oxygen concentration and glass fiber volume fraction on pyrolysate mass flux from samples of polypropylene/glass fiber composite (PP/G) is underway. The research is conducted as part of a larger project to develop a test methodology for flammability of materials, particularly composites, in the microgravity and variable oxygen concentration environment of spacecraft and space structures. Samples of PP/G sized at 30x30x10 mm are flush mounted in a flow tunnel, which provides a flow of oxidizer over the surface of the samples at a fixed value of 1 m/s and oxygen concentrations varying between 18 and 30%. Each sample is exposed to a constant external radiant heat flux at a given value, which varies between tests from 10 to 24 kW/m2. Continuous sample mass loss and surface temperature measurements are recorded for each test. Some tests are conducted with an igniter and some are not. In the former case, the research goal is to quantify the critical mass flux at ignition for the various environmental and material conditions described above. The later case generates a wider range of mass flux rates than those seen prior to ignition, providing an opportunity to examine the protective effects of blowing on oxidative pyrolysis and heating of the surface. Graphs of surface temperature and sample mass loss vs. time for samples of 30% PPG at oxygen concentrations of 18 and 21% are presented in the figures below. These figures give a clear indication of the lower pyrolysis rate and extended time to ignition that accompany a lower oxygen concentration. Analysis of the mass flux rate at the time of ignition gives good repeatability but requires further work to provide a clear indication of mass flux trends accompanying changes in environmental and material properties.

Rich, D. B.↗

Effects of Heat Flux, Oxygen Concentration and Glass Fiber Volume Fraction on Pyrolysate Mass Flux from Composite Solids

Experimental work on the effects of heat flux, oxygen concentration and glass fiber volume fraction on pyrolysate mass flux from samples of polypropylene/glass fiber composite (PP/G) is underway. The research is conducted as part of a larger project to develop a test methodology for flammability of materials, particularly composites, in the microgravity and variable oxygen concentration environment of spacecraft and space structures. Samples of PP/G sized at 30 x 30 x 10 mm are flush mounted in a flow tunnel, which provides a flow of oxidizer over the surface of the samples at a fixed value of 1 m/s and oxygen concentrations varying between 18 and 30%. Each sample is exposed to a constant external radiant heat flux at a given value, which varies between tests from 10 to 24 kW/sq m. Continuous sample mass loss and surface temperature measurements are recorded for each test. Some tests are conducted with an igniter and some are not. In the former case, the research goal is to quantify the critical mass flux at ignition for the various environmental and material conditions described above. The later case generates a wider range of mass flux rates than those seen prior to ignition, providing an opportunity to examine the protective effects of blowing on oxidative pyrolysis and heating of the surface. Graphs of surface temperature and sample mass loss vs. time for samples of 30% PPG at oxygen concentrations of 18 and 21% are presented in the figures below. These figures give a clear indication of the lower pyrolysis rate and extended time to ignition that accompany a lower oxygen concentration. Analysis of the mass flux rate at the time of ignition gives good repeatability but requires further work to provide a clear indication of mass flux trends accompanying changes in environmental and material properties.

Rich, D. B.↗

Automated Processing of ISIS Topside Ionograms into Electron Density Profiles

Modeling of the topside ionosphere has for the most part relied on just a few years of data from topside sounder satellites. The widely used Bent et al. (1972) model, for example, is based on only 50,000 Alouette 1 profiles. The International Reference Ionosphere (IRI) (Bilitza, 1990, 2001) uses an analytical description of the graphs and tables provided by Bent et al. (1972). The Alouette 1, 2 and ISIS 1, 2 topside sounder satellites of the sixties and seventies were ahead of their times in terms of the sheer volume of data obtained and in terms of the computer and software requirements for data analysis. As a result, only a small percentage of the collected topside ionograms was converted into electron density profiles. Recently, a NASA-funded data restoration project has undertaken and is continuing the process of digitizing the Alouette/ISIS ionograms from the analog 7-track tapes. Our project involves the automated processing of these digital ionograms into electron density profiles. The project accomplished a set of important goals that will have a major impact on understanding and modeling of the topside ionosphere: (1) The TOPside Ionogram Scaling and True height inversion (TOPIST) software was developed for the automated scaling and inversion of topside ionograms. (2) The TOPIST software was applied to the over 300,000 ISIS-2 topside ionograms that had been digitized in the fkamework of a separate AISRP project (PI: R.F. Benson). (3) The new TOPIST-produced database of global electron density profiles for the topside ionosphere were made publicly available through NASA s National Space Science Data Center (NSSDC) ftp archive at . (4) Earlier Alouette 1,2 and ISIS 1, 2 data sets of electron density profiles from manual scaling of selected sets of ionograms were converted fiom a highly-compressed binary format into a user-friendly ASCII format and made publicly available through nssdcftp.gsfc.nasa.gov. The new database for the topside ionosphere established as a result of this project, has stimulated a multitude of new studies directed towards a better description and prediction of the topside ionosphere. Marinov et al. (2004) developed a new model for the upper ion transition height (Oxygen to Hydrogen and Helium) and Bilitza (2004) deduced a correction term for the I N topside electron density model. Kutiev et al. (2005) used this data to develop a new model for the topside ionosphere scale height (TISH) as a function of month, local time, latitude, longitude and solar flux F10.7. Comparisons by Belehaki et al. (2005) show that TISH is in general agreement with scale heights deduced from ground ionosondes but the model predicts post-midnight and afternoon maxima whereas the ionosonde data show a noon maximum. Webb and Benson (2005) reported on their effort to deduce changes in the plasma temperature and ion composition from changes in the topside electron density profile as recorded by topside sounders. Limitations and possible improvements of the IRI topside model were discussed by Coisson et al. (2005) including also the possible use of the NeQuick model, Our project progressed in close collaboration and coordination with the GSFC team involved in the ISIS digitization effort. The digitization project was highly successful producing a large amount of digital topside ionograms. Several no-cost extensions of the TOPIST project were necessary to keep up with the pace and volume of the digitization effort.

Reinisch, bodo W.↗

Historical Pattern Analysis of Global Geothermal Power Capacity Development: Preprint

Between 1913 and 1958, Italy was the only country with an operational geothermal power plant until New Zealand installed its first plant in 1958. At present, 24 countries are involved in the geothermal power market, and they have a combined installed capacity of 16,127 GW. This study analyzes the historical patterns of geothermal power capacity in the world and in individual countries to investigate the ideal global geothermal development pattern by examining the annual cumulative capacity (ACC) and the annual capacity addition (ACA) graphs of the historical development of geothermal power capacity in 24 countries. First, the global patterns are analyzed using these graphs in five periods (1945-1957, 1958-1976, 1977-1991, 1992-2002, and 2003-2020) that are marked by a series of characteristics of ACA peaks separated by two major troughs. Then, five characteristic patterns are developed in five periods globally. These patterns correspond to the early-stage linear, the first acceleration, the first steady-state linear, the second acceleration, and the second steady-state linear developments. A positive relationship exists between global patterns and the 5-year shifted oil-price curve: two major factors influenced global development: 1) increasing oil prices and increasing awareness of global climate change, and 2) global development of geothermal power. Last, we investigate these patterns in each country. The top ten countries, which comprise 93.3% of the world's total installed capacity are separated into five groups based on the availability and characteristics of patterns globally developed in five periods. Group-1 (the United States) has an installed capacity of 3,794 MW, Group-2 (Mexico and Philippines) 963-1935 MW, Group-3 (New Zealand, Italy, Iceland, and Japan) 601-1,037 MW, and Group-4 (Indonesia, Kenya, and Turkiye) 944-2,356 MW. The remaining 14 countries (6.7%), which are called Group 5, are still in an immature stage and have installed capacities of 7-262 MW and are not involved in pattern analysis. Overall, geothermal power in the world is in its third stage of development, which had its peak development after 1977. A fourth development peak may be expected to occur after this through business-as-usual cases. The biggest barrier to the development of the global geothermal power market is the risk associated with exploration and drilling. If risk mitigation systems and funds are employed, the growth of geothermal power production projects could accelerate.

GEOTHERMAL ENERGY↗

Analysis of a Computational Framework for Bayesian Inverse Problems: Ensemble Kalman Updates and MAP Estimators under Mesh Refinement

This paper analyzes a popular computational framework to solve infinite-dimensional Bayesian inverse problems, discretizing the prior and the forward model in a finite-dimensional weighted inner product space. We demonstrate the benefit of working on a weighted space by establishing operator-norm bounds for finite element and graph-based discretizations of Matérn-type priors and deconvolution forward models. For linear-Gaussian inverse problems, we develop a general theory to characterize the error in the approximation to the posterior. We also embed the computational framework into ensemble Kalman methods and MAP estimators for nonlinear inverse problems. Furthermore, our operator-norm bounds for prior discretizations guarantee the scalability and accuracy of these algorithms under mesh refinement.

Bayesian inverse problem↗

Defect Diffusion Graph Neural Networks for Materials Discovery in High-Temperature Energy Applications

Here, the migration of crystallographic defects dictates material properties and performance for a plethora of technological applications. Density functional theory (DFT)-based nudged elastic band (NEB) calculations are a powerful computational technique for predicting defect migration activation energy barriers, yet they become prohibitively expensive for high-throughput screening of defect diffusivities. Without introducing hand-crafted (i.e., chemistry- or structure-specific) descriptors, we propose a generalized deep learning approach to train surrogate models for NEB energies of vacancy migration by hybridizing graph neural networks with transformer encoders and simply using pristine host structures as input. With sufficient training data, computationally efficient and simultaneous inference of vacancy defect thermodynamics and migration activation energies can be obtained to compute temperature-dependent vacancy diffusivities and to down-select candidates for more thorough DFT analysis or experiments. Thus, as we specifically demonstrate for potential water-splitting materials, candidates with desired defect thermodynamics, kinetics, and host stability properties can be more rapidly targeted from open-source databases of experimentally validated or hypothetical materials.

14 SOLAR ENERGY↗

BCSR on GPU: A Way Forward Extreme-scale Graph Processing on Accelerator-enabled Frontier Supercomputer

Handling large graphs in a distributed environment requires effective partitioning across processors and efficient management of local partitions. In 2D partitioning, local graphs often become too sparse, making memory-efficient data structures crucial. Using the Compressed Sparse Row (CSR) format wastes space, especially for > 83% of vertices with empty edges for the sparse graphs. This study explores bit-CSR (BCSR), a modified CSR representation, on GPUs to reduce memory usage in graph computations. We achieved 16.67% memory savings on a sparse rmat dataset with 268 million vertices and 357 million edges, without performance degradation, supported by both theoretical and experimental storage savings of 33%. However, we observed a 1.7× slowdown in degree lookup times due to bitwise operations on AMD CPUs. This analysis highlights the potential of BCSR on GPUs for improving Graph500 benchmark performance on GPU-accelerated systems, such as the Frontier supercomputer.

Sattar, Naw Safrin↗

Optimal payload ascent trajectories of winged vehicles

A trajectory analysis is presented which shows the advantage of lifting ascent and the penalties associated with increased wing mass due to increases in aerodynamic loading. The vehicle analyzed is made up of an orbiter, an expendable liquid fuel tank, and two solid rocket motors. Results of the tradeoff are illustrated in a graph of the total vehicle mass plotted against the trajectory constraint. Results indicate that even when structural mass penalties are considered lifting ascent provides a payload delivery.

Martin, J. A.↗