Explicit finite volume modeling of aircraft anti-icing and de-icing

In this paper, a runback water and ice prediction model is extended to anti-icing and thermal de-icing situations. The resulting coupled equations that govern thin-film flow, ice accretion, and heat conduction in the multilayered system substrate-ice-water are solved using an explicit finite volume approach. The procedure is implemented in the three-dimensional icing code ICECREMO2, and both structured and unstructured grids can be considered. Numerical results are presented to compare the present code simulations to some data provided by other ice prediction codes and to show the capabilities of the present numerical tool.

The Application Of Radioisotopes In The Study Of Estuarine Sedimentary Processes

The spatial and temporal distributions of some radionuclides in effluents originating from the British Nuclear Fuels Ltd (BNFL) reprocessing plant at Windscale, which are released into the Irish Sea, have been studied in sediments at 16 sites in the salt marsh region near Newbiggin on the Esk estuary Cumbria, England. The concentration of non-conservative radionuclides in surface sediments of the area cannot be described by a single parameter, but there is a high correlation with organic C, Cu, Al and the Si : Al ratio with particle size. The preservation of the historical record of the BNFL effluents in the Esk sediments is dependent on the hydrology of the area, as it effects such processes as accretion, erosion and remixing. From the 106Ru and 210Po concentrations and the 137Cs : 134Cs ratio in the sediment profiles with depth, we have identified these processes. Sedimentation rates at sites of accretion vary between 0·5 and 3 cm year−1. However, at some sites they appear to be much higher, approximately 6 cm year−1 in the top 10 cm, but they are not consistent throughout the depth profiles. This may be a true reflection of variable accretion related to sediment type, or one which is influenced by surficial mixing. Some cores showed evidence of continuous accretion but no significant radioactivity was detected at depths below 35–40 cm, indicating an overall sedimentation rate of approximately 1·5 cm year−1 for the 25–30-year period since BNFL effluents first entered the Irish Sea.

Molecular Hydrogen Emission from Protoplanetary Disks

We have modeled self-consistently the density and temperature profiles of gas and dust in protoplanetary disks, taking into account irradiation from a central star. Making use of this physical structure, we have calculated the level populations of molecular hydrogen and the line emission from the disks. As a result, we can reproduce the observed strong line spectra of molecular hydrogen from protoplanetary disks, both in the ultraviolet (UV) and the near-infrared, but only if the central star has a strong UV excess radiation.

Molecular Distributions in the Inner Regions of Protostellar Disks

The distributions of molecules in the inner regions of a protostellar disk are presented. These were calculated using an uncoupled chemical/dynamical model, with a numerical integration of the vertical disk structure. A comparison between models with and without the effects of X-ray ionisation is made, and molecules are identified which are good tracers of the ionisation level in this part of the disk, notably CN and C_2H. In the region considered in this paper (r

Molecular Hydrogen Emission from Protoplanetary Disks: Effects of X-ray Irradiation and Dust Evolution

Detailed models for the density and temperature profiles of gas and dust in protoplanetary disks are constructed by taking into account X-ray and UV irradiation from a central T Tauri star, as well as dust size growth and settling toward the disk midplane. The spatial and size distributions of dust grains are numerically computed by solving the coagulation equation for settling dust particles, with the result that the mass and total surface area of dust grains per unit volume of the gas in the disks are very small, except at the midplane. The H2 level populations and line emission are calculated using the derived physical structure of the disks. X-ray irradiation is the dominant heating source of the gas in the inner disk and in the surface layer, while the UV heating dominates otherwise. If the central star has strong X-ray and weak UV radiation, the H2 level populations are controlled by X-ray pumping, and the X-rayinduced transition lines could be observable. If the UV irradiation is strong, the level populations are controlled by thermal collisions or UV pumping, depending on the dust properties. As the dust particles evolve in the disks, the gas temperature at the disk surface drops because the grain photoelectric heating becomes less efficient. This makes the level populations change from LTE to non-LTE distributions, which results in changes to the line ratios. Our results suggest that dust evolution in protoplanetary disks could be observable through the H2 line ratios. The emission lines are strong from disks irradiated by strong UV and X-rays and possessing small dust grains; such disks will be good targets in which to observe H2 emission.

Spectroscopy of astrophysical plasmas in the laboratory

Most of the matter in the universe is in the few form of a plasma. Over the past years physicists have produced laboratory plasmas that can mimic those observed in astrophysics. The best known is probably the tokamak, which has similar physical conditions and plasma processes to those found in collisionally dominated solar and stellar transition regions and coronae. Spectroscopy of such laboratory plasmas, in, particular at, ultraviolet and X-ray wavelengths, has greatly aided our understanding of their astrophysical counterparts. More recently, experiments have been performed on the Z Machine at the Sandia National Laboratory in the USA with the aim of creating, for the first time, steady-state photoionization-dominated plasmas that recreate the conditions found in some accretion-powered X-ray sources, such as X-ray binaries. In the future, experiments are envisaged with laser-produced plasmas at AWE Aldermaston that may be able to mimic the steady-state conditions found in high-energy accretion-powered sources, including the central regions of active galaxies.

Imaging the cool stars in the interacting binaries AE Aqr, BV Cen and V426 Oph

It is well known that magnetic activity in late-type stars increases with increasing rotation rate. Using inversion techniques akin to medical imaging, the rotationally broadened profiles from such stars can be used to reconstruct `Doppler images' of the distribution of cool, dark starspots on their stellar surfaces. Interacting binaries, however, contain some of the most rapidly rotating late-type stars known and thus provide important tests of stellar dynamo models. Furthermore, magnetic activity is thought to play a key role in their evolution, behaviour and accretion dynamics. Despite this, we know comparatively little about the magnetic activity and its influence on such binaries. In this review we summarise the concepts behind indirect imaging of these systems, and present movies of the starspot distributions on the cool stars in some interacting binaries. We conclude with a look at the future opportunities that such studies may provide.

A Brown Dwarf Mass Donor in an Accreting Binary

A long-standing and unverified prediction of binary star evolution theory is the existence of a population of white dwarfs accreting from substellar donor stars. Such systems ought to be common, but the difficulty of finding them, combined with the challenge of detecting the donor against the light from accretion, means that no donor star to date has a measured mass below the hydrogen burning limit. We applied a technique that allowed us to reliably measure the mass of the unseen donor star in eclipsing systems. We were able to identify a brown dwarf donor star, with a mass of 0.052 ± 0.002 solar mass. The relatively high mass of the donor star for its orbital period suggests that current evolutionary models may underestimate the radii of brown dwarfs.

Time-resolved optical photometry of the ultracompact binary 4U 0614+091

We present a detailed optical study of the ultracompact X-ray binary 4U 0614+091. We have used 63 hr of time-resolved optical photometry taken with three different telescopes (IAC80, NOT, and SPM) to search for optical modulations. The power spectra of each data set reveals sinusoidal modulations with different periods, which are not always present. The strongest modulation has a period of 51.3 minutes, a semiamplitude of 4.6 mmag, and is present in the IAC80 data. The SPM and NOT data show periods of 42 minutes and 64 minutes, respectively, but with much weaker amplitudes, 2.6 mmag and 1.3 mmag, respectively. These modulations arise from either X-ray irradiation of the inner face of the secondary star and/or a superhump modulation from the accretion disk, or quasiperiodic modulations in the accretion disk. It is unclear whether these periods/quasi-periodic modulations are related to the orbital period; however, the strongest period of 51.3 minutes is close to earlier tentative orbital periods. Further observations taken over a long baseline are encouraged.

The first polarimetric signatures of infrared jets in X-ray binaries

We present near-infrared linear spectropolarimetry of a sample of persistent X-ray binaries, Sco X-1, Cyg X-2, and GRS 1915+105. The slopes of the spectra are shallower than what is expected from a standard steady state accretion disk, and can be explained if the near-infrared flux contains a contribution from an optically thin jet. For the neutron star systems, Sco X-1 and Cyg X-2, the polarization levels at 2.4 mu m are 1.3% +/- 0.10% and 5.4% +/- 0.7%, respectively, which is greater than the polarization level at 1.65 mu m. This cannot be explained by interstellar polarization or electron scattering in the anisotropic environment of the accretion flow. We propose that the most likely explanation is that this is the polarimetric signature of synchrotron emission arising from close to the base of the jets in these systems. In the black hole system GRS 1915+105 the observed polarization, although high (5.0% +/- 1.2% at 2.4 mu m), may be consistent with interstellar polarization. For Sco X-1 the position angle of the radio jet on the sky is approximately perpendicular to the near-infrared position angle (electric vector), suggesting that the magnetic field is aligned with the jet. These observations may be a first step toward probing the ordering, alignment, and variability of the outflow magnetic field in a region closer to the central accreting object than is observed in the radio band.

ULTRACAM photometry of the eclipsing cataclysmic variables GY Cnc, IR Com and HT Cas

We present high-speed, three-colour photometry of the eclipsing cataclysmic variables GY Cnc, IR Com and HT Cas. We find that the sharp eclipses in GY Cnc and IR Com are due to eclipses of the white dwarf. There is some evidence for a bright-spot on the edge of the accretion disc in GY Cnc, but not in IR Com. Eclipse mapping of HT Cas is presented which shows changes in the structure of the quiescent accretion disc. Observations in 2002 show the accretion disc to be invisible except for the presence of a bright-spot at the disc edge. 2003 observations, however, clearly show a bright inner disc and the bright-spot to be much fainter than in 2002. Although no outburst was associated with either set of quiescent observations, the system was similar to 0.6 mJy brighter in 2003, mainly due to the enhanced emission from the inner disc. We propose that these changes are due to variations in the mass-transfer rate from the secondary star and through the disc. The disc colours indicate that it is optically thin in both its inner and outer regions. We estimate the white dwarf temperature of HT Cas to be 15 000 +/- 1000 K in 2002 and 14 000 +/- 1000 K in 2003.

The component masses of the cataclysmic variable V347 Puppis

We present time-resolved spectroscopy and photometry of the double-lined eclipsing cataclysmic variable V347 Pup ( = LB 1800). There is evidence of irradiation on the inner hemisphere of the secondary star, which we correct for using a model to give a secondary-star radial velocity of K R = 198 +/- 5 km s(-1). The rotational velocity of the secondary star in V347 Pup is found to be v sin i = 131 +/- 5 km s(-1) and the system inclination is i = 84degrees.0 +/- 2degrees.3. From these parameters we obtain masses of M-1 = 0.63 +/- 0.04 M for the white dwarf primary and M-2 = 0.52 +/- 0.06 M for the M0.5V secondary star, giving a mass ratio of q = 0.83 +/- 0.05. On the basis of the component masses, and the spectral type and radius of the secondary star in V347 Pup, we find tentative evidence for an evolved companion. V347 Pup shows many of the characteristics of the SW Sex stars, exhibiting single-peaked emission lines, high-velocity S-wave components and phase-offsets in the radial velocity curve. We find spiral arms in the accretion disc of V347 Pup and measure the disc radius to be close to the maximum allowed in a pressureless disc.

Stokes imaging, Doppler mapping and Roche tomography of the AM Herculis system V834 Cen

We report on new simultaneous phase-resolved spectroscopic and polarimetric observations of the polar (AM Herculis star) V834 Cen during a high state of accretion. Strong emission lines and high levels of variable circular and linear polarization are observed over the orbital period. The polarization data are modelled using the Stokes imaging technique of Potter et al. The spectroscopic emission lines are investigated using the Doppler tomography technique of Marsh and Horne and the Roche tomography technique of Dhillon and Watson. Up to now, all three techniques have been used separately to investigate the geometry and accretion dynamics in cataclysmic variables. For the first time, we apply all three techniques to simultaneous data for a single system. This allows us to compare and test each of the techniques against each other and hence to derive a better understanding of the geometry, dynamics and system parameters of V834 Cen.

Roche tomography of cataclysmic variables - II. Images of the secondary stars in AM Her, QQ Vul, IP Peg and HU Aqr

We present a set of Roche tomography reconstructions of the secondary stars in the cataclysmic variables AM Her, QQ Vul, IP Peg and HU Aqr. The image reconstructions show distinct asymmetries in the irradiation pattern for all four systems that can be attributed to shielding of the secondary star by the accretion stream/column in AM Her, QQ Vul and HU Aqr, and increased irradiation by the bright-spot in IP Peg. We use the entropy landscape technique to derive accurate system parameters (M-1, M-2, i and gamma) for the four binaries. In principle, this technique should provide the most reliable mass determinations available, since the intensity distribution across the secondary star is known. We also find that the intensity distribution can systematically affect the value of gamma derived from circular orbit fits to radial velocity variations.

The role of H2D+ in the deuteration of interstellar molecules

Collisions between H-3(+) and HD in molecular clouds lead to the fractionation of deuterium in H2D+ at temperatures below 20 K. In this article, we describe the chemistry of H2D+ and discuss how variations in temperature and elemental abundances affect the level of fractionation in H2D+ and other species. We describe how accretion of gas-phase molecules on to cold dust grains enhances the deuteration in several molecules including doubly deuterated molecules. Mie show that the ion-neutral drift velocities attained in slow Alfven waves can destroy H2D+ in non-thermal reactions. As a result, the degree of fractionation can be reduced and we discuss observational consequences of such a model for the dark dust cloud TMC-1.