Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel data structures”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Towards Superior Software Portability with SHAD and HPX C++ Libraries

As hardware architectures and software stacks complexity grows, development productivity, performance and software portability, quickly evolve from desirable features to actual needs. SHAD, the Scalable High-performance Algorithms and Data-structures C++ library is designed to mitigate these issues: it provides general purpose building blocks as well as high-level custom utilities, and offers a shared-memory programming abstraction which facilitates the programming of complex systems, scaling up to High Performance Computing clusters. SHAD’s portability is achieved through an abstract runtime interface, which decouples the upper layers of the library and hides the low level details of the underlying architecture. This layer enables SHAD to interface with different runtime/threading systems, e.g. Intel TBB and Global Memory and Threading (GMT). However, current backends targeting distributed systems, rely on a centralized controller which may possibly limit scalability up to hundreds of nodes and creates a network hot spot due to all to one communication for synchronization, and possibly resulting in degraded performance at high process counts. In this research, we explore HPX, the C++ standard library for parallelism and concurrency, as an additional backend in support of the SHAD library, and present the methodologies in support of local and remote task executions in SHAD with respect to HPX. Finally, we evaluate the proposed system by comparing against existing backends of SHAD and analyzing their performance on C++ Standard Template Library algorithms.

Wu, Nanmiao↗
DOE OSTI2021

Understanding Lustre Internals. Second Edition

The Lustre file system has become a preferred storage resource for systems on the Top500 list, and it is often the file system of choice for small- to medium-sized HPC systems that require parallel shared access to data. Several resources exist to help users deploy and configure Lustre, but the same cannot be said for resources that explain the inner workings of the Lustre source code. A previous ORNL technical report entitled "Understanding Lustre Filesystem Internals" (ORNL/TM-2009/117) provided an excellent summary of Lustre subsystem operations. However, that report is over a decade old and is based on Lustre version 1.6. Since that report was published, Lustre has evolved significantly. Several subsystems underwent significant code changes and many new features have been added to the file system, bringing the current Lustre version up to 2.15.This report aims to document and explain the internal workings of the latest version of the Lustre file system. It will provide more complete and up-to-date information than the previous technical report and should serve as a foundational document for anyone interested in Lustre software development. Key data structures will be described along with the APIs used for interaction among the various Lustre subsystems. Although the Lustre software is constantly being developed, the details in this document should remain relevant for the forseeable future.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Nozomi Cis-Lunar Phase Orbit Determination

Japan's Institute of Space and Astronautical Science (ISAS) launched Nozomi, its first mission to the planet Mars using the newly developed M-V launch vehicle on July 3, 1998. Scientific objectives of the mission are to study the structure and dynamics of the Martian upper atmosphere and its interaction with the solar wind. Nozomi is a cooperative mission between ISAS and the National Aeronautics and Space Administration (NASA). The NASA contribution includes navigation and tracking services provided by the Jet Propulsion Laboratory (JPL). The spacecraft also serves as an engineering demonstration of basic technology for planetary exploration. One of the new technologies was a unique trajectory, developed by ISAS, which used solar gravitational perturbations at the weak stability boundary as an aid to achieve an Earth-Mars transfer orbit. This trajectory saves approximately 120 m/s of Delta V compared to direct hyperbolic insertion and is considered an enabling technology for the mission. Nozomi was the first spacecraft to employ this trajectory and provided on-orbit validation of the technique. The trajectory was achieved by initially placing the spacecraft in a highly elliptical cis-lunar phasing orbit. Six maneuvers were performed during this period to correct injection errors and target an outbound lunar swingby in September 1998. The gravity assist from the lunar swingby raised apogee to the vicinity of the weak stability boundary. After three more targeting maneuvers, Nozomi performed an inbound lunar swingby followed immediately by a powered Earth swingby in late December 1998. A 420 m/s Trans Mars Insertion (TMI) burn at the final Earth periapsis was intended to place the spacecraft on a heliocentric trajectory leading to Mars orbit insertion in October 1999. Orbit determination for Nozomi is performed in parallel by both ISAS and the Multi-Mission Navigation (MMNAV) group at JPL. This was an advantage for the mission because each group would generate solutions based on data collected from their respective tracking networks. Spacecraft events, such as sequence uplinks and maneuvers, were generally scheduled during passes at the Usuda tracking station in Japan. As a result, maneuver design and reconstruction was derived from MMNAV solutions based on JPL tracking data obtained immediately prior to or following maneuvers. Data was also exchanged between ISAS and MMNAV so orbit determination could be performed on joint data sets in support of critical targeting late in the cis-lunar phase. In this paper, information regarding the MMNAV orbit determination effort for the first six months of the mission is presented. The spacecraft trajectory is characterized first, followed by a discussion of the orbit determination estimation procedure and models. Results from selected orbit solutions are presented and compared against reconstructed trajectories. One area of emphasis in this paper is orbit determination in the vicinity of the weak stability boundary. Precise navigation was necessary to target the second lunar swingby and the powered Earth swingby. Delivery accuracy of 150 m was required for these critical encounters, but a number of factors contributed to the general degradation of orbit determination accuracy. This included the fact that the spacecraft was at apogee, at a range of 1.7 million km and moving at less than I km/sec perpendicular to the line of sight. Nozomi was also close to zero degrees declination where there are known limitations on orbit determination performance. Finally, S-band tracking data was acquired through the Nozomi backup low gain antenna. This antenna is offset from the axis of this spin stabilized spacecraft and superimposed large signatures in the Doppler and range data. These difficulties were overcome by combining long data arcs, spanning several maneuvers, with a high fidelity solar pressure model. The model included a physically accurate representation of the spacecraft structure and a high time resolution orientation model. Observation modeling included the removal of the spin induced Doppler bias, spin signature and per pass correction of range calibration errors applied for data leading up to critical events. As a result, all orbit determination goals were met. A second area of emphasis in this paper is the JPL tracking and orbit determination effort in support of the TMI maneuver. TMI occurred out of contact with ground stations and the JPL Goldstone tracking complex had the first pass following the bum. As a result, MMNAV had the responsibility to make a rapid assessment of the maneuver performance. MMNAV made the determination that a 100 m/s under bum had occurred and promptly informed ISAS via voice lines. ISAS immediately began preparations for a correction maneuver (TMIc), which had to be performed during the next Usuda pass. The near real time assessment by MMNAV provided accurate antenna frequency and pointing updates for the spacecraft acquisition at Usuda and the close coordination between the two agencies enabled the design and successful execution of the TMc maneuver. Propellant consumption during the correction burn dictated that the mission be redesigned. ISAS developed a new plan which adds 3 full solar orbits, two Earth swingbys and one lunar swingby with arrival at Mars in January 2004. The final Mars orbit will still enable the mission to achieve all of its science objectives.

Ryne, Mark↗
NASA NTRS2004

A Massively Parallel Hybrid Dusty-Gasdynamics and Kinetic Direct Simulation Monte Carlo Model for Planetary Applications

In order to understand the global structure, dynamics, and physical and chemical processes occurring in the upper atmospheres, exospheres, and ionospheres of the Earth, the other planets, comets and planetary satellites and their interactions with their outer particles and fields environs, it is often necessary to address the fundamentally non-equilibrium aspects of the physical environment. These are regions where complex chemistry, energetics, and electromagnetic field influences are important. Traditional approaches are based largely on hydrodynamic or magnetohydrodynamic (MHD) formulations and are very important and highly useful. However, these methods often have limitations in rarefied physical regimes where the molecular collision rates and ion gyrofrequencies are small and where interactions with ionospheres and upper neutral atmospheres are important. At the University of Michigan we have an established base of experience and expertise in numerical simulations based on particle codes which address these physical regimes. The Principal Investigator, Dr. Michael Combi, has over 20 years of experience in the development of particle-kinetic and hybrid kinetichydrodynamics models and their direct use in data analysis. He has also worked in ground-based and space-based remote observational work and on spacecraft instrument teams. His research has involved studies of cometary atmospheres and ionospheres and their interaction with the solar wind, the neutral gas clouds escaping from Jupiter s moon Io, the interaction of the atmospheres/ionospheres of Io and Europa with Jupiter s corotating magnetosphere, as well as Earth s ionosphere. This report describes our progress during the year. The contained in section 2 of this report will serve as the basis of a paper describing the method and its application to the cometary coma that will be continued under a research and analysis grant that supports various applications of theoretical comet models to understanding the inner comae of comets (grant NAGS- 13239 from the Planetary Atmospheres program).

Combi, Michael R.↗
NASA NTRS2015

Quantification of Plume-Soil Interaction and Excavation Due to the Sky Crane Descent Stage

The quantification of the particulate erosion that occurs as a result of a rocket exhaust plume impinging on soil during extraterrestrial landings is critical for future robotic and human lander mission design. The aerodynamic environment that results from the reflected plumes results in dust lifting, site alteration and saltation, all of which create a potentially erosive and contaminant heavy environment for the lander vehicle and any surrounding structures. The Mars Science Lab (MSL), weighing nearly one metric ton, required higher levels of thrust from its retro propulsive systems and an entirely new descent system to minimize these effects. In this work we seek to quantify plume soil interaction and its resultant soil erosion caused by the MSL's Sky Crane descent stage engines by performing three dimensional digital terrain and elevation mapping of the Curiosity rover's landing site. Analysis of plume soil interaction altitude and time was performed by detailed examination of the Mars Descent Imager (MARDI) still frames and reconstructed inertial measurement unit (IMU) sensor data. Results show initial plume soil interaction from the Sky Crane's eight engines began at ground elevations greater than 60 meters and more than 25 seconds before the rovers' touchdown event. During this time, viscous shear erosion (VSE) was dominant typically resulting in dusting of the surface with flow propagating nearly parallel to the surface. As the vehicle descended and decreased to four powered engines plume-plume and plume soil interaction increased the overall erosion rate at the surface. Visibility was greatly reduced at a height of roughly 20 meters above the surface and fell to zero ground visibility shortly after. The deployment phase of the Sky Crane descent stage hovering at nearly six meters above the surface showed the greatest amount of erosion with several large particles of soil being kicked up, recirculated, and impacting the bottom of the rover chassis. Image data obtained from MSL's navigation camera (NAVCAM) pairs on Sols 002, 003, and 016 were used to virtually recreate local surface topography and features around the rover by means of stereoscopic depth mapping. Images taken simultaneously by the left and right navigation cameras located on the rover's mast assembly spaced 42 centimeters were used to generate a three dimensional depth map from flat, two dimensional images of the same feature at slightly different angles. Image calibration with physical hardware on the rover and known terrain features were used to provide scaling information that accurately sizes features and regions of interest within the images. Digital terrain mapping analysis performed in this work describe the crater geometry (shape, radius, and depth), eroded volume, volumetric erosion rate, and estimated mass erosion rate of the Hepburn, Sleepy Dragon, Burnside, and Goulburn craters. Crater depths ranged from five to ten centimeters deep influencing an area as wide as two meters in some cases. The craters formed were highly asymmetrical and generally oblong primarily due to the underlying bedrock formations underneath the surface. Comparison with ground tests performed at the NASA AMES Planetary Aeolian Laboratory (PAL) by Mehta showed good agreement with volumetric erosion rates and crater sizes of large particle soil simulants, providing validation to Earth based ground tests of Martian regolith.

Vizcaino, Jeffrey↗
DOE OSTI2018

Characterization of the ePix100a and the FastCCD Semiconductor Detectors for the European XFEL

Here, the European X-ray Free Electron Laser (EuXFEL) is a research facility providing spatially coherent X-ray flashes in the energy range from 0 . 25 keV to 25 keV of unprecedented brilliance and with unique time structure: X-ray pulses with a 4 . 5 MHz repetition rate arranged in trains with 2700 pulses every 100 ms. The facility operates three photon beamlines called SASE 1, SASE 2 and SASE 3. Each of the beamlines is hosting two scientific experiments. The SASE 1 beamline started its user operation in September 2017, followed by successful first lasing at the SASE 2 beamline in May 2018. Early user experiments are planned to start in 2019 at this beamline, while early user experiments for the SASE 3 beamline are scheduled for the end of 2018. The quality of the experimental data will gain substantial benefits from an accurate characterization and calibration of the X-ray detectors. Supplementing high repetition rate detectors at MHz speeds, slower detectors such as the ePix100a and the FastCCD will be operated at the train repetition rate of 10 Hz. These 2D silicon pixelized detectors use fast parallel column-wise readout implemented as a CCD or as a hybrid pixel detector. In the following, characterization and analysis approaches for the FastCCD and the ePix100a detectors are discussed and the performance of the detectors is evaluated using appropriate state-of-the-art analysis techniques.

47 OTHER INSTRUMENTATION↗
NASA NTRS2022

X-ray Polarization Detection of Cassiopeia A with IXPE

We report on a ∼5σ detection of polarized 3–6 keV X-ray emission from the supernova remnant Cassiopeia A (Cas A) with the Imaging X-ray Polarimetry Explorer (IXPE). The overall polarization degree of 1.8% ± 0.3% is detected by summing over a large region, assuming circular symmetry for the polarization vectors. The measurements imply an average polarization degree for the synchrotron component of ∼2.5%, and close to 5% for the X-ray synchrotron-dominated forward shock region. These numbers are based on an assessment of the thermal and nonthermal radiation contributions, for which we used a detailed spatial-spectral model based on Chandra X-ray data. A pixel-by-pixel search for polarization provides a few tentative detections from discrete regions at the ∼ 3σ confidence level. Given the number of pixels, the significance is insufficient to claim a detection for individual pixels, but implies considerable turbulence on scales smaller than the angular resolution. Cas A's X-ray continuum emission is dominated by synchrotron radiation from regions within ≲1017 cm of the forward and reverse shocks. We find that (i) the measured polarization angle corresponds to a radially oriented magnetic field, similar to what has been inferred from radio observations; (ii) the X-ray polarization degree is lower than in the radio band (∼5%). Since shock compression should impose a tangential magnetic-field structure, the IXPE results imply that magnetic fields are reoriented within ∼1017 cm of the shock. If the magnetic-field alignment is due to locally enhanced acceleration near quasi-parallel shocks, the preferred X-ray polarization angle suggests a size of 3 × 1016 cm for cells with radial magnetic fields.

Polarimetry↗
DOE OSTI2023

Massively Parallel Bayesian Model Calibration and Uncertainty Quantification with Applications to Nuclear Fuels and Materials

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. Fuel behavior in a nuclear reactor is governed by a complex network of mechanisms interacting with various other physics aspects in the reactor system. Any model developed to represent the fuel behavior will likely be idealized resulting in uncertainties in their predictions compared to the observed data. As such, this report was motivated by the need to identify the sources of uncertainties and quantify and propagate them through the fuel model outputs. Such quantification of uncertainties will establish a level of model trustworthiness, identify approaches to improve the model trustworthiness, and even guide optimal experiment design for maximal information gain. To accomplish the uncertainty quantification for computational models, this report has relied on the Bayesian framework which provides probabilistic treatment of models their inputs and outputs. The current state-of-the-art on performing Bayesian Uncertainty Quantification (UQ) for nuclear engineering models using High Performance Computing (HPC) resources have been reviewed. Implementation of capabilities for massively parallel Bayesian UQ in Multiphysics Object-Oriented Simulation Environment (MOOSE) is discussed. Several verification cases are discussed to verify the accuracy of the quantified uncertainties using the developed computational capabilities in MOOSE. Then, the problem of quantifying the uncertainties in TRI-Structural isOtropic (TRISO) fuel silver release is addressed. For the first time, the uncertainties arising from the TRISO Fission Gas Release (FGR) model due to model inadequacy and experimental noise are quantified. Also, the Bayesian capabilities are applied to the calibration of the MATPRO creep model, a widely used model in several fuel assessment cases. The impact of the prediction uncertainties in the MATPRO model on the fuel cladding behavior as part of the TRIBULATION assessment case (which is an integral effects case) is investigated. This report concludes with a discussion on the future work for the UQ for computational models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2010

The OMPS Limb Profiler Instrument: Two-Dimensional Retrieval Algorithm

The upcoming Ozone Mapper and Profiler Suite (OMPS), which will be launched on the NPOESS Preparatory Project (NPP) platform in early 2011, will continue monitoring the global distribution of the Earth's middle atmosphere ozone and aerosol. OMPS is composed of three instruments, namely the Total Column Mapper (heritage: TOMS, OMI), the Nadir Profiler (heritage: SBUV) and the Limb Profiler (heritage: SOLSE/LORE, OSIRIS, SCIAMACHY, SAGE III). The ultimate goal of the mission is to better understand and quantify the rate of stratospheric ozone recovery. The focus of the paper will be on the Limb Profiler (LP) instrument. The LP instrument will measure the Earth's limb radiance (which is due to the scattering of solar photons by air molecules, aerosol and Earth surface) in the ultra-violet (UV), visible and near infrared, from 285 to 1000 nm. The LP simultaneously images the whole vertical extent of the Earth's limb through three vertical slits, each covering a vertical tangent height range of 100 km and each horizontally spaced by 250 km in the cross-track direction. Measurements are made every 19 seconds along the orbit track, which corresponds to a distance of about 150km. Several data analysis tools are presently being constructed and tested to retrieve ozone and aerosol vertical distribution from limb radiance measurements. The primary NASA algorithm is based on earlier algorithms developed for the SOLSE/LORE and SAGE III limb scatter missions. All the existing retrieval algorithms rely on a spherical symmetry assumption for the atmosphere structure. While this assumption is reasonable in most of the stratosphere, it is no longer valid in regions of prime scientific interest, such as polar vortex and UTLS regions. The paper will describe a two-dimensional retrieval algorithm whereby the ozone distribution is simultaneously retrieved vertically and horizontally for a whole orbit. The retrieval code relies on (1) a forward 2D Radiative Transfer code (to model limb radiances within a non-uniform atmosphere and evaluate 2D analytical partial derivatives) and (2) an optimal estimator inversion routine. The algorithm uses the typically sparse nature of the kernel matrices as well as fast matrix inversion techniques to allow for fast inversion of limb data with efficient memory management (as was done for MIPAS data processing). While the method has so far only been developed in the context of Single Scatter, the paper will show how the CPU intensive Multiple Scatter modeling can be implemented using parallel CPU processing. Initial results will be presented in terms of retrieved ozone profiles and code performance.

Rault, Didier F.↗
NASA NTRS1994

Observation of anomalous slow-mode shock and reconnection layer in the dayside magnetospause

Plasma and magnetic field data from the International Sun-Earth Explorer ISEE 2 spacecraft recorded during an outbound crossing of the dayside, northern hemisphere magnetopause in October 29, 1979, provide evidence for a slow shock (SS) in the observed reconnection layer. This layer is found to be bounded on the magnetosheath side by the SS; near the magnetospheric side of the layer, a second current sheet is found that may have been rotational discontinuity (RD). The direction of the accelerated plasma flow, the earthward sense of the normal magnetic filed across the SS and RD, and the relative orientation of the SS and the RD all indicate that the reconection site was located south of the spacecraft. Quantitative tests show that, allowing for experiemental uncertanties, data taken upstream and downstream of the SS are consistent with coplanarity and other Rankine-Hugoniot (RH) conditions. Examiniations of the flow parameters indicates two anomalous properties of the SS: the upstream flow, viewed in the deHoffmann-Teller frame, is superalfvenic and the downstream plasma is firehose unstable. In comparison to the long-wavelength slow-mode phase speed, however, the flow in the upstream region is super slow, while in the downstream region it is subslow, as requeired for a slow-mode shock. Further properties of the shock include a large decrease in total enthalpy across it, indicating the escape of a sizable heat flux from the shock structure, and the occurrence of a polarization reversal of the tangential magnetic field within the shock layer, a feature that is predicted by linear double-polytropic Hall-MHD and results from a large increase in pressure anistropy from the upstream region, where p(sub parallel) approximately equal to p(sub perp), to the downstream region, where p(sub parrallel) greater than p(sub perp). Quantitative tests of the RD-like discontinuity show that it satisfies the necessary RH conditions within experimental uncertainties and indicate that the flow across it is earthward, that is, it is the same as for the SS. Thus the RD is propagating behind the SS in the direction away from the Earth. The reversal in the order of the discontinuities is attributed to the reversal of the slow and intermediate phase speeds produced by the pressure anistropy in the region between the RD and the SS. The linear dispersion relation and polarization properties of dispersive MHD waves in an anisotropic plasma are examined in the appendix.

Walthour, D. W↗
NASA NTRS1979

Estimation of the ocean geoid near the Blake Escarpment using GEOS-3 satellite altimetry

The accuracy with which the local ocean geoid structure could be determined using satellite altimetry data was investigated. The undulation and along-track component of the vertical deflection for selected passes of GEOS-3 near the Blake Escarpment were estimated and compared with independent analogous estimates based on U. S. Navy surface gravimetric survey data. The results of these comparisons show agreement in the geoid undulation values generally to within one or two meters. The nature of the discrepancy in the undulation values was primarily that of a bias error believed to be due essentially to radial orbit uncertainties. The agreement between the vertical deflection estimates was not significantly affected by orbit uncertainties over the track lengths considered in this study (100 - 1500 km), and the comparisons show typical rms differences of between one and two arc secs. In addition, the capability of the altimeter to resolve short wavelength features of the geoid was determined. This analysis involved spectrum and cross spectrum analysis of sets of closely spaced parallel subtracks to determine statistically significant short wavelength geoid resolution capability. The results of this analysis show that resolution can be achieved down to wavelengths as short as 30 km - 80 km depending on regional geoid variations.

Brammer, R. F.↗
NASA NTRS1994

Analytic model of aurorally coupled magnetospheric and ionospheric electrostatic potentials

This paper describes modest but significant improvements on earlier studies of electrostatic potential structure in the auroral region using the adiabatic auroral arc model. This model has crucial nonlinearities (connected, for example. with aurorally produced ionization) which have hampered analysis; earlier work has either been linear, which I will show is a poor approximation or, if nonlinear, either numerical or too specialized to study parametric dependencies. With certain simplifying assumptions I find new analytic nonlinear solutions fully exhibiting the parametric dependence of potentials on magnetospheric (e.g.. cross-tail potential) and ionospheric (e.g., recombination rate) parameters. No purely phenomenological parameters are introduced. The results are in reasonable agreement with observed average auroral potential drops, inverted-V scale sizes, and dissipation rates. The dissipation rate is quite comparable to tail energization and transport rates and should have a major effect on tail and magnetospheric dynamics. This paper gives various relations between the cross-tail potential and auroral parameters (e.g., total parallel currents and potential drops) which can be studied with existing data sets.

Cornwall, J. M.↗
NASA NTRS1999

Development and Validation of a Fast, Accurate and Cost-Effective Aeroservoelastic Method on Advanced Parallel Computing Systems

Progress to date towards the development and validation of a fast, accurate and cost-effective aeroelastic method for advanced parallel computing platforms such as the IBM SP2 and the SGI Origin 2000 is presented in this paper. The ENSAERO code, developed at the NASA-Ames Research Center has been selected for this effort. The code allows for the computation of aeroelastic responses by simultaneously integrating the Euler or Navier-Stokes equations and the modal structural equations of motion. To assess the computational performance and accuracy of the ENSAERO code, this paper reports the results of the Navier-Stokes simulations of the transonic flow over a flexible aeroelastic wing body configuration. In addition, a forced harmonic oscillation analysis in the frequency domain and an analysis in the time domain are done on a wing undergoing a rigid pitch and plunge motion. Finally, to demonstrate the ENSAERO flutter-analysis capability, aeroelastic Euler and Navier-Stokes computations on an L-1011 wind tunnel model including pylon, nacelle and empennage are underway. All computational solutions are compared with experimental data to assess the level of accuracy of ENSAERO. As the computations described above are performed, a meticulous log of computational performance in terms of wall clock time, execution speed, memory and disk storage is kept. Code scalability is also demonstrated by studying the impact of varying the number of processors on computational performance on the IBM SP2 and the Origin 2000 systems.

Goodwin, Sabine A.↗
NASA NTRS2016

DeepSAT: A Deep Learning Approach to Tree-Cover Delineation in 1-m NAIP Imagery for the Continental United States

High resolution tree cover classification maps are needed to increase the accuracy of current land ecosystem and climate model outputs. Limited studies are in place that demonstrates the state-of-the-art in deriving very high resolution (VHR) tree cover products. In addition, most methods heavily rely on commercial softwares that are difficult to scale given the region of study (e.g. continents to globe). Complexities in present approaches relate to (a) scalability of the algorithm, (b) large image data processing (compute and memory intensive), (c) computational cost, (d) massively parallel architecture, and (e) machine learning automation. In addition, VHR satellite datasets are of the order of terabytes and features extracted from these datasets are of the order of petabytes. In our present study, we have acquired the National Agriculture Imagery Program (NAIP) dataset for the Continental United States at a spatial resolution of 1-m. This data comes as image tiles (a total of quarter million image scenes with ~60 million pixels) and has a total size of ~65 terabytes for a single acquisition. Features extracted from the entire dataset would amount to ~8-10 petabytes. In our proposed approach, we have implemented a novel semi-automated machine learning algorithm rooted on the principles of "deep learning" to delineate the percentage of tree cover. Using the NASA Earth Exchange (NEX) initiative, we have developed an end-to-end architecture by integrating a segmentation module based on Statistical Region Merging, a classification algorithm using Deep Belief Network and a structured prediction algorithm using Conditional Random Fields to integrate the results from the segmentation and classification modules to create per-pixel class labels. The training process is scaled up using the power of GPUs and the prediction is scaled to quarter million NAIP tiles spanning the whole of Continental United States using the NEX HPC supercomputing cluster. An initial pilot over the state of California spanning a total of 11,095 NAIP tiles covering a total geographical area of 163,696 sq. miles has produced true positive rates of around 88 percent for fragmented forests and 74 percent for urban tree cover areas, with false positive rates lower than 2 percent for both landscapes.

Imagery↗
NASA NTRS1989

Detection of edges using local geometry

Researchers described a new representation, the local geometry, for early visual processing which is motivated by results from biological vision. This representation is richer than is often used in image processing. It extracts more of the local structure available at each pixel in the image by using receptive fields that can be continuously rotated and that go to third order spatial variation. Early visual processing algorithms such as edge detectors and ridge detectors can be written in terms of various local geometries and are computationally tractable. For example, Canny's edge detector has been implemented in terms of a local geometry of order two, and a ridge detector in terms of a local geometry of order three. The edge detector in local geometry was applied to synthetic and real images and it was shown using simple interpolation schemes that sufficient information is available to locate edges with sub-pixel accuracy (to a resolution increase of at least a factor of five). This is reasonable even for noisy images because the local geometry fits a smooth surface - the Taylor series - to the discrete image data. Only local processing was used in the implementation so it can readily be implemented on parallel mesh machines such as the MPP. Researchers expect that other early visual algorithms, such as region growing, inflection point detection, and segmentation can also be implemented in terms of the local geometry and will provide sufficiently rich and robust representations for subsequent visual processing.

Gualtieri, J. A.↗
NASA NTRS1989

Incorporation of a progressive failure analysis method in the CSM testbed software system

Analysis of the postbuckling behavior of composite shell structures pose many difficult and challenging problems in the field of structural mechanics. Current analysis methods perform well for most cases in predicting the postbuckling response of undamaged components. To predict component behavior accurately at higher load levels, the analysis must include the effects of local material failures. The CSM testbed software system is a highly modular structural analysis system currently under development at Langley Research Center. One of the primary goals of the CSM testbed is to provide a software environment for the development of advanced structural analysis methods and modern numerical methods which will exploit advanced computer architecture such as parallel-vector processors. Development of a progressive failure analysis method consists of the design and implementation of a processor which will perform the ply-level progressive failure analysis and the development of a geometrically nonlinear analysis procedure which incorporates the progressive failure processor. Regarding the development of the progressive failure processor, two components are required: failure criteria and a degradation model. For the initial implementation, the failure criteria of Hashin will be used. For a matrix failure which typically indicates the development of transverse matrix cracks, the ply properties will be degraded. Work to date includes the design of the progressive failure analysis processor and initial plans for the controlling geometrically nonlinear analysis procedure. The implementation of the progressive failure analysis has begun. Access to the model database and the Hashin failure criteria are completed. Work is in progress on the input/output operations for the processor related data and the finite element model updating procedures. In total the progressive failure processor is approximately one-third complete.

Arenburg, Robert T.↗
DOE OSTI2024

Control And Optimization Modular Modeling Application For Nuclear Deployment

The purpose of the COMMAND code is to provide a flexible, scalable tool for use in developing, integrating, and testing the technologies necessary for achieving autonomous operations of advanced nuclear reactors. The code enables users to efficiently implement custom simulations and experiments by combining key methods from different software modules. These modules are focused on: modeling and simulation tools, such as nuclear simulation tools used for high-fidelity modeling (e.g., Reactor Excursion and Leak Analysis Program [RELAP5-3D] and Monte Carlo N-Particle [MCNP]); machine learning and optimization tools (e.g., anomaly detection and data-driven modeling techniques); advanced control in its digital, high-performance, and supervisory control forms (e.g., proportional integral derivative (PID) control and model predictive control (MPC); and integration with hardware through industrial communication protocols. To ensure flexibility and scalability, COMMAND was designed to be both modular—the software “pieces” all inherit from generic building blocks and can be combined and connected to create complicated simulations—and high performing—designed for parallel processing, enabling simulations and experiments to take advantage of multi-core computers, servers, and nodes. The code is written in the Python programming language due to the language's popularity, active community, and open-source and cross-platform nature. Maintaining consistency with other simulation tools used within the nuclear energy community, users implement simulations and experiments through text input files, which define components, parameters, connections, etc., through lines of text. Given that COMMAND is written in Python, these input files are native Python scripts, and so use the standard Python structure and formatting. This also enables users to take advantage of Python's extensive package library to develop custom capabilities for their specific use cases.

Faber, Jacob [Idaho National Laboratory (INL), Ida↗
DOE OSTI2021

An Overview of Non-LWR Vessel Cooling Systems for Passive Decay Heat Removal (Technical Letter Final Report)

The following report provides a technical review of various reactor vessel cooling system (VCS) concepts under consideration for decay heat removal (DHR) in non-light water advanced reactor designs. This review focuses on ex-vessel designs, including the Reactor Cavity Cooling System (RCCS), the Reactor Vessel Auxiliary Cooling System (RVACS), and hybrid iterations, using both air and water cooling to achieve their heat removal function. Based on a literature review of publicly available sources published between 1979 and 2021, a technical summary is presented detailing existing and planned VCS design options, their applicability to specific reactor type, and review of authored research and development studies. Following an assessment of the availability of data and modeling tools, an evaluation was performed analyzing their likely performance during normal, degraded, and accident conditions, including reliability, stability, and longevity. With renewed consideration of air-based DHR systems by some US vendors, there is a need to fully understand the complexities inherent to natural circulation systems, and more importantly, how they may influence safety-related heat removal functions and system performance. For DHR systems that rely on air-based mode of cooling, this entails quantification of the impacts of low flow conditions during start-up, relative elevations of inlet and outlet ducts for below-grade installations, effects of the use of multiple parallel chimneys for redundancy, and impact of weather conditions at the plant site. For systems relying on water-based mode cooling, reliance on limited coolant inventory in an event of pipe break or extended duration accident scenarios, sensitivities related to stability of boiling flow and heat transfer conditions, and quantification of full-scale structural vibrations on balance of plant structures are among the complicating factors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗