Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid 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 613 records · Page 34

Plasma physics in strong-field regimes: Theories and simulations

In strong electromagnetic fields, unique plasma phenomena and applications emerge whose description requires recently developed theories and simulations [Y. Shi, Ph.D. thesis, Princeton University (2018)]. In the classical regime, to quantify effects of strong magnetic fields on three-wave interactions, a convenient formula is derived by solving the fluid model to the second order in general geometry. As an application, magnetic resonances are exploited to mediate laser pulse compression, whereby higher intensity pulses can be produced in wider frequency ranges, as confirmed by particle-in-cell simulations. In even stronger fields, relativistic-quantum effects become important, and a plasma model based on scalar quantum electrodynamics (QED) is developed which unveils observable corrections to Faraday rotation and cyclotron absorption in strongly magnetized plasmas. Beyond the perturbative regime, lattice QED is extended as a numerical tool for plasma physics, using which the transition from wakefield acceleration to electron-positron pair production is captured when laser intensity exceeds the Schwinger threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations of the shock-driven Kelvin–Helmholtz instability in inclined gas curtains

In this paper, we present simulation results for the two-dimensional, shock-driven Kelvin–Helmholtz instability. Simulations are performed with a Mach 2.0 shock propagating through a finite-thickness curtain of gas inclined at an angle α0=30° with respect to the shock plane. After the passage of the shock, the gas curtain is accelerated along its axis. A perturbation develops due to shock reflection near the lower wall, and a Kelvin–Helmholtz instability forms near the vertical center of the curtain. This is the first known numerical reproduction of these phenomena that have previously been observed in experiments with an inclined cylindrical gas column. The effects of varying Mach number and column width were explored in detail to complement experimental data. Additionally, the dependence of the Kelvin–Helmholtz wavelength on Mach number closely matches the relationship observed in experiments. This supports the notion that the observed instability is effectively two-dimensional and inviscid (like classical Kelvin–Helmholtz). The growth rate of the perturbations in the gas curtain was also found to be similar for different Mach numbers. The perturbation at the curtain foot, previously unreported in experiments, was found to have a similar relationship to Mach number as the Kelvin–Helmholtz instability. Both perturbation wavelengths are found to be proportional to layer width. Simulations were performed with the fast interfaces and transport in the atmosphere, an exascale ready, graphics processing unit-accelerated compressible flow solver developed at the University of New Mexico.

42 ENGINEERING↗

Shock-driven dispersal of a corrugated finite-thickness particle layer

A research area emerging in the multiphase flow community is the study of shock-driven multiphase instability (SDMI), a gas–particle analog of the traditional fluid-fluid Richtmyer–Meshkov instability (RMI). In this work, we study the interaction of planar air shocks with corrugated glass particle curtains through the use of numerical simulations with an Eulerian–Lagrangian approach. One objective of this study is to compare the simulated particle curtains to a comparable set of shock tube experiments performed to analyze traditional RMI of a gas curtain. The simulations are set to match the experimental shock Mach numbers and perturbation wavelengths (3.6 and 7.2 mm) while also matching the Atwood number of the experiments to the multiphase Atwood number of the simulations. Varying particle diameters are tested in the simulations to explore the impact of particle diameter on the evolution of the particle curtain. This simulation setup allows for a one-to-one comparison between RMI and SDMI under comparable conditions while also allowing for a separate study into the validity of the use of the multiphase Atwood number to compare the single-phase and multiphase instabilities. In particular, we show that the comparison depends on the diameter of the particles (thus, dependent on the Stokes number of the flow). A second objective of this study is to analyze the effect of the initial particle volume fraction on the evolution of the curtain and the behavior of the instability. This is done through analyzing the effect of the multiphase terms of the vorticity evolution equation on the vorticity deposition in SDMI. Also discussed is the effect of the particle diameter on the multiphase generation terms as well as in the baroclinic vorticity generation term in SDMI as the shock passes over the curtain.

42 ENGINEERING↗

Mitigation of deceleration-phase Rayleigh–Taylor instability growth in inertial confinement fusion implosions

Rayleigh–Taylor growth during shell deceleration is one of the main limiting factors for target performance in inertial confinement fusion implosions. Here, using analytical scaling laws and hydrodynamic simulations, we show that such amplification can be mitigated by reducing the initial mass density in the central target region. The perturbation growth reduction is caused by a smaller hot-spot convergence ratio during deceleration, increased density scale length, and enhanced ablation stabilization. The required central density reduction can be achieved using the dynamic shell formation concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exact Coherent Structures and Phase Space Geometry of Preturbulent 2D Active Nematic Channel Flow

Confined active nematics exhibit rich dynamical behavior, including spontaneous flows, periodic defect dynamics, and chaotic “active turbulence.” Here, we study these phenomena using the framework of exact coherent structures, which has been successful in characterizing the routes to high Reynolds number turbulence of passive fluids. Exact coherent structures are stationary, periodic, quasiperiodic, or traveling wave solutions of the hydrodynamic equations that, together with their invariant manifolds, serve as an organizing template of the dynamics. We compute the dominant exact coherent structures and connecting orbits in a preturbulent active nematic channel flow, which enables a fully nonlinear but highly reduced-order description in terms of a directed graph. Using this reduced representation, we compute instantaneous perturbations that switch the system between disparate spatiotemporal states occupying distant regions of the infinite-dimensional phase space. Overall, our results lay the groundwork for a systematic means of understanding and controlling active nematic flows in the moderate- to high-activity regime.

36 MATERIALS SCIENCE↗

Viscosity Measurements in Extreme Conditions

Radiation-hydrodynamics simulations are critical to the NNSA complex for understanding high energy density physics and for addressing problems in national security science. However, current rad-hydro codes typically do not account for viscosity, which generally leads to large uncertainties in problems involving mixing of materials under shock loading. This report details progress made towards experimental measurements of fluid viscosity at high pressures and temperatures achieved by shock wave compression techniques, as well as complementary material model development and hydrodynamic simulations.

36 MATERIALS SCIENCE↗

Electron-ion relaxation in a dense plasma

The microscopic physics of the thermonuclear runaway in highly degenerate carbon-oxygen cores is investigated to determine if and how a detonation wave is generated. An expression for the electron-ion relaxation time is derived under the assumption of large degeneracy and extreme relativity of the electrons in a two-temperature plasma. Since the nuclear burning time proves to be several orders of magnitude shorter than the relaxation time, it is concluded that in studying the structure of the detonation wave the electrons and ions must be treated as separate fluids.

Littleton, J. E.↗

Vortical dissipation in two-dimensional shear flows

An exact expression is derived for the viscous dissipation function of a real homogeneous and isotropic fluid, which has terms associated with the square of vorticity, wave radiation, and dilatation. The implications of the principle of maximal dissipation rate, are explored by means of this equation for a parallel channel flow and a cylindrical vortex flow. The consequences of a condition of maximum dissipation rate on the growth of disturbances in an unsteady, laminar shear layer are apparently consistent with predictions and observations of maximum growth rate of vortical disturbances. Finally, estimates of the magnitudes of several dissipative components of an unsteady vortex flow are obtained from measurements of a periodic wall jet.

Horne, W. Clifton↗

Density fluctuation spectra in magnetohydrodynamic turbulence

It is shown that within the framework of nearly incompressible magnetohydrodynamics, but not within that of neutral-fluid hydrodynamics, a k exp -5/3 inertial-range wave number density fluctuation spectrum is to be expected at the same times that k exp -5/3 kinetic energy and magnetic energy cascade spectra are present. A previous discrepancy between theory and observation in the local interstellar medium and solar wind is thereby resolved.

Montgomery, D.↗

The onset of Alfvenic turbulence

An investigation is conducted on how low-frequency MHD oscillations in a warm plasma may undergo a transition from a coherent state to one of turbulence. A driven/dissipative derivative nonlinear Schroedinger equation is derived from the fluid equations. The time evolution of an arbitrary spectrum of waves is analyzed in the case where one k-mode is unstable, with the rest damped. It is found that the transition from order to chaos in the driven/dissipative system is correlated with the existence or absence of 'breathing' solitons in the associated conservative system.

Ghosh, S.↗

A general analysis of propagation along multiple-layer superconducting stripline and microstrip transmission lines

A rigorous spectral-domain formulation for a superconducting stripline or microstrip transmission line with a multiple-layer dielectric substrate is presented. The formulation models the strip conductor as a surface current with an equivalent surface impedance, where the surface impedance is approximated in closed form when the strip is either much thinner or much thicker than a penetration depth. In either case the surface impedance is related to the complex conductivity of the material, which is calculated from a two-fluid model. Results are presented to show the slow-wave propagation and attenuation along both microstrip and stripline packages in a realistic multiple-layer configuration, which accounts for the field penetration into the superconducting ground planes.

Nghiem, David↗

A general analysis of propagation along multiple-layer superconducting stripline and microstrip transmission lines

A rigorous spectral-domain formulation for a superconducting stripline or microstrip transmission line with a multiple-layer dielectric substrate is presented. The formulation models the strip conductor as a surface current with an equivalent surface impedance, where the surface impedance is approximated in closed form when the strip is either much thinner or much thicker than a penetration depth. In either case the surface impedance is related to the complex conductivity of the material, which is calculated from a two-fluid model. Results are presented to show the slow-wave propagation and attenuation along both microstrip and stripline packages in a realistic multiple-layer configuration, which accounts for the field penetration into the superconducting ground planes.

Nghiem, David↗

A second-order theory for transverse ion heating and momentum coupling due to electrostatic ion cyclotron waves

A second-order theory for electrostatic instabilities driven by counterstreaming ion beams is developed which describes momentum coupling and heating of the plasma via wave-particle interactions. Exchange rates between the waves and particles are derived, which are suitable for the fluid equations simulating microscopic effects on macroscopic scales. Using a fully kinetic simulation, the electrostatic ion cyclotron instability due to counterstreaming H(+) beams has been simulated. A power spectrum from the kinetic simulation is used to evaluate second-order exchange rates. The calculated heating and momentum loss from second-order theory is compared to the numerical simulation.

Miller, Ronald H.↗

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗