Results for “variational Monte Carlo” — Page 27

NASA NTRS2002

System Level Uncertainty Assessment for Collaborative RLV Design

A collaborative design process utilizing Probabilistic Data Assessment (PDA) is showcased. Given the limitation of financial resources by both the government and industry, strategic decision makers need more than just traditional point designs, they need to be aware of the likelihood of these future designs to meet their objectives. This uncertainty, an ever-present character in the design process, can be embraced through a probabilistic design environment. A conceptual design process is presented that encapsulates the major engineering disciplines for a Third Generation Reusable Launch Vehicle (RLV). Toolsets consist of aerospace industry standard tools in disciplines such as trajectory, propulsion, mass properties, cost, operations, safety, and economics. Variations of the design process are presented that use different fidelities of tools. The disciplinary engineering models are used in a collaborative engineering framework utilizing Phoenix Integration's ModelCenter and AnalysisServer environment. These tools allow the designer to join disparate models and simulations together in a unified environment wherein each discipline can interact with any other discipline. The design process also uses probabilistic methods to generate the system level output metrics of interest for a RLV conceptual design. The specific system being examined is the Advanced Concept Rocket Engine 92 (ACRE-92) RLV. Previous experience and knowledge (in terms of input uncertainty distributions from experts and modeling and simulation codes) can be coupled with Monte Carlo processes to best predict the chances of program success.

Charania, A. C.↗

NASA NTRS1993

Some effects of dust on photometry of high-z galaxies: Confounding the effects of evolution

Photometric observations of very distant galaxies--e.g., color vs. z or magnitude vs. z, have been used over the past decade or so in investigations into the evolution of the stellar component. Numerous studies have predicted significant color variations as a result of evolution, in addition to the shifting of different rest wavelengths into the band of observation. Although there is significant scatter, the data can be fit with relatively straightforward, plausible models for galaxian evolution. In very few cases are the effects of dust extinction included in the models. This is due in a large part to the uncertainty about the distribution and optical properties of the grains, and even whether or not they are present in significant numbers in some types of galaxies such as ellipticals. It is likely that the effects of dust on broadband observations are the greatest uncertainty in studies of very distant galaxies. We use a detailed Monte Carlo radiative transfer model within a spherical geometry for different star/dust distributions to examine the effects of dust on the broadband colors of galaxies as a function of redshift. The model fully accounts for absorption and angular redistribution in scattering. In this summary, we consider only the effects on color vs. redshift for three simple geometries each with the same total dust optical depth. Elsewhere at this conference, Capuano, Thronson, & Witt consider other effects of altering the relative dust/star distribution.

Thronson, H. A., Jr.↗

NASA NTRS1986

Monte Carlo models for the terrestrial exosphere over a solar cycle

Improved Monte Carlo exosphere simulations for the concentration and velocity distribution of hydrogen in the terrestrial exosphere were formulated for minimum, medium, and maximum solar cycle conditions. Both the classical exobase source of hot hydrogen and the plasmaspheric source were included, along with solar radiation pressure and photoionization. It was found that at solar minimum the hydrogen from the charge exchange of hot ions in the plasmasphere exceeds that from the exobase source not only for escape but for the population at geocentric distances greater than 2 earth radii. At about two earth radii, the equivalent temperature is about 50 percent greater than that of the exobase, a situation similar to that of the 'two-temperature' Venusian exosphere. Diurnal variations were calculated for altitudes of up to 30 earth radii; for all levels of solar activity, the inclusion of radiation pressure led to a semidiurnal variation at 3-8 earth radii and a high-altitude night/day asymmetry constituting the 'geotail'.

Tinsley, B. A.↗

DOE OSTI2020

Data Assimilative Optimization of WSA Source Surface and Interface Radii using Particle Filtering

TheWang-Sheeley-Arge (WSA) model estimates solar wind speed and interplanetary magnetic field polarity in the inner heliosphere using global photospheric magnetic field maps. WSA employs the Potential Field Source Surface (PFSS) and Schatten Current Sheet (SCS) models to determine the Sun's global coronal magnetic field. The PFSS and SCS models are connected through two radial parameters, the source surface and interface radii, which specify the overlap region between the inner SCS and outer PFSS models. Though both radii values are adjustable, they have typically been fixed to 2.5 solar radii. Our work highlights how solar wind predictions improve when the radii are allowed to vary over time. Data assimilation using particle filtering (sequential Monte Carlo) is used to infer optimal values over a fixed time window. Solar wind model predictions and satellite observations are compared with a newly developed quality-of-agreement prediction metric. The agreement metric between the model and observations is assumed to correspond to the probability of the two key WSA model parameters, the source surface and interface radii, where the highest metric value implies the optimal radii.We find that the optimal particle filter values of solar radii can perform twice as well as standard values for an exploratory period during Carrington Rotation 1901, with these values also reducing nonphysical kinking effects seen in solar magnetic field lines. Data assimilation choices of input realization and time frame have implications for variation in the solar wind over time. We present this work's theoretical context and practical applications for prediction accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRSDate not supplied

Computational Investigation of Oxidative Etch Pitting in FiberForm and Its Impact on Material Properties

Oxidation-driven carbon erosion does not occur uniformly but rather through the development of localized etch pits at active surface sites. These active sites form due to atomic defects on the carbon surface, making them significantly more reactive than the surrounding, non-defective areas. As a result, these sites are the first to react during ablation, leading to their removal. This process creates new defects in neighboring atoms, increasing their reactivity and causing localized carbon removal around these active sites. In this way, the highly reactive defective areas serve as nucleation points for the formation and growth of etch pits, which can have adverse effects on structural integrity of FiberForm. To better understand how these etch pits impact the material properties of carbon fiber microstructures, we have developed a new capability within the direct simulation Monte Carlo (DSMC) framework to capture the etch pit formation process. This capability, integrated into the DSMC code SPARTA (Stochastic Parallel Rarefied-gas Time-accurate Analyzer), models material removal in the presence of active sites, leading to the formation of etch pits. The current work focuses on studying the effects of these etch pits on the material properties of FiberForm, a widely used base material in thermal protection systems (TPS). The microstructure of virgin FiberForm, obtained via X-ray microtomography, is imported into SPARTA to generate the ablated geometries with etch pits. These modified microstructures are then analyzed using the Porous Microstructure Analysis (PuMA) software to compute various material properties, including elasticity, thermal conductivity, and permeability. We investigate the variation of these properties due to the complex surface topology changes caused by etch pit formation. Additionally, we compare the effects of pitting with the conventional model of shrinking fibers, traditionally used to simulate the ablation of carbon structures. Significant differences emerge between the two approaches. Consequently, this physically realistic model of material removal through etch pit formation offers improved accuracy in predicting the degradation of carbon-based TPS during oxidation. It also provides insights into other mechanisms, such as spallation, where chunks of material are removed into the flow due to etch pit growth. Ultimately, this model enhances our understanding of failure modes in these materials during ablation.

PuMA↗

NASA NTRS1987

Variational approach to probabilistic finite elements

Probabilistic finite element method (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties, and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗

NASA NTRS1991

Variational approach to probabilistic finite elements

Probabilistic finite element methods (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗

NASA NTRS1990

Near infrared radiances observed by the UK C130 multi-channel radiometer during the marine stratocumulus IFO and preliminary comparison with model calculations

A preliminary analysis of some of the narrow band radiance data measured on the U.K. Meteorological Office's C130 aircraft during the marine stratocumulus intensive field observation of First ISCCP Regional Experiment (FIRE), San Diego 29 June to 18 July 1987, is presented. The data are compared with Monte Carlo calculations of the reflectance and transmittance of the cloud based upon the observed droplet size distribution. The main scientific question being addressed is whether there is any evidence of anomalous absorption within the cloud which had been observed in similar measurements (Rozenberg et al., 1974; Twomey and Cocks, 1982; Foot, 1988). The measurements also indicate the potential for remotely sensing cloud properties. The data and method of presentation discussed here clearly separates out clouds in terms of the size of the cloud droplets. All of the daytime C130 FIRE flights have been studied and are consistent with the data presented here. There appears to be no peculiarities that might arise, for example if pollution were to be a significant factor in determining cloud absorption. Variation in the inferred size parameters, r sub e, along runs are also very small.

Foot, J. S.↗

NASA NTRS2008

Stochastic Transport Theory for Investigating the Three-Dimensional Canopy Structure from Space Measurements

Radiation reflected from vegetation canopies exhibits high spatial variation. Satellite-borne sensors measure the mean intensities emanating from heterogeneous vegetated pixels. The theory of radiative transfer in stochastic media provides the most logical linkage between satellite observations and the three-dimensional canopy structure through a closed system of simple equations which contains the mean intensity and higher statistical moments directly as its unknowns. Although this theory has been a highly active research field in recent years, its potential for satellite remote sensing of vegetated surfaces has not been fully realized because of the lack of models of a canopy pair-correlation function that the stochastic radiative transfer equations require. The pair correlation function is defined as the probability of finding simultaneously phytoelements at two points. This paper presents analytical and Monte Carlo generated pair correlation functions. Theoretical and numerical analyses show that the spatial correlation between phytoelements is primarily responsible for the effects of the three-dimensional canopy structure on canopy reflective and absorptive properties. The pair correlation function, therefore, is the most natural and physically meaningful measure of the canopy structure over a wide range of scales. The stochastic radiative transfer equations naturally admit this measure and thus provide a powerful means to investigate the three-dimensional canopy structure from space. Canopy reflectances predicted by the stochastic equations are assessed by comparisons with the PARABOLA measurements from coniferous and broadleaf forest stands in the BOREAS Southern Study Areas. The pair correlation functions are derived from data on tree structural parameters collected during field campaigns conducted at these sites. The simulated canopy reflectances compare well with the PARABOLA data.

Huang, Dong↗

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↗

DOE OSTI2021

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water–solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS ( i.e. polymer in 100% D 2 O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water–polymer interactions ( i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water–polymer interactions due to the increase in the number of forced water–polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water–polymer interactions and their impact on materials design.

36 MATERIALS SCIENCE↗