Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle dynamics 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 199 records · Page 11

A Massively Parallel Hybrid Dusty-Gasdynamics and Kinetic Direct Simulation Monte Carlo Model for Planetary Applications

In order to understand the global structure, dynamics, and physical and chemical processes occurring in the upper atmospheres, exospheres, and ionospheres of the Earth, the other planets, comets and planetary satellites and their interactions with their outer particles and fields environs, it is often necessary to address the fundamentally non-equilibrium aspects of the physical environment. These are regions where complex chemistry, energetics, and electromagnetic field influences are important. Traditional approaches are based largely on hydrodynamic or magnetohydrodynamic (MHD) formulations and are very important and highly useful. However, these methods often have limitations in rarefied physical regimes where the molecular collision rates and ion gyrofrequencies are small and where interactions with ionospheres and upper neutral atmospheres are important. At the University of Michigan we have an established base of experience and expertise in numerical simulations based on particle codes which address these physical regimes. The Principal Investigator, Dr. Michael Combi, has over 20 years of experience in the development of particle-kinetic and hybrid kinetichydrodynamics models and their direct use in data analysis. He has also worked in ground-based and space-based remote observational work and on spacecraft instrument teams. His research has involved studies of cometary atmospheres and ionospheres and their interaction with the solar wind, the neutral gas clouds escaping from Jupiter s moon Io, the interaction of the atmospheres/ionospheres of Io and Europa with Jupiter s corotating magnetosphere, as well as Earth s ionosphere. This report describes our progress during the year. The contained in section 2 of this report will serve as the basis of a paper describing the method and its application to the cometary coma that will be continued under a research and analysis grant that supports various applications of theoretical comet models to understanding the inner comae of comets (grant NAGS- 13239 from the Planetary Atmospheres program).

Combi, Michael R.↗

The neutral E region zonal winds during intense postmidnight diffuse aurora - Response to observed particle fluxes

Zonal winds associated with diffuse auroras were simulated using observed fluxes of electron precipitation and a high-resolution time-dependent numerical model described by Waltershcheid et al. (1985). A very strong dynamically unstable E region jet is simulated for an electric field of 50 mV/m, and the results suggest a connection with omega band. These zonal winds are not as strong as those reported by Lyons and Walterscheid (1985), but are much stronger (by a factor of 3 or greater) than the E region winds simulated by Fuller-Rowell (1985). The effects of large-scale cross-arc winds were simulated to evaluate their effect; it is shown that strong cross-arc winds are required to prevent the formation of a strong E zonal jet. The results agree qualitatively with radar observations of zonal winds by Johnson et al. (1987).

Walterscheid, R. L.↗

Overview of High Speed Close-Up Imaging in an Icing Environment

The Icing Branch and Imaging Technology Center at NASA Glenn Research Center have recently been involved in several projects where high speed close-up imaging was used to investigate water droplet impact/splash, and also ice particle impact/bounce in an icing wind tunnel. The combination of close-up and high speed imaging capabilities were required because the particles being studied were relatively small (d < 1 mm in diameter), and the impact process occurred in a very short time period (t(sub impact) << 1 sec). High speed close-up imaging was utilized to study the dynamics of droplet impact and splash in simulated Supercooled Large Droplet (SLD) icing conditions. The objective of this test was to evaluate the capability of a ultra high speed camera system to acquire quantitative information about the impact process (e.g., droplet size, velocity). Imaging data were obtained in an icing wind tunnel for spray cloud MVD > 50 m. High speed close-up imaging was also utilized to characterize the impact of ice particles on an airfoil with a thermally protected leading edge. The objective of this investigation was to determine whether ice particles tend to "stick" or "bounce" after impact. Imaging data were obtained for cases where the airfoil surface was heated and unheated. Based on the results from this test, follow on tests were conducted to investigate ice particle impact on the sensing elements of water content measurement devices. This paper will describe the use of the imaging systems to support these experimental investigations, present some representative results, and summarize what was learned about the use of these systems in an icing environment.

Miller, Dean R.↗

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↗

Consequences of the Breakout Model for Particle Acceleration in CMEs and Flares

The largest and most efficient particle accelerators in the solar system are the giant events consisting of a fast coronal mass ejection (CME) and an intense X-class solar flare. Both flares and CMEs can produce l0(exp 32) ergs or more in nonthermal particles. Two general processes are believed to be responsible: particle acceleration at the strong shock ahead of the CME, and reconnection-driven acceleration in the flare current sheet. Although shock acceleration is relatively well understood, the mechanism by which flare reconnection produces nonthermal particles is still an issue of great debate. We address the question of CME/flare particle acceleration in the context of the breakout model using 2.5D MHD simulations with adaptive mesh refinement (AMR). The AMR capability allows us to achieve ultra-high numerical resolution and, thereby, determine the detailed structure and dynamics of the flare reconnection region. Furthermore, we employ newly developed numerical analysis tools for identifying and characterizing magnetic nulls, so that we can quantify accurately the number and location of magnetic islands during reconnection. Our calculations show that flare reconnection is dominated by the formation of magnetic islands. In agreement with many other studies, we find that the number of islands scales with the effective Lundquist number. This result supports the recent work by Drake and co-workers that postulates particle acceleration by magnetic islands. On the other hand, our calculations also show that the flare reconnection region is populated by numerous shocks and other indicators of strong turbulence, which can also accelerate particles. We discuss the implications of our calculations for the flare particle acceleration mechanism and for observational tests of the models.

Antiochos, S. K.↗

Annual Research Briefs - 1996, Center for Turbulence Research

Partial contents include: (1) Simulation of jet impingement heat transfer with the k-e-V2 model; (2) A new methodology for turbulence modelers using Direct Numerical Simulation (DNS) database analysis; (3) A particle representation model for the deformation of homogeneous turbulence; (4) Modeling of inhomogeneous compressible turbulence using a two-scale statistical theory; (5) Lewis number effects on partially premixed flames; (6) Conditional moment closure of mixing and reaction in turbulent non-premixed combustion; (7) Dynamic models for LES of turbulent front propagation with a spectral method; (8) Jet noise using large eddy simulation; (9) Towards numerical simulations of trailing-edge aeroacoustics; (10) A model for the onset of vortex breakdown; (11) Compressing an elliptic vortex: transition to turbulence by tumble breakdown; (12) A new method for the adaptive control of vortex-wall interactions; (13) Symmetries in turbulent boundary layer flows; (14) Near-wall models in large eddy simulations of flow behind a backward-facing step; (15) Assessment of the subgrid-scale models at low and high Reynolds numbers; (16) Large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids; (17) Progress on LES of flow past a circular cylinder; (18) Experimental investigation of flow through an asymmetric plane diffuser; (19) Progress in the large-eddy simulation of an asymmetric plane diffuser; (20) A numerical study of self-similarity in a turbulent plane wake using large-eddy simulation; (21) On the generation of turbulent inflow conditions for boundary layer simulations, and (22) A new class of finite difference schemes.

COMPUTATIONAL FLUID DYNAMICS↗

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗

Numerical simulations of interacting disk galaxies

Galaxy-galaxy interactions have long attracted many extragalactic astronomers in various aspects. A number of computer simulations performed in the 1970s have successfully reproduced the peculiar morphologies observed in interacting disk galaxies and clarified that tidal deformation explains most of the observed global peculiarities. However, most of these simulations have used test particles in modelling the disk component. Tidal response of a self-gravitating disk remains to be further clarified. Another topic which is intensely discussed at present is the relation between galaxy-galaxy interactions and activity. Many observations suggest that interactions trigger strong starbursts and possibly active galactic nuclei (AGN). However, the detailed mechanism of triggering is not yet clear. It is vital here to understand the dynamics of interstellar gas. In order to understand various phenomena related to galaxy-galaxy interactions (mainly for disk galaxies), the author performed a series of numerical simulations on close galaxy encounters which includes both interstellar gas and self-gravitating disk components. In these simulations, the galaxy model to be perturbed (target galaxy) consists of a halo and a disk. The halo was treated as a rigid spherical gravitational field which is assumed to remain fixed during the interaction. The disk is composed of stars and gas. The stellar disk was constructed by 20000 collisionless particles of the same mass. Those particles move in the halo gravitational field, interacting with each other and with the perturber. Therefore, the self-graviy of the disk is properly taken into account. Stellar particles were initially given circular velocities with small random motions required to stabilize the disk against local axisymmetric disturbances. The gravitational field of the stellar disk was calculated by the particle-mesh scheme (e.g. Hockney and Eastwood 1981). The gaseous component was modelled by the cloud-particle scheme (e.g. Roberts and Hausman 1984). Here, the authors represent the gas as an ensemble of small spheres (i.e. clouds) and include the creation of an OB star in a cloud-cloud collision and subsequent velocity push on nearby clouds due to a supernova explosion.

Noguchi, Masafumi↗

Kinetic simulations of the Kelvin-Helmholtz instability at the magnetopause

Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the Kelvin-Helmholtz instability for a magnetopauselike configuration. The unidirectional magnetic field is essentially transverse to the plasma flow velocity, which is the most unstable case according to linear theory and models the flow dynamics in the subsolar region of the magnetopause for northward interplanetary magnetic field. We recover effects analogous to those found in MHD simulations, including a mode cascade to longer wavelengths. The boundary layer consists of coherent structure and is not well described by a diffusive process. Isolated structures on the order of the ion gyroradius are formed which can cross the boundary in either direction. We describe how the time evolution of these structures represents transport across boundary layers, and we consider the possible connection of these entities to flux transfer events and other structure seen in the low-latitude boundary layer at the Earth's magnetopause as well as to flux ropes commonly observed near the ionopause of Venus. We also discuss the relation of the hybrid calculations to previous MHD simulations and to observations.

Thomas, V. A.↗

Dynamics Of Ions In A Radio-Frequency Quadrupole Trap

Report describes computer-simulation study of motions of various numbers of ions in Paul trap. Study part of continuing effort to understand motions of trapped charged particles (atoms, ions, molecules, or dust particles). Motions characterized in terms of heating by radio-frequency fields, formation of crystallike structures in cold clouds of trapped particles, and other phenomena important in operation of radio-frequency traps in frequency standards.

Prestage, John D.↗

Observations of Deep Ionospheric F-Region Density Depletions with FPMU Instrumentation and Their Relationship with the Global Dynamics of the June 22-23, 2015 Geomagnetic Storm

The magnetic storm that commenced on June 22, 2015 was one of the largest storms in the current solar cycle. During this event, ionospheric F-region density measurements from the Floating Potential Measurement Unit (FPMU) on board the International Space Station (ISS) show dramatic depletions in the post-sunset (nighttime) local time sector at equatorial latitudes starting in the main phase of the storm and persisting on several subsequent orbits into the next day. Putting these low-latitude measurements in context with the global dynamics of the storm, we will present results from simulations and observations in our efforts to better understand the effects of this storm on the different regions of the coupled ionosphere-magnetosphere. The consequences of the magnetospheric penetration electric field and their role in the occurrence of these equatorial spread F observations will be investigated through the results of the SAMI3-RCM numerical model, a coupled ionosphere-magnetosphere model with self-consistent large-scale electrodynamics. Specifically, we will investigate the transient signatures of the interplanetary magnetic field component, Bz, and its role in driving the global convection electric field and ionospheric density redistribution. Lastly, measurements from the AMPERE Birkeland currents, DMSP drift velocities and the particle flux dropouts observed from the Magnetospheric Multiscale Mission (MMS) will be correlated with the FPMU density depletions and each other. Together these observations and simulation results will be assembled to provide each region’s context to the global dynamics and time evolution of the storm.

Coffey, Victoria↗

Observations of Deep Ionospheric F-Region Density Depletions with FPMU Instrumentation and Their Relationship with the Global Dynamics of the June 22-23, 2015 Geomagnetic Storm

The magnetic storm that commenced on June 22, 2015 was one of the largest storms in the current solar cycle. During this event, ionospheric F-region density measurements from the Floating Potential Measurement Unit (FPMU) on board the International Space Station (ISS) show dramatic depletions in the post-sunset (nighttime) local time sector at equatorial latitudes starting in the main phase of the storm and persisting on several subsequent orbits into the next day. Putting these low-latitude measurements in context with the global dynamics of the storm, we will present results from simulations and observations in our efforts to better understand the effects of this storm on the different regions of the coupled ionosphere-magnetosphere. The consequences of the magnetospheric penetration electric field and their role in the occurrence of these equatorial spread F observations will be investigated through the results of the SAMI3-RCM numerical model, a coupled ionosphere-magnetosphere model with self-consistent large-scale electrodynamics. Specifically, we will investigate the transient signatures of the interplanetary magnetic field component, Bz, and its role in driving the global convection electric field and ionospheric density redistribution. Lastly, measurements from the AMPERE Birkeland currents, DMSP drift velocities and the particle flux dropouts observed from the Magnetospheric Multiscale Mission (MMS) will be correlated with the FPMU density depletions and each other. Together these observations and simulation results will be assembled to provide each region's context to the global dynamics and time evolution of the storm.

Coffey, Victoria↗

Development of test methods for scale model simulation of aerial applications in the NASA Langley Vortex Research Facility

As part of basic research to improve aerial applications technology, methods were developed at the Langley Vortex Research Facility to simulate and measure deposition patterns of aerially-applied sprays and granular materials by means of tests with small-scale models of agricultural aircraft and dynamically-scaled test particles. Interactions between the aircraft wake and the dispersed particles are being studied with the objective of modifying wake characteristics and dispersal techniques to increase swath width, improve deposition pattern uniformity, and minimize drift. The particle scaling analysis, test methods for particle dispersal from the model aircraft, visualization of particle trajectories, and measurement and computer analysis of test deposition patterns are described. An experimental validation of the scaling analysis and test results that indicate improved control of chemical drift by use of winglets are presented to demonstrate test methods.

Jordan, F. L., Jr.↗

Formation and stability of the self-consistent one-dimensional tail current sheet

The paper investigates the formation, the structure, and the stability of self-consistent one-dimensional current sheets in which the ions carry most of the current and momentum (the occurrence of which was suggested by observations of Mitchell et al., 1990; and Sergeev et al., 1990). Results of the analysis showed that, for the case of a cold current sheet, the characteristic thickness lamba equals to about (Bz/B0) exp 4/3 c/omega(p0), where Bz is the normal field component, B0 is the asymptotic magnitude of the reversing field, and c/omega(p0)is the collisionless ion skin depth based on lobe density. A two-dimensional self-consistent dynamical simulation model is developed, which demonstrates that these idealized current sheets are unstable to kink perturbations driven by the anisotropic pressure distribution produced by the chaotic nature of the particle orbits in a field-reversal region.

Pritchett, P. L.↗

On the Interaction Between a Flow Instability and Inertial Particles in Supersonic Jet Plume-Surface Interaction

Extraterrestrial landings often require spacecraft to deploy supersonic jets directed at dusty planetary surfaces to ensure controlled descents. This interaction, termed plume-surface interaction (PSI), critically influences crater formation and the dynamics of ejected particles, which can obscure vision, damage onboard sensors, and impact nearby infrastructure with significant particle velocities. The physics of PSI remain poorly understood due to the complexities in simulating the involved multiphysics and experimentally replicating reduced atmospheric pressure conditions. This article discusses a collaborative effort to investigate PSI in reduced atmospheric pressure environments using a Mach 5.3 jet impinging on a granular bed. The results highlight observations of an azimuthal pattern in the ejecta distribution, characterized by alternating ejecta streaks of varying particle concentrations. This phenomenon suggests the presence of a low Reynolds number instability, a notable discovery given the supersonic nature of the jet.

PSI↗

Molecular Dynamics Calculations

The development of thermodynamics and statistical mechanics is very important in the history of physics, and it underlines the difficulty in dealing with systems involving many bodies, even if those bodies are identical. Macroscopic systems of atoms typically contain so many particles that it would be virtually impossible to follow the behavior of all of the particles involved. Therefore, the behavior of a complete system can only be described or predicted in statistical ways. Under a grant to the NASA Lewis Research Center, scientists at the Case Western Reserve University have been examining the use of modern computing techniques that may be able to investigate and find the behavior of complete systems that have a large number of particles by tracking each particle individually. This is the study of molecular dynamics. In contrast to Monte Carlo techniques, which incorporate uncertainty from the outset, molecular dynamics calculations are fully deterministic. Although it is still impossible to track, even on high-speed computers, each particle in a system of a trillion trillion particles, it has been found that such systems can be well simulated by calculating the trajectories of a few thousand particles. Modern computers and efficient computing strategies have been used to calculate the behavior of a few physical systems and are now being employed to study important problems such as supersonic flows in the laboratory and in space. In particular, an animated video (available in mpeg format--4.4 MB) was produced by Dr. M.J. Woo, now a National Research Council fellow at Lewis, and the G-VIS laboratory at Lewis. This video shows the behavior of supersonic shocks produced by pistons in enclosed cylinders by following exactly the behavior of thousands of particles. The major assumptions made were that the particles involved were hard spheres and that all collisions with the walls and with other particles were fully elastic. The animated video was voted one of two winning videos in a competition held at the meeting of the American Physical Society's Division of Fluid Dynamics, held in Atlanta, Georgia, in November 1994. Of great interest was the result that in every shock there were a few high-speed precursor particles racing ahead of the shock, carrying information about its impending arrival. Most recently, Dr. Woo has been applying molecular dynamics techniques to the problem of determining the drag produced by the space station truss structure as it flies through the thin residual atmosphere of low-Earth orbit. This problem is made difficult by the complex structure of the truss and by the extreme supersonic nature of the flow. A fully filled section of the truss has already been examined, and drag predictions have been made. Molecular dynamics techniques promise to make realistic drag calculations possible even for very complex partially filled truss segments flying at arbitrary angles.

Source record↗

Validation and Simulation of Ares I Scale Model Acoustic Test - 3 - Modeling and Evaluating the Effect of Rainbird Water Deluge Inclusion

The Ares I Scale Model Acoustics Test (ASMAT) is a series of live-fire tests of scaled rocket motors meant to simulate the conditions of the Ares I launch configuration. These tests have provided a well documented set of high fidelity measurements useful for validation including data taken over a range of test conditions and containing phenomena like Ignition Over-Pressure and water suppression of acoustics. Building on dry simulations of the ASMAT tests with the vehicle at 5 ft. elevation (100 ft. real vehicle elevation), wet simulations of the ASMAT test setup have been performed using the Loci/CHEM computational fluid dynamics software to explore the effect of rainbird water suppression inclusion on the launch platform deck. Two-phase water simulation has been performed using an energy and mass coupled lagrangian particle system module where liquid phase emissions are segregated into clouds of virtual particles and gas phase mass transfer is accomplished through simple Weber number controlled breakup and boiling models. Comparisons have been performed to the dry 5 ft. elevation cases, using configurations with and without launch mounts. These cases have been used to explore the interaction between rainbird spray patterns and launch mount geometry and evaluate the acoustic sound pressure level knockdown achieved through above-deck rainbird deluge inclusion. This comparison has been anchored with validation from live-fire test data which showed a reduction in rainbird effectiveness with the presence of a launch mount.

Strutzenberg, Louise L.↗

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗