Results for “Variational Monte Carlo” — Page 14

NASA NTRS2001

Reliability Analysis and Reliability-Based Design Optimization of Circular Composite Cylinders Under Axial Compression

This report describes the preliminary results of an investigation on component reliability analysis and reliability-based design optimization of thin-walled circular composite cylinders with average diameter and average length of 15 inches. Structural reliability is based on axial buckling strength of the cylinder. Both Monte Carlo simulation and First Order Reliability Method are considered for reliability analysis with the latter incorporated into the reliability-based structural optimization problem. To improve the efficiency of reliability sensitivity analysis and design optimization solution, the buckling strength of the cylinder is estimated using a second-order response surface model. The sensitivity of the reliability index with respect to the mean and standard deviation of each random variable is calculated and compared. The reliability index is found to be extremely sensitive to the applied load and elastic modulus of the material in the fiber direction. The cylinder diameter was found to have the third highest impact on the reliability index. Also the uncertainty in the applied load, captured by examining different values for its coefficient of variation, is found to have a large influence on cylinder reliability. The optimization problem for minimum weight is solved subject to a design constraint on element reliability index. The methodology, solution procedure and optimization results are included in this report.

Rais-Rohani, Masoud↗

NASA NTRS2004

Lunar and Planetary Science XXXV: Outer Solar System

The session 'Outer Solar System" inlcuded:Monte Carlo Modeling of [O I] 630 nm Auroral Emission on Io; The Detection of Iron Sulfide on Io; Io and Loki in 2003 as Seen from the Infrared Telescope Facility Using Mutual Satellite and Jupiter Occultations; Mapping of the Zamama-Thor Region of Io; First Solar System Results of the Spitzer Space Telescope; Mapping the Surface of Pluto with the Hubble Space Telescope; Experimental Study on Fischer-Tropsch Catalysis in the Circum-Saturnian Subnebula; New High-Pressure Phases of Ammonia Dihydrate; Gas Hydrate Stability at Low Temperatures and High Pressures with Applications to Mars and Europa; Laboratory UV Photolysis of Planetary Ice Analogs Containing H2O + CO2 (1:1); The OH Stretch Infrared Band of Water Ice and Its Temperature and Radiation Dependence; Band Position Variations in Reflectance Spectra of the Jovian Satellite Ganymede; Comparison of Porosity and Radar Models for Europa s Near Surface; Combined Effects of Diurnal and Nonsynchronous Surface Stresses on Europa; Europa s Northern Trailing Hemisphere: Lineament Stratigraphic Framework; Europa at the Highest Resolution: Implications for Surface Processes and Landing Sites; Comparison of Methods to Determine Furrow System Centers on Ganymede and Callisto; Resurfacing of Ganymede by Liquid-Water Volcanism; Layered Ejecta Craters on Ganymede: Comparisons with Martian Analogs; Evaluation of the Possible Presence of CO2-Clathrates in Europa s Icy Shell or Seafloor; Geosciences at Jupiter s Icy Moons: The Midas Touch; Planetary Remote Sensing Science Enabled by MIDAS (Multiple Instrument Distributed Aperture Sensor); and In Situ Surveying of Saturn s Rings.

Source record↗

NASA NTRS1960

An Exploratory Statistical Analysis of a Planet Approach-Phase Guidance Scheme Using Angular Measurements with Significant Error

An exploratory analysis of vehicle guidance during the approach to a target planet is presented. The objective of the guidance maneuver is to guide the vehicle to a specific perigee distance with a high degree of accuracy and minimum corrective velocity expenditure. The guidance maneuver is simulated by considering the random sampling of real measurements with significant error and reducing this information to prescribe appropriate corrective action. The instrumentation system assumed includes optical and/or infrared devices to indicate range and a reference angle in the trajectory plane. Statistical results are obtained by Monte-Carlo techniques and are shown as the expectation of guidance accuracy and velocity-increment requirements. Results are nondimensional and applicable to any planet within limits of two-body assumptions. The problem of determining how many corrections to make and when to make them is a consequence of the conflicting requirement of accurate trajectory determination and propulsion. Optimum values were found for a vehicle approaching a planet along a parabolic trajectory with an initial perigee distance of 5 radii and a target perigee of 1.02 radii. In this example measurement errors were less than i minute of arc. Results indicate that four corrections applied in the vicinity of 50, 16, 15, and 1.5 radii, respectively, yield minimum velocity-increment requirements. Thrust devices capable of producing a large variation of velocity-increment size are required. For a vehicle approaching the earth, miss distances within 32 miles are obtained with 90-percent probability. Total velocity increments used in guidance are less than 3300 feet per second with 90-percent probability. It is noted that the above representative results are valid only for the particular guidance scheme hypothesized in this analysis. A parametric study is presented which indicates the effects of measurement error size, initial perigee, and initial energy on the guidance requirements. Measurement error size significantly affects both guidance accuracy and velocity-increment expenditure. The initial trajectory, as given by its perigee and energy, affects the velocity-increment expenditure but not final guidance accuracy.

Friedlander, Alan L.↗

NASA NTRS1998

The Postshock Chemical Lifetimes of Outflow Tracers and a Possible New Mechanism to Produce Water Ice Mantles

We have used a coupled time-dependent chemical and dynamical model to investigate the lifetime of the chemical legacy in the wake of C-type shocks. We concentrate this study on the chemistry of H2O and O2, two molecules which are predicted to have abundances that are significantly affected in shock-heated gas. Two models are presented: (1) a three-stage model of preshock, shocked, and postshock gas; and (2) a Monte Carlo cloud simulation where we explore the effects of stochastic shock activity on molecular gas over a cloud lifetime. For both models we separately examine the pure gas-phase chemistry as well as the chemistry including the interactions of molecules with grain surfaces. In agreement with previous studies, we find that shock velocities in excess of 10 km/s are required to convert all of the oxygen not locked in CO into H2O before the gas has an opportunity to cool. For pure gas phase models the lifetime of the high water abundances, or "H2O legacy," in the postshock gas is approximately (4-7) x 10(exp 5) yr, independent of the gas density. A density dependence for the lifetime of H2O is found in gas-grain models as the water molecules deplete onto grains at the depletion timescale. Through the Monte Carlo cloud simulation we demonstrate that the time-average abundance of H2O, the weighted average of the amount of time gas spends in preshock, shock, and postshock stages, is a sensitive function of the frequency of shocks. Thus we predict that the abundance of H2O, and to a lesser extent O2, can be used to trace the history of shock activity in molecular gas. We use previous large-scale surveys of molecular outflows to constrain the frequency of 10 km/s shocks in regions with varying star formation properties and discuss the observations required to test these results. We discuss the postshock lifetimes for other possible outflow tracers (e.g., SiO and CH3OH) and show that the differences between the lifetimes for various tracers can produce potentially observable chemical variations between younger and older outflows. For gas-grain models we find that the abundance of water-ice on grain surfaces can be quite large and is comparable to that observed in molecular clouds. This offers a possible alternative method to create water mantles without resorting to grain surface chemistry: gas heating and chemical modification due to a C-type shock and subsequent depletion of the gas-phase species onto grain mantles.

Bergin, Edwin A.↗

NASA NTRSDate not supplied

Computational Study of Oxidative Etch Pitting in FiberForm and the Effect on Its Material Properties

Erosion of carbon due to oxidation does not occur uniformly but through the formation of localized etch pits because of active surface sites. These active sites are formed due to the presence of atomic defects on the carbon surface, and have much higher reactivity compared to average non-defective sites. Thus, these active sites are first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these “active” sites. In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. In order to understand the influence of these etch pits on the material properties of carbon fiber microstructures, we have developed a new capability within direct simulation Monte Carlo (DSMC) to capture the etch pit formation process. This capability is developed within the DSMC code SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) and can model the material removal in presence of active sites leading to the formation of etch pits. The focus of the current work will be to study the effect of etch pits on the material properties of FiberForm, a commonly used base material within many thermal protection system materials (TPS). The microstructure of virgin FiberForm obtained directly from X-ray microtomography experiments is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported within the Porous Microstructure Analysis (PuMA) software and various material properties such as elasticity, thermal conductivity, and permeability are computed. The variation of these properties as a result of the complex evolution of the surface topology due to etch pit formation is studied and analyzed. Furthermore, the effect of pitting is compared to the case of shrinking fibers, which has been the standard for modelling ablation of carbon structures; and significant differences are observed. Thus, such a physically realistic modeling of material removal through the formation of etch pits will be helpful in predicting the degradation of carbon-based TPS more accurately during oxidation; as well as other mechanisms such as spallation, which involves the removal of chunks of material into the flow due to etch pit growth. This will ultimately improve our understanding of the failure modes in these materials due to ablation.

Carbon Ablators↗

NASA NTRSDate not supplied

Computational Study of Oxidative Etch Pitting in FiberForm and the Effect on Its Material Properties

Erosion of carbon due to oxidation does not occur uniformly but through the formation of localized etch pits because of active surface sites. These active sites are formed due to the presence of atomic defects on the carbon surface, and have much higher reactivity compared to average non-defective sites. Thus, these active sites are first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these “active” sites. In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. In order to understand the influence of these etch pits on the material properties of carbon fiber microstructures, we have developed a new capability within direct simulation Monte Carlo (DSMC) to capture the etch pit formation process. This capability is developed within the DSMC code SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) and can model the material removal in presence of active sites leading to the formation of etch pits. The focus of the current work will be to study the effect of etch pits on the material properties of FiberForm, a commonly used base material within many thermal protection system materials (TPS). The microstructure of virgin FiberForm obtained directly from X-ray microtomography experiments is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported within the Porous Microstructure Analysis (PuMA) software and various material properties such as elasticity, thermal conductivity, and permeability are computed. The variation of these properties as a result of the complex evolution of the surface topology due to etch pit formation is studied and analyzed. Furthermore, the effect of pitting is compared to the case of shrinking fibers, which has been the standard for modelling ablation of carbon structures; and significant differences are observed. Thus, such a physically realistic modeling of material removal through the formation of etch pits will be helpful in predicting the degradation of carbon-based TPS more accurately during oxidation; as well as other mechanisms such as spallation, which involves the removal of chunks of material into the flow due to etch pit growth. This will ultimately improve our understanding of the failure modes in these materials due to ablation.

Carbon Ablators↗

NASA NTRSDate not supplied

Simulation of Oxidative Etch Pit Formation and Growth on FiberForm

Erosion of carbon surfaces due to oxidation does not occur uniformly but through the formation of localized etch pits because of active surface sites. These active sites are formed due to the presence of atomic defects on the carbon surface and have much higher reactivity compared to average non-defective sites. Thus, these active sites are the first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective, and increasing their reactivity, thus leading to localized carbon removal around these “active” sites. In this manner, these highly reactive defects serve as nucleation sites for the formation and growth of etch pits, with detrimental effects on the structural integrity. In order to understand the influence of these etch pits on the material properties of carbon fiber microstructures, we have developed a new capability within direct simulation Monte Carlo (DSMC) to capture the etch pit formation process. This capability is developed within the DSMC code SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) and can model the material removal in the presence of active sites leading to the formation of etch pits. The focus of the current work will be to study the effect of etch pits on the material properties of FiberForm, the precursor substrate of the PICA Thermal Protection System (TPS) material. The microstructure of virgin FiberForm obtained directly from X-ray microtomography scans is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported into the Porous Microstructure Analysis (PuMA) software and various material properties such as thermal conductivity, elasticity, and permeability are computed. The variation of these properties because of the complex evolution of the surface topology due to the formation of etch pits is studied and analyzed. Furthermore, the effect of pitting is compared to the case of uniform radial shrinking of fibers, which has been the standard for modelling ablation of carbon structures, and significant differences are observed. Thus, a physically realistic model of material removal through the formation of etch pits will be helpful in predicting the degradation of carbon-based TPS more accurately during oxidation; as well as other mechanisms such as spallation, which involves the removal of chunks of material into the flow due to the growth of etch pits. This will ultimately improve our understanding of the failure modes in these materials due to ablation.

DSMC↗

NASA NTRSDate not supplied

Simulation of Oxidative Etch Pit Formation and Growth on FiberForm

Erosion of carbon surfaces due to oxidation does not occur uniformly but through the formation of localized etch pits because of active surface sites. These active sites are formed due to the presence of atomic defects on the carbon surface and have much higher reactivity compared to average non-defective sites. Thus, these active sites are the first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective, and increasing their reactivity, thus leading to localized carbon removal around these “active” sites. In this manner, these highly reactive defects serve as nucleation sites for the formation and growth of etch pits, with detrimental effects on the structural integrity. In order to understand the influence of these etch pits on the material properties of carbon fiber microstructures, we have developed a new capability within direct simulation Monte Carlo (DSMC) to capture the etch pit formation process. This capability is developed within the DSMC code SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) and can model the material removal in the presence of active sites leading to the formation of etch pits. The focus of the current work will be to study the effect of etch pits on the material properties of FiberForm, the precursor substrate of the PICA Thermal Protection System (TPS) material. The microstructure of virgin FiberForm obtained directly from X-ray microtomography scans is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported into the Porous Microstructure Analysis (PuMA) software and various material properties such as thermal conductivity, elasticity, and permeability are computed. The variation of these properties because of the complex evolution of the surface topology due to the formation of etch pits is studied and analyzed. Furthermore, the effect of pitting is compared to the case of uniform radial shrinking of fibers, which has been the standard for modelling ablation of carbon structures, and significant differences are observed. Thus, a physically realistic model of material removal through the formation of etch pits will be helpful in predicting the degradation of carbon-based TPS more accurately during oxidation; as well as other mechanisms such as spallation, which involves the removal of chunks of material into the flow due to the growth of etch pits. This will ultimately improve our understanding of the failure modes in these materials due to ablation.

DSMC↗

NASA NTRS1994

The time-delay spectrum of GX 5-1 in its horizontal branch

Using a cross-spectral technique we investigate time delays between intensity variations of GX 5-1 in 10 X-ray spectral channels. The data were taken during a 1989 Ginga observation during which the source was in its horizontal-branch spectral state. We develope a new method to measure 'time-delay spectra' in fixed Fourier frequency ranges and use it to determine the energy and intensity dependence of time delays in the low-frequency noise (nu less than 2 Hz), the horizontal branch quasi-periodic oscillations (QPO), and the QPO second harmonic. These are the first time-delay spectra of a Z-source in its horizontal branch, and the first detection of time delays in the second harmonic. We consider two mechanisms for the production of the time lags: Comptonization and evolving shots. We perform Monte Carlo simulations of Compton scattering in a homogeneous, isotropic, central corona and show that it qualitatively explain the observed energy and time-delay spectra, but that it cannot explain the differences in the QPO first and second harmomnic time-delay spectra, nor the observed dependence of the QPO fractional rms variability upon energy. We consider implications of our results for millisecond pulsar searches in low-mass X-ray binaries.

Vaughan, B.↗

NASA NTRS2013

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant↗

NASA NTRS2004

Remote Sensing of Icy Galilean Moon Surface and Atmospheric Composition Using Low Energy (1 eV-4 keV) Neutral Atom Imaging

We describe a low energy neutral atom imager suitable for composition measurements Europa and other icy Galilean moons in the Jovian magnetosphere. This instrument employs conversion surface technology and is sensitive to either neutrals converted to negative ions, neutrals converted to positive ions and the positive ions themselves depending on the power supply. On a mission such as the Jupiter Icy Moons Orbiter (JIMO), two back-to-back sensors would be flown with separate power supplies fitted to the neutral atom and iodneutral atom sides. This will allow both remote imaging of 1 eV < E < 4 keV neutrals from icy moon surfaces and atmospheres, and in situ measurements of ions at similar energies in the moon ionospheres and Jovian magnetospheric plasma. The instrument provides composition measurements of the neutrals and ions that enter the spectrometer with a mass resolution dependent on the time-of-flight subsystem and capable of resolving molecules. The lower energy neutrals, up to tens of eV, arise from atoms and molecules sputtered off the moon surfaces and out of the moon atmospheres by impacts of more energetic (keV to MeV) ions from the magnetosphere. Direct Simulation Monte Carlo (DSMC) models are used to convert measured neutral abundances to compositional distributions of primary and trace species in the sputtered surfaces and atmospheres. The escaping neutrals can also be detected as ions after photo- or plasma-ionization and pickup. Higher energy, keV neutrals come from charge exchange of magnetospheric ions in the moon atmospheres and provide information on atmospheric structure. At the jovicentric orbits of the icy moons the presence of toroidal gas clouds, as detected at Europa's orbit, provide M e r opportunities to analyze both the composition of neutrals and ions originating from the moon surfaces, and the characteristics of magnetospheric ions interacting with neutral cloud material. Charge exchange of low energy ions near the moons, and directional distributions of the resultant neutrals, allow indirect global mapping of magnetic field structures around the moons. Temporal variation of the magnetic structures can be linked to induced magnetic fields associated with subsurface oceans.

Collier, M. R.↗

NASA NTRS2016

Influence of Population Variation of Physiological Parameters in Computational Models of Space Physiology

The earliest manifestations of Visual Impairment and Intracranial Pressure (VIIP) syndrome become evident after months of spaceflight and include a variety of ophthalmic changes, including posterior globe flattening and distension of the optic nerve sheath. Prevailing evidence links the occurrence of VIIP to the cephalic fluid shift induced by microgravity and the subsequent pressure changes around the optic nerve and eye. Deducing the etiology of VIIP is challenging due to the wide range of physiological parameters that may be influenced by spaceflight and are required to address a realistic spectrum of physiological responses. Here, we report on the application of an efficient approach to interrogating physiological parameter space through computational modeling. Specifically, we assess the influence of uncertainty in input parameters for two models of VIIP syndrome: a lumped-parameter model (LPM) of the cardiovascular and central nervous systems, and a finite-element model (FEM) of the posterior eye, optic nerve head (ONH) and optic nerve sheath. Methods: To investigate the parameter space in each model, we employed Latin hypercube sampling partial rank correlation coefficient (LHSPRCC) strategies. LHS techniques outperform Monte Carlo approaches by enforcing efficient sampling across the entire range of all parameters. The PRCC method estimates the sensitivity of model outputs to these parameters while adjusting for the linear effects of all other inputs. The LPM analysis addressed uncertainties in 42 physiological parameters, such as initial compartmental volume and nominal compartment percentage of total cardiac output in the supine state, while the FEM evaluated the effects on biomechanical strain from uncertainties in 23 material and pressure parameters for the ocular anatomy. Results and Conclusion: The LPM analysis identified several key factors including high sensitivity to the initial fluid distribution. The FEM study found that intraocular pressure and intracranial pressure had dominant impact on the peak strains in the ONH and retro-laminar optic nerve, respectively; optic nerve and lamina cribrosa stiffness were also important. This investigation illustrates the ability of LHSPRCC to identify the most influential physiological parameters, which must therefore be well-characterized to produce the most accurate numerical results.

biomechanics↗

NASA NTRS2016

Integrating Chlorophyll fapar and Nadir Photochemical Reflectance Index from EO-1/Hyperion to Predict Cornfield Daily Gross Primary Production

The concept of light use efficiency (Epsilon) and the concept of fraction of photosynthetically active ration (PAR) absorbed for vegetation photosynthesis (PSN), i.e., fAPAR (sub PSN), have been widely utilized to estimate vegetation gross primary productivity (GPP). It has been demonstrated that the photochemical reflectance index (PRI) is empirically related to ε. An experimental US Department of Agriculture (USDA) cornfield in Maryland was selected as our study field. We explored the potential of integrating fAPAR(sub chl) (defined as the fraction of PAR absorbed by chlorophyll) and nadir PRI (PRI(sub nadir)) to predict cornfield daily GPP. We acquired nadir or near-nadir EO-1/Hyperion satellite images that covered the cornfield and took nadir in-situ field spectral measurements. Those data were used to derive the PRI(sub nadir) and fAPAR (sub chl). The fAPAR (sub chl) is retrieved with the advanced radiative transfer model PROSAIL2 and the Metropolis approach, a type of Markov Chain Monte Carlo (MCMC) estimation procedure. We define chlorophyll light use efficiency Epsilon (sub chl) as the ratio of vegetation GPP as measured by eddy covariance techniques to PAR absorbed by chlorophyll (Epsilon(sub chl) = GPP/APAR (sub chl). Daily Epsilon (sub chl) retrieved with the EO-1 Hyperion images was regressed with a linear equation of PRI (sub nadir) Epsilon (sub chl) = Alpha × PRI (sub nadir) + Beta). The satellite Epsilon(sub chl- PRI (sub nadir) linear relationship for the cornfield was implemented to develop an integrated daily GPP model [GPP = (Alpha × PRI(sub nadir) + Beta) × fAPAR (sub chl) × PAR], which was evaluated with fAPAR (sub chl) and PRI (sub nadir) retrieved from field measurements. Daily GPP estimated with this fAPAR (sub chl-) PRI (nadir) integration model was strongly correlated with the observed tower in-situ daily GPP (R(sup 2) = 0.93); with a root mean square error (RMSE) of 1.71 g C mol−(sup -1) PPFD and coefficient of variation (CV) of 16.57%. Both seasonal Epsilon (sub chl) and PRI (sub nadir) were strongly correlated with fAPAR (sub chl ) retrieved from field measurements, which indicates that chlorophyll content strongly affects seasonal epsilon (sub chl) and PRI (sub nadir). We demonstrate the potential capacity to monitor GPP with space-based visible through shortwave infrared (VSWIR) imaging spectrometers such as NASA's soon to be decommissioned EO- 1/Hyperion and the future Hyperspectral Infrared Imager (HyspIRI).

HyspIRI↗

NASA NTRS2005

DSMC Simulations of Ballute Aerothermodynamics Under Hypersonic Rarefied Conditions

This paper presents computational results obtained with the direct simulation Monte Carlo (DSMC) method for towed ballute applications. A ballute is an inflatable drag device that can be used to create a large amount of drag at high altitudes. Consequently, ballutes provide a potential technology for achieving aerocapture when the primary spacecraft velocity reduction (ΔV) is achieved at much higher altitudes than with the conventional rigid aeroshell. Since the ΔV is achieved at relatively high altitudes, rarefaction can be significant and is the motivation for the current study with the DSMC method. Computed surface and flow-field results are presented for a toroidal ballute, isolated tethers when exposed to free-stream flow conditions, and the flow interactions resulting from a toroidal ballute when towed by a six meter diameter Mars Pathfinder shaped (without tethers) spacecraft. All results presented are for Earth entry at velocities of 14 to 7 km/s (primary focus is at 8.55 km/s, same as some previous Titan aerocapture studies) and altitudes of 200 to 100 km. Variations of drag and heating coefficients as a function of rarefaction are presented. A description of the flow structure is provided and also an explanation of how it is affected by shock interactions produced solely by the ballute and those resulting from the two body combination of towed ballute and spacecraft is also given.

James N Moss↗

NASA NTRSDate not supplied

In-Depth Modeling and Simulation Analysis of Artemis Missions Using the Impact Probabilistic Risk Assessment Tool

BACKGROUND The Artemis campaign is a Moon exploration program with a series of six planned missions, five of which will be crewed. These five crewed missions will contain a single mission segment (space flight), or multiple mission segments involving space flight (Orion), lunar landing (LTV) and/or space habitat (Gateway). Each crewed segment faces the risk of unique medical conditions, necessitating medical sets/kits tailored to those specificities. To support and enable a data-driven and evidence-based decision-making process through out a mission’s life cycle, a software tool called IMPACT was developed. Using probabilistic risk assessment (PRA) methodologies, IMPACT (Informing Mission Planning via Analysis of Complex Tradespaces) is a novel tool built for analyzing the possibility of encountering complex medical risks during space flight, and for identifying the medical resources and capabilities needed to treat those potential at-risk medical conditions. IMPACT achieves this by performing hundreds of thousands of Monte Carlo simulations of missions to build aggregate pictures of medical risk. During an extended simulation modeling phase, IMPACT generated analytical results for medical risks, and the medical resources and capabilities to address those risks, for every segment of every crewed Artemis mission. This presentation will highlight the reliability, consistency and validity of IMPACT’s computational modeling techniques and will showcase the library of analytical outcomes generated for the Artemis missions. OVERVIEW During the early stages of IMPACT’s design, architecture and technical requirements collection, “scenarios” (use cases) - achievement goals required for acceptance testing, were identified by stakeholders. IMPACT successfully completed the scenario testing requirements and undertook an extensive operational run phase utilizing a wide range of input combinations with a goal of delivering a cohesive, trustworthy, reliable, vast, and diverse body of evidence. The intent of these modeling runs was to validate consistency in output, ensure solidity of executable operations and to streamline processes by identifying areas requiring efficiency improvements. Using the many missions of Artemis, IMPACT ran variations of operational runs to assess the output for acceptable, as well as unusual characteristics. This rigorous long-term “shakedown” analysis was implemented to help build a collective body of evidence to aid in securing a high level of confidence, reliability, and validity in the output, whether from the applicational components of IMPACT, or the entirety of the operational process. ANTICIPATED ANALYSIS AND CONCLUSION This presentation will discuss the various categories of input criteria; the comparisons in the application of these input criteria to various Artemis missions; the preparation and collection of the body of evidence, and reliability of the computational modeling techniques. This paper serves as an initial analytical overview of IMPACT’s probabilistic risk assessment (PRA) medical risk outputs covering Artemis missions and is not intended to be deemed the official medical response for the Artemis campaign.

Crew Composition↗

NASA NTRS1998

A Monte Carlo Uncertainty Analysis of Ozone Trend Predictions in a Two Dimensional Model

We use Monte Carlo analysis to estimate the uncertainty in predictions of total O3 trends between 1979 and 1995 made by the Goddard Space Flight Center (GSFC) two-dimensional (2D) model of stratospheric photochemistry and dynamics. The uncertainty is caused by gas-phase chemical reaction rates, photolysis coefficients, and heterogeneous reaction parameters which are model inputs. The uncertainty represents a lower bound to the total model uncertainty assuming the input parameter uncertainties are characterized correctly. Each of the Monte Carlo runs was initialized in 1970 and integrated for 26 model years through the end of 1995. This was repeated 419 times using input parameter sets generated by Latin Hypercube Sampling. The standard deviation (a) of the Monte Carlo ensemble of total 03 trend predictions is used to quantify the model uncertainty. The 34% difference between the model trend in globally and annually averaged total O3 using nominal inputs and atmospheric trends calculated from Nimbus 7 and Meteor 3 total ozone mapping spectrometer (TOMS) version 7 data is less than the 46% calculated 1 (sigma), model uncertainty, so there is no significant difference between the modeled and observed trends. In the northern hemisphere midlatitude spring the modeled and observed total 03 trends differ by more than 1(sigma) but less than 2(sigma), which we refer to as marginal significance. We perform a multiple linear regression analysis of the runs which suggests that only a few of the model reactions contribute significantly to the variance in the model predictions. The lack of significance in these comparisons suggests that they are of questionable use as guides for continuing model development. Large model/measurement differences which are many multiples of the input parameter uncertainty are seen in the meridional gradients of the trend and the peak-to-peak variations in the trends over an annual cycle. These discrepancies unambiguously indicate model formulation problems and provide a measure of model performance which can be used in attempts to improve such models.

Considine, D. B.↗

NASA NTRS1996

A Burst Mode, Ultrahigh Temperature UF4 Vapor Core Reactor Rankine Cycle Space Power System Concept

Static and dynamic neutronic analyses have been performed on an innovative burst mode (100's of MW output for a few thousand seconds) Ulvahigh Temperature Vapor Core Reactor (UTVR) space nuclear power system. The NVTR employs multiple, neutronically-coupled fissioning cores and operates on a direct, closed Rankine cycle using a disk Magnetohydrodynamic (MHD) generater for energy conversion. The UTVR includes two types of fissioning core regions: (1) the central Ultrahigh Temperature Vapor Core (UTVC) which contains a vapor mixture of highly enriched UF4 fuel and a metal fluoride working fluid and (2) the UF4 boiler column cores located in the BeO moderator/reflector region. The gaseous nature of the fuel the fact that the fuel is circulating, the multiple coupled fissioning cores, and the use of a two phase fissioning fuel lead to unique static and dynamic neutronic characteristics. Static neutronic analysis was conducted using two-dimensional S sub n, transport theory calculations and three-dimensional Monte Carlo transport theory calculations. Circulating-fuel, coupled-core point reactor kinetics equations were used for analyzing the dynamic behavior of the UTVR. In addition to including reactivity feedback phenomena associated with the individual fissioning cores, the effects of core-to-core neutronic and mass flow coupling between the UTVC and the surrounding boiler cores were also included in the dynamic model The dynamic analysis of the UTVR reveals the existence of some very effectlve inherent reactivity feedback effects that are capable of quickly stabilizing this system, within a few seconds, even when large positive reactivity insertions are imposed. If the UTVC vapor fuel density feedback is suppressed, the UTVR is still inherently stable because of the boiler core liquid-fuel volume feedback; in contrast, suppression of the vapor fuel density feedback in 'conventional" gas core cavity reactors causes them to become inherently unstable. Due to the strength of the negative reactivity feedback in the UTVR, it is found that external reactivity insertions alone are inadequate for bringing about significant power level changes during normal reactor operations. Additional methods of reactivity control such as variations in the gaseous fuel mass flow rate, are needed to achieve the desired power level oontrol.

Dugan, E. T.↗

NASA NTRS2010

Updating Mars-GRAM to Increase the Accuracy of Sensitivity Studies at Large Optical Depths

The Mars Global Reference Atmospheric Model (Mars-GRAM) is an engineering-level atmospheric model widely used for diverse mission applications. Mars-GRAM s perturbation modeling capability is commonly used, in a Monte-Carlo mode, to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL). During the Mars Science Laboratory (MSL) site selection process, it was discovered that Mars-GRAM, when used for sensitivity studies for MapYear=0 and large optical depth values such as tau=3, is less than realistic. From the surface to 80 km altitude, Mars-GRAM is based on the NASA Ames Mars General Circulation Model (MGCM). MGCM results that were used for Mars-GRAM with MapYear set to 0 were from a MGCM run with a fixed value of tau=3 for the entire year at all locations. This has resulted in an imprecise atmospheric density at all altitudes. As a preliminary fix to this pressure-density problem, density factor values were determined for tau=0.3, 1 and 3 that will adjust the input values of MGCM MapYear 0 pressure and density to achieve a better match of Mars-GRAM MapYear 0 with Thermal Emission Spectrometer (TES) observations for MapYears 1 and 2 at comparable dust loading. Currently, these density factors are fixed values for all latitudes and Ls. Results will be presented from work being done to derive better multipliers by including variation with latitude and/or Ls by comparison of MapYear 0 output directly against TES limb data. The addition of these more precise density factors to Mars-GRAM 2005 Release 1.4 will improve the results of the sensitivity studies done for large optical depths.

Justh, Hiliary L.↗