Comparison of macro- and microscopic solutions of the Riemann problem II. Two-phase shock tube

The Riemann problem is one of the basic building blocks for numerical methods in computational fluid mechanics. Nonetheless, there are still open questions and gaps in theory and modelling for situations with complex thermodynamic behavior. In this series, we compare numerical solutions of the macroscopic flow equations with molecular dynamics simulation data. To enable molecular dynamics for sufficiently large scales in time and space, we selected the truncated and shifted Lennard-Jones potential for which also highly accurate equations of state are available. A comparison of a two-phase Riemann problem is shown, which involves a liquid and a vapor phase, with an undergoing phase transition. The loss of hyperbolicity allows for the occurrence of anomalous wave structures. We successfully compare the molecular dynamics data with two macroscopic numerical solutions obtained by either assuming local phase equilibrium or by imposing a kinetic relation and allowing for metastable states.

[1]  Gerald Warnecke,et al.  Exact solutions to the Riemann problem for compressible isothermal Euler equations for two-phase flows with and without phase transition , 2012 .

[2]  Claus-Dieter Munz,et al.  A contribution to the construction of diffusion fluxes for finite volume and discontinuous Galerkin schemes , 2007, J. Comput. Phys..

[3]  Grégoire Allaire,et al.  MODELLING AND SIMULATION OF LIQUID-VAPOR PHASE TRANSITION IN COMPRESSIBLE FLOWS BASED ON THERMODYNAMICAL EQUILIBRIUM ∗ , 2012 .

[4]  Siegfried Müller,et al.  The Riemann Problem for the Euler Equations with Nonconvex and Nonsmooth Equation of State: Construction of Wave Curves , 2006, SIAM J. Sci. Comput..

[5]  Claus-Dieter Munz,et al.  A sharp interface method for compressible liquid-vapor flow with phase transition and surface tension , 2015, J. Comput. Phys..

[6]  P. Raviart,et al.  Numerical Approximation of Hyperbolic Systems of Conservation Laws , 1996, Applied Mathematical Sciences.

[7]  Brij Lal,et al.  A textbook of sound , 1995 .

[8]  R. Saurel,et al.  Dynamic relaxation processes in compressible multiphase flows. Application to evaporation phenomena , 2013 .

[9]  R. Saurel,et al.  A simple and fast phase transition relaxation solver for compressible multicomponent two-phase flows , 2017 .

[10]  Stefan Becker,et al.  ls1 mardyn: The massively parallel molecular dynamics code for large systems , 2014, Journal of chemical theory and computation.

[11]  C. Rohde,et al.  On Riemann solvers and kinetic relations for isothermal two-phase flows with surface tension , 2016, Zeitschrift für angewandte Mathematik und Physik.

[12]  C. Munz,et al.  Approximate Riemann solver for compressible liquid vapor flow with phase transition and surface tension , 2017, Computers & Fluids.

[13]  Claus-Dieter Munz,et al.  Explicit Discontinuous Galerkin methods for unsteady problems , 2012 .

[14]  C. Rohde,et al.  THE SHARP-INTERFACE APPROACH FOR FLUIDS WITH PHASE CHANGE: RIEMANN PROBLEMS AND GHOST FLUID TECHNIQUES , 2007 .

[15]  Dominique Pelletier,et al.  Perspective on the geometric conservation law and finite element methods for ALE simulations of incompressible flow , 2009, J. Comput. Phys..

[16]  J. Peraire,et al.  Sub-Cell Shock Capturing for Discontinuous Galerkin Methods , 2006 .

[17]  Claus-Dieter Munz,et al.  Comparison of macro- and microscopic solutions of the Riemann problem I. Supercritical shock tube and expansion into vacuum , 2020, J. Comput. Phys..

[18]  D. Bedeaux,et al.  Interface film resistivities for heat and mass transfers-integral relations verified by non-equilibrium molecular dynamics. , 2006, The journal of physical chemistry. B.

[19]  Joachim Gross,et al.  Analysis of Interfacial Transport Resistivities of Pure Components and Mixtures Based on Density Functional Theory , 2015 .

[20]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[21]  M. Dumbser,et al.  On the use of tabulated equations of state for multi-phase simulations in the homogeneous equilibrium limit , 2019, Shock Waves.

[22]  Richard Saurel,et al.  Modelling evaporation fronts with reactive Riemann solvers , 2005 .

[23]  Thomas Ertl,et al.  Challenges and Opportunities using Software-defined Visualization in MegaMol , 2017 .

[24]  Rémi Abgrall,et al.  Modelling phase transition in metastable liquids: application to cavitating and flashing flows , 2008, Journal of Fluid Mechanics.

[25]  Guangwei Yuan,et al.  A robust and contact resolving Riemann solver on unstructured mesh, Part II, ALE method , 2014, J. Comput. Phys..

[26]  Ferdinand Thein,et al.  On the Impossibility of First-Order Phase Transitions in Systems Modeled by the Full Euler Equations , 2019, Entropy.

[27]  Dick Bedeaux,et al.  Nonequilibrium thermodynamics of interfaces using classical density functional theory. , 2008, The Journal of chemical physics.

[28]  Dominique Pelletier,et al.  Geometric Conservation Law and Finite Element Methods for 3D unsteady Simulations of Incompressible Flow , 2008 .

[29]  M. Ishii,et al.  Thermo-Fluid Dynamics of Two-Phase Flow , 2007 .

[30]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[31]  D. Bedeaux,et al.  Non-equilibrium Thermodynamics of Heterogeneous Systems , 2008, Series on Advances in Statistical Mechanics.

[32]  Claus-Dieter Munz,et al.  Shock Capturing for Discontinuous Galerkin Methods using Finite Volume Subcells , 2014 .

[33]  Hans Hasse,et al.  Vapor−Liquid Interface of the Lennard-Jones Truncated and Shifted Fluid: Comparison of Molecular Simulation, Density Gradient Theory, and Density Functional Theory , 2018, The Journal of Physical Chemistry C.

[34]  Wolfgang Dahmen,et al.  Exact Riemann solution for the Euler equations with nonconvex and nonsmooth equation of state , 2005 .

[35]  C. Munz,et al.  An Approximate Riemann Solver for Advection–Diffusion Based on the Generalized Riemann Problem , 2019, Communications on Applied Mathematics and Computation.

[36]  Phillip Colella,et al.  Efficient Solution Algorithms for the Riemann Problem for Real Gases , 1985 .

[37]  G. Meier,et al.  Complete adiabatic evaporation of highly superheated liquid jets , 1992, Journal of Fluid Mechanics.

[38]  G. Warnecke,et al.  Bubbles in liquids with phase transition , 2012 .

[39]  S. Fechter.  Compressible multi-phase simulation at extreme conditions using a discontinuous Galerkin scheme , 2015 .

[40]  David A. Kopriva,et al.  Discontinuous Galerkin spectral element approximations on moving meshes , 2011, J. Comput. Phys..

[41]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[42]  Thomas Ertl,et al.  MegaMol—A Prototyping Framework for Particle-Based Visualization , 2015, IEEE Transactions on Visualization and Computer Graphics.

[43]  J. K. Knowles,et al.  Kinetic relations and the propagation of phase boundaries in solids , 1991 .

[44]  G. Yadigaroglu,et al.  Explosive vaporization of superheated liquids by boiling fronts , 2001 .

[45]  Alexandre Joel Chorin,et al.  Riemann Problems for Reacting Gas, with Applications to Transition , 1982 .

[46]  H. Bruce Stewart,et al.  Two-phase flow: Models and methods , 1984 .

[47]  Christian Rohde,et al.  A relaxation Riemann solver for compressible two-phase flow with phase transition and surface tension , 2015 .

[48]  Claus-Dieter Munz,et al.  Efficient Parallelization of a Shock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells , 2016, Journal of Scientific Computing.

[49]  R. Lewis,et al.  Low-storage, Explicit Runge-Kutta Schemes for the Compressible Navier-Stokes Equations , 2000 .

[50]  A. B. Wood,et al.  A textbook of sound , 1930 .

[51]  Joseph E. Shepherd,et al.  Evaporation waves in superheated dodecane , 1999, Journal of Fluid Mechanics.

[52]  Roland Span,et al.  Equation of State for the Lennard-Jones Truncated and Shifted Model Fluid , 2015 .

[53]  Jadran Vrabec,et al.  Communication: Evaporation: Influence of heat transport in the liquid on the interface temperature and the particle flux. , 2016, The Journal of chemical physics.

[54]  James K. Knowles,et al.  On the driving traction acting on a surface of strain discontinuity in a continuum , 1990 .

[55]  Lev Truskinovsky,et al.  Kinks versus Shocks , 1993 .

[56]  C. Zeiler.  Liquid Vapor Phase Transitions: Modeling, Riemann Solvers and Computation , 2016 .

[57]  Jadran Vrabec,et al.  Evaporation sampled by stationary molecular dynamics simulation. , 2019, The Journal of chemical physics.

[58]  Claus-Dieter Munz,et al.  An explicit discontinuous Galerkin scheme with local time-stepping for general unsteady diffusion equations , 2008, J. Comput. Phys..

[59]  R. Menikoff,et al.  The Riemann problem for fluid flow of real materials , 1989 .

[60]  James K. Knowles,et al.  Evolution of Phase Transitions: A Continuum Theory , 2006 .

[61]  Michael Dumbser,et al.  High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: Viscous heat-conducting fluids and elastic solids , 2015, J. Comput. Phys..

[62]  Michael M. Resch,et al.  TweTriS: Twenty trillion-atom simulation , 2019, Int. J. High Perform. Comput. Appl..

[63]  Karsten Schatz,et al.  MegaMol – a comprehensive prototyping framework for visualizations , 2019, The European Physical Journal Special Topics.

[64]  H. Hasse,et al.  Equation of state for the Lennard-Jones truncated and shifted fluid with a cut-off radius of 2.5 σ based on perturbation theory and its applications to interfacial thermodynamics , 2018 .

[65]  Guangwei Yuan,et al.  A robust and contact resolving Riemann solver on unstructured mesh, Part I, Euler method , 2014, J. Comput. Phys..

[66]  Roland Span,et al.  How well does the Lennard-Jones potential represent the thermodynamic properties of noble gases? , 2017 .