We present results of full 3D hydrodynamical and radiative transfer simulations of the colliding stellar winds in the massive binary system Eta Carinae. We accomplish this by applying the SimpleX algorithm for 3D radiative transfer on an unstructured Voronoi-Delaunay grid to recent 3D smoothed particle hydrodynamics (SPH) simulations of the binary colliding winds. We use SimpleX to obtain detailed ionization fractions of hydrogen and helium, in 3D, at the resolution of the original SPH simulations. We investigate several computational domain sizes and Luminous Blue Variable primary star mass-loss rates. We furthermore present new methods of visualizing and interacting with output from complex 3D numerical simulations, including 3D interactive graphics and 3D printing. While we initially focus on Eta Car, the methods employed can be applied to numerous other colliding wind (WR 140, WR 137, WR 19) and dusty 'pinwheel' (WR 104, WR 98a) binary systems. Coupled with 3D hydrodynamical simulations, SimpleX simulations have the potential to help determine the regions where various observed time-variable emission and absorption lines form in these unique objects.

Eta Carinae is an ideal astrophysical laboratory for studying massive binary interactions and evolution, and stellar wind-wind collisions. Recent three-dimensional (3D) simulations set the stage for understanding the highly complex 3D flows in $\eta$ Car. Observations of different broad high- and low-ionization forbidden emission lines provide an excellent tool to constrain the orientation of the system, the primary's mass-loss rate, and the ionizing flux of the hot secondary. In this work we present the first steps towards generating synthetic observations to compare with available and future HST/STIS data. We present initial results from full 3D radiative transfer simulations of the interacting winds in $\eta$ Car. We use the SimpleX algorithm to post-process the output from 3D SPH simulations and obtain the ionization fractions of hydrogen and helium assuming three different mass-loss rates for the primary star. The resultant ionization maps of both species constrain the regions where the observed forbidden emission lines can form. Including collisional ionization is necessary to achieve a better description of the ionization states, especially in the areas shielded from the secondary's radiation. We find that reducing the primary's mass-loss rate increases the volume of ionized gas, creating larger areas where the forbidden emission lines can form. We conclude that post processing 3D SPH data with SimpleX is a viable tool to create ionization maps for $\eta$ Car.

(Abridged) Interacting galaxies are well-known for their high star formation rates and rich star cluster populations, but the rapidly changing tidal field can also efficiently destroy clusters. We use numerical simulations of merging disc galaxies to investigate which mechanism dominates. The simulations include a model for the formation and dynamical disruption of the entire star cluster population. We find that the dynamical heating of clusters by tidal shocks is about an order of magnitude higher in interacting galaxies than in isolated galaxies. This is driven by the increased gas density, and is sufficient to destroy star clusters at a higher rate than new clusters are formed: the total number of clusters in the merger remnant is 2-50% of the amount in the progenitor discs, with low-mass clusters being disrupted preferentially. By adopting observationally motivated selection criteria, we find that the observed surplus of star clusters in nearby merging galaxies is caused by the bias to detect young, massive clusters. We provide a general expression for the survival fraction of clusters, which increases with the gas depletion time-scale. Due to the preferential disruption of low-mass clusters, the mass distribution of the surviving star clusters in a merger remnant develops a peak at a mass of about 10^3 Msun, which evolves to higher masses at a rate of 0.3-0.4 dex per Gyr. The peak mass initially depends weakly on the galactocentric radius, but this correlation disappears as the system ages. We discuss the similarities between the cluster populations of the simulated merger remnants and (young) globular cluster systems. Our results suggest that the combination of cluster formation and destruction should be widespread in the dense star-forming environments at high redshifts, which could provide a natural origin to present-day globular cluster systems.

(Abridged) The formation and evolution of star cluster populations are related to the galactic environment. Cluster formation is governed by processes acting on galactic scales, and star cluster disruption is driven by the tidal field. In this paper, we present a self-consistent model for the formation and evolution of star cluster populations, for which we combine an N-body/SPH galaxy evolution code with semi-analytic models for star cluster evolution. The model includes star formation, feedback, stellar evolution, and star cluster disruption by two-body relaxation and tidal shocks. We apply the model by simulating a suite of 9 isolated disc galaxies and 24 galaxy mergers. The evolutionary histories of individual clusters in these simulations are discussed to illustrate how the environment of clusters changes in time and space. The resulting variability of the disruption rate with time and space affects the properties of star cluster populations. The combined effect of clusters escaping their dense formation sites (`cluster migration') and the preferential disruption of clusters residing in dense environments (`natural selection') implies that the mean disruption rate of the population decreases with cluster age. This affects the slope of the cluster age distribution, which becomes a function of the star formation rate density. The evolutionary histories of clusters in a galaxy merger vary widely and determine which clusters survive the merger. This impacts the age distributions and the locations of the surviving clusters at all times during a merger. We conclude that accounting for the interplay between the formation, disruption, and orbital histories of clusters enables a more sophisticated interpretation of observed properties of cluster populations, thereby extending the role of cluster populations as tracers of galaxy evolution.

We present an improved version of the SimpleX method for radiative transfer on an unstructured Delaunay grid. The grid samples the medium through which photons are transported in an optimal way for fast radiative transfer calculations. In this paper we study the detailed working of SimpleX in test problems and show improvements over the previous version of the method. We have applied a direction conserving transport scheme that correctly transports photons in the optically thin regime, a regime where the original SimpleX algorithm lost its accuracy. In addition, a scheme of dynamic grid updates is employed to ensure correct photon transport when the optical depth changes during a simulation. For the application to large data sets, the method is parallellised for distributed memory machines using MPI. To test the new method, we have performed standard tests for cosmological radiative transfer. We show excellent correspondence with both the analytical solution (when available) and to the results of other codes compared to the former version of SimpleX, without sacrificing the benefits of the high computational speed of the method.

Context. Analytical and numerical analysis of the SimpleX radiative transfer algorithm, which features transport on a Delaunay triangulation. Aims. Verify whether the SimpleX radiative transfer algorithm conforms to mathematical expectations, to develop error analysis and present improvements upon earlier versions of the code. Methods. Voronoi-Delaunay tessellation, classical Markov theory. Results. Quantitative description of the error properties of the SimpleX method. Numerical validation of the method and verification of the analytical results. Improvements in accuracy and speed of the method. Conclusions. It is possible to transport particles such as photons in a physically correct manner with the SimpleX algorithm. This requires the use of weighting schemes or the modification of the point process underlying the transport grid. We have explored and applied several possibilities.

In this paper, we present a new method for the solution of those linear transport processes that may be described by a Master Equation, such as electron, neutron and photon transport, and more exotic variants thereof. We base our algorithm on a Markov process on a Voronoi-Delaunay grid, a nonperiodic lattice which is derived from a random point process that is chosen to optimally represent certain properties of the medium through which the transport occurs. Our grid is locally translation and rotation invariant in the mean. We illustrate our approach by means of a particular example, in which the expectation value of the length of a grid line corresponds to the local mean free path. In this example, the lattice is a direct representation of the `free path space' of the medium. Subsequently, transport is defined as simply moving particles from one node to the next, interactions taking place at each point. We derive the statistical properties of such lattices, describe the limiting behavior, and show how interactions are incorporated as global coefficients. Two elementary linear transport problems are discussed: that of free ballistic transport, and the transport of particles through a scattering medium. We also mention a combination of these two. We discuss the efficiency of our method, showing that it is much faster than most other methods, because the operation count does not scale with the number of sources. (ABRIDGED)

We present a new numerical approach that is able to solve the multi-dimensional radiative transfer equations in all opacity regimes on a Lagrangian, unstructured network of characteristics based on a stochastic point process. Our method reverses the limiting procedure used to derive the transfer equations, by going back to the original Markov process. Thus, we reduce this highly complex system of coupled differential equations to a simple one-dimensional random walk on a graph, which is shown to be computationally very efficient. Specifically, we use a Delaunay graph, which makes it possible to combine our scheme with a new smoothed particle hydrodynamics (SPH) variant proposed by Pelupessy et al.(2003). We show that the results of applying a two-dimensional implementation of our method with various suitable test cases agree with the analytical results, and we point out the advantages of using our method with inhomogeneous point distributions, showing examples in the progress. Hereafter, we present a supplement to our method, which can be useful in cases where the medium is optically very thin, and we conclude by stating some anticipated properties of this method in three dimensions, and announce future extensions.

We present N-body/SPH simulations of the evolution of an isolated dwarf galaxy including a detailed model for the ISM, star formation and stellar feedback. Depending on the strength of the feedback, the modelled dwarf galaxy shows periodic or quasi-periodic bursts of star formation of moderate strength. The period of the variations is related to the dynamical timescale, of the order of $1.5~10^8$ yr. We show that the results of these simulations are in good agreement with recent detailed observations of dwarf irregulars (dIrr) and that the peculiar kinematic and morphological properties of these objects,as revealed by high resolution HI studies, are fully reproduced. We discuss these results in the context of recent surveys of dwarf galaxies and point out that if the star formation pattern of our model galaxy is typical for dwarf irregulars this could explain the scatter of observed properties of dwarf galaxies. Specifically, we show that the time sampled distribution of the ratio between the instanteneous star formation rate (SFR) and the mean SFR is similar to that distribution in observed sample of dwarf galaxies.

At the end of their lives low mass stars such as our Sun lose most of their mass. The resulting planetary nebulae show a wide variety of shapes, from spherical to highly bipolar. According to the generalized interacting stellar winds model, these shapes are due to an interaction between a very fast tenuous outflow, and a denser environment left over from an earlier slow phase of mass loss. Previous analytical and numerical work shows that this mechanism can explain cylindrically symmetric nebulae very well. However, many circumstellar nebulae have a multipolar or point-symmetric shape. With two-dimensional calculations, Icke showed that these seemingly enigmatic forms can be easily reproduced by a two-wind model in which the confining disk is warped, as is expected to occur in irradiated disks. Here, we present the extension to fully three-dimensional adaptive mesh refinement simulations of such an interaction.

Previous analytical and numerical work shows that the generalized interacting stellar winds model can explain the observed bipolar shapes of planetary nebulae very well. However, many circumstellar nebulae have a multipolar or point-symmetric shape. With two-dimensional calculations, Icke showed that these seemingly enigmatic forms can be easily reproduced by a two-wind model in which the confining disk is warped, as is expected to occur in irradiated disks. In this contribution we present the extension to fully three-dimensional adaptive mesh refinement simulations of such an interaction.

We report the discovery of a double ring structure in the waist of the nebula surrounding eta Carinae. The rings are detected in the mid-IR dust continuum at wavelengths of 7.9, 11.9, 12.9 and 20 micronmeter. The dust in the rings has a temperature of about 300 K. The orientation of the rings is inclined with respect to the axis of the homunculus by either 37 or 58 degrees. The central star is not in the projected centre of the structure defined by the two rings. This geometry is reminiscent of that seen in SN1987A and some planetary nebulae. We discuss several possible origins for this remarkable geometry and its orientation.

We have combined time dependent hydrodynamics with a two-fluid model for dust driven AGB winds. Our calculations include self-consistent gas chemistry, grain formation and growth, and a new implementation of the viscous momentum transfer between grains and gas. This allows us to perform calculations in which no assumptions about the completeness of momentum coupling are made. We derive new expressions to treat time dependent and non-equilibrium drift in a hydro code. Using a stationary state calculation for IRC +10216 as initial model, the time dependent integration leads to a quasi-periodic mass loss in the case where dust drift is taken into account. The time scale of the variation is of the order of a few hundred years, which corresponds to the time scale needed to explain the shell structure of the envelope of IRC +10216 and other AGB and post-AGB stars, which has been a puzzle since its discovery. No such periodicity is observed in comparison models without drift between dust and gas.

We have developed an implementation for the momentum transfer force in numerical two fluid hydrodynamics. This form of the frictional coupling between gas and grains is consistent with the microscopic interactions between the two components. The coupling force gives rise to a drift velocity of the grains with respect to the gas. We apply this mechanism to the outflow of (Post-) AGB objects. Our numerical hydrodynamics code calculates self consistently the dynamics of these outflows, as well as the nucleation and growth of grains and equilibrium chemistry of the gas. Grain nucleation and growth are processes that depend strongly on the rate of gas-grain collisions. Hence, the drift velocity becomes an important variable. The tight connection between grain chemistry and drift causes the system to become extremely sensitive to small changes in almost any parameter. This may be a cause for deviation from (spherical) symmetry and structure.