Generation of density inhomogeneities by magnetohydrodynamic waves in two dimensions

Using two dimensional simulations, we study the formation of structures with a high-density contrast by magnetohydrodynamic waves in regions in which the ratio of thermal to magnetic pressure is small. The initial state is a uniform background perturbed by fast-mode wave. Our most significant result is that dense structures persist for far longer in a two-dimensional simulation than in the one-dimensional case. Once formed, these structures persist as long as the fast-mode amplitude remains high.Comment: 6 pages, 7 figures, accepted by MNRA

Mass transfer characteristics in structured packing for CO2 emission reduction processes

Acid gas treating and CO2 capture from flue gas by absorption have gained wide importance over the past few decades. With the implementation of more stringent environmental regulations and the awareness of the greenhouse effect, the need for efficient removal of acid gases such as CO2 (carbon dioxide) has increased significantly. Therefore, additional effort for research in this field is inevitable. For flue gas processes the ratio of absorption solvent to gas throughput is very different compared to acid gas treating processes owing to the atmospheric pressures and the dilution effect of combustion air. Moreover, in flue gas applications pressure drop is a very important process parameter. Packing types are required that allow for low pressure drop in combination with high interfacial areas at low liquid loading per square meter. The determination of interfacial areas in gas-liquid contactors by means of the chemical method (Danckwerts, P. V. Gas-liquid reactions; McGraw-Hill: London, 1970) has been very frequently applied. Unfortunately, many of the model systems proposed in the literature are reversible and therefore this condition possibly is not met. Versteeg et al. (Versteeg, G. F.; Kuipers, J. A. M.; Beckum, F. P. H.; van Swaaij, W. P. M. Chem. Eng. Sci. 1989, 44, 2292) have demonstrated that for reversible reactions the conditions for the determination of the interfacial area by means of the chemical method are much more severe. In a study by Raynal et al. (Raynal, L.; Ballaguet, J. P.; Berrere-Tricca, C. Chem. Eng. Sci. 2004, 59, 5395), it has been shown that there is a dependency of the interfacial area on the packing height. Unfortunately, most model systems used, e.g., CO2-caustic soda (as used by Raynal et al.), are much more complex and consist of (a set of) reversible reaction(s). The natures of these systems make the conditions at which the interfacial area can be determined much more severe and put more limitations on the process conditions and experimental equipment than a priori can be expected. Therefore, an extended absorption model is required to determine the conditions at which the interfacial area can be measured without detailed knowledge of the values of the liquid-side mass transfer coefficient, k1, beforehand.

Evolving grain-size distributions embedded in gas flows

We present a numerical approach for accurately evolving a dust grain-size distribution undergoing number-conserving (such as sputtering) and/or mass-conserving (such as shattering) processes. As typically observed interstellar dust distributions follow a power-law, our method adopts a power-law discretisation and uses both the grain mass and number densities in each bin to determine the power-law parameters. This power-law method is complementary to piecewise-constant and linear methods in the literature. We find that the power-law method surpasses the other two approaches, especially for small bin numbers. In the sputtering tests the relative error in the total grain mass remains below 0.01% independent of the number of bins N, while the other methods only achieve this for N > 50 or higher. Likewise, the shattering test shows that the method also produces small relative errors in the total grain numbers while conserving mass. Not only does the power-law method conserve the global distribution properties, it also preserves the inter-bin characteristics so that the shape of the distribution is recovered to a high degree. This does not always happen for the constant and linear methods, especially not for small bin numbers. Implementing the power-law method in a hydrodynamical code thus minimises the numerical cost whilst maintaining high accuracy. The method is not limited to dust grain distributions, but can also be applied to the evolution of any distribution function, such as a cosmic-ray distribution affected by synchrotron radiation or inverse-Compton scattering

Shock-triggered formation of magnetically-dominated clouds

To understand the formation of a magnetically dominated molecular cloud out of an atomic cloud, we follow the dynamical evolution of the cloud with a time-dependent axisymmetric magnetohydrodynamic code. A thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium interacts with the cloud. As a fast-mode shock propagates through the cloud, the gas behind it becomes thermally unstable. The $\beta$ value of the gas also becomes much smaller than the initial value of order unity. These conditions are ideal for magnetohydrodynamic waves to produce high-density clumps embedded in a rarefied warm medium. A slow-mode shock follows the fast-mode shock. Behind this shock a dense shell forms, which subsequently fragments. This is a primary region for the formation of massive stars. Our simulations show that only weak and moderate-strength shocks can form cold clouds which have properties typical of giant molecular clouds.Comment: 7 pages, 6 figures, accepted by Astronomy and Astrophysic

Strong coupling of magnons in a YIG sphere to photons in a planar superconducting resonator in the quantum limit

We report measurements of a superconducting coplanar waveguide resonator (CPWR) coupled to a sphere of yttrium-iron garnet. The non-uniform CPWR field allows us to excite various magnon modes in the sphere. Mode frequencies and relative coupling strengths are consistent with theory. Strong coupling is observed to several modes even with, on average, less than one excitation present in the CPWR. The time response to square pulses shows oscillations at the mode splitting frequency. These results indicate the feasibility of combining magnonic and planar superconducting quantum devices.Comment: 5 pages, 4 figure

Kiloparsec-Scale Simulations of star formations in disk galaxies. III. Structure and dynamics of filaments and clumps in giant molecular clouds

We present hydrodynamic simulations of self-gravitating dense gas in a galactic disk, exploring scales ranging from 1 kpc down to ~0.1 pc. Our primary goal is to understand how dense filaments form in giant molecular clouds (GMCs). These structures, often observed as infrared dark clouds (IRDCs) in the Galactic plane, are thought to be the precursors to massive stars and star clusters, so their formation may be the rate-limiting step controlling global star formation rates in galactic systems as described by the Kennicutt–Schmidt relation. Our study follows on from Van Loo et al., which carried out simulations to 0.5 pc resolution and examined global aspects of the formation of dense gas clumps and the resulting star formation rate. Here, using our higher resolution, we examine the detailed structural, kinematic, and dynamical properties of dense filaments and clumps, including mass surface density (Σ) probability distribution functions, filament mass per unit length and its dispersion, lateral Σ profiles, filament fragmentation, filament velocity gradients and infall, and degree of filament and clump virialization. Where possible, these properties are compared to observations of IRDCs. By many metrics, especially too large mass fractions of high ${\Sigma }\gt 1\;{\rm g}\;{\rm c}{{{\rm m}}^{-2}}$ material, too high mass per unit length dispersion due to dense clump formation, too high velocity gradients, and too high velocity dispersion for a given mass per unit length, the simulated filaments differ from observed IRDCs. We thus conclude that IRDCs do not form from global fast collapse of GMCs. Rather, we expect that IRDC formation and collapse are slowed significantly by the influence of dynamically important magnetic fields, which may thus play a crucial role in regulating galactic star formation rates

GMC Collisions as Triggers of Star Formation. I. Parameter Space Exploration with 2D Simulations

We utilize magnetohydrodynamic (MHD) simulations to develop a numerical model for GMC-GMC collisions between nearly magnetically critical clouds. The goal is to determine if, and under what circumstances, cloud collisions can cause pre-existing magnetically subcritical clumps to become supercritical and undergo gravitational collapse. We first develop and implement new photodissociation region (PDR) based heating and cooling functions that span the atomic to molecular transition, creating a multiphase ISM and allowing modeling of non-equilibrium temperature structures. Then in 2D and with ideal MHD, we explore a wide parameter space of magnetic field strength, magnetic field geometry, collision velocity, and impact parameter, and compare isolated versus colliding clouds. We find factors of ~2-3 increase in mean clump density from typical collisions, with strong dependence on collision velocity and magnetic field strength, but ultimately limited by flux-freezing in 2D geometries. For geometries enabling flow along magnetic field lines, greater degrees of collapse are seen. We discuss observational diagnostics of cloud collisions, focussing on 13CO(J=2-1), 13CO(J=3-2), and 12CO(J=8-7) integrated intensity maps and spectra, which we synthesize from our simulation outputs. We find the ratio of J=8-7 to lower-J emission is a powerful diagnostic probe of GMC collisions

Non-thermal radio emission from O-type stars III. Is Cyg OB2 No. 9 a wind-colliding binary?

The star Cyg OB2 No. 9 is a well-known non-thermal radio emitter. Recent theoretical work suggests that all such O-stars should be in a binary or a multiple system. However, there is no spectroscopic evidence of a binary component. Re-analysis of radio observations from the VLA of this system over 25 years has revealed that the non-thermal emission varies with a period of 2.35+-0.02 yr. This is interpreted as a strong suggestion of a binary system, with the non-thermal emission arising in a wind-collision region. We derived some preliminary orbital parameters for this putative binary and revised the mass-loss rate of the primary star downward from previous estimates.Comment: 13 pages, 5 figures, includes online data, accepted by A&

Can single O stars produce non-thermal radio emission?

We present a model for the non-thermal radio emission from presumably single O stars, in terms of synchrotron emission from relativistic electrons accelerated in wind-embedded shocks. These shocks are associated with an unstable, chaotic wind. The main improvement with respect to earlier models is the inclusion of the radial dependence of the shock velocity jump and compression ratio, based on 1D hydrodynamical simulations. The decrease of the velocity jump and the compression ratio as a function of radius produces a rapidly decreasing synchrotron emissivity. This effectively prohibits the models from reproducing the spectral shape of the observed non-thermal radio emission. We investigate a number of "escape routes" by which the hydrodynamical predictions might be reconciled with the radio observations. Although these escape routes reproduce the observed spectral shape, none of these escape routes are physically plausible. In particular, re-acceleration by feeding an electron distribution through a number of shocks, is in contradiction with current hydrodynamical simulations. These hydrodynamical simulations have their limitations, most notably the use of 1D. At present, it is not feasible to perform 2D simulations of the wind out to the distances required for synchrotron-emission models. Based on the current hydrodynamic models, we suspect that the observed non-thermal radio emission from O stars cannot be explained by wind-embedded shocks associated with the instability of the line-driving mechanism. The most likely alternative mechanism is synchrotron emission from colliding winds. That would imply that all O stars with non-thermal radio emission should be members of binary or multiple systems.Comment: 10 pages, 8 figures, accepted for publication by A&