1. Introduction
It has been recognized for a long time that interstellar dust can be elevated outside galactic disks into their halos and circumgalactic gas (CGM), and farther away into extragalactic space either by radiation pressure or by supernova-driven winds (e.g., [
1,
2,
3,
4,
5]). An interstellar magnetic field seems to inhibit radiation-driven dust outflows as mentioned in [
6]; however, under certain conditions, large-scale magneto-gravitational (Parker) instability can enhance them [
7]. Observational evidence of dust’s presence in the galactic halo at heights
–5 kpc were reported first in the galaxy NGC 891 [
8]. Later on, clear traces of dust at 10–12 kpc above the disk plane were confirmed in [
9,
10]. Dusty halos extending into the circumgalactic environment at much larger scales,
60–200 kpc, have been found in two nearby galaxies, NGC 2835 and NGC 3521 [
11]. This result was confirmed in later studies [
12,
13,
14].
More recently, the presence of dust on megaparsec scales,
Mpc, with a considerable dust mass,
, and a dust-to-gas mass fraction of the order of the local ISM, extending from parent galaxies to the intergalactic medium, has been demonstrated in [
15,
16]. It is clear that the transport of unshielded dust particles over such long distances in a hot environment in a CGM is accompanied by the destruction of the grains via thermal or kinetic sputtering [
17]. Therefore, it is usually assumed that such transport is carried out via the expelling of dense and cold gas clouds and clumps containing dust (for a more recent discussion, see review [
18]).
Estimates of the dust mass fraction in galaxy clusters (GCs) are less certain. Measurements of extinction and reddening reveal the presence of dust up to radii ∼
Mpc with the total dust mass in the intracluster medium (ICM) of ∼
, i.e, the dust-to-gas mass ratio is around 0.4% of the local ISM, assuming the ICM mass ∼
[
19,
20]. On the other hand, estimates of the dust-to-gas mass ratio inferred from the IR observation of GCs only allow an upper limit to be put, which vary for different GC samples from <
% [
21] to <
% [
22,
23] of the local ISM. However, in both cases, the dust-to-gas mass ratio is considerably lower than in the local ISM, indicating a much stronger degree of dust destruction than is found in the galactic and circumgalactic haloes. The reasons for such differences are not very obvious because the distances between galaxies in a GC of the order of ∼
Mpc are smaller than the sizes of the CGM regions in [
15,
16]; here,
is the number of galaxies in the cluster. Numerical analyses of these reasons are important for understanding the life cycle of dust on intergalactic scales.
It is known that naked (unshielded) dust particles are destroyed in a rather dense (
cm
−3) and hot (
K) environment in a CGM and in an ICM gas while traveling on timescales close to the crossing times typical for a CGM and ICM [
18,
24]. Usually, dust expelled out of galaxies is supposed to be confined into denser gas clouds. Under such conditions, dust survives for lengths of time comparable to the clouds’ lifetimes that suffered destruction by hydrodynamical instabilities while moving outwards. Once clouds are destructed, dust particles of 0.1 μm can survive no longer than ∼100 Myr. This is, therefore, our primary motivation to consider hydrodynamical destruction of dusty clouds moving through a hot intracluster gas.
The problem we are pursuing here is equivalent to a gas cloud interacting with a wind, with an additional process connected with dust destruction when dust particles come into direct contact with a hot and dense ambient gas. This requires an adequate numerical procedure which would allow tracking an immediate environmental conditions in a gas volume carrying a dust particle.
To construct such models, one needs to track the cloud material as accurately as possible in order to diminish numerical diffusion of a dust particle into a hot environment. Complications come from different dynamics of gas and dust along with an unavoidable influence of numerical mixing of cold cloud and hot ambient gas while the cloud moves.
Our goal here is to minimize these effects. In several studies, the dynamics of dust has been described by using Lagrangian (trace) particles [
25]. This approach cannot be considered as adequate in general due to dust being inertial and suffering drag force by a gas [
26]; its mass can become smaller due to destruction during the evolution. In this paper, we can distinguish the dynamics of Lagrangian and dust particles, and compare them. In
Section 2, we describe our model. In
Section 3, we consider the tracking methods. In
Section 4, we describe the dynamics of the ‘cloud-wind’ interaction. In
Section 5, we summarize our results.
3. Cloud Tracking
Figure 1 shows the mass fraction of the colored gas for different levels of the tracer concentration
C:
(left panel). As defined, the concentration is
inside the cloud at the initial moment. During evolution, the tracer enclosed initially inside the cloud mixes with ambient noncolored gas due to numerical diffusion. This makes such gas belong to the cloud. Till
, the mass of the colored gas with level
remains close to the initial mass of the cloud. Whereas for
the mass of the colored gas is increased by less than 5% to the initial mass of the cloud, for a lower value of
, it is about 10% at this time. During further evolution, the ratio
for
stays around the same until
and begins to drop gradually until ∼
, and then even more substantially. The ratio for
varies weakly and increases not more than 20% till the final time
, while the ratio for the lowest level
increases monotonously at a factor of more than 2.6 at this time.
The increase in numerical resolution leads to better tracking of cloud material. The right panel in
Figure 1 presents the mass fraction of the colored gas for the runs with different spatial resolutions (the ratio of the cell size to the cloud radius) for
. As expected, the ratio becomes closer to being the same as that for a higher resolution, which requires larger computational resources.
Using trace particles, we follow the gas in a cloud. Lagrangian particles (lps) can trace the gas because their velocities are kept equal to that of the gas in a cell, where a particle is located, at each instance. It is obvious that the number of particles enclosed in a cell should be large enough for reasonable tracking of the gas in the cloud. Usually, this number is about ten particles per cell (e.g., [
42]). In our task, a cloud can be destructed and its material can be spread over a larger volume. Thus, we use more particles in our runs.
Each trace particle is ordered to trace a mass fraction of a cloud equal to
(where
is the total number of the particles). Thus, the total mass of these particles remains constant by definition, which can be seen in
Figure 1.
Further, we consider how such particles trace the thermal phase of a gas belonging initially to a cloud. We assume that a particle carries the information of a part of the cloud material in which it is located.
One should note that the increase in the colored gas mass can be considered as mixing of the ‘cloud’ material with the ambient gas. In this sense, such an increase is the result not only of partial mixing with the intercloud medium but also numerical diffusion. A cloud crushed by shock or wind passes through a mixing phase (e.g., [
48]); its duration is determined by the Kelvin–Helmholtz instability. However, cooling can start to play a role in the dynamics of a gas behind the shock front. Its influence depends on the metallicity, shock velocity, and external heating sources. Under such conditions, the cooling can become efficient and the development of the Kelvin–Helmholtz instability is suppressed. Thus, it is interesting to exclude these mixing effects as much as possible. Lagrangian particles provide sharper gradients and more intermittent distribution of density fields (e.g., [
45,
46]). Based on this, we apply trace particles for tracking cloud material. Another aim of using such particles is to compare their dynamics with dust particles.
4. Dynamics
Let us follow the evolution of the diffuse cloud that interacted with a hot postshock inflow. As mentioned above, after the interaction begins, the shock front moves into the cloud with the velocity of factor
lower. This corresponds to 30 km s
−1 for the model considered here. The gas behind the front in the cloud is heated to
K. The cooling time for the metallicity 0.1
is about 30 kyr, which is significantly shorter than the crushing time
Myr. Such an efficient cooling prevents the stripping of outer layers of the cloud. Indeed, the stripping operates under the action of the Kelvin–Helmholtz instability with the characteristic time
, where
with the maximum growth rate is comparable to the cloud radius
as noted by Klein et al. [
29]. It is seen that the time needed for developing the longest scale of the Kelvin–Helmholtz instability coincides with the crushing time. Therefore, one can assume that the cloud compressed by the shock front and the postshock flow does not experience a considerable mass loss due to stripping at timescales shorter than the crushing time
t tcc. One can see that the Kelvin–Helmholtz instability is suppressed within several crushing times under the conditions considered here. However, it can become significant again; certainly in the case where the compressed lump does not come to be gravitationally unstable. A short discussion about the importance of cooling can be found in
Appendix A, where the interaction of an adiabatic cloud with a shockwave is described.
Figure 2 presents the distributions of various values for gas and dust at
Myr
. Here, we describe in detail the slices for the fiducial model with
(the upper row of panels). Before we proceed to this description, we note that the morphology of the cloud material depends slightly on the spatial resolution. In the lower panels, one can find the maps for the model with the lowest resolution
. It is worth noting that the morphology of the cloud as a whole stays the same. The increase in resolution leads to more detailed structures. In the previous studies, the minimum ratio of
needed for converged mass growth and cloud entrainment was established [
30]. We use three times higher resolution for our fiducial model. In addition, the destruction of an adiabatic cloud is presented in
Appendix A. We carried out several runs for different spatial resolutions and numbers of trace/dust particles and found a good convergence.
One can assume that the extended structures on the edges of the compressed cloud seen in the gas density field for
Myr
(first panel of
Figure 2) are connected with the stripping. Similar structures can be found in the gas temperature field (second panel). However, the estimates presented above are not in favor of stripping. We argue that inflowing gas interacting with the cloud flows around it and forms such structures. During this process, the ambient gas mixes numerically with the cloud represented here by colored gas (third panel). As a result, up to
Myr
the cloud gas traced by color forms a cometary tail extended up to ∼
kpc or about ∼5 times larger than the initial radius of the cloud. Note that to follow the evolution of such structures till the final time
, we extend the grid size up to
kpc, which corresponds to the number of cells equaling 3000 in the direction of the flow. We present the distributions for the earliest time, and mainly, the cloud material extends highly along the flow propagation. Moreover, this time is of importance for understanding the relation between the cooling and Kelvin–Helmholtz timescales.
At the same time, the cloud material tracked by Lagrangian particles does not demonstrate any cometary tail (fourth panel): the cloud ‘represented’ by Lagrangian particles is compressed to a pancake-like structure. Their velocity field follows the gas by definition. Initially, they are attributed to the gas belonging to the cloud, and there are no trace particles outside it. Thus, under the suppression of the Kelvin–Helmholtz instability due to efficient cooling, the gas associated with the cloud is compressed into a pancake-like structure.
Dust particles show different dynamics; they represent a separate inertial ‘fluid’ interacting with gas by drag force, which is responsible for the gaining velocity of dust particles. The dynamics of the dust particles is determined by the stopping time
[
49] (
is the density of the grain material,
a is the grain size,
is the gas density,
is the relative velocity between gas and dust). The time needed to accelerate dust particles to the velocity of the postshock inflow is around 100 kyr. Within this timescale, dust swept out from the outer layers of the cloud gains the velocity close to the postshock inflow (fifth panel of
Figure 2). This time appears to be longer than the cooling time, during which the gas inside the cloud becomes colder and denser, which leads to an increase in the stopping rate of dust
, and as a consequence, the amount of swept-out dust is apparently low.
Figure 3 (left panel) presents the mass of a gas of different thermal phases (the labels correspond to the logarithm of gas temperature). These phases are tracked by colored fluid for
(thin solid lines),
(thin dashed lines), and Lagrangian particles (thick lines). The gas belonging to the cloud has a temperature within
–
K during the first period of ∼
–
. After that, the major part of the gaseous mass cools down and becomes enclosed in the cold phase with
K. This process is completed at
t ~ 3
tcc. It is worth noting that two methods for tracing cloud material give very close values for the mass fractions enclosed in the cold (
K) and diffuse (
K) phases. In addition, the mass fractions of the cloud material located in these phases by colored fluid for two levels
C > 0.1 and
C > 0.3 almost coincide with each other. A similar behavior can be seen for the averaged (mass-weighted) velocity of both cold and diffuse gas (right panel of
Figure 3). During the first ~0.5
tcc the velocities of these phases increase almost linearly up to ∼30 km/s and then remain around this value.
A significant discrepancy between values determined by using two methods for tracking cloud material is found for the hot ( K) phase: it reaches more than one order of magnitude for both mass fraction (compare red lines at left panel) and velocity (see right panel). According to the estimates of the cooling and Kelvin–Helmholtz timescales made above, the stripping is not efficient in the model considered here. Within , the mass fraction of the hot gas remains around 0.1 for and drops several times for . Apparently, the colored gas in the hot phase originates from numerical diffusion. This strong dependence on the level of color is in favor of numerical mixing between the gas of the cloud represented by colored fluid and the gas flowing around the cloud.
Lagrangian particles keep the total mass of the cloud exactly by definition (
Figure 1). One should note that, when we interpolate fluid velocity at the Lagrangian particle location, we use a shape function (here the “Cloud-In-Cell” function is chosen), so that each particle can give a nonzero contribution not only to the computational zone hosting the particle, but also to its left and right neighbors (see, e.g., [
37,
38,
39]). This can be considered as numerical mixing as well. However, the contributions are added only to the neighboring cells and are not spread to other ones. Under the conditions considered here, the temperature of the gas belonging to the cloud is below
K due to efficient cooling, and hence, a low amount of gas can be found in the hot phase due to the interpolating procedure: this value barely exceeds 0.01. Note that this mass fraction of the hot phase tracked by Lagrangian particles is much lower than that traced by colored fluid. The velocity of such hot gas is significantly lower than the unperturbed flow because in the interpolation only, the contributions from the neighboring cells are included. This picture takes place until the destruction of the cloud is inefficient.
At
–
, the mass fraction of the cold phase reaches its maximum. Afterwards, it begins to decrease apparently due to the destruction of the cold phase, which leads to the growth of the fraction in the diffuse phase. This comes from that the Kelvin–Helmholtz instability, which begins to play a role again. At this time, the largest size of the dense lump
1 decreases by ∼3 times, and the overdensity
reaches ∼100. Thus, the instability can develop at
t 3
tcc. One can note a remarkable discrepancy in both the cold and diffuse fractions after
t (2.5–3)
tcc between the two methods of tracking clouds: it reaches a factor of 1.5–2. It should be stressed that the fraction of the cold gas in the case of tracking by colored fluid decreases faster compared to that tracked by Lagrangian particles.
For the cloud tracked by colored fluid, the cold fraction begins to drop just after . This is accompanied by the growth of the diffuse fraction of the cloud from ∼30% at ∼ to more than 70% at . As the considerable fraction of the cloud mass transforms into the diffuse phase, it can give a wrong impression that the cloud is destroyed.
If the gas of the cloud is tracked by particles, one can see that the cloud keeps its identity longer. Only after t 5tcc does the cold fraction decrease with simultaneous growth of both the diffuse and hot fractions. The cold fraction is higher than that determined by using the tracking by colored fluid.
At
, almost 99% (by mass) of the gas belonging initially to the cloud moves with a velocity of ∼30 km s
−1; that is around 10% of the unperturbed postshock wind velocity. This part of the cloud is located in the cold and diffuse phases (
Figure 3). This result reflects a low entrainment efficiency when approaching this time. This conclusion is almost independent of the method of tracking cloud.
The mass fractions of dust particles contained within a gas in both the diffuse and cold phases evolve closely to these Lagrangian particles (compare left panels of
Figure 3 and
Figure 4). At the same time, in the hot phase, the dust particles demonstrate different dynamics (
Figure 4). This is because they are swept off by the shockwave from the sides of the cloud into the hot gas, where they become involved in the hot postshock flow (fifth panel of
Figure 2). One can see that the velocity of the dust associated with the hot phase rapidly increases up to ∼200 km/s during ∼
–
, and after
, reaches ∼250 km/s, which is ∼
of the wind velocity. Under such conditions, dust particles are efficiently destroyed by sputtering. This determines the mass loss of dust: the decrease in the total dust mass fraction is about 10% up to
(as shown by the right axis of the left panel in
Figure 4). During the later evolution, the mass-averaged dust velocity decreases due to destruction of the cloud and, as a consequence, replenishes the hot phase with still undestroyed dust particles stripped from the cloud and with lower velocities.
Above we consider the values averaged over a few coarse thermal bins. In order to take a look at more detailed distributions, we choose two time moments: 5 and
.
Figure 5 present the mass fraction distribution functions of temperature, number density, and velocity of a gas in a cloud for these moments. The two peaks clearly seen in all distributions correspond to the cold dense clump and to the hot rarefied flow, respectively. The fraction of a gas in the intermediate range is smaller than those enclosed in both the cold and hot phases, but it increases with time. The distributions of the cloud fractions tracked by colored fluid demonstrate a considerable excess over those tracked by Lagrangian particles; this becomes more remarkable for the diffuse gas with
T 5 × 10
3 K and
cm
−3 moving with
v 50 km/s. This discrepancy increases at 10
tcc in comparison with 5
tcc and will continue to grow during further evolution due to ongoing cloud destruction.
At both time moments, the mass distributions for Lagrangian and dust particles demonstrate similar behaviors for the cold phase only. This means that Lagrangian particles provide an efficient tracing of dust associated with only cloud material in the cold phase. The mass fraction of dust associated with hot gas is less than several percent. Such dust particles are of importance due to their efficient sputtering: at first, they can enrich the ambient gas with metals and, secondly, by becoming smaller in size they emit in a short infrared range. In the range
–
K, there is a significant gap in the distribution for the gas in which dust particles are enclosed at
(blue line, upper left panel of
Figure 5). More exactly, the mass fraction in this range is less than 0.01%. Due to sweeping, this fraction increases more than 100 times up to the time
(lower left panel). The same can be found for the intermediate values of the number density of the gas where the dust particles are located: one can see that the mass fraction of the dust associated with the gas of
–
cm
−3 is also less than 0.01%, and at
, the fraction increases by 30–100 times (see the two middle panels).
At , one can see the considerable differences in the velocity distributions for v 50 km/s, which become more significant for v 100 km/s (upper right panel). During further evolution, the discrepancy between the Lagrangian and dust particles decreases, whereas the distribution of the gas velocity tracked by colored fluid keeps the significant difference from other distributions.
5. Conclusions
In this work, we analyzed how a diffuse radiatively cooled cloud () interacts with the hot 300 km/c postshock wind. We utilized two methods to track the cloud material: the first is based on a ‘colored’ fluid approach, and the second one uses the Lagrangian (trace) particles. We followed the dynamics of the dust grains interacting with the ambient gas by the drag force up to their destruction after they are swept out of the ‘parent’ cloud by the hot flow. We find that the radiative losses can play a significant role in the dynamics of a gas behind the shock front. We find that under the conditions considered in this work, the cooling becomes efficient during the earliest times and the development of the Kelvin–Helmholtz instability is suppressed. In particular, the destruction of a warm cloud of size 100 pc with metallicity higher than 0.1 interacting with the hot flow of the ICM is expected to be delayed up to several crushing times for such a cloud.
For a gas tracked by colored fluid, we show that the total mass of a gas associated with the cloud strongly depends on the level of tracer concentration higher than that of the gas that is believed to belong to the cloud: for a low level, , the discrepancy reaches a factor of several times at ten crushing times of the cloud (); the increase in the level results in a smaller difference but leads to unpredictable variations in the total mass estimated by using colored fluid during the evolution. Certainly, the increase in spatial resolution a gives lower difference. Such an increase in the colored gas mass can be represented by mixing of the ‘cloud’ material with the ambient gas. This process includes unavoidable effects of numerical diffusion.
We find that when traced by Lagrangian particles, the cloud material has sharper gradients and more intermittent density field distribution than when traced by the colored fluid. However, it provides an alternative approach to trace the cloud material. There is a special interest in this method because trace particles are used for modeling dust [
25].
We found that at almost 99% (by mass) of the gas, initially attributed to the cloud, has at most 10% of the postshock wind velocity. We interpret this as being a result of the weak entrainment efficiency of the cloud. We find that, irrespective of the approach used to trace the cloud material, such gas is enclosed in the cold ( K) and diffuse ( K) phases. However, the gas fractions in the cold and diffuse phases differ depending on the method: the colored fluid approach results in about two times lower mass being found in the cold phase.
Our results suggest that the mass fraction of the cloud found in the hot phase ( K) depends on the tracer concentration level. This fraction appeared to be significantly higher for the fiducial level of than that determined in the case when the Lagrangian particles are used to trace the gas of the cloud.
Using tracking methods, one can infer the cloud entrainment efficiency, and the distribution of cloud material between various thermal phases. The two methods considered in this work allow the detection of differences in the dynamics of the cloud material. In particular, they reveal variations in the gas distribution in each of the assumed thermal phases, as well as distinctive mixing of the cloud material with the ambient gas. We note that these methods depend on spatial resolution: the accuracy with which the cloud material is traced increases for a smaller cell size and/or larger number of trace particles.
Dust particles can be swept out from their ‘parent’ cloud material; further, their dynamics in diffuse and hot rarefied gas differs from that of Lagrangian particles. In dense parts of a cloud, the velocity of the shockwave is significantly lower, so the dust dynamics as an ensemble of particles is close to that of Lagrangian ones. However, the difference becomes more prominent in the cloud parts mixed efficiently due to an irregular velocity field. In addition, the dynamics of dust particles is changed by their destruction, which becomes significant in the hot, high-velocity gas. Therefore, the dust particles cannot be correctly described as Lagrangian ones in a wide range of conditions that appear in the process of the interaction between a cloud and a hot wind (or strong shockwave).
The destruction of dusty clouds interacting with a hot flow is closely related to the accurate interpretation of the extinction measurements in the ICM. The large scatter of the extinction values can be connected not only with observational uncertainties (e.g., see [
18,
50] for a review), but also with erroneous conclusions about the destruction of dusty clouds moving in the hot gas.