Search NASA⌕ Search

SEARCH · Search NASA

Results for “Implementation difficulty”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Satellite Mapping of PM2.5 Episodes in the Wintertime San Joaquin Valley: A "Static" Model Using Column Water Vapor

The use of satellite Aerosol Optical Thickness (AOT) from imaging spectrometers has been successful in quantifying and mapping high PM2.5 (particulate matter mass <2.5µm diameter) episodes for pollution abatement and health studies. However, some regions have high PM2.5 but poor estimation success. The challenges in using Aerosol Optical Thickness (AOT) from imaging spectrometers to characterize PM2.5 worldwide was especially evident in the wintertime San Joaquin Valley (SJV). The SJV's attendant difficulties of high-albedo surfaces and very shallow, variable vertical mixing also occur in other significantly polluted regions around the world. We report on more accurate PM2.5 maps for the whole-winter period in the SJV, Nov 14, 2012?Dec 11, 2013. Intensive measurements by including NASA aircraft were made for several weeks in that winter, the DISCOVER-AQ California mission.We found success with a relatively simple method based on calibration and checking with surface monitors and a characterization of vertical mixing, and incorporating specific understandings of the region's climatology. We estimate PM2.5 to within ~7µgm?3 RMSE and with R values of ~0.9, based on remotely sensed MAIAC (Multi-Angle Implementation of Atmospheric Correction) observations, and that certain further work will improve that accuracy. Mapping is at 1km resolution. This allows a time sequence of mapped aerosols at 1km for cloud-free days. We describe our technique as a "static estimation". Estimation procedures like this one, not dependent on well-mapped source strengths or on transport error, should help full source-driven simulations by deconstructing processes. They also provide a rapid method to create a long-term climatology.Essential features of the technique are (a) daily calibration of the AOT to PM2.5 using available surface monitors, and (b) characterization of mixed-layer dilution using column water vapor (CWV, otherwise "precipitable water"). We noted that on multi-day timescales both water vapor and particles share near-surface sources and both fall to very low values with altitude; indeed, both are largely removed by precipitation. The existence of layers of H2O or aerosol not within the mixed layer adds complexity, but mixed-effects statistical regression captures essential proportionality of PM2.5 and the ratio variable (AOT/CWV). Accuracy is much higher than previous statistical models, and can be extended to the whole Aqua-satellite data record. The maps and time-series we show suggest a repeated pattern for large valleys like the SJV ? progressive stabilization of the mixing height after frontal passages: PM2.5 is somewhat more determined by day-by-day changes in mixing than it is by the progressive accumulation of pollutants (revealed as increasing AOT).

Chatfield, Robert B.↗
DOE OSTI2018

Algorithm for Establishing Safeguards Termination Concentration Limits

This paper presents an algorithm for estimating the maximum concentration limits for terminating safeguards on nuclear materials. These limits are important because such material is no longer subject to accounting, reporting or inspections after the International Atomic Energy Agency (IAEA) has approved the termination of safeguards on nuclear material declared by a State to be waste.Maximum concentration limits for different waste forms generated within nuclear fuel cycle facilities must be sufficiently conservative such that termination of safeguards does not become a weak link in the safeguards system. The concentration limits must also remain practical in order to ensure IAEA resources are used effectively and efficiently to implement safeguards in the State. These termination limits should be technically-based while remaining objective and reasonable to implement by the State. In support of the original IAEA guidance on concentration limits for termination of safeguards prepared in the early 1990s, the U.S. developed an algorithm to estimate concentrations that would make recovery of nuclear material from waste on which safeguards had been terminated at least as unattractive as undeclared production from uranium ore or diversion of similar nuclear material. In 2016, the IAEA sought technical advice from a meeting of experts from selected member States to support updating its internal guidance on termination of safeguards. Their recommendations included extending the guidance to consider additional waste forms. The experts in the 2016 meeting recommended that the IAEA limits should more clearly reflect the technical difficulty and level of effort required to recover nuclear material from the various waste forms. This paper identifies and assesses the most likely techniques a State might use to recover nuclear material from the waste forms identified by the experts and determines how the safeguards termination algorithm could be revised to more accurately reflect the difficulty of recovering one significant quantity of nuclear material from these waste forms.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Qubit Lattice Algorithms based on the Schrodinger-Dirac representation of Maxwell Equations and their Extensions

It is well known that Maxwell equations can be expressed in a unitary Schrodinger-Dirac representation for homogeneous media. However, difficulties arise when considering inhomoge- neous media. A Dyson map points to a unitary field qubit basis, but the standard qubit lattice algorithm of interleaved unitary collision-stream operators must be augmented by some sparse non-unitary potential operators that recover the derivatives on the refractive indices. The effect of the steepness of these derivatives on two dimensional scattering is examined with simulations showing quite complex wavefronts emitted due to transmissions/reflections within the dielectric objects. Maxwell equations are extended to handle dissipation using Kraus operators. Then, our theoretical algorithms are extended to these open quantum systems. A quantum circuit diagram is presented as well as estimates on the required number of quantum gates for implementation on a quantum computer.

Vahala, George↗
DOE OSTI2023

Qubit Lattice Algorithms Based on the Schrödinger-Dirac Representation of Maxwell Equations and Their Extensions

It is well known that Maxwell equations can be expressed in a unitary Schrodinger-Dirac representation for homogeneous media. However, difficulties arise when considering inhomogeneous media. A Dyson map points to a unitary field qubit basis, but the standard qubit lattice algorithm of interleaved unitary collision-stream operators must be augmented by some sparse non-unitary potential operators that recover the derivatives on the refractive indices. Here, the effect of the steepness of these derivatives on two-dimensional scattering is examined with simulations showing quite complex wavefronts emitted due to transmissions/reflections within the dielectric objects. Maxwell equations are extended to handle dissipation using Kraus operators. Then, our theoretical algorithms are extended to these open quantum systems. A quantum circuit diagram is presented as well as estimates on the required number of quantum gates for implementation on a quantum computer.

2D electromagnetic scattering↗
NASA NTRS2002

Alternative Metrics for Evaluating the Resilence of Advanced Life Support Systems

Ensuring the safety of the crew is a key performance requirement of a life support system. However, a number of conceptual and practical difficulties arise when devising metrics to concretely measure the ability of a life support system to maintain critical functions in the presence of anticipated and unanticipated faults. Resilience is a dynamic property of a life support system that depends on the complex interactions between faults, controls and system hardware. We review some of the approaches to understanding the robustness or resilience of complex systems being developed in diverse fields such as ecology, software engineering and cell biology and discuss their applicability to regenerative life support systems. We also consider how approaches to measuring resilience vary depending on system design choices such as the definition and choice of the nominal operating regime. Finally, we explore data collection and implementation issues such as the key differences between the instantaneous or conditional and average or overall measures of resilience. Extensive simulation of a hybrid computational model of a water revitalization subsystem (WRS) with probabilistic, component-level faults provides data about off-nominal behavior of the system. The data are used to consider alternative measures of resilience as predictors of the system's ability to recover from component-level faults.

Bell, Ann Maria↗
NASA NTRS2004

Enhanced Software for Scheduling Space-Shuttle Processing

The Ground Processing Scheduling System (GPSS) computer program is used to develop streamlined schedules for the inspection, repair, and refurbishment of space shuttles at Kennedy Space Center. A scheduling computer program is needed because space-shuttle processing is complex and it is frequently necessary to modify schedules to accommodate unanticipated events, unavailability of specialized personnel, unexpected delays, and the need to repair newly discovered defects. GPSS implements constraint-based scheduling algorithms and provides an interactive scheduling software environment. In response to inputs, GPSS can respond with schedules that are optimized in the sense that they contain minimal violations of constraints while supporting the most effective and efficient utilization of space-shuttle ground processing resources. The present version of GPSS is a product of re-engineering of a prototype version. While the prototype version proved to be valuable and versatile as a scheduling software tool during the first five years, it was characterized by design and algorithmic deficiencies that affected schedule revisions, query capability, task movement, report capability, and overall interface complexity. In addition, the lack of documentation gave rise to difficulties in maintenance and limited both enhanceability and portability. The goal of the GPSS re-engineering project was to upgrade the prototype into a flexible system that supports multiple- flow, multiple-site scheduling and that retains the strengths of the prototype while incorporating improvements in maintainability, enhanceability, and portability.

Barretta, Joseph A.↗
NASA NTRS1999

An Adaptive Kalman Filter Using a Simple Residual Tuning Method

One difficulty in using Kalman filters in real world situations is the selection of the correct process noise, measurement noise, and initial state estimate and covariance. These parameters are commonly referred to as tuning parameters. Multiple methods have been developed to estimate these parameters. Most of those methods such as maximum likelihood, subspace, and observer Kalman Identification require extensive offline processing and are not suitable for real time processing. One technique, which is suitable for real time processing, is the residual tuning method. Any mismodeling of the filter tuning parameters will result in a non-white sequence for the filter measurement residuals. The residual tuning technique uses this information to estimate corrections to those tuning parameters. The actual implementation results in a set of sequential equations that run in parallel with the Kalman filter. A. H. Jazwinski developed a specialized version of this technique for estimation of process noise. Equations for the estimation of the measurement noise have also been developed. These algorithms are used to estimate the process noise and measurement noise for the Wide Field Infrared Explorer star tracker and gyro.

Harman, Richard R.↗
NASA NTRS1999

An Adaptive Kalman Filter using a Simple Residual Tuning Method

One difficulty in using Kalman filters in real world situations is the selection of the correct process noise, measurement noise, and initial state estimate and covariance. These parameters are commonly referred to as tuning parameters. Multiple methods have been developed to estimate these parameters. Most of those methods such as maximum likelihood, subspace, and observer Kalman Identification require extensive offline processing and are not suitable for real time processing. One technique, which is suitable for real time processing, is the residual tuning method. Any mismodeling of the filter tuning parameters will result in a non-white sequence for the filter measurement residuals. The residual tuning technique uses this information to estimate corrections to those tuning parameters. The actual implementation results in a set of sequential equations that run in parallel with the Kalman filter. Equations for the estimation of the measurement noise have also been developed. These algorithms are used to estimate the process noise and measurement noise for the Wide Field Infrared Explorer star tracker and gyro.

Harman, Richard R.↗
DOE OSTI2021

Implicit highly-coupled single-ion Hall-MHD formulation for hybrid particle-in-cell codes

The rudiments of a particle-based single-fluid two-temperature magnetohydrodynamic (MHD) algorithm have been outlined in Thoma et al. (2013). The extension of this algorithm to include the effect of Hall physics is described in this paper. An implicit leapfrog version of the algorithm, which allows timesteps large compared to the resistive decay time and other relevant timescales, has recently been added to a hybrid particle-in-cell code. In standard MHD the Hall term in the generalized Ohm’s law can often be neglected when the Hall parameter is small. This term must, however, be retained in regimes where it is non-negligible. The retention of displacement current in Maxwell’s equations avoids the numerical difficulties associated with the whistler mode, which are encountered in standard explicit Hall-MHD codes, and allows the algorithm to be incorporated into hybrid particle-in-cell codes, for which particles may migrate from a kinetic to fluid to MHD description based upon local ambient plasma conditions. A highly-coupled implicit Hall-MHD formalism is presented, in which displacement current can either be retained or neglected. Even when displacement current is neglected, the highly-coupled implicit formalism avoids the restrictive timesteps for the whistler mode in explicit Hall-MHD codes. A comparison of numerical and analytic dispersion analysis demonstrates the feasibility of this approach and establishes relevant constraints to assure numerical stability. The implementation of the algorithm is described, and test simulation results in 1D and 2D in both linear and nonlinear regimes are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Assessing the difficulty of capturing the distribution function of neutrinos in neutron star merger simulations

The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Development and verification of local/global analysis techniques for laminated composites

Analysis and design methods for laminated composite materials have been the subject of considerable research over the past 20 years, and are currently well developed. In performing the detailed three-dimensional analyses which are often required in proximity to discontinuities, however, analysts often encounter difficulties due to large models. Even with the current availability of powerful computers, models which are too large to run, either from a resource or time standpoint, are often required. There are several approaches which can permit such analyses, including substructuring, use of superelements or transition elements, and the global/local approach. This effort is based on the so-called zoom technique to global/local analysis, where a global analysis is run, with the results of that analysis applied to a smaller region as boundary conditions, in as many iterations as is required to attain an analysis of the desired region. Before beginning the global/local analyses, it was necessary to evaluate the accuracy of the three-dimensional elements currently implemented in the Computational Structural Mechanics (CSM) Testbed. It was also desired to install, using the Experimental Element Capability, a number of displacement formulation elements which have well known behavior when used for analysis of laminated composites.

Griffin, O. Hayden, Jr.↗
NASA NTRS2006

Vivid Motor Imagery as an Adaptation Method for Head Turns on a Short-Arm Centrifuge

Artificial gravity (AG) has been proposed as a potential countermeasure to the debilitating physiological effects of long duration space flight. The most economical means of implementing AG may be through the use of a short-radius (2m or less) centrifuge. For such a device to produce gravitational forces comparable to those on earth requires rotation rates in excess of 20 revolutions per minute (rpm). Head turns made out of the plane of rotation at these rates, as may be necessary if exercise is combined with AG, result in cross-coupled stimuli (CCS) that cause adverse side effects including motion sickness, illusory sensations of motion, and inappropriate eye movements. Recent studies indicate that people can adapt to CCS and reduce these side effects by making multiple head turns during centrifuge sessions conducted over consecutive days. However, about 25% of the volunteers for these studies have difficulty tolerating the CCS adaptation paradigm and often drop out due to motion sickness symptoms. The goal of this investigation was to determine whether vivid motor imagery could be used as a pseudostimulus for adapting subjects to this unique environment. Twenty four healthy human subjects (14 males, 10 females), ranging in age from 21 to 48 years (mean 33, sd 7 years) took part in this study. The experimental stimuli were produced using the NASA JSC short-arm centrifuge (SAC). Subjects were oriented supinely on this device with the nose pointed toward the ceiling and head centered on the axis of rotation. Thus, centrifuge rotation was in the body roll plane. After ramp-up the SAC rotated clockwise at a constant rate of 23 rpm, producing a centrifugal force of approximately 1 g at the feet. Semicircular canal CCS were produced by having subjects make yaw head turns from the nose up (NU) position to the right ear down (RED) position and from RED to NU. Each head turn was completed in about one second, and a 30 second recovery period separated consecutive head movements. Participants were randomly assigned to one of three groups (n=8 per group): physical adapters (PA), mental adapters (MA), or a control group (CG). Each subject participated in a one hour test session on each of three consecutive days. Each test session consisted of an initial (preadaptation) period during which the subject performed six CCS maneuvers in the dark, followed by an adaptation period with internal lighting on the centrifuge, and a final (postadaptation) period during which six more CCS maneuvers were performed in the dark. For the PA group, the adaptation period consisted of performing 30 additional CCS maneuvers in the light. For the MA and CG group the centrifuge was ramped down to 0 rpm after the pre-adaptation period and ramped back up to 23 rpm before the post-adaptation period. For the both of these groups, the adaptation period consisted of making 30 CCS maneuvers in the light with the centrifuge stationary (so no cross-coupling occurred). MA group subjects were instructed to vividly imagine the provocative sensations produced by the preadaptation CCS maneuvers in terms of magnitude, duration, and direction of illusory body tilt, as well as any accompanying levels of motion sickness. CG group subjects were asked to answer low imagery content questions (trivial pursuit) during each adaptation period head turn. During the 30 second recovery following each head turn, psychophysical data were collected including self reports of motion sickness, magnitude and direction estimates of illusory body tilt, and the overall duration of these sensations. A multilevel mixed effects linear regression analysis performed on all response variables indicated that all three groups experienced some psychophysical adaptation across the three test sessions. For illusory tilt magnitude, the PA group exhibited the most overall adaptation, followed by the MA group, and the CG group. The slopes of these adaptation trajectories by group over day were significantly diffent from one another. For the perceived duration of sensations, the CG group again exhibited the least amount of adaptation. However, the rates of adaptation of the PA and the MA groups were indistinguishable, suggesting that the imagined pseudostimulus appeared to be just as effective a means of adaptation as the actual stimulus. The MA group's rate of adaptation to motion sickness symptoms was also comparable to the PA group. The use of vivid motor imagery may be an effective method for adapting to the illusory sensations and motion sickness symptoms produced by cross-coupled stimuli. For space-based AG applications, this technique may prove quite useful in retaining astronauts considered highly susceptible to motion sickness as it reduces the number of actual CCS required to attain adaptation.

Newby, N. J.↗
DOE OSTI2026

Protection System Validation Using Post-Event Anomaly Classification with Machine Learning

Power system protection devices have transitioned over the past few decades from mechanical to analog devices, then to solid state and finally digital. Relays and their associated critical network of equipment have significantly increased in complexity. Even internally, relays have gained significant intricacy, with relatively simple overcurrent or differential functions now being assisted by a myriad of other functions. This is necessary as the grid becomes more complex, but it brings increased difficulty in monitoring and upkeep. Misoperation caused by improper relay settings or malicious actions is a constant challenge faced by all utilities. These improper settings can be difficult to identify and may require exhaustive post-mortem analysis, typically after a major outage event has already occurred. A mechanism is needed for monitoring the behavior of protection systems to validate that they act and perform as expected. This work presents a concept for a machine learning (ML) system capable of validating the performance of protection systems by classifying anomalous events and characterizing protection system responses based solely on available current and voltage measurements. As a first step in its development, an experimental dataset is generated, and a random forest model is implemented with high accuracy in distinguishing four power system scenarios.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Protection System Validation with Machine Learning Anomaly Classification

A poster for the Early Career Poster Session. Power system protection devices have transitioned over the past few decades from mechanical to analog devices, then to solid state and finally digital. Relays and their associated critical network of equipment have significantly increased in complexity. Even internally, relays have gained significant intricacy, with relatively simple overcurrent or differential functions now being assisted by a myriad of other functions. This is necessary as the grid becomes more complex, but it brings increased difficulty in monitoring and upkeep. Misoperation caused by improper relay settings or malicious actions is a constant challenge faced by all utilities. These improper settings can be difficult to identify and may require exhaustive post-mortem analysis, typically after a major outage event has already occurred. A mechanism is needed for monitoring the behavior of protection systems to validate that they act and perform as expected. This work presents a concept for a machine learning (ML) system capable of validating the performance of protection systems by classifying anomalous events and characterizing protection system responses based solely on available current and voltage measurements. As a first step in its development, an experimental dataset is generated, and a random forest model is implemented with high accuracy in distinguishing four power system scenarios.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Lessons Learned From Public Outreach: Effective Incorporation of the Stakeholder Voice - 18447

The licensing process for any facility associated with nuclear power or aspects of the nuclear fuel cycle can be expected to encounter intense stakeholder interest, and often significant challenges. Over the last decade, staff of the Center for Nuclear Waste Regulatory Analyses (CNWRA{sup R}) at Southwest Research Institute{sup R} (SwRI{sup R}) provided comprehensive support to the U.S. Nuclear Regulatory Commission (NRC) for public outreach efforts across a wide variety of nuclear activities. Many of these projects included public involvement and consultation activities for compliance with the National Environmental Policy Act (NEPA), including scoping, public meetings, responding to public comments on draft documents, and holding small meetings with affected local, state, and tribal governments. In addition, we have supported NRC public outreach related to nuclear facility licensing activities (including the proposed repository at Yucca Mountain), rulemaking, and regulatory guidance development. This experience has offered us a broad view of the challenges faced by industry, regulators, and stakeholders. In the process, we have observed some best practices that lead to effective stakeholder communication. Communication about complex technical issues is challenging, particularly for societally controversial topics, such as civilian use of nuclear power and management of radioactive waste. Most commonly, we have observed a lack of mutual trust and respect among industry, regulators, and stakeholders that creates an atmosphere ripe for suspicion and misunderstanding. Stakeholders have, in public comments, criticized the lack of appropriate venues to raise concerns in which they feel they can get meaningful responses from the nuclear industry and regulators. Conversely, industry faces difficulty in reaching compromise with stakeholders, who sometimes misinterpret or misrepresent factual information. Within this contentious environment, regulators are constrained by their statutory framework, even as they are challenged to be more responsive, transparent, and fair to both stakeholders and industry. In spite of these obstacles, there are examples where feedback from stakeholders was heard and utilized. Because of successful collection of appropriate feedback, regulators became more adept at understanding the full range of perspectives from both industry and stakeholders, and subsequently did a better job at effectively incorporating feedback and communicating their decisions and findings in a clear and understandable way. Ultimately, the success of these examples was based on a measure of mutual respect, which was reached through effective communication and implementation of some best practices. These best practices include (i) outreach planning that is abundantly mindful of the various stakeholder points of view, site history, regional characteristics, and probable future and related actions, combined with advance identification of key stakeholders; (ii) development and distribution of information that facilitates public participation; and (iii) use of avenues by which stakeholders can provide meaningful and effective comments. Illustrative examples include (i) how concerns raised by stakeholders were heard by an applicant and resulted in a revised project proposal; (ii) use of a variety of interaction formats that incorporated input across a broad spectrum of stakeholders; (iii) utilizing stakeholder feedback received during scoping; and (iv) how a facility underestimated the impact of 'outsider' stakeholder voices that were in contrast to the known local community. By moving past the viewpoint that public outreach is a 'necessary evil' and engaging only to check a procedural box, each participant can benefit from mutually satisfactory and improved outcomes. (author)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2013

Final Report - Regulatory Considerations for Adaptive Systems

This report documents the findings of a preliminary research study into new approaches to the software design assurance of adaptive systems. We suggest a methodology to overcome the software validation and verification difficulties posed by the underlying assumption of non-adaptive software in the requirementsbased- testing verification methods in RTCA/DO-178B and C. An analysis of the relevant RTCA/DO-178B and C objectives is presented showing the reasons for the difficulties that arise in showing satisfaction of the objectives and suggested additional means by which they could be satisfied. We suggest that the software design assurance problem for adaptive systems is principally one of developing correct and complete high level requirements and system level constraints that define the necessary system functional and safety properties to assure the safe use of adaptive systems. We show how analytical techniques such as model based design, mathematical modeling and formal or formal-like methods can be used to both validate the high level functional and safety requirements, establish necessary constraints and provide the verification evidence for the satisfaction of requirements and constraints that supplements conventional testing. Finally the report identifies the follow-on research topics needed to implement this methodology.

Wilkinson, Chris↗
DOE OSTI2018

nbodykit: An Open-source, Massively Parallel Toolkit for Large-scale Structure

We present nbodykit, an open-source, massively parallel Python toolkit for analyzing large-scale structure (LSS) data. Using Python bindings of the Message Passing Interface, we provide parallel implementations of many commonly used algorithms in LSS. nbodykit is both an interactive and scalable piece of scientific software, performing well in a supercomputing environment while still taking advantage of the interactive tools provided by the Python ecosystem. Existing functionality includes estimators of the power spectrum, two- and three-point correlation functions, a friends-of-friends grouping algorithm, mock catalog creation via the halo occupation distribution technique, and approximate N-body simulations via the FastPM scheme. The package also provides a set of distributed data containers, insulated from the algorithms themselves, that enables nbodykit to provide a unified treatment of both simulation and observational data sets. nbodykit can be easily deployed in a high-performance computing environment, overcoming some of the traditional difficulties of using Python on supercomputers. We provide performance benchmarks illustrating the scalability of the software. The modular, component-based approach of nbodykit allows researchers to easily build complex applications using its tools. The package is extensively documented at http://nbodykit.readthedocs.io, which also includes an interactive set of example recipes for new users to explore. As open-source software, we hope nbodykit provides a common framework for the community to use and develop in confronting the analysis challenges of future LSS surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Validation of Loci-Stream for Autogenous Pressurization of Cryogenic Propellant Tank

Autogenous pressurization of cryogenic propellant tanks eliminates the need to have an additional pressurant tank on the space vehicle, which is highly advantageous due to reduced vehicle mass and design complexity. Autogenous pressurization therefore is one of the key technologies for deep space exploration and long-term space missions. The complex interaction of thermal gradients, turbulence and phase change near the interface make the problem a challenging one to model. Nodal analysis tools and reduced order models are unable to capture the necessary physics. 3-D CFD analyses are necessary to fully characterize autogenous pressurization. CFD analyses pose their own difficulties. The requisite CFD tool to tackle this problem need to be modular with the ability to incorporate various physics models, efficient, and computationally scalable for simulating flight size tanks. NASA MSFC's Loci-Stream CFD tool along with the VOF module is a great candidate to fit this mold. We demonstrate our modeling approach and validation of Loci-Stream for predicting autogenous pressurization of a flight scale propellant tank in order for the solver to serve as a reliable design and analysis tool for NASA's CFM application needs. Liquid hydrogen tank pressurization tests carried out at the MSFC test stand 300 facilities provide reliable validation data for this purpose. These tests were modeled using the Loci-Stream solver with a newly implemented two-phase sharp interface treatment. We show that our modeling approach and CFD solver predict the autogenous pressurization phenomena satisfactorily, and document challenging aspects of modeling this problem.

cryogenic fluid management↗