Research Page

Adviser: Dr. Daniel Meiron
Adviser: Dr. Oleg Schilling (LLNL)
Thesis Title: Simulation and Analysis of
Two- and Three-Dimensional Single-Mode Richtmyer-Meshkov
Instability using Weighted Essentially Non-Oscillatory and Vortex Methods


Research Overview  

  • A central objective of my research is to investigates the physics of the Richtmyer-Meshkov instability.
  • The Richtmyer-Meshkov instability is a fundamental fluid instability that occurs when perturbations on an interface separating gases with different properties grow following the passage of a shock. This instability is of great fundamental interest in fluid dynamics, as well as of interest to inertial confinement fusion, and to supernovae dynamics.
  • The instability derives its name from the linear instability analysis and numerical simulations of Richtmyer [Comm. Pure Appl Math. 8, 297 (1960)], and the shock tube experiments of Meshkov [Sov. Fluid Dyn. 4, 101 (1969)]
  • The goal of this research is two-fold;
    1. understand the evolution of the instability, model the growth of the mixing layer in the nonlinear phase and following reshock, as well as predict the statistical properties and dynamics of turbulent mixing induced by the instability.
    2. understand the large scale dynamics of the instability, including the bubble and spike dynamics, interface stretching, and large scale properties of the instability.
  • Two numerical methods are employed in the investigation as they offer complementary views of flow features
    1. the compressible Euler equation solver based on the Weighted Essentially Non-Oscillatory (WENO) shock capturing method, is used to determine the dynamics of the mixing layer, as well as the statistical and turbulent properties of mixing
    2. the incompressible vortex method is used to model the bubble and spike dynamics, and the interface evolution
    • Thesis Title: Investigation of the Richtmyer-Meshkov Instability in Complex Geometries Using the Weighted Essentially Non-Oscillatory method and Vortex Methods
      Research Goals
       
      A central objective of the present work is to establish a systematic procedure to investigate the dynamics of the mixing process induced by the Richtmyer-Meshkov instability, and more generally by complex hydrodynamic flows. The methods used are adapted from classical investigations of turbulence and turbulent mixing, and synthesize high-resolution numerical simulation data, theoretical models for instability growth, and available experimental data. This procedure results in:
      1. the application of a modern, high-resolution, flexible numerical method that has been validated against available experimental data;
      2. a numerical database that provides quantities that can be compared to model predictions and to experimental measurements, as well as quantities that have not been modeled (or are difficult to model) or are not available experimentally;
      3. numerical data for configurations extended to times beyond what is possible to achieve experimentally, or for configurations that are difficult to achieve experimentally;
      4. a systematic understanding of the important effects of spatial resolution and formal order of the method on quantities of interest to modeling the instability evolution and mixing.

    Program Goals

    This study is part of a larger, longer-term program aimed at:
    1. developing improved theoretical models for instability growth in the nonlinear regime, as well as for the evolution following reshock;
    2. investigating closure models for ensemble-averaged descriptions of turbulent transport and mixing, as well as for the development of subgrid-scale models for large-eddy simulations;
    3. improving the numerical methods used to simulate complex hydrodynamic flows induced by shocks;
    4. aiding the design of new experimental configurations and new experimental diagnostics.



    Publications

    • "High-resolution simulations and modeling of reshocked single-mode Richtmyer-Meshkov instability. I. Comparison to experimental data and to amplitude growth model predictions" by M. Latini, O. Schilling, and W.-S. Don, Physics of Fluids, 2006 (submitted)
    • "High-resolution simulations and modeling of reshocked single-mode Richtmyer-Meshkov instability. II. Physics of reshock and mixing" by O. Schilling, M. Latini, and W.-S. Don, Physics of Fluids, 2006 (submitted) 
    • "Effects of WENO flux reconstruction order and spatial resolution on reshocked two-dimensional Richtmyer-Meshkov instability" by M. Latini, O. Schilling, and W.-S. Don, Journal of Computational Physics, 2006 (submitted)
    • Weighted Essentially Non-Oscillatory Simulations and Modeling of Complex Hydrodynamic Flows. Part 1. Regular Shock Refraction  by M. Latini and O. Schilling, Tech. Report UCRLTR205132, Lawrence Livermore National Laboratory, 2004
    •  Weighted Essentially Non-Oscillatory Simulations and Modeling of Complex Hydrodynamic Flows. Part 2. Single-Mode Richtmyer-Meshkov Instability with Reshock  by M. Latini and O. Schilling, Tech. Report, Lawrence Livermore National Laboratory, 2004

    Conference and Poster Presentations

    • Talk: Quantitative analysis of three-dimensional reshocked Richtmyer-Meshkov instability-induced mixing using different orders of WENO flux reconstruction, 58 APS Division of Fluid Dynamics Meeting, Chicago, IL, November 2005 (with O. Schilling and W.-S. Don)
    • Talk:  Assessment of gradient-diffusion closures for modeling turbulent transport in three-dimensional Richtmyer-Meshkov instability-induced mixing with reshock, 58 APS Division of Fluid Dynamics Meeting, Chicago, IL, November 2005 (with O. Schilling and W.-S. Don)
    • Poster: Turbulent Transport Properties of 3D Richtmyer-Meshkov Instability with Reshock, HEDP Summer School, Berkeley, CA, August 2005, with O. Schilling and W.-S. Don
    • Talk:  Investigation of the Richtmyer-Meshkov Instability in Complex Geometries with vortex methods, 57 APS Division of Fluid Dynamics Meeting, Seattle, WA, November 2004 (with Dr. Daniel Meiron, Dr. Paul Dimotakis, Dr. Oleg Schilling)
    • Talk: Investigation of three-dimensional Richtmyer-Meshkov instability-induced mixing with reshock, 57APS Division of Fluid Dynamics Meeting, Seattle, WA, November 2004 (with Dr. Oleg Schilling, Dr. Wai-Sun Don)
    • Poster: Investigation of the Large-Scale and Statistical Properties of Richtmyer-Meshkov Instability-Induced Mixing Richtmyer-Meshkov Instability in Complex Geometries, Caltech ASC Center Review, October 2004 (with Dr. Oleg Schilling)
    • Talk:  Investigation of the Large-Scale and Statistical Properties of Richtmyer-Meshkov Instability-Induced Mixing Richtmyer-Meshkov Instability in Complex Geometries, 9th International Workshop on the Physics of Compressible Turbulence and Mixing, Cambridge, UK, July 2004 (with Dr. Oleg Schilling)
    • Invited Talk:  Investigation of Richtmyer-Meshkov Instability in Complex Geometries, Division of Applied Mathematics, Brown University, Providence, RI, March 2004
    • Talk: Investigation of Richtmyer-Meshkov Instability in Complex Geometries, 56 APS Division of Fluid Dynamics Meeting, Rutherford, NJ, November 2003 (with Dr. Oleg Schilling, Dr. Wai-Sun Don)
    • Talk: WENO Simulations of Single- and Multi-Mode Richtmyer-Meshkov Instability-Induced Mixing, APS Division of Fluid Dynamics Meeting, Rutherford, NJ, November 2003 (with Dr. Oleg Schilling, Dr. Wai-Sun Don)
    • Poster: Investigation of the Oblique Shock Richtmyer-Meshkov Instability, ASC Center Review, October 2003 (with Dr. Oleg Schilling, Dr. Wai-Sun Don)




    Other Research or Teaching Activities


    • Lecture Notes on Hyperbolic systems (AE 101,  10/2004, PDF)
    • QR Algorithm presentation (ACM student seminar 02/2004, PDF)
    • QR Algorithm report (02/2004 PDF)

    Main Page

    Last updated: May 12, 2006
    comments e-mail: mlatini@acm.caltech.edu