SEARCH · Search NASA

Results for “wave parameter”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Identification of stimulated Raman side scattering in near-spherical coronal plasmas on OMEGA EP

Recent spherical-target laser–plasma interaction experiments, performed on the OMEGA EP laser, have been analyzed for stimulated Raman scattering (SRS). This has been motivated by results obtained on the National Ignition Facility (NIF) that have demonstrated the importance of SRS, and in particular SRS side scatter, for directly driven inertial confinement fusion (ICF) conditions [Rosenberg et al. Phys. Rev. Lett. 120, 055001 (2018); Michel et al. Phys. Rev. E 99, 033203 (2019)]. The analysis, based on a generalized ray tracing approach, is described and is shown to explain the observed scattered light spectra: it identifies SRS convective scattering, from portions of each incident beam where the scattered electromagnetic wave is generated in the direction parallel to contours of constant density, as the dominant contribution. This result is novel, as SRS is mostly associated with plasmas of higher electron temperature (T e ≳ 3–4 keV) and longer density scale length (≳600 μ m) than those considered here (T e ~ 1–3 keV, Ln~150 μm) for the relevant overlapped irradiation intensities l ≲ 10 15 W/cm 2 . Furthermore, a combination of the high single-beam intensity on OMEGA EP and the side-scattering geometry, that has been identified, is responsible for this result. It is further suggested that the OMEGA EP platform could provide a good surrogate in which to develop SRS mitigation strategies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

DOE OSTI2023

Cross-beam energy transfer in conditions relevant to direct-drive implosions on OMEGA

In cross-beam energy transfer (CBET), the interference of two laser beams ponderomotively drives an ion-acoustic wave that coherently scatters light from one beam into the other. This redirection of laser beam energy can severely inhibit the performance of direct-drive inertial confinement fusion (ICF) implosions. To assess the role of nonlinear and kinetic processes in direct-drive-relevant CBET, the energy transfer between two laser beams in the plasma conditions of an ICF implosion at the OMEGA laser facility was modeled using particle-in-cell simulations. For typical laser beam intensities, the simulations are in excellent agreement with linear kinetic theory, indicating that nonlinear processes do not play a role in direct-drive implosions. At higher intensities, CBET can be modified by pump depletion, backward stimulated Raman scattering, or ion trapping, depending on the plasma density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Reconfigurable plasma-dielectric hybrid photonic crystal as a platform for electromagnetic wave manipulation and computing

The introduction of 3D printing has enabled fabrication of photonic crystal devices with complex crystal structures that would be challenging to construct using other fabrication methods. In this paper, we construct and characterize a photonic crystal consisting of two 4 x 8 x 8 cubic lattices composed of spherical silicon nitride elements straddling a layer of 8 plasma discharge tubes, creating a 3D hybrid plasma photonic crystal device. Integrating under-dense gaseous plasma elements provides a unique coupling dynamic between the dielectric spheres and the cylindrical plasma discharges, creating a monolithic hybrid photonic crystal with solid state and reconfigurable elements. The device has resonant modes that have attenuation peaks that are either switchable, tunable in amplitude, or tunable in frequency with variations in plasma density. The response of these bands seen with varying plasma density is confirmed through simulations when effects due to the heating of the photonic crystal from the gaseous plasma elements are accounted for in the experiments. We discuss how this reconfigurable device may be used and expanded upon for applications in photonic artificial neural networks and optical computing systems.

3D printing↗

NASA NTRS1991

Large scale motions of Neptune's bow shock: Evidence for control of the shock position by the rotation phase of Neptune's magnetic field

The Voyager 2 spacecraft observed high levels of Langmuir waves before the inbound crossing of Neptune's bow shock, thereby signifying magnetic connection of the bow shock. The Langmuir waves occurred in multiple bursts throughout two distinct periods separated by an 85 minute absence of wave activity. The times of onsets, peaks, and disappearances of the waves were used together with the magnetic field directions and spacecraft position, to perform a 'remote-sensing' analysis of the shape and location of Neptune's bow shock prior to the inbound bow shock crossing. The bow shock is assumed to have a parabolidal shape with a nose location and flaring parameter determined independently for each wave event. The remote-sensing analysis give a shock position consistent with the time of the inbound shock crossing. The flaring parameter of the shock remains approximately constant throughout each period of wave activity but differs by a factor of 10 between the two periods. The absence of waves between two periods of wave activity coincides with a large rotation of the magnetic field and a large increase in the solar wind ram pressure' both these effects lead to magnetic disconnection of the spacecraft from shock. The planetwards motion of the shock's nose from 38.5 R(sub N) to 34.5 R(sub N) during the second time period occurred while the solar wind ram pressure remained constant to within 15 percent. This second period of planetwards motion of the shock is therefore strong evidence for Neptune's bow shock moving in response to the rotation of Neptune's oblique, tilted magnetic dipole. Normalizing the ram pressure, the remotely-sensed shock moves sunwards during the first wave period and planetwards in the second wave period. The maximum standoff distance occurs while the dipole axis is close to being perpendicular to the Sun-Neptune direction. The remote-sensing analysis provides strong evidence that the location of Neptune's bow shock is controlled by Neptune's rotation phase.

Cairns, Iver H.↗

DOE OSTI2022

Tin removal by an annular surface wave plasma antenna in an extreme ultraviolet lithography source

Tin contamination of the collector mirror surface remains one of the crucial issues of EUV (Extreme Ultraviolet) sources, directly impacting the availability of the tool. Hydrogen plasma-based tin removal processes employ hydrogen radicals and ions to interact with tin deposits to form gaseous tin hydride (SnH 4 ), which can be removed through pumping. An annular surface wave plasma (SWP) source developed at the University of Illinois—Urbana Champaign is integrated into the cone and perimeter of the collection mirror for in situ tin removal. The SWP is characterized by high ion and radical densities, low electron temperature, and local generation where etching is needed. This method has the potential to significantly reduce downtime and increase mirror lifetime. Radical probe measurements show hydrogen radical densities in the order of 10 19 m -3 , while Langmuir probe measurements show electron temperatures of up to 6 eV and plasma densities on the order of 10 17–18 m -3 . The generated ions are essential to the tin cleaning and have sufficiently low energy to cause no damage to the collector capping layer. Tin etch rates of up to 270 nm/min were observed in a variety of experimental conditions, including various powers, pressures, flowrates, and temperatures. In conclusion, the high etch rates demonstrated in this study exceed the expected contamination rate of the EUV source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2006

Modeling the Observed QBO and Inter-Annual Variations of the Diurnal Tide in the Mesosphere

In the current version of the Numerical Spectral Model (NSM), the Quasi-biennial Oscillation (QBO) is generated primarily by small-scale gravity waves (GW) from Hines' Doppler Spread Parameterization (DSP). The model does not have topography, and the planetary waves are solely generated by instabilities. We discuss a 3D modeling study that describes the QBO extending from the stratosphere into the upper mesosphere, where the oscillation produces significant inter-annual variations in the diurnal tide. The numerical results are compared with temperature measurements from the SABER (TIMED) and MLS (UARS) instruments obtained by Huang et al. (2006). With a GW source that peaks at the Equator and is taken to be isotropic and independent of season, the NSM generates a QBO with variable periods around 26 months and zonal wind amplitudes of almost 25 m/s at 30 km. As reported earlier, the NSM reproduces the observed equinoctial maxima in the diurnal tide at altitudes around 95 km. The modeled QBO modulates the tide such that the seasonal amplitude maxima can vary from one year to another by as much as 30%. To shed light on the underlying mechanisms, the relative importance of the advection terms are discussed, and they are shown to be important in the stratosphere. At altitudes above 80 km, however, the QBO-related inter-annual variations of the tide are generated primarily by GW momentum deposition. In qualitative agreement with the SABER measurements, the model generates distinct zonal-mean QBO temperature variations in the stratosphere and mesosphere. In the stratosphere, the computed amplitudes are not much smaller than those observed, and the rate of downward propagation at the Equator is reproduced. The modeled temperature amplitudes in the mesosphere, however, are much smaller than those observed. The observed and computed temperature variations of the QBO peak at the Equator but extend with phase reversals to high latitudes, in contrast to the zonal winds that are confined to equatorial latitudes. Hemispherical asymmetries also appear in both the model results and the observations. The temperature amplitudes outside the equatorial region however tend to occur at lower latitudes in the model results. While there is qualitative agreement between the TIMED measurements and the model prediction, there are some areas of significant disagreement that require us to reexamine the present version of the NSM. The numerical results critically depend on the chosen parameters that determine the wave forcing, and there are a number of avenues to improve the performance of the model that had not been tuned to fit the observations. The GW spectrum and its latitude dependence in the troposphere are not well known, and numerical experiments are discussed that describe the related model response. While it appears that eastward propagating Kelvin waves and westward propagating Rossby gravity waves are not the primary source to generate the QBO, the GW forcing can seed the oscillation and act as a catalyst to enhance effectiveness of these planetary waves.

Mayr, Hans G.↗

NASA NTRS2015

Parameter Estimation for Compact Binaries with Ground-Based Gravitational-Wave Observations Using the LALInference

The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star (BNS), a neutron star - black hole binary (NSBH) and a binary black hole (BBH), where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with non-spinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence (CBC) parameter space.

gravitational wave↗

DOE OSTI2022

Multi-harmonic Rutherford island theory

Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1980

Unsteady aerodynamics of conventional and supercritical airfoils

The unsteady aerodynamics of a conventional and a supercritical airfoil are compared by examining measured chordwise unsteady pressure time-histories from four selected flow conditions. Although an oscillating supercritical airfoil excites more harmonics, the strength of the airfoil's shock wave is the more important parameter governing the complexity of the unsteady flow. Whether they are conventional or supercritical, airfoils that support weak shock waves induce unsteady loads that are qualitatively predictable with classical theories; flows with strong shock waves are sensitive to details of the shock-wave and boundary-layer interaction and cannot be adequately predicted.

Davis, S. S.↗

NASA NTRS1986

Effects of cone surface waviness and freestream noise on transition in supersonic flow

A comparison of transition on wavy-wall and smooth-wall cones in a Mach 3.5 wind tunnel is made under conditions of either low freestream noise (quiet flow) or high freestream noise (noisy flow). The noisy flow compares to that found in conventional wind tunnels while the quiet flow gives transitional Reynolds numbers on smooth sharp cones comparable to those found in flight. The waves were found to have a much smaller effect on transition than similar sized trip wires. A satisfatory correlating parameter for the effect of waves on transition was simply the wave height-to-length ratio. A given value of this ratio was found to cause the same percentage change in transition location in quiet and noisy flows.

Morrisette, E. L.↗

DOE OSTI2025

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Adaptive Variational Quantum Imaginary Time Evolution Approach for Ground State Preparation

Abstract An adaptive variational quantum imaginary time evolution (AVQITE) approach is introduced that yields efficient representations of ground states for interacting Hamiltonians on near‐term quantum computers. It is based on McLachlan's variational principle applied to imaginary time evolution of variational wave functions. The variational parameters evolve deterministically according to equations of motions that minimize the difference to the exact imaginary time evolution, which is quantified by the McLachlan distance. Rather than working with a fixed variational ansatz, where the McLachlan distance is constrained by the quality of the ansatz, the AVQITE method iteratively expands the ansatz along the dynamical path to keep the McLachlan distance below a chosen threshold. This ensures the state is able to follow the quantum imaginary time evolution path in the system Hilbert space rather than in a restricted variational manifold set by a predefined fixed ansatz. AVQITE is used to prepare ground states of H 4 , H 2 O, and BeH 2 molecules, where it yields compact variational ansätze and ground state energies within chemical accuracy. Polynomial scaling of circuit depth with system size is shown through a set of AVQITE calculations of quantum spin models. Finally, quantum Lanczos calculations are demonstrated alongside AVQITE without additional quantum resource costs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Generalized slab universal instability and its appearance in pair plasma

In this work, a generalized linear dispersion relation of electromagnetic slab universal modes is derived, taking into account arbitrary ion charge state, electron finite Larmor radius (FLR) effects, and Debye shielding effects. As a consequence, it provides more accurate predictions and is applicable to a wider range of plasmas. We find that electron FLR effects have a weakly stabilizing effect on the slab universal instability in electron–ion plasma, while Debye shielding strongly stabilizes this instability when λ D approaches ρ i (λ D is the Debye length and ρi is the ion gyroradius). In particular, we examine the stability criterion for this instability in electron–positron pair plasmas and find that the instability persists in this simplest plasma system as long as the pair plasma number density exceeds the critical value n c = B 2 /(8πm e c 2 ).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Laser propagation in a subcritical foam: Subgrid model

Here, we present a subgrid model for laser propagation in a subcritical foam. Our model describes the expansion of laser-irradiated foam elements that are below the resolution of the simulation grid and predicts the plasma conditions that result from burning down the foam. Our model can be included as a module within a larger multiphysics code, and we have implemented it within the code pF3D, which is used for simulating a laser-plasma interaction. The model predicts a reduced propagation velocity for a laser through a subcritical foam compared to simulating that foam as a homogeneous gas. This is attributed to the laser energy that goes into burning down the foam microstructure. We compare our model against experimental data by simulating a 2 mg/cc SiO 2 foam shot performed at the Janus laser facility at the Lawrence Livermore National Laboratory. pF3D simulations with the foam model predict hot ion temperatures. This leads to a reduction in the level of stimulated Brillouin scattering (SBS), bringing the simulated level of SBS into agreement with the data. Intensity fluctuations at the foam front due to laser speckles and refraction result in ion temperature fluctuations when the foam burns down. These drive long-lived electron density fluctuations on scales that are large compared to the pore size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Laser-driven plasma sources of intense, ultrafast, and coherent radiation

High-power lasers can deliver extreme light intensities, but avoiding damage in optical components requires large beam sizes, hindering further advances. The use of plasma as a medium for generating and manipulating light avoids the damage thresholds of solid materials and can support extraordinarily bright radiation. In this work, we discuss here how parametric plasma amplification and relativistic high-order harmonic generation offer paths to the development of light sources with peak powers beyond the capabilities of solid-state optics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

High frequency impedance characteristics of a tunable microplasma device

Computational studies on high frequency impedance characteristics of a microplasma device are reported. While microplasma is ignited using a primary excitation signal, frequency response of plasma impedance is determined by a secondary high frequency probe signal with significantly lower voltage amplitude such that it does not influence the plasma parameters. Here, the computational model utilizing the drift–diffusion approximation is first validated by comparing with experimental data for microplasmas ignited at pressures ranging from 1 to 5 Torr. In spite of quantitative discrepancies, good overall agreement is obtained between the measured frequency response of impedance of the discharge. Comparisons are also presented for various plasma parameters including mean electron number density, sheath thickness, mean electron temperature, and collision frequency that were inferred from the impedance measurements. The computational model is then used to perform simulations of near-atmospheric pressure microplasmas with the probe signal frequency ranging from 3 to 20 GHz. The simulations demonstrate the presence of a resonance frequency at which the impedance vanishes. More importantly, it is shown that this resonant frequency can be tuned effectively by suitably modifying the operating parameters (gap size, pressure, and excitation voltage). The simulated impedance characteristics are used to determine the effective plasma inductance and capacitance using a non-linear fitting approach, thereby showing the dependence of these electrical parameters on the plasma operating conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Statistical analysis of non-Maxwellian electron distribution functions measured with angularly resolved Thomson scattering

Angularly resolved Thomson scattering is a novel extension of Thomson scattering, enabling the measurement of the electron velocity distribution function over many orders of magnitude. Here, details of the theoretical basis of the technique and the instrument designed for this measurement are described. Angularly resolved Thomson-scattering data from several experiments are shown with descriptions of the corresponding distribution functions. A reduced model describing the distribution function is given and used to perform a Monte Carlo analysis of the uncertainty in the measurements. The electron density and temperature were determined to a precision of 12% and 21%, respectively, on average, while all other parameters defining the distribution function were generally determined to better than 20%. It was found that these uncertainties were primarily due to limited signal to noise and instrumental effects. Furthermore, measurements with this level of precision were sufficient to distinguish between Maxwellian and non-Maxwellian distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗