Influence of adaptive mesh refinement and the hydro solver on shear-induced mass stripping in a minor-merger scenario
We compare two different codes for simulations of cosmological structure - formation to investigate the sensitivity of hydrodynamical instabilities - to numerics, in particular, the hydro solver and the application of - adaptive mesh refinement (AMR). As a simple test problem, we consider an - init...
Saved in:
| Main Authors: | , , , , |
|---|---|
| Format: | Article (Journal) |
| Language: | English |
| Published: |
2015
|
| In: |
Astronomy and computing
Year: 2014, Volume: 9, Pages: 49-63 |
| ISSN: | 2213-1345 |
| DOI: | 10.1016/j.ascom.2014.11.003 |
| Online Access: | Verlag, lizenzpflichtig, Volltext: https://doi.org/10.1016/j.ascom.2014.11.003 Verlag, lizenzpflichtig, Volltext: http://adsabs.harvard.edu/abs/2015A%26C.....9...49S 
|
| Author Notes: | W. Schmidt, J. Schulz, L. Iapichino, F. Vazza, A.S. Almgren |
| Summary: | We compare two different codes for simulations of cosmological structure - formation to investigate the sensitivity of hydrodynamical instabilities - to numerics, in particular, the hydro solver and the application of - adaptive mesh refinement (AMR). As a simple test problem, we consider an - initially spherical gas cloud in a wind, which is an idealized model for - the merger of a subcluster or galaxy with a big cluster. Based on an - entropy criterion, we calculate the mass stripping from the subcluster - as a function of time. Moreover, the turbulent velocity field is - analyzed with a multi-scale filtering technique. We find remarkable - differences between the commonly used PPM solver with directional - splitting in the ENZO code and an unsplit variant of PPM in the NYX - code, which demonstrates that different codes can converge to - systematically different solutions even when using uniform grids. For - the test case of an unbound cloud, AMR simulations reproduce - uniform-grid results for the mass stripping quite well, although the - flow realizations can differ substantially. If the cloud is bound by a - static gravitational potential, however, we find strong sensitivity to - spurious fluctuations which are induced at the cutoff radius of the - potential and amplified by the bow shock. This gives rise to substantial - deviations between uniform-grid and AMR runs performed with ENZO, while - the mass stripping in NYX simulations of the subcluster is nearly - independent of numerical resolution and AMR. Although many factors - related to numerics are involved, our study indicates that unsplit - solvers with advanced flux limiters help to reduce grid effects and to - keep numerical noise under control, which is important for - hydrodynamical instabilities and turbulent flows. |
|---|---|
| Item Description: | Available online 21 November 2014 Gesehen am 06.07.2020 |
| Physical Description: | Online Resource |
| ISSN: | 2213-1345 |
| DOI: | 10.1016/j.ascom.2014.11.003 |