Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluids”

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 1,441 records · Page 80

The Prediction and Analysis of Jet Flows and Scattered Turbulent Mixing Noise about Flight Vehicle Airframes

Jet flows interacting with nearby surfaces exhibit a complex behavior in which acoustic and aerodynamic characteristics are altered. The physical understanding and prediction of these characteristics are essential to designing future low noise aircraft. A new approach is created for predicting scattered jet mixing noise that utilizes an acoustic analogy and steady Reynolds-averaged Navier-Stokes solutions. A tailored Green's function accounts for the propagation of mixing noise about the airframe and is calculated numerically using a newly developed ray tracing method. The steady aerodynamic statistics, associated unsteady sound source, and acoustic intensity are examined as jet conditions are varied about a large flat plate. A non-dimensional number is proposed to estimate the effect of the aerodynamic noise source relative to jet operating condition and airframe position.The steady Reynolds-averaged Navier-Stokes solutions, acoustic analogy, tailored Green's function, non-dimensional number, and predicted noise are validated with a wide variety of measurements. The combination of the developed theory, ray tracing method, and careful implementation in a stand-alone computer program result in an approach that is more first principles oriented than alternatives, computationally efficient, and captures the relevant physics of fluid-structure interaction.

Miller, Steven A. E.↗

Robe Development for Electrical Conductivity Analysis in an Electron Gun Produced Helium Plasma

The use of magnetohydrodynamic (MHD) power conversion systems, potentially coupled with a fission power source, is currently being investigated as a driver for an advanced propulsion system, such as a plasma thruster. The efficiency of a MHD generator is strongly dependent on the electrical conductivity of the fluid that passes through the generator; power density increases as fluid conductivity increases. Although traditional MHD flows depend on thermal ionization to enhance the electrical conductivity, ionization due to nuclear interactions may achieve a comparable or improved conductivity enhancement while avoiding many of the limitations inherent to thermal ionization. Calculations suggest that nuclear-enhanced electrical conductivity increases as the neutron flux increases; conductivity of pure He-3 greater than 10 mho/m may be achievable if exposed to a flux greater than 10(exp 12) neutrons/cm2/s.) However, this remains to be demonstrated experimentally. An experimental facility has been constructed at the Propulsion Research Center at the NASA Marshall Space Flight Center, using helium as the test fluid. High energy electrons will be used to simulate the effects of neutron-induced ionization of helium gas to produce a plasma. These experiments will be focused on diagnosis of the plasma in a virtually static system; results will be applied to future tests with a MHD system. Initial experiments will utilize a 50 keV electron gun that can operate at up to a current of 200 micro A. Spreading the electron beam over a four inch diameter window results in an electron flux of 1.5x 10(exp 13) e/sq cm/s. The equivalent neutron flux that would produce the same ionization fraction in helium is 1x10(exp 12) n/sq cm/s. Experiments will simulate the neutron generated plasma modeled by Bitteker, which takes into account the products of thermal neutron absorption in He-3, and includes various ion species in estimating the conductivity of the resulting plasma. Several different probes will be designed and implemented to verify the plasma kinetics model. System parameters and estimated operating ranges are summarized. The predicted ionization fraction, electron density, and conductivity levels are provided in for an equivalent neutron flux of 1x10(exp 12) n/cm2/s. Understanding the complex plasma kinetics throughout a MHD channel is necessary to design an optimal power conversion system for space propulsion applications. The proposed experiments seek to fully characterize the helium plasma and to determine the reliability of each measurement technique, such that they may be applied to more advanced MHD studies. The expected value of each plasma parameter determined from theoretical models will be verified experimentally by several independent techniques to determine the most reliable method of obtaining each parameter. The results of these experiments will be presented in the final paper.

Bragg-Sitton, Shannon M.↗

The Prediction and Analysis of Jet Flows and Scattered Turbulent Mixing Noise About Flight Vehicle Airframes

Jet flows interacting with nearby surfaces exhibit a complex behavior in which acoustic and aerodynamic characteristics are altered. The physical understanding and prediction of these characteristics are essential to designing future low noise aircraft. A new approach is created for predicting scattered jet mixing noise that utilizes an acoustic analogy and steady Reynolds-averaged Navier-Stokes solutions. A tailored Green's function accounts for the propagation of mixing noise about the air-frame and is calculated numerically using a newly developed ray tracing method. The steady aerodynamic statistics, associated unsteady sound source, and acoustic intensity are examined as jet conditions are varied about a large at plate. A non-dimensional number is proposed to estimate the effect of the aerodynamic noise source relative to jet operating condition and airframe position. The steady Reynolds-averaged Navier-Stokes solutions, acoustic analogy, tailored Green's function, non- dimensional number, and predicted noise are validated with a wide variety of measurements. The combination of the developed theory, ray tracing method, and careful implementation in a stand-alone computer program result in an approach that is more first principles oriented than alternatives, computationally efficient, and captures the relevant physics of fluid-structure interaction.

Miller, Steven A.↗

Laplace's equation and the Dirichlet-Neumann map in multiply connected domains

A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.

Greenbaum, A.↗

Development of High Speed Imaging and Analysis Techniques Compressible Dynamics Stall

Dynamic stall has limited the flight envelope of helicopters for many years. The problem has been studied in the laboratory as well as in flight, but most research, even in the laboratory, has been restricted to surface measurement techniques such as pressure transducers or skin friction gauges, except at low speed. From this research, it became apparent that flow visualization tests performed at Mach numbers representing actual flight conditions were needed if the complex physics associated with dynamic stall was to be properly understood. However, visualization of the flow field during compressible conditions required carefully aligned and meticulously reconstructed holographic interferometry. As part of a long-range effort focused on exposing of the physics of compressible dynamic stall, a research wind tunnel was developed at NASA Ames Research Center which permits visual access to the full flow field surrounding an oscillating airfoil during compressible dynamic stall. Initially, a stroboscopic schlieren technique was used for visualization of the stall process, but the primary research tool has been point diffraction interferometry(PDI), a technique carefully optimized for use in th is project. A review of the process of development of PDI will be presented in the full paper. One of the most valuable aspects of PDI is the fact that interferograms are produced in real time on a continuous basis. The use of a rapidly-pulsed laser makes this practical; a discussion of this approach will be presented in the full paper. This rapid pulsing(up to 40,000 pulses/sec) produces interferograms of the rapidly developing dynamic stall field in sufficient resolution(both in space and time) that the fluid physics of the compressible dynamic stall flowfield can be quantitatively determined, including the gradients of pressure in space and time. This permits analysis of the influence of the effect of pitch rate, Mach number, Reynolds number, amplitude of oscillation, and other parameters on the dynamic stall process. When interferograms can be captured in real time, the potential for real-time mapping of a developing unsteady flow such as dynamic stall becomes a possibility. This has been achieved in the present case through the use of a high-speed drum camera combined with electronic circuitry which has resulted in a series of interferograms obtained during a single cycle of dynamic stall; images obtained at the rate of 20 KHz will be presented as a part of the formal presentation. Interferometry has been available for a long time; however, most of its use has been limited to visualization. The present research has focused on use of interferograms for quantitative mapping of the flow over oscillating airfoils. Instantaneous pressure distributions can now be obtained semi-automatically, making practical the analysis of the thousands of interferograms that are produced in this research. A review of the techniques that have been developed as part of this research effort will be presented in the final paper.

Chandrasekhara, M. S.↗

Laboratory simulated hydrothermal alteration of sedimentary organic matter from Guaymas Basin, Gulf of California

High temperature alteration of sedimentary organic matter associated with marine hydrothermal systems involves complex physical and chemical processes that are not easily measured in most natural systems. Many of these processes can be evaluated indirectly by examining the geochemistry of the hydrothermal system in the laboratory. In this investigation, an experimental organic geochemical approach to studying pyrolysis of sedimentary organic matter is applied to the hydrothermal system in the Guaymas Basin, Gulf of California. A general survey of hydrothermal oils and extractable organic matter (bitumen) in hydrothermally altered sediments identified several homologous series of alkanones associated with a high temperature hydrothermal origin. The alkanones range in carbon number from C11 to C30 with no carbon number preference. Alkan-2-ones are in highest concentrations, with lower amounts of 3-, 4-, 5- (and higher) homologs. The alkanones appear to be pyrolysis products synthesized under extreme hydrothermal conditions. Hydrous pyrolysis and confinement pyrolysis experiments were performed to simulate thermally enhanced diagenetic and catagenetic changes in the immature sedimentary organic matter. The extent of alteration was measured by monitoring the n-alkanes, acyclic isoprenoids, steroid and triterpenoid biomarkers, polycyclic aromatic hydrocarbons and alkanones. The results were compared to bitumen extracts from sediments which have been naturally altered by a sill intrusion and accompanied hydrothermal fluid flow. These pyrolysis experiments duplicated many of the organic matter transformations observed in the natural system. Full hopane and sterane maturation occurred after 48 hr in experiments at 330 deg C with low water/rock mass ratios (0.29). A variety of radical and ionic reactions are responsible for the organic compound conversions which occur under extreme hydrothermal conditions. Short duration pyrolysis experiments revealed that a portion of the hydrocarbons generated from kerogen was observed to go through alkene intermediates, and the rate of alkene isomerization was influenced by the ionic strength and catalytic mineral phases. Confinement of the organic pyrolysate to the bulk sediment accelerated the rates of the biomarker epimerization reactions, suggesting that these reactions are influenced strongly by the association of the inorganic matrix, and that the relative rates of some ionic and radical reactions can be influenced by the water/rock ratio during the pyrolysis experiments.

Leif, Roald N.↗

MESSENGER: Exploring Mercury's Magnetosphere

The MESSENGER mission to Mercury offers our first opportunity to explore this planet s miniature magnetosphere since the brief flybys of Mariner 10. Mercury s magnetosphere is unique in many respects. The magnetosphere of Mercury is among the smallest in the solar system; its magnetic field typically stands off the solar wind only - 1000 to 2000 km above the surface. For this reason there are no closed drift paths for energetic particles and, hence, no radiation belts. The characteristic time scales for wave propagation and convective transport are short and kinetic and fluid modes may be coupled. Magnetic reconnection at the dayside magnetopause may erode the subsolar magnetosphere allowing solar wind ions to impact directly the regolith. Inductive currents in Mercury s interior may act to modify the solar wind interaction by resisting changes due to solar wind pressure variations. Indeed, observations of these induction effects may be an important source of information on the state of Mercury s interior. In addition, Mercury s magnetosphere is the only one with its defining magnetic flux tubes rooted in a planetary regolith as opposed to an atmosphere with a conductive ionospheric layer. This lack of an ionosphere is probably the underlying reason for the brevity of the very intense, but short-lived, - 1-2 min, substorm-like energetic particle events observed by Mariner 10 during its first traversal of Mercury s magnetic tail. Because of Mercury s proximity to the sun, 0.3 - 0.5 AU, this magnetosphere experiences the most extreme driving forces in the solar system. All of these factors are expected to produce complicated interactions involving the exchange and re-cycling of neutrals and ions between the solar wind, magnetosphere, and regolith. The electrodynamics of Mercury s magnetosphere are expected to be equally complex, with strong forcing by the solar wind, magnetic reconnection at the magnetopause and in the tail, and the pick-up of planetary ions all driving field-aligned electric currents. However, these field-aligned currents do not close in an ionosphere, but in some other manner. In addition to the insights- into magnetospheric physics offered by study of the solar wind - Mercury system, quantitative specification of the "external" magnetic field generated by magnetospheric currents is necessary for accurate determination of the strength and multi-polar decomposition of Mercury s intrinsic magnetic field. MESSENGER S highly capable instrumentation and broad orbital coverage will greatly advance our understanding of both the origin of Mercury s magnetic field and the acceleration of charged particles in small magnetospheres. In. this article, we review what is known about Mercury s magnetosphere and describe the MESSENGER science team s strategy for obtaining answers to the outstanding science questions surrounding the interaction of the solar wind with Mercury and its small, but dynamic, magnetosphere.

Slavin, James A.↗

A Model for Macromolecular Crystallization

Macromolecular crystallization is a complex process. involving a system which typically has 5 or more components (macromolecule, water, buffer + counter ion, and precipitant). Whereas small molecules have only several well defined contacts in the crystal lattice, macromolecules generally have 10's or even 100's of contacts between molecules. These can range from hydrogen bonds (direct or water-mediated), through van der Waals, hydrophobic, salt bridges, and ion-mediated contacts. The latter interactions are stronger and require some specificity in the molecular alignment, while the others are weaker, more prevalent, and more promiscuous, i.e., can often be readily broken and reformed between other sites. Formation of a consistent, ordered, 3D structure may be impossible in the absence of any or presence of too many strong interactions. Further complicating the process is the inherent structural asymmetry of monomeric single chain macromolecules. The process of crystal nucleation and growth involves the ordered assembly of growth units into a defined 3D lattice. We suggest that for many macromolecules, particularly those that are monomeric, this involves a preliminary solution-phase assembly process into a growth unit having some symmetry prior to addition to the lattice, recapitulating the initial stages of the nucleation process. If this model is correct then fluids and crystal growth models assuming a strictly monodisperse nutrient solution need to be revised. Experimental evidence, based upon face growth rate, AFM, and fluorescence energy transfer data, for a postulated model of the nucleation of tetragonal lysozyme crystals and how it transitions into crystal growth will be presented.

Pusey, Marc L.↗

A Model for Tetragonal Lysozyme Crystal Nucleation and Growth

Macromolecular crystallization is a complex process, involving a system that typically has 5 or more components (macromolecule, water, buffer + counter ion, and precipitant). Whereas small molecules have only a few contacts in the crystal lattice, macromolecules generally have 10's or even 100's of contacts between molecules. These can range from hydrogen bonds (direct or water-mediated), through van der Waals, hydrophobic, salt bridges, and ion-mediated contacts. The latter interactions are stronger and require some specificity in the molecular alignment, while the others are weaker, more prevalent, and more promiscuous, i.e., can be readily broken and reformed between other sites. Formation of a consistent, ordered, 3D structure may be difficult or impossible in the absence of any or presence of too many strong interactions. Further complicating the process is the inherent structural asymmetry of monomeric (single chain) macromolecules. The process of crystal nucleation and growth involves the ordered assembly of growth units into a defined 3D lattice. We suggest that for many macromolecules, particularly those that are monomeric, this involves a preliminary solution-phase assembly process into a growth unit having some symmetry prior to addition to the lattice, recapitulating the initial stages of the nucleation process. If this model is correct then fluids and crystal growth models assuming a strictly monodisperse nutrient solution need to be revised. This model has been developed from experimental evidence based upon face growth rate, AFM, and fluorescence energy transfer data for the nucleation and growth of tetragonal lysozyme crystals.

Pusey, Marc L.↗

Microgravity Flow Boiling Experiments With Liquid-Vapor Mixture Inlet Onboard the International Space Station

This article is part of a series of studies culminating from the multi-objective Flow Boiling and Condensation Experiment (FBCE) onboard the International Space Station, which utilized the Flow Boiling Module (FBM) for experiments between February and July 2022. This study investigates microgravity flow boiling of n-Perfluorohexane with liquid-vapor mixture (two-phase) inlet conditions to the FBM with either one or two opposite walls heated. The FBM’s channel has a rectangular cross-sectional area of 5.0 × 2.5 mm 2 and a heated length of 114.6 mm. Key parameters of interest include mass velocity (180 – 2400 kg/m 2 s), inlet quality (-0.01 – 0.87), inlet pressure (120 – 200 kPa), and heat flux (1.8 W/cm 2 to critical heat flux), and a large database is amassed. The flow is visualized via a high-speed video camera and photographs are recorded at each heating increment to assess the periodic flow patterns within the channel and the near-wall interfacial behavior. Flow patterns are complex and mainly characterized by high- and low-density fronts alternately traversing the channel to yield high- and low-density-dominant periods of boiling. At all operating conditions, high-density fronts are faster during high-density-dominant periods. At low inlet qualities, the flow is annular near the channel inlet with a central vapor core surrounded by an annular liquid layer. Each high-density front having a high liquid fraction leaves a thin liquid layer sheared onto the heated walls. Boiling occurs within the liquid layer and a vapor layer is formed next to the heated wall. Inlet quality and mass velocity most dictate the overall flow patterns followed by heating configuration, and to a much lesser extent, heat flux and inlet pressure. Heat transfer characteristics are assessed via averaged boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and parametric curves of local and averaged heat transfer coefficients. Inlet pressure has an insignificant effect on heat transfer. At similar operating conditions, both the heating configurations yield similar trends and values of heat transfer coefficient and critical heat flux (CHF, slightly higher for single-sided) even though double-sided heating adds twice the heat to the fluid and doubly raises local quality. The heat fluxes required for both onset of nucleate boiling degradation and CHF are larger at high mass velocities and low inlet qualities. For a fixed inlet quality, high mass velocities yield higher average heat transfer coefficients at both lower and higher heat fluxes, while the nucleate boiling regime at intermediate heat fluxes is unaffected. For a fixed mass velocity, higher inlet qualities yield higher and lower average heat transfer coefficients at lower and higher heat fluxes, respectively.

Microgravity↗

The Impact of ARM on Climate Modeling

Climate models are among humanity's most ambitious and elaborate creations. They are designed to simulate the interactions of the atmosphere, ocean, land surface, and cryosphere on time scales far beyond the limits of deterministic predictability, and including the effects of time-dependent external forcings. The processes involved include radiative transfer, fluid dynamics, microphysics, and some aspects of geochemistry, biology, and ecology. The models explicitly simulate processes on spatial scales ranging from the circumference of the Earth down to one hundred kilometers or smaller, and implicitly include the effects of processes on even smaller scales down to a micron or so. The atmospheric component of a climate model can be called an atmospheric global circulation model (AGCM). In an AGCM, calculations are done on a three-dimensional grid, which in some of today's climate models consists of several million grid cells. For each grid cell, about a dozen variables are time-stepped as the model integrates forward from its initial conditions. These so-called prognostic variables have special importance because they are the only things that a model remembers from one time step to the next; everything else is recreated on each time step by starting from the prognostic variables and the boundary conditions. The prognostic variables typically include information about the mass of dry air, the temperature, the wind components, water vapor, various condensed-water species, and at least a few chemical species such as ozone. A good way to understand how climate models work is to consider the lengthy and complex process used to develop one. Lets imagine that a new AGCM is to be created, starting from a blank piece of paper. The model may be intended for a particular class of applications, e.g., high-resolution simulations on time scales of a few decades. Before a single line of code is written, the conceptual foundation of the model must be designed through a creative envisioning that starts from the intended application and is based on current understanding of how the atmosphere works and the inventory of mathematical methods available.

atmospheric circulation↗

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗

Broken Ergodicity in MHD Turbulence in a Spherical Domain

Broken ergodicity (BE) occurs in Fourier method numerical simulations of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence. Although naive statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that low-wave-number coefficients have non-zero mean values that can be very large compared to the associated standard deviation. In other words, large-scale coherent structure (i.e., broken ergodicity) in homogeneous MHD turbulence can spontaneously grow out of random initial conditions. Eigenanalysis of the modal covariance matrices in the probability density functions of ideal statistical theory leads to a theoretical explanation of observed BE in homogeneous MHD turbulence. Since dissipation is minimal at the largest scales, BE is also relevant for resistive magnetofluids, as evidenced in numerical simulations. Here, we move beyond model magnetofluids confined by periodic boxes to examine BE in rotating magnetofluids in spherical domains using spherical harmonic expansions along with suitable boundary conditions. We present theoretical results for 3-D and 2-D spherical models and also present computational results from dynamical simulations of 2-D MHD turbulence on a rotating spherical surface. MHD turbulence on a 2-D sphere is affected by Coriolus forces, while MHD turbulence on a 2-D plane is not, so that 2-D spherical models are a useful (and simpler) intermediate stage on the path to understanding the much more complex 3-D spherical case.

Shebalin, John V.↗

Diversity in C-Xanes Spectra Obtained from Carbonaceous Solid Inclusions from Monahans Halite

Monahans meteorite (H5) contains fluid inclusion- bearing halite (NaCl) crystals [1]. Microthermometry and Raman spectroscopy showed that the fluid in the inclusions is an aqueous brine and they were trapped near 25degC [1]. Their continued presence in the halite grains requires that their incorporation into the H chondrite asteroid was post metamorphism [2]. Abundant solid inclusions are also present in the halites. The solid inclusions include abundant and widely variable organics [2]. Analyses by Raman microprobe, SEM/EDX, synchrotron X-ray diffraction and TEM reveal that these grains include macromolecular carbon similar in structure to CV3 chondrite matrix carbon, aliphatic carbon compounds, olivine (Fo99−59), high- and low-Ca pyroxene, feldspars, magnetite, sulfides, lepidocrocite, carbonates, diamond, apatite and possibly the zeolite phillipsite [3]. Here we report organic analyses of these carbonaceous residues in Monahans halite using C-, N-, and O- X-ray absorption near edge structure (XANES). Samples and Methods: Approximately 100 nm-thick sections were extracted with a focused ion beam (FIB) at JSC from solid inclusions from Monahans halite. The sections were analyzed using the scanning transmission X-ray microscope (STXM) on beamline 5.3.2.2 at the Advanced Light Source, Lawrence Berkeley National Laboratory for XANES spectroscopy. Results and Discussion: C-XANES spectra of the solid inclusions show micrometer-scale heterogeneity, indicating that the macromolecular carbon in the inclusions have complex chemical variations. C-XANES features include 284.7 eV assigned to aromatic C=C, 288.4-288.8 eV assigned to carboxyl, and 290.6 eV assigned to carbonate. The carbonyl features obtained by CXANES might have been caused by the FIB used in sample preparation. No specific N-XANES features are observed. The CXANES spectra obtained from several areas in the FIB sections include type 1&2 chondritic IOM like, type 3 chondritic IOM like, and none of the above. The natures of the macromolecular carbon in the solid inclusions observed by C-XANES are consistent with the previous studies showing that the carbonaceous solid inclusions have not originated from Monahans parent body [1-3], and have various origins, including various chondritic meteorite parent bodies as well as other unknown source(s).

Kebukawa, Y.↗

Mars Ascent Vehicle Hybrid Propulsion Configuration

As part of a Mars Sample Return (MSR) campaign, two Mars Ascent Vehicle (MAV) configurations have been designed in parallel. Each ascent vehicle configuration has a different propulsion system which ultimately leads to two unique vehicle designs. As part of a Preliminary Architecture Assessment (PAA), these vehicle designs were developed to the same level of maturity in order to inform the selection of one of the vehicles as the point of departure design for the campaign. The selection will be made in November 2019. The initial MSR architecture called for a hybrid-based propulsion MAV. This type of propulsion system calls for a solid wax motor that would utilize liquid MON-25 as an oxidizer. Hybrid rocket propulsion allows for more flexibility than traditional solid or liquid propulsion options, and typically benefits from the advantages of both. A hybrid motor can be throttled and shut down easily, and avoids significant risk in manufacturing and handling. On a theoretical level, hybrid motors perform at a higher specific impulse (Isp) than solid motors. The primary disadvantage of hybrid motors comes from additional complexity and significantly less flight heritage and low Technology Readiness Level (TRL). This paper describes the design of the hybrid propulsion configuration. An additional paper will be published describing the design of the solid propulsion configuration1. The hybrid propulsion configuration MAV was developed in 2019 by NASA Marshall Space Flight Center (MSFC) in association with NASA Jet Propulsion Laboratory (JPL). It features a Single Stage to Orbit (SSTO) design with an SP7A solid wax fuel and MON-25 liquid oxidizer. The liquid portion of the vehicle allows for a Liquid Injection Thrust Vector Controller (LITVC) as well as hypergolic propellant additives for ignition. The vehicle was designed to deliver approximately 0.31kg of Martian geological samples to a circular orbit at Mars of 343km at a 25o inclination. Although hybrid propulsion in general has been used on launch vehicles in the past, the integrated vehicle subsystems that operate in conjunction with these propulsion elements do not typically operate in a Martian environment, which in this application can get as cold as -40oC. The PAA advanced the maturity of these subsystems by performing detailed design and analysis on the vehicle with respect to structures and mechanisms, Guidance/Navigation/Control (GNC) systems, avionics, Reaction Control System (RCS), LITVC, thermal environments, and advanced Computational Fluid Dynamics (CFD). This paper will summarize the results of these studies.

Yaghoubi, Darius↗

Finite Element Models for Electron Beam Freeform Fabrication Process

Electron beam freeform fabrication (EBF3) is a member of an emerging class of direct manufacturing processes known as solid freeform fabrication (SFF); another member of the class is the laser deposition process. Successful application of the EBF3 process requires precise control of a number of process parameters such as the EB power, speed, and metal feed rate in order to ensure thermal management; good fusion between the substrate and the first layer and between successive layers; minimize part distortion and residual stresses; and control the microstructure of the finished product. This is the only effort thus far that has addressed computer simulation of the EBF3 process. The models developed in this effort can assist in reducing the number of trials in the laboratory or on the shop floor while making high-quality parts. With some modifications, their use can be further extended to the simulation of laser, TIG (tungsten inert gas), and other deposition processes. A solid mechanics-based finite element code, ABAQUS, was chosen as the primary engine in developing these models whereas a computational fluid dynamics (CFD) code, Fluent, was used in a support role. Several innovative concepts were developed, some of which are highlighted below. These concepts were implemented in a number of new computer models either in the form of stand-alone programs or as user subroutines for ABAQUS and Fluent codes. A database of thermo-physical, mechanical, fluid, and metallurgical properties of stainless steel 304 was developed. Computing models for Gaussian and raster modes of the electron beam heat input were developed. Also, new schemes were devised to account for the heat sink effect during the deposition process. These innovations, and others, lead to improved models for thermal management and prediction of transient/residual stresses and distortions. Two approaches for the prediction of microstructure were pursued. The first was an empirical approach involving the computation of thermal gradient, solidification rate, and velocity (G,R,V) coupled with the use of a solidification map that should be known a priori. The second approach relies completely on computer simulation. For this purpose a criterion for the prediction of morphology was proposed, which was combined with three alternative models for the prediction of microstructure; one based on solidification kinetics, the second on phase diagram, and the third on differential scanning calorimetry data. The last was found to be the simplest and the most versatile; it can be used with multicomponent alloys and rapid solidification without any additional difficulty. For the purpose of (limited) experimental validation, finite element models developed in this effort were applied to three different shapes made of stainless steel 304 material, designed expressly for this effort with an increasing level of complexity. These finite element models require large computation time, especially when applied to deposits with multiple adjacent beads and layers. This problem can be overcome, to some extent, by the use of fast, multi-core computers. Also, due to their numerical nature coupled with the fact that solid mechanics- based models are being used to represent the material behavior in liquid and vapor phases as well, the models have some inherent approximations that become more pronounced when dealing with multi-bead and multi-layer deposits.

Chandra, Umesh↗

Validation of High-Fidelity CFD Simulations for Rocket Injector Design

Computational fluid dynamics (CFD) has the potential to improve the historical rocket injector design process by evaluating the sensitivity of performance and injector-driven thermal environments to the details of the injector geometry and key operational parameters. Methodical verification and validation efforts on a range of coaxial injector elements have shown the current production CFD capability must be improved in order to quantitatively impact the injector design process. This paper documents the status of a focused effort to compare and understand the predictive capabilities and computational requirements of a range of CFD methodologies on a set of single element injector model problems. The steady Reynolds-Average Navier-Stokes (RANS), unsteady Reynolds-Average Navier-Stokes (URANS) and three different approaches using the Large Eddy Simulation (LES) technique were used to simulate the initial model problem, a single element coaxial injector using gaseous oxygen and gaseous hydrogen propellants. While one high-fidelity LES result matches the experimental combustion chamber wall heat flux very well, there is no monotonic convergence to the data with increasing computational tool fidelity. Systematic evaluation of key flow field regions such as the flame zone, the head end recirculation zone and the downstream near wall zone has shed significant, though as of yet incomplete, light on the complex, underlying causes for the performance level of each technique. 1 Aerospace Engineer and Combustion CFD Team Leader, MS ER42, NASA MSFC, AL 35812, Senior Member, AIAA. 2 Professor and Director, Computational Combustion Laboratory, School of Aerospace Engineering, 270 Ferst Dr., Atlanta, GA 30332, Associate Fellow, AIAA. 3 Reilly Professor of Engineering, School of Mechanical Engineering, 585 Purdue Mall, West Lafayette, IN 47907, Fellow, AIAA. 4 Principal Member of Technical Staff, Combustion Research Facility, 7011 East Avenue, MS9051, Livermore, CA 94550, Associate Fellow, AIAA. 5 J. L. and G. H. McCain Endowed Chair, Mechanical Engineering, 104 Research Building East, University Park, PA 16802, Fellow, AIAA. American Institute of Aeronautics and Astronautics 1

Tucker, P. Kevin↗

Kestrel Results at Liftoff Conditions for a Space Launch System Configuration in Proximity to the Launch Tower

Aerodynamic data books for Space Launch System vehicles require databases for the integrated forces and moments and section loads during liftoff and transition to the ascent phase of flight. While the force and moment database can be generated from wind tunnel results, computational analyses are necessary to provide the extensive surface information required to generate proper lineloads. Of the two flight regimes, the liftoff problem is the more costly and complex situation to simulate, as it requires modeling of the vehicle in proximity to the launch tower. The effects of massive separation on the leeward pressure fields of both the tower and vehicle are not well captured with RANS methods, necessitating the use of more advanced methods, such as Delayed Detached Eddy Simulation, in conjunction with computational grids sufficiently refined to resolve the wakes. Details on the computational setup for employing the Kestrel flow solver to address the liftoff problem are presented. The methodology involves the use of independent unstructured near-body grids for the vehicle and the tower, overset by a solution adaptive Cartesian off-body grid. Results from the simulations are compared to experimental results from a test in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel.

Computational fluid dynamics↗