Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modeling kinetic effects of charged vacancies on electromechanical responses of ferroelectrics: Rayleighian approach

Understanding the time-dependent effects of charged vacancies on the electromechanical responses of materials is at the forefront of research for designing materials exhibiting metal-insulator transitions and memristive behavior. A Rayleighian approach is used to develop a model for studying the nonlinear kinetics of the reaction leading to generation of vacancies and electrons via the dissociation of vacancy-electron pairs. Also, diffusion and elastic effects of charged vacancies are considered to model polarization-electric potential and strain-electric potential hysteresis loops. The model captures multiphysics phenomena by introducing couplings among polarization, the electric potential, stress, strain, and concentrations of charged (multivalent) vacancies and electrons (treated as classical negatively charged particles), where the concentrations can vary due to association-dissociation reactions. A derivation of coupled time-dependent equations based on the Rayleighian approach is presented. Three limiting cases of the governing equations are considered, highlighting the effects of (1) nonlinear reaction kinetics on the generation of charged vacancies and electrons, (2) Vegard's law (i.e., the concentration-dependent local strain) on asymmetric strain-electric potential relations, and (3) coupling between a fast component and the slow component of the net polarization on the polarization-electric-field relations. The Rayleighian approach discussed in this work should pave the way for developing a multiscale modeling framework in a thermodynamically consistent manner while capturing multiphysics phenomena in ferroelectric materials. Published by the American Physical Society 2025

Kumar, Rajeev (ORCID:0000000194943488)↗
NASA NTRS1982

The kinetic chemistry of dense interstellar clouds

A model of the time-dependent chemistry of dense interstellar clouds is formulated to study the dominant chemical processes in carbon and oxygen isotope fractionation, the formation of nitrogen-containing molecules, and the evolution of product molecules as a function of cloud density and temperature. The abundances of the dominant isotopes of the carbon- and oxygen-bearing molecules are calculated. The chemical abundances are found to be quite sensitive to electron concentration since the electron concentration determines the ratio of H3(+) to He(+), and the electron density is strongly influenced by the metals abundance. For typical metal abundances and for H2 cloud density not less than 10,000 molecules/cu cm, nearly all carbon exists as CO at late cloud ages. At high cloud density, many aspects of the chemistry are strongly time dependent. Finally, model calculations agree well with abundances deduced from observations of molecular line emission in cold dense clouds.

Graedel, T. E.↗
NASA NTRS2001

High Temperature, Slow Strain Rate Forging of Advanced Disk Alloy ME3

The advanced disk alloy ME3 was designed in the HSR/EPM disk program to have extended durability at 1150 to 1250 F in large disks. This was achieved by designing a disk alloy and process producing balanced monotonic, cyclic, and time-dependent mechanical properties. combined with robust processing and manufacturing characteristics. The resulting baseline alloy, processing, and supersolvus heat treatment produces a uniform, relatively fine mean grain size of about ASTM 7, with as-large-as (ALA) grain size of about ASTM 3. There is a long term need for disks with higher rim temperature capabilities than 1250 F. This would allow higher compressor exit (T3) temperatures and allow the full utilization of advanced combustor and airfoil concepts under development. Several approaches are being studied that modify the processing and chemistry of ME3, to possibly improve high temperature properties. Promising approaches would be applied to subscale material, for screening the resulting mechanical properties at these high temperatures. n obvious path traditionally employed to improve the high temperature and time-dependent capabilities of disk alloys is to coarsen the grain size. A coarser grain size than ASTM 7 could potentially be achieved by varying the forging conditions and supersolvus heat treatment. The objective of this study was to perform forging and heat treatment experiments ("thermomechanical processing experiments") on small compression test specimens of the baseline ME3 composition, to identify a viable forging process allowing significantly coarser grain size targeted at ASTM 3-5, than that of the baseline, ASTM 7.

Gabb, Timothy P.↗
DOE OSTI2024

Chiral kinematic theory and converse vortical effects

Response theories in condensed matter typically describe the response of an electron fluid to external electromagnetic fields, while perturbations on neutral particles are often designed to mimic such fields. Here, we study the response of fermions to a space-time-dependent velocity field, thereby sidestepping the issue of gauge charge. First, we use a semiclassical chiral kinematic theory to obtain the local density of current and extract the orbital magnetization. The theory immediately predicts a "converse vortical effect," defined as an orbital magnetization driven by linear velocity. It receives contributions from magnetic moments on the Fermi surface and the Berry curvature of the occupied bands. Then, transcending semiclassics via a complementary Kubo formalism reveals that the uniform limit of a clean system receives only the Berry curvature contribution while other limits sense the Fermi surface magnetic moments too. We propose CoSi as a candidate material and suggest magnetometry of a sample under a thermal gradient to detect the effect. Overall, our study sheds light on the effects of a space-time-dependent velocity field on electron fluids and paves the way for exploring quantum materials using new probes and perturbations.

Chen, Kai↗
NASA NTRS2003

CARES/LIFE Ceramics Analysis and Reliability Evaluation of Structures Life Prediction Program

This manual describes the Ceramics Analysis and Reliability Evaluation of Structures Life Prediction (CARES/LIFE) computer program. The program calculates the time-dependent reliability of monolithic ceramic components subjected to thermomechanical and/or proof test loading. CARES/LIFE is an extension of the CARES (Ceramic Analysis and Reliability Evaluation of Structures) computer program. The program uses results from MSC/NASTRAN, ABAQUS, and ANSYS finite element analysis programs to evaluate component reliability due to inherent surface and/or volume type flaws. CARES/LIFE accounts for the phenomenon of subcritical crack growth (SCG) by utilizing the power law, Paris law, or Walker law. The two-parameter Weibull cumulative distribution function is used to characterize the variation in component strength. The effects of multiaxial stresses are modeled by using either the principle of independent action (PIA), the Weibull normal stress averaging method (NSA), or the Batdorf theory. Inert strength and fatigue parameters are estimated from rupture strength data of naturally flawed specimens loaded in static, dynamic, or cyclic fatigue. The probabilistic time-dependent theories used in CARES/LIFE, along with the input and output for CARES/LIFE, are described. Example problems to demonstrate various features of the program are also included.

Nemeth, Noel N.↗
NASA NTRS2010

Winds from Luminous Late-Type Stars: II. Broadband Frequency Distribution of Alfven Waves

We present the numerical simulations of winds from evolved giant stars using a fully non-linear, time dependent 2.5-dimensional magnetohydrodynamic (MHD) code. This study extends our previous fully non-linear MHD wind simulations to include a broadband frequency spectrum of Alfven waves that drive winds from red giant stars. We calculated four Alfven wind models that cover the whole range of Alfven wave frequency spectrum to characterize the role of freely propagated and reflected Alfven waves in the gravitationally stratified atmosphere of a late-type giant star. Our simulations demonstrate that, unlike linear Alfven wave-driven wind models, a stellar wind model based on plasma acceleration due to broadband non-linear Alfven waves, can consistently reproduce the wide range of observed radial velocity profiles of the winds, their terminal velocities and the observed mass loss rates. Comparison of the calculated mass loss rates with the empirically determined mass loss rate for alpha Tau suggests an anisotropic and time-dependent nature of stellar winds from evolved giants.

Airapetian, V.↗
NASA NTRS1998

New Approach to Diagnosing Properties of Protoplanetary Disks

In this paper we suggest that subjecting the observationally derived properties of protoplanetary disks to the evolutionary interpretation yields new insights into the working of those disks, and offers valuable constraints on their models. We propose that the global properties of individual disks, such as their accretion rates and disk masses, sorted by the mass of the central star, can be indexed by the age of the star to simulate the evolution of a single disk. Using data from published surveys of T Tauri stars, we show that accretion rate data, and disk mass data for the lowest mass stars, form well-defined evolutionary tracks. The higher mass stars show a definitive negative correlation between accretion rates and star ages. We use the time-dependent alpha-disk model of the viscous protoplanetary disk to link the theory to observations. The data are consistent with the standard theoretical paradigm, but not with the layered accretion model. The best fits to the data are obtained for the standard models that start with disks that are about one-third of the mass of the central star and have their angular momenta, j, and alpha-coefficients linked by the relationship j varies as Solar mass(exp 3/2)alpha(exp 1/3). The proportionality constant in this relationship, when derived from the accretion rate data, differs from the constant derived from the disk mass data. We argue that the accretion rate data are more reliable. Taking into account typical values of the specific angular momentum of disk-forming matter, we obtain alpha is greater than or equal to 10(exp -2). A complete time-dependent standard disk model, built on the parameters determined from the best-fit procedure, is presented. Such a model constitutes a good point of departure for various theoretical studies aimed at the issue of formation of planetary systems and the character of protoplanetary disks.

Stepinski, Tomasz F.↗
DOE OSTI2026

Performing Numerical Analysis of Cybersecurity Options Using Dynamic Risk Analysis Tool EMRALD

Cyberattacks can have many different attack paths, durations, and goals. There are also many different mitigation options involving hardware, software, and/or humans. Considering a cyber threat should involve defense-in-depth methods and a quantitative or numerical evaluation of overall effectiveness against dynamic, time-dependent attacks to make cost and risk-informed decisions. Typical cyberattack modeling methods only provide a qualitative evaluation. The main areas of cybersecurity are confidentiality, integrity, and availability. For companies with cyber-physical systems such as advanced nuclear reactors, cyber-related safety is a requirement set by North American Electric Reliability and the U.S. Nuclear Regulatory Commission. They are also concerned about availability or reliability as a business case. As cyber threats are evolving to a business-for-hire structure, more attacks may focus on disrupting business success and reliability, causing financial and economic stability risk. Companies want to know business reliability and recovery from those threats, and that requires modeling physical behavior of the targets. Dynamic-state-based and Markov-based modeling provides a method for better cyber scenario modeling with different tools having issues such as state-base explosion. Dynamic modeling enables time and conditional features not found in other numerical evaluation methods. EMRALD (Event Modeling Risk Assessment using Lined Diagrams) is a dynamic risk analysis modeling and simulation tool and has features that reduce modeling issues. It has been used to model different time-dependent events including plant behavior and operator procedures. As a general modeling tool, EMRALD can also be used to model cyberattack scenarios with varying mitigation options and quantify effectiveness, producing numerical data for risk-informed decisions. This paper uses EMRALD to demonstrate that dynamic numerical risk analysis can be used for cyber threat modeling to provide insights for design decision-making and optimize defense strategies. Keywords: cyber modeling; cyber-physical systems; numerical cyber modeling

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS2014

Terminal Area Simulation System User's Guide - Version 10.0

The Terminal Area Simulation System (TASS) is a three-dimensional, time-dependent, large eddy simulation model that has been developed for studies of wake vortex and weather hazards to aviation, along with other atmospheric turbulence, and cloud-scale weather phenomenology. This document describes the source code for TASS version 10.0 and provides users with needed documentation to run the model. The source code is programed in Fortran language and is formulated to take advantage of vector and efficient multi-processor scaling for execution on massively-parallel supercomputer clusters. The code contains different initialization modules allowing the study of aircraft wake vortex interaction with the atmosphere and ground, atmospheric turbulence, atmospheric boundary layers, precipitating convective clouds, hail storms, gust fronts, microburst windshear, supercell and mesoscale convective systems, tornadic storms, and ring vortices. The model is able to operate in either two- or three-dimensions with equations numerically formulated on a Cartesian grid. The primary output from the TASS is time-dependent domain fields generated by the prognostic equations and diagnosed variables. This document will enable a user to understand the general logic of TASS, and will show how to configure and initialize the model domain. Also described are the formats of the input and output files, as well as the parameters that control the input and output.

Switzer, George F.↗
NASA NTRS2019

Exploration of Failure and Potential Damage Markers in Ti-6Al-4V

To improve prognostics for predicting component life, coupon tests were conducted on Ti-6Al-4V to investigate the development of damage. Various test types were employed (tension, creep, stress relaxation, and fatigue) under a wide range of temperatures and loading rates to engage various amounts of time-dependent and damage behavior. Both stiffness degradation and Poisson’s ratio were monitored to see if they were useful in indicating the initiation and accumulation of damage and subsequently the onset of failure. Stiffness exhibited large decreases with increasing strain when calculated using engineering stress and strain. However, minimal stiffness change was observed when using the true-stress and true- strain values. Detailed microscopy and image analysis were performed to document the actual physical damage that occurred within the various specimens. Two types of damage were observed: shallow surface cracks and internal pores. The severity of damage was greater in the region of localized necking and diminished further away from the neck. The area fraction of pores was higher in the necked region, yet even there it did not exceed more than a few percent and, in most cases, was significantly lower. It was observed that conditions that induced more time-dependent deformation (e.g., slower strain rate testing) resulted in a larger number of pores and larger pore sizes. The pores were always found to occur in the softer β-phase. The severity of observed damage was inconsistent with the large drop in engineering stiffness, but was more aligned with stiffness changes from the true-stress and -strain values. Poisson’s ratio was monitored on every test and showed promise as a damage indicator. It was found to be more sensitive than either the ultimate tensile strength or the inception of tertiary creep. However, calculation of Poisson’s ratio was found to be more involved than expected because of anisotropy and the presence of strain gradients along the specimen at large applied strains. First, the material had a moderate texture and resulted in anisotropy of the strains. Second, the evolution of Poisson’s ratio as a function of axial strain showed strange and unexpected behaviors, which necessitated more extensive analyses of its development. This led to an investigation on the use of small versus large strain approaches and the importance of making consistent measurements; that is, axial and transverse strain measurements over the same volume of material. Additionally, finite element analysis of a tensile sample indicated that non- uniform strain occurred in the sample prior to attaining the maximum strength. This, along with the location at which transverse strain was measured, had a significant influence on the resulting values of Poisson’s ratio. Based on this, lessons learned are offered for future testing.

Ti-6Al-4V↗
DOE OSTI2021

Damage development in neutron-irradiated concrete in a test reactor: Hygro-thermal and mechanical simulations

This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Higher Ion Selectivity with Lower Energy Usage Promoted by Electro-osmotic Flow in the Transport through Conical Nanopores

The tradeoff between selectivity and throughput presents fundamental challenges to improve desalination and charge storage, salinity gradient-based energy harvesting, memory device/circuit development, and so forth. The well-known ion-current rectification and the recently resolved time-dependent transport hysteresis in conical nanopores or asymmetric nanointerfaces offer new opportunities for the selective transport of matter. This report shows that electro-osmotic flow (EOF) is an overlooked factor that increases ion selectivity while maintaining enhanced transport throughput in rectified nanoscale electrokinetic transport. The increased selectivity originates primarily from the suppression of anion flux by the fluid flow in the opposite direction under the applied electrical field. By solving the Poisson and Nernst–Planck (PNP) equations without and with coupled Navier–Stokes (PNP–NS) equations, the EOF effects on cation and anion transport are unequivocally revealed in asymmetric nanopipettes. The flux of cations and anions as well as the transference number and flow velocity are elucidated using the models and boundary conditions validated by previous experiments. A dimensionless parameter, radius over the Debye length, reveals optimal ion selectivity and energy cost at intermediate ion concentrations and nanopore sizes, up to hundreds of millimolars and tens of nanometers under time-dependent potential stimulus. Furthermore, the fundamental insights into EOF at nanointerfaces suggest new routes/strategies for better separation, analysis, and energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2015

Generalized Fluid System Simulation Program (GFSSP) - Version 6

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, flow control valves and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. Users can introduce new physics, non-linear and time-dependent boundary conditions through user-subroutine.

Majumdar, Alok↗
NASA NTRS1998

The Impact of FU Orionis Outbursts and the Solar Nebula

Protostellar systems are variable on many timescales. One of the most dramatic forms of variability known to occur in low mass stellar systems is the FU Orionis outburst (Herbig 1977). Throughout a typical outburst lasting several decades, system luminosities may be a hundred times what is typical of the quiesent state. FU Orionis outburst events are thought to have significant impact on the thermal structure of the protosolar nebula; their existence has been used to explain features in the meteoritic record from thermally induced homogenization to the formation of chondrules. Until recently, the magnitude of the likely effect from such outbursts has been largely speculative due to the lack of a detailed understanding of the outburst mechanism. Recent numerical models (Bell\& Lin 1994) have demonstrated the viability of the observational hypothesis (Hartmann\& Kenyon 1985) that the radiation observed during outburst is emitted by a luminous circumstellar disk transporting mass at a thousand times the quiesent rate. Light curves and color and line width evolution observed in FU Orionis systems are naturally explained by time dependent outbursting model disks (Bell et al. 1995). The radial temperature structure and shape of the disk during outburst derived from these models may be used to calculate the outburst's expected impact on primitive material at various radii throughout the disk. In this review, we will begin by discussing what is known about the FU Orionis outburst phenomenon from recent observations and theory including statistically deduced outburst timescales and observed peak temperatures. Unless covered by another author, we will discuss the evidence which suggests that outburst radiation is emitted by a circumstellar disk rather than by the star and will briefly review the thermal instability as a mechanism for outburst. We will then report on recent work which investigates the likely heating of solar nebula material due to FU Orionis outbursts including the following effects: (1) heating of the planet forming region by direct radiation from the hot inner nebula; (2) heating by the diffuse radiation field of a coccooning envelope; and (3) time-dependent penetration of the increased luminosity from the above sources into the optically thick nebula. Some of this work is currently in progress. The potential effects on condensation and migration in the nebula and the thermal processing of solids will be evaluated.

Bell, Robbins↗
DOE OSTI2024

Modern chemical graph theory

Abstract Graph theory has a long history in chemistry. Yet as the breadth and variety of chemical data is rapidly changing, so too do graph encoding methods and analyses that yield qualitative and quantitative insights. Using illustrative cases within a basic mathematical framework, we showcase modern chemical graph theory's utility in Chemists' analysis and model development toolkit. The encoding of both experimental and simulation data is discussed at various levels of granularity of information. This is followed by a discussion of the two major classes of graph theoretical analyses: identifying connectivity patterns and partitioning methods. Measures, metrics, descriptors, and topological indices are then introduced with an emphasis upon enhancing interpretability and incorporation into physical models. Challenging data cases are described that include strategies for studying time dependence. Throughout, we incorporate recent advancements in computer science and applied mathematics that are propelling chemical graph theory into new domains of chemical study. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Structure and Mechanism > Computational Materials Science Structure and Mechanism > Molecular Structures

Leite, Leonardo S. G.↗
DOE OSTI2019

The Belle II Physics Book

The Belle II-Theory interface Platform (B2TiP) was created as a physics prospects working group of the Belle II collaboration in June 2014. It offered a platform where theorists and experimentalists work together to elucidate the potential impacts of the Belle II program, which includes a wide scope of physics topics: B physics, charm, τ, quarkonium physics, electroweak precision measurements and dark sector searches. It is composed of nine working groups (WGs), which are coordinated by teams of theory and experiment conveners: WG1) Semileptonic and leptonic B decays, WG2) Radiative and Electroweak penguins, WG3) φ1 and φ2 (time-dependent CP violation) measurements, WG4) φ3 measurement, WG5) Charmless hadronic B decay, WG6) Charm, WG7) Quarkonium(like), WG8) τ and lowmultiplicity processes, WG9) New Physics. We organised workshops twice a year from 2014 until 2016, which moved from KEK in Japan to Europe and the Americas, gathering experts in the respective fields to discuss with Belle II members. One of the goals for B2TiP was to propose so-called “golden- and silver-channels”: we asked each working group to choose among numerous possible measurements, those that would have the highest potential impact and to focus on them for the writeup. Theorists scrutinised the role of those measurements in terms of understanding the theory behind them, and estimated the theoretical uncertainties, now achievable as well as prospects for the future. For flavour physics, having tight control of hadronic uncertainties is one of the most crucial aspects in the field, and this is considered as an important criteria to determine the golden or silver-channels. Experimentalists, on the other hand, investigated the expected improvements with data from Belle II. For the channels where the errors are dominated by statistical uncertainties, or where systematic errors are reducible, the errors can decrease rapidly as more data becomes available. The impact of the upgraded performance from Belle II is a crucial element for reducing the uncertainties: we therefore include the latest available studies of the detector efficiency using Monte Carlo simulated events. We list the golden(silver)-channel table in this first chapter, as a guide for the chapters that follow. The book is not a collection of reports based on talks given at the workshops. The working group conveners endeavoured to construct a coherent document that can be used by Belle II collaborators, and others in the field of flavour physics, as a reference. There were two books of a similar type written in the past, “The BaBar book” [1] and “The Physics of the B factories” [2]. In order to avoid too much repetition with respect to those references, we refer to them as possible for introductory material. We would like to thank the section editors and contributing authors for the many stimulating discussions and their tremendous effort to bring the book together.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1995

Propagation of nuclear burning fronts on accreting neutron stars: X-ray bursts and sub-hertz noise

We identify a new regime of time dependent helium burning for high accretion rate neutron stars and suggest that this burning is the origin of the low-level luminosity variations (on timescales of 10-10(exp 4) s, designated the 'very low-frequency noise'(VLFN) by van der Klis and collaborators) always detected in the brightest accreting X-ray sources. Only two nuclear burning regimes were previously recognized. At accretion rates in excess of the Eddington limit (dot-M approximately greater than (1-3) x 10(exp -8) solar mass/yr), the accreted matter fuses steadily. At very low dot-M, the star's entire surface is rapidly (approximately less than 10 s) burned by a fast propagating convective burning front at regular intervals, giving quasi-periodic Type I X-ray bursts. We show that for the observationally interesting range of 5 x 10(exp -10) solar mass/yr approximately less than dot-M approximately less than 10(exp -8) solar mass/yr, parts of the stellar surface burn slowly. At these accretion rates, a local thermonuclear instability starts a fire which propagates horizontally at v approximately 300 cm/s. The fire propagates around the flammable surface in roughly the same time it takes to accrete enough fuel for the next instability (approximately 10(exp 3)-10(exp 4), so that only a few fires are burning at once, giving rise to large luminosity flares. Nuclear burning is always time dependent for sub-Eddington local accretion rates: a local patch undergoes a recurrent cycle, accumulation fuel for hours until it becomes thermally unstable or is 'ignited' by a nearby burning region. The global pattern of burning and the resulting luminosity are thus very dependent on how fast nuclear fires spread around the star. The nuclear burning luminosity is not uniform over the stellar surface and so may provide a handle on measuring, or constraining, the spin periods of these neutron stars.

Bildsten, Lars↗
DOE OSTI2026

Characterizing skyrmion flow phases with principal component analysis

Principal component analysis (PCA) is a powerful method that can identify patterns in large, complex data sets by constructing low-dimensional order parameters from higher-dimensional feature vectors. There are increasing efforts to use space-and-time-dependent PCA to detect transitions in nonequilibrium systems that are difficult to characterize with equilibrium methods. Here, we demonstrate that feature vectors incorporating the position and velocity information of driven skyrmions moving through random disorder permit PCA to resolve different types of disordered skyrmion motion as a function of driving force and the ratio of the Magnus force to the dissipation. Since the Magnus force creates gyroscopic motion and a finite Hall angle, skyrmions can exhibit a greater range of flow phases than what is observed in overdamped driven systems with quenched disorder. We show that in addition to identifying previously known skyrmion flow phases, PCA detects several additional phases, including different types of channel flow, moving fluids, and partially ordered states. Guided by the PCA analysis, we further characterize the disordered flow phases to elucidate the different microscopic dynamics and show that the changes in the PCA-derived order parameters can be connected to features in bulk transport measures, including the transverse and longitudinal velocity-force curves, differential conductivity, topological defect density, and changes in the skyrmion Hall angle as a function of drive. We discuss how asymmetric feature vectors can be used to improve the resolution of the PCA analysis, and how this technique can be extended to find disordered phases in other nonequilibrium systems with time-dependent dynamics.

36 MATERIALS SCIENCE↗