Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic material flow analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

Dislocation density-based plasticity model from massive discrete dislocation dynamics database

In this work, we present a dislocation density-based strain hardening model for single crystal copper through a systematic coarse-graining analysis of more than 200 discrete dislocation dynamics (DDD) simulations of plastic deformation under uniaxial tension. The proposed constitutive model has two components: a generalized Taylor relation connecting resolved shear stresses to dislocation densities on individual slip systems, and a generalized Kocks-Mecking model for dislocation multiplication. The DDD data strongly suggests a logarithmic dependence of flow stress on the plastic shear strain rate on each slip system, and, equivalently, an exponential dependence of the plastic shear strain rate on the resolved shear stress. Hence the proposed generalized Taylor relation subsumes the Orowan relation for plastic flow. The DDD data also calls for a correction to the Kocks-Mecking model of dislocation multiplication to account for the increase of dislocation density on slip systems with negligible plastic shear strain rate. This is accomplished by allowing the multiplication rate on each slip system to include contributions from the plastic strain rates of the two coplanar slip systems. The resulting constitutive model successfully captures the strain hardening rate dependence on the loading orientation as predicted by the DDD simulations, which is also consistent with existing experiments.

36 MATERIALS SCIENCE↗

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamics and chemical mode analysis of plasma thermal-chemical instability

The stability of the weakly ionized plasma and the transition from a stable homogeneous discharge to unstable filaments play an important role in gas laser physics, plasma-assisted combustion, chemical reforming, and material synthesis. Here, theoretical stability analysis and thermal-chemical mode analysis were performed to understand the mechanism of plasma thermal-chemical instability by using a zero-dimensional plasma system with both simplified and detailed chemical kinetics of H 2 /O 2 /N 2 mixtures. The plasma dynamic and kinetic models accounted for multiple physical mechanisms in the chemically-reactive weakly ionized plasma, including ionization, attachment/detachment, recombination, vibrational and electronic energy relaxation, convective and diffusive species/heat removal, Joule heating, and detailed chemical kinetics. An analytical criterion and the explosive mode species/temperature pointers were formulated while the representative active species were identified for different thermal-chemical modes. The results showed that in addition to the classical thermal-ionization mechanism, various chemical modes from chemical heat imbalance and elementary kinetics significantly modified the time dynamics and the stability of the weakly ionized plasma. The present analysis provides insights and guidance to control plasma instability using chemical kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generalized Analysis of Dynamic Flow Fouling on Heat Transfer Surfaces

Several energy applications involve flow of untreated heat transfer fluids that contain trace amounts of mineral salts, such as calcium carbonate and calcium sulfate, which show reverse solubility with temperature. The gradual accretion of these salts onto the heat transfer surfaces causes fouling that degrades thermal performance over time. While fouling of bare surfaces has been extensively studied, fouling of nonwetting superhydrophobic or liquid-infused surfaces is less understood, especially under dynamic flow conditions. This study presents a systematic experimental analysis to elucidate dynamic flow fouling of calcium carbonate and calcium sulfate on smooth, superhydrophobic and liquid-infused surfaces in a generalized manner. Here, copper tubes with modified inner surface wettability are fabricated and subjected to forced convection at different flow rates and fouling salt concentrations. Fouling characteristics of the two salts on the different surfaces are examined through microstructural analysis, and quantified in terms of asymptotic fouling resistance that measures the total extent of fouling, and a threshold time, which correlates to the onset of significant fouling. Using Taguchi design of experiments, the study reports, for the first time, closed-form analytical relationships for the asymptotic fouling thermal resistance and threshold time on the Reynolds number, dimensionless concentration of fouling agent and dimensionless infused liquid viscosity that represents the different surface types in a unified manner. It is shown that the analytical model accurately predicts the fouling characteristics in up to about 97% of the data, based on which contour maps of optimum surface designs for minimizing fouling resistance are presented. The results provide information on designing surfaces for reduced fouling or for estimating the fouling characteristics for a given surface in applications.

36 MATERIALS SCIENCE↗

A Novel Experimental Approach to Understand the Transport of Nanodrugs

Nanoparticle-based drugs offer attractive advantages like targeted delivery to the diseased site and size and shape-controlled properties. Therefore, understanding the particulate flow of the nanodrugs is important for effective delivery, accurate prediction of required dosage, and developing efficient drug delivery platforms for nanodrugs. In this study, the transport of nanodrugs including flow velocity and deposition is investigated using three model metal oxide nanodrugs of different sizes including iron oxide, zinc oxide, and combined Cu-Zn-Fe oxide synthesized via a modified polyol approach. The hydrodynamic size, size, morphology, chemical composition, crystal phase, and surface functional groups of the water-soluble nanodrugs were characterized via dynamic light scattering, transmission electron microscopy, scanning electron microscopy-energy dispersive X-ray, X-ray diffraction, and fourier transform infrared spectroscopy, respectively. Two different biomimetic flow channels with customized surfaces are developed via 3D printing to experimentally monitor the velocity and deposition of the different nanodrugs. A diffusion dominated mechanism of flow is seen in size ranges 92 nm to 110 nm of the nanodrugs, from the experimental velocity and mass loss profiles. The flow velocity analysis also shows that the transport of nanodrugs is controlled by sedimentation processes in the larger size ranges of 110–302 nm. However, the combined overview from experimental mass loss and velocity trends indicates presence of both diffusive and sedimentation forces in the 110–302 nm size ranges. It is also discovered that the nanodrugs with higher positive surface charges are transported faster through the two test channels, which also leads to lower deposition of these nanodrugs on the walls of the flow channels. The results from this study will be valuable in realizing reliable and cost-effective in vitro experimental approaches that can support in vivo methods to predict the flow of new nanodrugs.

36 MATERIALS SCIENCE↗

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗

Weather Research and Forecasting Model Sensitivity Comparisons for Warm Season Convective Initiation

Mesoscale weather conditions can significantly affect the space launch and landing operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). During the summer months, land-sea interactions that occur across KSC and CCAFS lead to the formation of a sea breeze, which can then spawn deep convection. These convective processes often last 60 minutes or less and pose a significant challenge to the forecasters at the National Weather Service (NWS) Spaceflight Meteorology Group (SMG). The main challenge is that a "GO" forecast for thunderstorms and precipitation is required at the 90 minute deorbit decision for End Of Mission (EOM) and at the 30 minute Return To Launch Site (RTLS) decision at the Shuttle Landing Facility. Convective initiation, timing, and mode also present a forecast challenge for the NWS in Melbourne, FL (MLB). The NWS MLB issues such tactical forecast information as Terminal Aerodrome Forecasts (TAFs), Spot Forecasts for fire weather and hazardous materials incident support, and severe/hazardous weather Watches, Warnings, and Advisories. Lastly, these forecasting challenges can also affect the 45th Weather Squadron (45 WS), which provides comprehensive weather forecasts for shuttle launch, as well as ground operations, at KSC and CCAFS. The need for accurate mesoscale model forecasts to aid in their decision making is crucial. Both the SMG and the MLB are currently implementing the Weather Research and Forecasting Environmental Modeling System (WRF EMS) software into their operations. The WRF EMS software allows users to employ both dynamical cores - the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM). There are also data assimilation analysis packages available for the initialization of the WRF model- the Local Analysis and Prediction System (LAPS) and the Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS). Having a series of initialization options and WRF cores, as well as many options within each core, provides SMG and NWS MLB with a lot of flexibility. It also creates challenges, such as determining which configuration options are best to address specific forecast concerns. The goal of this project is to assess the different configurations available and to determine which configuration will best predict warm season convective initiation in East-Central Florida. Four different combinations of WRF initializations will be run (ADAS-ARW, ADAS-NMM, LAPS-ARW, and LAPS-NMM) at a 4-km resolution over the Florida peninsula and adjacent coastal waters. Five candidate convective initiation days using three different flow regimes over East-Central Florida will be examined, as well as two null cases (non-convection days). Each model run will be integrated 12 hours with three runs per day, at 0900, 1200, and 1500 UTe. ADAS analyses will be generated every 30 minutes using Level II Weather Surveillance Radar-1988 Doppler (WSR-88D) data from all Florida radars to verify the convection forecast. These analyses will be run on the same domain as the four model configurations. To quantify model performance, model output will be subjectively compared to the ADAS analyses of convection to determine forecast accuracy. In addition, a subjective comparison of the performance of the ARW using a high-resolution local grid with 2-way nesting, I-way nesting, and no nesting will be made for select convective initiation cases. The inner grid will cover the East-Central Florida region at a resolution of 1.33 km. The authors will summarize the relative skill of the various WRF configurations and how each configuration behaves relative to the others, as well as determine the best model configuration for predicting warm season convective initiation over East-Central Florida.

Watson, Leela R.↗

Aerosol Particle Deposition on a Spent Nuclear Fuel Assembly Spacer Grid

The flow and particle deposition patterns on surfaces in an idealized spacer grid for a 17x17 pressurized water reactor (PWR) assembly in a spent fuel canister are modeled using computational fluid dynamics (CFD) with laminar flow. The effects of gravitational settling, non-Stokesian flow, and particle slip are first rigorously analyzed. From the analysis, non-Stokesian effects and slip may be neglected for the particle sizes and conditions expected in a canister. For particles that do not settle out, a swirling flow pattern at the corners of a spacer grid channel directs particles to the leeward side of the flow vanes where much of the deposition occurs. Particle deposition increases with increasing particle diameter. Deposition also increases with decreasing flow velocity as this provides more time for particles to settle and deposit on the leeward side of the flow vanes. The fraction of particles that are transmitted through a spacer grid is determined as a function of inlet gas velocity and particle diameter by running the CFD calculation for each set of conditions and for each particle diameter. Curve fits of the transmission curve as a function of particle diameter for a specified spacer grid and flow velocity are applied to a lognormal particle mass density function for the inlet particles. The resulting mass density function and aerosol mass fraction that passes through the spacer grid can be determined analytically without resorting to numerical iteration. A sample calculation of the analytical solution is demonstrated for a lognormal particle mass density function.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Novel Applications of Magnetic Fields for Fluid Flow Control and for Simulating Variable Gravity Conditions

Static and dynamic magnetic fields have been used to control convection in many materials processing applications. In most of the applications, convection control (damping or enhancement) is achieved through the Lorentz force that can be tailored to counteract/assist dominant system flows. This technique has been successfully applied to liquids that are electrically conducting, such as high temperature melts of semiconductors, metals and alloys, etc. In liquids with low electrical conductivity such as ionic solutions of salts in water, the Lorentz force is weak and hence not very effective and alternate ways of flow control are necessary. If the salt in solution is paramagnetic then the variation of magnetic susceptibility with temperature and/or concentration can be used for flow control. For thermal buoyancy driven flows this can be accomplished in a temperature range below the Curie point of the salt. The magnetic force is proportional to the magnetic susceptibility and the product of the magnetic field and its gradient. By suitably positioning the experiment cell in the magnet, system flows can be assisted or countered, as desired. A similar approach can be extended to diamagnetic substances and fluids but the required magnetic force is considerably larger than that required for paramagnetic substances. The presentation will provide an overview of work to date on a NASA fluid physics sponsored project that aims to test the hypothesis of convective flow control using strong magnetic fields in protein crystal growth. The objective is to understand the nature of the various forces that come into play, delineate causative factors for fluid flow and to quantify them through experiments, analysis, and numerical modeling. The seminar will report specifically on the experimental results using paramagnetic salts and solutions in magnetic fields and compare them to analytical predictions. Applications of the concept to protein crystallization studies will be discussed. The use of strong magnetic fields for terrestrially simulating variable gravity environments and applications supporting the NASA Exploration Initiative will also be briefly discussed.

Ramachandran, N.↗

Corrosion Prevention in Copper Combustion Chamber Liners of Liquid Oxygen/Methane Booster Engines

Prior research conducted by United Technologies Research Center and Rockwell International Rocketdyne Division encountered severe copper corrosion while flowing hydrocarbon fuels through copper cooling tubes. These experimental results have very important implications for the development of low-cost, regeneratively cooled, oxygen/hydrocarbon bipropellant booster engines. Aerojet TechSystems undertook a program sponsored by NASA/Lewis Research Center in 1986 with two objectives: (1) to define the corrosive interaction process that occurs between hydrocarbon fuels and copper combustion chamber liner materials, using both static and dynamic test methods, and (2) to identify and demonstrate protective measures against this corrosive process. The results of the first program task demonstrated that the most damaging corrosive process was caused by trace amounts of sulfur-containing impurities in the hydrocarbon fuel which react with the copper chamber liner material to form cuprous sulfide (Cu2S), as reported in Paper No. AlAA 88-3215. Experiments were conducted which demonstrated this corrosive process with Mil-Spec RP-1, propane, and methane fuels in contact with OFHC, NASA-Z, and ZrCu coppers. The preliminary results of the second program task demonstrated that electrodeposited gold and platinum coatings held promise of preventing corrosion of copper chamber liners by methane contaminated with sulfur-containing impurities, as reported in Paper No. AIAA 89-2738. Experiments were conducted which demonstrated greatly reduced corrosion of copper by hydrogen sulfide (H2S) and methyl mercaptan (CH3SH) as measured impurities in methane. This paper presents the concluding, confirmatory results from the second program task. The efficacy of a protective gold coating on NASA-Z copper cooling channel walls was conclusively demonstrated. Tests were conducted in the Aerojet Carbothermal Test Facility which provide realistic simulations of booster engine cooling channel conditions such as temperature, pressure, flow velocity, and heat flux. Dynamic test series were performed with methane containing 5 and 10 ppm (by volume) hydrogen sulfide. Gold-plated and unplated copper alloy specimens were evaluated. The tests demonstrated that gold coatings were effective in preventing corrosion of the copper. Posttest metallographic examinations of the unplated specimens showed severe corrosion as a result of reaction with the sulfur containing contaminant in the fuel. In contrast, posttest metallographic examination of the gold-coated specimens showed no corrosion under similar operating conditions. The findings from metallurgical and chemical analyses of the copper specimens, chemical analysis of the methane, and description of all test conditions are provided. Conclusions are reached as to the efficacy of the plating methods and the effectiveness of the metal coatings in the prevention of corrosion in hydrocarbon-fueled booster engine combustion chamber liners.

S D Rosenberg↗

Rayleigh-Taylor Instability: A Status Review of Experimental Designs and Measurement Diagnostics

The focus of experiments and the sophistication of diagnostics employed in Rayleigh-Taylor instability (RTI) induced mixing studies have evolved considerably over the past seven decades. The first theoretical analysis by Taylor and the two-dimensional experimental results by Lewis on RTI in 1950 examined single-mode RTI using conventional imaging techniques. Over the next 70 years, several experimental designs have been used to creating an RTI unstable interface between two materials of different densities. These early experiments though innovative, were arduous to diagnose and provided little information on the internal, turbulent structure and initial conditions of the RT mixing layer. Coupled with the availability of high-fidelity diagnostics, the experiments designed and developed in the last three decades allow detailed measurements of various turbulence statistics that have allowed broadly to validate and verify late-time nonlinear models and mix-models for buoyancy-driven flows. Besides, they have provided valuable insights to solve several long-standing disagreements in the field. Herein, this review serves as an opportunity to discuss the understanding of the RTI problem and highlight valuable insights gained into the RTI driven mixing process with a focus on low to high Atwood number (> 0.1) experiments.

42 ENGINEERING↗

Crystallization history of the 1984 Mauna Loa lava flow

During a 3-week eruption in 1984, Mauna Loa produced vent lavas that increased in crystallinity from less than 1 to 30%, and 27-km-long flows that increased in crystallinity as they moved downstream. We examined the crystallization history of these lavas using crystal size distribution (CSD) analysis to study the rates of crystallization, viscosity increase, and latent heating. Typical average growth and nucleation rates were 5 x 10(exp -9) cm/s and 5/cu cm/s for microphenocrysts (20- to 500-micron size crystals nucleated in the rift zone) and 5 x 10(exp -8) cm/s and 5 x 10(exp 4)/cu cm/s for microlites (1- to 20-micron size crystals nucleated in the channel). These crystallization rates are high compared with those found in other CSD studies of igneous rocks, probably due to highly nonequilibrium conditions brought on by rapid degassing in the rift zone and cooling in the lava channel. Growth and nucleation rates decreased with time at the vent and with distance downstream. The maximum downstream total crystallinity measured is 39% (25% microlites, 14% microphenocrysts) in a quenched sample 14 km from the vent. Growth and nucleation rates cannot be calculated for postemplacement samples, but they place upper limits of 53-58% on the amount of crystallization in the channel 9-20 km from the vent. Crystallization could have been mostly responsible for the 10(exp 5)-fold downstream increase in apparent viscosity, although degassing and increasing incorporation of solid lava fragments also contributed. Another effect of crystallization on the lava flow was the sizeable latent heating (0.01 J/g/s over the first half of the flow length, if the crystallinity of downstream quench samples is representative of the hot fluid core), which may have been counteracted by entrainment of cooler material. Measurements of crystallization are shown to be crucial in the study of lava flow emplacement dynamics.

Crisp, Joe↗

Synoptic and chemical evolution of the Antarctic vortex in winter and spring, 1987

The dynamical evolution of the vortex at least up to 50 mb is dominated by synoptic scales in the troposphere. In particular, there is a clear response when poleward extension of tropospheric anticyclones from latitudes of 40 and 50 S to 70 and 80 S occurs. This response is evident in isentropic potential vorticity maps, TOMS ozone fields and SAM II polar stratospheric clouds. An important feature of the high latitude Southern Hemisphere lower stratosphere is a transition at potential temperatures in the 390 to 400 K range. This transition, the vortopause, is clearly marked in aircraft profiles of O3, H2O, N2O and ClO at latitudes 68 to 72 S near the Antarctic peninsula, and also cross-sections of potential vorticity and potential temperature; above it, the isopleths are more closely spaced than below it. The aircraft measurements of H2O, O3, NO sub y, N2O, ClO, and whole air data are examined in material coordinates, theta and P sub theta, backed up by trajectory analysis. The evolution of the chemical mixing ratios is examined in these coordinates as a function of time from mid-August to late September. Conclusions are drawn about the rates of change and their causes. The meteorological and aircraft data are examined for evidence of the following kinds of motion with respect to the vortex: ingress of air aloft, subsidence, peeling off of air to lower latitudes, and folding of the vortopause. Conclusions are presented regarding the evidence for a chemical sink of ozone above and below theta = 400 K, and whether the vortex has a mass flow through it, or if the chemical sink operates on a fixed mass of air. Implications for mid-latitudes are briefly considered in the light of the conclusions.

Tuck, A. F.↗

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and an efficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad midwall temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. There are many development shortcomings in DASSH detailed in this document, but this version is functional for modern analysis needs. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Effects of Traveling Magnetic Field on Dynamics of Solidification

The Lorentz body force induced in electrically conducting fluids can be utilized for a number of materials processing technologies. An application of strong static magnetic fields can be beneficial for damping convection present during solidification. On the other hand, alternating magnetic fields can be used to reduce as well as to enhance convection. However, only special types of time dependent magnetic fields can induce a non-zero time averaged Lorentz force needed for convection control. One example is the rotating magnetic field. This field configuration induces a swirling flow in circular containers. Another example of a magnetic field configuration is the traveling magnetic field (TMF). It utilizes axisymmetric magnetostatic waves. This type of field induces an axial recirculating flow that can be advantageous for controlling axial mass transport, such as during solidification in long cylindrical tubes. Incidentally, this is the common geometry for crystal growth research. The Lorentz force induced by TMF can potentially counter-balance the buoyancy force, diminishing natural convection, or even setting up the flow in reverse direction. Crystal growth process in presence of TMF can be then significantly modified. Such properties as the growth rate, interface shape and macro segregation can be affected and optimized. Melt homogenization is the other potential application of TMF. It is a necessary step prior to solidification. TMF can be attractive for this purpose, as it induces a basic flow along the axis of the ampoule. TMF can be a practical alloy mixing method especially suited for solidification research in space. In the theoretical part of this work, calculations of the induced Lorentz force in the whole frequency range have been completed. The basic flow characteristics for the finite cylinder geometry are completed and first results on stability analysis for higher Reynolds numbers are obtained. A theoretical model for TMF mixing is also developed. In the experimental part, measurements of flow induced by TMF in a column of mercury (Hg) are presented. Also, an alloy mixing of Bi-Sn of the eutectic composition is demonstrated. A traveling magnetic field of 4mT at 3kHz applied for 120 minutes is found to be sufficient to homogenize an alloy enclosed in a 1cm diameter and 12 cm long tube.

Source record↗

EVA Swab Kit: Tools and Techniques for Collecting Aseptic Samples from Crewed Space Missions

Introduction: When we send humans to search for life on other planets, we'll need to know what we brought with us versus what may already be there. To ensure our crewed spacecraft meet planetary protection requirements—and to protect our science from human contamination—we'll need to assess and verify whether micro-organisms may be leaking/venting from our spacesuits. This requires collecting samples under Extravehicular Activity (EVA) conditions. Detailed, systematic research on forward contamination from robotic spacecraft has been steadily progressing since the Viking missions, but systematic studies of contamination from space suits has not been conducted in many years. The modern EMU (Extravehicular Mobility Unit) suit used by NASA is designed to leak at rates as high as 100 cc/min. Before humans land on Mars there is a critical need to understand the types and quantities of microbes that could be introduced via space suits. The Human Forward Contamination Assessment team at NASA’s Johnson Space Center (JSC) has developed a prototype EVA swab tool [1,2,3,4] designed for use in space to sample cleaned and uncleaned space suits to determine the present day microbial load and eventually the rate of leakage. The ability to assess microbial leakage early in advanced space suit and life support system design cycles will help avoid costly hardware redesign later. Test Objectives: The primary objective of EMU testing was to characterize the type of micro-organisms typically found on or near selected suit pressure joints under suit differential pressure conditions. Most human-borne microbes can fit through a 0.5 to 1.0 µm gap. Knowing which joints are more likely to leak will inform hardware design decisions. Knowing which types of micro-organisms may leak from EVA suits provides a basis for subsequent studies to characterize the viability of those organisms under destination conditions, as well as how far they might spread through natural or human-influenced processes. That data, in turn, will inform exploration mission operations and hardware design. The secondary objective of testing was to evaluate the interface between a fully suited test subject and the EVA swab tool at vacuum. Bulky EVA suits can restrict movement and limit visibility through the helmet visor. Fully suited testing is important for identifying tool design issues prior to flight. At exploration destinations, such as Mars, suited crew may be required to periodically sample their suits as part of an environmental monitoring protocol. Suit Microbial Sampling Results: This report details results of microbial swabs collected from current flight suit configurations worn by crew members assigned to upcoming ISS expedition missions as well as swabs collected from prototype suits intended for use on the Orion spacecraft. These tests were intended to characterize the types of contaminants found on flight suits under current, typical handling conditions. No attempt was made to change suit handling procedures, provide additional sterilization, or to limit typical potential contaminant sources. Using culture based techniques, we cultivated 235 CFU (colony forming units) comprised of 26 bacterial species and one fungal species on the outside of the suits. The fungal species and 14 of the bacterial species were unique to the suit surfaces and were not detected in any of the background samples collected within the chambers. We sequenced 755,434 ribosomal fragments on all of the suit surfaces from swab samples. 557,016 of these sequences represent DNA that survived at least 4 hours at vacuum. These sequences formed 2,464 OTU's (Operational Taxonomic Units, 97% similarity) showing low diversity in the samples. The most abundant sequences that survived vacuum belong to the genera Staphyloccocus, Ralstona, Bacillus and Rhodobacter all of which are common to the human microbiome. [5] See Danko et al., (2021) for more complete details of these first analyses. Further analysis of EVA suit materials with respect to the efficacy of various cleaning protocols and engineered containment solutions is planned to inform suit design for NASA’s Artemis Moon to Mars program crew testing. Swab Tool Function Results: The kit was demonstrated for fit and function in suited subject vacuum tests to determine how well the tool worked as an aseptic microbial sampling device as well as to identify any design elements that could be upgraded for EVA task specific improvement. It was found that sample acquisition efficacy could be enhanced by redesign of the sample canister to end-effector interface. Several modifications of the sample caddy assemblies to optimize EVA safety and functionality were also identified. Consequently, fabrication of the redesigned sample canister to end-effector assembly interfaces and and the sample caddy assemblies are required. Fabrication of sixteen flight sample canister assemblies (8 per each of two EVA Swab Kits) and two sample caddy assemblies are in process to be followed by hardware testing and certification to produce two flight-certified EVA Swab Kits for transport to ISS no earlier than summer of 2022. Sampling Strategy: The International Space Station is an ideal testbed for systematic studies of contamination from crewed vehicles since it has been continuously occupied for 20 years and exposed to non-terrestrial conditions. We will sample the exterior of the ISS during EVA using a purpose-built swab tool capable of maintaining sterility while undergoing temperature changes from -151 to +121°C under hard vacuum. Prior to each EVA, the project team will work with ISS mission managers to identify precise sampling locations, which will vary by EVA based on the translation paths and worksites scheduled for that particular EVA. Ideally, translation path handrails and areas near ECLSS (Environmental Control and Life Support System) external vent openings on a spacecraft would be assessed. There are currently more than a dozen ECLSS external vents on the ISS. Some are connected to systems that vent waste products, while others are intended to equalize cabin pressure. As EVA opportunity allows, microbial samples from any of these external vents would provide a valuable data point, though some will be more useful than others. Four criteria have been identified to help prioritize sampling sites near vents: • EVA Accessibility: To minimize cost, it is desired to piggy-back onto a planned EVA. Therefore, the sampling location must be readily accessible by an EVA crew • Type of Vented Products: Vent products that have been in direct contact with crew, such as cabin air, are more likely to contain microorganisms than vent products associated with isolated systems, such as experiment module combustion products. • Mass of Vented Products: Higher-flow vents are more likely to contain detectible levels of microbial contaminants than lower-flow vents. • Local Environment: Sample locations with relatively benign local conditions, such as warm surfaces shielded from direct ultraviolet (UV) radiation exposure, may be more likely to support microbial growth than locations with harsher local environmental conditions. Because EVA accessibility is the most important criteria, the proposal team worked with an astronaut and flight controllers using the Dynamic Onboard Ubiquitous Graphics (DOUG) tool. The DOUG virtual environment allows an operator to “fly” around the current ISS vehicle configuration to assess EVA translation paths, attach points, and keep-out zones. While analysis on station or rapid return to Earth would be preferable, samples collected from the exterior of the ISS have already been exposed to temperature variations between -157 and +121 °C as well as hard vacuum. Therefore, they should be fairly stable and robust. We hypothesize that samples collected from the ISS exterior could be stored for up to 6 months at -80°C without degradation. Sample canisters will be returned to Earth while frozen at -80°C for analysis, and sterilized canisters can be re-flown back to ISS to support additional sampling opportunities Relevance to NASA Exploration Objectives: These data will allow us to identify new or improved methods, technologies, and procedures for spacecraft sterilization and leakage mitigation to minimize the amount of contamination introduced to the environment by human explorers. This work is funded by NASA research grant: NNH18ZDA001N-PPR References: [1] Bell, M.S. et al. (2015) LPS XLVI, Abst. #1832 [2] Rucker et al. (2018) 42nd COSPAR (PPP.3) [3] Bell, M.S. et al. (2019) Mars Extant Life Conference, Abst. #5096.[4] Bell, M.S. et al., (2020) 43rd COSPAR (BO.2).[5] Danko D, et.al.,(2021)Front.Microbiol.12:608478.

Mary Suzanne Bell↗

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗