Search NASA⌕ Search

SEARCH · Search NASA

Results for “KNUDSEN NUMBER”

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 127 records · Page 7

Comparison of DSMC and CFD Solutions of Fire II Including Radiative Heating

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. These flows may also contain significant radiative heating. To prepare for these missions, NASA is developing the capability to simulate rarefied, ionized flows and to then calculate the resulting radiative heating to the vehicle's surface. In this study, the DSMC codes DAC and DS2V are used to obtain charge-neutral ionization solutions. NASA s direct simulation Monte Carlo code DAC is currently being updated to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced Quantum-Kinetic chemistry model, and to include electronic energy levels as an additional internal energy mode. The Fire II flight test is used in this study to assess these new capabilities. The 1634 second data point was chosen for comparisons to be made in order to include comparisons to computational fluid dynamics solutions. 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. It is shown that there can be quite a bit of variability in the vibrational temperature inferred from DSMC solutions and that, from how radiative heating is computed, the electronic temperature is much better suited for radiative calculations. To include the radiative portion of heating, the flow-field solutions are post-processed by the non-equilibrium radiation code HARA. Acceptable agreement between CFD and DSMC flow field solutions is demonstrated and the progress of the updates to DAC, along with an appropriate radiative heating solution, are discussed. In addition, future plans to generate more high fidelity radiative heat transfer solutions are discussed.

Liechty, Derek S.↗

Thermal Analysis of Low Layer Density Multilayer Insulation Test Results

Investigation of the thermal performance of low layer density multilayer insulations is important for designing long-duration space exploration missions involving the storage of cryogenic propellants. Theoretical calculations show an analytical optimal layer density, as widely reported in the literature. However, the appropriate test data by which to evaluate these calculations have been only recently obtained. As part of a recent research project, NASA procured several multilayer insulation test coupons for calorimeter testing. These coupons were configured to allow for the layer density to be varied from 0.5 to 2.6 layer/mm. The coupon testing was completed using the cylindrical Cryostat-l00 apparatus by the Cryogenics Test Laboratory at Kennedy Space Center. The results show the properties of the insulation as a function of layer density for multiple points. Overlaying these new results with data from the literature reveals a minimum layer density; however, the value is higher than predicted. Additionally, the data show that the transition region between high vacuum and no vacuum is dependent on the spacing of the reflective layers. Historically this spacing has not been taken into account as thermal performance was calculated as a function of pressure and temperature only; however the recent testing shows that the data is dependent on the Knudsen number which takes into account pressure, temperature, and layer spacing. These results aid in the understanding of the performance parameters of MLI and help to complete the body of literature on the topic.

Johnson, Wesley L.↗

Extension of the Viscous Collision Limiting Direct Simulation Monte Carlo Technique to Multiple Species

There are many flows fields that span a wide range of length scales where regions of both rarefied and continuum flow exist and neither direct simulation Monte Carlo (DSMC) nor computational fluid dynamics (CFD) provide the appropriate solution everywhere. Recently, a new viscous collision limited (VCL) DSMC technique was proposed to incorporate effects of physical diffusion into collision limiter calculations to make the low Knudsen number regime normally limited to CFD more tractable for an all-particle technique. This original work had been derived for a single species gas. The current work extends the VCL-DSMC technique to gases with multiple species. Similar derivations were performed to equate numerical and physical transport coefficients. However, a more rigorous treatment of determining the mixture viscosity is applied. In the original work, consideration was given to internal energy non-equilibrium, and this is also extended in the current work to chemical non-equilibrium.

Liechty, Derek S.↗

Significance of DSMC Computed Aerothermal Environments in the Rarefied Regime for Atmospheric Entry Material Response

During Mars atmospheric entry, the Mars Science Laboratory (MSL) was protected by a 4.5 meters diameter ablative heatshield assembled in 113 tiles. The heatshield was made of NASA's flagship ablative material, the Phenolic Impregnated Carbon Ablator (PICA). Prior work compared the traditional one-dimensional and three-dimensional material response models at different locations in the heatshield. It was observed that the flow was basically one-dimensional in the nose and flank regions, but three-dimensional flow effects were observed in the outer flank. The objective of this work is to study the effects of the aerothermal environment on the material response. We extend prior work by computing aerothermal environments using the direct simulation Monte Carlo (DSMC) code SPARTA and the CFD code Data Parallel Line Relaxation (DPLR). SPARTA is used to compute environment in the rarefied regime prior to 48.4s of entry where the Knudsen number is such that the Navier-Stokes equations can be inaccurate. Similarly to previous work, the DPLR software is used to compute the hypersonic environment for laminar then turbulent boundary layer assumptions from 48.4 s up to 100 s after Entry Interface (EI) along the MSL 08-TPS-02/01a trajectory. We observe that extending the aerothermal environments to times prior to 48.4 s modifies the thermal response of the heat shield at the surface and in-depth; however the effects on the recession are minimal. Additionally, using the assumption of a turbulent boundary layer versus a laminar one leads to higher surface and in-depth temperatures, larger recession, and a displacement of the peak heating and peak recession location.

Aerothermal↗

Comparison of CFD and DSMC Using Calibrated Transport Parameters

Hypersonic re-entry flows span a wide range of length scales where regions of both rarefied and continuum flow exist. Traditional computational fluid dynamics (CFD) techniques do not provide an accurate solution for the rarefied regions of such ‘mixed’ flow fields. Although direct simulation Monte Carlo (DSMC) can be used to accurately capture both the continuum and rarefied features of ‘mixed’ flow fields, they are computationally expensive when employed to simulate the low Knudsen number continuum regimes. Thus, a hybrid framework for seamlessly combining the two methodologies, CFD and DSMC, continues to be a topic of significant research effort. Ensuring consistency in the reaction kinetics and transport models employed within CFD and DSMC is a crucial requirement for obtaining a reliable solution from a hybrid framework for combined continuum/rarefied high speed flows. This paper represents one of the first studies to utilize the calibrated transport parameters developed to ensure consistency between CFD and DSMC solvers. The new variable soft sphere (VSS) parameters are compared to both previous “standard” variable hard sphere (VHS) parameters and also to solutions from the CFD transport properties that the new parameters were developed to reproduce.

Liechty, D. S.↗

Comparing Particle Flow Regimes in the L2K Arcjet with Martian Entry Conditions

In this study, the US3D-DUST code will be used to simulate particles in the L2K dusty-flow test environments and compare them to conditions in the shock layer of the Schiaparelli capsule during its entry into the Martian atmosphere in 2016. The particle Reynolds, Mach, and Knudsen numbers, key non-dimensional numbers that characterize the dusty-flow environments, will be compared between the L2K and Martian shock layers. This study will provide a preliminary assessment of how well-suited L2K dusty-flow conditions are for particle drag model validation and the development of particle impact surface damage correlations. The results of this study may help to guide the design of future dusty flow environments in the L2K facility

Dusty flow↗

On the Equivalence of Gas and Liquid Leak Flow Rates for Application to Propulsion and Propellant Transfer Systems

The equivalence of leak rates for gas flow and liquid flow through an arbitrary small leak size in pressurized systems has been investigated theoretically and experimentally. Measurements were taken for gaseous helium (GHe) rates and water flow issuing from the same leak configuration. The predictions for the correlation of GHe leak rates to water leak rates have been satisfactorily compared with present test data and those of other investigators over a wide range leak flow rates (many orders of magnitude) corresponding to a wide range of leak sizes (0.1 m to 0.1 mm). The results show that at extremely small leak sizes (sub-micron range), the effects of Knudsen number in the free molecular regime (rarefied gases) become important. By a proper scaling procedure, the measured flow rate for water could be applied to equivalent propellant leaks in flight operations. In the present tests the micro-channel range of the leak size is satisfactorily simulated by repeatable precision set needle position valve. Test results show that in general the theoretical model predicts the data reasonably well for GHe flow rates in the range of 0.001 to 1 sccs. Outside this range of GHe flow rates some adjustments (improvements) to the analytical model are required along with further testing.

Max Kandula↗

The Effect of Drag Model on Heatshield Recession due to Particle Impacts for Martian Spacecraft

A spacecraft entering the Martian atmosphere during a dust storm may experience recession to the heatshield due to dust particle impacts. Aerodynamic drag is the primary force that determines the trajectory of the dust particles through the shock layer. This paper examines the effect of particle drag model on the heatshield recession. Three particle drag models are assessed including two that are intended to be applicable over a wide range of particle flow conditions. Particle trajectories are computed in conditions measured during the 2007 major global dust storm. It was found that accounting for Knudsen number and compressibility effects made a large difference in the estimated particle impact velocity. The two drag models that were valid for transitional, compressible particle flow environments predicted only slightly-different amounts of heatshield recession due to dust particle impacts. A brief description of the effects of non-spherical particles on drag coefficients is provided.

Heatshield↗

The Effect of Drag Model on Heatshield Recession due to Particle Impacts for Martian Spacecraft

A spacecraft entering the Martian atmosphere during a dust storm may experience recession to the heatshield due to dust particle impacts. Aerodynamic drag is the primary force that determines the trajectory of the dust particles through the shock layer. This paper examines the effect of particle drag model on the heatshield recession. Three particle drag models are assessed including two that are intended to be applicable over a wide range of particle flow conditions. Particle trajectories are computed in conditions measured during the 2007 major global dust storm. It was found that accounting for Knudsen number and compressibility effects made a large difference in the estimated particle impact velocity. The two drag models that were valid for transitional, compressible particle flow environments predicted only slightly-different amounts of heatshield recession due to dust particle impacts. A brief description of the effects of non-spherical particles on drag coefficients is provided.

Heatshield↗

Adaptive simulations enable computational design of electron beam processing of nanomaterials with supersonic micro-jet precursor

Focused Electron Beam Induced Processing (FEBIP) is a powerful tool for the “direct-write” of nanomaterials with the possibility of atomistic control on suspended 2D material substrates. FEBIP capabilities have been significantly expanded by using a localized jet-based delivery of precursors, and especially when a thermally energized supersonic micro-jet enhances the delivery of mass flux along with controlling the far-from-equilibrium thermodynamic state of adsorbed adatoms. The possibilities of growing nanomaterials with “dialed-in” composition with ultra-high growth rates and aspect ratios up to 100:1 have been demonstrated using the supersonic micro-jet-FEBIP. Bringing this scientific discovery to the level of maturity required for practical applications in additive nanomanufacturing requires simulation tools that are capable of capturing the complex flow physics of micro-jet-substrate interactions that bridge a wide range of flow regimes from the high-density gas micro-jet expanding into the vacuum environment of FEBIP. To address this significant computational challenge, a new approach has been developed and described in this work for a multiscale adaptive DSMC (Direct Simulation Monte Carlo) algorithm to predict the micro-jet gas dynamics in an FEBIP environment, spanning the full range of flow regimes from low Knudsen (Kn) number O(0.01) continuum flow to high Kn of O(10) for the molecular flow within a unified computational framework. Here, the fundamental principles of the adaptive DSMC algorithm are described, its viability as a DSMC technique is demonstrated, and computational improvements for the benchmark cases are discussed in the context of advancing 3D nanofabrication with FEBIP. Ultimately, combining the first principle simulations via adaptive DSMC with the complementary experimental data will enable the creation of the powerful CAD tools for in silico design and optimal operation of micro-jet-FEBIP.

36 MATERIALS SCIENCE↗

One-dimensional outgassing problem

A predominantly numerical study, using the direct-simulation Monte Carlo method, of the flow from outgassing surfaces with plane, cylindrical, or spherical geometry is presented. The numerical results for the plane flow case are compared with existing theoretical and experimental results and also with a new formulation of the jump conditions across the Knudsen layer between the surface and the continuum outflow. This outflow is sonic for the free expansion case but subsonic in the presence of back-pressure effects. The self-scattered return flux is 19% of the outgassed flux in the sonic outflow case and increases at lower Mach numbers. Thermal accomodation and reflection of this return flux produces an outgassing cooling effect. The cylindrical and spherical geometries introduce an additional Knudsen number defined by the ratio of the mean free path in the effusing gas to the radius. The plane flow results apply when this Knudsen number is less than 0.001 and, when it is of order unity or above, the return flux is inversely proportional to it.

Tran Cong, T.↗

Modeling Ablation of Fibrous Materials from Bulk to Knudsen Regime

Material-environment interactions are analyzed at microscopic scale to explain the lower than expected density observed by post-flight analysis of the char layer on the Stardust shield. Mass transfer, ablation (oxidation), and surface recession of fibrous material is simulated in 3D using a Monte-Carlo simulation tool. Ablation is found to occur either at the surface or in volume depending on Knudsen and Thiele number values. This study supports the idea of volume ablation followed by possible carbon fiber spallation that may explain post-flight analyses.

Lachaud, Jean↗

Heat Transfer from Cylinders in Transition From Slip Flow to Free-Molecule Flow

Over 600 measured heat-transfer coefficients in the transition from slip to free-molecule flow have been correlated by using the Nusselt number Nu as a function of the Knudsen Kn and Reynolds Re (or Mach M) numbers. The experimental range for these heat-transfer data from transverse cylinders in air corresponds to the following dimensionless groups: M, 0.10 to 0.90; Re, 0.03 to 11.5; Kn, 0.10 to 5.0. The total air temperature T(sub t) was maintained constant at 80 F, but wire temperature was Varied from 150 to 580 F. At Kn=0.10, Nu extrapolates smoothly into slip-flow empirical curves that show Nu as a function of Re and M or Kn. The correlation gradually changes from the square root of Re(sub t) dependence characteristic of continuum flow to first-power Re dependence as Kn increases (decreasing Re). At the experimental limit Kn ft 5.0, the Nu data correlate with a mean fractional error of 413 percent by the prediction of free-molecule-flow theory. In comparing experimental results with theory, an accommodation coefficient of 0.57+/-0.07 was inferred from the heat-transfer data, which were obtained with etched tungsten wire in air. The wire recovery temperature T(sub e) was measured and compared with existing data and theory in terms of a ratio eta(equivalent to T(sub e)/T(sub t). The results can be divided into three groups by Kn criteria: For Kn less than 2.01, eta is independent of Kn, and eta decreases from 1.0 to 0.97 as M increases from 0 to 0.90; for 2.0 less than Kn less than 5.0, eta is a function of both Kn and M in this transition region to fully developed free-molecule flow; and for Kn greater than 5.0, eta predicted by free-molecule-flow theory is observed and increases from 1.0 to 1.08 as M increases from 0 to 0.90, again independent of Kn. Therefore, these T(sub e) data provide a guide to the boundary of fully developed free-molecule flow, which is.inferred from this research to exist for Kn greater than 5.0. This boundary criterion is substantiated by other published data on T(sub e) at supersonic speeds.

Cybulski, Ronald J.↗

Generalized hydrodynamic transport in lattice-gas automata

The generalized hydrodynamics of two-dimensional lattice-gas automata is solved analytically in the linearized Boltzmann approximation. The dependence of the transport coefficients (kinematic viscosity, bulk viscosity, and sound speed) upon wave number k is obtained analytically. Anisotropy of these coefficients due to the lattice symmetry is studied for the entire range of wave number, k. Boundary effects due to a finite mean free path (Knudsen layer) are analyzed, and accurate comparisons are made with lattice-gas simulations.

Luo, Li-Shi↗

Investigation of hypersonic rarefied flow on a spherical nose of the AOTV

The Navier-Stokes (NS) equations were integrated numerically for investigating the flow characteristics on the forepart of the spherical nose of a space vehicle such as the AOTV or AFE by a modified Accelerated Successive Replacement (ASR) scheme under hypersonic rarefied conditions. Technical feasibility of the mathematical approach was demonstrated by computing the flowfield on a spherical nose under conditions that the AFE encounters at times t = 15 and 20 seconds after its reentry into the atmosphere. Local similar solutions for the merged layer flow along the stagnation line of the sphere were developed. These are correct to the same degree of accuracy as the NS equations. These solutions provided stagnation line boundary conditions for the domain of integration on the spherical noise. Also, a parametric study of the stagnation line solution was made with a view to understand the flow characteristics in tunnels with different ambient fluids. Analytical expressions for surface slip temperature, jump conditions, and concentration level in the presence of the real gas effects at the top of the Knudsen layer were derived and used to calculate the stagnation line flowfield with nonequilibrium dissociation and ionization. A number of graphics were drawn to illustrate the basic physics of the flowfields. The present analysis can be extended to include real gas effects and to bodies of arbitrary shapes. It can further provide boundary conditions for integrating the NS equations in the near wake region.

Jain, Amolak C.↗

Interactions of Water Vapor with Oxides at Elevated Temperatures

Many volatile metal hydroxides form by reaction of the corresponding metal oxide with water vapor. These reactions are important in a number of high temperature corrosion processes. Experimental methods for studying the thermodynamics of metal hydroxides include: gas leak Knudsen cell mass spectrometry, free jet sampling mass spectrometry, transpiration and hydrogen-oxygen flame studies. The available experimental information is reviewed and the most stable metal hydroxide species are correlated with position in the periodic table. Current studies in our laboratory on the Si-O-H system are discussed.

Jacobson, Nathan↗

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗

Photophoretic Propulsion Enabling Mesosphere Exploration NIAC Phase I Final Report

This Phase I report presents a comprehensive study on photophoretic flyers—innovative, ultralight, solar-powered vehicles that harness photophoretic forces generated via Knudsen pumping to achieve sustained flight in the mesosphere (50–80 km altitude). By integrating advanced materials such as nanocardboard— characterized by its extremely low areal density (~1 g/m²) and high bending stiffness—with ultrathin light-absorbing coatings, the project converts incident solar radiation directly into a directed thrust. Extensive experimental investigations, coupled with high-fidelity computational fluid dynamics (CFD) simulations using ANSYS Fluent, validate the concept across various three-dimensional geometries, including spherical, conical, and rocket-shaped configurations. These simulations bridge the gap between free-molecular and continuum flow regimes, demonstrating that optimized designs can generate lift forces sufficient to support kilogram-scale payloads even in low-pressure environments. At the heart of this innovation is the use of Knudsen pumping, where temperature gradients across porous surfaces induce directional gas flow, creating a modest overpressure that provides lift. The report introduces an analytical framework that interpolates between the well-known low-Reynolds number drag regime and the high-Reynolds number momentum theory. This model accurately predicts lift based on design parameters such as microchannel dimensions, porous wall geometry, areal density, and nozzle exit area. For instance, simulations indicate that 10-meter-scale structures with carefully engineered porous walls can achieve the necessary pressure differential to support scientifically significant payloads (~1 kg). The study also explores a hybrid propulsion approach that combines solar buoyancy with photophoretic lift. Initially, solar heating creates a buoyant force that elevates the flyer to mesospheric altitudes. Once in the optimal pressure range, the photophoretic mechanism—powered by Knudsen pumping—takes over as the primary source of lift, ensuring stable, long-duration flight. This dual-mode operation not only facilitates the deployment of photophoretic flyers but also broadens the potential applications for mesospheric exploration. In addition to propulsion, the report investigates the integration of photophoretic thrusters for trajectory control of existing research balloons in the upper stratosphere. By dynamically adjusting the nozzle orientation and controlling the flow-through velocity, these thrusters provide precise maneuverability, enabling the flyers to counteract atmospheric disturbances and adjust their flight paths in real time. For example, a photophoretic thruster approximately 7.5 by 7.5 meters in size could be unfolded below a payload gondola of a 60 million-cubic-foot zero-pressure balloon. Such a thruster can provide horizontal speed control of approximately 1 m/s using only sunlight and no moving parts (except those needed to track the Sun and control the jet direction). Importantly, photophoretic thrusters operate more efficiently at higher altitudes, which is complementary to known trajectory control techniques, such as propellers and tethered wings, which are more effective at lower altitudes. Finally, the report identifies three scientific research thrusts where mesospheric aircraft technology can have a profound impact: atmospheric tides, characterization of gravity waves, and investigation of mesospheric instabilities. Overall, the findings of this Phase I project represent a significant advancement in photophoretic propulsion technology. By demonstrating that large-scale, ultralight structures can be powered solely by solar radiation—via carefully engineered Knudsen pumping—this work lays a robust foundation for scalable, near-space flight architectures. Future refinements in material fabrication, structural optimization, and integrated trajectory control are expected to further enhance performance, paving the way for operational demonstrations that could revolutionize atmospheric science, remote sensing, and communication networks.

Knudsen Pump↗