SEARCH · Search NASA

Results for “Parallel particle tracing”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Mapping the magnetosheath field between the magnetopause and the bow shock - Implications for magnetospheric particle leakage

An approximate picture of the volumes occupied by particles that originate in the vicinity of the magnetopause is obtained by mapping magnetosheath magnetic field lines which drape over the magnetopause through the bow shock. Subsets of these field lines that connect to potential sites of magnetic merging on the magnetopause are also traced in the event that the particle leakage occurs preferentially where normal components of the field are present across that boundary. The results of this modeling exercise suggest that energetic magnetospheric particles which are not scattered by magnetosheath magnetic fluctuations are likely to exit the magnetosheath in the region of the quasi-parallel shock.

Luhmann, J. G.↗

NASA NTRS2023

Electron Acceleration and Heating during Magnetic Reconnection in the Earth's Quasi-parallel Bow Shock

We perform a 2.5-dimensional particle-in-cell simulation of a quasi-parallel shock, using parameters for the Earth's bow shock, to examine electron acceleration and heating due to magnetic reconnection. The shock transition region evolves from the ion-coupled reconnection dominant stage to the electron-only reconnection dominant stage, as time elapses. The electron temperature enhances locally in each reconnection site, and ion-scale magnetic islands generated by ion-coupled reconnection show the most significant enhancement of the electron temperature. The electron energy spectrum shows a power law, with a power-law index around 6. We perform electron trajectory tracing to understand how they are energized. Some electrons interact with multiple electron-only reconnection sties, and Fermi acceleration occurs during multiple reflections. Electrons trapped in ion-scale magnetic islands can be accelerated in another mechanism. Islands move in the shock transition region, and electrons can obtain larger energy from the in-plane electric field than the electric potential in those islands. These newly found energization mechanisms in magnetic islands in the shock can accelerate electrons to energies larger than the achievable energies by the conventional energization due to the parallel electric field and shock drift acceleration. This study based on the selected particle analysis indicates that the maximum energy in the nonthermal electrons is achieved through acceleration in ion-scale islands, and electron-only reconnection accounts for no more than half of the maximum energy, as the lifetime of sub-ion-scale islands produced by electron-only reconnection is several times shorter than that of ion-scale islands.

Solar magnetic reconnection↗

DOE OSTI2015

Electron acceleration by parallel and perpendicular electric fields during magnetic reconnection without guide field

In this work, electron acceleration due to the electric field parallel to the background magnetic field during magnetic reconnection with no guide field is investigated by theory and two-dimensional electromagnetic particle-in-cell simulations and compared with acceleration due to the electric field perpendicular to the magnetic field. The magnitude of the parallel electric potential shows dependence on the ratio of the plasma frequency to the electron cyclotron frequency as (ω pe /Ω e ) -2 and on the background plasma density as N$^{-½}_{b}$ . In the Earth's magnetotail, the parameter ω pe /Ω e ~9 and the background (lobe) density can be of the order of 0.01 cm -3 , and it is expected that the parallel electric potential is not large enough to accelerate electrons up to 100 keV. Therefore, we must consider the effect of the perpendicular electric field to account for electron energization in excess of 100 keV in the Earth's magnetotail. Trajectories for high-energy electrons are traced in a simulation to demonstrate that acceleration due to the perpendicular electric field in the diffusion region is the dominant acceleration mechanism, rather than acceleration due to the parallel electric fields in the exhaust regions. For energetic electrons accelerated near the X line due to the perpendicular electric field, pitch angle scattering converts the perpendicular momentum to the parallel momentum. On the other hand, for passing electrons that are mainly accelerated by the parallel electric field, pitch angle scattering converting the parallel momentum to the perpendicular momentum occurs. In this way, particle acceleration and pitch angle scattering will generate heated electrons in the exhaust regions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Collaborative Project. 3D Radiative Transfer Parameterization Over Mountains/Snow for High-Resolution Climate Models. Fast physics and Applications

Under the support of the aforementioned DOE Grant, we have made two fundamental contributions to atmospheric and climate sciences: (1) Develop an efficient 3-D radiative transfer parameterization for application to intense and intricate inhomogeneous mountain/snow regions. (2) Innovate a stochastic parameterization for light absorption by internally mixed black carbon and dust particles in snow grains for understanding and physical insight into snow albedo reduction in climate models. With reference to item (1), we divided solar fluxes reaching mountain surfaces into five components: direct and diffuse fluxes, direct- and diffuse-reflected fluxes, and coupled mountain-mountain flux. “Exact” 3D Monte Carlo photon tracing computations can then be performed for these solar flux components to compare with those calculated from the conventional plane-parallel (PP) radiative transfer program readily available in climate models. Subsequently, Parameterizations of the deviations of 3D from PP results for five flux components are carried out by means of the multiple linear regression analysis associated with topographic information, including elevation, solar incident angle, sky view factor, and terrain configuration factor. We derived five regression equations with high statistical correlations for flux deviations and successfully incorporated this efficient parameterization into WRF model, which was used as the testbed in connection with the Fu-Liou-Gu PP radiation scheme that has been included in the WRF physics package. Incorporating this 3D parameterization program, we conducted simulations of WRF and CCSM4 to understand and evaluate the mountain/snow effect on snow albedo reduction during seasonal transition and the interannual variability for snowmelt, cloud cover, and precipitation over the Western United States presented in the final report. With reference to item (2), we developed in our previous research a geometric-optics surface-wave approach (GOS) for the computation of light absorption and scattering by complex and inhomogeneous particles for application to aggregates and snow grains with external and internal mixing structures. We demonstrated that a small black (BC) particle on the order of 1 μm internally mixed with snow grains could effectively reduce visible snow albedo by as much as 5–10%. Following this work and within the context of DOE support, we have made two key accomplishments presented in the attached final report.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1997

Visualization of Unsteady Computational Fluid Dynamics

The current compute environment that most researchers are using for the calculation of 3D unsteady Computational Fluid Dynamic (CFD) results is a super-computer class machine. The Massively Parallel Processors (MPP's) such as the 160 node IBM SP2 at NAS and clusters of workstations acting as a single MPP (like NAS's SGI Power-Challenge array and the J90 cluster) provide the required computation bandwidth for CFD calculations of transient problems. If we follow the traditional computational analysis steps for CFD (and we wish to construct an interactive visualizer) we need to be aware of the following: (1) Disk space requirements. A single snap-shot must contain at least the values (primitive variables) stored at the appropriate locations within the mesh. For most simple 3D Euler solvers that means 5 floating point words. Navier-Stokes solutions with turbulence models may contain 7 state-variables. (2) Disk speed vs. Computational speeds. The time required to read the complete solution of a saved time frame from disk is now longer than the compute time for a set number of iterations from an explicit solver. Depending, on the hardware and solver an iteration of an implicit code may also take less time than reading the solution from disk. If one examines the performance improvements in the last decade or two, it is easy to see that depending on disk performance (vs. CPU improvement) may not be the best method for enhancing interactivity. (3) Cluster and Parallel Machine I/O problems. Disk access time is much worse within current parallel machines and cluster of workstations that are acting in concert to solve a single problem. In this case we are not trying to read the volume of data, but are running the solver and the solver outputs the solution. These traditional network interfaces must be used for the file system. (4) Numerics of particle traces. Most visualization tools can work upon a single snap shot of the data but some visualization tools for transient problems require dealing with time.

Haimes, Robert↗

NASA NTRS1994

Effect of B(sub y) on neutral line ridges and dynamical source ordering

In this paper we study the effect of a uniform magnetic field B(sub y) in the cross-tail y direction on the dynamics and distribution of energetic ions in a current sheet model including a neutral line. Martin and Speiser (1988) have previously shown that a 'ridge' in the velocity space distribution function is a remote neutral line signature in a two-dimensional field without B(sub y). Our results show that as B(sub y) is increased, using nominal tail parameters, there is very little change in the ridge signature for small B(sub y) (up to about 4 times B(sub z)). For intermediate values (up to about 10 times times B(sub Y)) the ridge becomes observable further from the neutral line, while close to the X line the ridge is strongly modified. For large B(sub y) (of the order of B(sub x)) the ridge evolves into the new neutral line signatures which depend strongly on whether the observations point is above or below the current sheet. These results are used to estimate B(sub y) less than 2 nT in the Active Magentospheric Particle Traces Explorers/Ion Release Module (AMPTE/IRM) event modeled by Speiser and Martin (1994), which is consistent with onboard magnetometer measurements. We further show that the cause of the new structures appearing at large B(sub Y) is the same source seperation effect that produced the original ridge and that all these signatures are relatively insensitive to a small parallel electric field. Finally, we point our a dynamical ordering which occurs at large B(sub y): Using high-resolution simulations, we find a common boundary separating regions of positive and negative initial x, y, and z position, as well as initial pitch angle. This boundary is purely dynamical in nature and is independent of modeled source distributions.

Martin, R. F., Jr.↗

NASA NTRS2009

Particle Acceleration and Associated Emission from Relativistic Shocks

Five talks consist of a research program consisting of numerical simulations and theoretical development designed to provide an understanding of the emission from accelerated particles in relativistic shocks. The goal of this lecture is to discuss the particle acceleration, magnetic field generation, and radiation along with the microphysics of the shock process in a self-consistent manner. The discussion involves the collisionless shocks that produce emission from gamma-ray bursts and their afterglows, and producing emission from supernova remnants and AGN relativistic jets. Recent particle-in-cell simulation studies have shown that the Weibel (mixed mode two-stream filamentation) instability is responsible for particle (electron, positron, and ion) acceleration and magnetic field generation in relativistic collisionless shocks. 3-D RPIC code parallelized with MPI has been used to investigate the dynamics of collisionless shocks in electron-ion and electron-positron plasmas with and without initial ambient magnetic fields. In this lecture we will present brief tutorials of RPIC simulations and RMHD simulations, a brief summary of recent RPIC simulations, mechanisms of particle acceleration in relativistic shocks, and calculation of synchrotron radiation by tracing particles. We will discuss on emission from the collisionless shocks, which will be calculated during the simulation by tracing particle acceleration self-consistently in the inhomogeneous magnetic fields generated in the shocks. In particular, we will discuss the differences between standard synchrotron radiation and the jitter radiation that arises in turbulent magnetic fields.

Nishkawa, Ken-Ichi↗

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 OSTI2017

SCIDAC Center for simulation of wave particle interactions CompX participation

Harnessing the energy that is released in fusion reactions would provide a safe and abundant source of power to meet the growing energy needs of the world population. The next step toward the development of fusion as a practical energy source is the construction of ITER, a device capable of producing and controlling the high performance plasma required for self-sustaining fusion reactions, or “burning” plasma. The input power required to drive the ITER plasma into the burning regime will be supplied primarily with a combination of external power from radio frequency waves in the ion cyclotron range of frequencies and energetic ions from neutral beam injection sources, in addition to internally generated Ohmic heating from the induced plasma current that also serves to create the magnetic equilibrium for the discharge. The ITER project is a large multi-billion dollar international project in which the US participates. The success of the ITER project depends critically on the ability to create and maintain burning plasma conditions, it is absolutely necessary to have physics-based models that can accurately simulate the RF processes that affect the dynamical evolution of the ITER discharge. The Center for Simulation of WavePlasma Interactions (CSWPI), also known as RF-SciDAC, is a multi-institutional collaboration that has conducted ongoing research aimed at developing: (1) Coupled core-to-edge simulations that will lead to an increased understanding of parasitic losses of the applied RF power in the boundary plasma between the RF antenna and the core plasma; (2) Development of models for core interactions of RF waves with energetic electrons and ions (including fusion alpha particles and fast neutral beam ions) that include a more accurate representation of the particle dynamics in the combined equilibrium and wave fields; and (3) Development of improved algorithms that will take advantage of massively parallel computing platforms at the petascale level and beyond to achieve the needed physics, resolution, and/or statistics to address these issues. CompX provides computer codes and analysis for the calculation of the electron and ion distributions in velocity-space and plasma radius which are necessary for reliable calculations of power deposition and toroidal current drive due to combined radiofrequency and neutral beam at high injected powers. It has also contributed to ray tracing modeling of injected radiofrequency powers, and to coupling between full-wave radiofrequency wave models and the distribution function calculations. In the course of this research, the Fokker-Planck distribution function calculation was made substantially more realistic by inclusion of finite-width drift-orbit effects (FOW). FOW effects were also implemented in a calculation of the phase-space diffusion resulting from radiofrequency full-wave models. Average level of funding for CompX was approximately three man-months per year.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1996

Non-Linear Cosmological Power Spectra in Real and Redshift Space

We present an expression for the non-linear evolution of the cosmological power spectrum based on Lagrangian trajectories. This is simplified using the Zel'dovich approximation to trace particle displacements, assuming Gaussian initial conditions. The model is found to exhibit the transfer of power from large to small scales expected in self-gravitating fields. Some exact solutions are found for power-law initial spectra. We have extended this analysis into red-shift space and found a solution for the non-linear, anisotropic redshift-space power spectrum in the limit of plane-parallel redshift distortions. The quadrupole-to-monopole ratio is calculated for the case of power-law initial spectra. We find that the shape of this ratio depends on the shape of the initial spectrum, but when scaled to linear theory depends only weakly on the redshift-space distortion parameter, beta. The point of zero-crossing of the quadrupole, kappa(sub o), is found to obey a simple scaling relation and we calculate this scale in the Zel'dovich approximation. This model is found to be in good agreement with a series of N-body simulations on scales down to the zero-crossing of the quadrupole, although the wavenumber at zero-crossing is underestimated. These results are applied to the quadrupole-to-monopole ratio found in the merged QDOT plus 1.2-Jy-IRAS redshift survey. Using a likelihood technique we have estimated that the distortion parameter is constrained to be beta greater than 0.5 at the 95 percent level. Our results are fairly insensitive to the local primordial spectral slope, but the likelihood analysis suggests n = -2 un the translinear regime. The zero-crossing scale of the quadrupole is k(sub 0) = 0.5 +/- 0.1 h Mpc(exp -1) and from this we infer that the amplitude of clustering is sigma(sub 8) = 0.7 +/- 0.05. We suggest that the success of this model is due to non-linear redshift-space effects arising from infall on to caustic and is not dominated by virialized cluster cores. The latter should start to dominate on scales below the zero-crossing of the quadrupole, where our model breaks down.

Taylor, A. N.↗

NASA NTRS2006

The Development Of Enabling Technologies For Submillimeter-Wave Remote Sensing of Ice Clouds From Space

This paper presents the chronological development of technologies and techniques that have led to a satellite mission concept aimed at quantifying the temporal and spatial distributions of upper tropospheric ice clouds. The Submillimeter-wave and Infrared Ice Cloud Experiment (SIRICE) is an Earth System Science Pathfinder mission concept designed to improve our understanding of the upper tropospheric water cycle and its coupling to the Earth s radiation budget. Ice outflow from convective storm systems is known to play an important role in regional energy budgets; however, ice generation and subsequent precipitation and sublimation are poorly quantified. SIRICE will provide measurements of ice cloud distributions and microphysical properties which are needed for understanding the crucial link between the hydrologic and energy cycles. The SIRICE measurement platform is comprised of two integrated instruments, the Submillimeter/millimeter-wave radiometer (SM4) and the Infrared Cloud Ice Radiometer (IRCIR). The primary instrument is the SM4, a conical scanner that provides a 1600 km swath of the Earth's surface at 53 degree incidence. The SM4 has 6 linearly polarized receivers measuring 12 spectral bands centered at 183 GHz, 325 GHz, 448 GHz, 643 GHz and 874 GHz; two receivers at 643 GHz measure horizontal and vertical polarizations. Submillimeter-wavelengths are well suited to the remote sensing of ice clouds due to the relative size of the wavelengths to particle sizes. Upwelling emission from lower tropospheric water vapor is scattered by the ice clouds thus causing a brightness temperature depression at submillimeter wavelengths. The IRCIR is a push broom imager with approximately 1500 km swath and spectral channels at 11 and 12 micrometers. This combination of coincident infrared and submillimeter-wavelength measurements were chosen because of its ability to provide retrieval of ice water path and median particle size for a wide range of ice clouds from thin cirrus to thick anvil structures. Over the past decade there has been a parallel development of submillimeter-wave technologies, demonstration instruments, and remote sensing techniques that have led to the present SIRICE mission concept. Mapping of these developmental paths reveals the origins, rational and maturity of features of the SIRICE payload such as its channel selection, compact design, and multipoint calibration. This presentation traces the evolution of the SIRICE mission concept from the early 1990's to its present status.

Racette, Paul↗

DOE OSTI2020

Threadsafe Dynamic Neighbor Lists for Monte Carlo Ray Tracing

Monte Carlo (MC) transport codes offer high-fidelity modeling of particle transport physics, but their high computational cost makes them impractical for many applications. For some applications such as multiphysics and depletion that use finely discretized geometries, a large portion of this computational cost is attributable to ray tracing. Neighbor lists are a well-known method for accelerating ray-tracing calculations in a MC code, but despite their prevalence, little work has been published on the details of their implementation. The fine details can have a significant impact on performance, particularly when using shared-memory parallelism. This paper addresses these details of implementation with a discussion of different neighbor list schemes and their impact on software runtime. Performance tests were run by using OpenMC on a pin-cell problem discretized with up to 200 axial regions. The results demonstrate that switching from surface-based to cell-based neighbor lists leads to a 10 faster calculation rate for the most fine discretization. Finally, using a threadsafe shared-memory data structure results in a 20% faster calculation rate versus simple threadprivate neighbor lists. Results here show that a data structure that is contiguous in memory improves performance by only 1% to 2% over noncontiguous linked lists.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1992

Simultaneous mapping of the unsteady flow fields by Particle Displacement Velocimetry (PDV)

Current experimental and computational techniques must be improved in order to advance the prediction capability of the longitudinal vortical flows shed by underwater vehicles. The generation, development, and breakdown mechanisms of the shed vortices at high Reynolds numbers are not fully understood. The ability to measure hull separated vortices associated with vehicle maneuvering does not exist at present. The existing point-by-point measurement techniques can only capture approximately the large 'mean' eddies but fail to meet the dynamics of small vortices during the initial stage of generation. A new technique, which offers a previously unavailable capability to measure the unsteady cross-flow distribution in the plane of the laser light sheet, is called Particle Displacement Velocimetry (PDV). PDV consists of illuminating a thin section of the flowfield with a pulsed laser. The water is seeded with microscopic, neutrally buoyant particles containing imbedded fluorescing dye which responds with intense spontaneous fluorescence with the illuminated section. The seeded particles in the vortical flow structure shed by the underwater vehicle are illuminated by the pulse laser and the corresponding particle traces are recorded in a single photographic frame. Two distinct approaches were utilized for determining the velocity distribution from the particle traces. The first method is based on matching the traces of the same particle and measuring the distance between them. The direction of the flow can be identified by keeping one of the pulses longer than the other. The second method is based on selecting a small window within the image and finding the mean shift of all the particles within that region. The computation of the auto-correlation of the intensity distribution within the selected sample window is used to determine the mean displacement of particles. The direction of the flow is identified by varying the intensity of the laser light between pulses. Considerable computational resources are required to compute the auto-correction of the intensity distribution. Parallel processing will be employed to speed up the data reduction. A few examples of measured unsteady vortical flow structures shed by the underwater vehicles will be presented.

Huang, Thomas T.↗

DOE OSTI2024

A closer look in the mirror: reflections on the matter/dark matter coincidence

We argue that the striking similarity between the cosmic abundances of baryons and dark matter, despite their very different astrophysical behavior, strongly motivates the scenario in which dark matter resides within a rich dark sector parallel in structure to that of the standard model. The near cosmic coincidence is then explained by an approximate ℤ$_{2}$ exchange symmetry between the two sectors, where dark matter consists of stable dark neutrons, with matter and dark matter asymmetries arising via parallel WIMP baryogenesis mechanisms. Taking a top-down perspective, we point out that an adequate ℤ$_{2}$ symmetry necessitates solving the electroweak hierarchy problem in each sector, without our committing to a specific implementation. A higher-dimensional realization in the far UV is presented, in which the hierarchical couplings of the two sectors and the requisite ℤ$_{2}$-breaking structure arise naturally from extra-dimensional localization and gauge symmetries. We trace the cosmic history, paying attention to potential pitfalls not fully considered in previous literature. Residual ℤ$_{2}$-breaking can very plausibly give rise to the asymmetric reheating of the two sectors, needed to keep the cosmological abundance of relativistic dark particles below tight bounds. We show that, despite the need to keep inter-sector couplings highly suppressed after asymmetric reheating, there can naturally be order-one couplings mediated by TeV scale particles which can allow experimental probes of the dark sector at high energy colliders. Massive mediators can also induce dark matter direct detection signals, but likely at or below the neutrino floor.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1986

A test of Lee's quasi-linear theory of ion acceleration by interplanetary traveling shocks

Lee's (1983) quasi-linear theory of ion acceleration is tested using ISEE-3 measurements of the November 12, 1978 quasi-parallel interplanetary shock. His theory accounts with varying degrees of precision for the energetic proton spatial profiles; the dependence of the spectral index of the power law proton velocity distribution upon the shock compression ratio; the power law dependence of the upstream proton scalelength upon energy; the absolute magnitude of the upstream proton scale length; the behavior of the energetic proton anisotropy upstream and downstream of the shock; the behavior of the alpha-particle proton ratio upstream; the equality of the spatial scale lengths at the shock of the upstream waves and of the protons that resonate with them; and the dependence of the integrated wave energy density upon the proton energy density at the shock. However, the trace magnetic field frequency spectra disagree with his theory in two ways. The part of the spectrum that can resonate with the observed protons via first-order cyclotron resonance is flat, whereas Lee's theory predicts an f exp - 7/4 frequency dependence for the November 12 shock. Higher frequency waves, which could not resonate with the observed upstream protons, increased in amplitude as the shock approached, suggesting that they too were generated by the shock.

Kennel, C. F.↗

DOE OSTI2017

Generation of Highly Oblique Lower Band Chorus Via Nonlinear Three-Wave Resonance

Chorus in the inner magnetosphere has been observed frequently at geomagnetically active times, typically exhibiting a two-band structure with a quasi-parallel lower band and an upper band with a broad range of wave normal angles. But recent observations by Van Allen Probes confirm another type of lower band chorus, which has a large wave normal angle close to the resonance cone angle. It has been proposed that these waves could be generated by a low-energy beam-like electron component or by temperature anisotropy of keV electrons in the presence of a low-energy plateau-like electron component. This paper, however, presents an alternative mechanism for generation of this highly oblique lower band chorus. Through a nonlinear three-wave resonance, a quasi-parallel lower band chorus wave can interact with a mildly oblique upper band chorus wave, producing a highly oblique quasi-electrostatic lower band chorus wave. This theoretical analysis is confirmed by 2-D electromagnetic particle-in-cell simulations. Furthermore, as the newly generated waves propagate away from the equator, their wave normal angle can further increase and they are able to scatter low-energy electrons to form a plateau-like structure in the parallel velocity distribution. As a result, the three-wave resonance mechanism may also explain the generation of quasi-parallel upper band chorus which has also been observed in the magnetosphere.

58 GEOSCIENCES↗

NASA NTRS2022

Non-Protein Amino Acids Identified in Carbon-Rich Hayabusa Particles

Amino acid abundances in acid-hydrolyzed hot water extracts of gold foils containing five Category 3 (carbon-rich) Hayabusa particles were studied using liquid chromatography with tandem fluorescence and accurate mass detection. Initial particle analyses using field emission scanning electron microscopy with energy-dispersive X-ray spectrometry indicated that the particles were composed mainly of carbon. Prior to amino acid analysis, infrared and Raman microspectroscopy showed some grains possessed primitive organic carbon. Although trace terrestrial contamination, namely L-protein amino acids, was observed in all Hayabusa extracts, several terrestrially uncommon non-protein amino acids were also identified. Some Hayabusa particles contained racemic (D≈L) mixtures of the non-protein amino acids β-aminoisobutyric acid (β-AIB) and β-amino-n-butyric acid (β-ABA) at low abundances ranging from 0.09 to 0.31 nmol g −1 . Larger abundances of the non-protein amino acid β-alanine (9.2 nmol g −1 , ≈4.5 times greater than background levels) were measured in an extract of three Hayabusa particles. This β-alanine abundance was ≈6 times higher than that measured in an extract of a CM2 Murchison grain processed in parallel. The comparatively high β-alanine abundance is surprising as asteroid Itokawa is similar to amino acid-poor LL ordinary chondrites. Elevated β-alanine abundances and racemic β-AIB and β-ABA in Hayabusa grains suggested these compounds have non-biological and plausibly non-terrestrial origins. These results are the first evidence of plausibly extraterrestrial amino acids in asteroid material from a sample-return mission and demonstrate the capabilities of the analytical protocols used to study asteroid Ryugu and Bennu samples returned by the JAXA Hayabusa2 and NASA OSIRIS-REx missions, respectively.

Eric T. Parker↗

DOE OSTI2026

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗