Geophysical and Astrophysical Fluid Dynamics

Abstract:

We present the results of non-linear, hydrodynamic simulations, in three dimensions, of the tidal perturbation of accretion discs in binary systems where the orbit is circular and not necessarily coplanar with the disc mid-plane. The accretion discs are assumed to be geometrically thin, and of low mass relative to the stellar mass so that they are governed by thermal pressure and viscosity, but not self-gravity. The parameters that we consider in our models are the ratio of the orbital distance to the disc radius, D/R, the binary mass ratio, Ms/Mp, the initial inclination angle between the orbit and disc planes, δ, and the Mach number in the outer parts of the unperturbed disc, ℳ. Since we consider non-self-gravitating discs, these calculations are relevant to protostellar binaries with separations below a few hundred au. For binary mass ratios of around unity and D/R in the range 3 to 4, we find that the global evolution of the discs is governed primarily by the value of ℳ. For relatively low Mach numbers (i.e. ℳ = 10 to 20) we find that the discs develop a mildly warped structure, are tidally truncated, and undergo a near rigid body precession at a rate which is in close agreement with analytical arguments. For higher Mach numbers (ℳ ≈ 30), the evolution is towards a considerably more warped structure, but the disc none the less maintains itself as a long-lived, coherent entity. A further increase in Mach number to ℳ = 50 leads to a dramatic disruption of the disc as a result of differential precession, since the sound speed is too low to allow efficient communication between constituent parts of the disc. Additionally, it is found that the inclination angle between the disc and the orbital angular momentum vectors evolves on a longer time-scale, which is probably the viscous evolution time-scale of the disc. The calculations are relevant to a number of observed astrophysical phenomena, including the precession of jets associated with young stars, the high spectral index of some T Tauri stars, and the light curves of X-ray binaries such as Hercules X-1 which suggest the presence of precessing accretion discs.

Abstract:

We compute here the spectral energy distribution (SED) of warped T Tauri discs in a general way. In a previous paper we analytically calculated, in a linear approximation, the response of a circumstellar disc to tidal forces owing to a stellar companion in a non-coplanar young binary system. Here, we consider tidally induced warps of larger amplitude, and we use these previous results to parametrize the disc deformation. We then compute the energy emitted in a given direction by the system of the warped disc and central star, taking into account shadowing effects. We find that the parametrized warp model produces a broad variety of synthetic SEDs. Some of them are comparable to those of T Tauri stars with infrared excess (Class II sources), whereas others resemble Class I protostellar sources. By comparing models with actual observations of both a T Tauri star with high spectral index and a Class I source, we find that the derived warp and disc parameters are not unrealistic, and we conclude that tidal interactions in T Tauri binary systems with intermediate separations must play a role in shaping the SEDs of these stars.

Abstract:

A meteorological data assimilation scheme has been developed for Mars, based on techniques used for the current operational scheme at the UK Meteorological Office (UKMO). The scheme has been designed to interface with a range of models of varying complexity, from a simple primitive equations model to a full martian general circulation model with detailed parameterizations of the most important physical processes. The scheme was originally intended to be used primarily for analysis of temperature profiles from the Pressure Modulator Infrared Radiometer (PMIRR) onboard the Mars Observer (MO) spacecraft. It is, however, capable of analysing any asynoptic, randomly-distributed dataset, which may include surface pressure, temperature, velocity or constituent data originating either from other past or future spacecraft missions or from models. It is proposed to use the scheme not only for interpretation of forthcoming data from future spacecraft but also for model inter-comparisons, for hind-casting earlier martian data for model validation and for data impact studies in the planning of new missions. Initial results have demonstrated successful assimilations of trial orbiter and lander data into models with different initial conditions and physics, although a large number of landers might be required if orbiter data were not available. An optimum approach would combine orbiter data with simultaneous measurements from a small number of landers.