Results for “balanced clustering” — Page 4

NASA NTRS2018

The Diversity of Refractory Organic Material in Comets

Organic matter exists in comets (most notably in 81P/Wild 2 [Stardust], 67P/Churyomov- Gerasimenko (67P/C-G) [Rosetta], chondritic porous IDPS, and UCAMMs) and in primitive carbonaceous chondrites that likely retain some chemistry that reflects an origin in the prenatal cold molecular cloud (Alexander+2017). Heavy isotopic enrichments, 15N/14N and possibly D/H, signify preserved molecular cloud organics. In the cold outer disk, if grains are lofted above the disk mid-plane then organics likely experience significant UV processing (Ciesla+2012). In remote sensing of comet comae, organics in the dust are considered refractory or semirefractory. Semi-refractory organics have limited comae lifetimes and produce distributed sources of molecules (H2CO and CO). Rosetta's close passes of 67P/C-G's nucleus (10-15 km) reveals a distributed source of glycine, methyl amine and ethylamine (Altwegg+2016). Cometary samples and primitive meteorites have two types of organic matter: (1) acid-insoluble organic matter (IOM), which is a macromolecular polymer with a mixture of aromatic and aliphatic moieties, and (2) labile, soluble organics, which includes the amino acids, such as glycine (Stardust, Elsila+2009). Meteoritic IOM is robust, withstanding experimental temperatures of 1200 K (Dobrica+2011, Cody+2008). Nanoglobules are a type of IOM; they have a distinct physical structure, but often share the same chemistry as the other IOM from the same meteorite. Moderate-sized PAHs (20 C-atoms) are detected in Stardust samples (Clemett+2010). Refractory organic IOM is ubiquitous yet has a great diversity of abundances between cometary samples. IOM is in primitive chondrites, 67P/C-G (Rosetta), 81P/Wild 2 (Stardust), 1P/Halley, 26P/Grigg-Skjellerup, UCAMMs, anhydrous IDPs, and in chondritic porous IDPs (CP IDPs) and larger cluster IDPs (e.g., Fray+2016, Fomenkova+94, Busemann+ 09, Dobrica+2011, Dobrica+2012). 81P's refractory organic matter is of two types (De Gregorio+2011): nanoglobules of highly aromatic refractory organic matter and polyaromatic carbonyl-containing organic matter, which is similar to IOM in primitive meteorites and IDPs. Fray+2017 estimate that 50% of carbon in 67P/C-G is in IOM. 67P/C-G's organics appear to lack the soluble organic matter, aliphatic carbon, amino acids, and PAHs (Fray+2016). Other notable aspects of the diversity in IOM in cometary samples are the ranges of atomic ratios of N/C, O/C, and H/C, and the range of isotopic enrichments of 15N/14N and D/H. Aqueous and thermal processing on asteroids changes the balance of soluble to insoluble organics, and may be important for diversifying the range of OM delivered to Earth.

Wooden, Diane↗

NASA NTRSDate not supplied

Communicating Metrics of Land Surface Temperature Variability Using Multi-sensor Machine Learning

Land surface temperature (LST) is a key climate observable used to detect changes in the Earth’s surface energy budget that influence carbon and water cycles. Land surface temperature exhibits strong diurnal variability, which geostationary satellites can observe at scale thanks to their temporal resolution. Due to anthropogenic climate and land use changes, the surface energy balance has been considerably modified and may be described by changes in diurnal temperature range and extremes. Using high performance computing and datasets from the NASA Earth Exchange, we exploit co-located, co-temporal observations from low-earth orbit (LEO) and geostationary (GEO) sensors to develop a deep learning-based method for LEO-to-GEO algorithm emulation. Our model is trained to predict MODIS Terra LST from GOES-16 thermal bands and achieves validation error <2K. Application of the model to unseen times of day (observed by MODIS Aqua) and a new GEO sensor (Himawari-8) observing an unseen spatial domain, demonstrate the generalization of the deep learning model across space, time and spectra. Further, time series clustering approaches are examined with the objective of identifying key indicators of change in diurnal cycling and extremes on a continental scale. Communicating LST variability observed by geostationary satellites can have impacts in multiple disciplines, from understanding of snow, vegetation and soil dynamics, to recognizing trends in heat events relevant to human health.

Kate Duffy↗

NASA NTRS2017

Orion Multi-Purpose Crew Vehicle Solving and Mitigating the Two Main Cluster Pendulum Problem

The Orion Multi-purpose Crew Vehicle (MPCV) Orion spacecraft will return humans from beyond earth's orbit, including Mars and will be required to land 20,000 pounds of mass safely in the ocean. The parachute system nominally lands under 3 main parachutes, but the system is designed to be fault tolerant and land under 2 main parachutes. During several of the parachute development tests, it was observed that a pendulum, or swinging, motion could develop while the Crew Module (CM) was descending under two parachutes. This pendulum effect had not been previously predicted by modeling. Landing impact analysis showed that the landing loads would double in some places across the spacecraft. The CM structural design limits would be exceeded upon landing if this pendulum motion were to occur. The Orion descent and landing team was faced with potentially millions of dollars in structural modifications and a severe mass increase. A multidisciplinary team was formed to determine root cause, model the pendulum motion, study alternate canopy planforms and assess alternate operational vehicle controls & operations providing mitigation options resulting in a reliability level deemed safe for human spaceflight. The problem and solution is a balance of risk to a known solution versus a chance to improve the landing performance for the next human-rated spacecraft.

Ali, Yasmin↗

NASA NTRS2018

Improved Hydrology over Peatlands in a Global Land Modeling System

Peatlands of the Northern Hemisphere represent an important carbon pool that mainly accumulated since the last ice age under permanently wet conditions in specific geological and climatic settings. The carbon balance of peatlands is closely coupled to water table dynamics. Consequently, the future carbon balance over peatlands is strongly dependent on how hydrology in peatlands will react to changing boundary conditions, e.g. due to climate change or regional water level drawdown of connected aquifers or streams. Global land surface modeling over organic-rich regions can provide valuable global-scale insights on where and how peatlands are in transition due to changing boundary conditions. However, the current global land surface models are not able to reproduce typical hydrological dynamics in peatlands well. We implemented specific structural and parametric changes to account for key hydrological characteristics of peatlands into NASA's GEOS-5 Catchment Land Surface Model (CLSM, Koster et al. 2000). The main modifications pertain to the modeling of partial inundation, and the definition of peatland-specific runoff and evapotranspiration schemes. We ran a set of simulations on a high performance cluster using different CLSM configurations and validated the results with a newly compiled global in-situ dataset of water table depths in peatlands. The results demonstrate that an update of soil hydraulic properties for peat soils alone does not improve the performance of CLSM over peatlands. However, structural model changes for peatlands are able to improve the skill metrics for water table depth. The validation results for the water table depth indicate a reduction of the bias from 2.5 to 0.2 m, and an improvement of the temporal correlation coefficient from 0.5 to 0.65, and from 0.4 to 0.55 for the anomalies. Our validation data set includes both bogs (rain-fed) and fens (ground and/or surface water influence) and reveals that the metrics improved less for fens. In addition, a comparison of evapotranspiration and soil moisture estimates over peatlands will be presented, albeit only with limited ground-based validation data. We will discuss strengths and weaknesses of the new model by focusing on time series of specific validation sites.

Bechtold, M.↗

NASA NTRSDate not supplied

Effect of Engine Thrust and Isp Tradeoffs and Alternate Propellants on ΔV Budget and Architecture Mass for 1st Generation Nuclear Thermal Propulsion Flight Test Systems

Following the first Nuclear Thermal Propulsion (NTP) system test, also known as DRACO, the next NTP system to be developed would be the 1st Generation NTP system. An analysis was conducted to determine the performance of different vehicle configurations utilizing hydrogen (H-NTP), ammonia (A-NTP), and methane (M-NTP) as propellants launched onboard commercial launch vehicles. This was enabled by a quasi-steady-state power balance engine model and vehicle component physics that sized the vehicle system using Master Equipment List (MEL) parameters. The analysis considered configuration cases outlined by a matrix of different mission classes and vehicle configurations that covered the design space of the 1st Generation NTP system to explore various propellant options, engine architectures, and mission scopes. Four mission classes were evaluated which included single burn missions performing maneuvers having a ΔV of 1 km/s and 2.5 km/s followed by 2-burn and 4-burn missions that aimed to exhaust the launch vehicles volume or mass limitations. In all cases, the NTP systems with the lowest thrust class had the longest burn time of which A-NTP and M-NTP systems provided the longest and shortest burn times depending on the launch vehicle used while H-NTP systems tended to cluster together in the middle. Longer burn times could be useful as a testing platform to increase the time for data accumulation. Across the multi-burn cases, H-NTP systems were found to be volume limited while A-NTP and M-NTP systems were mass limited. Both single burn missions showed that A-NTP configurations provided the lowest dry mass given that ammonia had the highest density with comparable performance to M-NTP systems and no requirement for cryocoolers. The results showed that beyond the propellant selection type, the launch vehicle selection, which included Starship, New Glenn, Vulcan, Falcon 9 (recoverable), and Falcon Heavy (recoverable), was a primary driving factor in the test vehicle’s capabilities. Trends were determined based on a set of dimensionless parameters that included the mass ratio of inert mass to initial wetted mass, ratio of specific impulse to burn time, and a dimensionless engine parameter (ratio of impulse to engine momentum). These relationships found the “knees-in-the-curves” that could be a significant point of reference for the designer as they indicate a change in the trend which is located at a specific impulse to burn time ratio of 1, a mass ratio of 0.6, and an engine performance parameter of 4. This study did not attempt to make a recommendation rather provide the tools for the reader to select their own configuration based on their needs.

Propellant↗

NASA NTRSDate not supplied

Effect of Engine Thrust and Isp Tradeoffs and Alternate Propellants on ΔV Budget and Architecture Mass for 1st Generation Nuclear Thermal Propulsion Flight Test Systems

Following the first Nuclear Thermal Propulsion (NTP) system test, also known as DRACO, the next NTP system to be developed would be the 1st Generation NTP system. An analysis was conducted to determine the performance of different vehicle configurations utilizing hydrogen (H-NTP), ammonia (A-NTP), and methane (M-NTP) as propellants launched onboard commercial launch vehicles. This was enabled by a quasi-steady-state power balance engine model and vehicle component physics that sized the vehicle system using Master Equipment List (MEL) parameters. The analysis considered configuration cases outlined by a matrix of different mission classes and vehicle configurations that covered the design space of the 1st Generation NTP system to explore various propellant options, engine architectures, and mission scopes. Four mission classes were evaluated which included single burn missions performing maneuvers having a ΔV of 1 km/s and 2.5 km/s followed by 2-burn and 4-burn missions that aimed to exhaust the launch vehicles volume or mass limitations. In all cases, the NTP systems with the lowest thrust class had the longest burn time of which A-NTP and M-NTP systems provided the longest and shortest burn times depending on the launch vehicle used while H-NTP systems tended to cluster together in the middle. Longer burn times could be useful as a testing platform to increase the time for data accumulation. Across the multi-burn cases, H-NTP systems were found to be volume limited while A-NTP and M-NTP systems were mass limited. Both single burn missions showed that A-NTP configurations provided the lowest dry mass given that ammonia had the highest density with comparable performance to M-NTP systems and no requirement for cryocoolers. The results showed that beyond the propellant selection type, the launch vehicle selection, which included Starship, New Glenn, Vulcan, Falcon 9 (recoverable), and Falcon Heavy (recoverable), was a primary driving factor in the test vehicle’s capabilities. Trends were determined based on a set of dimensionless parameters that included the mass ratio of inert mass to initial wetted mass, ratio of specific impulse to burn time, and a dimensionless engine parameter (ratio of impulse to engine momentum). These relationships found the “knees-in-the-curves” that could be a significant point of reference for the designer as they indicate a change in the trend which is located at a specific impulse to burn time ratio of 1, a mass ratio of 0.6, and an engine performance parameter of 4. This study did not attempt to make a recommendation rather provide the tools for the reader to select their own configuration based on their needs.

Propellant↗

NASA NTRSDate not supplied

Effect of Engine Thrust and Isp Tradeoffs and Alternate Propellants on ΔV Budget and Architecture Mass for 1st Generation Nuclear Thermal Propulsion Flight Test Systems

Following the first Nuclear Thermal Propulsion (NTP) system test, also known as DRACO, the next NTP system to be developed would be the 1st Generation NTP system. An analysis was conducted to determine the performance of different vehicle configurations utilizing hydrogen (H-NTP), ammonia (A-NTP), and methane (M-NTP) as propellants launched onboard commercial launch vehicles. This was enabled by a quasi-steady-state power balance engine model and vehicle component physics that sized the vehicle system using Master Equipment List (MEL) parameters. The analysis considered configuration cases outlined by a matrix of different mission classes and vehicle configurations that covered the design space of the 1st Generation NTP system to explore various propellant options, engine architectures, and mission scopes. Four mission classes were evaluated which included single burn missions performing maneuvers having a ΔV of 1 km/s and 2.5 km/s followed by 2-burn and 4-burn missions that aimed to exhaust the launch vehicles volume or mass limitations. In all cases, the NTP systems with the lowest thrust class had the longest burn time of which A-NTP and M-NTP systems provided the longest and shortest burn times depending on the launch vehicle used while H-NTP systems tended to cluster together in the middle. Longer burn times could be useful as a testing platform to increase the time for data accumulation. Across the multi-burn cases, H-NTP systems were found to be volume limited while A-NTP and M-NTP systems were mass limited. Both single burn missions showed that A-NTP configurations provided the lowest dry mass given that ammonia had the highest density with comparable performance to M-NTP systems and no requirement for cryocoolers. The results showed that beyond the propellant selection type, the launch vehicle selection, which included Starship, New Glenn, Vulcan, Falcon 9 (recoverable), and Falcon Heavy (recoverable), was a primary driving factor in the test vehicle’s capabilities. Trends were determined based on a set of dimensionless parameters that included the mass ratio of inert mass to initial wetted mass, ratio of specific impulse to burn time, and a dimensionless engine parameter (ratio of impulse to engine momentum). These relationships found the “knees-in-the-curves” that could be a significant point of reference for the designer as they indicate a change in the trend which is located at a specific impulse to burn time ratio of 1, a mass ratio of 0.6, and an engine performance parameter of 4. This study did not attempt to make a recommendation rather provide the tools for the reader to select their own configuration based on their needs.

Propellant↗

NASA NTRS2009

Molecular Cooling as a Probe of Star Formation: Spitzer Looking Forward to Herschel

We explore here the question of how cloud physics can be more directly probed when one observes the majority of cooling emissions from molecular gas. For this purpose we use results from a recent Spitzer Space Telescope study of the young cluster of embedded objects in NGC1333. For this study we mapped the emission from eight pure H2 rotational lines, from S(0) to S(7). The H2 emission appears to be associated with the warm gas shocked by the multiple outflows present in the region. The H2 lines are found to contribute to 25 - 50% of the total outflow luminosity, and can be used to more directly ascertain the importance of star formation feedback on the natal cloud. From these lines, we determine the outflow mass loss rate and, indirectly, the stellar infall rate, the outflow momentum and the kinetic energy injected into the cloud over the embedded phase. The latter is found to exceed the binding energy of individual cores, suggesting that outflows could be the main mechanism for cores disruption. Given the recent launch of Herschel and the upcoming operational lifetime of SOFIA we discuss how studies of molecular cooling can take a step beyond understanding thermal balance to exploring the origin, receipt, and transfer of energy in atomic and molecular gas in a wide range of physical situations.

Bergin, Edwin A.↗

NASA NTRS2019

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗

NASA NTRS2017

Measuring Shock Stage of ltokawa Regolith Grains by Electron Back-Scattered Diffraction and Synchrotron X-Ray Diffraction

We have been analyzing Itokawa samples in order to definitively establish the degree of shock experienced by the regolith of asteroid Itokawa, and to devise a bridge between shock determinations by standard light optical petrography, crystal structures as determined by electron and X-ray diffraction techniques. We are making measurements of olivine crystal structures and using these to elucidate critical regolith impact processes. We use electron back-scattered diffraction (EBSD) and synchrotron X-ray diffraction (SXRD). We are comparing the Itokawa samples to L and LL chondrite meteorites chosen to span the shock scale experienced by Itokawa, specifically Chainpur (LL3.4, Shock Stage 1), Semarkona (LL3.00, S2), Kilabo (LL6, S3), NWA100 (L6, S4) and Chelyabinsk (LL5, S4). In SXRD we measure the line broadening of olivine reflections as a measure of shock stage. In this presentation we concentrate on the EBSD work. We employed JSC's Supra 55 variable pressure FEG-SEM and Bruker EBSD system. We are not seeking actual strain values, but rather indirect strain-related measurements such as extent of intra-grain lattice rotation, and determining whether shock state "standards" (meteorite samples of accepted shock state, and appropriate small grain size) show strain measurements that may be statistically differentiated, using a sampling of particles (number and size range) typical of asteroid regoliths. Using our system we determined that a column pressure of 9 Pa and no C-coating on the sample was optimal. We varied camera exposure time and gain to optimize mapping performance, concluding that 320x240 pattern pixilation, frame averaging of 3, 15 kV, and low extractor voltage yielded an acceptable balance of hit rate (>90%), speed (11 fps) and map quality using an exposure time of 30 ms (gain 650). We found that there was no strong effect of step size on Grain Orientation Spread (GOS) and Grain Reference Orientation Deviation angle (GROD-a) distribution; there was some effect on grain average Kernel Average Misorientation (KAM) (reduced with smaller step size for the same grain), as expected. We monitored GOS, Maximum Orientation Spread (MOS) and GROD-a differences between whole olivine grains and sub-sampled areas, and found that there were significant differences between the whole grain dataset and subsets, as well as between subsets, likely due to sampling-related "noise". Also, in general (and logically) whole grains exhibit greater degrees of cumulative lattice rotation. Sampling size affects the apparent strain character of the grain, at least as measured by GOS, MOS and GROD-a. There were differences in the distribution frequencies of GOS and MOS between shock stages, and in plots of MOS and GOS vs. grain diameter. These results are generally consistent with those reported this year. However, it is unknown whether the differences between samples of different shock states exceeds the clustering of these values to the extent that shock stage determinations can still be made with confidence. We are investigating this by examination of meteorites with higher shock stage 4 to 5. Our research will improve our understanding of how small, primitive solar system bodies formed and evolved, and improve understanding of the processes that determine the history and future of habitability of environments on other solar system bodies. The results will directly enrich the ongoing asteroid and comet exploration missions by NASA and JAXA, and broaden our understanding of the origin and evolution of small bodies in the early solar system, and elucidate the nature of asteroid and comet regolith.

Zolensky, Michael↗

NASA NTRS2007

Automated Knowledge Discovery From Simulators

A computational method, SimLearn, has been devised to facilitate efficient knowledge discovery from simulators. Simulators are complex computer programs used in science and engineering to model diverse phenomena such as fluid flow, gravitational interactions, coupled mechanical systems, and nuclear, chemical, and biological processes. SimLearn uses active-learning techniques to efficiently address the "landscape characterization problem." In particular, SimLearn tries to determine which regions in "input space" lead to a given output from the simulator, where "input space" refers to an abstraction of all the variables going into the simulator, e.g., initial conditions, parameters, and interaction equations. Landscape characterization can be viewed as an attempt to invert the forward mapping of the simulator and recover the inputs that produce a particular output. Given that a single simulation run can take days or weeks to complete even on a large computing cluster, SimLearn attempts to reduce costs by reducing the number of simulations needed to effect discoveries. Unlike conventional data-mining methods that are applied to static predefined datasets, SimLearn involves an iterative process in which a most informative dataset is constructed dynamically by using the simulator as an oracle. On each iteration, the algorithm models the knowledge it has gained through previous simulation trials and then chooses which simulation trials to run next. Running these trials through the simulator produces new data in the form of input-output pairs. The overall process is embodied in an algorithm that combines support vector machines (SVMs) with active learning. SVMs use learning from examples (the examples are the input-output pairs generated by running the simulator) and a principle called maximum margin to derive predictors that generalize well to new inputs. In SimLearn, the SVM plays the role of modeling the knowledge that has been gained through previous simulation trials. Active learning is used to determine which new input points would be most informative if their output were known. The selected input points are run through the simulator to generate new information that can be used to refine the SVM. The process is then repeated. SimLearn carefully balances exploration (semi-randomly searching around the input space) versus exploitation (using the current state of knowledge to conduct a tightly focused search). During each iteration, SimLearn uses not one, but an ensemble of SVMs. Each SVM in the ensemble is characterized by different hyper-parameters that control various aspects of the learned predictor - for example, whether the predictor is constrained to be very smooth (nearby points in input space lead to similar output predictions) or whether the predictor is allowed to be "bumpy." The various SVMs will have different preferences about which input points they would like to run through the simulator next. SimLearn includes a formal mechanism for balancing the ensemble SVM preferences so that a single choice can be made for the next set of trials.

Burl, Michael↗

NASA NTRS1997

Infrared Line Emission from Molecular Gas Heated by X-Rays and Energetic Electrons

"I propose to carry out a detailed study using infrared observations (and in some cases, optical and ultraviolet observations) of dense interstellar gas exposed to intense fluxes of X-rays and/or energetic electrons. This is undoubtedly the dominant source of line emission for clouds exposed to X-rays from active galactic nuclei, supernova shocks, or embedded X-ray sources (e.g., X-ray binaries), or to high-temperature or relativistic electrons in galaxy clusters, near powerful radio sources, or supernova remnants. Detailed physical and chemical models of such clouds will be used to analyze infrared observations of the Great Annihilator X-ray source in the Galactic Center, cD galaxies in massive cooling flows, and the nuclei of Seyfert galaxies which will be obtained with the Infrared Space Observatory (ISO), UV and optical observations of the Crab Nebula obtained with the Hubble Space Telescope, and ground-based near-infrared observations of Seyfert nuclei. Results from this work will also be of great relevance to observations obtained with the Submillimeter Wave Astronomical Satellite (SWAS). In the first year of funding of this proposal, my chief collaborators (D.J. Hollenbach and A.G.G.M. Tielens, both of NASA Ames Research Center) and I concentrated on completing our models of the physical conditions in, and the resulting line emission from, dense gas irradiated by X-rays. As noted in the original proposal, some important physical processes were not yet thoroughly incorporated into our models at the time of submission. We completed our modeling of the physical conditions and line emission for essentially the entire range of parameter space (five orders of magnitude in X-ray flux to gas density ratio) occupied by typical dense interstellar clouds in which the gas is mostly neutral and X-rays are important for the ionization, chemistry, and thermal balance.

Maloney, Philip R.↗

NASA NTRS2013

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric↗

NASA NTRS2000

A Study of Bubble and Slug Gas-Liquid Flow in a Microgravity Environment

The influence of gravity on the two-phase flow dynamics is obvious.As the gravity level is reduced,there is a new balance between inertial and interfacial forces, altering the behavior of the flow. In bubbly flow,the absence of drift velocity leads to spherical-shaped bubbles with a rectilinear trajectory.Slug flow is a succession of long bubbles and liquid slug carrying a few bubbles. There is no flow reversal in the thin liquid film as the long bubble and liquid slug pass over the film. Although the flow structure seems to be simpler than in normal gravity conditions,the models developed for the prediction of flow behavior in normal gravity and extended to reduced gravity flow are unable to predict the flow behavior correctly.An additional benefit of conducting studies in microgravity flows is that these studies aide the development of understanding for normal gravity flow behavior by removing the effects of buoyancy on the shape of the interface and density driven shear flows between the gas and the liquid phases. The proposal calls to study specifically the following: 1) The dynamics of isolated bubbles in microgravity liquid flows will be analyzed: Both the dynamics of spherical isolated bubbles and their dispersion by turbulence, their interaction with the pipe wall,the behavior of the bubbles in accelerated or decelerated flows,and the dynamics of isolated cylindrical bubbles, their deformation in accelerated/decelerated flows (in converging or diverging channels), and bubble/bubble interaction. Experiments will consist of the use of Particle Image Velocimetry (PIV) and Laser Doppler Velocimeters (LDV) to study single spherical bubble and single and two cylindrical bubble behavior with respect to their influence on the turbulence of the surrounding liquid and on the wall 2) The dynamics of bubbly and slug flow in microgravity will be analyzed especially for the role of the coalescence in the transition from bubbly to slug flow (effect of fluid properties and surfactant), to identify clusters that promote coalescence and transition the void fraction distribution in bubbly and slug flow,to measure the wall friction in bubbly flow. These experiments will consist of multiple bubbles type flows and will utilize hot wire and film anemometers to measure liquid velocity and wall shear stress respectively and double fiber optic probes to measure bubble size and velocity as a function of tube radius and axial location.

McQuillen, J.↗

NASA NTRSDate not supplied

Metallic Powder Core Tubular Wire (PCTW) Development for Additive Manufacturing

NASA Langley Research Center (LaRC) in Hampton, VA, and the Center for Welding, Joining, and Coating Research (CWJCR) in the Department of Metallurgical and Materials Engineering at the Colorado School of Mines (CSM) in Golden, CO, have collaborated for nearly fifteen years to develop metallic powder core tubular wire (PCTW) feedstock for additive manufacturing (AM) using electron beam deposition. The motivation for this development was to offset dealloying observed during deposition of titanium and aluminum alloy solid wire feedstock in the LaRC electron beam freeform fabrication (EBF3) system and to explore novel alloy compositions by tailoring the powder fill. EBF3 is a fusion-based AM process that is performed under vacuum and is consequently prone to vaporization loss of low vapor pressure alloying elements. Research spanning four graduate student research projects succeeded in developing the basic methodology to fabricate PCTW, the mass balance equations used to define the powder fill, microstructure control through particulate inoculation, and the feasibility of producing metal matrix composite materials. Metallic PCTW is comprised of a thin foil sheath wrapped around a powder blend. Fabrication of PCTW at CSM was accomplished using a tubular wire mill system to first bend a flat strip into a u-shaped cross section, deposit a powder fill, create a lap joint to seal the wire, and draw to final diameter. Powder management systems were developed to ensure uniform mixing of powder constituents and identify feed rates for optimum fill ratios. Processing parameters were successfully developed for titanium and aluminum alloy sheaths to account for material work hardening during drawing. Mass balance equations were developed to account for compositional contributions of the sheath and powder fill. PCTW composition can be adjusted through selection of the sheath or powder materials. PCTWs were successfully developed to offset Al loss in deposits of Ti-6-4 and Mg loss in Al 6061, resulting in deposits that were within the composition limits for each material. Ti-6-4 deposits were fabricated using PCTW comprised of commercially pure (CP) titanium sheath and a fill of pre-alloyed Ti-6-4 powder with additions of elemental Al and V powders. Additions of both elements were required to account for the CP Ti from the sheath, and Al was further enhanced to offset documented losses in EBF3 deposits. For Al 6061 deposits, composition was modulated through the choice of sheath material. The PCTW was made using an Al 5052 sheath to supplement Mg content and a fill composition of pre-alloyed 6061 powder with additions of Si and Cu.The slower solidification rate of the EBF3 process compared with laser powder bed fusion processes can lead to coarse-grained microstructures and mechanical property anisotropy. EBF3deposits of Ti-6-4 exhibited large epitaxially grown β grains, which result in degraded mechanical properties in the deposition direction. Microstructure refinement was achieved through a combination of electron beam modulation and the use of a PCTW designed for Ti-6-4 deposits and modified by the addition of powder made from Fe and B. Resulting deposits exhibited reductions in β grain size, α phase colony intercept length and lath width, and an increased volume fraction of β phase. The mechanism of microstructure refinement was related to precipitation of boride particles. Hardness testing confirmed that mechanical property anisotropy was reduced and an overall increase in hardness was achieved.Metal matrix composite (MMC) materials offer strength and stiffness benefits over unreinforced materials. Two methods were investigated using PCTW technology for producing Al 6061 basedMMC materials; the direct addition of SiC reinforcing particles and the in-situ generation of reinforcements through reaction synthesis of precursor powders. Al 6061 based MMCs reinforced with SiC exhibited particle clustering and formation of aluminum carbides during deposition. The application of a Ni coating to the SiC particles prior to incorporation into the powder fill eliminated aluminum carbide formation and aided particle distribution during deposition. Modulation of beam focus and power further reduced aluminum carbide formation and aided particle distribution. In reaction synthesis, precursor additions undergo an exothermic reaction during deposition to form new ceramic products, generally of finer scale than the initial precursor components. The product phases provide reinforcement in the MMC and act as nucleation sites during solidification, resulting in grain refinement, improved mechanical properties, and reduced susceptibility to solidification cracking. Higher tensile strength and stiffness were demonstrated in Al 6061 based MMCs with 2 vol% reinforcement. These combined research efforts demonstrate the potential of PCTW technology for customizing of alloy composition and achieving microstructure control. A manufacturing process for titanium- and aluminum-based PCTWs was developed based on iterative design, theoretical calculations, and empirical data. PCTW technology successfully mitigated Al loss from Ti-6-4 and Mg loss in Al 6061, resulted in grain refinement in Ti-6-4, and demonstrated fabrication of Al based MMCs

Marcia S Domack↗