Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level Of Conservation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Site 300 B-851 Permit: Radiological Assumptions and Source Term - Air Dispersion Modeling for Radiological Dose Assessment

The following modeling assumptions and source term that is provided in the attachments to this document are meant to serve as a guide to the presentation delivered on April 23rd, 2020, to the San Joaquin Air District Permit Services slides. This analysis is in response to the District’s “Project Incompleteness Letter” (Facility Number: N-472 and Project Number: N-1173492) regarding the potential for depleted uranium (DU) to be resuspended from areas of contaminated soil at Site 300 as a result of high explosive (HE) experiments increased to a requested permitted maximum of 1,000 pound detonations with an annual HE throughput not to exceed 7,500 pounds per year. LLNL has taken every step in the radiological health risk concerns in the modeling assumptions and final source term to be at a level of maximum conservativeness. As an example, forty soil samples were taken in the B-851 area (LLNL-AR-738708 – 2018) and assayed for uranium-234, uranium-235, and uranium-238 concentrations. The maximum uranium soil concentrations were used in the assumptions (potential for resuspended dust occurring from detonation experiments) and not the use of median uranium soil concentrations – a more realistic outcome. Additionally, these uranium soil concentrations were considered to be at 100% uptake in the vegetation for prescribed burns; vegetation uptake factors were not used as that would lower the final source term. This level of conservativeness was maintained throughout the modeling assumptions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates, bilaterian animals, and filamentous fungi). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ~400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization. We, in this study, found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ~100 million years, the microsynteny of many genes—especially genes in metabolic clusters (e.g., in the GAL gene cluster)—was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A Conservative Finite-Volume Based Interface-Tracking Algorithm Using the Signed Distance Function

Methods for tracking an interface between two fluid phases are developed to ensure desired fluid properties, conservation, and stability are preserved in a finitevolume (FV) discretization. Here, the interface is tracked using a level set method where the signed distance function implicitly defines the interface. Marching methods are used to evaluate the value of the signed distance function, including a novel initialization method to initialize any implicit function to the signed distance function around sharp corners in the level set. Global conservation and consistency with a set of governing equations is enforced by a compression coefficient that measures the volumetric compression or expansion due to inaccuracies in the level set evaluation. A redistribution method is integrated into the volume correction to eliminate the small-cell instability while maintaining global conservation. This suite of methods is implemented and tested using static uniform velocity, and potential flow cases with multiple interface geometries. Results show these methods achieve up to second order accuracy, and are conservative. The application for these methods is intended to track the interface of a 3D printing filament in a finite-volume discretization of the all-speed Navier-Stokes equations.

42 ENGINEERING↗
DOE OSTI2023

Techno-economic analysis of the water, energy, and greenhouse gas emissions impacts from the adoption of water efficiency practices in the U.S. manufacturing sector

Water is a critical resource for all of manufacturing. As such, mitigating risks to the supply of water to manufacturing plants facilitates maintaining a vital manufacturing sector. Reductions in water use through water efficiency is one method for mitigating risk, however little is known about its potential in the United States (U.S.) manufacturing sector. Herein, we address this gap and calculate the water reduction potential in U.S. manufacturing from commercially available water efficiency opportunities specific to the manufacturing sector. A list of water efficiency measures specific to the manufacturing sector is compiled and applied to estimates of U.S. manufacturing water withdrawals. The energy and carbon dioxide equivalent (CO 2 e) emission impacts from implementation of these measures are also quantified to better inform the resource, environmental, and economic trade-offs from implementation of these measures. We find that the water savings potential from the evaluated measures is up to 60 percent of U.S. manufacturing water withdrawal for certain subsectors, with the energy and emissions impact dependent on the measure. Further, we calculate the levelized cost of conserved water for several water efficiency measures (-$\$4,000$ to $\$85$/thousand m 3 , where negative values indicate net cost savings and positive values indicate net cost expenditures) and compare them to the levelized cost of alternative water sources ($\$600$ to $\$3,400$/thousand m 3 ). The results show significant opportunities for water and energy use reductions at levelized costs at least one order of magnitude lower than alternative water supplies, with some being revenue-generating.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Vulnerability of Wetlands Due to Projected Sea-Level Rise in the Coastal Plains of the South and Southeast United States

Coastal wetlands are vulnerable to accelerated sea-level rise, yet knowledge about their extent and distribution is often limited. We developed a land cover classification of wetlands in the coastal plains of the southern United States along the Gulf of Mexico (Texas, Louisiana, Mississippi, Alabama, and Florida) using 6161 very-high (2 m per pixel) resolution WorldView-2 and WorldView-3 satellite images from 2012 to 2015. Area extent estimations were obtained for the following vegetated classes: marsh, scrub, grass, forested upland, and forested wetland, located in elevation brackets between 0 and 10 m above sea level at 0.1 m intervals. Sea-level trends were estimated for each coastal state using tide gauge data collected over the period 1983–2021 and projected for 2100 using the trend estimated over that period. These trends were considered conservative, as sea level rise in the region accelerated between 2010 and 2021. Estimated losses in vegetation area due to sea level rise by 2100 are projected to be at least 12,587 km 2 , of which 3224 km 2 would be coastal wetlands. Louisiana is expected to suffer the largest losses in vegetation (80%) and coastal wetlands (75%) by 2100. Such high-resolution coastal mapping products help to guide adaptation plans in the region, including planning for wetland conservation and coastal development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Worldwide Population Genomics Reveal Long-Term Stability of the Mitochondrial Genome Architecture in a Keystone Marine Plant

Mitochondrial genomes (mitogenomes) of flowering plants are composed of multiple chromosomes. Recombination within and between the mitochondrial chromosomes may generate diverse DNA molecules termed isoforms. The isoform copy number and composition can be dynamic within and among individual plants due to uneven replication and homologous recombination. Nonetheless, despite their functional importance, the level of mitogenome conservation within species remains understudied. Whether the ontogenetic variation translates to evolution of mitogenome composition over generations is currently unknown. Here we show that the mitogenome composition of the seagrass Zostera marina is conserved among worldwide populations that diverged ca. 350,000 years ago. Using long-read sequencing, we characterized the Z. marina mitochondrial genome and inferred the repertoire of recombination-induced configurations. To characterize the mitochondrial genome architecture worldwide and study its evolution, we examined the mitogenome in Z. marina meristematic region sampled in 16 populations from the Pacific and Atlantic oceans. Our results reveal a striking similarity in the isoform relative copy number, indicating a high conservation of the mitogenome composition among distantly related populations and within the plant germline, despite a notable variability during individual ontogenesis. Our study supplies a link between observations of dynamic mitogenomes at the level of plant individuals and long-term mitochondrial evolution.

(sub)stoichiometric shifting↗
DOE OSTI2023

EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations

Here, we present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as “EZ” as shorthand for “Esirkepov meets ZigZag”. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Structural insights into strigolactone catabolism by carboxylesterases reveal a conserved conformational regulation

Phytohormone levels are regulated through specialized enzymes, participating not only in their biosynthesis but also in post-signaling processes for signal inactivation and cue depletion. Arabidopsis thaliana (At) carboxylesterase 15 (CXE15) and carboxylesterase 20 (CXE20) have been shown to deplete strigolactones (SLs) that coordinate various growth and developmental processes and function as signaling molecules in the rhizosphere. Here, we elucidate the X-ray crystal structures of AtCXE15 (both apo and SL intermediate bound) and AtCXE20, revealing insights into the mechanisms of SL binding and catabolism. The N-terminal regions of CXE15 and CXE20 exhibit distinct secondary structures, with CXE15 characterized by an alpha helix and CXE20 by an alpha/beta fold. These structural differences play pivotal roles in regulating variable SL hydrolysis rates. Our findings, both in vitro and in planta, indicate that a transition of the N-terminal helix domain of CXE15 between open and closed forms facilitates robust SL hydrolysis. The results not only illuminate the distinctive process of phytohormone breakdown but also uncover a molecular architecture and mode of plasticity within a specific class of carboxylesterases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Discovery of the new isotope 251 Lr and the α-decay fine structure in 253 L

The products of the 203, 205 Tl ( 50 Ti, 2n) fusion-evaporation reactions were studied using the recently commissioned Argonne Gas-Filled Analyzer at Argonne National Laboratory. Two α-decay activities with energies of 9210(19) and 9246(19) keV and half-lives of $42^{+42}_{-14}$ and $24.4^{+7.0}_{-4.5}$ ms were observed which were followed by the known α decays of 247 Md and 243 Es . They are interpreted as originating from the 1/2 ⁻ [521] and 7/2 ⁻ [514] single-proton Nilsson states in the hitherto unknown isotope 251 Lr . From the measured Q α values the 1/2 ⁻ level was placed 117(27) keV above the 7/2 - level in 251 Lr in contrast to 255 Lr where the 1/2 ⁻ level is the lowest. Also, the α decay of 253 Lr was studied in more detail and a new α line at 8660(20) keV was found and a new half-life value of 2.46(32) s for an isomeric state in 253 Lr was measured. The 251, 253, 255 Lr Q α values were compared with predictions of various mass models. The relative energies of the 1/2 ⁻ [521] and 7/2 ⁻ [514] single-proton Nilsson states in 251, 253, 255 Lr isotopes were compared with results of the cranking shell model with pairing treated using the particle-number-conserving method. The level separation and, in particular, the level order change between 251 Lr and 255 Lr was reproduced only when the hexacontetrapole deformation ϵ 6 was included in the calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Enrichment of root-associated Streptomyces strains in response to drought is driven by diverse functional traits and does not predict beneficial effects on plant growth

The genus Streptomyces has consistently been found enriched in drought-stressed plant root microbiomes, yet the ecological basis and functional variation underlying this enrichment at the strain and isolate level remain unclear. Using two 16S rRNA sequencing methods with different levels of taxonomic resolution, we confirmed drought-associated enrichment (DE) of Streptomyces in field-grown sorghum roots and identified five closely related but distinct amplicon sequence variants (ASVs) belonging to the genus with variable drought enrichment patterns. From a culture collection of sorghum root endophytes, we selected 12 Streptomyces isolates representing these ASVs for phenotypic and genomic characterization. Whole-genome sequencing revealed substantial variation in gene content, even among closely related isolates, and exometabolomic profiling showed distinct metabolic responses to media supplemented with drought- versus well-watered root tissue. Traits linked to drought survival, including osmotic stress tolerance, siderophore production, and carbon utilization, varied widely among isolates and were not phylogenetically conserved. Using a broader panel of 48 Streptomyces, we demonstrate that DE scores, determined through mono-association experiments in gnotobiotic sorghum systems, showed high variability and lacked correlation with plant growth promotion. Pangenome-wide association identified orthogroups involved in osmolyte transport (e.g., proP) and membrane biosynthesis (e.g., fabG) as positively associated with DE, though most associations lacked phylogenetic signal. Collectively, these results demonstrate that Streptomyces DE is not a conserved genus-level trait but is instead strain-specific and functionally heterogeneous. Furthermore, DE in the root microbiome was shown not to predict beneficial effects on plant growth. This work underscores the need to resolve functional traits at the strain level and highlights the complexity of microbe-host-environment interactions under abiotic stress.

Fonseca-Garcia, Citlali↗
DOE OSTI2026

Photogalvanic Effects in Surface States of Topological Insulators under Perpendicular Magnetic Fields

We present a theoretical study of the nonlinear magneto-optical shift conductivity in the surface states of the prototypical topological insulator Bi$_2$Se$_3$ under a perpendicular quantizing magnetic field. By describing the electronic states as Landau levels and using a perturbative approach, we derive the microscopic expression for the shift conductivity $\sigma^{(2);\alpha\beta\gamma}(-\omega,\omega)$, where $\alpha,\beta,\gamma=\pm$ stand for the circular polarization of light and $\omega$ is the light frequency; the spectra are further decomposed into contributions from the interband and intraband optical transitions, for which the selection rules are identified. Considering that the system possesses $C_3$ point group of symmetry, the nonzero components of the conductivity tensor are $\sigma^{(2);-++}=[\sigma^{(2);+--}]^\ast$. Therefore, a pure circularly polarized light generates zero shift current. In the clean limit, the conductivities are nonzero only for discrete photon energies because of the discrete Landau levels and energy conservation, and they become Lorentzian lineshapes with the inclusion of damping, which relaxes the condition of energy conservation. The dependence of the spectra on the damping parameters, the magnetic fields, and the chemical potentials is investigated in detail. Our results reveal that the shift current is highly tunable by the chemical potential and the magnetic field. Furthermore, these results underscore the potential of topological insulators for tunable, strong nonlinear magneto-optical applications.

Landau level↗
DOE OSTI2020

Harvest Intensity Effects on Carbon Stocks and Biodiversity Are Dependent on Regional Climate in Douglas-Fir Forests of British Columbia

Temperate forests provide crucial ecosystems services as living sinks for atmospheric carbon (C) and repositories of biodiversity. Applying harvesting at intensities that minimize losses offers one means for mitigating global change. However, little is known of overstory retention levels that best conserve ecosystem services in different regional climates, and likewise as climate changes. To quantify the effect of harvest intensity on C stocks and biodiversity, we compared five harvesting intensities (clearcutting, seedtree retention, 30% patch retention, 60% patch retention, and uncut controls) across a climatic aridity gradient that ranged from humid to semi-arid in the Douglas-fir (Pseudotsuga menziesii) forests of British Columbia. We found that increased harvesting intensity reduced total ecosystem, above ground, and live tree C stocks 1 year post-harvest, and the magnitude of these losses were negatively correlated with climatic aridity. In humid forests, total ecosystem C ranged from 50% loss following clearcut harvest, to 30% loss following large patch retention harvest. In arid forests this range was 60 to 8% loss, respectively. Where lower retention harvests are sought, the small patch retention treatment protected both C stocks and biodiversity in the arid forests, whereas the seedtree method performed as well or better in the humid forests. Below ground C stocks declined by an average of 29% after harvesting, with almost all of the loss from the forest floor and none from the mineral soil. Of the secondary pools, standing and coarse deadwood declined in all harvesting treatments regardless of cutting intensity or aridity, while C stocks in fine fuels and stumps increased. The understory plant C pool declined across all harvesting intensities in the humid forests, but increased in arid forests. Shannon’s diversity and richness of tree and bryoid species declined with harvesting intensity, where tree species losses were greatest in the humid forests and bryoid losses greatest in arid forests. Shrub and herb species were unaffected. This study showed that the highest retention level was best at reducing losses in C stocks and biodiversity, and clearcutting the poorest, and while partial retention of canopy trees can reduce losses in these ecosystem services, outcomes will vary with climatic aridity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Computationally Efficient Linear Semi-Lagrangian Scheme for the Advection of Microphysical Variables in Cloud-Resolving Models

An efficient semi-Lagrangian (SL) numerical scheme for the advection of positive-definite microphysical variables in cloud-resolving models is proposed. The scheme is a linear combination of SL schemes of the first and second order. The proposed scheme is compared with two high-order, monotonic, nonlinear Eulerian schemes, in two 2D tests. In the first test, advection of a positive scalar field with large gradients is performed within complicated idealized steering flows. In the second test, a hail storm is simulated using the 2D Hebrew University Cloud Model, describing the interactions between 15 size distributions of different microphysical quantities (each distribution is described by 43 mass bins). The proposed scheme is computationally efficient, has a low numerical diffusion and a high level of mass conservation accuracy, and preserves the sum of multiple advected variables as well as the linear relationships between the advected scalar variables. Although the proposed SL scheme is not exactly positive definite, the negative values and their number, which may appear only in the case of very strong gradients, are negligible. The scheme produces results similar to those of nonlinear Eulerian advection schemes while reducing computation time by approximately a factor of 10.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multi-Agent Hierarchical Deep Reinforcement Learning for HVAC Control With Flexible DERs

As electricity consumption in commercial and residential buildings continues to rise, reducing energy costs presents an increasing challenge. Heating, ventilating, and air-conditioning (HVAC) systems, which typically account for 40%-50% of a building's energy use, are prime targets for energy savings. Intelligent control of HVAC temperature through the exploitation of HVAC load flexibility brings significant potential to reduce energy consumption and electricity expenses. The nonlinear models of HVAC systems challenge traditional control methods, while the uncertainty introduced by HVAC load flexibility complicates distributed energy resource (DER) management using conventional optimal dispatch techniques. In response to these challenges, we propose a hierarchical multi-agent deep reinforcement learning (DRL) approach. The lower-level agents focus on balancing comfort and energy conservation, while the upper-level DRL agents optimize the use of DERs to reduce peak demand based on the control outcomes of the HVAC by the lower-level agents. Here, in the upper-level agents, we incorporate a multi-agent structure based on ensemble learning, which acts based on historical and current data without relying on precise load forecasting to address the delayed rewarding issue in DRL. This allows for the effective reduction of energy costs. The proposed method is tested using a real-world microgrid comprising 413 buildings in Southern California, and the results demonstrate that our approach can significantly reduce overall electricity bills while ensuring the comfort of consumers and residents.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Mechanisms underlying leaf photosynthetic acclimation to warming and elevated CO 2 as inferred from least‐cost optimality theory

Abstract The mechanisms responsible for photosynthetic acclimation are not well understood, effectively limiting predictability under future conditions. Least‐cost optimality theory can be used to predict the acclimation of photosynthetic capacity based on the assumption that plants maximize carbon uptake while minimizing the associated costs. Here, we use this theory as a null model in combination with multiple datasets of C 3 plant photosynthetic traits to elucidate the mechanisms underlying photosynthetic acclimation to elevated temperature and carbon dioxide (CO 2 ). The model‐data comparison showed that leaves decrease the ratio of the maximum rate of electron transport to the maximum rate of Rubisco carboxylation ( J max / V cmax ) under higher temperatures. The comparison also indicated that resources used for Rubisco and electron transport are reduced under both elevated temperature and CO 2 . Finally, our analysis suggested that plants underinvest in electron transport relative to carboxylation under elevated CO 2 , limiting potential leaf‐level photosynthesis under future CO 2 concentrations. Altogether, our results show that acclimation to temperature and CO 2 is primarily related to resource conservation at the leaf level. Under future, warmer, high CO 2 conditions, plants are therefore likely to use less nutrients for leaf‐level photosynthesis, which may impact whole‐plant to ecosystem functioning.

Smith, Nicholas G.↗
DOE OSTI2020

Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation

Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Phenomena-based graph representations and applications to chemical process simulation

Rapid and robust simulation of chemical processes is critical to conduct process design, optimization, techno-economic analysis, and sustainability analysis. Yet, efficiently solving simulation models remains a challenge due to the highly coupled and nonlinear nature of the underlying algebraic equations that capture the physical phenomena taking place in the process (e.g., material and energy conservation, phase equilibrium, reactions). In this work, we show that graph-theoretic representations of the physical phenomena within unit operations can help navigate and decompose equations to systematically identify alternative approaches for fast and robust numerical solutions. Specifically, we present a graph-theoretic abstraction that captures the connectivity between the model variables/equations and use this abstraction to group variables/equations into fundamental phenomena. We show that phenomena-based decomposition of the underlying equations can help decouple nonlinearities and enforce material/energy conservation at the process level to accelerate convergence. The proposed decomposition approach differs from the more traditional sequential modular simulation approach, in which equations are grouped and decomposed by unit operations. We implemented the phenomena-based decomposition in BioSTEAM—an open-source process simulation platform in Python—and demonstrated that this approach can converge a variety of separation process models. Compared to sequential modular simulation, the phenomena-based approach can converge idealized systems faster, but it can be slower for (or even fail to converge) highly coupled and nonideal process systems.

Convergence↗
DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗