Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle simulations”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Impacts of Anisotropic Porosity on Heat Transfer and Off-Gassing during Biomass Pyrolysis

The pore structure of biogenic materials imbues the ability to deliver water and nutrients through a plant from root to leaf. This anisotropic pore granularity can also play a significant role in processes such as biomass pyrolysis that are used to convert these materials into useful products like heat, fuel, and chemicals. Evolutions in modeling of biomass pyrolysis as well as imaging of pore structures allow for further insights into the concerted physics of phase change-induced off-gassing, heat transfer, and chemical reactions. In this work, we report a biomass single particle model which incorporates these physics to explore the impact of implementing anisotropic permeability and diffusivity on the conversion time and yields predicted for pyrolysis of oak and pine particles. Simulation results showed that anisotropic permeability impacts predicted conversion time more than 2 times when the Biot number is above 0.1 and pyrolysis numbers (Py 1 , Py 2 ) are less than 20. Pore structure significantly impacts predicted pyrolytic conversion time (>8 times) when the Biot number is above 1 and the pyrolysis number is below 1, i.e., the “conduction controlled” regime. Therefore, these nondimensional numbers reflect that when internal heat conduction limits pyrolysis performance, internal pyrolysis off-gassing further retards effective heat transfer rates as a closely coupled phenomenon. Overall, this study highlights physically meaningful opportunities to improve particle-scale pyrolysis modeling and experimental validation relevant to a variety of feedstock identities and preparations, guiding the future design of pyrolyzers for efficient biomass conversion.

09 BIOMASS FUELS↗
NASA NTRS2010

Depth Dose Distribution Study within a Phantom Torso after Irradiation with a Simulated Solar Particle Event at NSRL

The adequate knowledge of the radiation environment and the doses incurred during a space mission is essential for estimating an astronaut's health risk. The space radiation environment is complex and variable, and exposures inside the spacecraft and the astronaut's body are compounded by the interactions of the primary particles with the atoms of the structural materials and with the body itself Astronauts' radiation exposures are measured by means of personal dosimetry, but there remains substantial uncertainty associated with the computational extrapolation of skin dose to organ dose, which can lead to over- or underestimation of the health risk. Comparisons of models to data showed that the astronaut's Effective dose (E) can be predicted to within about a +10% accuracy using space radiation transport models for galactic cosmic rays (GCR) and trapped radiation behind shielding. However for solar particle event (SPE) with steep energy spectra and for extra-vehicular activities on the surface of the moon where only tissue shielding is present, transport models predict that there are large differences in model assumptions in projecting organ doses. Therefore experimental verification of SPE induced organ doses may be crucial for the design of lunar missions. In the research experiment "Depth dose distribution study within a phantom torso" at the NASA Space Radiation Laboratory (NSRL) at BNL, Brookhaven, USA the large 1972 SPE spectrum was simulated using seven different proton energies from 50 up to 450 MeV. A phantom torso constructed of natural bones and realistic distributions of human tissue equivalent materials, which is comparable to the torso of the MATROSHKA phantom currently on the ISS, was equipped with a comprehensive set of thermoluminescence detectors and human cells. The detectors are applied to assess the depth dose distribution and radiation transport codes (e.g. GEANT4) are used to assess the radiation field and interactions of the radiation field with the phantom torso. Lymphocyte cells are strategically embedded at selected locations at the skin and internal organs and are processed after irradiation to assess the effects of shielding on the yield of chromosome damage. The initial focus of the present experiment is to correlate biological results with physical dosimetry measurements in the phantom torso. Further on, the results of the passive dosimetry within the anthropomorphic phantoms represent the best tool to generate reliable data to benchmark computational radiation transport models in a radiation field of interest. The presentation will give first results of the physical dose distribution, the comparison with GEANT4 computer simulations based on a Voxel model of the phantom, and a comparison with the data from the chromosome aberration study.

Berger, Thomas↗
DOE OSTI2025

Overview of Geant4 Physics

Geant4 is a powerful toolkit for simulating particle interactions across a wide energy range (0 to ~TeV). It offers comprehensive physics processes for various particle types, allowing users to define custom processes using the G4VProcess interface. Major categories include electromagnetic, hadronic, decay, parameterized processes, and transportation. A critical feature is the PhysicsList, which specifies particles and processes for simulations. Users can choose between G4VUserPhysicsList for simpler environments or G4VModularPhysicsList for more complex scenarios, leveraging pre-defined or custom G4VPhysicsConstructors. Additionally, Geant4 provides reference physics lists and the G4PhysicsListFactory for well-established applications, facilitating simulation development in diverse domains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1986

A stochastic model for particle impingements on orbiting spacecraft

A general methodology for simulating particle impingements on orbiting spacecraft is developed. Major steps in the modeling process are presented as (1) modeling objective, (2) construction of the spacecraft geometrical model, (3) simulation of the particles in the space environment, (4) particle impact and subsequent events of interest, and (5) results of the simulation. A simulation of the expected meteoroid impingements on the Hubble Space Telescope and the resulting angular momentum transfers which can cause telescope pointing disturbances is given to illustrate these methods.

Howell, L. W., Jr.↗
NASA NTRS1986

A stochastic model for particle impingements on orbiting spacecraft

A general methodology for simulating particle impingements on orbiting spacecraft is developed. Major steps in the modeling process are presented as (1) modeling objective, (2) construction of the spacecraft geometrical model, (3) simulation of the particles in the space environment, (4) particle impact and subsequent events of interest, and (5) results of the simulation. A simulation of the expected meteoroid impingements on the Hubble Space Telescope and the resulting angular momentum transfers which can cause telescope pointing disturbances is given to illustrate these methods.

Howell, Leonard W., Jr.↗
NASA NTRS1990

Collisional electrostatic ion cyclotron waves as a possible source of energetic heavy ions in the magnetosphere

A new mechanism is proposed for the source of energetic heavy ions (NO/+/, O2/+/, and O/+/) found in the magnetosphere. Simulations using a multispecies particle simulation code for resistive current-driven electrostatic ion cyclotron waves show transverse and parallel bulk heating of bottomside ionospheric heavy ion populations. The dominant mechanism for the transverse bulk heating is resonant ion heating by wave-particle ion trapping. Using a linear kinetic dispersion relation for a magnetized, collisional, homogenous, and multiion plasma, it is found that collisional electrostatic ion cyclotron waves near the NO(+), O2(+), and O(+) gyrofrequencies are unstable to field-aligned currents of 50 microA/sq m for a typical bottomside ionosphere.

Providakes, Jason↗
NASA NTRS1991

Interrelation of soft and hard X-ray emissions during solar flares. II - Simulation model

Two-dimensional electrostatic particle simulations are presented which incorporate the effect of quasi-static electric fields on particle dynamics as well as effects associated with wave-particle interactions induced by the accelerated particles. The properties of the soft and hard X-ray and microwave emissions from such systems are examined. In particular, it is shown that acceleration by quasi-static electric fields and heating via wave-particle interactions produces electron distributions with a broken-power law, similar to those inferred from hard X-ray spectra. Also, heating of the ambient plasma gives rise to a region of hot plasma propagating down to the chromosphere at about the ion sound speed.

Winglee, R. M.↗
NASA NTRS2011

Process to create simulated lunar agglutinate particles

A method of creating simulated agglutinate particles by applying a heat source sufficient to partially melt a raw material is provided. The raw material is preferably any lunar soil simulant, crushed mineral, mixture of crushed minerals, or similar material, and the heat source creates localized heating of the raw material.

Gustafson, Robert J.↗
DOE OSTI2020

Simulating sticky particles: A Monte Carlo method to sample a stratification

Many problems in materials science and biology involve particles interacting with strong, short-ranged bonds that can break and form on experimental timescales. Treating such bonds as constraints can significantly speed up sampling their equilibrium distribution, and there are several methods to sample probability distributions subject to fixed constraints. We introduce a Monte Carlo method to handle the case when constraints can break and form. More generally, the method samples a probability distribution on a stratification: a collection of manifolds of different dimensions, where the lower-dimensional manifolds lie on the boundaries of the higher-dimensional manifolds. We show several applications of the method in polymer physics, self-assembly of colloids, and volume calculation in high dimensions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Discrete ordinates analysis of the forced-flight variance reduction technique in Monte Carlo neutral particle transport simulations

This paper presents mathematical formulations and methods to predict the effect of forced-flight variance reduction on Monte Carlo tally variance and calculation time. This includes deducing biasing operators that are then used to construct a history-score probability density function (HSPDF), which represents all possible Monte Carlo random walks and gives the probability of a Monte Carlo history scoring in a tally from a particular phase-space position. The history-score moment equations (HSMEs), the statistical moments of the HSPDF, are then derived to calculate the statistical behavior of the Monte Carlo tally when forced-flight variance reduction is applied. In addition, the future-time equation (FTE) is derived to predict the Monte Carlo computational time as a result of applying forced-flight variance reduction. The solutions of the HSMEs and FTE can be used to predict Monte Carlo computational cost. Furthermore, this work also describes a discrete ordinates method to solve the forced-flight HSMEs and FTE. Several 1-D and 2-D test problems verify that the derivations are performed and implemented correctly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗
DOE OSTI2022

Simulating energetic particle losses in JET plasmas with a reverse integration biasing scheme

We report that an integrated energetic particle transport model has been constructed in JET plasmas constrained by experimental fast ion loss measurements. The model incorporates a synthetic fast ion loss detector identical to JET’s thin-foil Faraday cup fast ion loss detector array. The loss model combines analyses from the TRANSP and ORBIT-kick codes with enhanced features for producing the synthetic diagnostic. Extensions to the ORBIT code framework allow a full-orbit representation within the vacuum region that can map particles directly to an installed detector geometry. Since synthetic fast ion loss detectors are plagued by weak loss statistics, a novel reverse integration biasing scheme has been implemented to boost computational efficiency. The model is validated against experimental loss measurements induced by long-lived kink modes and is found to be in good agreement. This confirms the development of a fully integrated transport/loss model which can be quantitatively verified against experiment allowing for future validation and predictive studies. The model is particularly useful for more complicated plasma scenarios that involve multiple fast ion species such as JET’s 2021 DT-campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

INTEGRATED WORKFLOW MANAGEMENT FOR PARTICLE ACCELERATOR SIMULATION

Supercomputing systems are used for a wide range of computationally demanding tasks in many fields of science and engineering. They play a key role in numerical simulation, in which mathematical models are computed in order to simulate the behavior of physical systems. Scientists and engineers that use supercomputers for numerical simulation often have their productivity limited by the need to manually organize and manage extremely large amounts of data that are often produced and consumed by the software programs run on these systems. Recognizing these limitations, Kitware Inc. (Clifton Park, NY) and SLAC National Accelerator Laboratory (Menlo Park, CA) are developing an advanced software platform that can reduce the cognitive overhead required by knowledge workers when using supercomputers for numerical simulation. Phase I of the project is complete and includes the development of new capabilities for organizing simulation project files, improvements to the user interface and overall usability, and deployment of a “middle tier” server to sit between user desktop machines and supercomputers to offload much of the data management workload. The project also developed prototype software for executing sequences of numerical simulations, and a prototype for migrating supercomputing software to cloud-based computing systems to provide a potential alternative to supercomputers with different logistical and price-to-performance tradeoffs.

Tourtellott, John↗