skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: RICH: OPEN-SOURCE HYDRODYNAMIC SIMULATION ON A MOVING VORONOI MESH

Abstract

We present here RICH, a state-of-the-art two-dimensional hydrodynamic code based on Godunov’s method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time-advancement schemeS as well as a novel parallelization scheme based on Voronoi tessellation. Using our code, we study the pros and cons of a moving mesh (in comparison to a static mesh). We also compare its accuracy to other codes. Specifically, we show that our implementation of external sources and time-advancement scheme is more accurate and robust than is AREPO when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Lloyd iterations) and its effects. We find that in most cases a moving mesh gives better results than a static mesh, but it is not universally true. In the case where matter moves in one way and a sound wave is traveling in the other way (such that relative to the grid the wave is not moving) a static mesh gives better results than a moving mesh. We perform an analytic analysis for finite difference schemes that reveals that a Lagrangianmore » simulation is better than a Eulerian simulation in the case of a highly supersonic flow. Moreover, we show that Voronoi-based moving mesh schemes suffer from an error, which is resolution independent, due to inconsistencies between the flux calculation and the change in the area of a cell. Our code is publicly available as open source and designed in an object-oriented, user-friendly way that facilitates incorporation of new algorithms and physical processes.« less

Authors:
; ;  [1]
  1. Racah Institute of Physics, The Hebrew University, 91904, Jerusalem (Israel)
Publication Date:
OSTI Identifier:
22340087
Resource Type:
Journal Article
Resource Relation:
Journal Name: Astrophysical Journal, Supplement Series; Journal Volume: 216; Journal Issue: 2; Other Information: Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; ACCURACY; ALGORITHMS; ANALYTIC FUNCTIONS; COMPARATIVE EVALUATIONS; COMPUTERIZED SIMULATION; HYDRODYNAMICS; INTERPOLATION; LAGRANGIAN FUNCTION; SOUND WAVES; SUPERSONIC FLOW; TWO-DIMENSIONAL CALCULATIONS

Citation Formats

Yalinewich, Almog, Steinberg, Elad, and Sari, Re’em. RICH: OPEN-SOURCE HYDRODYNAMIC SIMULATION ON A MOVING VORONOI MESH. United States: N. p., 2015. Web. doi:10.1088/0067-0049/216/2/35.
Yalinewich, Almog, Steinberg, Elad, & Sari, Re’em. RICH: OPEN-SOURCE HYDRODYNAMIC SIMULATION ON A MOVING VORONOI MESH. United States. doi:10.1088/0067-0049/216/2/35.
Yalinewich, Almog, Steinberg, Elad, and Sari, Re’em. Sun . "RICH: OPEN-SOURCE HYDRODYNAMIC SIMULATION ON A MOVING VORONOI MESH". United States. doi:10.1088/0067-0049/216/2/35.
@article{osti_22340087,
title = {RICH: OPEN-SOURCE HYDRODYNAMIC SIMULATION ON A MOVING VORONOI MESH},
author = {Yalinewich, Almog and Steinberg, Elad and Sari, Re’em},
abstractNote = {We present here RICH, a state-of-the-art two-dimensional hydrodynamic code based on Godunov’s method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time-advancement schemeS as well as a novel parallelization scheme based on Voronoi tessellation. Using our code, we study the pros and cons of a moving mesh (in comparison to a static mesh). We also compare its accuracy to other codes. Specifically, we show that our implementation of external sources and time-advancement scheme is more accurate and robust than is AREPO when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Lloyd iterations) and its effects. We find that in most cases a moving mesh gives better results than a static mesh, but it is not universally true. In the case where matter moves in one way and a sound wave is traveling in the other way (such that relative to the grid the wave is not moving) a static mesh gives better results than a moving mesh. We perform an analytic analysis for finite difference schemes that reveals that a Lagrangian simulation is better than a Eulerian simulation in the case of a highly supersonic flow. Moreover, we show that Voronoi-based moving mesh schemes suffer from an error, which is resolution independent, due to inconsistencies between the flux calculation and the change in the area of a cell. Our code is publicly available as open source and designed in an object-oriented, user-friendly way that facilitates incorporation of new algorithms and physical processes.},
doi = {10.1088/0067-0049/216/2/35},
journal = {Astrophysical Journal, Supplement Series},
number = 2,
volume = 216,
place = {United States},
year = {Sun Feb 01 00:00:00 EST 2015},
month = {Sun Feb 01 00:00:00 EST 2015}
}