Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Periodic systems”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Method for Projecting Cloud Shadows Onto a Central Receiver Field to Predict Receiver Damage

This work demonstrates methods of mapping high-spatial-resolution direct normal irradiance (DNI) data from satellites, Total Sky Imagers (TSIs), and analogous data sources onto a heliostat field for characterizing the spatial and temporal variation of the incident flux on a central receiver tower during cloud transient events. The mapping methods are incorporated into an optical software module that interfaces with CoPylot–SolarPILOT’s python API– to provide computationally efficient optical simulation of the heliostat field and the solar power tower. Eventually, this optical model will be incorporated into optimization models whereby a plant operator can understand the effects of cloud transient events on overall power production and receiver lifetime due to creep-fatigue damage and therefore make better informed decisions about receiver shutdown events. By more accurately modelling the effects of cloud events on receiver flux maps, this work may determine the magnitude and frequency of thermal cycling on receiver tubes and panels using actual or realistic cloud shapes instead of averaged DNI values–which may undercount the total cycle number. This work may also prevent unnecessary plant shutdowns due to overly precautionary control strategies and characterize the relative impact of various cloud types on receiver life. We plan to eventually integrate this methodology into the System Advisor Model (SAM) to improve performance model accuracy during periods of cloudiness. In this paper, we demonstrate generating DNI maps and mapping them to a solar field in CoPylot using 10 m resolution data from publicly available Sentinel-2 satellite data over the Crescent Dunes plant.

Mullin, Matthew↗
DOE OSTI2020

Geothermal Deep Direct Use for Turbine Inlet Cooling in East Texas

The National Renewable Energy Laboratory (NREL), the Southern Methodist University Geothermal Laboratory (SMU), Eastman Chemical (Longview, TX), and TAS (Houston, TX) evaluated the feasibility of using geothermal heat to improve the performance of a natural-gas power plant in East Texas. The area of interest is the Eastman Chemical plant in Longview, Texas, which is on the northwestern margin of a geologic region known as the Sabine Uplift. The feasibility study focused on determining the potential for accessing a subsurface hot-water geothermal resource within a 10-km radius of the site to provide thermal energy for absorption chillers. Wells within a 20-km radius are included for broader geological comparison to determine the heat flow, temperature-at-depth, field porosity and permeability. The lithologies of most interest are the Lower Cretaceous Trinity Group and Upper Jurassic Cotton Valley Group. The deeper Cotton Valley formations are hotter (averaging 117 to 130°C), yet permeability and porosity are low. The shallower Trinity Group contains more variability in permeability and porosity and lower temperatures averaging about 98 to 117°C. The shallower formations are considered despite the lower temperature because of increased ability to produce larger volumes of water and extract enough heat before reinjection. The complete SMU analysis is available in the National Geothermal Data System (NGDS). Tapping such deep geothermal sources for direct heating (as opposed to power generation) is known as geothermal deep direct use (DDU). Geothermal DDU has potential across a wide swath of the United States but is underutilized due to challenging project economics associated with developing a deep geothermal resource for what are typically small-scale, variable-demand projects. This project examines the feasibility of geothermal energy integration in a natural-gas combined cycle power station in East Texas. The DDU resource is tapped to drive absorption chillers (24/7) for production of chilled water at 5-10°C (41-50°F). This chilled water is stored until needed, which allows for continuous operation with a relatively small-capacity geothermal/absorption chiller system. When conditions are favorable, the chilled water is dispatched to cool the air entering the compressor stage of a gas combustion turbine. This process, known as turbine inlet cooling (TIC), boosts power production during periods of high temperature and high-power demand. Such systems can enhance grid reliability and reduce the cost for peak-demand power. A simulation model of the power plant was developed in IPSEpro software and validated against operational data from the plant. This model allowed the team to estimate the additional power that could be produced by applying TIC under different operating and ambient conditions. Absorption chiller performance was estimated from vendor sources to determine the production rate of chilled water from the geothermal resource. Geothermal drilling and development costs were estimated using NREL's GEOPHIRES 2.0. The expected lower drilling costs in this region led to an estimated cost of geothermal heat of about $4/MMBtu (1.4 cents/kWh t ). The estimated cost for the absorption chillers and TIC hardware were obtained from literature sources and project partners. Hourly data were obtained for weather, natural gas and electricity prices, and plant operating state for 2017, which served as a representative year. NREL estimated the capital cost, operating cost, and additional electricity production and revenue for different combinations of geothermal capacity, chiller capacity, and water storage-tank size. The analysis drove toward smaller geothermal and chiller systems to reduce equipment cost. A relatively low-cost water storage tank accumulated the near-continuous chilled water output for later use when TIC was most valued.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

The Rapid X-Ray and UV Evolution of ASASSN-14ko

Abstract ASASSN-14ko is a recently discovered periodically flaring transient at the center of the active galactic nucleus (AGN) ESO 253−G003 with a slowly decreasing period. Here, we show that the flares originate from the northern, brighter nucleus in this dual-AGN, post-merger system. The light curves for the two flares that occurred in 2020 May and September are nearly identical over all wavelengths. For both events, Swift observations showed that the UV and optical wavelengths brightened in unison. The effective temperature of the UV/optical emission rises and falls with the increase and subsequent decline in the luminosity. The X-ray flux, by contrast, first rapidly drops over ∼2.6 days, rises for ∼5.8 days, drops again over ∼4.3 days, and then recovers. The X-ray spectral evolution of the two flares differ, however. During the 2020 May peak the spectrum softened with increases in the X-ray luminosity, while we observed the reverse for the 2020 September peak. We found a small change in the period derivative, which seems to indicate that the system does not have a static period derivative and there is some stochasticity in its evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Control Oriented Model of Cabin-HVAC System in a Long-Haul Trucks for Energy Management Applications

Super Truck II is a 48V mild hybrid class 8 truck with an all auxiliary loads powered purely by the battery pack. Electric Heating Ventilation and Air Conditioning (HVAC) load is the most prominent battery load during the hotel period, when the truck driver is resting inside the sleeper. For the PACCAR Super Truck II (ST-II) project a 48 V battery system provides the required power during the hotel period. A cabin-HVAC model estimates the electric load on the 48V battery system, allowing the control system to implement an efficient energy management strategy that avoids engine idling during the hotel period. The thermal model accounts for the sun load due to the time of day and the geographic location of the truck during the hotel period. The cabin-HVAC model has two parts. First, a grey box model with two heat exchangers (Condenser and Evaporator) working in unison with refrigerant mass flow rate as an input and HVAC load as an output. Second, a two-node cabin model formulated to estimate the cabin temperature as a function of the Global Horizontal Irradiance (GHI), HVAC load and ambient temperature. The models are calibrated using experimental cabin-HVAC system data as for long-haul class 8 truck (e.g. ST-II). Here, the model simulations show that the overall Root Mean Square Error (RMSE) value of 0.4°C between the experimental and simulated cabin temperature.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Dominant Wave Energy Systems and Conditional Wave Resource Characterization for Coastal Waters of the United States

Opportunities and constraints for wave energy conversion technologies and projects are evaluated by identifying and characterizing the dominant wave energy systems for United States (US) coastal waters using marginal and joint distributions of the wave energy in terms of the peak period, wave direction, and month. These distributions are computed using partitioned wave parameters generated from a 30 year WaveWatch III model hindcast, and regionally averaged to identify the dominant wave systems contributing to the total annual available energy ( A A E ) for eleven distinct US wave energy climate regions. These dominant wave systems are linked to the wind systems driving their generation and propagation. In addition, conditional resource parameters characterizing peak period spread, directional spread, and seasonal variability, which consider dependencies of the peak period, direction, and month, are introduced to augment characterization methods recommended by international standards. These conditional resource parameters reveal information that supports project planning, conceptual design, and operation and maintenance. The present study shows that wave energy resources for the United States are dominated by long-period North Pacific swells (Alaska, West Coast, Hawaii), short-period trade winds and nor’easter swells (East Coast, Puerto Rico), and wind seas (Gulf of Mexico). Seasonality, peak period spread, and directional spread of these dominant wave systems are characterized to assess regional opportunities and constraints for wave energy conversion technologies targeting the dominant wave systems.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

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

Strong Correspondence in Evapotranspiration and Carbon Dioxide Fluxes Between Different Eddy Covariance Systems Enables Quantification of Landscape Heterogeneity in Dryland Fluxes

Abstract The eddy covariance method is widely used to investigate fluxes of energy, water, and carbon dioxide at landscape scales, providing important information on how ecological systems function. Flux measurements quantify ecosystem responses to environmental perturbations and management strategies, including nature‐based climate‐change mitigation measures. However, due to the high cost of conventional instrumentation, most eddy covariance studies employ a single system, limiting spatial representation to the flux footprint. Insufficient replication may be limiting our understanding of ecosystem behavior. To address this limitation, we deployed eight lower‐cost eddy covariance systems in two clusters around two conventional eddy covariance systems in the Chihuahuan Desert of North America for a period of 2 years. These dryland settings characterized by large temperature variations and relatively low carbon dioxide fluxes represented a challenging setting for eddy covariance. We found very good closure of energy and water balance across all systems (within ±9% of unity). We found very good correspondence between the lower‐cost and conventional systems' fluxes of sensible heat (with concordance correlation coefficient (CCC) of ≥0.87), latent energy (evapotranspiration; CCC ≥ 0.89), and useful correspondence in the net ecosystem exchange ((NEE); with CCC ≥ 0.4) at the daily temporal resolution. Relative to the conventional systems, the low‐frequency systems were characterized by a higher level of random error, particularly in the NEE fluxes. Lower‐cost systems can enable wider deployment affording better replication and sampling of spatiotemporal variability at the expense of greater measurement noise that might be limiting for certain applications. Replicated eddy covariance observations may be useful when addressing gaps in the existing monitoring of critical and underrepresented ecosystems and for measuring areas larger than a single flux footprint.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Automatic Search of Cataclysmic Variables Based on LightGBM in LAMOST-DR7

The search for special and rare celestial objects has always played an important role in astronomy. Cataclysmic Variables (CVs) are special and rare binary systems with accretion disks. Most CVs are in the quiescent period, and their spectra have the emission lines of Balmer series, HeI, and HeII. A few CVs in the outburst period have the absorption lines of Balmer series. Owing to the scarcity of numbers, expanding the spectral data of CVs is of positive significance for studying the formation of accretion disks and the evolution of binary star system models. At present, the research for astronomical spectra has entered the era of Big Data. The Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST) has produced more than tens of millions of spectral data. the latest released LAMOST-DR7 includes 10.6 million low-resolution spectral data in 4926 sky regions, providing ideal data support for searching CV candidates. To process and analyze the massive amounts of spectral data, this study employed the Light Gradient Boosting Machine (LightGBM) algorithm, which is based on the ensemble tree model to automatically conduct the search in LAMOST-DR7. Finally, 225 CV candidates were found and four new CV candidates were verified by SIMBAD and published catalogs. This study also built the Gradient Boosting Decision Tree (GBDT), Adaptive Boosting (AdaBoost), and eXtreme Gradient Boosting (XGBoost) models and used Accuracy, Precision, Recall, the F1-score, and the ROC curve to compare the four models comprehensively. Experimental results showed that LightGBM is more efficient. The search for CVs based on LightGBM not only enriches the existing CV spectral library, but also provides a reference for the data mining of other rare celestial objects in massive spectral data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Iterated projected position algorithm for constructing exponentially localized generalized Wannier functions for periodic and nonperiodic insulators in two dimensions and higher

Localized bases play an important role in understanding electronic structure. In periodic insulators, a natural choice of localized basis is given by the Wannier functions which depend on a choice of unitary transform known as a gauge transformation. Over the past few decades, there have been many works that have focused on optimizing the choice of the gauge so that the corresponding Wannier functions are maximally localized or reflect some symmetry of the underlying system. In this work, we consider fully nonperiodic materials where the usual Wannier functions are not well defined and gauge optimization is impractical. To tackle the problem of calculating exponentially localized generalized Wannier functions in both periodic and nonperiodic systems, we discuss the ‘iterated projected position (IPP)” algorithm. Here, the IPP algorithm is based on matrix diagonalization and therefore unlike optimization-based approaches, it does not require initialization and cannot get stuck at a local minimum. Furthermore, the IPP algorithm is guaranteed by a rigorous analysis to produce exponentially localized functions under certain mild assumptions. We numerically demonstrate that the IPP algorithm can be used to calculate exponentially localized bases for the Haldane model, the Kane-Mele model (in both Z 2 invariant even and Z 2 invariant odd phases), and the p x + ip y model on a quasicrystal lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Two Young Planetary Systems around Field Stars with Ages between 20 and 320 Myr from TESS

Planets around young stars trace the early evolution of planetary systems. We report the discovery and validation of two planetary systems with ages ≲300 Myr from observations by the Transiting Exoplanet Survey Satellite (TESS). The 40--320 Myr old G star TOI-251 hosts a 2.74{sub −0.18}{sup +0.18} R{sub ⊕} mini-Neptune with a 4.94 day period. The 20--160 Myr old K star TOI-942 hosts a system of inflated Neptune-sized planets, with TOI-942b orbiting in a period of 4.32 days with a radius of 4.81{sub −0.20}{sup +0.20} R{sub ⊕} and TOI-942c orbiting in a period of 10.16 days with a radius of 5.79{sub −0.18}{sup +0.19} R{sub ⊕}. Though we cannot place either host star into a known stellar association or cluster, we can estimate their ages via their photometric and spectroscopic properties. Both stars exhibit significant photometric variability due to spot modulation, with measured rotation periods of ∼3.5 days. These stars also exhibit significant chromospheric activity, with age estimates from the chromospheric calcium emission lines and X-ray fluxes matching that estimated from gyrochronology. Both stars also exhibit significant lithium absorption, similar in equivalent width to well-characterized young cluster members. TESS has the potential to deliver a population of young planet-bearing field stars, contributing significantly to tracing the properties of planets as a function of their age.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Digital Real-Time Simulation and Power Quality Analysis of a Hydrogen-Generating Nuclear-Renewable Integrated Energy System

This paper investigates the challenges and solutions associated with integrating a hydrogen-generating nuclear-renewable integrated energy system (NR-IES) under a transactive energy framework. The proposed system directs excess nuclear power to hydrogen production during periods of low grid demand while utilizing renewables to maintain grid stability. Using digital real-time simulation (DRTS) in the Typhoon HIL 404 model, the dynamic interactions between nuclear power plants, electrolyzers, and power grids are analyzed to mitigate issues such as harmonic distortion, power quality degradation, and low power factor caused by large non-linear loads. A three-phase power conversion system is modeled using the Typhoon HIL 404 model and includes a generator, a variable load, an electrolyzer, and power filters. Active harmonic filters (AHFs) and hybrid active power filters (HAPFs) are implemented to address harmonic mitigation and reactive power compensation. The results reveal that the HAPF topology effectively balances cost efficiency and performance and significantly reduces active filter current requirements compared to AHF-only systems. During maximum electrolyzer operation at 4 MW, the grid frequency dropped below 59.3 Hz without filtering; however, the implementation of power filters successfully restored the frequency to 59.9 Hz, demonstrating its effectiveness in maintaining grid stability. Future work will focus on integrating a deep reinforcement learning (DRL) framework with real-time simulation and optimizing real-time power dispatch, thus enabling a scalable, efficient NR-IES for sustainable energy markets.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Real-time monitoring and prediction of water quality parameters and algae concentrations using microbial potentiometric sensor signals and machine learning tools

We report the overarching hypothesis of this study was that temporal microbial potentiometric sensor (MPS) signal patterns could be used to predict changes in commonly monitored water quality parameters by using artificial intelligence/machine learning tools. To test this hypothesis, the study first examines a proof of concept by correlating between MPS's signals and high algae concentrations in an algal cultivation pond. Then, the study expanded upon these findings and examined if multiple water quality parameters could be predicted in real surface waters, like irrigation canals. Signals generated between the MPS sensors and other water quality sensors maintained by an Arizona utility company, including algae and chlorophyll, were collected in real time at time intervals of 30 min over a period of 9 months. Data from the MPS system and data collected by the utility company were used to train the ML/AI algorithms and compare the predicted with actual water quality parameters and algae concentrations. Based on the composite signal obtained from the MPS, the ML/AI was used to predict the canal surface water's turbidity, conductivity, chlorophyll, and blue-green algae (BGA), dissolved oxygen (DO), and pH, and predicted values were compared to the measured values. Initial testing in the algal cultivation pond revealed a strong linear correlation (R 2 = 0.87) between mixed liquor suspended solids (MLSS) and the MPSs' composite signals. The Normalized Root Mean Square Error (NRMSE) between the predicted values and measured values were <6.5%, except for the DO, which was 10.45%. The results demonstrate the usefulness of MPSs to predict key surface water quality parameters through a single composite signal, when the ML/AI tools are used conjunctively to disaggregate these signal components. The maintenance-free MPS offers a novel and cost-effective approach to monitor numerous water quality parameters at once with relatively high accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

In situ nanoscale imaging of moiré superlattices in twisted van der Waals heterostructures

Abstract Direct visualization of nanometer-scale properties of moiré superlattices in van der Waals heterostructure devices is a critically needed diagnostic tool for study of the electronic and optical phenomena induced by the periodic variation of atomic structure in these complex systems. Conventional imaging methods are destructive and insensitive to the buried device geometries, preventing practical inspection. Here we report a versatile scanning probe microscopy employing infrared light for imaging moiré superlattices of twisted bilayers graphene encapsulated by hexagonal boron nitride. We map the pattern using the scattering dynamics of phonon polaritons launched in hexagonal boron nitride capping layers via its interaction with the buried moiré superlattices. We explore the origin of the double-line features imaged and show the mechanism of the underlying effective phase change of the phonon polariton reflectance at domain walls. The nano-imaging tool developed provides a non-destructive analytical approach to elucidate the complex physics of moiré engineered heterostructures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Activated vortex lattice transition in a superconductor with combined sixfold and twelvefold anisotropic interactions

Abstract Numerical simulations are used to examine the transition dynamics between metastable and ground state vortex lattice phases in a system with combined sixfold and twelvefold contributions to the vortex–vortex interactions. The system is initially annealed using a twelvefold anisotropy, yielding domains of two different orientations and separated by grain boundaries. The vortex–vortex interaction is then suddenly changed to a sixfold anisotropy, rendering the twelvefold state metastable. Applying a drive that mimics an oscillating magnetic field causes the metastable state to decay, indicated by the structure factor that evolves from twelve to six peaks. The results fit the behavior seen in recent small-angle neutron scattering studies of the vortex lattice in MgB 2 . At higher drive amplitudes, the decay exhibits a two step process in which the initial fast decrease is followed by a slower regime where avalanches or bursts are correlated with dislocation annihilation events. The results are compared to other types of metastable systems with quenched disorder that decay under a periodic external drive.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Size, shade, or shape? The contribution of galaxies of different types to the star formation history of the Universe from SDSS-IV MaNGA

ABSTRACT By fitting stellar populations to the fourth generation of the Sloan Digital Sky Survey (SDSS-IV) Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey observations of ∼7000 suitably weighted individual galaxies, we reconstruct the star formation history of the Universe, which we find to be in reasonable agreement with previous studies. Dividing the galaxies by their present-day stellar mass, we demonstrate the downsizing phenomenon, whereby the more massive galaxies hosted the most star formation at earlier times. Further dividing the galaxy sample by colour and morphology, we find that a galaxy’s present-day colour tells us more about its historical contribution to the cosmic star formation history than its current morphology. We show that downsizing effects are greatest among galaxies currently in the blue cloud, but that the level of downsizing in galaxies of different morphologies depends quite sensitively on the morphological classification used, due largely to the difficulty in classifying the smaller low-mass galaxies from their ground-based images. Nevertheless, we find agreement that among galaxies with stellar masses $M_{\star } \gt 6\times 10^{9}\, \mathrm{ M}_{\odot }$, downsizing is most significant in spirals. However, there are complicating factors. For example, for more massive galaxies, we find that colour and morphology are predictors of the past star formation over a longer time-scale than in less massive systems. Presumably this effect is reflecting the longer period of evolution required to alter these larger galaxies’ physical properties, but shows that conclusions based on any single property do not tell the full story.

Peterken, Thomas↗
DOE OSTI2022

Four new deeply eclipsing white dwarfs in Zwicky Transient Facility

ABSTRACT We present the results of a search for deeply eclipsing white dwarfs in the Zwicky Transient Facility (ZTF) Data Release 4 (DR4). We identify nine deeply eclipsing white dwarf candidates, four of which we followed up with high-cadence photometry and spectroscopy. Three of these systems show total eclipses in the ZTF data and our follow-up Apache Point Observatory 3.5 m telescope observations. Even though the eclipse duration is consistent with sub-stellar companions, our analysis shows that all four systems contain a white dwarf with low-mass stellar companions of ∼0.1 M⊙. We provide mass and radius constraints for both stars in each system based on our photometric and spectroscopic fitting. Finally, we present a list of 41 additional eclipsing WD+M candidates identified in a preliminary search of ZTF DR7, including 12 previously studied systems. We identify two new candidate short-period, eclipsing, white dwarf–brown dwarf binaries within our sample of 41 WD+M candidates based on Pan-STARRS colours.

Kosakowski, A. (ORCID:0000000298781647)↗
DOE OSTI2022

Commensuration effects on skyrmion Hall angle and drag for manipulation of skyrmions on two-dimensional periodic substrates

Here, we examine the dynamics of an individually driven skyrmion moving through a background lattice of skyrmions coupled to a 2D periodic substrate as we vary the ratio of the number of skyrmions to the number of pinning sites across commensurate and incommensurate conditions. As the skyrmion density increases, the skyrmion Hall angle is nonmonotonic, dropping to low or zero values in commensurate states and rising to an enhanced value in incommensurate states. Under commensuration, the driven skyrmion is channeled by a symmetry direction of the pinning array and exhibits an increased velocity. At fillings for which the skyrmion Hall angle is zero, the velocity has a narrow band noise signature, while for incommensurate fillings, the skyrmion motion is disordered and the velocity noise is broad band. Under commensurate conditions, multi-step depinning transitions appear and the skyrmion Hall angle is zero at low drives but becomes finite at higher drives, while at incommensurate fillings there is only a single depinning transition. As the gyrotropic component of the skyrmion dynamics, called the Magnus force, increases, peaks in the velocity that appear in commensurate regimes cross over to dips, and new types of directional locking effects can arise in which the skyrmion travels along other symmetry directions of the background lattice. At large Magnus forces, and particularly at commensurate fillings, the driven skyrmion can experience a velocity boost in which the skyrmion moves faster than the applied drive due to the alignment of the Magnus-induced velocity with the driving direction. In some cases, an increase of the Magnus force can produce regimes of enhanced pinning when the skyrmion is forced to move along a nonsymmetry direction of the periodic pinning array. This is in contrast to systems with random pinning, where increasing the Magnus force generally reduces the pinning effect. We demonstrate these dynamics for both square and triangular substrates and map out the different regimes as a function of filling fraction, pinning force, and the strength of the Magnus force in a series of dynamic phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Numerical heating in particle-in-cell simulations with Monte Carlo binary collisions

The binary Monte Carlo (MC) collision algorithm is a standard and robust method to include binary Coulomb collision effects in particle-in-cell (PIC) simulations of plasmas. Here we show that the coupling between PIC and MC algorithms can give rise to (nonphysical) numerical heating of the system that significantly exceeds that observed when these algorithms operate independently. We argue that this deleterious effect results from an inconsistency between the particle motion associated with MC collisions and the work performed by the collective electromagnetic field on the PIC grid. This inconsistency manifests as the (artificial) stochastic production of electromagnetic energy, which ultimately heats the plasma particles. Here, the MC-induced numerical heating can significantly impact the evolution of the simulated system for long simulation times (≳10 3 collision periods, for typical numerical parameters). We describe the source of the MC-induced numerical heating analytically and discuss strategies to minimize it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗