Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Timing Analysis”

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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Forward sensitivity analysis and mode dependent control for closure modeling of Galerkin systems

Model reduction by projection-based approaches is often associated with losing some of the important features that contribute towards the dynamics of the retained scales. As a result, a mismatch occurs between the predicted trajectories of the original system and the truncated one. We put forth a framework to apply a continuous time control signal in the latent space of the reduced order model (ROM) to account for the effect of truncation. We set the control input using parameterized models by following energy transfer principles. Our methodology relies on observing the system behavior in the physical space and using the projection operator to restrict the feedback signal into the latent space. Then, we leverage the forward sensitivity method (FSM) to derive relationships between the feedback and the desired mode-dependent control. We test the performance of the proposed approach using two test cases, corresponding to viscous Burgers and vortex merger problems at high Reynolds number. Results show that the ROM trajectory with the applied FSM control closely matches its target values in both the data-dense and data-sparse regimes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A new framework for X-ray photon correlation spectroscopy analysis from polycrystalline materials

Here, we report a new analytical framework for interpreting data from X-ray photon correlation spectroscopy (XPCS) measurements on polycrystalline materials characterized by strong scattering intensity variations at fixed wavevector magnitude (i.e., anisotropic scattering). Currently, no analytical method exists for the interpretation of the time-dependent anisotropic scattering from such materials. The framework is applied to interrogate the dynamics of a spherical micelle-forming diblock copolymer melt below the order-disorder transition, wherein finite size grains of micellar body-centered cubic structure produce anisotropic scattering. A wealth of analytical information is recovered from a simple measurement, including distributions of relaxation times and speeds associated with micelles within grains. The findings of this study demonstrate the efficacy of this new analytical method, which may be readily adapted for application to a variety of materials and systems.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Stochastic Dynamic Analysis for Power Systems Under Uncertain Variability

This paper proposes a novel method to analyze impacts of uncertain variability on power system dynamics. There is considerable interest in integrating intermittent renewable energy and plug-in electric vehicles into power systems. Therefore, power systems are in an environment with increasing uncertain variability. In this paper, the uncertain variability is described as a continuous-time stochastic process, and thus, the power system under uncertain variability is modeled by stochastic differential equations. To quantify impacts of uncertain variability on power system dynamics, an intra-region probability index is presented. Based on the stochastic averaging method, an analytical method is proposed to calculate the intra-region probability. Compared to Monte Carlo simulation, the proposed method is many orders of magnitude faster without sacrificing the result accuracy. As a result, several insights are given to improve the power system dynamics under uncertain variability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Long-Term Dynamics of Atmospheric Sulfur Dioxide in Urban and Rural Regions of China: Urbanization and Policy Impacts

High levels of sulfur dioxide (SO 2 ) due to human activities pose a serious air pollution issue in China, especially in urban agglomerations. However, limited research has investigated the impact of anthropogenic emissions on higher SO 2 concentrations in urban regions compared to rural areas in China. Here, we analyzed the trends in SO 2 concentrations from 1980 to 2021 in China using the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) dataset. SO 2 column concentrations from the Copernicus Atmosphere Monitoring Service (CAMS) and the Ozone Monitoring Instrument (OMI) during the years 2007–2021 were also examined for validation and comparison purposes. Eight representative areas, including four urban regions (Pearl River Delta [PRD], Beijing-Tianjin-Hebei [BTH], Yangtze River Delta [YRD], and Sichuan Basin [SCB]) and four rural regions (Northeast Region [NER], Mongolian Region [MR], West Region [WR], and Tibetan Plateau Region [TR]) were selected for the analysis. Overall, a significant but fluctuating increase in SO 2 concentrations over China was observed during 1980–2021. During 1980–1997 and 2000–2010, there was an increase in SO 2 concentration, while during 1997–2000 and 2010–2021, a decreasing trend was observed. The average increase in SO 2 concentration was approximately 16 times higher in urban regions than in the rural background. We also found that SO 2 dynamics were highly associated with expansion of urban areas, population density, and gross domestic product. Nonetheless, since 2007, SO 2 concentrations have exhibited a downward trend, which is mainly attributed to the air pollution policies implemented by the Chinese government. Our findings highlight the need for further studies on the impact of SO 2 on regional climate change in China.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

PySFD: comprehensive molecular insights from significant feature differences detected among many simulated ensembles

Motivation Many modeling analyses of molecular dynamics (MD) simulations are based on a definition of states that can be (groups of) clusters of simulation frames in a feature space composed of molecular coordinates. With increasing dimension of this feature space (due to the increasing size or complexity of a simulated molecule), it becomes very difficult to cluster the underlying MD data and estimate a statistically robust model. To mitigate this “curse of dimensionality”, one can reduce the feature space, e.g., with principal component or time-lagged independent component analysis transformations, focusing the analysis on the most important modes of transitions. In practice, however, all these reduction strategies may neglect important molecular details that are susceptible to experimental verification. Results To recover such molecular details, I have developed PySFD (Significant Feature Differences analyzer for Python), a multi-processing software package that efficiently selects significantly different features of any user-defined feature type among potentially many different simulated state ensembles, such as meta-stable states of a Markov State Model (MSM). Applying PySFD on MSMs of an aggregate of 300 microseconds MD simulations recently performed on the major histocompatibility complex class II (MHCII) protein, I demonstrate how this toolkit can extract and visualize valuable mechanistic information from big MD simulation data, e.g., in form of networks of dynamic interaction changes connecting functionally relevant sites of a protein complex.

Biochemistry & Molecular Biology↗
DOE OSTI2026

Multiprobe constraints on early and late time dark energy

We perform a multiprobe analysis combining cosmic microwave background (CMB) data from Planck and the Atacama Cosmology Telescope (ACT), ACT CMB lensing, and large-scale structure (LSS) measurements from the Dark Energy Spectroscopic Instrument (DESI), including DESI Legacy Imaging Survey (LS) galaxies and baryon acoustic oscillations (BAOs). We present the first 5×2 pt analysis of ACT DR6 lensing, DESI LS, and Planck Integrated Sachs-Wolfe (ISW). Within Λ⁢CDM, this yields 𝑆 8 =𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 =0.819 ±0.016, in good agreement with primary CMB inferences and provides a sound-horizon-free Hubble constant constraint of 𝐻 0 =70.0±4.4 km s −1 Mpc −1 . Then, combining with CMB primary and BAO, we reconfirm a CMB–BAO discrepancy in the Ω 𝑚 –$\frac{𝐷𝑣}{𝑟𝑑}$ plane, which is heightened when combining BAO with the 5 ×2 pt data vector. We explore two dark-energy extensions that may reconcile this: an early-time modification, early dark energy (EDE), and late-time dynamical dark energy (DDE) parametrized by 𝑤 0 ⁢𝑤 𝑎 . For CMB primary +BAO +5 ×2 pt, we find a 3.3⁢𝜎 preference for DDE over Λ⁢CDM, while EDE is modestly favored at 2.3⁢𝜎. The models address different shortcomings of Λ⁢CDM: DDE relaxes the neutrino mass bound (𝑀 𝜈 <0.17 eV vs <0.050 eV under Λ⁢CDM), making it compatible with neutrino oscillation measurements, while EDE raises the Hubble constant to 𝐻 0 =70.5±1.2 km s −1 Mpc −1 , easing the discrepancy with SH0ES. However, neither model resolves both issues simultaneously. Our analysis indicates that both DDE and EDE remain viable extensions of Λ⁢CDM within current uncertainties and demonstrates the capacity of combined probes to place increasingly stringent constraints on cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Hybrid methods for cybersecurity analysis :

Early 2010 saw a signi cant change in adversarial techniques aimed at network intrusion: a shift from malware delivered via email attachments toward the use of hidden, embedded hyperlinks to initiate sequences of downloads and interactions with web sites and network servers containing malicious software. Enterprise security groups were well poised and experienced in defending the former attacks, but the new types of attacks were larger in number, more challenging to detect, dynamic in nature, and required the development of new technologies and analytic capabilities. The Hybrid LDRD project was aimed at delivering new capabilities in large-scale data modeling and analysis to enterprise security operators and analysts and understanding the challenges of detection and prevention of emerging cybersecurity threats. Leveraging previous LDRD research e orts and capabilities in large-scale relational data analysis, large-scale discrete data analysis and visualization, and streaming data analysis, new modeling and analysis capabilities were quickly brought to bear on the problems in email phishing and spear phishing attacks in the Sandia enterprise security operational groups at the onset of the Hybrid project. As part of this project, a software development and deployment framework was created within the security analyst work ow tool sets to facilitate the delivery and testing of new capabilities as they became available, and machine learning algorithms were developed to address the challenge of dynamic threats. Furthermore, researchers from the Hybrid project were embedded in the security analyst groups for almost a full year, engaged in daily operational activities and routines, creating an atmosphere of trust and collaboration between the researchers and security personnel. The Hybrid project has altered the way that research ideas can be incorporated into the production environments of Sandias enterprise security groups, reducing time to deployment from months and years to hours and days for the application of new modeling and analysis capabilities to emerging threats. The development and deployment framework has been generalized into the Hybrid Framework and incor- porated into several LDRD, WFO, and DOE/CSL projects and proposals. And most importantly, the Hybrid project has provided Sandia security analysts with new, scalable, extensible analytic capabilities that have resulted in alerts not detectable using their previous work ow tool sets.

Dunlavy, Daniel M.↗
DOE OSTI2012

Numerical Analysis of the Dynamics of Two- and Three-Dimensional Fluidized Bed Reactors using an Euler-Lagrange Approach

Biomass thermochemical conversion, often done in fluidized beds, recently gained a lot of attention due to its potential to efficiently produce renewable liquid fuels. Optimization of reactor design and operating conditions, however, requires a fundamental understanding of bed dynamics. In this work, a numerical framework based on an Euler-Lagrange approach is developed and used to perform and analyze large-scale simulations of two- and three-dimensional periodic fluidized beds. Collisions are handled using a soft-sphere model. An efficient parallel implementation allows one to explicitly track over 30 million particles, which is representative of the number of particles found in lab-scale reactor, therefore demonstrating the capability of Lagrangian approaches to simulate realistic systems at that scale. An on-the-fly bubble identification and tracking algorithm is used to characterize bubble dynamics for inlet velocities up to 9 times the minimum fluidization velocity. Statistics for gas volume fraction, gas and particle velocities, bed expansion, and bubble size and velocity, is compared across the two- and three-dimensional configurations, and comparison with literature data generally shows good agreement. The wide distribution of gas residence times observed in the simulations is linked to the different gas hold-up characteristics of the gas-solid system.

09 BIOMASS FUELS↗
DOE OSTI2019

Interaction of the Westerlies with the Tibetan Plateau in Determining the Mei-Yu Termination

This article explores how the termination of the mei-yu is dynamically linked to the westerlies impinging on the Tibetan Plateau. It is found that the mei-yu stage terminates when the maximum upper-tropospheric westerlies shift beyond the northern edge of the plateau, around 40°N. This termination is accompanied by the disappearance of tropospheric northerlies over northeastern China. The link between the transit of the jet axis across the northern edge of the plateau, the disappearance of northerlies, and termination of the mei-yu holds on a range of time scales from interannual through seasonal and pentad. Diagnostic analysis indicates that the weakening of the meridional moisture contrast and meridional wind convergence, mainly resulting from the disappearance of northerlies, causes the demise of the mei-yu front. The authors propose that the westerlies migrating north of the plateau and consequent weakening of the extratropical northerlies triggers the mei-yu termination. Model simulations are employed to test the causality between the jet and the orographic downstream northerlies by repositioning the northern edge of the plateau. As the plateau edge extends northward, orographic forcing on the westerlies strengthens, leading to persistent strong downstream northerlies and a prolonged mei-yu. Idealized simulations with a dry dynamical core further demonstrate the dynamical link between the weakening of orographically forced downstream northerlies with the positioning of the jet from south to north of the plateau. Changes in the magnitude of orographically forced stationary waves are proposed to explain why the downstream northerlies disappear when the jet axis migrates beyond the northern edge of the plateau.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Probing Surface Plasmon Dynamics in Periodic Nanostructures through Ultrafast Electron Microscopy

Surface plasmon polaritons (SPPs) can be manipulated to localize and guide light in subwavelength distances, enabling them to find applications in a wide range of areas, from sensing to quantum computing. Among several methods of SPP excitation, periodic arrays of nano- and microstructures are of particular interest, as they enable engineering SPP properties through structural parameters. Here, in this study, using the photon-induced near-field electron microscopy (PINEM) technique, we investigated the mode formation, coupling, interference, and decay of SPPs in square and hexagonal arrays of circular nanoholes under both visible and near-infrared excitation. Polarization-resolved analysis revealed the key factors governing SPP localization and interference patterns, showing that the periodicity and symmetry of the array primarily determine the SPP interference patterns and their orientation, while pump polarization mainly modulates their intensity. Time-resolved PINEM measurements demonstrated the spatial dependence of the SPP temporal characteristics. In addition, cathodoluminescence (CL) spectroscopy was employed to examine the intrinsic plasmonic characteristics of the structure. Finite difference time domain (FDTD) simulations showed strong agreement with both PINEM and CL measurements on the spatial and spectral behavior of SPPs. Understanding the spatiotemporal dynamics of SPPs on nanostructures beyond the diffraction limit is crucial for optimizing plasmonic structures for advanced photonic and quantum technologies.

Plasmonics↗
DOE OSTI2013

New Methods of Uncertainty Quantification for Mixed Discrete-Continuous Variable Models

The scale and complexity of problems such as designing power grids or planning for climate change is growing rapidly, driving the development of complicated computer models. More complex models have longer run times and incorporate larger numbers of inputs, both continuous and discrete. For example, a detailed physics model may have continuous variables such as temperature, height or pressure along with discrete variables that indicate the choice of a material for a particular piece or the model to be used to calculate air flow. A power grid design model may have continuous variables such as generation capacity, power flow or demand along with discrete variables such as number of generators, number of transmission lines or binary variables to indicate whether or not a node is chosen for generation expansion. A growing awareness of uncertainty and the desire to make risk-informed decisions is causing uncertainty quantification (UQ) to be more routine and often required. UQ provides the underpinnings necessary to establish confidence in models and their use; therefore, much time and effort is being invested in creating efficient approaches for UQ. However, these efforts have been focused on models that take continuous variables as inputs. When discrete inputs are thrown into the mix, the basic approach is to repeat the UQ analysis for each combination of discrete inputs or some subset thereof; this rapidly becomes intractable. Because of the computational complexity inherent in mixed discrete-continuous models, researchers will focus on the uncertainty in their particular problem finding ways to take advantage of symmetries, simplifications or structures. For example, uncertainty propagation in certain dynamical systems can be efficiently carried out after various decomposition steps or uncertainty propagation in stochastic programming is confined to scenario generation. Unfortunately models are not always available for such machinations: models may be embedded in legacy codes, may utilize commercial off the shelf codes or may be created by stringing a series of codes together. It is also time consuming to start each problem from scratch; worse there may not be any simplifications or symmetries to take advantage of. For these situations a UQ method developed for any black box function is necessary. This report documents a new conceptual model for performing UQ for mixed discrete-continuous models which not only applies to any simulator function, but allows the use of the efficient UQ methods that have been developed for continuous inputs only. The conceptual model is presented and an estimation procedure is fleshed out for one class of problems. This is applied to variations of a mixed discrete-continuous optimization test problem. This procedure provides comparable results to a benchmark solution with fewer function evaluations.

Bauman, Lara E.↗
DOE OSTI2020

Evolution of AISI 304L stainless steel part properties due to powder recycling in laser powder-bed fusion

Selective laser melting (SLM) is a powder-bed fusion process in Additive Manufacturing (AM) that uses a laser beam to selectively fuse layers of powder into near net-shape components with little porosity. However, inconsistencies in the part properties due to the presence of defects in as-built components has hindered the widespread adoption of SLM for industrial applications forcing researchers to study the sources of variation for quality control purposes. A critical area suspected of creating variation in the part properties is the feedstock, where batch-to-batch differences as well as changes in the powder properties with reuse have the potential to affect performance. During processing, laser spatter and condensate eject from the melt pool and deposit into the powder-bed surrounding the parts. These particulates, more commonly known as ejecta, differ morphologically and chemically from the virgin powder compromising its reusability. In this study, 304L stainless steel powder was recycled for a total of 5 times through a systematic approach aimed at accelerating powder degradation to reveal the influence on both the tensile and impact toughness properties. Through analysis of variance (ANOVA), it was found that tensile properties did not change with reuse, while the impact toughness showed a steady decline illustrating the differences in static and dynamic part properties due to powder recycling.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Lyapunov spectrum of the separated flow around the NACA 0012 airfoil and its dependence on numerical discretization

We investigate the impact of numerical discretization on the Lyapunov spectrum of separated flow simulations. The two-dimensional chaotic flow around the NACA 0012 airfoil at a low Reynolds number and large angle of attack is considered to that end. Time, space and accuracy-order refinement studies are performed to examine each of these effects separately. Numerical results show that the time discretization has a small impact on the dynamics of the system, whereas the spatial discretization can dramatically change them. Also, the finite-time Lyapunov exponents associated to unstable modes are shown to be positively skewed, and quasi-homoclinic tangencies are observed in the attractor of the system. The implications of these results on flow physics and sensitivity analysis of chaotic flows are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Origin of Stabilization of Ligand-Centered Mixed Valence Ruthenium Azopyridine Complexes: DFT Insights for Neuromorphic Applications

Redox-driven conductance changes are critical processes in molecular- and coordination-complex-based memristive thin films and devices that are envisioned for neuromorphic technologies, but fundamental mechanisms of conductance switching are not fully understood. Here, we explore charge disproportionation (CD) processes in [Ru II L 2 ](PF 6 ) 2 molecular systems that intrinsically involve interfragment charge transfer (IFCT). Using a combination of ab initio molecular dynamics simulation (AIMD), time-dependent density functional theory (TD-DFT), and density functional theory (DFT) calculations, we investigate the electron transfer mechanisms and the roles of temperature and cell volumetric expansion in facilitating the counterion movements and electronic transitions required for low-cost IFCT and charge redistribution. A detailed analysis of the density of states and TD-DFT calculations highlights that unpaired electrons play a crucial role in low-energy transitions, with the azo (N=N) groups of the ligand serving as the primary sites for electronic transport between molecular fragments, further stabilizing the asymmetric state. Localization of added electrons on azo ligands occurs with negligible change at the Ru centers, supported by atomic volume expansions up to +4.74 bohr 3 , and goes along with a progressive reduction of the HOMO−LUMO gap across redox states, suggesting enhanced conductivity. The TD-DFT analysis reveals a dominant IFCT excitation at 2082.76 nm in the doubly reduced (22) state, while a stabilization energy of 1.20 eV of the asymmetric (13) state relative to the symmetric (22) state is predicted by constrained DFT. Periodic DFT and AIMD simulations emulating a molecular film show that the stabilization of the asymmetric state, relative to a symmetric one, translates in net charge separation values (order of ∼0.33 e) that are strongly linked to increased counterion mobility (average counterion displacements exceeding 0.7 Å per atom during CD events) and the involvement of azo groups in electron redistribution. These findings, which align with previously reported experimental and computational data, provide key insights into the IFCT mechanisms and electronic transport facilitated by azo groups, with important implications for redox-driven memristive and neuromorphic technologies.

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

Improving Time Step Convergence in an Atmosphere Model With Simplified Physics: The Impacts of Closure Assumption and Process Coupling

Convergence testing is a common practice in the development of dynamical cores of atmospheric models but is not as often exercised for the parameterization of sub-grid physics. An earlier study revealed that the stratiform cloud parameterizations in several predecessors of the Energy Exascale Earth System Model (E3SM) showed strong time-step sensitivity and slower-than-expected convergence when the model's time step was systematically refined. In this work, a simplified atmosphere model is configured that consists of the spectral-element dynamical core of the E3SM atmosphere model coupled with a large-scale condensation parameterization based on commonly used assumptions. This simplified model also resembles E3SM and its predecessors in the numerical implementation of process coupling and shows poor time-step convergence in short ensemble tests. We present a formal error analysis to reveal the expected time-step convergence rate and the conditions for obtaining such convergence. Numerical experiments are conducted to investigate the root causes of convergence problems. We show that revisions in the process coupling and closure assumption help to improve convergence in short simulations using the simplified model; the same revisions applied to a full atmosphere model lead to significant changes in the simulated long-term climate. This work demonstrates that causes of convergence issues in atmospheric simulations can be understood by combining analyses from physical and mathematical perspectives. Addressing convergence issues can help to obtain a discrete model that is more consistent with the intended representation of the physical phenomena.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Stochastic journeys of cell progenies through compartments and the role of self-renewal, symmetric and asymmetric division

Abstract Division and differentiation events by which cell populations with specific functions are generated often take place as part of a developmental programme, which can be represented by a sequence of compartments. A compartment is the set of cells with common characteristics; sharing, for instance, a spatial location or a phenotype. Differentiation events are transitions from one compartment to the next. Cells may also die or divide. We consider three different types of division events: (i) where both daughter cells inherit the mother’s phenotype (self-renewal), (ii) where only one of the daughters changes phenotype (asymmetric division), and (iii) where both daughters change phenotype (symmetric division). The self-renewal probability in each compartment determines whether the progeny of a single cell, moving through the sequence of compartments, is finite or grows without bound. We analyse the progeny stochastic dynamics with probability generating functions. In the case of self-renewal, by following one of the daughters after any division event, we may construct lifelines containing only one cell at any time. We analyse the number of divisions along such lines, and the compartment where lines terminate with a death event. Analysis and numerical simulations are applied to a five-compartment model of the gradual differentiation of hematopoietic stem cells and to a model of thymocyte development: from pre-double positive to single positive (SP) cells with a bifurcation to either SP4 or SP8 in the last compartment of the sequence.

Science & Technology - Other Topics↗
DOE OSTI2023

Automated instant labeling chemistry workflow for real-time monitoring of monoclonal antibody N -glycosylation

With the transition toward continuous bioprocessing, process analytical technology (PAT) is becoming necessary for rapid and reliable in-process monitoring during biotherapeutics manufacturing. Bioprocess 4.0 is looking to build end-to-end bioprocesses that include PAT-enabled real-time process control. This is especially important for drug product quality attributes that can change during bioprocessing, such as protein N-glycosylation, a critical quality attribute for most monoclonal antibody (mAb) therapeutics. Glycosylation of mAbs is known to influence their efficacy as therapeutics and is regulated for a majority of mAb products on the market today. Currently, there is no method to truly measure N-glycosylation using on-line PAT, hence making it impractical to design upstream process control strategies. We recently described the N-GLYcanyzer workflow: an integrated PAT unit that measures mAb N-glycosylation within 3 hours of automated sampling from a bioreactor. Here, we integrated Agilent's Instant Procainamide (InstantPC) based chemistry workflow into the N-GLYcanyzer PAT unit to allow for nearly 10× faster near real-time analysis of mAb glycoforms. Furthermore, our methodology is explained in detail to allow for replication of the PAT workflow as well as present a case study demonstrating the use of this PAT to autonomously monitor a mammalian cell perfusion process at the bench scale to gain increased knowledge of mAb glycosylation dynamics during continuous biologics manufacturing using Chinese hamster ovary (CHO) cells.

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

Self-consistent nonlinear kinetic simulations of the anomalous Doppler instability of suprathermal electrons in plasmas

Suprathermal tails in the distributions of electron velocities parallel to the magnetic field are found in many areas of plasma physics, from magnetic confinement fusion to solar system plasmas. Parallel electron kinetic energy can be transferred into plasma waves and perpendicular gyration energy of particles through the anomalous Doppler instability (ADI), provided that energetic electrons with parallel velocities v{sub ||}≥(ω+Ω{sub ce})/k{sub ||} are present; here Ω{sub ce} denotes electron cyclotron frequency, ω the wave angular frequency, and k{sub ||} the component of wavenumber parallel to the magnetic field. This phenomenon is widely observed in tokamak plasmas. Here, we present the first fully self-consistent relativistic particle-in-cell simulations of the ADI, spanning the linear and nonlinear regimes of the ADI. We test the robustness of the analytical theory in the linear regime and follow the ADI through to the steady state. By directly evaluating the parallel and perpendicular dynamical contributions to j·E in the simulations, we follow the energy transfer between the excited waves and the bulk and tail electron populations for the first time. We find that the ratio Ω{sub ce}/(ω{sub pe}+Ω{sub ce}) of energy transfer between parallel and perpendicular, obtained from linear analysis, does not apply when damping is fully included, when we find it to be ω{sub pe}/(ω{sub pe}+Ω{sub ce}); here ω{sub pe} denotes the electron plasma frequency. We also find that the ADI can arise beyond the previously expected range of plasma parameters, in particular when Ω{sub ce}>ω{sub pe}. The simulations also exhibit a spectral feature which may correspond to the observations of suprathermal narrowband emission at ω{sub pe} detected from low density tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗