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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Simulation of the current-driven electrostatic ion cyclotron instability

Results are presented of a self-consistent particle simulation of the growth and saturation of electrostatic ion cyclotron (EIC) waves along auroral field lines. The driving mechanism for the waves is an initial drifting Maxwellian distribution for the electrons. The magnetic field is taken to lie in the simulation plane, so that modes such as the ion acoustic instability and the EIC instability can be excited, thus making it possible to study the behavior of EIC waves in the presence of ion acoustic turbulence. A dc electric field, which tends to keep the current constant along the field lines, is included. It is shown that the EIC instability grows to a moderate level before saturating due to plateau formation. The perpendicular ion temperature increases by typically 40%, with an increase arising from bulk heating only. The anomalous resistivity, associated with the EIC instability is determined and compared with previous calculations and with that due to the ion acoustic instability. Implications of the results for observation of ion cyclotron waves on auroral field lines are discussed.

Pritchett, P. L.↗
NASA NTRS1982

Nonlinear longitudinal resonance interaction of energetic charged particles and VLF waves in the magnetosphere

The longitudinal resonance of waves and energetic electrons in the Earth's magnetosphere, and the possible role this resonance may play in generating various magnetospheric phenomena are studied. The derivation of time-averaged nonlinear equations of motion for energetic particles longitudinally resonant with a whistler mode wave propagating with nonzero wave normal is considered. It is shown that the wave magnetic forces can be neglected at lower particle pitch angles, while they become equal to or larger than the wave electric forces for alpha 20 deg. The time-averaged equations of motion were used in test particle simulation which were done for a wide range of wave amplitudes, wave normals, particle pitch angles, particle parallel velocities, and in an inhomogeneous medium such as the magnetosphere. It was found that there are two classes of particles, trapped and untrapped, and that the scattering and energy exchange for those two groups exhibit significantly different behavior.

Tkalcevic, S.↗
NASA NTRS1997

Current Collection in a Magnetic Field

It is found that the upper-bound limit for current collection in the case of strong magnetic field from the current is close to that given by the Parker-Murphy formula. This conclusion is consistent with the results obtained in laboratory experiments. This limit weakly depends on the shape of the wire. The adiabatic limit in this case will be easily surpassed due to strong magnetic field gradients near the separatrix. The calculations can be done using the kinetic equation in the drift approximation. Analytical results are obtained for the region where the Earth's magnetic field is dominant. The current collection can be calculated (neglecting scattering) using a particle simulation code. Dr. Singh has agreed to collaborate, allowing the use of his particle code. The code can be adapted for the case when the current magnetic field is strong. The needed dm for these modifications is 3-4 months. The analytical description and essential part of the program is prepared for the calculation of the current in the region where the adiabatic description can be used. This was completed with the collaboration of Drs. Khazanov and Liemohn. A scheme of measuring the end body position is also proposed. The scheme was discussed in the laboratory (with Dr. Stone) and it was concluded that it can be proposed for engineering analysis.

Krivorutsky, E. N.↗
DOE OSTI2021

Investigation of theoretical scaling laws using large eddy simulations for airborne spreading of viral contagion from sneezing and coughing

Using a set of large eddy point-particle simulations, we explore the fluid dynamics of an ejected puff resulting from a cough/sneeze. The ejection contains over 61,000 potentially virus-laden droplets at an injection Reynolds number of about 46,000, comparable to an actual cough/ sneeze. We observe that global puff properties, such as centroid, puff volume, momentum, and buoyancy vary little across realizations. Other properties, such as maximum extent, shape, and edge velocity of the puff, may exhibit substantial variation. In many realizations, a portion of the puff splits off and advances along a random direction, while keeping airborne droplet nuclei afloat. This peeled-off portion provides a mechanism for virus-laden droplets to travel over large distances in a short amount of time. We also observe that the vast majority of droplets remain suspended within the puff after all liquid has evaporated. The main objectives of the study are to (i) evaluate assumptions of Balachandar’s et al. theory [Int. J. Multiphase Flow 132, 103439 (2020)], which include buoyancy effects, shape of the puff, and droplet evaporation rate, (ii) obtain values of closure parameters, which include location and time of the virtual origin, and puff entrainment and drag coefficients, and (iii) evaluate the accuracy of the theory in predicting the shape, size, and location of the puff, as well as droplet number density long after ejection. The theory adequately predicts global puff properties including size, velocity, and distance traveled, the largest size of droplets that exit the puff due to settling, and the droplet size distribution within the puff long after ejection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1983

Simulations of beam excited minor species gyroharmonics in the Porcupine experiment

An active experiment for the study of wave particle interaction in the presence of a perpendicular beam type distribution is considered. The experiment involves the injection of a beam perpendicular to the magnetic field. Beam parameters can then be varied while local plasma parameters and associated wave phenomena are measured. Such an experiment was conducted as part of the Porcupine project. The present investigation has the objective to explain theoretically and by means of plasma particle simulations, some of the wave observations obtained. In the experiment, a main rocket payload and four subpayloads were employed. The subpayloads were ejected radially in different directions at an altitude of 240 km. One of the spinning subpayloads contained a plasma gun which emitted a beam of 200 eV xenon ions toward the main payload. The data recorded on the main payload revealed the existence of electrostatic waves with peaks at hydrogen gyroharmonics. A series of computer simulations approximating the xenon experiment was performed. The interpretation of the data is discussed.

Roth, I.↗
DOE OSTI2025

Assessment of Buffer-IPyC Thermomechanical Debonding Behavior Using New Experiment Strength Data in BISON

TRIstructural ISOtropic (TRISO) single sub-millimeter-diameter nuclear fuel particle consists of a spherical fuel kernel surrounded by four coating layers. During reactor operation, the TRISO coating layers are subjected to irradiation-induced dimensional changes and the associated thermomechanical behavior of each layer. One of the observed behaviors is gap formation between the buffer and IPyC layer due to the porous buffer’s irradiation-induced shrinkage exceeding that of the IPyC layer. To predict TRISO thermomechanical performance, multi-physics models have been built that are continually updated and refined. The BISON code is a finite element-based nuclear fuel performance code that may be used for TRISO particle simulations. One of the recent additions to the BISON code is the ability to model the process of layer debonding. This paper will focus on the simulation results of the improved BISON debonding model that will utilize updated strengths measured from irradiated TRISO fuel particles. The new experimental strength data from micromechanical tests of irradiated TRISO fuel samples were exercised in the BISON simulations and compared to baseline strength data to assess their applicability in the models. This also includes updated buffer-IPyC bond strengths to simulate layer delamination. Based on current experimental observations it is noted that the buffer-IPyC separation occurs not exactly at the junction of these two layers, but more on the side of the buffer layer. This improved modeling approach using experimental strength data to characterize buffer-IPyC debonding and its potential subsequent cracking will be presented in the paper along with comparisons to available experimental observations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Validation of the TOPAS Monte Carlo toolkit for HDR brachytherapy simulations

The goal of this work is to validate the user-friendly Geant4-based Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for brachytherapy applications. Brachytherapy simulations performed with TOPAS were systematically compared with published TG-186 reference data. The photon emission energy spectrum, the air-kerma strength, and the dose-rate constant of the model-based dose calculation algorithm (MBDCA)-WG generic Ir-192 source were extracted. For dose calculations, a track-length estimator was implemented. Additionally, the four Joint AAPM/ESTRO/ABG MBDCA-WG test cases were evaluated through histograms of the local and global dose difference volumes. A prostate, a palliative lung, and a breast case were simulated. For each case, the dose ratio map, the histogram of the global dose difference volume, and cumulative dose-volume histograms were calculated. The air-kerma strength was (9.772 ± 0.001) × 10 -8 U Bq -1 (within 0.3% of the reference value). The dose-rate constant was 1.1107 ± 0.0005 cGy h -1 U -1 (within 0.01% of the reference value). For all cases, at least 96.9% of voxels had a local dose difference within [-1%, 1%] and at least 99.9% of voxels had a global dose difference within [-0.1%, 0.1%]. The implemented track-length estimator scorer was more efficient than the default analog dose scorer by a factor of 237. For all clinical cases, at least 97.5% of voxels had a global dose difference within [-1%, 1%]. Dose-volume histograms were consistent with the reference data. TOPAS was validated for high-dose-rate brachytherapy simulations following the TG-186 recommended approach for MBDCAs. Built on top of Geant4, TOPAS provides broad access to a state-of-the-art Monte Carlo code for brachytherapy simulations.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1990

Simulation of electrostatic turbulence due to sheared flows parallel and transverse to the magnetic field

A spatially two-dimensional electrostatic particle simulation code is used to examine the stability of a plasma equilibrium characterized by a localized transverse dc electric field and a magnetic-field-aligned electron drift for L much less than Lx, where Lx is the simulation length in the x direction and L is the scale length associated with the dc electric field. It is found that the dc electric field and the field-aligned current can together play a synergistic role to enable the excitation of electrostatic waves even when the threshold values of the field-aligned drift and the E x B drift are individually subcritical. The simulation results indicate that a broadband turbulence is associated with such an equilibrium.

Nishikawa, K.-I.↗
NASA NTRS1988

Simulations of active ion injection experiments on ARCS 3

The ARCS 3 sounding rocket experiment used two argon guns oriented parallel and perpendicular to the geomagnetic field. The parameters of this experiment have been used in a series of particle simulations in order to interpret ion beam interaction with the multispecies ionospheric plasma. Along with the lower hybrid mode between the electron and the hydrogen gyrofrequencies, ion-ion hybrid modes between the hydrogen and helium gyrofrequencies and between the helium and oxygen gyrofrequencies are identified in the simulations. Besides these cold plasma modes, warm Bernstein modes of the respective ion species are identified. The parallel argon beam injections excite lower hybrid waves perpendicular to the magnetic field because of the finite thermal spread of the parallel beam. Perpendicular ion heating of background hydrogen, helium, and oxygen is comparable in the simulations. Comparisons are made with the wave and particle data from the ARCS 3 rocket experiment.

Hudson, M. K.↗
NASA NTRS1984

Windblown sand on Venus - Preliminary results of laboratory simulations

Small particles and winds of sufficient strength to move them have been detected from Venera and Pioneer-Venus data and suggest the existence of aeolian processes on Venus. The Venus wind tunnel (VWT) was fabricated in order to investigate the behavior of windblown particles in a simulated Venusian environment. Preliminary results show that sand-size material is readily entrained at the wind speeds detected on Venus and that saltating grains achieve velocities closely matching those of the wind. Measurements of saltation threshold and particle flux for various particle sizes have been compared with theoretical models which were developed by extrapolation of findings from Martian and terrestrial simulations. Results are in general agreement with theory, although certain discrepancies are apparent which may be attributed to experimental and/or theoretical-modeling procedures. Present findings enable a better understanding of Venusian surface processes and suggest that aeolian processes are important in the geological evolution of Venus.

Greeley, R.↗
DOE OSTI2022

On the perturbation effect and LET dependence of beam quality correction factors in carbon ion beams

Background: In a recent study, we reported beam quality correction factors, f Q , in carbon ion beams using Monte Carlo (MC) methods for a cylindrical and a parallel-plate ionization chamber (IC). A non-negligible perturbation effect was observed; however, the magnitude of the perturbation correction due to the specific IC subcomponents was not included. Furthermore, the stopping power data presented in the International Commission on Radiation Units and Measurements (ICRU) report 73 were used, whereas the latest stopping power data have been reported in the ICRU report 90. Purpose: The aim of this study was to extend our previous work by computing f Q correction factors using the ICRU 90 stopping power data and by reporting IC-specific perturbation correction factors. Possible energy or linear energy transfer (LET) dependence of the f Q correction factor was investigated by simulating both pristine beams and spread-out Bragg peaks (SOBPs). Methods: The TOol for PArticle Simulation (TOPAS)/GEANT4 MC code was used in this study. A 30 × 30 × 50 cm 3 water phantom was simulated with a uniform 10 × 10 cm 2 parallel beam incident on the surface. A Farmer-type cylindrical IC (Exradin A12) and two parallel-plate ICs (Exradin P11 and A11) were simulated in TOPAS using the manufacturer-provided geometrical drawings. The f Q correction factor was calculated in pristine carbon ion beams in the 150–450 MeV/u energy range at 2 cm depth and in the middle of the flat region of four SOBPs. The k Q correction factor was calculated by simulating the f Qo correction factor in a 60 Co beam at 5 cm depth. The perturbation correction factors due to the presence of the individual IC subcomponents, such as the displacement effect in the air cavity, collecting electrode, chamber wall, and chamber stem, were calculated at 2 cm depth for monoenergetic beams only. Additionally, the mean dose-averaged and track-averaged LET was calculated at the depths at which the f Q was calculated. Results: The ICRU 90 f Q correction factors were reported. The p dis correction factor was found to be significant for the cylindrical IC with magnitudes up to 1.70%. The individual perturbation corrections for the parallel-plate ICs were <1.0% except for the A11 p cel correction at the lowest energy. The f Q correction for the P11 IC exhibited an energy dependence of >1.00% and displayed differences up to 0.87% between pristine beams and SOBPs. Conversely, the f Q for A11 and A12 displayed a minimal energy dependence of <0.50%. The energy dependence was found to manifest in the LET dependence for the P11 IC. A statistically significant LET dependence was found only for the P11 IC in pristine beams only with a magnitude of <1.10%. Conclusions: The perturbation and k Q correction factor should be calculated for the specific IC to be used in carbon ion beam reference dosimetry as a function of beam quality.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

Assessment of buffer-IPyC thermomechanical debonding behavior using new experimental strength data in BISON

TRIstructural ISOtropic (TRISO) fuel is a nuclear fuel commonly used in High Temperature Gas-cooled Reactors (HTGRs). A single sub-millimeter-diameter TRISO fuel particle consists of a spherical fuel kernel surrounded by four coating layers: a low-density pyrocarbon buffer layer, an inner pyrolytic carbon (IPyC) layer, a silicon carbide (SiC) layer, and an outer pyrolytic carbon (OPyC) layer. The kernel is commonly made of UO2 or a mixture of uranium carbide and uranium oxide (UCO). During reactor operation, the TRISO coating layers are subjected to irradiation-induced dimensional changes and the associated thermomechanical behavior of each layer. One of the observed behaviors is gap formation between the buffer and IPyC layer due to the porous buffer’s irradiation-induced shrinkage exceeding that of the IPyC layer. Not all irradiated particles will experience buffer-IPyC gap formation. The debonding may be partial, or it may be nearly total. However, from post-irradiation examination of UCO TRISO fuels irradiated as part of the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program, it was concluded that partial buffer-IPyC debonding was the most common type of buffer-IPyC interaction. To predict TRISO thermomechanical performance, multi-physics models have been built that are being continually updated and refined. The BISON code is a finite element-based nuclear fuel performance code that may be used for 1D, 2D, and 3D TRISO particle simulations. This code is used to calculate fuel temperature, kernel swelling, buffer densification, thermal and irradiation creep, fracture, and fission gas production and release. One of the recent additions to the BISON code is the ability to model the process of layer debonding. This paper will focus on the simulation results of the improved BISON debonding model that will utilize updated strengths measured from irradiated AGR TRISO fuel particles. The new experimental strength data from micromechanical tests of irradiated TRISO fuel samples were exercised in the BISON simulations and compared to baseline strength data to assess their applicability in the models. This also includes updated buffer-IPyC bond strengths to simulate layer delamination. Based on current experimental observations it is noted that the buffer-IPyC separation occurs not exactly at the junction of these two layers, but more on the side of the buffer layer. That observation is also implemented in the TRISO interface debonding model. This improved modeling approach using experimental strength data to characterize buffer-IPyC debonding and its potential subsequent cracking will be presented in the paper along with comparisons to available experimental observations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Simulation studies of plasma waves in the electron foreshock - The generation of Langmuir waves by a gentle bump-on-tail electron distribution

Particle simulation experiments were used to study the basic physical ingredients needed for building a global model of foreshock wave phenomena. In particular, the generation of Langmuir waves by a gentle bump-on-tail electron distribution is analyzed. It is shown that, with appropriately designed simulations experiments, quasi-linear theory can be quantitatively verified for parameters corresponding to the electron foreshock.

Dum, C. T.↗
NASA NTRS1991

Simulation of whistler waves excited in the presence of a cold plasma cloud - Implications for the CRRES mission

A one-dimensional electromagnetic particle simulation model is constructed to study the excitation of whistler waves in the presence of a cold plasma cloud for conditions representative of those after the release of lithium in the inner plasma sheet during the Combined Release and Radiation Effect Satellite mission. The results indicate that a standing-wave pattern with discrete wave frequencies is formed within the cloud. The magnetic wave amplitude inside the cloud, which is limited by quasi-linear diffusion, is of the order of several nanoteslas. Assuming a magnetospheric loss cone of 5 deg, the observed pitch angle diffusion produced by the whistler waves is sufficient to put the electrons on strong diffusion.

Pritchett, P. L.↗
NASA NTRSDate not supplied

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↗
NASA NTRS1985

A simulation study of kilometric radiation generation along an auroral field line

A model for the high-altitude auroral zone is combined with local particle simulations in order to assess the quantitative implications of relativistic dispersion for the generation of auroral kilometric radiation. The auroral zone model is described, including the distribution functions for the primary and secondary electron populations and the variations of the total number density with altitude. The analysis shows that the most intense cyclotron maser emission should occur in the 1.75-2.0 R(E) altitude range. At lower altitudes, the secondary electrons are dominant, and the instability saturates at a low level. At higher altitudes, the energy of the primary electron is lower, and again the radiation level drops. The simulations indicate that linear growth rates of the order of 0.002 Omega(e) occur in this region and that the conversion of primary electron energy into AKR is about one percent.

Pritchett, P. L.↗
NASA NTRS1986

A simulation study of the critical ionization velocity process

The critical ionization velocity process is studied by first investigating a coupled system of equations describing the production of several ion species and electrons by impact ionization, their collisions with neutrals, and the heating of electrons. Analytic relations derived from this were tested with the help of a particle simulation, including collisional processes between neutrals and plasma particles. It was found that resistive heating of electrons plays an important role when the density of the neutrals is high, and that electron heating due to lower hybrid waves is significant when the neutral density is low. In both cases, the control of the plasma production rate by the ratio of the beam velocity to the critical velocity was verified.

Machida, S.↗