Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular simulation”

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 667 records · Page 37

A study of molecular contamination

Contamination kinetics of typical space materials were studied using several techniques, including ellipsometry. Data were obtained on contaminant deposition rates, sticking coefficients, and reevaporation of the contamination. Materials studied included silicone sealants, a thermal control paint, and a graphite laminate as contaminant sources, and gold coated surfaces, a thermal control paint and a second surface mirror material as the contaminees. Both deposition and reevaporation rates were found to be nonlinear with time, and followed an exponential decay. A hypothesis is presented that the nonlinear behavior was a result of disproportionation of the molecular weight fractions of the polymers. The measured sticking coefficients were found to vary substantially with temperature over narrow ranges. Simulated space radiation was found to reduce the evaporation rates of contaminants by several orders of magnitude.

Allen, T. H.↗

A comparative study of prebiotic and present day translational models

It is generally recognized that the understanding of the molecular basis of primitive translation is a fundamental step in developing a theory of the origin of life. However, even in modern molecular biology, the mechanism for the decoding of messenger RNA triplet codons into an amino acid sequence of a protein on the ribosome is understood incompletely. Most of the proposed models for prebiotic translation lack, not only experimental support, but also a careful theoretical scrutiny of their compatibility with well understood stereochemical and energetic principles of nucleic acid structure, molecular recognition principles, and the chemistry of peptide bond formation. Present studies are concerned with comparative structural modelling and mechanistic simulation of the decoding apparatus ranging from those proposed for prebiotic conditions to the ones involved in modern biology. Any primitive decoding machinery based on nucleic acids and proteins, and most likely the modern day system, has to satisfy certain geometrical constraints. The charged amino acyl and the peptidyl termini of successive adaptors have to be adjacent in space in order to satisfy the stereochemical requirements for amide bond formation. Simultaneously, the same adaptors have to recognize successive codons on the messenger. This translational complex has to be realized by components that obey nucleic acid conformational principles, stabilities, and specificities. This generalized condition greatly restricts the number of acceptable adaptor structures.

Rein, R.↗

Research in Computational Astrobiology

We present results from several projects in the new field of computational astrobiology, which is devoted to advancing our understanding of the origin, evolution and distribution of life in the Universe using theoretical and computational tools. We have developed a procedure for calculating long-range effects in molecular dynamics using a plane wave expansion of the electrostatic potential. This method is expected to be highly efficient for simulating biological systems on massively parallel supercomputers. We have perform genomics analysis on a family of actin binding proteins. We have performed quantum mechanical calculations on carbon nanotubes and nucleic acids, which simulations will allow us to investigate possible sources of organic material on the early earth. Finally, we have developed a model of protobiological chemistry using neural networks.

Chaban, Galina↗

Comparisons of CFD Simulations of Icing Wind Tunnel Clouds with Experiments Conducted at the NASA Propulsions Systems Laboratory

This paper evaluates simulation predictions against experimental test data of icing clouds that were produced during 2018 ice crystal icing physics tests conducted at the NASA Propulsion Systems Laboratory icing wind tunnel. Aero-thermal and cloud parameters are set and known upstream at the tunnel inlet and spray system, but change as the cloud and air thermodynamically interact as the flowing masses reach the tunnel test section. Utilizing the ANSYS Fluent Discrete Phase Model function, 3D computational fluid dynamics (CFD) simulations were performed, capturing the thermodynamic interactions between the test parameters, and providing predictions of the aero thermal and cloud conditions at the tunnel test section. Simulation predictions were compared with test data that were measured at the tunnel exit plane. Evaluations focused on the cloud concentration (total water content), humidity content, and air temperature. CFD simulation predictions showed areas of agreement and disagreement. Simulations showed that cloud concentration profiles at the test section are strongly related to the initial spray nozzle pattern used at the tunnel inlet. Experimental data suggest that greater dispersion of the cloud occurred as the simulated cloud predicted areas of high and low cloud concentration compared to test data profiles. Simulations, however, captured the magnitude and location of the change in humidity content and the change in air temperature due to the presence of the cloud reasonably well, when compared to test data. This result would suggest that while the ANSYS simulation did not fully predict the spreading of the cloud as measured during experiment, it did capture evaporation and the molecular movements of air and vapor relatively well.

CFD↗

Comparisons of CFD Simulations of Icing Wind Tunnel Clouds with Experiments Conducted at the NASA Propulsions Systems Laboratory

This paper evaluates simulation predictions against experimental test data of icing clouds that were produced during 2018 ice crystal icing physics tests conducted at the NASA Propulsion Systems Laboratory icing wind tunnel. Aero-thermal and cloud parameters are set and known upstream at the tunnel inlet and spray system, but change as the cloud and air thermodynamically interact as the flowing masses reach the tunnel test section. Utilizing the ANSYS Fluent Discrete Phase Model function, 3D computational fluid dynamics (CFD) simulations were performed, capturing the thermodynamic interactions between the test parameters, and providing predictions of the aero thermal and cloud conditions at the tunnel test section. Simulation predictions were compared with test data that were measured at the tunnel exit plane. Evaluations focused on the cloud concentration (total water content), humidity content, and air temperature. CFD simulation predictions showed areas of agreement and disagreement. Simulations showed that cloud concentration profiles at the test section are strongly related to the initial spray nozzle pattern used at the tunnel inlet. Experimental data suggest that greater dispersion of the cloud occurred as the simulated cloud predicted areas of high and low cloud concentration compared to test data profiles. Simulations, however, captured the magnitude and location of the change in humidity content and the change in air temperature due to the presence of the cloud reasonably well, when compared to test data. This result would suggest that while the ANSYS simulation did not fully predict the spreading of the cloud as measured during experiment, it did capture evaporation and the molecular movements of air and vapor relatively well.

CFD↗

Ambiguities in the identification of giant molecular cloud complexes from longitude-velocity diagrams

Techniques which use longitude-velocity diagrams to identify molecular cloud complexes in the disk of the Galaxy are investigated by means of model Galactic disks generated from N-body cloud-particle simulations. A procedure similar to the method used to reduce the low-level emission in Galactic l-v diagrams is employed to isolate complexes of emission in the model l-v diagram (LVCs) from the 'background'clouds. The LVCs produced in this manner yield a size-line-width relationship with a slope of 0.58 and a mass spectrum with a slope of 1.55, consistent with Galactic observations. It is demonstrated that associations identified as LVCs are often chance superpositions of clouds spread out along the line of sight in the disk of the model system. This indicates that the l-v diagram cannot be used to unambiguously determine the location of molecular cloud complexes in the model Galactic disk. The modeling results also indicate that the existence of a size-line-width relationship is not a reliable indicator of the physical nature of cloud complexes, in particular, whether the complexes are gravitationally bound objects.

Adler, David S.↗

Time-dependent Computational Studies of Premixed Flames in Microgravity

This report describes the research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of NASA Microgravity Science and Applications Program. The primary focus of this research is on investigating fundamental questions concerning the propagation and extinction of premixed flames in earth gravity and in microgravity environments. Our approach is to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. Recently more advanced wall boundary conditions such as isothermal and no-slip have been added to the model. This enables the simulation of flames propagating in more practical systems than before. We have used the numerical simulations to investigate the effects of heat losses and buoyancy forces on the structure and stability of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the interaction between the various processes leading to flame instabilities and extinguishment, and to study the dynamics of cell formation and splitting. Our studies have been able to bring out the differences between upward- and downward-propagating flames and predict the zero-gravity behavior of these flames. The simulations have also highlighted the dominant role of wall heat losses in the case of downward-propagating flames. The simulations have been able to qualitatively predict the formation of multiple cells and the cessation of cell-splitting. Our studies have also shown that some flames in a microgravity environment can be extinguished due to a chemical instability and without any external losses. However, further simulations are needed to more completely understand upward-propagating and zero-gravity flames as well as to understand the potential effect of radiative heat losses.

Kailasanath, K.↗

Simulations of Li+ in Ionic Liquids: Structure, Transport, and Electrochemical Windows

Ionic liquids have been proposed as candidate electrolytes for a number of electrochemical applications. The Li+ solvation structure in these liquids is of central importance to electrolyte properties, like ionic conductivity and electrochemical stability. To this point, we employ simulations at three different size scales to better understand various aspects of the interplay between Li+ solvation structure and dynamics. The smallest systems are Li(Anion)n clusters that are treated with high-accuracy density functional theory (DFT) techniques to provide insight into solvation shell structure through energetics and comparisons to experimental IRRaman spectra. Mid-range sized liquid-phase systems (12-24 ion pairs) are treated with DFT molecular dynamics (MD) to provide temperature-dependent insight into Li+ solvation structure, diffusion, and electrochemical window. The largest systems (144-216 ion pairs) are treated with polarizable MD simulations to evaluate the influence of Li-networks on structure and provide size independent values of transport properties. We perform this procedure on three technologically important ionic liquids and comment on property correlations with solvation structure.

batteries↗

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species With Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for nearequilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion of the heating and is then compared to the total heating measured in flight.

Liechty, Derek S.↗

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species with Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for near-equilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion.

Liechty, Derek S.↗

Tidal coupling with the lower atmosphere (invited review)

The various ways are reviewed in which propagating tidal components excited in the mesophere and below affect the structure of the thermosphere and ionosphere above 100 km. Dynamo effects are not treated here. The physical processes affecting the propagation of upward propagating tides are examined and how they are interrelated in the context of a numerical model. Propagating diurnal and semidiurnal tides which reach thermospheric heights are excited primarily by insolation absorption by tropospheric water vapor (0 to 5 km) and stratospheric/mesospheric ozone (40 to 60 km), respectively. Simulation of these oscillations requires consideration of mean zonal winds and meridional temperature gradients, and the damping effects of turbulent and molecular dissipation, radiative cooling, and ion drag. These effects must be considered on a spherical rotating atmosphere extending from the ground to above 300 km, as they are in the model developed by Forbes depicted schematically.

Forbes, J. M.↗

Role of Oxidative Damage in Radiation-Induced Bone Loss

During prolonged spaceflight, astronauts are exposed to both microgravity and space radiation, and are at risk for increased skeletal fragility due to bone loss. Evidence from rodent experiments demonstrates that both microgravity and ionizing radiation can cause bone loss due to increased bone-resorbing osteoclasts and decreased bone-forming osteoblasts, although the underlying molecular mechanisms for these changes are not fully understood. We hypothesized that excess reactive oxidative species (ROS), produced by conditions that simulate spaceflight, alter the tight balance between osteoclast and osteoblast activities, leading to accelerated skeletal remodeling and culminating in bone loss. To test this, we used the MCAT mouse model; these transgenic mice over-express the human catalase gene targeted to mitochondria, the major organelle contributing free radicals. Catalase is an anti-oxidant that converts reactive species, hydrogen peroxide into water and oxygen. This animal model was selected as it displays extended lifespan, reduced cardiovascular disease and reduced central nervous system radio-sensitivity, consistent with elevated anti-oxidant activity conferred by the transgene. We reasoned that mice overexpressing catalase in mitochondria of osteoblast and osteoclast lineage cells would be protected from the bone loss caused by simulated spaceflight. Over-expression of human catalase localized to mitochondria caused various skeletal phenotypic changes compared to WT mice; this includes greater bone length, decreased cortical bone area and moment of inertia, and indications of altered microarchitecture. These findings indicate mitochondrial ROS are important for normal bone-remodeling and skeletal integrity. Catalase over-expression did not fully protect skeletal tissue from structural decrements caused by simulated spaceflight; however there was significant protection in terms of cellular oxidative damage (MDA levels) to the skeletal tissue. Furthermore, we used an array of countermeasures (Antioxidant diets and injections) to prevent the radiation-induced bone loss, although these did not prevent bone loss, analysis is ongoing to determine if these countermeasure protected radiation-induced damage to other tissues.

oxidative damage↗

High-Performance Computing Optimization for Aladyn – Adaptive Neural Network Molecular Dynamics Mini-Application

This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.

Yamakov, Vesselin I.↗

Direct Simulation Monte Carlo Studies of the Gas Sampling for the VATMOS-SR Mission Concept

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept which would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere, and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus, with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a full picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of atmospheric volatile signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process in order to determine how measured isotopic ratios of noble gases in the sample compare to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept, and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed at Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield a comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters as well as the modeling parameters will be performed.

direct simulation Monte Carlo↗

Ionic Polyimides: New High Performance Polymers for Additive Manufacturing

There is currently a very limited set of engineering polymers that have been demonstrated as viable for use in 3-D printing. Additive manufacturing of custom components will require a much larger array of polymers, especially those with physical, thermal, chemical, and mechanical properties that can be tailor-made. The development of ‘Ionic Polyimides’ offers a solution to this shortage by combining the well understood and widely accepted properties of conventional polyimides, with a new approach to polymer synthesis. Polyimides and polymeric ionic liquids (poly(ILs)) are at the forefront of advanced polymer materials, each with their own set of advantages and disadvantages. While it is clear that more types of polymer materials are needed for fused deposition modeling (FDM) additive manufacturing, there is a need to explore these classes of materials. The synthesis process developed by the Bara Research Group at the University of Alabama allows full control over polymer structure, nanostructure, thermal, electrical, and physical properties making them a prime candidate for use in the additive manufacturing process. Furthermore, the new process allows us to tailor-make a high strength polymer that can be used to fabricate filament feedstock instead of pellets for 3D printing. The primary objective of this proposal is to determine the relationship between molecular structure, physical properties, and performance of ionic polyimides. Further, we seek to determine their utility as materials suitable for additive manufacturing of components used in aerospace vehicles, with an emphasis on characterizing and simulating their thermal behaviors and properties. This proposal addresses the need for fundamental research on a customizable polymer filament feedstock for 3-D printing with tailor-made properties potentially making it superior to the commercial blends offered in industry today. The deliverables for this project are the creation of a database that will detail the relationships between the molecular structure and physical properties for the ionic polyimide of interest (e.g. Tg/Tm (Glass Transition Temperature divided by Melting Point)) relative to different ionic polyimide structures). This new database will provide a “road map” to the development of the first generation of materials and ultimately proof-of-concept.

Jackson, Enrique↗

Spacecraft contamination investigation by direct simulation Monte Carlo - Contamination on UARS/HALOE

A three-dimensional version of the direct simulation Monte Carlo method is adapted to study the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite, with the ultimate goal of estimating the cumulative level of molecular deposits on critical optical surfaces of the Halogen Occultation Experiment. Emphasis is placed on describing the development of preprocessing codes that simplify the incorporation of geometry and outgassing data into the main run code. A novel use of parallel processing software allows use of multiple processors to solve different regions of the overall flowfield. Initial results relating to contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the Halogen Occultation telescope's aperture plane are presented, along with data related to code performance.

Rault, Didier F. G.↗

Synthesis of Large Molecules in Cometary Ice Analogs: Physical Properties Related to Self-Assembly Processes

The combination of realistic laboratory simulations and infrared observations have revolutionized our understanding of interstellar dust and ice-the main component of comets. Since comets and carbonaceous micrometeorites may have been important sources of volatiles and carbon compounds on the early Earth, their organic composition may be related to the origin of life. Ices on grains in molecular clouds contain a variety of simple molecules. The D/H ratios of the comets Hale-Bopp and Hyakutake are consistent with a primarily interstellar ice mixture. Within the cloud and especially in the presolar nebula through the early solar system, these icy grains would have been photoprocessed by the ultraviolet producing more complex species such as hexamethylenetetramine, polyoxymethylenes, and simple keones. We reported at the 1999 Bioastronomy meeting laboratory simulations studied to identify the types of molecules which could have been generated in pre-cometary ices. Experiments were conducted by forming a realistic interstellar mixed-molecular ice (H2O, CH3OH, NH3 and CO) at approximately 10 K under high vacuum irradiated with UV light from a hydrogen plasma lamp. The gas mixture was typically 100:50:1:1, however when different ratios were used material with similar characteristics was still produced. The residue that remained after warming to room temperature was analyzed by HPLC, and by several mass spectrometric methods. This material contains a rich mixture of complex compounds with mass spectral profiles resembling those found in IDPs and meteorites. Surface tension measurements show that an amphiphilic component is also present. These species do not appear in various controls or in unphotolyzed samples. Residues from the simulations were also dispersed in aqueous media for microscopy. The organic material forms 10-40 gm diameter droplets that fluoresce at 300-450 nm under UV excitation. These droplets have a morphology and internal structure which appear strikingly similar to those produced by extracts of the Murchison meteorite. Together, these results suggest a link between organic material photochemically synthesized on the cold grains in dense, interstellar molecular clouds and compounds that may have contributed to the organic inventory of the primitive Earth. For example, the amphiphilic properties of such compounds permit self-assembly into the membranous boundary structures that required for the first forms of cellular life.

Dworkin, Jason P.↗

The VATMOS-SR Mission Concept: DSMC Studies of the Gas Sampling

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept that would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a complete picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O, and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of volatile atmospheric signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process to determine how measured isotopic ratios of noble gases in the sample compared to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open-source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed on Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters, as well as the modeling parameters, will be performed.

direct simulation Monte Carlo↗