1,623 research outputs found
Artificial viscosity model to mitigate numerical artefacts at fluid interfaces with surface tension
The numerical onset of parasitic and spurious artefacts in the vicinity of uid interfaces with surface tension is an important and well-recognised problem with respect to the accuracy and numerical stability of interfacial ow simulations. Issues of particular interest are spurious capillary waves, which are spatially underresolved by the computational mesh yet impose very restrictive time-step requirements, as well as parasitic currents, typically the result of a numerically unbalanced curvature evaluation. We present an arti cial viscosity model to mitigate numerical artefacts at surface-tension-dominated interfaces without adversely a ecting the accuracy of the physical solution. The proposed methodology computes an additional interfacial shear stress term, including an interface viscosity, based on the local ow data and uid properties that reduces the impact of numerical artefacts and dissipates underresolved small scale interface movements. Furthermore, the presented methodology can be readily applied to model surface shear viscosity, for instance to simulate the dissipative e ect of surface-active substances adsorbed at the interface. The presented analysis of numerical test cases demonstrates the e cacy of the proposed methodology in diminishing the adverse impact of parasitic and spurious interfacial artefacts on the convergence and stability of the numerical solution algorithm as well as on the overall accuracy of the simulation results
Galaxy Tracers and Velocity Bias
This paper examines several methods of tracing galaxies in N-body simulations
and their effects on the derived galaxy statistics, especially measurements of
velocity bias. Using two simulations with identical initial conditions, one
following dark matter only and the other following dark matter and baryons,
both collisionless and collisional methods of tracing galaxies are compared to
one another and against a set of idealized criteria. None of the collisionless
methods proves satisfactory, including an elaborate scheme developed here to
circumvent previously known problems. The main problem is that galactic
overdensities are both secularly and impulsively disrupted while orbiting in
cluster potentials. With dissipation, the baryonic tracers have much higher
density contrasts and much smaller cross sections, allowing them to remain
distinct within the cluster potential. The question remains whether the
incomplete physical model introduces systematic biases. Statistical measures
determined from simulations can vary significantly based solely on the galaxy
tracing method utilized. The two point correlation function differs most on
sub-cluster scales with generally good agreement on larger scales. Pairwise
velocity dispersions show less uniformity on all scales addressed here. All
tracing methods show a velocity bias to varying degrees, but the predictions
are not firm: either the tracing method is not robust or the statistical
significance has not been demonstrated. Though theoretical arguments suggest
that a mild velocity bias should exist, simulation results are not yet
conclusive.Comment: ApJ, in press, 23 pages, plain TeX, 8 of 13 figures included, all
PostScript figures (4.8 MB) available via anonymous ftp from
ftp://astro.princeton.edu/summers/tracers . Also available as POPe-616 on
http://astro.princeton.edu/~library/prep.htm
Mass Estimates of X-Ray Clusters
We use cosmological gas dynamic simulations to investigate the accuracy of
galaxy cluster mass estimates based on X-ray observations. The experiments
follow the formation of clusters in different cosmological models and include
the effects of gravity, pressure gradients, and hydrodynamical shocks. A subset
of our ensemble also allows for feedback of mass and energy from galactic winds
into the intracluster medium. We find that mass estimates based on the
hydrostatic, isothermal beta-model are remarkably accurate when evaluated at
radii where the cluster mean density is between 500-2500 times the critical
density. Applied to 174 artificial ROSAT images constructed from the
simulations, the distribution of the estimated-to-true mass ratio is nearly
unbiased and has a standard deviation of 14-29%. The scatter can be
considerably reduced (to 8-15%) by using an alternative mass estimator that
exploits the tightness of the mass-temperature relation found in the
simulations. The improvement over beta-model estimates is due to the
elimination of the variance contributed by the gas outer slope parameter. We
discuss these findings and their implications for recent measurements of
cluster baryon fractions.Comment: TeX, 24p; 11 Postscript figs. Submitted to the Astrophysical Journa
Galaxy tracers in N-body simulations
Using the method of smoothed particle hydrodynamics, we have modeled the formation of a compact group of galaxies with sufficient resolution to trace galaxies. Radiative cooling allows the baryons to dissipate their thermal energy and collapse to overdensities characteristic of real galaxies. With their cross section greatly reduced, these galaxy tracers remain distinct during cluster formation while their dark matter halos merge. In addition, the number density, the mass distribution function, and even the morphology of these objects are similar to those of observed galaxies. A viable population of galaxy tracers can be unambiguously defined
The Ultimate Halo Mass in a LCDM Universe
In the far future of an accelerating LCDM cosmology, the cosmic web of
large-scale structure consists of a set of increasingly isolated halos in
dynamical equilibrium. We examine the approach of collisionless dark matter to
hydrostatic equilibrium using a large N-body simulation evolved to scale factor
a = 100, well beyond the vacuum--matter equality epoch, a_eq ~ 0.75, and 53/h
Gyr into the future for a concordance model universe (Omega_m ~ 0.3,
Omega_Lambda ~ 0.7). The radial phase-space structure of halos -- characterized
at a < a_eq by a pair of zero-velocity surfaces that bracket a dynamically
active accretion region -- simplifies at a > 10 a_eq when these surfaces merge
to create a single zero-velocity surface, clearly defining the halo outer
boundary, rhalo, and its enclosed mass, mhalo. This boundary approaches a fixed
physical size encompassing a mean interior density ~ 5 times the critical
density, similar to the turnaround value in a classical Einstein-deSitter
model. We relate mhalo to other scales currently used to define halo mass
(m200, mvir, m180b) and find that m200 is approximately half of the total
asymptotic cluster mass, while m180b follows the evolution of the inner zero
velocity surface for a < 2 but becomes much larger than the total bound mass
for a > 3. The radial density profile of all bound halo material is well fit by
a truncated Hernquist profile. An NFW profile provides a somewhat better fit
interior to r200 but is much too shallow in the range r200 < r < rhalo.Comment: 5 pages, 3 figures, submitted to MNRAS letter
Estimation of curvature from volume fractions using parabolic reconstruction on two-dimensional unstructured meshes
This paper proposes a method to estimate the curvature of an interface represented implicitly by discrete volume fractions on an unstructured two-dimensional mesh. The method relies on the computation of local parabolic reconstructions of the interface. The parabolic reconstruction of the interface in a given computational cell is obtained by solving a local non-linear minimisation problem, and only requires additional information from two neighbouring cells. This compactness ensures a robust behaviour on poorly-resolved interfaces. The proposed method is proven to be analogous to the height-function method for Cartesian configurations with consistent heights, and can be interpreted as a generalisation of the height-function method to meshes of any type. Tests are conducted on a range of interfaces with known curvature. The method is shown to converge with mesh refinement with the same order of accuracy as the height-function method for all three types of meshes tested, i.e. Cartesian, triangular, and polygonal
Quantitative Estimates of Environmental Effects on the Star Formation Rate of Disk Galaxies in Clusters of Galaxies
A simple model is constructed to evaluate the change of star formation rate
of a disk galaxy due to environmental effects in clusters of galaxies. Three
effects, (1) tidal force from the potential well of the cluster, (2) increase
of external pressure when the galaxy plows into the intracluster medium, (3)
high-speed encounters between galaxies, are investigated. General analysis
indicates that the star formation rate increases significantly when the
pressure of molecular clouds rises above in yr. The tidal force from the potential well of the cluster increases
pressures of molecular clouds in a disk galaxy infalling towards the cluster
center. Before the galaxy reaches the cluster center, the star formation rate
reaches a maximum. The peak is three to four times larger than the initial
value. If this is the main mechanism of the Butcher-Oemler effect, blue
galaxies are expected to be located within kpc from the center of
the cluster. However this prediction is inconsistent with the recent
observations. The increase of external pressure when the galaxy plows into the
intracluster medium does not change star formation rate of a disk galaxy
significantly. The velocity perturbation induced by a single high-speed
encounter between galaxies is too small to affect star formation rate of a disk
galaxy, while successive high-speed encounters (galaxy harassment) trigger star
formation activity because of the accumulation of gas in the galaxy center.
Therefore, the galaxy harassment remains as the candidate for a mechanism of
the Butcher-Oemler effect.Comment: 12 pages, 13 figures. To be published in Ap
Putting Randomized Compiler Testing into Production (Artifact)
This artifact accompanies our experience report for our compiler testing technology transfer project: taking the GraphicsFuzz research project on randomized metamorphic testing of graphics shader compilers, and building the necessary tooling around it to provide a highly automated process for improving the Khronos Vulkan Conformance Test Suite (CTS) with test cases that expose fuzzer-found compiler bugs, or that plug gaps in test coverage. The artifact consists of two Dockerfiles and associated files that can be used to build two Docker containers. The containers include our main tool for performing fuzzing: gfauto. The containers allow the user to fuzz SwiftShader, a software Vulkan implementation, finding 4 bugs. The user will also perform some line coverage analysis of SwiftShader using our tools to synthesize a small test that increases line coverage. Ubuntu, gfauto, SwiftShader, and other dependencies inside the Docker containers are fixed at specific versions, and all random seeds are set to specific values. Thus, all examples should reproduce faithfully on any machine
- …