Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear waves”

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 523 records · Page 29

The benefits of CMB delensing

The effects of gravitational lensing of the cosmic microwave background (CMB) have been measured at high significance with existing data and will be measured even more precisely in future surveys. Reversing the effects of lensing on the observed CMB temperature and polarization maps provides a variety of benefits. Delensed CMB spectra have sharper acoustic peaks and more prominent damping tails, allowing for improved inferences of cosmological parameters that impact those features. Delensing reduces B-mode power, aiding the search for primordial gravitational waves and allowing for lower variance reconstruction of lensing and other sources of secondary CMB anisotropies. Lensing-induced power spectrum covariances are reduced by delensing, simplifying analyses and improving constraints on primordial non-Gaussianities. Biases that result from incorrectly modeling nonlinear and baryonic feedback effects on the lensing power spectrum are mitigated by delensing. All of these benefits are possible without any changes to experimental or survey design. Here, we develop a self-consistent, iterative, all-orders treatment of CMB delensing on the curved sky and demonstrate the impact that delensing will have with future surveys.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electrostatic approximation

The plasma weak turbulence theory is a perturbative nonlinear theory, which has been proven to be quite valid in a number of applications. However, the standard weak turbulence theory found in the literature is fully developed for highly idealized unmagnetized plasmas. As many plasmas found in nature and laboratory are immersed in a background static magnetic field, it is necessary to extend the existing discussions to include the effects of ambient magnetic field. Such a task is quite formidable, however, which has prevented fundamental and significant progresses in the subject matter. The central difficulty lies in the formulation of the complete nonlinear response functions for magnetized plasmas. The present paper derives the nonlinear susceptibilities for weakly turbulent magnetized plasmas up to the third order nonlinearity, but in doing so, a substantial reduction in mathematical complexity is achieved by the use of Bessel function addition theorem (or sum rule). The present paper also constructs the weak turbulence wave kinetic equation in a formal sense. For the sake of simplicity, however, the present paper assumes the electrostatic interaction among plasma particles. Fully electromagnetic generalization is a subject of a subsequent paper.

Physics↗

Deep Koopman Neural Network for Analyzing High-Energy-Density Simulations of Electrical Wire Explosions

Megaampere-scale electrical wire experiments (EWEs) provide a platform for studying magnetohydrodynamic (MHD) instability growth in magneto-inertial fusion (MIF) devices. Even when nonlinear simulations of these experiments can digitally reproduce much of the experimentally observed instability growth, interpreting the results and understanding mode growth and evolution can be non-trivial. As a first step toward providing better interpretation of these simulation features, this work investigates the use of a deep neural network that uses Koopman operator theory to analyze the dynamics of pulsed-power-driven explosions of EWEs. This deep neural network is trained on 1-D resistive MHD simulations of EWEs. This neural network learns to transform the nonlinear data into a lower-dimensional representation where the time dynamics are linear. Layers of this neural network are shown to learn features of the simulations, including the locations of shock waves and different physical regimes of the simulation. Using the learned features, the network can compress a time state of the simulation consisting of 5120 data point into a 36-parameter lower-dimensional latent space embedding. Furthermore, these embeddings are shown to be clustered in the latent space by initial radius and time state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

OC6 phase I: Improvements to the OpenFAST predictions of nonlinear, low-frequency responses of a floating offshore wind turbine platform

We report in the OC5 and OC6 projects, the authors observed a persistent underprediction of the nonlinear, low-frequency responses of an offshore wind semisubmersible with many mid-fidelity engineering models, including the OpenFAST tool developed by the National Renewable Energy Laboratory. Both the low-frequency wave excitation in surge and pitch and the resulting resonance motions were severely underpredicted. In response, we developed several modifications to the OpenFAST model from the OC5/6 projects to improve the predictions of the low-frequency wave loads and responses. All modifications are in the modeling of the viscous drag forces. Efforts were made to provide physical justifications to the changes and to limit the number of additional parameters requiring tuning, so that the modified model can be applied to other floating wind systems. With the proposed modifications, the predictions of the low-frequency surge and pitch wave loads on a fixed floater and the resonance responses of a floating structure are both significantly improved with a single set of model coefficients, which leads to good agreement with the measurements from the OC6 wave-basin experimental campaign.

17 WIND ENERGY↗

Single-particle dynamics of microbunching

The mechanism of steady-state microbunching (SSMB) has been proposed [D. F. Ratner and A. W. Chao, Phys. Rev. Lett. 105, 154801 (2010)] to generate high-power coherent radiation at a high repetition rate or in continuous-wave mode using electron storage rings. In this paper, the related single-particle dynamics are theoretically and numerically studied, and important results are presented. The investigated effects are longitudinal quantum radiation excitation, nonlinear momentum compaction, and linear and nonlinear coupling of the transverse and longitudinal motion. Although this analysis is oriented toward SSMB, some of the analyzed effects are also crucial in cases such as coherent harmonic generation, bunch slicing, bunch compression, free-electron laser beam transport lines, and quasi-isochronous rings, which involve precise longitudinal phase-space manipulations.

43 PARTICLE ACCELERATORS↗

Impact of Biofouling on Point Absorber Wave Energy Converter Performance and Control

Biofouling is a well-documented problem in naval engineering, but little is known about its effect on wave energy converter (WEC) performance. In this study, the software WEC-Sim is used to simulate the performance of a point absorber WEC that has been biofouled by “hard” species (e.g., mussels, barnacles) to varying degrees. Specifically, biofouling is assumed to change the nonlinear drag forces acting on the WEC, which have quantifiable effects on key performance characteristics such as optimal damping conditions, power, peak displacement, and peak velocity. The results of this analysis are then used to discuss strategies for WEC control as it relates to biofouling. Furthermore, the results show that average power production can decrease by as much as 15% with heavy biofouling and require an adjustment of the optimal control law by up to 20%.

Biofouling↗

Evidence for a thermally driven charge-density-wave transition in 1T-TaS 2 thin-film devices: Prospects for GHz switching speed

We report on the room-temperature switching of 1T-TaS 2 thin-film charge-density-wave devices, using nanosecond-duration electrical pulsing to construct their time-resolved current–voltage characteristics. The switching action is based upon the nearly commensurate to incommensurate charge-density-wave phase transition in this material, which has a characteristic temperature of 350 K at thermal equilibrium. For sufficiently short pulses, with rise times in the nanosecond range, self-heating of the devices is suppressed, and their current–voltage characteristics are weakly nonlinear and free of hysteresis. This changes as the pulse duration is increased to ~200 ns, where the current develops pronounced hysteresis that evolves nonmonotonically with the pulse duration. By combining the results of our experiments with a numerical analysis of transient heat diffusion in these devices, we clearly reveal the thermal origins of their switching. In spite of this thermal character, our modeling suggests that suitable reduction of the size of these devices should allow their operation at GHz frequencies.

36 MATERIALS SCIENCE↗

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulating plasma wave propagation on a superconducting quantum chip

Quantum computers may one day enable the efficient simulation of strongly coupled plasmas that lie beyond the reach of classical computation in regimes where quantum effects are important and the scale separation is large. Here, in this article, we take a first step toward efficient simulation of quantum plasmas by demonstrating linear plasma wave propagation on a superconducting quantum chip. Using high-fidelity and highly expressive device-native gates, combined with an error-mitigation technique, we simulate the scattering of laser pulses from inhomogeneous plasmas. Our approach is made feasible by the identification of a suitable local spin model whose excitations mimic plasma waves, and whose circuit implementation requires a lower gate count than other proposed approaches that would require a future fault-tolerant quantum computer. This work opens avenues to study more complicated phenomena that cannot be simulated efficiently on classical computers, such as nonlinear quantum dynamics when strongly coupled plasmas are driven out of equilibrium.

general physics↗

On applications of quantum computing to plasma simulations

Quantum computing is gaining increased attention as a potential way to speed up simulations of physical systems, and it is also of interest to apply it to simulations of classical plasmas. However, quantum information science is traditionally aimed at modeling linear Hamiltonian systems of a particular form that is found in quantum mechanics, so extending the existing results to plasma applications remains a challenge. Here, we report a preliminary exploration of the long-term opportunities and likely obstacles in this area. First, we show that many plasma-wave problems are naturally representable in a quantumlike form and thus are naturally fit for quantum computers. Second, we consider more general plasma problems that include non-Hermitian dynamics (instabilities, irreversible dissipation) and nonlinearities. We show that by extending the configuration space, such systems can also be represented in a quantumlike form and thus can be simulated with quantum computers too, albeit that requires more computational resources compared to the first case. Third, we outline potential applications of hybrid quantum–classical computers, which include analysis of global eigenmodes and also an alternative approach to nonlinear simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimization of Superconducting Linac for Proton Improvement Plan-II (PIP-II)

PIP-II is an essential upgrade of the Fermilab complex that will enable the worlds most intense high-energy beam of neutrinos for the international Deep Underground Neutrino Experiment at LBNF and support a broad physics program at Fermilab. Ultimately, the PIP-II superconducting linac will be capable of accelerating the $H-$ CW beam to 800 MeV with an average power of 1.6 MW. To operate the linac with such high power, beam losses and beam emittance growth must be tightly controlled. In this paper, we present the results of global optimization of the Linac options towards a robust and efficient physics design for the superconducting section of the PIP-II linac. We also investigate the impact of the nonlinear field of the dipole correctors on the beam quality and derive the requirement on the field quality using statistical analysis. Finally, we assess the need to correct the quadrupole focusing produced by Half Wave, and Single Spoke accelerating cavities. We assess the feasibility of controlling the beam coupling in the machine by changing the polarity of the field of linac focusing solenoids

43 PARTICLE ACCELERATORS↗

Modulation of the Eastern Equatorial Pacific Seasonal Cycle by Tropical Instability Waves

Feedbacks from tropical instability waves (TIWs) on the seasonal cycle of the eastern Pacific Ocean are studied using two eddy-rich ocean simulations, with and without TIWs. By warming the equatorial waters by up to 0.4°C through nonlinear advection in boreal summer and fall, TIWs reduce the amplitude of the seasonal cycle in upper ocean temperatures. In addition, TIWs stabilize the upper part of the Equatorial Undercurrent (EUC) through enhanced barotropic energy conversion, leading to a year-round weakening by -0.15 m s -1 and preventing an unrealistic re-intensification in boreal fall usually found in non-eddy resolving models. A coarser simulation at 1-degree horizontal resolution fails to reproduce the TIW-induced nonlinear warming of equatorial waters, but succeeds in inhibiting the EUC re-intensification. This suggests a threshold effect in TIW strength, associated with the model's ability to simulate eddies, which may be responsible for long-standing biases displayed by global climate models in this region.

58 GEOSCIENCES↗

Pulsed RF schemes for tearing mode stabilization

The RF stabilization of tearing modes with current condensation has the potential to increase stabilization efficiency and loosen power localization requirements. Such benefits stem from the cooperative feedback between the RF deposition and resulting island temperature perturbation governed by diffusion. A self consistent treatment of the damping of an rf ray as it traverses the island shows that low damping scenarios can require unfavorably high powers to overcome initial power leakage and effectively capitalize on the nonlinear effect. In this work it is demonstrated that for such regimes,modulated stabilization schemes can achieve significant improvements in heating and current drive contributions to stabilization for the same average power as a continuous wave scheme. The impact of modulation frequency and duty cycle on the performance is explored, the results of which suggest modulation strategies in which the pulsing periods are kept on the order of a diffusive time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Giant nonlinear optical responses from photon-avalanching nanoparticles

Avalanche phenomena use steeply nonlinear dynamics to generate disproportionately large responses from small perturbations, and are found in a multitude of events and materials. Photon avalanching enables technologies such as optical phase-conjugate imaging, infrared quantum counting and efficient upconverted lasing. However, the photon-avalanching mechanism underlying these optical applications has been observed only in bulk materials and aggregates, limiting its utility and impact. In this study we report the realization of photon avalanching at room temperature in single nanostructures-small, Tm 3+ -doped upconverting nanocrystals-and demonstrate their use in super-resolution imaging in near-infrared spectral windows of maximal biological transparency. Avalanching nanoparticles (ANPs) can be pumped by continuous-wave lasers, and exhibit all of the defining features of photon avalanching, including clear excitation-power thresholds, exceptionally long rise time at threshold, and a dominant excited-state absorption that is more than 10,000 times larger than ground-state absorption. Beyond the avalanching threshold, ANP emission scales nonlinearly with the 26th power of the pump intensity, owing to induced positive optical feedback in each nanocrystal. This enables the experimental realization of photon-avalanche single-beam super-resolution imaging with sub-70-nanometre spatial resolution, achieved by using only simple scanning confocal microscopy and without any computational analysis. Pairing their steep nonlinearity with existing super-resolution techniques and computational methods, ANPs enable imaging with higher resolution and at excitation intensities about 100 times lower than other probes. The low photon-avalanching threshold and excellent photostability of ANPs also suggest their utility in a diverse array of applications, including sub-wavelength imaging and optical and environmental sensing.

36 MATERIALS SCIENCE↗

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗

Nonlinear post-compression in multi-pass cells in the mid-IR region using bulk materials

We numerically investigate the regime of nonlinear pulse compression at mid-IR wavelengths in a multi-pass cell (MPC) containing a dielectric plate. This post-compression setup allows for ionization-free spectral broadening and self-compression while mitigating self-focusing effects. We find that self-compression occurs for a wide range of MPC and pulse parameters and derive scaling rules that enable its optimization. We also reveal the solitonic dynamics of the pulse propagation in the MPC and its limitations and show that spatiotemporal/spectral couplings can be mitigated for appropriately chosen parameters. In addition, we reveal the formation of spectral features akin to quasi-phase matched degenerate four-wave mixing. Finally, we present two case studies of self-compression at 3-μm and 6-μm wavelengths using pulse parameters compatible with driving high-field physics experiments. The simulations presented in this paper set a framework for future experimental work using few-cycle pulses at mid-IR wavelengths.

Carlson, D. (ORCID:0000000346836463)↗

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↗