Search NASA⌕ Search

SEARCH · Search NASA

Results for “PICS”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Systematic Influence of Electronic Modification of Ligands on the Catalytic Rate of Water Oxidation by a Single-Site Ru-Based Catalyst

Catalytic water oxidation is an important process for the development of clean energy solutions and energy storage. Despite the significant number of reports on active catalysts, systematic control of the catalytic activity remains elusive. In this study, descriptors are explored that can be correlated with catalytic activity. [Ru(tpy)(pic) 2 (H 2 O)](NO 3 ) 2 and [Ru(EtO-tpy)(pic) 2 (H 2 O)](NO 3 ) 2 (where tpy=2,2': 6',2“-terpyridine, EtO-tpy=4'-(ethoxy)-2,':6',2”-terpyridine, pic=4-picoline) are synthesized and characterized by NMR, UV/Vis, EPR, resonance Raman, and X-ray absorption spectroscopy, and electrochemical analysis. Addition of the ethoxy group increases the catalytic activity in chemically driven and photocatalytic water oxidation. Thus, the effect of the electron-donating group known for the [Ru(tpy)(bpy)(H 2 O)] 2+ family is transferable to architectures with a tpy ligand trans to the Ru-oxo unit. Under catalytic conditions, [Ru(EtO-tpy)(pic) 2 (H 2 O)](NO 3 ) 2 displays new spectroscopic signals tentatively assigned to a peroxo intermediate. In conclusion, reaction pathways were analyzed by using DFT calculations. [Ru(EtO-tpy)(pic) 2 (H 2 O)](NO 3 ) 2 is found to be one of the most active catalysts functioning by a water nucleophilic attack mechanism.

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

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

PUMIPic: A mesh-based approach to unstructured mesh Particle-In-Cell on GPUs

Unstructured mesh particle-in-cell, PIC, simulations executing on the current and next generation of massively parallel systems require new methods for both the mesh and particles to achieve performance and scalability on GPUs. The traditional approach to implementing PIC simulations defines data structures and algorithms in terms of particles with a full copy of the unstructured mesh on every process. To effectively scale the unstructured mesh and particles, mesh-based PIC uses the unstructured mesh as the predominant data structure with the particles stored in terms of the mesh entities. Here, this paper details the PUMIPic library, a framework for developing efficient and performance-portable mesh-based PIC simulations on GPU systems. A pseudo physics simulation based on a five-dimensional gyro-kinetic code for modeling plasma physics is used to examine the performance of PUMIPic. Scaling studies of the unstructured mesh partition and number of particles are performed up to 4096 nodes of the Summit system at Oak Ridge National Laboratory. The studies show that mesh-based PIC can utilize a partitioned mesh and maintain scaling up to system limitations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Geometric electrostatic particle-in-cell algorithm on unstructured meshes

We present a geometric particle-in-cell (PIC) algorithm on unstructured meshes for studying electrostatic perturbations with frequency lower than electron gyrofrequency in magnetized plasmas. In this method, ions are treated as fully kinetic particles and electrons are described by the adiabatic response. The PIC method is derived from a discrete variational principle on unstructured meshes. To preserve the geometric structure of the system, the discrete variational principle requires that the electric field is interpolated using Whitney 1-forms, the charge is deposited using Whitney 0-forms and the electric field is computed by discrete exterior calculus. The algorithm has been applied to study the ion Bernstein wave (IBW) in two-dimensional magnetized plasmas. The simulated dispersion relations of the IBW in a rectangular region agree well with theoretical results. In a two-dimensional circular region with fixed boundary condition, the spectrum and eigenmode structures of the IBW are obtained from simulations. We compare the energy conservation property of the geometric PIC algorithm derived from the discrete variational principle with that of previous PIC methods on unstructured meshes. The comparison shows that the new PIC algorithm significantly improves the energy conservation property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Acceleration of the particle-in-cell code Osiris with graphics processing units

Fully relativistic particle-in-cell (PIC) simulations are crucial for advancing our knowledge of plasma physics. Modern supercomputers based on graphics processing units (GPUs) offer the potential to perform PIC simulations of unprecedented scale, but require robust and feature-rich codes that can fully leverage their computational resources. In this work, this demand is addressed by adding GPU acceleration to the PIC code Osiris. An overview of the algorithm, which features a CUDA extension to the underlying Fortran architecture, is given. Detailed performance benchmarks for thermal plasmas are presented, which demonstrate excellent weak scaling on NERSC's Perlmutter supercomputer and high levels of absolute performance. The robustness of the code to model a variety of physical systems is demonstrated via simulations of Weibel filamentation and laser-wakefield acceleration run with dynamic load balancing. Finally, measurements and analysis of energy consumption are provided that indicate that the GPU algorithm is up to ~14 times faster and ~7 times more energy efficient than the optimized CPU algorithm on a node-to-node basis. The described development addresses the PIC simulation community's computational demands both by contributing a robust and performant GPU-accelerated PIC code and by providing insight into efficient use of GPU hardware.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical study of magnetic island coalescence using magnetohydrodynamics with adaptively embedded particle-in-cell model

Collisionless magnetic reconnection typically requires kinetic treatment that is, in general, computationally expensive compared to fluid-based models. In this study, we use the magnetohydrodynamics with an adaptively embedded particle-in-cell (MHD-AEPIC) model to study the interaction of two magnetic flux ropes. This innovative model embeds one or more adaptive PIC regions into a global MHD simulation domain such that the kinetic treatment is only applied in regions where the kinetic physics is prominent. We compare the simulation results among three cases: (1) MHD with adaptively embedded PIC regions, (2) MHD with statically (or fixed) embedded PIC regions, and (3) a full PIC simulation. The comparison yields good agreement when analyzing their reconnection rates and magnetic island separations as well as the ion pressure tensor elements and ion agyrotropy. In order to reach good agreement among the three cases, large adaptive PIC regions are needed within the MHD domain, which indicates that the magnetic island coalescence problem is highly kinetic in nature, where the coupling between the macro-scale MHD and micro-scale kinetic physics is important.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Intermediate-Temperature Proton Exchange Membranes Based on Cerium Ultraphosphate Composited with Polybenzimidazole

This paper reports the rational fabrication and structural, thermal, mechanical and electrochemical characterization of a new type of intermediate-temperature (IT) polymer-inorganic composite (PIC) proton exchange membranes (PEMs) that are made of cerium ultraphosphate (CeP 5 O 14 —CUP) as the solid-state proton conductor composited with a high-temperature (HT) polybenzimidazole (PBI) as the polymeric binder. Here, flexible PBI-CUP PIC membranes with the thickness of ~135 μ m and CUP mass fraction of up to 75% were prepared by solution-casting without additional acid-doping (e.g., phosphoric acid). The proton conductivity of the fabricated IT-PIC-PEMs was up to 5.80 × 10 –2 S cm –1 as measured from a prototype IT PEM fuel cell (PEMFC) operated at 200 °C in the humidified hydrogen and air environment. This type of IT-PIC-PEMs also demonstrated sufficient mechanical strength and flexibility, excellent thermal stability (up to 350 °C), and very good durability of the proton conductivity (within the test duration of 500 h). The present experimental study shows the promising future of the IT-PIC-PEMs for applications in various IT electrochemical processes including IT-PEMFCs, IT-electrolyzers, etc.

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

Dynamical Mass Estimates of the β Pictoris Planetary System through Gaussian Process Stellar Activity Modeling

Nearly 15 yr of radial velocity (RV) monitoring and direct imaging enabled the detection of two giant planets orbiting the young, nearby star β Pictoris. The δ Scuti pulsations of the star, which overwhelm planetary signals, need to be carefully suppressed. In this work, we independently revisit the analysis of the RV data following a different approach than available in the literature to model the activity of the star. We show that a Gaussian process (GP) with a stochastically driven damped harmonic oscillator kernel can model the δ Scuti pulsations. It provides similar results to parametric models but with a simpler framework, using only three hyperparameters. It also enables us to model poorly sampled RV data that were excluded from previous analyses, hence extending the RV baseline by nearly five years. Altogether, the orbit and mass of both planets can be constrained from RV only, which was not possible with the parametric modeling. To characterize the system more accurately, we also perform a joint fit of all available relative astrometry and RV data. Our orbital solutions for β Pic b favor a low eccentricity of 0.029$_{−0.024}^{+0.061}$ and a relatively short period of 21.1$_{−0.8}^{+2.0}$ yr. The orbit of β Pic c is eccentric with 0.206$_{−0.063}^{+0.074}$ with a period of 3.36 ± 0.03 yr. We find model-independent masses of 11.7 ± 1.4 and 8.5 ± 0.5 M Jup for β Pic b and c, respectively, assuming coplanarity. The mass of β Pic b is consistent with the hottest start evolutionary models, at an age of 25 ± 3 Myr. A direct detection of β Pic c would provide a second calibration measurement in a coeval system.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Image retrieval and processing system version 2.0 development work

The Image Retrieval and Processing System (IRPS) is a software package developed at Washington University and used by the NASA Regional Planetary Image Facilities (RPIF's). The IRPS combines data base management and image processing components to allow the user to examine catalogs of image data, locate the data of interest, and perform radiometric and geometric calibration of the data in preparation for analysis. Version 1.0 of IRPS was completed in Aug. 1989 and was installed at several IRPS's. Other RPIF's use remote logins via NASA Science Internet to access IRPS at Washington University. Work was begun on designing and population a catalog of Magellan image products that will be part of IRPS Version 2.0, planned for release by the end of calendar year 1991. With this catalog, a user will be able to search by orbit and by location for Magellan Basic Image Data Records (BIDR's), Mosaicked Image Data Records (MIDR's), and Altimetry-Radiometry Composite Data Records (ARCDR's). The catalog will include the Magellan CD-ROM volume, director, and file name for each data product. The image processing component of IRPS is based on the Planetary Image Cartography Software (PICS) developed by the U.S. Geological Survey, Flagstaff, Arizona. To augment PICS capabilities, a set of image processing programs were developed that are compatible with PICS-format images. This software includes general-purpose functions that PICS does not have, analysis and utility programs for specific data sets, and programs from other sources that were modified to work with PICS images. Some of the software will be integrated into the Version 2.0 release of IRPS. A table is presented that lists the programs with a brief functional description of each.

Slavney, Susan H.↗
NASA NTRS1996

The Beta Pictoris Phenomenon in A-Shell Stars: Detection of Accreting Gas

We present the results of an expanded survey of A-shell stars using IUE high-dispersion spectra and find accreting, circumstellar gas in the line of sight to nine stars, in addition to the previously identified beta Pic, HR 10, and 131 Tau, which can be followed to between +70 and 100 km/s relative to the star. Two of the program stars, HD 88195 and HD 148283, show variable high-velocity gas. Given the small number of IUE spectra for our program stars, detection of high-velocity, accreting gas in 2/3 of the A-shell stars sampled indicates that accretion is an intrinsic part of the A-shell phenomenon and that beta Pic is not unique among main-sequence A stars in exhibiting such activity. Our program stars, as a group, have smaller column densities of high-velocity gas and smaller near-IR excesses compared with beta Pic. These features are consistent with greater central clearing of a remnant debris disk, compared with beta Pic, and suggest that the majority of field A-shell stars are older than beta Pic.

Grady, C. A.↗
NASA NTRS2014

Low-Mass Planar Photonic Imaging Sensor

Continuing on the successful progress of NIAC Phase I, this report summarizes the technical progress achieved under NIAC Phase II during the performance period September 19, 2014-June 18, 2017. During this period, the research team has made the following accomplishments: designed and layout a silica photonic integrated circuit (PIC) as a two baseline interferometric imager; constructed an experiment to utilize the two baselines for complex visibility measurement on a point source and a variable width slit; analyzed and studied the testbed results. (in collaboration with Lockheed Martin); designed and layout Si3N4 PICs for the low-resolution and high-resolution SPIDER telescope; fabricated the multi-layer Si3N4 PIC for low and high resolution SPIDER telescope; characterize the optical throughput and heater response for Si3N4 PIC for low and high resolution SPIDER telescopes; carried out imaging experiments using the Si3N4 PIC low-resolution version (in collaboration with Lockheed Martin); investigated signal-to-noise (SNR) ratio of SPIDER imager compared to the conventional panchromatic imager (in collaboration with Lockheed Martin); fulfilled the SNR simulation upon SPIDER imager (in collaboration with Lockheed Martin).

Yoo, S. J. Ben↗
NASA NTRS2017

Low-Mass Planar Photonic Imaging Sensor

Continuing on the successful progress of NIAC (NASA Innovative Advanced Concepts) Phase I, this report summarizes the technical progress achieved under NIAC Phase II during the performance period September 19, 2014 to June 18, 2017. During this period, the research team has made the following accomplishments: designed and layout a silica photonic integrated circuit (PIC) as a two baselineinterferometric imager; constructed an experiment to utilize the two baselines for complex visibility measurementon a point source and a variable width slit; analyzed and studied the testbed results (in collaboration with Lockheed Martin); designed and layout Si3N4 PICs for the low-resolution and high-resolution SPIDER (Segmented Planar Imaging Detector for Electro-Optical Reconnaissance) telescope; fabricated the multi-layer Si3N4 PIC for low and high resolution SPIDER telescope; characterized the optical throughput and heater response for Si3N4 PIC for low and highresolution SPIDER telescopes; carried out imaging experiments using the Si3N4 PIC low-resolution version (in collaboration with Lockheed Martin); investigated signal-to-noise (SNR) ratio of SPIDER imager compared to the conventional panchromatic imager (in collaboration with Lockheed Martin); fulfilled the SNR simulation upon SPIDER imager (in collaboration with Lockheed Martin).

Yoo, S. J. Ben↗
DOE OSTI2022

On numerical energy conservation for an implicit particle-in-cell method coupled with a binary Monte-Carlo algorithm for Coulomb collisions

We report conventional particle-in-cell (PIC) methods suffer from enhanced numerical heating (explicit PIC) or cooling (semi-implicit PIC) when coupled with a binary Monte-Carlo algorithm for Coulomb collisions. In this work, a fully-implicit θ-PIC scheme (with adjustable time-biasing parameter 1/2 ≤ θ ≤) is considered. The discrete change in energy of a closed system after a time step for this scheme scales with (1/2 - θ)C θ , where C θ is a positive definite quantity that depends on the frequency spectrum of the energy in the fields. Collisions lead to additional energy in the field fluctuations associated with high-frequency light waves produced by a numerical Bremsstrahlung process, which can result in a large increase in the numerical cooling rate for θ > 1/2. However, for θ = 1/2, energy is exactly conserved. The energy in the field fluctuations on long time scales agrees with that calculated using the equipartition theorem for a classical system in thermodynamic equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Assessment of model parameters in MFiX particle-in-cell approach

The limitations in numerical treatment of solids-phase in conventional methods like Discrete Element Model and Two-Fluid Model have facilitated the development of alternative techniques such as Particle-In-Cell (PIC). However, a number of parameters are involved in PIC due to its empiricism. In this work, global sensitivity analysis of PIC model parameters is performed under three distinct operating regimes common in chemical engineering applications, viz. settling bed, bubbling fluidized bed and circulating fluidized bed. Simulations were performed using the PIC method in Multiphase Flow with Interphase eXchanges (MFiX) developed by National Energy Technology Laboratory (NETL). A non-intrusive uncertainty quantification (UQ) based approach is applied using Nodeworks to first construct an adequate surrogate model and then identify the most influential parameters in each case. This knowledge will aid in developing an effective design of experiments and determine optimal parameters through techniques such as deterministic or statistical calibration.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Quantitative proton radiography and shadowgraphy for arbitrary intensities

Charged-particle radiography and shadowgraphy data can be directly inverted to obtain a line-integrated transverse Lorentz force or a line-integrated transverse refractive index gradient if intensity modulations due to scattering and absorption are negligible, and angular deflections are small. We develop a new direct-inversion algorithm based on plasma physics and compare it to a new Monge–Ampère code and an existing power diagram code. The measured or source intensity is represented by electrons subject to drag, and the other intensity by fixed ions. The decrease in kinetic plus electrostatic energy determines convergence. The displacement of the electrons from their initial to their equilibrium positions determines the line-integrated force or refractive index gradient. We have implemented two approaches: PIC (particle in cell) and Lagrangian fluid, in 1-D and 2-D. The PIC code works for arbitrary intensities, can work efficiently in parallel, and can make use of existing codes. The Lagrangian code requires less memory and is faster than the PIC code without massively parallel processing, but fails in 2-D for large intensity modulations. The Monge–Ampère code is by far the fastest in 2-D, without massively parallel processing, but fails for intensities with large voids, high contrast ratios and large deflections across the boundaries, and could not obtain the degree of convergence possible with the PIC code. As a result, the power diagram code was by far the slowest and most memory intensive, and failed for large peaks in the measured intensity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Comparison of atmospheric radionuclide dispersion models for a risk-informed consequence-driven advanced reactor licensing framework

Current nuclear facility emergency planning zones (EPZs) are based on outdated distance-based criteria, predating comprehensive dose and risk-informed frameworks. Recent advancements in simulation tools have permitted the development of site-specific, dose, and risk-based consequence-driven assessment frameworks. This study investigated the computation of advanced reactor (AR) EPZs using two atmospheric dispersion models: a straight-line Gaussian plume model (GPM) and a semi-Lagrangian Particle in Cell (PIC). Two case studies were conducted: (1) benchmarking the NRC SOARCA study for the Peach Bottom Nuclear Generating Station and (2) analyzing an advanced INL Heat Pipe Design A microreactor's end-of-cycle inventory. The dose criteria for both cases were 10 mSv at mean weather conditions and 50 mSv at 95th percentile weather conditions at 96 h post-release. Results demonstrated that GPM and PIC estimated similar mean peak dose levels for large boiling water reactors in the farfield case, placing EPZ limits beyond current regulations. For ARs with source terms remaining in the nearfield, PIC modeling without specific nearfield considerations could result in excessively high doses and inaccurate EPZ designations. PIC dispersion demonstrated an order of magnitude higher estimate of nearfield inhalation dose contribution when compared to GPM results. Furthermore, both models significantly reduced EPZ sizing within the nearfield. Thus, reductions in the AR source term may eliminate the need for a separate EPZ.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2020

The geometric theory of charge conservation in particle-in-cell simulations

In recent years, several gauge-symmetric particle-in-cell (PIC) methods have been developed whose simulations of particles and electromagnetic fields exactly conserve charge. While it is rightly observed that these methods’ gauge symmetry gives rise to their charge conservation, this causal relationship has generally been asserted via ad hoc derivations of the associated conservation laws. In this work, we develop a comprehensive theoretical grounding for charge conservation in gauge-symmetric Lagrangian and Hamiltonian PIC algorithms. For Lagrangian variational PIC methods, we apply Noether’s second theorem to demonstrate that gauge symmetry gives rise to a local charge conservation law as an off-shell identity. For Hamiltonian splitting methods, we show that the momentum map establishes their charge conservation laws. We define a new class of algorithms – gauge-compatible splitting methods – that exactly preserve the momentum map associated with a Hamiltonian system’s gauge symmetry – even after time discretization. This class of algorithms affords splitting schemes a decided advantage over alternative Hamiltonian integrators. We apply this general technique to design a novel, explicit, symplectic, gauge-compatible splitting PIC method, whose momentum map yields an exact local charge conservation law. Finally, our study clarifies the appropriate initial conditions for such schemes and examines their symplectic reduction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Carbamoylation as an Effective Tool in the Analysis of the Soman Nerve Agent Marker Pinacolyl Alcohol in Soil Matrices by EI-GC-MS and LC-HRMS

Pinacolyl alcohol (PA) is a Schedule 2 chemical commonly featured in most proficiency tests (PTs) administered by the Organisation for the Prohibition of Chemical Weapons (OPCW) due to its direct as well as forensic link to the nerve agent Soman. Therefore, its detection by Chemical Weapons Convention (CWC) inspection teams during on-site investigations is a strong indicator of the past or latent presence of Soman in the environment. Its small molecular weight (102), early elution time, and poor ionization profile make PA a challenging analyte to detect particularly at low concentrations (∼1–10 μg/g). In this work, 1,1′-carbonyldimidazole (CDI) has been used to effectively modify PA for the first time in two different soil matrices (Virginia type A soil and silt sediment) at two separate concentrations (1 and 10 μg/g) for its subsequent detection by EI-GC-MS and LC-HRMS methods. For the EI-GC-MS analysis, the PA carbamate derivative (PIC) exhibits improved chromatography relative to PA such as improved peak shape, increased molecular weight (196 for PIC and 102 for PA) and increased retention time (16.9 min for PIC and ∼4.1 min for PA). In addition, the derivatization also improves the detection of PA by LC-HRMS as the PIC product possesses protonation sites (i.e., imidazole ring) relative to none exhibited by PA. More importantly, the carbamoylation proceeds under mild conditions (55 °C, no base) and rapidly (3 h), characteristics that make it an appealing protocol for the analysis of PA during OPCW PTs or real case scenarios particularly in instances where it is present at low concentrations. It is anticipated that the protocol can be applied to the forensic analysis of this important Soman marker in various environmental matrices.

Chemistry↗