Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave data”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

VLBI Data Processing on Coronal Radio-sounding Experiments of Mars Express

The ESA’s Mars Express solar corona experiments were performed at two solar conjunctions in the years 2015 and 2017 by a number of radio telescopes in the European VLBI Network. This paper presents the methods to measure the frequency and phase fluctuations of the spacecraft radio signal, and the applications to study the characteristics of the plasma turbulence effects on the signal at a single station and at multiple stations via cross correlation. The power spectra of the frequency fluctuations observed between 4.9 and 76.3 R {sub s} have a power-law shape close to a Kolmogorov spectrum over the frequency interval ν {sub lo} < ν < ν {sub up}, where the nominal value of ν {sub lo} is set to 3 mHz and ν {sub up} is in the range of 0.03–0.15 Hz. The rms of the frequency fluctuations is presented as a function of the heliocentric distance. Furthermore, we analyze the variations of the electron column density fluctuations at solar offsets 4.9 R {sub s} and 9.9 R {sub s} and the cross-correlation products between the VLBI stations. The power density of the differential fluctuations between different stations decreases at ν < 0.01 Hz. Finally, the fast flow speeds of solar wind >700 km s{sup -1} are derived from the cross correlation of frequency fluctuations at ν < 0.01 Hz. The fast flow speeds of solar wind correspond to the high heliolatitude of the coronal region that the radio rays passed. The VLBI observations and analysis methods can be used to study the electron column density fluctuations and the turbulence at multiple spatial points in the inner solar wind by providing multiple lines of sight between the Earth and the spacecraft.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Radiation Detection System Integration Using Multiple-Input Multiple-Output Radios

This study describes the application development of multiple-input multiple-output radios to provide persistent mobile ad hoc network (MANET) for the Department of Homeland Security. By using Man Portable Unit (MPU5) fifth generation radios (manufactured by Persistent Systems) with the Android Team Awareness Kit (ATAK), an Android smartphone geospatial infrastructure and military situational awareness application, the Remote Sensing Laboratory has developed a MANET connectivity to monitor deployed nuclear/radiological search operation assets.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Possible $ν$ Source Class: 3-sigma Detection of High-Energy Neutrinos from Supermassive Black Hole Binary Candidates

Identifying the sources of high-energy (TeV-PeV) astrophysical neutrinos is crucial for studies in both astrophysics and particle physics. Despite extensive searches for more than a decade, which revealed several individual potential sources and only one potential source class, the origins of these neutrinos remain largely unresolved; thus, more source classes should be investigated. In this work, we conduct the first search for high-energy neutrino emission from a new source class, supermassive black hole binaries (SMBHBs), which are also theoretically motivated. We perform an unbinned maximum-likelihood-ratio analysis on our constructed catalog of 693 SMBHB candidates and 10 years of IceCube public data. Our results show positive correlations, with higher significance in more physically motivated scenarios and the highest significance at 3.0$σ$. In addition, we also study potential connections between SMBHBs' high-energy neutrino and nano-Hz gravitational-wave emissions, the latter being the main target of pulsar timing arrays. Our results provide the first evidence of SMBHBs being high-energy neutrino emitters.

Pugazhendhi A. D, Pugazhendhi D [Bangalore, Indian↗
DOE OSTI2020

Characterization and Scaling of Forced Convective Swirl in Sinusoidal Wavy-Plate-Fin Cores of Compact Heat Exchangers

The characterization and scaling of the thermal-hydraulic performance in wavy plate-fin compact heat exchanger cores, based on the understanding of the physical phenomena and heat transfer enhancement mechanism is delineated. Experimental data are presented for forced convection in air (Pr = 0.71) with flow rates in the range 50 ≤ Re ≤ 4000. A variety of wavy-fin cores that span viable applications, with geometrical attributes described by the cross section aspect ratio α (ratio of fin spacing to height), fin corrugation aspect ratio γ (ratio of 2× corrugation amplitude to wave pitch), and fin spacing ratio ζ (ratio of fin spacing to wave pitch), are considered. To characterize and correlate the vortex-flow mixing in interfin spaces, a Swirl number Sw is introduced from the balance of viscous, inertial and centrifugal forces. It is shown that the laminar, transitional and turbulent flow regimes can be explicitly identified by this Swirl number. Based on the experimental results and extended analysis, new correlations for Fanning friction factor f and Colburn factor j are developed with Sw, α, γ, and ζ as scaling parameters. Requisite expressions are devised by a superposition of both enhancement components due to the corrugated surface area enlargement and induced swirl flow field, and they are combined to cover the laminar, transitional and turbulent regimes by the method of asymptotic matching. The resulting generalized correlations are further shown to predict all available experimental data for f and j factors to within ±20% and ±15%, respectively.

42 ENGINEERING↗
DOE OSTI2022

A mass–momentum consistent coupling for mesh-adaptive two-phase flow simulations

Here, we present a novel mass-momentum consistent coupling between a geometric volume-of-fluid scheme and an incompressible flow solver with differing directional-splitting approaches. The advection of the volume fraction is performed using a direction-split algorithm, whereas the momentum advection algorithm uses a traditional unsplit, fractional-step approach. Both algorithms employ finite-volume discretizations based on Cartesian meshes. In solving the mass-momentum consistency problem, momentum fluxes at the cell faces are weighted by the density fluxes based on the already advected volume fraction. The success of our approach lies on introducing a Favre-averaged velocity interpolation at the liquid/gas interface along with a minmod slope limiter. Mesh-convergence studies show that when the minmod slope limiter is used, the two-phase solver retains an accuracy between first and second order, but when a purely upwind scheme is considered, its accuracy drops to first order. Finally, after considering several validation problems, the solver is shown to agree well with reference numerical and experimental data while retaining its robustness and efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quantum annealing algorithms for Boolean tensor networks

Abstract Quantum annealers manufactured by D-Wave Systems, Inc., are computational devices capable of finding high-quality heuristic solutions of NP-hard problems. In this contribution, we explore the potential and effectiveness of such quantum annealers for computing Boolean tensor networks. Tensors offer a natural way to model high-dimensional data commonplace in many scientific fields, and representing a binary tensor as a Boolean tensor network is the task of expressing a tensor containing categorical (i.e., $$\{0, 1\}$$ { 0 , 1 } ) values as a product of low dimensional binary tensors. A Boolean tensor network is computed by Boolean tensor decomposition, and it is usually not exact. The aim of such decomposition is to minimize the given distance measure between the high-dimensional input tensor and the product of lower-dimensional (usually three-dimensional) tensors and matrices representing the tensor network. In this paper, we introduce and analyze three general algorithms for Boolean tensor networks: Tucker, Tensor Train, and Hierarchical Tucker networks. The computation of a Boolean tensor network is reduced to a sequence of Boolean matrix factorizations, which we show can be expressed as a quadratic unconstrained binary optimization problem suitable for solving on a quantum annealer. By using a novel method we introduce called parallel quantum annealing, we demonstrate that Boolean tensor’s with up to millions of elements can be decomposed efficiently using a DWave 2000Q quantum annealer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗
DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗
DOE OSTI2022

Probing physics beyond the standard model: limits from BBN and the CMB independently and combined

We present new Big Bang Nucleosynthesis (BBN) limits on the cosmic expansion rate or relativistic energy density, quantified via the number N ν of equivalent neutrino species. We use the latest light element observations, neutron mean lifetime, and update our evaluation for the nuclear rates d + d → 3 He + n and d + d → 3 H+p. Combining this result with the independent constraints from the cosmic microwave background (CMB) yields tight limits on new physics that perturbs N ν and η prior to cosmic nucleosynthesis: a joint BBN+CMB analysis gives N ν = 2.898 ± 0.141, resulting in N ν < 3.180 at 2σ. We apply these limits to a wide variety of new physics scenarios including right-handed neutrinos, dark radiation, and a stochastic gravitational wave background. The strength of the independent BBN and CMB constraints now opens a new window: we can search for limits on potential changes in N ν and/or the baryon-to-photon ratio η between the two epochs. The present data place strong constraints on the allowed changes in N ν between BBN and CMB decoupling; for example, we find -0.708 < N ν CMB - N ν BBN < 0.328 in the case where η and the primordial helium mass fraction Yp are unchanged between the two epochs; we also give limits on the allowed variations in η or in (η, N ν ) jointly. We discuss scenarios in which such changes could occur, and show that BBN+CMB results combine to place important constraints on some early dark energy models to explain the H 0 tension. Looking to the future, we forecast the tightened precision for N ν arising from both CMB Stage 4 measurements as well as improvements in astronomical 4 He measurements. Here, we find that CMB-S4 combined with present BBN and light element observation precision can give σ(N ν ) ≃ 0.03. Such future precision would reveal the expected effect of neutrino heating (N eff -3 = 0.044) of the CMB during BBN, and would be near the level to reveal any particle species ever in thermal equilibrium with the standard model. Improved Y p measurements can push this precision even further.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Triaxiality and the nature of low-energy excitations in Ge 76

The deformation properties of the low-lying states in 76 Ge have been investigated following a safe energy Coulomb excitation measurement with the GRETINA tracking array and CHICO2 heavy-ion counter at the ATLAS accelerator facility at Argonne National Laboratory. A comprehensive set of transition and static E2 matrix elements were extracted from the measured differential Coulomb cross-sections, and compared with results of configuration interaction shell-model calculations and computations carried out within the framework of the generalized triaxial rotor model. The remarkable agreement between the calculated and experimental data supports a near-maximum triaxial deformation for the ground state of 76 Ge. Additionally, the degree of softness of the asymmetry in 76 Ge and 76 Se was investigated using rotational invariants generated from configuration interaction shell-model wave functions computed with the jj44b and JUN45 effective interactions. The resulting invariants are shown to be consistent with a stiff triaxial deformation in 76 Ge and a predominantly soft triaxial potential for 76 Se, in agreement with the conclusions of recent works by this collaboration.

59 ≤ A ≤ 89↗
DOE OSTI2020

Forward sensitivity approach for estimating eddy viscosity closures in nonlinear model reduction

In this paper, we propose a variational approach to estimate eddy viscosity using forward sensitivity method (FSM) for closure modeling in nonlinear reduced order models. FSM is a data assimilation technique that blends model's predictions with noisy observations to correct initial state and/or model parameters. We apply this approach on a projection based reduced order model (ROM) of the one-dimensional viscous Burgers equation with a square wave defining a moving shock, and the two-dimensional vorticity transport equation formulating a decay of Kraichnan turbulence. We investigate the capability of the approach to approximate an optimal value for eddy viscosity with different measurement configurations. Specifically, we show that our approach can sufficiently assimilate information either through full field or sparse noisy measurements to estimate eddy viscosity closure to cure standard Galerkin reduced order model (GROM) predictions. Therefore, our approach provides a modular framework to correct forecasting error from a sparse observational network on a latent space. We highlight that the proposed GROM-FSM framework is promising for emerging digital twin applications, where real-time sensor measurements can be used to update and optimize surrogate model's parameters.

42 ENGINEERING↗
DOE OSTI2020

Observation of Radar Echoes from High-Energy Particle Cascades

We report here the observation of radar echoes from the ionization trails of high-energy particle cascades. Data were taken at the SLAC National Accelerator Laboratory, where the full electron beam ( ~ 10 9 e - at ~ 10 GeV / e - ) was directed into a plastic target to simulate an ultrahigh-energy neutrino interaction. The target was interrogated with radio waves, and coherent radio reflections from the cascades were detected with properties consistent with theoretical expectations. This is the first definitive observation of radar echoes from high-energy particle cascades, which could lead to a viable neutrino detection technology for energies ≳ 10 16 eV .

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

A Study on Modeled Wind Speed Errors Using the U.S. Department of Energy Buoys

Pacific Northwest National Laboratory (PNNL) operates two AXYS WindSentinel lidar buoys for the U.S. Department of Energy’s Wind Energy Technologies Office. The purpose of these buoys is to collect hub-height winds and supporting meteorological and oceanographic information to facilitate the development of wind energy in the U.S. waters. The first deployment for one buoy was off the coast of Virginia from December 2014 to May 2016, and the first deployment for the other buoy was off the coast of New Jersey from November 2015 until February 2017. This report describes recent analysis of data collected during these first two deployments. Specifically, we compare hub-height wind speed estimates using Monin-Obukhov Similarity Theory (MOST) to the lidar measurements, and examine how those errors are affected by wind direction, atmospheric stability, wind-wave direction differences, and various measures of the wave-state. The comparisons are done using standard similarity functions based on MOST; including the Businger - Dyer, the Beljaars & Holtslag and the Vickers & Mahrt similarity functions. All models produce large errors over the range of atmospheric stabilities that were observed, with the largest errors occurring for stable flows. The Vickers & Mahrt function resulted in the largest overall bias and standard deviation, while Beljaars & Holtslag function gave the smallest bias and standard deviation due to its better performance under stable conditions. The models perform best under unstable conditions, but even in this regime there is a consistent overestimation of the wind speed of between roughly 0 to 1 ms -1 compared to the lidar measurements. We identify specific metocean conditions (i.e. stability and wind and wave directions) at each of the deployment locations that lead to large errors in MOST predictions. Finally, a coupled ocean-atmosphere model framework was investigated to simulate large errors in weather research forecasting (WRF).

17 WIND ENERGY↗
DOE OSTI2022

Deep Learning Based Superconducting Radio-Frequency Cavity Fault Classification at Jefferson Laboratory

This work investigates the efficacy of deep learning (DL) for classifying C100 superconducting radio-frequency (SRF) cavity faults in the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab. CEBAF is a large, high-power continuous wave recirculating linac that utilizes 418 SRF cavities to accelerate electrons up to 12 GeV. Recent upgrades to CEBAF include installation of 11 new cryomodules (88 cavities) equipped with a low-level RF system that records RF time-series data from each cavity at the onset of an RF failure. Typically, subject matter experts (SME) analyze this data to determine the fault type and identify the cavity of origin. This information is subsequently utilized to identify failure trends and to implement corrective measures on the offending cavity. Manual inspection of large-scale, time-series data, generated by frequent system failures is tedious and time consuming, and thereby motivates the use of machine learning (ML) to automate the task. This study extends work on a previously developed system based on traditional ML methods (Tennant and Carpenter and Powers and Shabalina Solopova and Vidyaratne and Iftekharuddin, Phys. Rev. Accel. Beams, 2020, 23, 114601), and investigates the effectiveness of deep learning approaches. The transition to a DL model is driven by the goal of developing a system with sufficiently fast inference that it could be used to predict a fault event and take actionable information before the onset (on the order of a few hundred milliseconds). Because features are learned, rather than explicitly computed, DL offers a potential advantage over traditional ML. Specifically, two seminal DL architecture types are explored: deep recurrent neural networks (RNN) and deep convolutional neural networks (CNN). We provide a detailed analysis on the performance of individual models using an RF waveform dataset built from past operational runs of CEBAF. In particular, the performance of RNN models incorporating long short-term memory (LSTM) are analyzed along with the CNN performance. Furthermore, comparing these DL models with a state-of-the-art fault ML model shows that DL architectures obtain similar performance for cavity identification, do not perform quite as well for fault classification, but provide an advantage in inference speed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Ambient field seismology in critical zone hydrological sciences

Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.

58 GEOSCIENCES↗
DOE OSTI2019

Equation of state constraints from multi-messenger observations of neutron star mergers

The very first detection of gravitational waves from a neutron star binary merger, GW170817, exceeded all expectations. The event was relatively nearby, which may point to a relatively high merger rate. It was possible to extract finite-size effects from the gravitational-wave signal, which constrains the nuclear equation of state. Also, an electromagnetic counterpart was detected at many wavebands from radio to gamma rays marking the begin of a new multi-messenger era involving gravitational waves. We describe how multi-messenger observations of GW170817 are employed to constrain the nuclear equation of state. Combining the information from the optical emission and the mass measurement through gravitational waves leads to a lower limit on neutron star radii. According to this conservative analysis, which employs a minimum set of assumptions, the radii of neutron stars with typical masses should be larger than about 10.7 km. This implies a lower limit on the tidal deformability of about 210, while much stronger lower bounds are not supported by the data of GW170817. The multi-messenger interpretation of GW170817 rules out very soft nuclear matter and complements the upper bounds on NS radii which are derived from the measurement of finite-size effects during the pre-merger phase. We highlight the future potential of multi-messenger observations and of GW measurements of the postmerger phase for constraining the nuclear equation of state. Finally, we propose an observing strategy to maximize the scientific yield of future multi-messenger observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The São Paulo Potential and the 3 He Breakup Reaction at 130 MeV on 93 Nb and 197 Au

Nuclear reactions induced by weakly bound nuclei often result in the emission of light particles as products of the projectile breakup. The incident particle is expected to fragment before reaching the interior of the target nuclei, a characteristic that tends to favor forward angular emission. Although the first models developed for such reactions were proposed many decades ago by Udagawa and Tamura, Ichimura, Austern, Vincent, and Kasano, and Hussein and McVoy, to this date most of the attention has been dedicated to the breakup of the deuteron and other very weakly bound systems, while more strongly bound projectiles have not been fully explored. Here, we describe the production of protons and deuterons from the breakup of 3He. The breakup of 3 He on two heavy targets, 93 Nb and 197 Au, is analyzed in this study. However, the incoming energy involved is much larger than that of the standard phenomenological optical potentials available in the literature for the 3He entrance channel. We have overcome this difficulty by using the São Paulo potential and adjusting the nuclear diffusivity. The deuteron inclusive spectra were calculated at several angles and compared well with the experimental data. Theoretical proton emission predictions are also given for future reference, since no inclusive measurements were performed for the targets under study. Finally, one of our goals is to verify the description of deuteron emission from fast projectiles, for which many partial waves contribute to the scattering process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Benchmarking magnetised three-wave coupling for laser backscattering: analytic solutions and kinetic simulations

Understanding magnetised laser–plasma interactions is important for controlling magneto-inertial fusion experiments and developing magnetically assisted radiation and particle sources. For nanosecond pulses at non-relativistic intensities, interactions are dominated by coherent three-wave interactions, whose nonlinear coupling coefficients became known only recently when waves propagate at oblique angles with the magnetic field. In this paper, backscattering coupling coefficients predicted by warm-fluid theory are benchmarked using particle-in-cell simulations in one spatial dimension, and excellent agreements are found for a wide range of plasma temperatures, magnetic field strengths and laser propagation angles, when the interactions are mediated by electron-dominant hybrid waves. Systematic comparisons between theory and simulations are made possible by a rigorous protocol. On the theory side, the initial boundary value problem of linearised three-wave equations is solved, and the transient-time solutions allow the effects of growth and damping to be distinguished. On the simulation side, parameters are carefully chosen and calibration runs are performed to ensure that comparisons are well controlled. Fitting simulation data to analytical solutions yields numerical growth rates that match theory predictions within error bars. Although warm-fluid theory is found to be valid for a wide parameter range, genuine kinetic effects have also been observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗