0
Research Papers: Fundamental Issues and Canonical Flows

Evaporation Effects in Shock-Driven Multiphase Instabilities

[+] Author and Article Information
Wolfgang J. Black

Department of Mechanical Engineering,
University of Missouri,
Columbia, MO 65202
e-mail: wjbvg5@missouri.com

Nicholas A. Denissen

Department of Computational Physics,
Los Alamos National Laboratory,
Los Alamos, NM 87545
e-mail: Denissen@lanl.gov

Jacob A. McFarland

Department of Mechanical Engineering,
University of Missouri,
Columbia, MO 65202
e-mail: mcfarlandja@missouri.edu

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received August 19, 2016; final manuscript received February 13, 2017; published online April 26, 2017. Assoc. Editor: Daniel Livescu.The United States Government retains, and by accepting the article for publication, the publisher acknowledges that the United States Government retains, a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for United States Government purposes.

J. Fluids Eng 139(7), 071204 (Apr 26, 2017) (15 pages) Paper No: FE-16-1528; doi: 10.1115/1.4036162 History: Received August 19, 2016; Revised February 13, 2017

This paper considers the effects of multiphase parameters on a shock-driven particle-laden hydrodynamic instability using simulations performed with the hydrocode FLAG, developed at Los Alamos National Laboratory. The classic sinusoidal interface common in instability literature is created using water particles seeded in a nitrogen–water vapor mixture. The simulations model a shock tube environment as the computational domain, to guide future experimentation. Multiphase physics in FLAG include momentum and energy coupling, with this paper discussing the addition of mass coupling through evaporation. The multiphase effects are compared to a dusty gas approximation, which ignores multiphase components, as well as to a multiphase case which ignores evaporation. Evaporation is then further explored by artificially changing parameters which effect the rate of evaporation as well as the amount of total evaporation. Among all these experiments, the driving force of the hydrodynamic instability is a shock wave with a Mach number of 1.5 and a system Atwood number of 0.11 across the interface. The analysis is continued into late time for select cases to highlight the effects of evaporation during complex accelerations, presented here as a reshock phenomenon. It was found that evaporation increases the circulation over nonevaporating particles postshock. Evaporation was also shown to change the postshock Atwood number. Reshock showed that the multiphase instabilities exhibited additional circulation deposition over the dusty gas approximation. Mixing measures were found to be affected by evaporation, with the most significant effects occurring after reshock.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.

References

Figures

Grahic Jump Location
Fig. 1

Analytical solutions versus FLAG simulations for nonevaporating particles: particle velocity (left) and particle temperature (right)

Grahic Jump Location
Fig. 2

Two-axis plot validation of the saturation temperature and vapor concentration for the evaporation model

Grahic Jump Location
Fig. 3

(a) The initial particle field and (b) the initial gas field. Note: This shows how the seeded and nonseeded phase are initialized separately despite being the same gas. This also can be used to show the DG case initial conditions.

Grahic Jump Location
Fig. 4

Evolution of the dusty gas and evaporation particles instabilities with respect to nondimensional time, τ

Grahic Jump Location
Fig. 5

The qualitative scalar fields for the physical cases at τ = 200

Grahic Jump Location
Fig. 6

Circulation over time for the physical cases

Grahic Jump Location
Fig. 7

Integral mixed width over time for the physical cases

Grahic Jump Location
Fig. 8

The qualitative scalar fields for the diffusion modifier cases at τ = 200

Grahic Jump Location
Fig. 9

The circulation over time for the diffusion modifier cases

Grahic Jump Location
Fig. 10

Integral mixed width over time for the diffusion modifier cases

Grahic Jump Location
Fig. 11

The qualitative scalar fields for the saturation pressure modifier cases at τ = 200

Grahic Jump Location
Fig. 12

Circulation over time for the saturation pressure modifier cases

Grahic Jump Location
Fig. 13

Integral mixed width over time for the saturation pressure modifier cases

Grahic Jump Location
Fig. 14

Evolution of the post-reshock interface with respect to nondimensional time, τR. Like Fig. 4, the particle cases are represented by the exemplar case.

Grahic Jump Location
Fig. 15

Density of the reshock cases with respect to τR

Grahic Jump Location
Fig. 16

Temperatures of the reshock cases with respect to τR

Grahic Jump Location
Fig. 17

Particle temperatures of the reshock cases with respect to τR

Grahic Jump Location
Fig. 18

Vorticity fields of the reshock cases with respect to τR

Grahic Jump Location
Fig. 19

Circulation over time for the reshock cases

Grahic Jump Location
Fig. 20

Integral mixed width over time for the reshock cases

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In