Search NASA⌕ Search

SEARCH · Search NASA

Results for “heat wave”

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

Scattering of radio frequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radio frequency (RF) electromagnetic waves that propagate through this region. It is important to know, whether used for diagnostics or for heating and current drive, the spectral properties of these scattered RF waves. The spectral changes influences the interpretation of the diagnostic-data obtained and the current and heating profiles. A full-wave, 3D electromagnetic code ScaRF (see Papadopoulos et al. 2019) has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma, anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can be for the study of scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U, DIII-D. For the case of medium-sized tokamaks, there's experimental evidence that drift waves and rippling modes are present in the edge region (see Ritz et al. 1984). Hence, we study the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes (see Papadopoulos et al. 2019). The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion (see Xiu & Karniadakis 2002) method in conjunction with the Smolyak sparse grid integration (see Papadopoulos et al. 2018) (PCE-SG). The PCE-SG method is proven accurate and much more efficient (roughly 2-orders of magnitude shorter execution time) compared to alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Higher-order particle representation for particle-in-cell simulations

In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r 0 . A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r 0 . As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. In this work, good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Predictions of core plasma performance for the SPARC tokamak

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with Ion Cyclotron Range-of-Frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach Q~11, well above its mission objective of Q>2. To build confidence that SPARC will be successful, physics-based integrated modeling has also been performed. The TRANSP code coupled with the theory-based TGLF turbulence model and EPED predictions for pedestal stability find that Q~9 is attainable in standard H-mode operation and confirms Q>2 operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modeled with the Monte-Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Predictions of core plasma performance for the SPARC tokamak

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with ion cyclotron range-of-frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach $\textit{Q}$ ≈ 11, which is well above its mission objective of $\textit{Q}$ > 2. To build confidence that SPARC will be successful, physics-based integrated modelling has also been performed. The TRANSP code coupled with the theory-based trapped gyro-Landau fluid (TGLF) turbulence model and EPED predictions for pedestal stability find that $\textit{Q}$ ≈ 9 is attainable in standard H-mode operation and confirms $\textit{Q}$ > 2 operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modelled with the Monte–Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Mechanisms in Regulating the Quasi‐Biennial Oscillation in Exascale Earth System Model Version 2

Abstract The update to the deep convection parameterization in the energy exascale Earth system model version 2 (E3SMv2) makes the simulated quasi‐biennial oscillation (QBO) have a shorter period than without the update when other tunable parameters the same since the update makes convection more intense but less frequent while leaving the time‐mean convective heating tendency almost unchanged. Amplitudes of momentum fluxes of parameterized gravity waves (GWs) are intensified since they are determined by the square of convective heating tendency. In addition, stronger planetary waves are simulated in E3SMv2 with the convection scheme update, partially due to the tropospheric precipitation change. Furthermore, there is evidence that planetary waves are intensified in the stratosphere due to the dissipation of enhanced parameterized GWs. These factors are found to be responsible for modulating the QBO simulated in E3SMv2.

Li, Yuanpu [Climate and Global Dynamics Division U↗
DOE OSTI2021

Development of an Ultrahigh-bandwidth Phase Contrast Imaging System for detection of electron scale turbulence and Gigahertz Radio-Frequency Waves

The study of waves and turbulence is vital to the development of future reactor-grade plasma devices developed in the quest for fusion energy. These fluctuations are responsible for moving heat and particles across the magnetic field, and a predictive understanding of them is needed to achieve the density and temperature required to sustain a plasma fusion reaction. While many techniques have been developed for measuring waves and fluctuations, every measurement method has limitations. There are relatively few techniques for measuring very high frequency fluctuations, such as radio frequency waves injected to heat the plasma, unstable waves driven by suprathermal particles, or short wavelength electrostatic waves driven by electron temperature or density gradients. The present project builds upon the proven phase contrast imaging (PCI) technique to extend the response of the diagnostic by orders of magnitude in frequency and almost a factor of ten in spatial resolution. PCI provides a measurement of electron density based on small angle scattering of a CO-2 laser beam by using optical techniques to render a phase shift as an intensity change on a detector. Due to the telecommunications revolution, technological development by manufacturers has focused on components in the near -infrared, so that high- quality lasers and detectors at 1.55 µm are readily available. Shifting PCI design to a new, shorter wavelength has numerous advantages and challenges. Similar detector performance is available with room-temperature arrays with GHz bandwidth, while the detector arrays for 10.6 µm required liquid nitrogen cooling and were therefore limited to a bandwidth of about 1 MHz. Shifting to a shorter wavelength reduces the angle at which the laser beam scatters off of plasma waves, allowing more such scattered components to pass through the aperture of the vacuum vessel port and be collected by the PCI, which increases the spatial resolution. Concomitant with these benefits, various questions of performance arise. At shorter laser wave- length, the sensitivity to mirror and lens quality is increased, the contribution of the laser to the overall system noise is increased, and the sensitivity to vibrations is increased. The custom optical components at the heart of the PCI technique were required to be properly scaled for the shorter wavelength, so fabrication technologies needed to be explored. This project was designed to show that a low noise, high response PCI system at 1.55 µm could be constructed and operated, and then to quantify potential issues to allow extrapolation to a full-size production system providing physics measurements on a large plasma device. The first stage, producing the custom optical component called a Phase Plate, was successfully achieved using two methods. First, an easily reproduced masking and coating technique was able to produce good phase plates with the required parameters. Second, a nanofabrication technique was found to produce extremely high quality phase plates at a competitive price. The PCI constructed with the new phase plates and 1.55 µm laser was found to provide excellent wavelength measurements with the theoretically expected response. The sensitivity to optical surface quality was found to be in line with previous measurements at 10.6 µm. The observed signal-to-noise ratio was similar to the theoretically expected value. The effect of vibrations on PCI was studied with the first measurement of the effect of beam motion on PCI response and comparison to theory, allowing for a quantitative prediction of the effect of vibrations on a production PCI system and the requirements for improved beam stabilization. PCI is an extremely cost-effective method to provide a low noise, absolutely calibrated measurement of plasma fluctuations across a wide spatial scale. This project has shown that a 1.55 µm PCI using modern techniques and components is less expensive than the 10.6 µm PCI of a few years ago, with the largest savings in phase plate fabrication and the infrared detector array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Diffraction of radio-frequency waves by plasma turbulence in the edge of a tokamak

The Mie-Lorenz-Debye theory of scattering of radio frequency waves (RF) by filamentary structures in the edge of a tokamak plasma has been established in prior publications. These results have been verified. Theory and computations show that RF can be reflected, refracted, and diffracted by turbulence. Furthermore, the filamentary structures lead to side scattering of waves and excitation of wave modes different from the launched RF. It is important to quantify the overall effect of turbulence on the propagation of RF since efficiency of heating and current drive by RF waves depends upon the wave characteristics. The commonly used RF ranges vary from the low frequency ion-cyclotron (IC) waves, to medium frequency lower-hybrid (LH) waves, to high-frequency electron cyclotron (EC) waves with wavelengths respectively, large, comparable, or smaller than the radial dimensions of the filaments. Currently, we discuss the scattering as well as the forces exerted by different RF in the edge plasma. Also, we present a code for 3D full wave scattering analysis of arbitrary tensor permittivity and plasma density profiles. The study is applicable to ITER-like plasmas, as well as to plasmas in medium sized tokamaks such as TCV, ASDEX-U and DIII-D and high field concepts like SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Thermal Conductivity Degradation in High Burnup U-Pu-Zr Fuel

Recent advancements in the characterization of irradiated U-Pu-Zr fuels have revealed complexities that challenge existing understanding of constituent redistribution. Traditionally, models have proposed three concentric regions within the fuel, each characterized by distinctive phases and porosity. However, through detailed analysis of high burnup U-Pu-Zr, we discovered the presence of four distinct constituent redistribution regions. Particularly novel is the observation of significant Pu redistribution, a previously unreported phenomenon that necessitates a reevaluation of current models. This work aims to delve deeper into these findings, seeking to correlate mesoscale measurements of thermal diffusivity and respective thermal conductivity with the phases present in each redistribution region. To achieve this objective, we employed mesoscale thermoreflectance methods using the unique, Idaho National Laboratory (INL) developed, Thermal Conductivity Microscope (TCM) at INL’s Irradiated Materials Characterization Laboratory. The TCM employs two tightly focused lasers: one for heating to generate periodic thermal waves in the substrate, and another spatially separated probe laser to detect changes in the optical reflectivity of the gold-coated substrate resulting from thermal wave diffusion. We conducted several thermal diffusivity measurements within each region of constituent redistribution of a U-19Pu-10Zr fuel pin cross section irradiated to 11 at. % burnup. The TCM measurement positions strategically aligned with transmission electron microscopy (TEM) lift-out locations previously collected from the fuel sample. Complementary microstructural analysis techniques such as optical and scanning electron microscopy (OM/SEM), electron probe microanalysis for chemical compositions, and TEM-based selective area electron diffraction (SAED) analysis for crystallographic insights into each phase were also utilized. This comprehensive approach allowed us to correlate local thermal diffusivity data with microstructural characteristics, enabling the computation of local thermal conductivity at each position. The significance of this contribution lies in its pioneering use of the TCM for ternary fuel mesoscale examination, shedding light on the previously overlooked effects of Pu redistribution on local thermal conductivity. By informing current models capturing constituent redistribution and heat transfer, our findings pave the way for more accurate predictions of metallic fuel performance. Moreover, this work sets the stage for future comparisons with similar TCM examinations on U-19Pu-10Zr fuels at ultra-low burnup, facilitating a comprehensive understanding of thermal property changes across different burnup levels. Ultimately, our study not only enriches our understanding of the thermophysical properties of individual redistribution regions within U-Pu-Zr fuel but also offers valuable insights for the design and operational parameters of proposed next-generation fast reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Microscopic Mechanisms of Glasslike Lattice Thermal Transport in Cubic Cu 12 Sb 4 S 13 Tetrahedrites

Materials based on cubic tetrahedrites (Cu 12 Sb 4 S 13 ) are useful thermoelectrics with unusual thermal and electrical transport properties, such as very low and nearly temperature-independent lattice thermal conductivity (κ L ). We explain the microscopic origin of the glasslike κ L in Cu 12 Sb 4 S 13 by explicitly treating anharmonicity up to quartic terms for both phonon energies and phonon scattering rates. We show that the strongly unstable phonon modes associated with trigonally coordinated Cu atoms are anharmonically stabilized above approximately 100 K and continue hardening with increasing temperature in accord with experimental data. This temperature-induced hardening effect reduces scattering of heat carrying acoustic modes by reducing the available phase space for three-phonon processes, thereby balancing the conventional ∝ T increase in scattering due to phonon population and yielding nearly temperature independent κ L . Furthermore, we find that very strong phonon broadening leads to a qualitative breakdown of the conventional phonon-gas model and modify the dominant heat transport mechanism from the particlelike phonon wave packet propagation to incoherent contributions described by the off-diagonal terms in the heat-flux operator, which are typically prevailing in glasses and disordered crystals. Furthermore, our work paves the way to a deeper understanding of glasslike thermal conductivity in complex crystals with strong anharmonicity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Parametric raytracing modeling for NSTX-U scenario development with high harmonic fast waves and neutral beam injection

High harmonic fast waves (HHFW) are a versatile heating and current drive tool for scenario development. Extensive modeling scans were performed to find optimal parameters for different uses of HHFW in National Spherical Tokamak Experiment (NSTX-U). Scans of plasma density, temperature, magnetic field, and antenna phasing were performed both with and without neutral beam injection. For speed of calculation, the ray-tracing code GENRAY coupled to the quasilinear Fokker–Planck code CQL3D was used. CQL3D allows for a more accurate description of the fast ion population, as well as for quasilinear effects such as HHFW-induced modifications of the distribution function. Best current drive results are obtained at elevated electron temperatures and with the lowest k φ phasing. Adding neutral beams however typically strongly reduces the HHFW current drive efficiency at the low density cases due to HHFW absorption on beam ions. Results of this parametric study will feed into scenario development and predict-first whole-shot modeling of NSTX-U discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Improved Simulation of East Asian Summer Monsoon in the High-resolution CESM1 and Its Causes

Based on the high- and low-resolution Community Earth System Model, version 1 (CESM1), and corresponding simulations from phase 6 of the Coupled Model Intercomparison Project (CMIP6), we compare the interannual variability of the East Asian summer monsoon (EASM). The EASM interannual variability is characterized by the anomalous western North Pacific anticyclone (WNPAC) circulation and the dipole rainfall pattern with a negative southern lobe over the western North Pacific and a positive northern lobe along the Meiyu–Baiu region, which is better reproduced by the high-resolution models. Here, the reason for the improvement in the high-resolution models has been attributed to the better simulation of the warm temperature advection from the wind anomalies on the climatological temperature gradient. Positive sea surface temperature (SST) anomalies over the tropical Indian Ocean are the key to the improved wind anomalies featuring a WNPAC in the high-resolution models. The warm SST anomalies over the tropical Indian Ocean strengthen the WNPAC by triggering a Kelvin-wave response to the enhanced heat release induced by the increased precipitation. Based on the mixed-layer heat budget analysis, the warm SST anomalies over the western Indian Ocean in the high-resolution CESM1 are tied to the anomalous easterly wind along the equator, which reduces surface evaporation and upwelling. Therefore, the better simulations of air-sea feedback and the oceanic mesoscale eddy over the western Indian Ocean are the key for the improved simulation of the EASM interannual variations in the high-resolution CESM1.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Physical, numerical, and computational challenges of modeling neutrino transport in core-collapse supernovae

The proposal that core collapse supernovae are neutrino driven is still the subject of active investigation more than 50 years after the seminal paper by Colgate and White. The modern version of this paradigm, which we owe to Wilson, proposes that the supernova shock wave is powered by neutrino heating, mediated by the absorption of electron-flavor neutrinos and antineutrinos emanating from the proto-neutron star surface, or neutrinosphere. Neutrino weak interactions with the stellar core fluid, the theory of which is still evolving, are flavor and energy dependent. The associated neutrino mean free paths extend over many orders of magnitude and are never always small relative to the stellar core radius. Thus, neutrinos are never always fluid like. Instead, a kinetic description of them in terms of distribution functions that determine the number density of neutrinos in the six-dimensional phase space of position, direction, and energy, for both neutrinos and antineutrinos of each flavor, or in terms of angular moments of these neutrino distributions that instead provide neutrino number densities in the four-dimensional phase-space subspace of position and energy, is needed. In turn, the computational challenge is twofold: (i) to map the kinetic equations governing the evolution of these distributions or moments onto discrete representations that are stable, accurate, and, perhaps most important, respect physical laws such as conservation of lepton number and energy and the Fermi–Dirac nature of neutrinos and (ii) to develop efficient, supercomputer-architecture-aware solution methods for the resultant nonlinear algebraic equations. In this review, we present the current state of the art in attempts to meet this challenge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Thermal Energy Transport in Oxide Nuclear Fuel

To efficiently capture the energy of the nuclear bond, advanced nuclear reactor concepts seek solid fuels that must withstand unprecedented temperature and radiation extremes. In these advanced fuels, thermal energy transport under irradiation is directly related to reactor performance as well as reactor safety. The science of thermal transport in nuclear fuel is a grand challenge as a result of both computational and experimental complexities. Here we provide a comprehensive review of thermal transport research on two actinide oxides: one currently in use in commercial nuclear reactors, uranium dioxide (UO 2 ), and one advanced fuel candidate material, thorium dioxide (ThO 2 ). In both materials, heat is carried by lattice waves or phonons. Crystalline defects caused by fission events effectively scatter phonons and lead to a degradation in fuel performance over time. Bolstered by new computational and experimental tools, researchers are now developing the foundational work necessary to accurately model and ultimately control thermal transport in advanced nuclear fuels. We begin by reviewing research aimed at understanding thermal transport in perfect single crystals. The absence of defects enables studies that focus on the fundamental aspects of phonon transport. Next, we review research that targets defect generation and evolution. Here the focus is on ion irradiation studies used as surrogates for damage caused by fission products. We end this review with a discussion of modeling and experimental efforts directed at predicting and validating mesoscale thermal transport in the presence of irradiation defects. While efforts in these research areas have been robust, challenging work remains in developing holistic tools to capture and predict thermal energy transport across widely varying environmental conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Effect of Laser Beam Profile on Rotating Lattice Single Crystal Growth in Sb2S3 Model Glass

Laser heating of chalcogenide glasses has successfully produced rotating lattice single crystals through a solid-solid transformation. To better understand the nature of complex, orientation-dependent lattice rotation, we designed heat profiles of the continuous wave laser by beam shaping, fabricated larger Sb2S3 crystal dots in Sb2S3 glass, and investigated the lattice rotation where the crystal could grow in all radial directions under a circular thermal gradient. The results show that the rate of lattice rotation is highly anisotropic and depends on crystallographic direction. The nature of this rotation is the same in crystals of different orientation relative to the surface. The growth directions that align with the slip planes show the highest rate of rotation and the rotation rate gradually decreases away from this direction. Additionally, the presence of multiple growth directions results in a complicated rotation system. We suggest that the growth front influences the density of dislocations introduced during growth under confinement and thus affects the lattice rotation rate in these crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reflection and transmission of electromagnetic pulses at a planar dielectric interface -- theory and quantum lattice simulations

There is considerable interest in the application of quantum information science to advance computations in plasma physics. A particular point of curiosity is whether it is possible to take advantage of quantum computers to speed up numerical simulations relative to conventional computers. Many of the topics in fusion plasma physics are classical in nature. In order to implement them on quantum computers it will require couching a classical problem in the language of quantum mechanics. Electromagnetic waves are routinely used in fusion experiments to heat a plasma or to generate currents in the plasma. The propagation of electromagnetic waves is described by Maxwell equations with an appropriate description of the plasma as a dielectric medium. Before advancing to the tensor dielectric of a magnetized plasma, this paper considers electromagnetic wave propagation in a one-dimensional inhomogeneous scalar dielectric. The classic theory of scattering of plane electromagnetic waves at a planar interface, separating two different dielectric media, leads to Fresnel equations for reflection and transmission coefficients. In contrast to plane waves, this paper is on the reflection and transmission of a spatially confined electromagnetic pulse. Following an analytical formulation for the scattering of a Gaussian pulse, it is deduced that the maximum transmission coefficient for a pulse is $\sqrt{n_2/n_1}$ times that for a plane wave; the incident and transmitted pulses propagating in dielectric media with refractive indices $n_1$ and $n_2$, respectively. The analytical theory is complemented by numerical simulations using a quantum lattice algorithm for Maxwell equations. The algorithm, based on the Riemann-Silberstein-Weber representation of the electromagnetic fields and expressed in term of qubits, is an interleaved sequence of entangling operators at each lattice site and unitary streaming operators which transmit information from one site to an adjacent lattice site. Besides substantiating results from the theory for Gaussian pulses, numerical simulations show their validity for non-Gaussian pulses. Apart from their time-asymptotic forms, the simulations display an interplay between the incident, reflected, and transmitted pulses in the vicinity of the transition region between two dielectric media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Real-time laser ultrasonic monitoring of laser-induced thermal processes

Abstract Intra- and inter-layer integrity of components fabricated with advanced manufacturing techniques, such as laser powder bed fusion, is dependent upon rapid heating, melting, and solidification processes. There is a need for new techniques to provide in situ feedback of these processes. Here a laser-based ultrasonic technique to probe thermal effects induced by a high-power continuous wave laser in titanium samples is described. Numerical simulations were performed to show that, for a spatially uniform heating beam, laser-induced surface acoustic waves are strongly influenced by surface heating conditions, are dispersive in the case of rapid heating, and that an abrupt velocity reduction happens upon the onset of surface melting. Furthermore, laser-based ultrasound experimental results which monitor the transient change of surface wave travel time associated with high power laser surface heating are provided. A pulsed laser is used to generate high frequency surface acoustic waves that propagate through the laser-heated region and are detected using a photorefractive crystal-based interferometer. Qualitative agreement is observed between theory and experiment with both showing a rapid reduction in the surface wave velocity at the onset of illumination and further decrease in surface wave velocity associated with melting. It is demonstrated that changes in the surface wave velocity can be used to track local heating and detect the onset of surface melting in real time.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Model predictive current profile control in tokamaks by exploiting spatially moving electron cyclotron current drives

The plasma control systems in next-generation tokamaks like ITER will balance competing control objectives to achieve the desired level of performance in advanced scenarios while preventing magnetohydrodynamic instabilities and disruptions. During normal tokamak operation, the points of incidence of the electromagnetic waves generated by the electron cyclotron heating and current drives (EC H&CDs) are usually fixed in space. However, the points of incidence can be modified in real-time by changing the angles of the mirrors that reflect the EC H&CD waves. Altering the points of incidence, in turn, varies the ability of the plasma control system to regulate a plasma property. For instance, changing the EC H&CD wave incidence location may place the power demands necessary to achieve a particular plasma target within saturation limits. Therefore, using the EC H&CD deposition location, which is related to the EC H&CD mirror angle, as a supplementary controllable variable may facilitate access to a given target scenario. However, active scenario-control algorithms have not been designed so far to fully exploit this capability in real time. In this work, a model predictive controller that can handle actuation locations as control inputs is developed. In particular, the controller is designed to regulate both the auxiliary powers and the EC H&CD deposition locations in a pre-defined optimal sense to achieve the control objective of attaining and sustaining a target current profile. Further, the proposed controller is tested for a DIII-D tokamak scenario in nonlinear simulations using the Control Oriented Transport SIMulator (COTSIM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗