Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Meshkov”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulations of three-layer Richtmyer–Meshkov mixing in a shock tube

Here, the Richtmyer–Meshkov instability causes perturbations to grow after a shock traverses a fluid density interface. This increases the mixing rate between fluid from either side of the interface. We use the Flash Eulerian hydrodynamic code to investigate alterations when a thin third layer of intermediate density is placed along the interface, effectively creating two adjacent unstable interfaces. This is a common occurrence in engineering applications where a thin barrier initially separates two materials. We find that the width of the mixing layer is similar or slightly reduced; however, the total mass of mixed material can actually increase. The mixing layer becomes more compact and efficient. However, the normalized mixed mass decreases, meaning that finger entrainment becomes more important than in the simple two-layer case. The effect of adding the central layer appears to decrease when the Atwood number is decreased. The Flash results are also benchmarked against two-layer experimental data from a shock tube at the University of Arizona.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Analysis of single-mode Richtmyer–Meshkov instability using high-order incompressible vorticity-streamfunction and shock-capturing simulations

In this study, two- and three-dimensional simulation results obtained using a new high-order incompressible, variable-density vorticity–streamfunction (VS) method and data from previous ninth-order weighted essentially nonoscillatory (WENO) shock-capturing simulations are used to investigate the nonlinear dynamics of single-mode Richtmyer–Meshkov instability using a model of a Mach 1.3 air(acetone)/SF6 shock tube experiment. A comparison of the density fields from both simulations with the experimental images demonstrates very good agreement in the large-scale structure with both methods but differences in the small-scale structure. The WENO method captures the small-scale disordered structure observed in the experiment, while the VS method partially captures such structure and yields a strong rotating core. The perturbation amplitude growth from the simulations generally agrees well with the experiment. The simulation bubble and spike amplitudes agree well at early times. At later times, the WENO bubble amplitude is smaller than the VS amplitude and vice versa for the spike amplitude. The predictions of nonlinear single-mode instability growth models are shown to agree with the simulation amplitudes at early-to-intermediate times but underpredict the amplitudes at later times in the nonlinear regime. Visualizations of the mass fraction and enstrophy isosurfaces, velocity and vorticity fields, and baroclinic vorticity production and vortex stretching terms from the three-dimensional simulations indicate that, with the exception of the small-scale structure within the rollups, the VS and WENO results are in good agreement.

42 ENGINEERING↗
DOE OSTI2024

Shock tube experiments on the three-layer Richtmyer–Meshkov instability

Here, a vertical shock tube is used for experiments on the three-layer Richtmyer–Meshkov instability. Two closely spaced membrane-less interfaces are formed by the flow of two different sects of three gases: one with air above CO 2 above SF 6 and the other with helium above air above SF 6 . The lightest of the three gases enters the shock tube at the top of the driven section and flows downward. Conversely, the heaviest gas enters at the bottom of the shock tube and flows upward while the intermediate density gas enters at the middle through porous plates. All three gases are allowed to escape through holes at the layer location, leaving an approximately 30-mm layer of intermediate-density gas suspended between the lightest gas from above and the heaviest gas from below. A single-mode, two-dimensional initial perturbation is then imposed on the lower interface by oscillating the shock tube in the horizontal direction. The flow is visualized by seeding the intermediate gas with particles and illuminating it with a pulsed laser. Image sequences are then captured using high-speed video cameras. Perturbation amplitude measurements are made from the three-layer system and compared with measurements from 2, two-layer systems. It is observed that the presence of the upper, initially flat interface produces a decrease in growth of instability amplitude in the nonlinear phase over an equivalent single-interface configuration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Temporal evolution of scalar modes in Richtmyer–Meshkov instability of inclined interface using high-speed PIV and PLIF measurements at 60 kHz

Abstract The current work presents simultaneous, high-speed measurements at 60,000 fields per second of velocity and mole fraction using particle image velocimetry (PIV) and planar laser induced acetone-fluorescence in a Richtmyer–Meshkov instability of an inclined interface (Atwood number, At = 0.22). Specifically, around 2 ms of temporal evolution of the vortex structures and their associated scalar modes immediately following the interface-reshock interaction is presented. Two initial interface conditions are discussed—(a) a sharp, inclined ‘single mode’ interface and (b) a ‘multi-mode’ interface where small perturbations are imposed on the single mode case. A 2D wavelet decomposition of the scalar flow field shows a highly intermittent distribution of small-scale variance throughout the interface even at late times. These are correlated strongly with the vortex structures and local turbulence intensity, where each small-scale scalar mode is sandwiched between two co-rotating vortex structures. This indicates that the interstitial regions between the vortices are significant hotspots of entrainment, which is then dispersed by the induced, counter-flow velocity fields. The multimode case demonstrates similar organization at large scales, while the scalar field is much more homogeneous at smaller scales. These observations highlight the importance of capturing the early time vortex evolution to accurately estimate any late time intermittency, especially where deposition of intense vorticity on sharp interfaces is present.

Engineering↗
DOE OSTI2026

Passive Freeze-Out of the Richtmyer-Meshkov Instability

The Richtmyer-Meshkov instability (RMI) poses a major challenge in inertial confinement fusion (ICF) due to its role in mixing and performance degradation. We report the first experimental observation of passive freeze-out of RMI in a low-pressure surrogate regime, an instability stagnation effect induced without modifying the driving pressure pulse or the target surface geometry. Using additively manufactured subsurface voids in a sinusoidal target, we convert a single shock into a sequence of weaker shocks that suppress instability growth upstream of the surface by over 70%. High-speed x-ray imaging and hydrodynamic simulations suggest that this suppression arises primarily from temporal shaping, with lesser contributions from spatial curvature and shock weakening. Our results demonstrate a driver-independent pathway for controlling shock-driven hydrodynamic instabilities relevant to ICF and other high energy density systems.

Materials science↗
DOE OSTI2021

Benchmarking near-term quantum devices with the variational quantum eigensolver and the Lipkin-Meshkov-Glick model

The variational quantum eigensolver is a promising algorithm for noisy intermediate scale quantum (NISQ) computation. Verification and validation of NISQ algorithms' performance on NISQ devices is an important task. Here, we consider the exactly diagonalizable Lipkin-Meshkov-Glick (LMG) model as a candidate for benchmarking NISQ computers. We use the Bethe Ansatz to construct eigenstates of the trigonometric LMG model using quantum circuits inspired by the LMG's underlying algebraic structure. We construct circuits with depth $\mathcal{O}$(N) and $\mathcal{O}$(log 2 N) that can prepare any trigonometric LMG eigenstate of N particles. The number of gates required for both circuits is $\mathcal{O}$(N). The energies of the eigenstates can then be measured and compared to the exactly known answers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Low viscosity of solid MgO at high pressures and strain rates measured using the laser-driven Richtmyer-Meshkov instability

Solids are often assumed to behave as viscous fluids under high-strain rates. This behavior has been studied experimentally in metals but largely unexplored in brittle ceramic materials. Here, in this study, we present a technique for measuring the viscosity of MgO using time-resolved velocimetry to track the growth rate of a perturbation caused by the Richtmyer-Meshkov instability at the OMEGA EP laser facility. To interpret the results, we use an in-house Eulerian hydrocode to simulate our experiments and model the plastic deformation of solid MgO as a viscous fluid. Results indicate that MgO has a surprisingly low upper bound to its effective viscosity of ∼10 2 ⁢ Pa s at 175 ± 15⁢ GPa, ∼3500⁢K, and 10 6 −10 7 s −1 strain rate.

Perez, Tyler M. [Johns Hopkins Univ., Baltimore, M↗
DOE OSTI2023

Quantum simulations in effective model spaces: Hamiltonian-learning variational quantum eigensolver using digital quantum computers and application to the Lipkin-Meshkov-Glick model

Quantum simulations offer the potential to predict the structure and dynamics of nuclear many-body systems that are beyond the capabilities of classical computing. Generally, preparing the ground state of strongly-interacting many-body systems relevant to nuclear physics is however inefficient, even using ideal quantum computers. In addition, currently available NISQ-era quantum devices possess modest numbers of qubits, limiting the size of quantum many-body systems that can be simulated. In this context, a reformulation of the quantum many-body problems using truncated model spaces and Hamiltonians is desirable to make them more amenable to near-term quantum computers. The importance of symmetries in low-energy theories, including effective field theories (EFTs), lattice quantum chromodynamics (QCD), and effective model spaces for nuclear systems, in particular their interplay with the reduction of active Hilbert spaces, is well known. Lesser known is the fact that the non-commutivity of some symmetries and truncations of the model space can be profitably combined with variational calculations to rearrange the entanglement into localized structures and enable more efficient simulations. Here, the goal of the present study is to explore and utilize the non-commutivity of symmetries and model-space truncations of quantum many-body systems important to nuclear physics, particularly in combination with variational algorithms for quantum simulations and effective Hamiltonian learning. We introduce an iterative hybrid classical-quantum algorithm, Hamiltonian learning variational quantum eigensolver (HL-VQE), that simultaneously optimizes an effective Hamiltonian, thereby rearranging entanglement into the effective model space, and the associated ground-state wavefunction. Quantum simulations, using classical computers and IBM's superconducting-qubit quantum computers, are performed to demonstrate the HL-VQE algorithm, in the context of the Lipkin-Meshkov-Glick (LMG) model of interacting fermions, where the Hamiltonian transformation corresponds to an orbital rotation. We use a mapping where the number of qubits scales with the $\log$ of the size of the effective model space, rather than the particle number. HL-VQE is found to provide an exponential improvement in LMG-model calculations of the ground-state energy and wavefunction, compared to naive truncations without Hamiltonian learning, throughout a significant fraction of the Hilbert space. In the context of EFT, this corresponds to counterterms scaling exponentially with the cut-off as opposed to power law. Implementations on IBM's QExperience quantum computers and simulators for 1- and 2-qubit effective model spaces are shown to provide accurate and precise results, reproducing classical predictions. For a range of parameters defining the LMG model, the HL-VQE algorithm is found to have better scaling of quantum resources requirements than previously explored algorithms. In particular, the HL-VQE scales efficiently over a large fraction of the model space, in contrast to VQE alone. This work constitutes a step in the development of entanglement-driven quantum algorithms for descriptions of nuclear many-body systems. This, in part, leverages the potential of noisy intermediate-scale quantum (NISQ) devices. The exponential scaling of counterterms observed in this study suggests the possibility of more general applicability to other non-perturbative EFTs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Exact, approximate, and hybrid treatments of viscous Rayleigh-Taylor and Richtmyer-Meshkov instabilities

We consider Rayleigh-Taylor and Richtmyer-Meshkov instabilities at the interface between two fluids, one or both of which may be viscous. We derive exact analytic expressions for the amplitude η(t) in the linear regime when only one of the fluids is viscous. The more general case is solved numerically using Laplace transforms. We compare the exact solutions of the initial-value problem with the approximate solutions of the eigenvalue problem used in a simple expression for η(t) in terms of two growth rates, γ+ and γ-. We then propose a hybrid model as an improvement on the approximate model. The hybrid model is based on the same expression for η(t) in terms of γ± but uses exact eigenvalues for γ+, the larger growth rate, and a relationship between γ- and γ+. We also discuss two concepts: isogrowth wave number pairs and asymptotic decay. The first relies on viscosity in one or both fluids to identify perturbations of two different wavelengths having the same γ+. The second concept, which is more general, can be found in viscous as well as inviscid fluids and requires only a specific initial growth rate $\dot{η}$$critical\atop{0}$ to force η(t) → 0 as t → ∞. We present several examples illustrating these two concepts and comparing exact, approximate, and hybrid treatments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

High-resolution Navier-Stokes simulations of Richtmyer-Meshkov instability with reshock

The interaction of a Mach-1.45 shock wave with a perturbed planar interface between sulfur hexafluoride and air is studied through high-resolution two-dimensional (2D) and three-dimensional (3D) shock-capturing adaptive mesh refinement simulations of multispecies Navier-Stokes equations. The sensitivities of time-dependent statistics on grid resolution for 2D and 3D simulations are evaluated to ensure the accuracy of the results. Here, the numerical results are used to examine the differences between the development of 2D and 3D Richtmyer-Meshkov instability during two different stages: (1) initial growth of hydrodynamic instability from multimode perturbations after first shock and (2) transition to chaotic or turbulent state after reshock. The effects of the Reynolds number on the mixing in 3D simulations are also studied through varying the transport coefficients.

42 ENGINEERING↗
DOE OSTI2022

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗
DOE OSTI2024

Suppression of Richtmyer-Meshkov Instability via Special Pairs of Shocks and Phase Transitions

The classical Richtmyer-Meshkov instability (RMI) is a hydrodynamic instability characterizing the evolution of an interface following shock loading. In contrast to other hydrodynamic instabilities such as Rayleigh-Taylor, it is known for being unconditionally unstable: regardless of the direction of shock passage, any deviations from a flat interface will be amplified. In this article, we show that for negative Atwood numbers, there exist special sequences of shocks which result in a nearly perfectly suppressed instability growth. Here, we demonstrate this principle computationally and experimentally with stepped fliers and phase transition materials. A fascinating immediate corollary is that in specific instances, a phase-transitioning material may self-suppress RMI.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Variational Simulation of the Lipkin-Meshkov-Glick Model on a Neutral Atom Quantum Computer

We simulate the Lipkin-Meshkov-Glick model using the variational-quantum-eigensolver algorithm on a neutral atom quantum computer. We test the ground-state energy of spin systems with up to 15 spins. Two different encoding schemes are used: an individual spin encoding where each spin is represented by one qubit, and an efficient Gray code encoding scheme that only requires a number of qubits that scales with the logarithm of the number of spins. This more efficient encoding, together with zero-noise extrapolation techniques, is shown to improve the fidelity of the simulated energies with respect to exact solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗
DOE OSTI2020

Dependence of Enstrophy Transport and Mixed Mass on Dimensionality and Initial Conditions in the Richtmyer–Meshkov Instability Induced Flows

Here we present a comparative study of the enstrophy budget and mixed mass between two- and three-dimensional flows induced by Richtmyer–Meshkov instability (RMI). Specifically, the individual contributions to the enstrophy budget due to the production from baroclinicity and from vortex stretching (which vanishes in two-dimensional (2D) flow) are delineated. This is enabled by a set of two- and three-dimensional computations at Atwood 0.5 having both narrow- and broad-band perturbations. A further three-dimensional (3D) computation is conducted at Atwood 0.9 using an identical narrowband perturbation to the Atwood 0.5 case to examine the sensitivity to density ratio. The mixed mass is also considered with the goal to obtain insight on how faithfully a simplified calculation performed in two dimensions can capture the mixed mass for an inertial confinement fusion (ICF) or other practical application. It is shown that the late time power law decay of variable density enstrophy is substantially different in two and three dimensions for the narrowband initial perturbation. The baroclinic production term is negligible in three dimensions (aside from the initial shock interaction), as vortex stretching is larger by two orders of magnitude. The lack of vortex stretching considerably reduces the decay rate in both narrowband and broadband perturbations in two dimensions. In terms of mixed mass, the lack of vortex stretching reduces the mixed mass in two dimensions compared to three in all cases. In the broadband cases, the spectral bandwidth in the 2D case is wider; hence, there is a longer time period of sustained linear growth which reduces the normalized mixed mass further.

42 ENGINEERING↗
DOE OSTI2023

Multi-body entanglement and information rearrangement in nuclear many-body systems: a study of the Lipkin–Meshkov–Glick model

Here, we examine how effective-model-space (EMS) calculations of nuclear many-body systems rearrange and converge multi-particle entanglement. The generalized Lipkin-Meshkov-Glick (LMG) model is used to motivate and provide insight for future developments of entanglement-driven descriptions of nuclei. The effective approach is based on a truncation of the Hilbert space together with a variational rotation of the qubits (spins), which constitute the relevant elementary degrees of freedom. The non-commutivity of the rotation and truncation allows for an exponential improvement of the energy convergence throughout much of the model space. Our analysis examines measures of correlations and entanglement, and quantifies their convergence with increasing cut-off. We focus on one- and two-spin entanglement entropies, mutual information, and $n$-tangles for $n=2,4$ to estimate multi-body entanglement. The effective description strongly suppresses entropies and mutual information of the rotated spins, while being able to recover the exact results to a large extent with low cut-offs. Naive truncations of the bare Hamiltonian, on the other hand, artificially underestimate these measures. The $n$-tangles in the present model provide a basis-independent measures of $n$-particle entanglement. While these are more difficult to capture with the EMS description, the improvement in convergence, compared to truncations of the bare Hamiltonian, is significantly more dramatic. We conclude that the low-energy EMS techniques, that successfully provide predictive capabilities for low-lying observables in many-body systems, exhibit analogous efficacy for quantum correlations and multi-body entanglement in the LMG model, motivating future studies in nuclear many-body systems and effective field theories relevant to high-energy physics and nuclear physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Mshock, a thin layer Richtmyer-Meshkov instability experiment

The Mshock campaign is studying the feed through of the Richtmyer-Meshkov instability (RMI) in a thin layer in the high energy density (HED) regime. Experiments utilize a beryllium shock tube filled with low density foam and a high density plastic to study the instability. Initial experiments were performed at the OMEGA laser facility and recently at the National Ignition facility. While the basic design of the experiment is the same, critical components were changed to optimize the data quantity and clarity for the more energetic NIF.

42 ENGINEERING↗