Search NASA⌕ Search

SEARCH · Search NASA

Results for “method of cells”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

PIC simulation of compressive and rarefactive dust ion-acoustic solitary waves

The nonlinear propagations of dust ion-acoustic solitary waves in a collisionless four-component unmagnetized dusty plasma system containing nonextensive electrons, inertial negative ions, Maxwellian positive ions, and negatively charged static dust grains have been investigated by the particle-in-cell method. By comparing the simulation results with those obtained from the traditional reductive perturbation method, it is observed that the rarefactive KdV solitons propagate stably at a low amplitude, and when the amplitude is increased, the prime wave form evolves and then gradually breaks into several small amplitude solitary waves near the tail of soliton structure. The compressive KdV solitons propagate unstably and oscillation arises near the tail of soliton structure. The finite amplitude rarefactive and compressive Gardner solitons seem to propagate stably.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Electrostatic and whistler instabilities excited by an electron beam

The electron beam-plasma system is ubiquitous in the space plasma environment. Here in this paper, using a Darwin particle-in-cell method, the excitation of electrostatic and whistler instabilities by a gyrating electron beam is studied in support of recent laboratory experiments. It is assumed that the total plasma frequency (ω pe ) is larger than the electron cyclotron frequency (Ω e ). The fast-growing electrostatic beam-mode waves saturate in a few plasma oscillations by slowing down and relaxing the electron beam parallel to the background magnetic field. Upon their saturation, the finite amplitude electrostatic beam-mode waves can resonate with the tail of the background thermal electrons and accelerate them to the beam parallel velocity. The slower-growing whistler waves are excited in primarily two resonance modes: (a) through Landau resonance due to the inverted slope of the beam electrons in the parallel velocity and (b) through cyclotron resonance by scattering electrons to both lower pitch angles and smaller energies. It is demonstrated that, for a field-aligned beam, the whistler instability can be suppressed by the electrostatic instability due to a faster energy transfer rate between the beam electrons and electrostatic waves. Such a competition of growth between whistler and electrostatic waves depends on the ratio of ω pe /Ω e . In terms of wave propagation, beam-generated electrostatic waves are confined to the beam region, whereas beam-generated whistler waves transport energy away from the beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Interaction of hydrogen with actinide dioxide (111) surfaces

The interaction of atomic and molecular hydrogen with actinide dioxide (AnO 2 , An = U, Np, Pu) (111) surfaces has been studied by DFT+U, where noncollinear 3k antiferromagnetic behaviour and spin-orbit interactions are considered. The adsorption of atomic hydrogen forms a hydroxide group, coupled to the reduction of an actinide ion. The energy of atomic hydrogen adsorption on the UO 2 (0.82 eV), NpO 2 (-0.10 eV), and PuO 2 (-1.25 eV) surfaces has been calculated. The dissociation of molecular hydrogen is not observed, demonstrated to be due to kinetic rather than thermodynamic factors. As a barrier to the formation of a second hydroxyl group, an unusual charge distribution has been shown. This could be a limitation of a (1·1) unit cell method or an artefact of the systems. The recombination of hydrogen ions on the AnO 2 (111) surfaces is preferred over hydroxide formation.

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

Vortex splitting in two-dimensional fluids and non-neutral electron plasmas with smooth vorticity profiles

Initially, elliptical, quasi-two-dimensional (2D) fluid vortices can split into multiple pieces if the aspect ratio is sufficiently large due to the growth and saturation of perturbations known as Love modes on the vortex edge. Presented here are experiments and numerical simulations, showing that the aspect ratio threshold for vortex splitting is significantly higher for vortices with realistic, smooth edges than that predicted by a simple “vortex patch” model, where the vorticity is treated as piecewise constant inside a deformable boundary. The experiments are conducted by exploiting the E × B drift dynamics of collisionless, pure electron plasmas in a Penning–Malmberg trap, which closely model 2D vortex dynamics due to an isomorphism between the Drift–Poisson equations describing the plasmas and the Euler equations describing ideal fluids. The simulations use a particle-in-cell method to model the evolution of a set of point vortices. The aspect ratio splitting threshold ranges up to about twice as large as the vortex patch prediction and depends on the edge vorticity gradient. This is thought to be due to spatial Landau damping, which decreases the vortex aspect ratio over time and, thus, stabilizes the Love modes. Near the threshold, asymmetric splitting events are observed in which one of the split products contains much less circulation than the other. These results are relevant to a wide range of quasi-2D fluid systems, including geophysical fluids, astrophysical disks, and drift-wave eddies in tokamak plasmas.

Physics↗
DOE OSTI2024

Conceptual design of a Doppler spectrometer for 102 m/s cross-field flows in tokamak divertors

The cross-field transport in the scrape-off-layers (SOLs) and divertors in tokamaks is of a similar size to the poloidal component of the parallel flow, thereby significantly impacting the plasma transport there. However, its direct observation has been challenging because the drift velocity (102–103 m/s) is significantly below the detection limit of conventional diagnostics. To realize cross-field ion flow measurement, a variety of systematic uncertainties in the system must be narrowed down. Here, we develop a conceptual design of the Doppler spectrometry that enables us to measure the impurity flows with 102-m/s accuracy based on an in situ wavelength-calibration technique developed in the astrophysics field, the iodine-cell method. We discuss its properties and applicability. In particular, the scaling relation between wavelength accuracy and various spectroscopic parameters is newly presented, which suggests the high importance of the wavelength resolution of the system. Based on transport simulations for the JT-60SA divertor, the feasibility of the system is assessed.

Fujii, Keisuke↗
DOE OSTI2015

Laser vacuum acceleration of a relativistic electron bunch

With regard to the problem of laser acceleration of a relativistic electron bunch we present a scheme of its vacuum acceleration directly by a relativistic intensity laser pulse. The energy of the electron bunch injected into the laser pulse leading edge increases during its coaxial movement to a thin, pulse-reflecting target. The laser-accelerated electrons continue to move free forward, passing through the target. The study of this acceleration scheme in the three-dimensional geometry is verified in a numerical simulation by the particle-in-cell method, which showed that the energy of a part of the electrons can increase significantly compared to the initial one. Restrictions are discussed, which impose limiting values of energy and total charge of accelerated electrons. (superstrong light fields)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗
DOE OSTI2015

Plasma physical parameters along CME-driven shocks. II. Observation–simulation comparison

In this work, we compare the spatial distribution of the plasma parameters along the 1999 June 11 coronal mass ejection (CME)-driven shock front with the results obtained from a CME-like event simulated with the FLIPMHD3D code, based on the FLIP-MHD particle-in-cell method. The observational data are retrieved from the combination of white-light coronagraphic data (for the upstream values) and the application of the Rankine–Hugoniot equations (for the downstream values). The comparison shows a higher compression ratio X and Alfvénic Mach number M{sub A} at the shock nose, and a stronger magnetic field deflection d toward the flanks, in agreement with observations. Then, we compare the spatial distribution of M{sub A} with the profiles obtained from the solutions of the shock adiabatic equation relating M{sub A}, X, and θ{sub Bn} (the angle between the upstream magnetic field and the shock front normal) for the special cases of parallel and perpendicular shock, and with a semi-empirical expression for a generically oblique shock. The semi-empirical curve approximates the actual values of M{sub A} very well, if the effects of a non-negligible shock thickness δ{sub sh} and plasma-to magnetic pressure ratio β{sub u} are taken into account throughout the computation. Moreover, the simulated shock turns out to be supercritical at the nose and sub-critical at the flanks. Finally, we develop a new one-dimensional Lagrangian ideal MHD method based on the GrAALE code, to simulate the ion-electron temperature decoupling due to the shock transit. Two models are used, a simple solar wind model and a variable-γ model. Both produce results in agreement with observations, the second one being capable of introducing the physics responsible for the additional electron heating due to secondary effects (collisions, Alfvén waves, etc.).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

ELECTRODYNAMICS OF AXISYMMETRIC PULSAR MAGNETOSPHERE WITH ELECTRON-POSITRON DISCHARGE: A NUMERICAL EXPERIMENT

We present the first self-consistent global simulations of pulsar magnetospheres with operating e {sup ±} discharge. We focus on the simple configuration of an aligned or anti-aligned rotator. The star is spun up from a zero (vacuum) state to a high angular velocity, and we follow the coupled evolution of its external electromagnetic field and plasma particles using the ''particle-in-cell'' method. A plasma magnetosphere begins to form through the extraction of particles from the star; these particles are accelerated by the rotation-induced electric field, producing curvature radiation and igniting e {sup ±} discharge. We follow the system evolution for several revolution periods, longer than required to reach a quasi-steady state. Our numerical experiment puts to test previous ideas for the plasma flow and gaps in the pulsar magnetosphere. We first consider rotators capable of producing pairs out to the light cylinder through photon-photon collisions. We find that their magnetospheres are similar to the previously obtained force-free solutions with a Y-shaped current sheet. The magnetosphere continually ejects e {sup ±} pairs and ions. Pair creation is sustained by a strong electric field along the current sheet. We observe powerful curvature and synchrotron emission from the current sheet, consistent with Fermi observations of gamma-ray pulsars. We then study pulsars that can only create pairs in the strong-field region near the neutron star, well inside the light cylinder. We find that both aligned and anti-aligned rotators relax to the ''dead'' state with suppressed pair creation and electric currents, regardless of the discharge voltage.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Special Issue: To Commemorate the Legacy of Malcolm J. Andrews—Scientist, Educator, and Leader

This special issue celebrates the scientific legacy of Dr. Malcolm J. Andrews, former Editor of the ASME Journal of Fluids Engineering (2010–2018). Malcolm devoted his career to the development and applications of computational, experimental, and analytical tools to investigate and solve diverse, complex fluid dynamics, and heat transfer problems. Malcolm’s most notable research accomplishments include: (i) contributions to innovative experiments, simulations, and theory to understand buoyancydriven turbulent mixing due to Rayleigh–Taylor instability; (ii) development of the multiphase particle-in-cell method for understanding fragmentation in spray jets, coal-water slurry, and other dense particulate flows; (iii) mixing in nano- and biomaterials in applications ranging from protein glues to carbon nanotube composites; (iv) determining dispersion in proteins to improve therapeutic drug delivery via micro-encapsulation; (v) innovative methods for storing information using nanopinned dipole magnetic vortices in superconducting materials; and (vi) modeling and design modifications for shell-and-tube heat exchangers, as well as heat transfer in rotating systems such as turbine blades. The results of his research have been published in a long series of trend-setting journal and conference papers, book chapters, and patents.

42 ENGINEERING↗
DOE OSTI2015

TU-F-CAMPUS-T-01: Potential of Using Cerium Oxide Nanoparticles (CONP) for Protecting Healthy Tissue During Accelerated Partial Breast Irradiation (APBI)

Purpose: The purpose of this research is to investigate the feasibility of using targeted cerium oxide nanoparticles (CONP) during APBI to protect healthy cells. Methods: In one approach, CONP are assumed to be incorporated in a micrometer-thick polymer film on the surface of routinely used mammosite balloon applicators for sustained in-situ release of the CONP. In case two, CONPs are administered directly into the lumpectomy cavity. The concentration of H{sub 2}O{sub 2} produced by ionizing radiation was estimated from previously published work using short range linear extrapolation. The assessment of CONPs concentration required to absorb corresponding H{sub 2}O{sub 2} to protect healthy tissue was calculated. Fick’s Second law of diffusion was employed to determine the initial concentration of CONP needed to achieve the minimum concentration for radioprotection at distance 1 cm and 2 cm from the lumpectomy cavity during APBI. The study was carried out for different nanoparticle sizes. Results: The initial concentration of CONPs required to get desired radioprotection concentration at 1 cm and 2 cm after 7 days was found to be 0.4089 mg per Kg and 59.7605 g per Kg respectively for 2 nm size nanoparticles. Using concentrations of 5 mg per kg of CONP that have been shown to be used to confer radioprotection (for about 7.97 Gy) in experiments it was observed that 4.5631, 8.5286, 10.9247, 22.0408, 43.6796 and 65.5618 number of days are required to achieve radioprotection at 1 cm for CONP of sizes 2 nm, 3.8 nm, 5 nm, 10 nm, 20 nm, 30 nm, respectively. Conclusion: Our preliminary results show that smaller size (2 nm and 3.8 nm) CONP would be suitable as radio protectant during APBI because they took a reasonable number of days, i.e. less than 10 days to reach tissues of 1 cm or 2 cm thickness.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

In Vivo Assay Reveals Microbial OleA Thiolases Initiating Hydrocarbon and β-Lactone Biosynthesis

Microbially produced β-lactones are found in antibiotic, antitumor, and antiobesity drugs. Long-chain olefinic membrane hydrocarbons have potential utility as fuels and specialty chemicals. The metabolic pathway to both end products share bacterial enzymes denoted as OleA, OleC, and OleD that transform acyl-CoA cellular intermediates into β-lactones. Bacteria producing membrane hydrocarbons via the Ole pathway additionally express a β-lactone decarboxylase, OleB. Both β-lactone and olefin biosynthesis pathways are initiated by OleA enzymes that define the overall structure of the final product. There is currently very limited information on OleA enzymes apart from the single representative from Xanthomonas campestris . In this study, bioinformatic analysis identified hundreds of new, putative OleA proteins, 74 proteins were screened via a rapid whole-cell method, leading to the identification of 25 stably expressed OleA proteins representing seven bacteria phyla.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Exploring life’s hidden majority: microbial dark matter symposium highlights

The Microbial Dark Matter Symposium held on August 28–29, 2025, in Laguna Beach, Orange County, CA, convened a multidisciplinary group of scientists to address the vast unknowns in microbial life—from uncultured taxa and uncharacterized proteins to elusive viruses and spacefaring microbes. Set against a scenic coastal backdrop, the symposium highlighted advances in single-cell genomics, proximity ligation sequencing, and artificial intelligence-ready bioinformatics, while also probing the limits of microbial persistence, metabolism, and ecological distribution. Sessions explored microbial dark matter from multiple dimensions: cultivability, where new strategies are enabling recovery of elusive microbes; functional ambiguity, where metagenomic dark zones are illuminated by computational annotation; and genomic representation, where single-cell methods bridge gaps left by shotgun community sequencing. Researchers shared breakthroughs in identifying atmospheric microbiomes, “dark oxygen” production in groundwater ecosystems, and microbial survival on the International Space Station. The symposium emphasized integration of methods, disciplines, and ecosystems, advancing a collective push to illuminate the microbial dark matter on Earth and beyond. By highlighting emerging tools, pressing questions, and cross-domain insights, the symposium underscored the need for collaborative, open, and adaptive approaches to study the microbial unknown. The meeting marks a pivotal moment in microbiology, where cultivating knowledge of the uncultivated promises transformative understanding of life, everywhere.

Podar, Mircea [ORNL] (ORCID:0000000327760205)↗
DOE OSTI2016

Properties and parameters of the electron beam injected into the mirror magnetic trap of a plasma accelerator

The parameters of the injector of an axial plasma beam injected into a plasma accelerator operating on the basis of gyroresonance acceleration of electrons in the reverse magnetic field are determined. The trapping of the beam electrons into the regime of gyroresonance acceleration is numerically simulated by the particle- in-cell method. The optimal time of axial injection of the beam into a magnetic mirror trap is determined. The beam parameters satisfying the condition of efficient particle trapping into the gyromagnetic autoresonance regime are found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Quantum Corrections and Magnetotransport in 3D Dirac Semimetal Cd{sub 3 –x}Mn{sub x}As{sub 2} Films

Thin films of solid solutions based on the three-dimensional Dirac semimetal Cd{sub 3}As{sub 2} with the addition of manganese are investigated. Cd{sub 3 –x}Mn{sub x}As{sub 2} films (x = 0, 0.05, and 0.1) 50–70 nm in thickness are formed on a glassceramic substrate using vacuum-thermal deposition from cadmium arsenide ingots doped by Mn and fabricated by direct alloying elements by the vacuum-cell method. The temperature and magnetic-field dependences of the resistance are measured and the transport parameters of the films under study are determined. Positive magnetoresistance of the characteristic shape corresponding to the contribution of the weak antilocalization effect is observed for films with x = 0 and 0.05. The contribution from the weak localization effect is observed at a higher Mn content (x = 0.1). This change in the quantum correction type as applied to topological semimetals points to reconstruction of the band structure and transition from the Dirac semimetal state into a trivial semiconductor phase, which corresponds to the critical Mn content x{sub c} ~ 0.07 in this case.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

picassodev/picasso

Picasso is a C++ template library for the development of exascale-class particle-in-cell methods. It contains data structures for particle and grid management using the ECP Kokkos and Cabana libraries, abstractions for particle and grid discretizations, and tools for interface tracking and reconstruction with particles.

Slattery, Stuart↗
DOE OSTI2018

Performance Optimization of a Beam Dynamics PIC Code On Hybrid Computer Architectures

The self-consistent multi-particle tracking based on particle-in-cell method (PIC) has been widely used in particle accelerator beam dynamics study. However, the PIC simulation is time-consuming and needs to use modern parallel computers for high resolution applications. In this paper, we implemented and optimized a parallel beam dynamics PIC code on two types of hybrid parallel computer architectures: one is the GPU and GPU cluster, while the other is the "Knight Landing" CPU cluster.

Liu, Zhicong↗
DOE OSTI2012

P24 Plasma Physics Summer School 2012 Los Alamos National Laboratory Summer lecture series for students [Slides]

This report covers the 2012 LANL summer lecture series for students. The lectures were: (1) Tom Intrator, P24 LANL: Kick off, Introduction - What is a plasma; (2) Bruno Bauer, Univ. Nevada-Reno: Derivation of plasma fluid equations; (3) Juan Fernandez, P24 LANL Overview of research being done in p-24; (4) Tom Intrator, P24 LANL: Intro to dynamo, reconnection, shocks; (5) Bill Daughton X-CP6 LANL: Intro to computational particle in cell methods; (6) Kirk Flippo, P24 LANL: High energy density plasmas; (7) Thom Weber, P24 LANL: Energy crisis, fission, fusion, non carbon fuel cycles; (8) Tom Awe, Sandia National Laboratory: Magneto Inertial Fusion; and (9) Yongho Kim, P24 LANL: Industrial technologies.

Astronomy & Astrophysics(79)↗