Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct simulation Monte Carlo”

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 433 records · Page 24

Simulation of Low-density Nozzle Plumes in Non-zero Ambient Pressures

The direct simulation Monte-Carlo (DSMC) method was applied to the analysis of low-density nitrogen plumes exhausting from a small converging-diverging nozzle into finite ambient pressures. Two cases were considered that simulated actual test conditions in a vacuum facility. The numerical simulations readily captured the complicated flow structure of the overexpanded plumes adjusting to the finite ambient pressures, including Mach disks and barrel shaped shocks. The numerical simulations compared well to experimental data of Rothe.

Chung, Chan-Hong↗

Meteors: A Delivery Mechanism of Organic Matter to The Early Earth

All potential exogenous pre-biotic matter arrived to Earth by ways of our atmosphere, where much material was ablated during a luminous phase called 1. meteors" in rarefied flows of high (up to 270) Mach number. The recent Leonid showers offered a first glimpse into the elusive physical conditions of the ablation process and atmospheric chemistry associated with high-speed meteors. Molecular emissions were detected that trace a meteor's brilliant light to a 4,300 K warm wake rather than to the meteor's head. A new theoretical approach using the direct simulation by Monte Carlo technique identified the source-region and demonstrated that the ablation process is critical in the heating of the meteor's wake. In the head of the meteor, organic carbon appears to survive flash heating and rapid cooling. The temperatures in the wake of the meteor are just right for dissociation of CO and the formation of more complex organic compounds. The resulting materials could account for the bulk of pre-biotic organic carbon on the early Earth at the time of the origin of life.

Jenniskens, Peter↗

Testing of Error-Correcting Sparse Permutation Channel Codes

A computer program performs Monte Carlo direct numerical simulations for testing sparse permutation channel codes, which offer strong error-correction capabilities at high code rates and are considered especially suitable for storage of digital data in holographic and volume memories. A word in a code of this type is characterized by, among other things, a sparseness parameter (M) and a fixed number (K) of 1 or "on" bits in a channel block length of N.

Shcheglov, Kirill, V.↗

Monte Carlo Simulation to Estimate Likelihood of Direct Lightning Strikes

A software tool has been designed to quantify the lightning exposure at launch sites of the stack at the pads under different configurations. In order to predict lightning strikes to generic structures, this model uses leaders whose origins (in the x-y plane) are obtained from a 2D random, normal distribution.

Mata, Carlos↗

Applications of Low Density Flow Techniques and Catalytic Recombination at the Johnson Space Center

The talk presents a brief background on defInitions of catalysis and effects associated with chemically nonequilibrium and low-density flows of aerospace interest. Applications of catalytic recombination on surfaces in dissociated flow are given, including aero heating on reentry spacecraft thermal protection surfaces and reflection of plume flow on pressure distributions associated with the space station. Examples include aero heating predictions for the X-38 test vehicle, the inlet of a proposed gas-sampling probe used in high enthalpy test facilities, and a parabolic body at angle of attack. The effect of accommodation coefficients on thruster induced pressure distributions is also included. Examples of tools used include simple aero heating formulas based on boundary layer solutions, an engineering approximation that uses axisymmetric viscous shock layer flow to simulate full three dimensional flow, full computational fluid dynamics, and direct simulation Monte-Carlo calculations. Methods of determining catalytic recombination rates in arc jet flow are discus ed. An area of catalysis not fully understood is the formation of single-wall carbon nanotubes (SWNT) with gas phase or nano-size metal particles. The Johnson Space Center is making SWNTs using both a laser ablation technique and an electric arc vaporization technique.

Scott, Carl D.↗

Towards Plume Impingement Modeling in Space Environments

After 30 years human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to act as a home base for Lunar exploration. A large space station must be assembled while in orbit, using a “piecemeal” approach. In this context, it means that the different modules will arrive at different times and attach to what is already in service. The ISS provides an excellent example of this approach, and the proposed Lunar Gateway will undergo a similar assembly process. This assembly is achieved via “docking” maneuvers between modules, which are made possible by sequential firings of the onboard reaction control system (RCS) thrusters. They work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking approach. The issue is that the gas from these thrusters' forms flow structures described as “plumes” and can impinge onto the outer surfaces of the space station, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion. All mechanisms that can damage the space station and must be avoided. Both permanent and visiting modules will have these RCS thruster exhaust impingement problems. Accurately and efficiently modeling these plume is an involved multi-physics calculation but also an important tool when designing the control algorithms of approaching modules. This poster presents progress towards this simulation on two fronts. First is the verification of OpenFOAM for rarefied plume impingement calculations by direct comparison to published DAC cases. Second is the estimation of plume impingement strikes over the time scale of an entire docking event. This is done by using a simple plume source flow model and a prescribed motion visualizer—coded in Python. Together these tools help push NASA's capabilities for simulating these plume impingement effects.

Rarefied Flows↗

The Venus Atmospheric Sample Return (VATMOS-SR) Mission Concept

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small space- craft mission concept that would return several gas samples from the upper atmosphere of Venus (< 110 km) to Earth for scientific analysis. This could be the first sample return mission for an extraterrestrial 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 ac- quire 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 (e.g. trajectory shown in Fig. 1).

direct simulation Monte Carlo↗

The Venus Atmospheric Sample Return (VATMOS-SR) Mission Concept

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small space- craft mission concept that would return several gas samples from the upper atmosphere of Venus (< 110 km) to Earth for scientific analysis. This could be the first sample return mission for an extraterrestrial 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 ac- quire 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.

direct simulation Monte Carlo↗

Realistic Covariance Generation for the GPM Spacecraft

We present several different methods for generating realistic predictive covariance, including Monte-Carlo simulations and more direct linear methods which require the addition of process noise. The Monte-Carlo simulation starts with an epoch uncertainty sample basis and propagates each trial to a time of interest in the future. The variance-covariance of the state elements as well as other higher order sample statistics can be readily computed from the propagated sample. While this method preserves the nonlinear effects on the propagated uncertainty, it is computationally intensive as a statistically significant sample size must be considered in the propagation process. Moreover, if the epoch covariance is optimistic, this can result in an underestimation of the prediction error. Another method is to propagate a state sensitivity matrix simultaneously with the satellite state, which allows the epoch covariances to be propagated forward in a linear fashion. This method does not preserve the non-linearity of the satellite state uncertainty but is much less computationally intensive. The propagated state covariance is the scaled to represent the realistic level of GPM state uncertainty via a "e-tuning process." The tuning process generates an inflation factor based on the observed error statistics of the predictive satellite trajectories when compared to the definitive ones. Difference tuning strategies are considered and compared via Goodness-of-Fit method testing for the Gaussian properties of the scaled covariance.

prediction error↗

Enabling the broader adoption of fusion simulation on complex geometry

This project addressed a key barrier to advanced fusion and nuclear simulation: the difficulty of performing high-fidelity Monte Carlo neutronics directly on complex, real-world CAD geometry. Traditional workflows require engineers to rebuild CAD models as simplified constructive solid geometry, a time-consuming and error-prone process that limits design iteration and broader adoption of simulation tools. The goal of this Phase I SBIR was to make CAD-based neutronics practical, accessible, and robust for industrial and research users. During the project, Coreform significantly enhanced the Direct Accelerated Geometry Monte Carlo (DAGMC) workflow and fully integrated it into Coreform Cubit as a first-class capability. Major achievements include optimized material assignment and surface meshing workflows, substantial performance improvements to geometry imprinting and preparation, native export of DAGMC models, and new visualization tools to support OpenMC source definition and lost-particle debugging. Coreform also expanded Cubit’s capabilities as a full OpenMC preprocessor, including the ability to convert OpenMC constructive solid geometry models back into CAD for visualization, multiphysics coupling, and debugging. In collaboration with Argonne National Laboratory, the project delivered comprehensive new DAGMC documentation and training materials, transforming DAGMC from a research-oriented tool into a production-ready workflow. Results were disseminated through tutorials, conference training, and multiple well-attended webinars demonstrating integrated CAD-based neutronics and multiphysics workflows. Overall, this project demonstrated that high-fidelity Monte Carlo simulations can be performed directly on complex CAD geometry, reducing setup time, improving usability, and enabling faster, more informed design decisions for fusion and nuclear energy systems.

42 ENGINEERING↗

An Application of the Direct Coulomb Electron Pair Production Process to the Energy Measurement of the "VH-Group" in the "Knee" Region of the "All-Particle" Energy Spectrum

The "all-particle" cosmic ray energy spectrum appears to be exhibiting a significant change in the spectral index just above approximately 3000 TeV. This could indicate (1) a change in the propagation of the cosmic rays in the galactic medium, and/or (2) the upper limit of the supernova shock wave acceleration mechanism, and/or (3) a new source of high-energy cosmic rays. Air shower and JACEE data indicate the spectral change is associated with a composition change to a heavier element mixture whereas DICE does not indicate this. A detector concept will be presented that utilizes the energy dependence of the production of direct Coulomb electron-positron pairs by energetic heavy ions. Monte Carlo simulations of a direct electron pair detector consisting of Pb target foils interleaved with planes of 1-mm square scintillating optical fibers will be discussed. The goal is to design a large area, non-saturating instrument to measure the energy spectrum of the individual cosmic ray elements in the "VH-group" for energies greater than 10 TeV/nucleon.

Derrickson, J. H.↗

Simulations of particle acceleration in parallel shocks: Direct comparison between Monte Carlo and one-dimensional hybrid codes

We have made a direct comparison between two different computer simulations of a plane, parallel, collisionless shock including particle acceleration to energies typical of those of diffuse ions observed at the earth bow shock. Despite the fact that the one-dimensional hybrid and Monte Carlo techniques employ entirely different algorithms, they give surprisingly close agreement in the overall shapes of the complete distribution functions for protons as well as heavier ions. Both methods show that energetic ions emerge smoothly from the background thermal plasma with approximately the same relative injection rate and that the fraction of the incoming plasma's energy flux that is converted into downstream enthalpy flux of the accelerated population (i.e., the acceleration efficiency) is similar in the two cases. The fraction of the downstream proton distribution made up of superthermal particles is quite large, with at least 10% of the energy flux going into protons with energies above 10 keV. In addition, an upstream precursor, produced by backstreaming energetic particles, is present in both shocks, although the Monte Carlo precursor is considerably longer than that produced in the hybrid shock. These results offer convincing evidence that, at least in these ways, the two simulations are consistent in their description of parallel shock structure and particle acceleration, and they lay the groundwork for development of shock models employing a combination of both methods.

Ellison, Donald C.↗

Probalistic Finite Elements (PFEM) structural dynamics and fracture mechanics

The purpose of this work is to develop computationally efficient methodologies for assessing the effects of randomness in loads, material properties, and other aspects of a problem by a finite element analysis. The resulting group of methods is called probabilistic finite elements (PFEM). The overall objective of this work is to develop methodologies whereby the lifetime of a component can be predicted, accounting for the variability in the material and geometry of the component, the loads, and other aspects of the environment; and the range of response expected in a particular scenario can be presented to the analyst in addition to the response itself. Emphasis has been placed on methods which are not statistical in character; that is, they do not involve Monte Carlo simulations. The reason for this choice of direction is that Monte Carlo simulations of complex nonlinear response require a tremendous amount of computation. The focus of efforts so far has been on nonlinear structural dynamics. However, in the continuation of this project, emphasis will be shifted to probabilistic fracture mechanics so that the effect of randomness in crack geometry and material properties can be studied interactively with the effect of random load and environment.

Liu, Wing-Kam↗

A physics-informed deep learning model of the hot tail runaway electron seed

A challenging aspect of the description of a tokamak disruption is evaluating the hot tail runaway electron seed that emerges during the thermal quench. This problem is made challenging due to the requirement of describing a strongly non-thermal electron distribution, together with the need to incorporate a diverse range of multiphysics processes, including magnetohydrodynamic instabilities, impurity transport, and radiative losses. Here this work develops a physics-informed neural network (PINN) tailored to the solution of the hot tail seed during an axisymmetric thermal quench. Here, a PINN is developed to identify solutions to the adjoint relativistic Fokker–Planck equation in the presence of a rapid quench of the plasma's thermal energy. It is shown that the PINN is able to accurately predict the hot tail seed across a range of parameters, including the thermal quench timescale, initial plasma temperature, and local current density, in the absence of experimental or simulation data. The hot tail PINN is verified by comparison with direct Monte Carlo simulations, with excellent agreement found across a broad range of thermal quench conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Single molecule tracking of bacterial cell surface cytochromes reveals dynamics that impact long-distance electron transport

Using a series of multiheme cytochromes, the metal-reducing bacterium Shewanella oneidensis MR-1 can perform extracellular electron transfer (EET) to respire redox-active surfaces, including minerals and electrodes outside the cell. While the role of multiheme cytochromes in transporting electrons across the cell wall is well established, these cytochromes were also recently found to facilitate long-distance (micrometer-scale) redox conduction along outer membranes and across multiple cells bridging electrodes. Recent studies proposed that long-distance conduction arises from the interplay of electron hopping and cytochrome diffusion, which allows collisions and electron exchange between cytochromes along membranes. However, the diffusive dynamics of the multiheme cytochromes have never been observed or quantified in vivo, making it difficult to assess their hypothesized contribution to the collision-exchange mechanism. Here, we use quantum dot labeling, total internal reflection fluorescence microscopy, and single-particle tracking to quantify the lateral diffusive dynamics of the outer membrane-associated decaheme cytochromes MtrC and OmcA, two key components of EET in S. oneidensis. We observe confined diffusion behavior for both quantum dot-labeled MtrC and OmcA along cell surfaces (diffusion coefficients D MtrC = 0.0192 ± 0.0018 µm 2 /s, D OmcA = 0.0125 ± 0.0024 µm 2 /s) and the membrane extensions thought to function as bacterial nanowires. We find that these dynamics can trace a path for electron transport via overlap of cytochrome trajectories, consistent with the long-distance conduction mechanism. The measured dynamics inform kinetic Monte Carlo simulations that combine direct electron hopping and redox molecule diffusion, revealing significant electron transport rates along cells and membrane nanowires.

59 BASIC BIOLOGICAL SCIENCES↗

Understanding Reentrance in Frustrated Magnets: The Case of the Er 2 Sn 2 O 7 Pyrochlore

Reentrance, the return of a system from an ordered phase to a previously encountered less-ordered one as a controlled parameter is continuously varied, is a recurring theme found in disparate physical systems, yet its microscopic cause is often not investigated thoroughly. Here, through detailed characterization and theoretical modeling, we uncover the microscopic mechanism behind reentrance in the strongly frustrated pyrochlore antiferromagnet Er 2 Sn 2 O 7 . We use single crystal heat capacity measurements to expose that Er 2 Sn 2 O 7 exhibits multiple instances of reentrance in its magnetic field B vs temperature T phase diagram for magnetic fields along three cubic high symmetry directions. Through classical Monte Carlo simulations, mean field theory, and classical linear spin-wave expansions, we argue that the origins of the multiple occurrences of reentrance observed in Er 2 Sn 2 O 7 are linked to soft modes. These soft modes arise from phase competition and enhance thermal fluctuations that entropically stabilize a specific ordered phase, resulting in an increased transition temperature for certain field values and thus the reentrant behavior. Our work represents a detailed examination into the mechanisms responsible for reentrance in a frustrated magnet and may serve as a template for the interpretation of reentrant phenomena in other physical systems.

36 MATERIALS SCIENCE↗

Lunar Outpost Life Support Architecture Study Based on a High Mobility Exploration Scenario

As scenarios for lunar surface exploration and habitation continue to evolve within NASA s Constellation program, so must studies of optimal life support system architectures and technologies. This paper presents results of a life support architecture study based on a 2009 NASA scenario known as Scenario 12. Scenario 12 represents a consolidation of ideas from earlier NASA scenarios and includes an outpost near the Lunar South Pole comprised of three larger fixed surface elements and four attached pressurized rovers. The scenario places a high emphasis on surface mobility, with planning assuming that all four crewmembers spend roughly 50% of the time away from the outpost on 3-14 day excursions in two of the pressurized rovers. Some of the larger elements can also be mobilized for longer duration excursions. This emphasis on mobility poses a significant challenge for a regenerative life support system in terms of cost-effective waste collection and resource recovery across multiple elements, including rovers with very constrained infrastructure resources. The current study considers pressurized rovers as part of a distributed outpost life support architecture in both stand-alone and integrated configurations. A range of architectures are examined reflecting different levels of closure and distributed functionality. Different lander propellant scavenging options are also considered involving either initial conversion of residual oxygen and hydrogen propellants to water or initial direct oxygen scavenging. Monte Carlo simulations are used to assess the sensitivity of results to volatile high-impact mission variables, including the quantity of residual lander propellants available for scavenging, the fraction of crew time away from the outpost on excursions, total extravehicular activity hours, and habitat leakage. Architectures are evaluated by estimating surpluses or deficits of water and oxygen per 180-day mission and differences in fixed and 10-year-total equivalent system mass (ESM) relative to a reference case. Results are presented based on current assumptions for Scenario 12 and based on Monte Carlo simulations with assumed probability distributions for the high-impact mission variables. The calculated probability of no water or oxygen resupply from Monte Carlo simulations provides a quantitative measure of system robustness that can be used for cost/benefit analyses to identify leading architecture candidates. Areas of technology improvement that are likely to have a significant impact are also suggested.

Lange, Kevin E.↗