Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hydrodynamic Instability”

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 361 records · Page 20

Gravitational instabilities in protostellar disks

The nonaxisymmetric stability of self-gravitating, geometrically thick accretion disks has been studied for protostellar systems having a wide range of disk-to-central object mass ratios. Global eigenmodes with four distinctly different characters were identified using numerical, nonlinear hydrodynamic techniques. The mode that appears most likely to arise in normal star formation settings, however, resembles the 'eccentric instability' that was identified earlier in thin, nearly Keplerian disks: It presents an open, one-armed spiral pattern that sweeps continuously in a trailing direction through more than 2-pi radians, smoothly connecting the inner and outer edges of the disk, and requires cooperative motion of the point mass for effective amplification. This particular instability promotes the development of a single, self-gravitating clump of material in orbit about the point mass, so its routine appearance in our simulations supports the conjecture that the eccentric instability provides a primary route to the formation of short-period binaries in protostellar systems.

Tohline, J. E.↗

Nonspherical supernova remnants. IV - Sequential explosions in OB associations

Multisupernova remnants, driven by sequential supernova explosions in OB associations, are modelled by means of two-dimensional hydrodynamical calculations. It is shown that due to the Rayleigh-Taylor instability the remnants quickly evolve into highly irregular structures. A critical evaluation of the multisupernova model as an explanation for supershells is given.

Tenorio-Tagle, G.↗

Gamma-Ray Bursts and Relativistic MHD around Coalescing Neutron Stars

In this proposal we performed the first relativistic hydrodynamic simulations of the binary neutron stars near coalescence. Calculations were performed in three spatial dimensions and were supplemented with magnetohydrodynamic simulations of the magnetic reconnection as the neutron star fluid responds to the relativistic forces of the binary. The ultimate goal of the proposed work was to develop a theoretical model for the temporal and spectral evolution of cosmological gamma-ray bursts produced during the magnetic reconnection. This is the first time such hydrodynamic calculations were performed and some surprising results were obtained. The neutrons stars were observed to compress and heat in the simulations. They could even collapse to black holes many orbits before coalescence. The question then arose as to whether this strange phenomenon could power a gamma-ray burst. The bulk of computing effort in this proposal was therefore devoted to an analysis of the implications of this possible new gamma-ray burst paradigm. The initial focus of the work performed under this proposal was therefore to utilize (3+1) dimensional and spherical numerical general relativistic hydrodynamics to study the origin, evolution, and parametric sensitivity of the collapse instability. We derived modified conditions of hydrostatic equilibrium for stars in the curved space of quasi-static orbits. We followed the magnetic field evolution for approximately one orbit (10 msec). The magnetic field energy was observed to exponentiate with an e-folding time of about 1 msec. Thus, the field should build up very quickly to a magnitude such that reconnection and back reaction of the fluid inhibits further growth.

Mathews, Grant J.↗

Optimal Disturbances in Boundary Layers Subject to Streamwise Pressure Gradient

Laminar-turbulent transition in shear flows is still an enigma in the area of fluid mechanics. The conventional explanation of the phenomenon is based on the instability of the shear flow with respect to infinitesimal disturbances. The conventional hydrodynamic stability theory deals with the analysis of normal modes that might be unstable. The latter circumstance is accompanied by an exponential growth of the disturbances that might lead to laminar-turbulent transition. Nevertheless, in many cases, the transition scenario bypasses the exponential growth stage associated with the normal modes. This type of transition is called bypass transition. An understanding of the phenomenon has eluded us to this day. One possibility is that bypass transition is associated with so-called algebraic (non-modal) growth of disturbances in shear flows. In the present work, an analysis of the optimal disturbances/streamwise vortices associated with the transient growth mechanism is performed for boundary layers in the presence of a streamwise pressure gradient. The theory will provide the optimal spacing of the control elements in the spanwise direction and their placement in the streamwise direction.

Tumin, Anatoli↗

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational study of instability and fill tube mitigation strategies for double shell implosions

Double shell capsules are an attractive alternative scheme for achieving robust alpha-heating at the National Ignition Facility due to their low convergence ratio and low predicted temperature for achieving volume ignition. Nevertheless, simulations suggest that double shell targets are more susceptible to the fill tube, used to fill the inner shell with liquid DT, than typical single-shell ignition capsule designs, due to the higher density gradient between the shell and the fill tube hole, a lower outer shell velocity, which prevents the implosion from catching up to the initial fill tube jet, and the absence of a rebounding shock through the foam to slow this jet. Double shells are also highly susceptible to the Rayleigh-Taylor instability at both interfaces with the high density inner shell. Combined, these effects are predicted by radiation-hydrodynamics simulations to reduce fuel confinement and temperature, resulting in reduced performance by a factor of ≈20–45, depending on design details, compared to idealized one-dimensional (1D) simulations. We discuss herein a mitigation strategy for both the interfacial instabilities and the fill tube that is predicted by simulations to decrease the yield degradation to a factor of ≈4. The mitigation strategy involves a modification of the capsule geometry as well as the use of a multishock pulse shape. The multishock pulse is required for the fill tube mitigation strategy and has the added benefit of stabilizing perturbations at the foam/pusher interface without decreasing 1D yield. In order to experimentally verify these predictions, we discuss the potential use of a hydrogrowth radiography platform that could be applied to test the proposed mitigation strategies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗

The fragmentation of proto-globular clusters. I - Thermal instabilities

The metal abundances among the stars within a typical globular cluster are remarkably homogeneous. This indicates that star formation in these systems was a globally coordinated event which occurred over a time span less than or comparable to the collapse time scale of the cluster. This issue is addressed by assuming that the fragmentation of a proto-globular cluster cloud proceeded in two steps. In the first step, thermal instability led to the rapid growth of initial fluctuations. This led to a large contrast in the dynamical time scales between the perturbations and the parent cloud, and the perturbations then underwent gravitational instabilities on short time scales. This process is modeled using one-dimensional hydrodynamic simulations of clouds both with and without external heat sources and self-gravity. The models include the effects of a non-equilibrium H2 abundance. The results indicate that fragmentation can occur on time scales significantly less than the dynamical time scale of the parent cloud.

Murray, Stephen D.↗

Detrimental effects and mitigation of the joint feature in double shell implosion simulations

Double shell capsules provide an attractive option in inertial confinement fusion experiments due to their potential for achieving a low-convergence, robust burn. However, these designs suffer from symmetry degradation and accompanying reduced fuel confinement due to the currently necessary joint between the two hemispheres of the outer shell. Here, the gap widens as a result of the excess ablation pressure produced by x rays that penetrate the joint during the drive phase, and this perturbation grows and imprints onto the inner shell during the collision. xRAGE Eulerian radiation-hydrodynamic simulations predict significant reductions in deuterium–tritium fusion yields compared to joint-less simulations when the depth of the outer joint is increased, whereas the performance is less sensitive to the depth of the inner gap. Here we examine the technique of plating the insides of the outer gap with a high-Z material to mitigate the impact of this feature. Gold-plating in quantities comparable to or exceeding the “missing” outer shell mass shows promise toward restoring both implosion symmetry and yield closer to the joint-less levels, and synthetic diagnostics suggest that high-energy x-ray radiographs can capture this shape retention of the inner and outer shells in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Automated Detection of Instability-Inducing Channel Geometry Transitions in Saint-Venant Simulation of Large-Scale River Networks

A new sweep-search algorithm (SSA) is developed and tested to identify the channel geometry transitions responsible for numerical convergence failure in a Saint-Venant equation (SVE) simulation of a large-scale open-channel network. Numerical instabilities are known to occur at “sharp” transitions in discrete geometry, but the identification of problem locations has been a matter of modeler’s art and a roadblock to implementing large-scale SVE simulations. The new method implements techniques from graph theory applied to a steady-state 1D shallow-water equation solver to recursively examine the numerical stability of each flowpath through the channel network. The SSA is validated with a short river reach and tested by the simulation of ten complete river systems of the Texas–Gulf Coast region by using the extreme hydrological conditions recorded during hurricane Harvey. The SSA successfully identified the problematic channel sections in all tested river systems. Subsequent modification of the problem sections allowed stable solution by an unsteady SVE numerical solver. The new SSA approach permits automated and consistent identification of problem channel geometry in large open-channel network data sets, which is necessary to effectively apply the fully dynamic Saint-Venant equations to large-scale river networks or for city-wide stormwater networks.

54 ENVIRONMENTAL SCIENCES↗

Constraining computational modeling of indirect drive double shell capsule implosions using experiments

Double shell capsule implosions are an alternative approach to achieving alpha heating on the National Ignition Facility. Current machining techniques construct the outer shell as two hemispheres that are glued together, and the deuterium and tritium (DT) liquid inside the inner shell will be injected by a fill tube. These features introduce asymmetries and jetting that may disrupt the confinement of the DT fuel if not carefully controlled. Simulations indicate that in order to achieve high yields in the laboratory, these features as well as susceptibility to the Rayleigh–Taylor instability (RTI) must be mitigated. Furthermore, due to uncertainties in computational models and the expense of using the best physics models at adequate resolution in three dimensions, our computational modeling must be constrained by experiments. We report on the results of recent hydrogrowth radiography and dual-axis keyhole experiments with double shell targets that have been used to evaluate our modeling of the outer shell joint as well as the impacts of high-energy x-ray preheat that strongly impacts RTI growth. Our simulations show good agreement with the experimental data and inform several important modeling choices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Asymptotic scaling laws for the stagnation conditions of Z-pinch implosions

Implosions of magnetically driven annular shells (Z pinches) are studied in the laboratory to produce high-energy-density plasmas. Such plasmas have a wide-range of applications including x-ray generation, controlled thermonuclear fusion, and astrophysics studies. In this work, we theoretically investigate the in-flight dynamics of a magnetically driven, imploding cylindrical shell that stagnates onto itself upon collision on axis. The converging flow of the Z-pinch is analyzed by considering the implosion trajectory in the (A, M) parametric plane, where A is the in-flight aspect ratio and M is the implosion Mach number. For an ideal implosion in the absence of instabilities and in the limit of A ≫ 1, we derive asymptotic scaling laws for hydrodynamic quantities evaluated at stagnation (e.g., density, temperature, and pressure) and for performance metrics (e.g., soft x-ray emission, K-shell x-ray emission, and neutron yield) as functions of target-design parameters.

Ruiz, D. E. [Sandia National Laboratories (SNL-NM)↗

A Parametric Study of the SASI Comparing General Relativistic and Nonrelativistic Treatments

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

The Stability of Radiatively Cooling Jets. 2: Nonlinear Evolution

We use two-dimensional time-dependent hydrodynamical simulations to follow the growth of the Kelvin-Helmholtz (K-H) instability in cooling jets into the nonlinear regime. We focus primarily on asymmetric modes that give rise to transverse displacements of the jet beam. A variety of Mach numbers and two different cooling curves are studied. The growth rates of waves in the linear regime measured from the numerical simulations are in excellent agreement with the predictions of the linear stability analysis presented in the first paper in this series. In the nonlinear regime, the simulations show that asymmetric modes of the K-H instability can affect the structure and evolution of cooling jets in a number of ways. We find that jets in which the growth rate of the sinusoidal surface wave has a maximum at a so-called resonant frequency can be dominated by large-amplitude sinusoidal oscillations near this frequency. Eventually, growth of this wave can disrupt the jet. On the other hand, nonlinear body waves tend to produce low-amplitude wiggles in the shape of the jet but can result in strong shocks in the jet beam. In cooling jets, these shocks can produce dense knots and filaments of cooling gas within the jet. Ripples in the surface of the jet beam caused by both surface and body waves generate oblique shock "spurs" driven into the ambient gas. Our simulations show these shock "spurs" can accelerate ambient gas at large distances from the jet beam to low velocities, which represents a new mechanism by which low-velocity bipolar outflows may be driven by high-velocity jets. Rapid entrainment and acceleration of ambient gas may also occur if the jet is disrupted. For parameters typical of protostellar jets, the frequency at which K-H growth is a maximum (or highest frequency to which the entire jet can respond dynamically) will be associated with perturbations with a period of - 200 yr. Higher frequency (shorter period) perturbations excite waves associated with body modes that produce internal shocks and only small-amplitude wiggles within the jet. The fact that most observed systems show no evidence for large-amplitude sinusoidal oscillation leading to disruption is indicative that the perturbation frequencies are generally large, consistent with the suggestion that pro- tostellar jets arise from the inner regions (r less than 1 AU) of accretion disks.

Stone, James M.↗

The Evolution of a Circumplanetary Disc with a Dead Zone

We investigate whether the regular Galilean satellites could have formed in the dead zone of a circumplanetary disc. A deadzone is a region of weak turbulence in which the magnetorotational instability is suppressed, potentially an ideal environmentfor satellite formation. With the grid-based hydrodynamic codeFAR GO3D, we examine the evolution of a circumplanetary discmodel with a dead zone. Material accumulates in the dead zone of the disc leading to a higher total mass and but a similartemperature profile compared to a fully turbulent disc model. The tidal torque increases the rate of mass transport through thedead zone leading to a steady-state disc with a dead zone that does not undergo accretion outbursts. We explore a range of disc,dead zone, and mass inflow parameters and find that the maximum mass of the disc is around 0.001MJ. Since the total solidmass of such a disc is much lower, we find that there is not sufficient material in the disc forin situformation of the Galileansatellites and that external supplement is required.

Cheng Chen↗

Hall interchange instability as a seed for helical magneto-Rayleigh–Taylor instabilities in magnetized liner inertial fusion Z-Pinches scaled from Z-Machine parameters to a next generation pulsed power facility

Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion concept that is studied on the 20-MA, 100-ns rise time Z Pulsed Power Facility at Sandia National Laboratories. Given the relative success of the platform, there is a wide interest in studying the scaled performance of this concept at a next-generation pulsed-power facility that may produce peak currents upward of 60 MA. An important aspect that requires more research is the instability dynamics of the imploding MagLIF liner, specifically how instabilities are initially seeded. It has been shown in magnetized 1-MA thin-foil liner Z-pinch implosion simulations that a Hall interchange instability (HII) effect can provide an independent seeding mechanism for helical magneto-Rayleigh–Taylor instabilities. Here in this paper, we explore this instability at higher peak currents for MagLIF using 2D discontinuous Galerkin PERSEUS simulations, an extended magneto-hydrodynamics code, which includes Hall physics. Our simulations of scaled MagLIF loads show that the growth rate of the HII is invariant to the peak current, suggesting that studies at 20-MA are directly relevant to 60-MA class machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiation and heat transport in divergent shock–bubble interactions

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrodynamic interactions in directional solidification

A binary liquid undergoes unidirectional solidification. The one-dimensional steady state is susceptible to morphological instability that causes the solid/liquid interface to change from a planar state to a cellular pattern. This paper examines the effects on this transition of volume-change convection, buoyancy-driven convection or forced flows. It emphasizes how flows alter stability limits, create scale and pattern changes in morphology, and create, through coupling, new instabilities. Emphasis is placed on the physical mechanisms of the interactions.

Davis, Stephen H.↗