Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupled mode theory”

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

Interaction of radio frequency waves with cylindrical density filaments: scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas , vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci. , vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

Physics↗
DOE OSTI2022

Interaction of radio frequency waves with cylindrical density filaments -- scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas, vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci., vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Plastic shrinkage of mortars cured with a paraffin-based compound – Bimodal neutron/X-ray tomography study

Early-age drying (immediately after casting) of mortars and the corresponding plastic shrinkage were studied using bimodal neutron/X-ray computed tomography. This novel, correlative 3D imaging mode enabled studying simultaneously and without any source of spurious perturbation the water migration and loss processes together with the corresponding deformations due to plastic shrinkage. Bimodal imaging opens up new possibilities for studying dynamic processes of coupled water transport and deformations in porous solids. The measurements were carried out on model systems (cylindrical mortar specimens with height of 19 mm). The study focused on the effect of a paraffin-based curing compound. Our results confirm that when the curing compound was applied directly onto the drying surface in a sufficient amount, both the evaporation rate and the rate of vertical displacement (settlement) were substantially reduced. The results shed a new light on the mechanisms of plastic shrinkage and the action of curing compounds.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Excitable dynamics in high-Lewis number premixed gas combustion at normal and microgravity

Freely-propagating, premixed gas flames in high-Lewis (Le) number, quiescent mixtures are studied experimentally in tubes of various diameter at normal (lg) and microgravity (mu g). A premixture of lean butane and oxygen diluted with helium, argon, neon, nitrogen or a mixture of multiple diluents is examined such that the thermal diffusivity of the mixture (and to a lesser extent, the mass diffusivity of the rate-limiting component) is systematically varied. In effect, different diluents allow variation of the Le without changing the chemistry. The flames are recorded with high speed cinematography and their stability is visually assessed. Different modes of propagation were observed depending on the diameter of the tubes (different conductive heat loss), the composition of the mixture and the g-level. At 1g, four modes of propagation were observed in small and intermediate diameter tubes (large conductive heat loss): (1) steadily propagating flames, (2) radial and longitudinal pulsating flames, (3) 'wavering' flames, and (4) rotating spiral flames. As the diameter of the tube increases, the radial modes become more pronounced while the longitudinal modes systematically disappear. Also, multiple, simultaneous, spatially-separated 'pacemaker' sites are observed in intermediate and large diameter tubes. Each site starts as a small region of high luminosity and develops into a flamelet which assumes the form of one of the fore mentioned modes. These flamelets eventually interact, annihilate each other in their regions of intersection and merge at their newly created free-ends. For very large tubes, radially-propagating wave-trains (believed to be 'trigger waves') are observed. These are analogous to the radial pulsations observed in the smaller diameter tubes. At mu g, three modes of propagation have been observed: (1) steadily propagating flames, (2) radial and longitudinal pulsating flames, and (3) multi-armed, rotating flames. Since the pulsating mode exists at mu g and 1g, buoyant flicker is not the mechanism which drives the pulsations. Moreover, all of the instabilities at 1g and mu g have characteristic frequencies on the O(100Hz). This value is lower than the fundamental, longitudinal acoustic frequencies of the tubes which suggests that the instabilities are not acoustically driven. The patterns formed by this reaction bear remarkable similarities with the patterns formed in most excitable media when the behavior of the system is driven by couplings between chemical reaction and diffusion (e.g., Belousov-Zhabotinsky reaction, Patterns in slime molds, spiral waves in the retina of a bird's eye). While it is recognized that the chemical mechanism associated with this premixed gas reaction is exponentially sensitive to temperature and undoubtedly different from those which govern previously observed excitable media (most are isothermal, or weakly exothermic, liquid phase reactions), similar spatial and temporal patterns should not come as a complete surprise considering heat and mass diffusion are self similar. It is concluded that this premixed gas system is a definitive example of a diffusive-thermal, gas-phase oscillator based on these experimental results and their favorable comparison with theory.

Pearlman, Howard↗
DOE OSTI2024

Linear response of molecular polaritons

In this article, we show that the collective light–matter strong coupling regime where N molecular emitters couple to the photon mode of an optical cavity can be mapped to a quantum impurity model where the photon is the impurity that is coupled to a bath of anharmonic transitions. In the thermodynamic limit where N ≫ 1, we argue that the bath can be replaced with an effective harmonic bath, leading to a dramatic simplification of the problem into one of the coupled harmonic oscillators. We derive simple analytical expressions for linear optical spectra (transmission, reflection, and absorption) where the only molecular input required is the molecular linear susceptibility. This formalism is applied to a series of illustrative examples, showing the role of temperature, disorder, vibronic coupling, and optical saturation of the molecular ensemble, explaining that it is useful even when describing an important class of nonlinear optical experiments. For completeness, we provide Appendixes A–C that include a self-contained derivation of the relevant spectroscopic observables for arbitrary anharmonic systems (for both large and small N) within the rotating-wave approximation. While some of the presented results herein have already been reported in the literature, we provide a unified presentation of the results as well as new interpretations that connect powerful concepts in open quantum systems and linear response theory with molecular polaritonics.

Chemistry↗
DOE OSTI2023

Fast particles in drift wave turbulence

This study aims to incorporate the effects of fast particles into our present fluid model for tokamak transport. The parameter ε f = ω / ω f, where ω is the mode frequency and ω f is the typical frequency of the fast particles, which enters as a factor in front of the fast particle response. Thus, for trapped fast particles, where ω f = ω pres the precession frequency of the fast particles, this parameter is of order 10 – 2 for drift waves, and thus, the fast particle response can be neglected. However, ε f will be of order 1 for fast particle modes such as in the fishbone instability. An important turbulence property, affecting both these limits, is resonance broadening. Effects of resonance broadening have recently been considered for fast particle instabilities, often coupled directly to the linear growth rate, while we here consider the original Dupree formulation where the turbulence directly drives a nonlinear frequency shift. Resonance broadening has a general tendency to counteract dissipative wave particle resonances. This has been observed for fast particle instabilities. Here, there is a resonant external source for the fast particles, so the instability survives if this source is dominant over the resonance broadening. For drift waves, however, external sources are not resonant since ε f << 1. Furthermore, the resonance broadening is able to remove the dissipative wave particle resonance completely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Polariton-assisted incoherent to coherent excitation energy transfer between colloidal nanocrystal quantum dots

Here, we explore the dynamics of energy transfer between two nanocrystal quantum dots placed within an optical microcavity. By adjusting the coupling strength between the cavity photon mode and the quantum dots, we have the capacity to fine-tune the effective coupling between the donor and acceptor. Introducing a non-adiabatic parameter, γ, governed by the coupling to the cavity mode, we demonstrate the system’s capability to shift from the overdamped Förster regime (γ $\ll$ 1) to an underdamped coherent regime (γ $\gg$ 1). In the latter regime, characterized by swift energy transfer rates, the dynamics are influenced by decoherence time. To illustrate this, we study the exciton energy transfer dynamics between two closely positioned CdSe/CdS core/shell quantum dots with sizes and separations relevant to experimental conditions. Employing an atomistic approach, we calculate the excitonic level arrangement, exciton–phonon interactions, and transition dipole moments of the quantum dots within the microcavity. These parameters are then utilized to define a model Hamiltonian. Subsequently, we apply a generalized non-Markovian quantum Redfield equation to delineate the dynamics within the polaritonic framework.

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

Catching the Killer: Dynamic Disorder Design Rules for Small–Molecule Organic Semiconductors

Small-molecule organic semiconductors are promising materials for applications ranging from solar cells to medical sensors. Their nearly infinite design space means that, in theory, it is possible to tailor the material to the exact specifications of a particular application. In reality; however, design rules to improve mobility (μ) remain elusive because of the complex calculations required to understand its limiters. Transient localization theory posits that charge carriers are slowed down by the collective phonon motions, called dynamic disorder, which localize charge carriers temporarily. Traditionally, a mode-by-mode analysis is performed to try and identify “killer” modes. Herein, the flaws are demonstrated with this analysis and present a streamlined simulation workflow that simplifies simulation of dynamic disorder and enables engineering of new molecular structures based on phonon dynamics. This workflow is combined with a novel visualization technique that enables per-atom insights into how μ is limited. Furthermore, this workflow is applied and analyzed to a series of -acenes and a series of BTBT-based molecules. Then design rules are identified in each series that identify possible ways to improve the μ beyond current experimental limits using phonon engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Interface velocity dependent solute trapping and phase selection during rapid solidification of laser melted hypo-eutectic Al-11at.%Cu alloy

Microstructure formation of a hypo-eutectic Al-11Cu (atom percent) alloy during solidification after laser melting has been studied by transmission electron microscopy (TEM). The evolution of the solid-liquid-interface velocity, V SL , during the solidification process has been determined from direct observation by in-situ TEM. This enabled correlating V SL with development of four distinct microstructure zones. Crystal growth mode transitions from planar to cellular, cellular to dendritic, dendritic to cellular, and cellular to planar, have been observed for the accelerating solid-liquid-interface. The transition from coupled two-phase growth to the single-phase growth occurred for V SL = V a = (0.80±0.05) m/s, where V a is the velocity of absolute stability, at the onset of banded morphology grain formation. The in-situ and post-mortem TEM uniquely permitted determination of the non-equilibrium solidus for the rapidly solidifying Al-11Cu alloy. Experimental evidence for solute clustering and chemical ordering tendencies at length scales on the order of ≤ 5nm has been detected in the single-phase regions of the banded grains. The structural features of the single-phase bands have been interpreted as signatures of ‘frozen in’ configurations present in the liquid boundary layer adjacent to the growing crystal, which has a width of about 3nm. Finally, the nano-scale spatiotemporal resolution experimental TEM studies performed here provided quantitative metrics, e.g. the solidification interface velocity dependent solute concentration of the α-Al phase and the near-atomic scale structure in the single-phase bands, that are uniquely suitable for comparison with theory and model predictions for solidification microstructure development in multicomponent alloys after laser melting.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Single Entity Electrochemistry in Nanopore Electrode Arrays: Ion Transport Meets Electron Transfer in Confined Geometries

Electrochemical measurements conducted in confined volumes provide a powerful and direct means to address scientific questions at the nexus of nanoscience, biotechnology, and chemical analysis. How are electron transfer and ion transport coupled in confined volumes and how does understanding them require moving beyond macroscopic theories? Also, how do these coupled processes impact electrochemical detection and processing? We address these questions by studying a special type of confined-volume architecture, the nanopore electrode array, or NEA, which is designed to be commensurate in size with physical scaling lengths, such as the Debye length, a concordance that offers performance characteristics not available in larger scale structures. The experiments described here depend critically on carefully constructed nanoscale architectures that can usefully control molecular transport and electrochemical reactivity. We begin by considering the experimental constraints that guide the design and fabrication of zero-dimensional nanopore arrays with multiple embedded electrodes. These zero-dimensional structures are nearly ideal for exploring how permselectivity and unscreened ion migration can be combined to amplify signals and improve selectivity by enabling highly efficient redox cycling. Our studies also highlight the benefits of arrays, in that molecules escaping from a single nanopore are efficiently captured by neighboring pores and returned to the population of active redox species being measured, benefits that arise from coupling ion accumulation and migration. These tools for manipulating redox species are well-positioned to explore single molecule and single particle electron transfer events through spectroelectrochemistry, studies which are enabled by the electrochemical zero-mode waveguide (ZMW), a special hybrid nanophotonic/nanoelectronic architecture in which the lower ring electrode of an NEA nanopore functions both as a working electrode to initiate electron transfer reactions and as the optical cladding layer of a ZMW. Furthermore, while the work described here is largely exploratory and fundamental, we believe that the development of NEAs will enable important applications that emerge directly from the unique coupled transport and electron-transfer capabilities of NEAs, including in situ molecular separation and detection with external stimuli, redox-based electrochemical rectification in individually encapsulated nanopores, and coupled sorters and analyzers for nanoparticles.

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

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗
NASA NTRSDate not supplied

Intermittency Model for Coupled CFD-Stability Transition Analysis in Hypersonic Flow

Accurate prediction of aerothermodynamic loads on hypersonic vehicles requires precise modeling of surface quantities across the boundary layer transition zone. The peak heating loads and total heat transfer are determined by parameters such as the transition zone length and the magnitude of potential overshoots in heat flux and skin friction beyond their respective values in fully turbulent flows. While previous studies on CFD integrated transition modeling using linear stability correlations have shown promise in modeling these features, they did not develop an intermittency model for hypersonic flows. This paper presents a data-driven approach to develop a model of this type by correlating experimental transition data with the relevant flow parameters. The resulting model demonstrates significant improvements over previous low-speed models in terms of predicting heat transfer distributions during the transition process associated with first and second mode instabilities in axisymmetric high-speed flows. Separate correlations for flight and ground test conditions are developed, and the potential to combine these correlations is discussed.

Transition↗
NASA NTRS2023

CFD Integrated Transition Modeling for High-Speed Flows via Coupled OVERFLOW-LASTRAC Analysis

This work details ongoing efforts at the NASA Langley Research Center to develop and validate a general-use CFD tool that includes built-in predictions of boundary-layer transition in high-speed flows. Existing tools for the coupling of the NASA OVERFLOW structured overset RANS solver and LASTRAC stability analysis code have been extended to capture boundary-layer transition in high-speed flows driven by either Mack’s first-mode or secondmode instabilities and/or crossflow instabilities. The efficacy of this coupled approach has been demonstrated by examining a variety of supersonic and hypersonic test cases for which experimental validation data is available. The intermittency prescription parameters inherited from prior low-speed applications are found to be suitable for some, but not all high-speed flow scenarios. The method is found to be robust for straight cone configurations and improvements necessary to accurately capture transition on geometries with spatially complex amplification factor envelopes with intermediate regions of slow N-factor variation are examined. The automated, coupled analysis is demonstrated for a 3d supersonic test case and found to perform well within the limits of the linear stability analysis on which it relies. Finally, a preliminary investigation of the method’s robustness to non-ideal CFD meshes is conducted.

High-speed flow↗
DOE OSTI2021

Charged lepton flavor violation at the EIC

We present a comprehensive analysis of the potential sensitivity of the Electron-Ion Collider (EIC) to charged lepton flavor violation (CLFV) in the channel ep → τX, within the model-independent framework of the Standard Model Effective Field Theory (SMEFT). We compute the relevant cross sections to leading order in QCD and electroweak corrections and perform simulations of signal and SM background events in various τ decay channels, suggesting simple cuts to enhance the associated estimated efficiencies. To assess the discovery potential of the EIC in τ-e transitions, we study the sensitivity of other probes of this physics across a broad range of energy scales, from pp → eτX at the Large Hadron Collider to decays of B mesons and τ leptons, such as τ → eγ, τ → eℓ+ℓ–, and crucially the hadronic modes τ → eY with Yϵ {π, K, ππ, Kπ, …}. We find that electroweak dipole and four-fermion semi-leptonic operators involving light quarks are already strongly constrained by τ decays, while operators involving the c and b quarks present more promising discovery potential for the EIC. An analysis of three models of leptoquarks confirms the expectations based on the SMEFT results. We also identify future directions needed to maximize the reach of the EIC in CLFV searches: these include an optimization of the τ tagger in hadronic channels, an exploration of background suppression through tagging b and c jets in the final state, and a global fit by turning on all SMEFT couplings, which will likely reveal new discovery windows for the EIC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Stability of Kuramoto networks subject to large and small fluctuations from heterogeneous and spatially correlated noise

Oscillatory networks subjected to noise are broadly used to model physical and technological systems. Due to their nonlinear coupling, such networks typically have multiple stable and unstable states that a network might visit due to noise. In this article, we focus on the assessment of fluctuations resulting from heterogeneous and spatially correlated noise inputs on Kuramoto model networks. We evaluate the typical, small fluctuations near synchronized states and connect the network variance to the overlap between stable modes of synchronization and the input noise covariance. Going beyond small to large fluctuations, we introduce the indicator mode approximation that projects the dynamics onto a single amplitude dimension. Such an approximation allows for estimating rates of fluctuations to saddle instabilities, resulting in phase slips between connected oscillators. Statistics for both regimes are quantified in terms of effective noise amplitudes that are compared and contrasted for several noise models. In conclusion, bridging the gap between small and large fluctuations, we show that a larger network variance does not necessarily lead to higher rates of large fluctuations.

42 ENGINEERING↗
NASA NTRS1993

Ignition and subsequent flame spread over a thin cellulosic material

Both ignition and flame spread on solid fuels are processes that not only are of considerable scientific interest but that also have important fire safety applications. Both types of processes, ignition and flame spread, are complicated by strong coupling between chemical reactions and transport processes, not only in the gas phase but also in the condensed phase. In most previous studies, ignition and flame spread were studied separately with the result that there has been little understanding of the transition from ignition to flame spread. In fire safety applications this transition is crucial to determine whether a fire will be limited to a localized, temporary burn or will transition into a growth mode with a potential to become a large fire. In order to understand this transition, the transient mechanisms of ignition and subsequent flame spread must be studied. However, there have been no definitive experimental or modeling studies, because of the complexity of the flow motion generated by buoyancy near the heated sample surface. One must solve the full Navier-Stokes equations over an extended region to represent accurately the highly unstable buoyant plume and entrainment of surrounding gas from far away. In order to avoid the complicated nature of the starting plume problem under normal gravity, previous detailed radiative ignition models were assumed to be one-dimensional or were applied at a stagnation point. Thus, these models cannot be extended to include the transition to flame spread. The mismatch between experimental and calculated geometries means that theories cannot be compared directly with experimental results in normal gravity. To overcome the above difficulty, theoretical results obtained without buoyancy can be directly compared with experimental data measured in a microgravity environment. Thus, the objective of this study is to develop a theoretical model for ignition and the transition to flame spread and to make predictions using the thermal and chemical characteristics of a cellulosic material which are measured in normal gravity.

Nakabe, Kazuyoshi↗
DOE OSTI2020

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving the linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the resistive DCON code [A. H. Glasser, Z. R. Wang, and J.-K. Park, Phys. Plasmas 23, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Additionally, the solution of the ideal outer region is critical to the asymptotic matching and is challenging in toroidal geometry due to the singular power series solution at the resonant surfaces. Thus, a systematic verification of the outer region Δ' matrix is made by reproducing the well-known analytical Δ' result in Furth et al. [Phys. Fluids 16, 1054–1063 (1073)] and by making a quantitative benchmark with the PEST3 code [A. Pletzer and R. L. Dewar, J. Plasma Phys. 45, 427–451 (1991)]. Finally, the reconstructed numerical solution of the resistive plasma response from the toroidal matching matrix is presented. Compared with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

From Atom to Engine: Understanding Fundamental Effects of Structure on Combustion Using Tandem Experiment and Computation

The ability to tie structural features to a fuel candidate’s ignition properties provides a path for rational design of advantaged fuels that facilitate higher-efficiency, lower-emitting combustion in engines. The time evolution of a fuel’s radical population directly correlates with global combustion parameters, and because the generated radicals depend on the initial structure and reaction pathways of the fuel molecule, it follows that the molecular structure of a fuel has a direct impact on this cascade, and therefore salient ignition-tied fuel properties. This work aims to understand the impact of structural features in the context of global combustion properties for direct application to emerging engine technologies. Specifically, we highlight the ability to identify potential fuel blendstocks for specific engine strategies by examining and elucidating the radical cascade tied to ignition. A representative subset of alcohols possessing a variety of structural features (e.g. branching, chain length) was selected for study using a robust combination of experimental and computational methods. The approach begins with identifying short-lived radicals and likely abstraction sites upon pyrolysis in a microreactor coupled with photoionization mass spectrometry for direct detection of product species. In tandem, electronic structure calculations were performed to calculate the relevant potential energy surfaces for fuel pyrolysis as well as the subsequent O2 addition-isomerization reactions, which govern ignition behavior. Furthermore, the theoretical energetic and ro-vibrational data are used to calculate rate constants for direct ignition modeling to explore how structural effects impact NTC behavior. The microreactor data combined with electronic structure theory calculations provide a well-informed picture of the reactions relevant to ignition for each fuel in this study, as well as the most likely radicals generated under engine conditions. This thorough fundamental picture of fuel decomposition and ignition chemistry is used to correlate the observed ignition behavior to molecular structural effects.

advanced biofuels↗