Search NASA⌕ Search

SEARCH · Search NASA

Results for “SPECIAL RELATIVITY THEORY”

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.

268 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Greedy emulators for nuclear two-body scattering

Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. Here, in this paper, we develop, implement, and test two model-driven emulators based on the (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples (“snapshots”) for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires many fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.

Bayesian methods↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

Origin, Bulk Chemical Composition and Physical Structure of the Galilean Satellites of Jupiter: A Post-Galileo Analysis

The origin of Jupiter and the Galilean satellite system is examined in the light of the new data that has been obtained by the NASA Galileo Project. In particular, special attention is given to a theory of satellite origin which was put forward at the start of the Galileo Mission and on the basis of which several predictions have now been proven successful. These predictions concern the chemical composition of Jupiter's atmosphere and the physical structure of the satellites. According to the proposed theory of satellite origin, each of the Galilean satellites formed by chemical condensation and gravitational accumulation of solid grains within a concentric family of orbiting gas rings. These rings were cast off equatorially by the rotating proto-Jovian cloud (PJC which contracted gravitationally to form Jupiter some 4 1/2 billion years ago. The PJC formed from the gas and grains left over from the gas ring that had been shed at Jupiter's orbit by the contracting proto-solar cloud (PSC Supersonic turbulent convection provides the means for shedding discrete gas rings. The temperatures T (sub n) of the system of gas rings shed by the PSC and PJC vary with their respective mean orbital radii R (sub n) (n = 0, 1,2,...) according as T (sub n) proportional to R (sub n) (exp -0.9). If the planet Mercury condenses at 1640 K, so accounting for the high density of that planet via a process of chemical fractionation between iron and silicates, then T (sub n) at Jupiter's orbit is 158 K. Only 35% of the water vapour condenses out. Thus fractionation between rock and ice, together with an enhancement in the abundance of solids relative to gas which takes place through gravitational sedimentation of solids onto the mean orbit of the gas ring, ensures nearly equal proportions of rock and ice in each of Ganymede and Callisto. Io and Europa condense above the H20 ice point and consist solely of hydrated rock (h-rock). The Ganymedan condensate consists of h-rock and H20 ice. For Callisto, NH3 ice makes up -5% of the condensate mass next to h-rock (approximately 50%) and H2O ice (approximately 45%). Detailed thermal and structural models for each of Europa, Ganymede and Callisto are constructed on the basis of the above initial bulk chemical compositions. For Europa (E), a predicted 2-zone model consisting of a dehydrated rock core of mass 0.912 M (sub E) and a 150 km thick frozen mantle of salty H2O yields a moment-of-inertia coefficient which matches the Galileo Orbiter gravity measurement. For Ganymede (G), a 3-zone model possessing an inner core of solid FeS and mass approximately 0.116 M (sub G), and an outer H2O ice mantle of mass approximately 0.502 M (sub G) is needed to explain the gravity data. Ganymede's native magnetic field was formed by thermoremanent magnetization of Fe3O4. A new Callisto (C) model is proposed consisting of a core of mass 0.826 M (sub C) containing a uniform mixture of h-rock (60% by mass) and H2O and NH3 ices, and capped by a mantle of pure ice. This model may have the capacity to yield a thin layer of liquid NH3 (raised dot) 2H2O at the core boundary, in line with Galileo's discovery of an induced magnetic field.

Prentice, A. J. R.↗
DOE OSTI2024

Foundations of Molecular 'Isotomics'

The naturally occurring rare isotopes are versions of common elements, such as hydrogen, carbon and oxygen, that contain a larger than usual number of neutrons in their atomic nuclei and therefore are higher in mass than the common atoms of that element. Isotopes exist for most elements and are found in most natural and synthetic materials, but are uneven in their distribution because chemical and physical processes are isotope-selective (e.g., a chemical reaction may proceed more rapidly for one isotope than for another). For this reason, abundances of isotopes in a material of interest can provide a record, or ‘signature’ of various features of that material’s origin and history. These signatures have been used in the geo, life, chemical and physical sciences in a wide variety of ways over close to 8 decades. However, many such applications struggle to reach unique interpretations of isotopic data because multiple factors combine to control a given sample’s overall isotopic content. That is, the factors controlling isotopic content are too numerous and complex to fully constrain from a simple measurement of a material’s isotope abundances. However, the distribution of isotopes within materials, at molecular scales potentially provides a vastly larger number and diversity of constraints on the chemical and physical processes that comprise a material’s history. The rare isotopes may be concentrated into one atomic position in a molecule relative to another, some proportion of molecules in a sample may contain two or more rare isotopes, and those multiply-isotope-substituted forms of molecules may also have uneven distributions of those isotopes across individual atomic sites. For these reasons, even small, seemingly simple molecules, such as sugars, amino acids or drug compounds, actually exist in a vast number of isotopically unique forms (often millions or more), and each one of those forms is in some sense an independent ‘vote’ on that sample’s history. This project has focused on opening this rich archive of information by enabling the creation of routinely and widely applicable ways of measuring and interpreting isotopic structures of molecules. This work has included the development of core technologies and analytical methods, advancing fundamental understanding of the physical and chemical properties of isotopic versions of molecules, and conducting proof of concept studies of illustrative geochemical, cosmochemical and forensic problems in order to show how these technologies, methods and principles come together to solve problems in new ways. A key to the success of this project was the adaptation of ‘Fourier transform mass spectrometry’ (FTMS) to the task of precisely measuring proportions of the rare, naturally occurring isotopic forms of molecules. FTMS is a highly specialized form of mass spectrometry that traps ions within magnetic or electrostatic cavities and, effectively, ‘listens’ (through registering of subtle electrical signals) to the harmonic signals they make while rapidly orbiting within those cavities. These signals have periods that are a function of their mass and strength (or ‘loudness’) that is proportional to their abundances. Thus, these signals constrain relative amounts of molecules that differ in their mass due to various isotopic substitutions. This technology has been essential to the identification of organic molecules in the life, chemical and environmental sciences for over 4 decades, but generally has lacked the control, stability and precision to meaningfully measure rare isotope forms of molecules. This project’s most fundamental contribution has been to modify FTMS, both in terms of hardware and methods, to enable such measurements. The raw data of molecular isotopic structure is tremendously voluminous and complex, so another important activity of this project has been developing the theoretical and data-science tools needed to interpret the data generated by this new form of isotopic measurement. A particularly challenging part of this task has been predicting molecular isotopic structure, as only through the comparison of measurements with predictions can we make progress on hypothesis driven research questions. We have attacked this this prediction task through a combination of first-principles chemical-physics models of the effects of isotope substitution on molecule properties and data-science models that permit us to generalize that chemical physics to cases that have not yet been studied by detailed chemical physics theory. The proof of concept applications we have pursued over the course of this study include biological reactions of amino acids and other biomolecules, non-biological synthesis of organic molecules in extra-terrestrial settings such as meteorites, petroleum geoscience questions concerning the origin and evolution of natural gas, oil and kerogen compounds, and forensic questions such as the sourcing of chemical weapons. The successes of these applications have laid the groundwork for the next phase of this field’s development, which will include larger scale and more ambitious studies of molecular isotopic structure as a means of diagnosing human diseases, such as cancer, and reconstructing detailed interpretations of the origin and evolution of organic molecules in modern and geological environments.

Cesar, Jaime↗
DOE OSTI2023

Synthesis and Observation of Emergent Phenomena in Epitaxial Heusler Compound Heterostructures

The proposal was on the synthesis and observation of emergent phenomena in epitaxial Heusler compound heterostructures. The large range of properties and number of Heusler compounds opens up a wide number of potential compounds that will exhibit emergent phenomena. The similarity, large range of relatively inexpensive, large area, high crystal quality, III-V bulk substrates, lattice parameters and the ability to tune the lattice parameters through ternary or quaternary III-V compound semiconductor epitaxial growth, makes III-V semiconductors an ideal choice for substrates for epitaxial growth of Heusler compounds. A number of half Heusler compounds have been predicted to exhibit band inversion, making them topological and are therefore expected to exhibit spin-momentum locked topological surface states with linear dispersion. Others are predicted to be semimetals with Weyl points and others semiconducting and magnetic. During the course of this grant, emphasis has been on investigating Heusler compounds with emergent phenomena and demonstrating the ability to tune their properties through alloying and strain. We have grown toplogical semimetal (PtLuSb, PtMnBi), Weyl (Co 2 MnAl, Co 2 TiGe), half metal (PtMnSb, Co 2 MnSi, Co 2 MnAl x Si 1-x , Co2FeAl), and semiconducting (CoTiSb, NiTiSn) and tuned their properties through alloying and epitaxial strain. We also investigated the closely related materials of rare-earth monopnictide, some of which have also been predicted to be topological. During the attempts to grow the PtMnBi, it was discovered that Bi, another predicted topological material when ultrathin, could be grown epitaxially on InSb, results for which are also reported here. The main focus for this effort has been on using variable photon energy and spin-dependent angle resolved photoemission (ARPES) to determine bulk band structure and surface states of pristine epitaxial films grown on III-V semiconductor and MgO substrates and correlate results with theory and transport measurements. Theory has been critical to interpretation of experimental results and has been essential in guiding experiments. The research benefited from several strong collaborations between the PIs and the beamline scientists at the Advanced Light Source at Lawrence Berkeley Laboratory, the Stanford Linear Accelerator Center (SLAC) at Stanford and at the Max Lab at Lund University in Sweden. The strong experiment - theory collaboration between the PI’s groups, the Palmstrøm group at UCSB and the Janotti group at the University of Delaware, has been critical for interpreting the experimental ARPES and magnetotransport measurements results and making predictions to guide experiments. Weekly interactive Zoom meetings made this work well. A collaboration between the Palmstrøm group and Dr. Alexei Fedorov at the Advanced Light Source (ALS) resulted in significant modifications to his end chamber to accommodate the vacuum suitcase that was designed and constructed in the Palmstrøm group at UCSB. In collaboration with beamline scientists, Drs. Makoto Hashimoto and Donghui Lu at SLAC, Palmstrøm made modifications to the vacuum suitcase and developed special sample holders that allowed samples to be grown in the Palmstrøm MBE systems at UCSB and transported in the UHV vacuum suitcase to SLAC for ARPES measurements. The development of the vacuum suitcase was essential for this grant as it has allowed variable photon energies to be used to identify surface versus bulk states on samples that could not be capped and decapped using As- or Sb-capping layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Recurrent features of amplitudes in planar $\mathcal{N}$ = 4 super Yang-Mills theory

The planar three-gluon form factor for the chiral stress tensor operator in planar maximally supersymmetric Yang-Mills theory is an analog of the Higgs-to-three-gluon scattering amplitude in QCD. The amplitude (symbol) bootstrap program has provided a wealth of high-loop perturbative data about this form factor, with results up to eight loops available. The symbol of the form factor at L loops is given by words of length 2L in six letters with associated integer coefficients. In this paper, we analyze this data, describing patterns of zero coefficients and relations between coefficients. We find many sequences of words whose coefficients are given by closed-form expressions which we expect to be valid at any loop order. Moreover, motivated by our previous machine-learning analysis, we identify simple recursion relations that relate the coefficient of a word to the coefficients of particular lower-loop words. These results open an exciting door for understanding scattering amplitudes at all loop orders.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Enhancement of EM Pump Performance through Modeling and Testing

Liquid sodium has seen development as a nuclear reactor coolant in the United States for over 50 years. While sodium has many desirable characteristics as a heat transfer fluid, its reactivity poses numerous engineering challenges for even the most basic plant components like pumps. Despite the challenges, mechanical pumps have seen success in sodium service albeit at the expense of simplicity and flexibility. In fact, significant failures of these pumps have been recorded which can be attributed to their complexity. Notwithstanding the complexity, mechanical pumps are still widely proposed for the main circuit pump of sodium cooled reactors. However, an alternative to mechanical pumps are Electromagnetic Induction Pumps (EMIPs) which use the principle of the Lorentz Force. EMIPs avoid mechanically moving parts in contact with sodium which simplifies the design by avoiding mechanical wear and seals. Therefore, these pumps are generally safer and more reliable than their mechanical counterparts. That being said, EMIPs suffer from relatively low efficiencies that range from 5% to 45%. Some of the inefficiencies are inherent to the EMIPs pumps; namely, inefficiencies arise in the form of resistive heating losses in the fluid itself. However, other inefficiencies arise from Edge Effects which include Finite Length and Finite Width effects. Due to the significant safety advantages EMIPs have compared to mechanical pumps, there is a strong motivation to integrate these technologies into the design of Generation IV liquid-metal-cooled nuclear power plants. However, EMIPs have low efficiencies partly attributed to Edge Effects. Much of the literature focusing on Edge Effects is dated or simply does not exist in the public domain. More recently, numerical models of EMIPs have shown promising results. However, these models are often custom solvers or use complex software suits which limit their application. This report describes the work completed to understand Edge Effects in two EMIPs designs. These designs are the Annular Linear Induction Pump (ALIP) and Permanent Magnet Induction Pump (PMIP) also known as a Moving Magnet Pump (MMP). The methodology developed during this project aim at measuring key performance parameters (pressure-flowrate and efficiency-slip performance) using high temperature liquid-sodium. Baseline experimental measurements were then used to verify the accuracy of numerical models developed in the software suit ANSYS/FLUENT. On the experimental side, instrumentation was developed to characterize the pressure-flowrate and efficiency-slip performance of both the ALIP and PMIP in high-temperature liquid-sodium. Additionally, instrumentation was developed to measure potential output pressure pulsations and the magnetic field distributions. Moreover, a specialized ALIP was constructed to study the impacts magnetic field tapering, also known as coil grading, on the Finite Length Effect. Lastly, experimental performance curves of the PMIP were collected at 300 oC and 600 oC and compared to analytical theory which accounted for Finite Width Effects. These results showed that the Finite Width Effect and non-uniform magnetic field in a PMIP impact the performance. On the numerical side, models were developed using the software suit FLUENT/ANSYS. These methods were initially compared against a simple MHD flow in a Permanent Magnet Electromagnet Flowmeter (PMEMFM). The results were significant in that they showed agreeance between the numerical model and demonstrated that PMEMFM’s could be calibrated numerically. These modeling methods were then applied to a prototypic PMIP. Results showed that this simple model provided reasonably accurate results of the PMIP p-Q performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRSDate not supplied

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗
NASA NTRS2015

A Vehicle Management End-to-End Testing and Analysis Platform for Validation of Mission and Fault Management Algorithms to Reduce Risk for NASA's Space Launch System

The engineering development of the new Space Launch System (SLS) launch vehicle requires cross discipline teams with extensive knowledge of launch vehicle subsystems, information theory, and autonomous algorithms dealing with all operations from pre-launch through on orbit operations. The characteristics of these spacecraft systems must be matched with the autonomous algorithm monitoring and mitigation capabilities for accurate control and response to abnormal conditions throughout all vehicle mission flight phases, including precipitating safing actions and crew aborts. This presents a large and complex system engineering challenge, which is being addressed in part by focusing on the specific subsystems involved in the handling of off-nominal mission and fault tolerance with response management. Using traditional model based system and software engineering design principles from the Unified Modeling Language (UML) and Systems Modeling Language (SysML), the Mission and Fault Management (M&FM) algorithms for the vehicle are crafted and vetted in specialized Integrated Development Teams (IDTs) composed of multiple development disciplines such as Systems Engineering (SE), Flight Software (FSW), Safety and Mission Assurance (S&MA) and the major subsystems and vehicle elements such as Main Propulsion Systems (MPS), boosters, avionics, Guidance, Navigation, and Control (GNC), Thrust Vector Control (TVC), and liquid engines. These model based algorithms and their development lifecycle from inception through Flight Software certification are an important focus of this development effort to further insure reliable detection and response to off-nominal vehicle states during all phases of vehicle operation from pre-launch through end of flight. NASA formed a dedicated M&FM team for addressing fault management early in the development lifecycle for the SLS initiative. As part of the development of the M&FM capabilities, this team has developed a dedicated testbed that integrates specific M&FM algorithms, specialized nominal and off-nominal test cases, and vendor-supplied physics-based launch vehicle subsystem models. Additionally, the team has developed processes for implementing and validating these algorithms for concept validation and risk reduction for the SLS program. The flexibility of the Vehicle Management End-to-end Testbed (VMET) enables thorough testing of the M&FM algorithms by providing configurable suites of both nominal and off-nominal test cases to validate the developed algorithms utilizing actual subsystem models such as MPS. The intent of VMET is to validate the M&FM algorithms and substantiate them with performance baselines for each of the target vehicle subsystems in an independent platform exterior to the flight software development infrastructure and its related testing entities. In any software development process there is inherent risk in the interpretation and implementation of concepts into software through requirements and test cases into flight software compounded with potential human errors throughout the development lifecycle. Risk reduction is addressed by the M&FM analysis group working with other organizations such as S&MA, Structures and Environments, GNC, Orion, the Crew Office, Flight Operations, and Ground Operations by assessing performance of the M&FM algorithms in terms of their ability to reduce Loss of Mission and Loss of Crew probabilities. In addition, through state machine and diagnostic modeling, analysis efforts investigate a broader suite of failure effects and associated detection and responses that can be tested in VMET to ensure that failures can be detected, and confirm that responses do not create additional risks or cause undesired states through interactive dynamic effects with other algorithms and systems. VMET further contributes to risk reduction by prototyping and exercising the M&FM algorithms early in their implementation and without any inherent hindrances such as meeting FSW processor scheduling constraints due to their target platform - ARINC 653 partitioned OS, resource limitations, and other factors related to integration with other subsystems not directly involved with M&FM such as telemetry packing and processing. The baseline plan for use of VMET encompasses testing the original M&FM algorithms coded in the same C++ language and state machine architectural concepts as that used by Flight Software. This enables the development of performance standards and test cases to characterize the M&FM algorithms and sets a benchmark from which to measure the effectiveness of M&FM algorithms performance in the FSW development and test processes.

Trevino, Luis↗
NASA NTRS2015

A Vehicle Management End-to-End Testing and Analysis Platform for Validation of Mission and Fault Management Algorithms to Reduce Risk for NASA's Space Launch System

The development of the Space Launch System (SLS) launch vehicle requires cross discipline teams with extensive knowledge of launch vehicle subsystems, information theory, and autonomous algorithms dealing with all operations from pre-launch through on orbit operations. The characteristics of these systems must be matched with the autonomous algorithm monitoring and mitigation capabilities for accurate control and response to abnormal conditions throughout all vehicle mission flight phases, including precipitating safing actions and crew aborts. This presents a large complex systems engineering challenge being addressed in part by focusing on the specific subsystems handling of off-nominal mission and fault tolerance. Using traditional model based system and software engineering design principles from the Unified Modeling Language (UML), the Mission and Fault Management (M&FM) algorithms are crafted and vetted in specialized Integrated Development Teams composed of multiple development disciplines. NASA also has formed an M&FM team for addressing fault management early in the development lifecycle. This team has developed a dedicated Vehicle Management End-to-End Testbed (VMET) that integrates specific M&FM algorithms, specialized nominal and off-nominal test cases, and vendor-supplied physics-based launch vehicle subsystem models. The flexibility of VMET enables thorough testing of the M&FM algorithms by providing configurable suites of both nominal and off-nominal test cases to validate the algorithms utilizing actual subsystem models. The intent is to validate the algorithms and substantiate them with performance baselines for each of the vehicle subsystems in an independent platform exterior to flight software test processes. In any software development process there is inherent risk in the interpretation and implementation of concepts into software through requirements and test processes. Risk reduction is addressed by working with other organizations such as S&MA, Structures and Environments, GNC, Orion, the Crew Office, Flight Operations, and Ground Operations by assessing performance of the M&FM algorithms in terms of their ability to reduce Loss of Mission and Loss of Crew probabilities. In addition, through state machine and diagnostic modeling, analysis efforts investigate a broader suite of failure effects and detection and responses that can be tested in VMET and confirm that responses do not create additional risks or cause undesired states through interactive dynamic effects with other algorithms and systems. VMET further contributes to risk reduction by prototyping and exercising the M&FM algorithms early in their implementation and without any inherent hindrances such as meeting FSW processor scheduling constraints due to their target platform - ARINC 653 partitioned OS, resource limitations, and other factors related to integration with other subsystems not directly involved with M&FM. The plan for VMET encompasses testing the original M&FM algorithms coded in the same C++ language and state machine architectural concepts as that used by Flight Software. This enables the development of performance standards and test cases to characterize the M&FM algorithms and sets a benchmark from which to measure the effectiveness of M&FM algorithms performance in the FSW development and test processes. This paper is outlined in a systematic fashion analogous to a lifecycle process flow for engineering development of algorithms into software and testing. Section I describes the NASA SLS M&FM context, presenting the current infrastructure, leading principles, methods, and participants. Section II defines the testing philosophy of the M&FM algorithms as related to VMET followed by section III, which presents the modeling methods of the algorithms to be tested and validated in VMET. Its details are then further presented in section IV followed by Section V presenting integration, test status, and state analysis. Finally, section VI addresses the summary and forward directions followed by the appendices presenting relevant information on terminology and documentation.

Trevino, Luis↗
DOE OSTI2024

Off-shell Partition Functions in 3d Gravity

We explore three-dimensional gravity with negative cosmological constant via canonical quantization. We focus on chiral gravity which is related to a single copy of PSL(2,R) Chern-Simons theory and is simpler to treat in canonical quantization. Its phase space for an initial value surface Σ is given by the appropriate moduli space of Riemann surfaces. We use geometric quantization to compute partition functions of chiral gravity on three-manifolds of the form Σ×S 1 , where Σ can have asymptotic boundaries. Most of these topologies do not admit a classical solution and are thus not amenable to a direct semiclassical path integral computation. We use an index theorem that expresses the partition function as an integral of characteristic classes over phase space. In the presence of n asymptotic boundaries, we use techniques from equivariant cohomology to localize the integral to a finite-dimensional integral over $\overline{M}$ g,n , which we evaluate in low genus cases. Higher genus partition functions quickly become complicated since they depend in an oscillatory way on Newton's constant. There is a precise sense in which one can isolate the non-oscillatory part which we call the fake partition function. We establish that there is a topological recursion that computes the fake partition functions for arbitrary Riemann surfaces Σ. As a result, there is a scaling limit in which the model reduces to JT gravity and our methods give a novel way to compute JT partition functions via equivariant localization.

Classical and Quantum Gravity↗
NASA NTRS1995

Kinetic structure of intermediate shocks: Implications for the magnetopause

A general study of the structure and stability of intermediate shocks (IS) in an isotropic plasma is presented using a hybrid as well as a resistive Hall MHD code. Special emphasis is put on the question of whether the rotational layers observed at the magnetopause can be intermediate shocks. The shocks are formed dynamically by the interaction between a flowing plasma and a stationary piston. Coplanar ISs (both strong and weak) are found to be stable in a collisionless plasma. The existence of slow shocks in a high beta plasma is also established for the first time. Noncoplanar ISs are found to be time-dependent, evolving toward a rotational discontinuity (RD) after some characteristic time tau which can be quite long (1000 Omega(exp -1), where Omega is the ion gyrofrequency). The value tau is larger the closer the rotation angle is to 180 deg. Rotations larger than 180 deg are found to be unstable, decaying into a state of minimum shear (i.e., rotation angle less than 180 deg). There are various length scales associated with an IS in the kinetic regime. The shortest scale is found to be the length scale over which rotation of the transverse component of the magnetic field takes place. This scale can have a half width as small as one ion inertial length (c/omega(sub p)) for electron sense rotations and 3c/omega(sub p) for ion sense rotations, for an upstream ion beta of unity. Both of these scales are consistent with the observed thickness at the magnetopause and identical to the corresponding RD scales. A detailed study of the mode conversion of the Alfven ion cyclotron waves (A/IC) waves across both slow and intermediate shocks and the resulting downstream wave spectrum are presented. The possibility that the large number of relfected ions observed at the magnetopause may be due to the presence of strong ISs is considered. The identification of strong ISs and their distinction from RDs should be possible in observations due to significant differences that exist between jump conditions and overall structure of the two discontinuities. The jumps in the plasma parameters across a weak IS are typically small. This together with the fact that the weak ISs and RDs have very similar thickness and other overall properties makes the distinction between weak ISs and RDs in the observations largely inconsequential. However, at large noncoplanarity angles the weak IS approaches the RD limit in a relatively short time (approximately less than 100 Omega (exp -1)). Thus, magnetopause rotations with large noncoplanarity angles are most likely either RDs or strong ISs. Finally, direct comparisons between fluid (resistive Hall MHD) and kinetic simulations show that fluid theory is not applicable to study of ISs in a collisionless plasma.

Karimabadi, H.↗
DOE OSTI2021

A trait-based framework for linking microbial communities with carbon transformations under precipitation change

Droughts are common throughout the world. As the climate changes, droughts may become more frequent and intense. Still, scientists are uncertain about how drought will affect the natural world, particularly the bacteria, fungi, and other microbes that live in soils. These tiny life forms are crucial because they control the Earth’s flows of carbon and essential nutrients. Researchers at the University of California, Irvine, and Lawrence Berkeley National Laboratory teamed up to study how the microbiome, or collection of bacteria and fungi in the soil, deals with drought. Since 2007, the researchers have used shelters with retractable roofs to prevent nearly half of normal rainfall from reaching grass and shrub ecosystems, and their soil microbiomes, in Southern California. The study team discovered that microbes have some clever tricks up their sleeve for surviving drought. When growing on dead grass as a food source, microbes ramp up production of specialized chemicals called osmolytes that keep their cells from drying out. But microbes growing on dead shrubs face another problem. Unlike grass, shrub leaves are a lousy food source. To digest and consume shrub leaves, microbes have to exude more enzymes, which act like biochemical chef knives that chop complex leaf molecules into bite-sized pieces. Carbon is the coin of the microbial realm, earned via enzyme action or slurping up dead plant juices. Microbes growing under normal conditions on tasty dead grass have it easy: they can spend most of their carbon currency on growth. With drought, life gets harder as microbes need to pay up for osmolytes and ramp down their growth. It gets worse with shrub leaves because microbes have to multi-task among growth, enzyme secretion, and osmolyte production. When drought hits, microbes on shrub leaves forgo the osmolytes, probably because losing the carbon revenue from enzyme investment would be a deal-breaker for survival. The next question tackled by the researchers asked how the genes controlling microbial lifestyles sort out across the tree of life. Most microbiologists thought these lifestyles would correspond to rather large branches on the tree. But the study team found that in fact, very closely related microbes differ in important and interesting ways. For example, bacteria that have nearly identical housekeeping genes respond distinctively to warming, rainfall, and plant chemistry. As a result, soil microbiomes are much richer in diversity than originally thought. And studying microbial diversity in much greater detail could open up many more possibilities for how microbial life deals with changing environmental conditions. Trying to understand microbial life without the fine details of genetic diversity is like trying to stream Netflix on a dial-up connection. By looking across the landscape, the researchers revealed that microbial diversity is absolutely critical for maintaining the planet’s flows of carbon and nutrients. The study team designed a new technology—microbial cages—for transplanting intact microbiomes. With the cages, the researchers could move microbiomes to novel environments and compare their ability to cycle carbon and nutrients. In some cases, performance tailed off when microbiomes found themselves in a new environment, but in other cases, performance rivaled or even exceeded that of the resident microbiome. And even the low performers eventually caught up to the native microbiomes if given sufficient time. These findings mean that microbiomes—at least in Southern California—may be resilient to climate shifts due to a high diversity of microbial lifestyles. Coping with heat and drought may literally be in their DNA. The last piece of the research puzzle, and a “Holy Grail” for microbial ecologists, is to forecast the behavior of diverse microbiomes. To meet this challenge, the study team developed new theory and computer models. For the first time, these models account for hundreds of different microbes and how their intricate lifestyles cope with drought. The models are starting to connect the tiniest microbes with the global cycles that sustain Earth’s farms, fields, and forests. With that connection, it will be easier for society to plan for a world with more droughts and other climate disruptions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Uranium-Containing and Thorium-Containing Anions Studied by Photoelectron Spectroscopy

An in-depth knowledge of actinide chemistry is fundamental to many aspects of nuclear science and technology, including the synthesis and processing of materials and the remediation of waste disposal sites. Among the actinides, the chemical bonding behaviors of actinium and thorium resemble those of the transition metals; the 5f-electrons of protactinium, uranium, neptunium, and plutonium often play important roles in their bonding; and among the still heavier elements, their bonding tends to mimic the lanthanide elements in terms of electron shielding and their f-electron contributions. Bonding that involves 5f-electrons, however, is especially important, in part because of the significance of uranium and plutonium, but also because these elements are among the few where f-electron participation in bonding is relatively common. This work focused on studying uranium-containing and thorium-containing anions in the gas phase using negative ion photoelectron spectroscopy. Since this technique directly probed valence electrons, it was uniquely positioned to address open questions regarding molecular bonding and electron configurations. A particularly important issue concerned how bonding in actinide-containing molecules was affected by modifications to their actinide atoms’ environment, i.e., due to their interaction with ligands. A closely related question was how actinide atoms’ suborbitals were qualitatively reordered and their energies quantitatively shifted as a result of their ligated environments. These were especially relevant issues in regard to uranium due to it having multiple possible oxidation states (OS) and the potential for 5f electron participation in bonding. The effects of ligands on oxidation states and 5f-orbital energies in uranium bonding was expected to be pronounced. Both ligands and excess electrons were seen as probes of actinide atoms within actinide-containing molecules. Our strategy for advancing knowledge of chemical bonding in the actinide-containing species utilized the synergy between experiments and theory, where in some cases experimental results validated theory and where in others computational results assisted in interpreting experiments. Calculations on actinide systems are terrifically challenging due to large spin-orbit interactions, relativistic effects, and just the sheer number of electrons involved. Even in the simplest species, e.g., U and U2, the most sophisticated, modern calculations carried out by the most experienced theorists often only approximate experimentally-measured values, such as electron affinities. For theory to provide confident predictions that can be used to solve real problems it needed an iterative and ultimately corrective mechanism by which its methods can develop further. Experiments can be used to identify when theory has failed; whereupon the subsequent process of using the experiment-theory interplay can be used to find the cause of the failure. Upon fixing it in one case, different test species can be proposed and studied by the experiment-theory combination to determine whether the problem has been corrected. Thus, experiments not only measure the values of molecular properties, they also provide navigational 3 beacons that keep computations off the reefs in an otherwise dark sea with few reference points. Experimental measurements in the actinide field are not only important, they are in actuality essential to computational progress. While it was not always possible to compare the theoreticallydetermined quantity of interest directly with the same experimentally-measured observable, it was usually possible to compare consequential properties that are both calculable and measurable. In the work completed here electron affinities and electronic state spacings were often sensitive consequential parameters. Reasonable agreement between measured and computational values signaled that a calculation that was very likely to be on-track. We had established collaborative relationships with five computational groups, all of which have expertise in computational actinide chemistry. Their PI’s are L. Cheng, D. Dixon, L. Gagliardi, K. Peterson, and B. Vlaisavljevich. Our close interaction with our theory partners led to us suggesting systems to them and them to us. This reciprocal interaction between our experimental and their computational results was among the most important strengths of this work and was a thread woven throughout. Even though anion photoelectron spectroscopic studies are conducted on anions, much of the information that they provide, pertains to the electronic structure of the neutral counterparts of those anions; among these are electron affinities and electronically excited state spacings. Our experimental tools included several specialized ion sources for forming the anionic species of interest, a mass spectrometer for identifying and mass-selecting them, and an anion photoelectron spectrometer for determining their electron affinities (EA) and characterizing the electronic states of the selected anions’ neutral counterparts. Anion photoelectron spectroscopy is conducted by crossing a mass-selected beam of anions with a fixed-frequency laser beam and energy-analyzing the resultant photodetached electrons. The photodetachment process is governed by the energyconserving relationship: hν = EBE + EKE, where hν is the photon’s energy, EBE is the electron binding (photodetachment transition) energy, and EKE is the electron’s kinetic energy. In our apparatus mass-selection is accomplished via time-of-flight mass spectrometry (TOF-MS), electron energy analysis is achieved with either a magnetic bottle or by velocity mapped imaging. Photodetachment of electrons from anions is implemented via either Nd:YAG or excimer lasers. The photodetachment transition energy, i.e., the EBE, between the ground vibrational and electronic state of an anion and the ground vibrational and electronic state of that anion’s neutral counterpart is the adiabatic electron affinity (EA) of that neutral molecule. Likewise, photodetachment transitions between the ground vibrational and electronic state of an anion and the various electronically-excited states of that anion’s corresponding neutral map the electronic spectrum of that neutral species, i.e., the spectral spacings in the photoelectron spectrum are a mirror image of the neutral’s electronic spectrum. It was, of course, crucial to be able to form the anionic species of interest. There, we had a particularly broad field of anion sources from which to choose. These included several variants of pulsed laser vaporization (LV), laser photoemission, infrared desorption plus photoemission, pulsed arc discharge (PACIS), electrospray ionization (ESI), and Rydberg electron transfer (RET). Each of these anion sources were readily combined with, i.e., connected to, the anion photoelectron spectroscopic portion of our apparatus as described above. Among the sources that utilize lasers, visible light for LV sources as well as IR for desorption sources are provided by Nd:YAG lasers. Ultraviolet photons are provided by both Nd:YAG and excimer lasers, whereas the excitation wavelengths for RET experiments come from two Nd:YAG-pumped dye lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

LA17-V-InverseRadTransportProblem-PD2Kb. Final Project Report

In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗