Catalogue Search | MBRL
Search Results Heading
Explore the vast range of titles available.
MBRLSearchResults
-
DisciplineDiscipline
-
Is Peer ReviewedIs Peer Reviewed
-
Item TypeItem Type
-
SubjectSubject
-
YearFrom:-To:
-
More FiltersMore FiltersSourceLanguage
Done
Filters
Reset
11
result(s) for
"Groth, Frederick"
Sort by:
N-body simulations of dark matter-baryon interactions
2025
Dark matter (DM) particles can interact with particles of the standard model. Although there are a number of constraints derived from direct and indirect detection experiments, the evolution of astrophysical objects could offer a promising probe. Obtaining predictions is challenging and primarily limited by our ability to simulate scattering between DM and baryonic particles within N-body and hydrodynamics simulations. We have developed the first scheme allowing for the simulation of these interacting dark matter (IDM) models, accurately accounting for their angular and velocity dependence, as well as the mass ratio between the DM and baryonic scattering partners. To describe DM-baryon interactions, we used an N-body code together with its implementation of smoothed-particle hydrodynamics and meshless finite mass. The interaction is realised in a pairwise fashion by creating a virtual scattering partner from the baryonic particle and allowing it to interact with a DM particle using a scattering routine initially developed for self-interacting dark matter (SIDM). After the interaction, the virtual particle is rejoined with the baryonic particle, fulfilling the requirements of energy and momentum conservation. Through several test problems, we demonstrated that we are able to reproduce the analytic solutions with our IDM scheme. This includes a test for scattering with a physical mass ratio of 1:1000, which is beyond the limits of SIDM simulations. We comment on various numerical aspects and challenges, and we describe the limitations of our numerical scheme. Furthermore, we study the impact of IDM on halo formation with a collapsing over-density. We find that it is possible to accurately model IDM within N-body and hydrodynamics simulations commonly used in astrophysics. Finally, our scheme allows for novel predictions to be made and new constraints on DM-baryon scattering to be set.
The Cosmological Simulation Code OpenGadget3 -- Implementation of Meshless Finite Mass
2023
Subsonic turbulence plays a major role in determining properties of the intra cluster medium (ICM). We introduce a new Meshless Finite Mass (MFM) implementation in OpenGadget3 and apply it to this specific problem. To this end, we present a set of test cases to validate our implementation of the MFM framework in our code. These include but are not limited to: the soundwave and Kepler disk as smooth situations to probe the stability, a Rayleigh-Taylor and Kelvin-Helmholtz instability as popular mixing instabilities, a blob test as more complex example including both mixing and shocks, shock tubes with various Mach numbers, a Sedov blast wave, different tests including self-gravity such as gravitational freefall, a hydrostatic sphere, the Zeldovich-pancake, and a \\(10^{15}M_{\\odot}\\) galaxy cluster as cosmological application. Advantages over SPH include increased mixing and a better convergence behavior. We demonstrate that the MFM-solver is robust, also in a cosmological context. We show evidence that the solver performs extraordinarily well when applied to decaying subsonic turbulence, a problem very difficult to handle for many methods. MFM captures the expected velocity power spectrum with high accuracy and shows a good convergence behavior. Using MFM or SPH within OpenGadget3 leads to a comparable decay in turbulent energy due to numerical dissipation. When studying the energy decay for different initial turbulent energy fractions, we find that MFM performs well down to Mach numbers \\(\\mathcal{M}\\approx 0.01\\). Finally, we show how important the slope limiter and the energy-entropy switch are to control the behavior and the evolution of the fluids.
Turbulence in Simulated Local Cluster Analogs: One-to-one comparisons between SLOW and XRISM/Hitomi
by
Seidel, Benjamin A
,
Vladutescu-Zopp, Stephan
,
Sorce, Jenny G
in
Galactic clusters
,
Simulation
,
Turbulence
2026
The XRISM Resolve X-ray spectrometer allows to gain detailed insight into gas motions of the intra cluster medium (ICM) of galaxy clusters. Current simulation studies focus mainly on statistical comparisons, making the comparison to the currently still small number of clusters difficult due to unknown selection effects. This study aims to bridge this gap, using simulated counterparts of Coma, Virgo, and Perseus from the SLOW constrained simulations. These clusters show excellent agreement in their properties and dynamical state with observations, thus providing an ideal testbed to understand the processes shaping the properties of the ICM. We find that the simulations match the order of the amount of turbulence for the three considered clusters, Coma being the most active, followed by Perseus, while Virgo is very relaxed. Typical turbulent velocities are a few \\(\\approx100\\) km s\\(^{-1}\\), very close to observed values. The resulting turbulent pressure support is \\(\\approx1\\%\\) for Virgo, \\(\\approx 6\\%\\) for Perseus, and \\(\\approx 8\\%\\) for Coma within the central \\(1-2\\%\\) of \\(R_{200}\\). Compared to previous simulations and observations, measured velocities and turbulent pressure support are on average lower, in line with XRISM findings, thus indicating the importance of selection effects.
Vortex-p: a Helmholtz-Hodge and Reynolds decomposition algorithm for particle-based simulations
by
Vallés-Pérez, David
,
Planelles, Susana
,
Groth, Frederick
in
Algorithms
,
Computer simulation
,
Decomposition
2024
Astrophysical turbulent flows display an intrinsically multi-scale nature, making their numerical simulation and the subsequent analyses of simulated data a complex problem. In particular, two fundamental steps in the study of turbulent velocity fields are the Helmholtz-Hodge decomposition (compressive+solenoidal; HHD) and the Reynolds decomposition (bulk+turbulent; RD). These problems are relatively simple to perform numerically for uniformly-sampled data, such as the one emerging from Eulerian, fix-grid simulations; but their computation is remarkably more complex in the case of non-uniformly sampled data, such as the one stemming from particle-based or meshless simulations. In this paper, we describe, implement and test vortex-p, a publicly available tool evolved from the vortex code, to perform both these decompositions upon the velocity fields of particle-based simulations, either from smoothed particle hydrodynamics (SPH), moving-mesh or meshless codes. The algorithm relies on the creation of an ad-hoc adaptive mesh refinement (AMR) set of grids, on which the input velocity field is represented. HHD is then addressed by means of elliptic solvers, while for the RD we adapt an iterative, multi-scale filter. We perform a series of idealised tests to assess the accuracy, convergence and scaling of the code. Finally, we present some applications of the code to various SPH and meshless finite-mass (MFM) simulations of galaxy clusters performed with OpenGadget3, with different resolutions and physics, to showcase the capabilities of the code.
Turbulent Pressure Support in Galaxy Clusters -- Impact of the Hydrodynamical Solver
by
Steinwandel, Ulrich P
,
Vallés-Pérez, David
,
Groth, Frederick
in
Galactic clusters
,
Red shift
,
Simulation
2024
The amount of turbulent pressure in galaxy clusters is still debated, especially as for the impact of the dynamical state and the hydro-method used for simulations. We study the turbulent pressure fraction in the intra cluster medium of massive galaxy clusters. We aim to understand the impact of the hydrodynamical scheme, analysis method, and dynamical state on the final properties of galaxy clusters from cosmological simulations. We perform non-radiative simulations of a set of zoom-in regions of seven galaxy clusters with Meshless Finite Mass (MFM) and Smoothed Particle Hydrodynamics (SPH). We use three different analysis methods based on: \\((i)\\) the deviation from hydrostatic equilibrium, \\((ii)\\) the solenoidal velocity component obtained by a Helmholtz-Hodge decomposition, and \\((iii)\\) the small-scale velocity obtained through a multi-scale filtering approach. We split the sample of simulated clusters into active and relaxed clusters. Our simulations predict an increased turbulent pressure fraction for active compared to relaxed clusters. This is especially visible for the velocity-based methods. For these, we also find increased turbulence for the MFM simulations compared to SPH, consistent with findings from more idealized simulations. The predicted non-thermal pressure fraction varies between a few percent for relaxed clusters and \\(\\approx13\\%\\) for active ones within the cluster center and increases towards the outskirts. No clear trend with redshift is visible. Our analysis quantitatively assesses the importance played by the hydrodynamical scheme and the analysis method to determine the non-thermal/turbulent pressure fraction. While our setup is relatively simple (non-radiative runs), our simulations show agreement with previous, more idealized simulations, and make a step further toward the understanding of turbulence.
Turbulence in Simulated Local Cluster Analogs: one-to-one comparisons between SLOW and XRISM/Hitomi
by
Seidel, Benjamin A
,
Vladutescu-Zopp, Stephan
,
Sorce, Jenny G
in
Galactic clusters
,
Simulation
,
Turbulence
2025
The XRISM Resolve X-ray spectrometer allows to gain detailed insight into gas motions of the intra cluster medium (ICM) of galaxy clusters. Current simulation studies focus mainly on statistical comparisons, making the comparison to the currently still small number of clusters difficult due to unknown selection effects. This study aims to bridge this gap, using simulated counterparts of Coma, Virgo, and Perseus from the SLOW constrained simulations. These clusters show excellent agreement in their properties and dynamical state with observations, thus providing an ideal testbed to understand the processes shaping the properties of the ICM. We find that the simulations match the order of the amount of turbulence for the three considered clusters, Coma being the most active, followed by Perseus, while Virgo is very relaxed. Typical turbulent velocities are a few \\(100\\) km s\\(^-1\\), very close to observed values. The resulting turbulent pressure support is \\(1\\%\\) for Virgo and \\( 3-4\\%\\) for Perseus and Coma within the central \\(1-2\\%\\) of \\(R_200\\). Compared to previous simulations and observations, measured velocities and turbulent pressure support are consistently lower, in line with XRISM findings, thus indicating the importance of selection effects.
Turbulent Pressure Support in Galaxy Clusters -- Impact of the Hydrodynamical Solver
by
Steinwandel, Ulrich P
,
Vallés-Pérez, David
,
Groth, Frederick
in
Cluster analysis
,
Galactic clusters
,
Impact analysis
2024
The amount of turbulent pressure in galaxy clusters is still debated, especially as for the impact of the dynamical state and the hydro-method used for simulations. We study the turbulent pressure fraction in the intra cluster medium of massive galaxy clusters. We aim to understand the impact of the hydrodynamical scheme, analysis method, and dynamical state on the final properties of galaxy clusters from cosmological simulations. We perform non-radiative simulations of a set of zoom-in regions of seven galaxy clusters with Meshless Finite Mass (MFM) and Smoothed Particle Hydrodynamics (SPH). We use three different analysis methods based on: \\((i)\\) the deviation from hydrostatic equilibrium, \\((ii)\\) the solenoidal velocity component obtained by a Helmholtz-Hodge decomposition, and \\((iii)\\) the small-scale velocity obtained through a multi-scale filtering approach. We split the sample of simulated clusters into active and relaxed clusters. Our simulations predict an increased turbulent pressure fraction for active compared to relaxed clusters. This is especially visible for the velocity-based methods. For these, we also find increased turbulence for the MFM simulations compared to SPH, consistent with findings from more idealized simulations. The predicted non-thermal pressure fraction varies between a few percent for relaxed clusters and \\(\\approx13\\%\\) for active ones within the cluster center and increases towards the outskirts. No clear trend with redshift is visible. Our analysis quantitatively assesses the importance played by the hydrodynamical scheme and the analysis method to determine the non-thermal/turbulent pressure fraction. While our setup is relatively simple (non-radiative runs), our simulations show agreement with previous, more idealized simulations, and make a step further toward the understanding of turbulence.
Stratigraphy of the Chinle Formation in the San Rafael Swell, Utah
by
Groth, Frederick A
in
Geology
1955
The unit formerly mapped as Shinarump conglomerate in the San Rafael Swell, Utah, is in reality a medial unit of the Chinle formation, and the name Mossback sandstone is applied in accordance with the present U.S. Geological Survey nomenclature. The overlying interval, formerly mapped as Chinle formation, is renamed Upper Chinle. Because of its prominence and consistency, the Mossback sandstone is raised to member status in the Chinle formation, and the name Monitor Butte is applied to the Chinle interval below the Mossback, which is also considered to be a member of the Chinle formation. The lithology of the Upper Chinle, Mossback sandstone member and the Monitor Butte member is described in detail. A correlation stereogram, showing the correlation of the various member units of the Chinle formation around the San Rafael Swell is included, and is discussed. Stokes' pediment concept of deposition for the Shinarump conglomerate is described, and applied to the Mossback sandstone. Generalizations concerning the climate and environment of deposition are presented.
Dissertation