Early-type stars observed in the ESO UVES Paranal Observatory Project - III. Sub-parsec and AU-scale structure in the interstellar medium

UVES interstellar observations from the Paranal Observatory Project are presented for early-type stars located in the line of sight to the nearby open clusters IC 2391 (Omni Vel) and NGC 6475 (M7), with spectroscopic resolution R similar to 80 000 and signal-to-noise ratios in the Ti II (3383 angstrom), Ca II K, CH+ (4232 angstrom), Na I D and K I (7698 angstrom) lines of several hundred. The sightlines are a mixture of cluster and non-cluster objects. A total of 22 early-type stars (A and B type) are present in our sample towards IC 2391, with 21 towards NGC 6475/M7, and enable us to probe for differences in column density on scales from similar to 0.07 to 7.3 and similar to 0.05 to 4.9 pc in the respective clusters. Additionally, towards Praesepe the Na I D interstellar variation only is probed towards 13 sightlines and transverse scales of similar to 0.16-10.7 pc at R = 70 000. Towards IC 2391 variations are found in Ti II, Ca II K and Na I D column density in different sightlines of up to 0.7, 1.0 and 1.8 dex (excluding one star), respectively. This kind of variability correlates well with the Hipparcos parallax of the objects, and probes structure within the Local Bubble. For cluster-only objects the variations are 0.3, 0.3 and 0.5 dex, respectively. For the field of view towards NGC6475 the corresponding maximum variations are somewhat smaller, being 0.5, 0.3, 0.8 and 1.0 dex for Ti II, Ca II K, Na I and K I, respectively, for all objects and 0.4, 0.2, 0.6 and 0.7 dex for the cluster-only objects. These are uncorrelated with parallax, and again demonstrate that Ca II K tends to be more smoothly distributed than Na I D. A few likely cluster sightlines show evidence for CH+ and variations in this molecular species of a factor of 10 in equivalent width over sub-pc scales. Towards Praesepe variation in interstellar Na I D is small, being a maximum of only similar to 0.4 dex (including measurement errors), but with fewer sightlines studied. Overall, the scatter in the data is similar for the singly ionized species Ti II and Ca II, lending more support to the hypothesis that these two species sample similar parts of the interstellar medium (ISM). This also appears to be the case for the neutral species Na I D and K I in the one cluster studied. Finally, multiple-epoch observations from a variety of archive sources are used to search for astronomical unit (au) scale structure in the ISM towards 46 sightlines. There are tentative indications of structure on scales of tens to thousands of au for three sightlines. Future observations will confirm the veracity or otherwise of the time-variable components and others presented.

Recherche d’étoiles jeunes de faible masse dans le voisinage solaire

Formées lors de l’effondrement gravitationnel d’un nuage de gaz moléculaire, les étoiles naissantes auront différentes masses variant entre 0.08 et environ 100M . La majorité de la population stellaire de la Galaxie est constituée d’étoiles dont la masse est inférieure à environ 0.6 M . Le dernier évènement de formation stellaire dans le voisinage solaire s’est produit dans la bulle locale il y a au plus 100 millions d’années, vraisemblablement provoqué par le passage d’une onde de choc dans le bras local de la Galaxie. C’est ainsi que se formèrent de jeunes associations d’étoiles dont les membres se caractérisent en particulier par une vitesse spatiale et une position commune dans la Galaxie. Les associations jeunes étant peu densément peuplées et relativement proches du Soleil, leurs membres se font plutôt rares et dispersés sur toute la voûte céleste. Jusqu’à présent, surtout les étoiles les plus massives (brillantes) ont été répertoriées. Les étoiles jeunes de faible masse, constituant la majorité de la population, restent pour la plupart à être identifiées. Les étoiles jeunes de faible masse représentent une population clef pour contraindre les modèles évolutifs des étoiles M et des naines brunes. Elles sont également d’excellentes candidates pour chercher des exoplanètes via les techniques d’imagerie directe. Ce mémoire présente une nouvelle méthode utilisant un modèle cinématique enrichi d’une analyse statistique Bayesienne pour identifier des étoiles jeunes de faible masse dans les associations beta Pictoris, Tucana-Horologium et AB Doradus. À partir d’un échantillon de 1080 étoiles K et M, toutes comportant des indicateurs de jeunesse tels l’émission Halpha et une forte luminosité dans les rayons X, leurs propriétés cinématiques (mouvement propre) et photométriques sont analysées pour en extraire 98 candidates hautement probables membres d’une des trois associations. Une confirmation de leur statut comme membre nécessitera en particulier une mesure de leur vitesse radiale (prédit par notre analyse) et une mesure de la largeur équivalente du lithium à 6708 Å pour mieux contraindre leur âge.

COMPARISONS OF THE INTERSTELLAR MAGNETIC FIELD DIRECTIONS OBTAINED FROM THE IBEX RIBBON AND INTERSTELLAR POLARIZATIONS

Variations in the spatial configuration of the interstellar magnetic field (ISMF) near the Sun can be constrained by comparing the ISMF direction at the heliosphere found from the Interstellar Boundary Explorer (IBEX) spacecraft observations of a ""Ribbon"" of energetic neutral atoms (ENAs), with the ISMF direction derived from optical polarization data for stars within similar to 40 pc. Using interstellar polarization observations toward similar to 30 nearby stars within similar to 90 degrees of the heliosphere nose, we find that the best fits to the polarization position angles are obtained for a magnetic pole directed toward ecliptic coordinates of lambda, beta similar to 263 degrees, 37 degrees (or galactic coordinates of l, b similar to 38 degrees, 23 degrees), with uncertainties of +/- 35 degrees based on the broad minimum of the best fits and the range of data quality. This magnetic pole is 33 degrees from the magnetic pole that is defined by the center of the arc of the ENA Ribbon. The IBEX ENA ribbon is seen in sight lines that are perpendicular to the ISMF as it drapes over the heliosphere. The similarity of the polarization and Ribbon directions for the local ISMF suggests that the local field is coherent over scale sizes of tens of parsecs. The ISMF vector direction is nearly perpendicular to the flow of local interstellar material (ISM) through the local standard of rest, supporting a possible local ISM origin related to an evolved expanding magnetized shell. The local ISMF direction is found to have a curious geometry with respect to the cosmic microwave background dipole moment.

THE INTERSTELLAR MAGNETIC FIELD CLOSE TO THE SUN. II.

The magnetic field in the local interstellar medium (ISM) provides a key indicator of the galactic environment of the Sun and influences the shape of the heliosphere. We have studied the interstellar magnetic field (ISMF) in the solar vicinity using polarized starlight for stars within 40 pc of the Sun and 90 degrees of the heliosphere nose. In Frisch et al. (Paper I), we developed a method for determining the local ISMF direction by finding the best match to a group of interstellar polarization position angles obtained toward nearby stars, based on the assumption that the polarization is parallel to the ISMF. In this paper, we extend the analysis by utilizing weighted fits to the position angles and by including new observations acquired for this study. We find that the local ISMF is pointed toward the galactic coordinates l, b = 47 degrees +/- 20 degrees, 25 degrees +/- 20 degrees. This direction is close to the direction of the ISMF that shapes the heliosphere, l, b = 33 degrees +/- 4 degrees, 55 degrees +/- 4 degrees, as traced by the center of the "Ribbon" of energetic neutral atoms discovered by the Interstellar Boundary Explorer (IBEX) mission. Both the magnetic field direction and the kinematics of the local ISM are consistent with a scenario where the local ISM is a fragment of the Loop I superbubble. A nearby ordered component of the local ISMF has been identified in the region l approximate to 0 degrees -> 80 degrees and b approximate to 0 degrees -> 30 degrees, where PlanetPol data show a distance-dependent increase of polarization strength. The ordered component extends to within 8 pc of the Sun and implies a weak curvature in the nearby ISMF of +/- 0 degrees.25 pc(-1). This conclusion is conditioned on the small sample of stars available for defining this rotation. Variations from the ordered component suggest a turbulent component of +/- 23 degrees. The ordered component and standard relations between polarization, color excess, and H-o column density predict a reasonable increase of N(H) with distance in the local ISM. The similarity of the ISMF directions traced by the polarizations, the IBEX Ribbon, and pulsars inside the Local Bubble in the third galactic quadrant suggest that the ISMF is relatively uniform over spatial scales of 8-200 pc and is more similar to interarm than spiral-arm magnetic fields. The ISMF direction from the polarization data is also consistent with small-scale spatial asymmetries detected in GeV-TeV cosmic rays with a galactic origin. The peculiar geometrical relation found earlier between the cosmic microwave background dipole moment, the heliosphere nose, and the ISMF direction is supported by this study. The interstellar radiation field at +/- 975 angstrom does not appear to play a role in grain alignment for the low-density ISM studied here.

The locations of recent supernovae near the Sun from modelling 60Fe transport

The signature of 60Fe in deep-sea crusts indicates that one or more supernovae exploded in the solar neighbourhood about 2.2 million years ago1–4. Recent isotopic analysis is consistent with a core-collapse or electron-capture supernova that occurred 60 to 130 parsecs from the Sun5. Moreover, peculiarities in the cosmic ray spectrum point to a nearby supernova about two million years ago6. The Local Bubble of hot, diffuse plasma, in which the Solar System is embedded, originated from 14 to 20 supernovae within a moving group, whose surviving members are now in the Scorpius– Centaurus stellar association7,8. Here we report calculations of the most probable trajectories and masses of the supernova progenitors, and hence their explosion times and sites. The 60Fe signal arises from two supernovae at distances between 90 and 100 parsecs. The closest occurred 2.3 million years ago at present-day galactic coordinates l = 327°, b = 11°, and the second-closest exploded about 1.5 million years ago at l = 343°, b = 25°, with masses of 9.2 and 8.8 times the solar mass, respectively. The remaining supernovae, which formed the Local Bubble, contribute to a smaller extent because they happened at larger distances and longer ago (60Fe has a half- life of 2.6 million years9,10). There are uncertainties relating to the nucleosynthesis yields and the loss of 60Fe during transport, but they do not influence the relative distribution of 60Fe in the crust layers, and therefore our model reproduces the measured relative abundances very well.

Adaptive finite element analysis of elliptic problems based on bubble-type local mesh generation

A new mesh adaptivity algorithm that combines a posteriori error estimation with bubble-type local mesh generation (BLMG) strategy for elliptic differential equations is proposed. The size function used in the BLMG is defined on each vertex during the adaptive process based on the obtained error estimator. In order to avoid the excessive coarsening and refining in each iterative step, two factor thresholds are introduced in the size function. The advantages of the BLMG-based adaptive finite element method, compared with other known methods, are given as follows: the refining and coarsening are obtained fluently in the same framework; the local a posteriori error estimation is easy to implement through the adjacency list of the BLMG method; at all levels of refinement, the updated triangles remain very well shaped, even if the mesh size at any particular refinement level varies by several orders of magnitude. Several numerical examples with singularities for the elliptic problems, where the explicit error estimators are used, verify the efficiency of the algorithm. The analysis for the parameters introduced in the size function shows that the algorithm has good flexibility.

Relevance of local parallel theory to the linear stability of laminar separation bubbles

Laminar separation bubbles are thought to be highly non-parallel, and hence global stability studies start from this premise. However, experimentalists have always realized that the flow is more parallel than is commonly believed, for pressure-gradient-induced bubbles, and this is why linear parallel stability theory has been successful in describing their early stages of transition. The present experimental/numerical study re-examines this important issue and finds that the base flow in such a separation bubble becomes nearly parallel due to a strong-interaction process between the separated boundary layer and the outer potential flow. The so-called dead-air region or the region of constant pressure is a simple consequence of this strong interaction. We use triple-deck theory to qualitatively explain these features. Next, the implications of global analysis for the linear stability of separation bubbles are considered. In particular we show that in the initial portion of the bubble, where the flow is nearly parallel, local stability analysis is sufficient to capture the essential physics. It appears that the real utility of the global analysis is perhaps in the rear portion of the bubble, where the flow is highly non-parallel, and where the secondary/nonlinear instability stages are likely to dominate the dynamics.

Unsteady separation and vortex shedding from a laminar separation bubble over a bluff body

Boundary layers are subject to favorable and adverse pressure gradients because of both the temporal and spatial components of the pressure gradient. The adverse pressure gradient may cause the flow to separate. In a closed loop unsteady tunnel we have studied the initiation of separation in unsteady flow past a constriction (bluff body) in a channel. We have proposed two important scalings for the time when boundary layer separates. One is based on the local pressure gradient and the other is a convective time scale based on boundary layer parameters. The flow visualization using a dye injection technique shows the flow structure past the body. Nondimensional shedding frequency (Strouhal number) is calculated based on boundary layer and momentum thicknesses. Strouhal number based on the momentum thickness shows a close agreement with that for flat plate and circular cylinder.

Bubble Behavior and Heat Transfer in Quasi-Steady Pool Boiling in Microgravity

Pool boiling of degassed FC-72 on a plane plate heater has been studied experimentally in microgravity. A quasi-steady heating method is adopted, in which the heating voltage is controlled to increase exponentially with time. Compared with terrestrial experiments, bubble behaviors are very different, and have direct effect on heat transfer. Small, primary bubbles attached on the surface seem to be able to suppress the activation of the cavities in the neighborhoods, resulting in a slow increase of the wall temperature with the heat flux. For the high subcooling, the coalesced bubble has a smooth surface and a small size. It is difficult to cover the whole heater surface, resulting in a special region of gradual transitional boiling in which nucleate boiling and local dry area can co-exist. No turning point corresponding to the transition from nucleate boiling to film boiling can be observed. On the contrary, the surface oscillation of the coalesced bubble at low subcooling may cause more activated nucleate sites, and then the surface temperature may keep constant or even fall down with the increasing heat flux. Furthermore, an abrupt transition to film boiling can also be observed. It is shown that heat transfer coefficient and CHF increase with the subcooling or pressure in microgravity, as observed in normal gravity.

Horizontal Bloch line motion in magnetic bubble materials

The purpose of this work is to extend experimental and theoretical understanding of horizontal Bloch line (HBL) motion in magnetic bubble materials. The present theory of HBL motion is reviewed, and then extended to include transient effects in which the internal domain wall structure changes with time. This is accomplished by numerically solving the equations of motion for the internal azimuthal angle ɸ and the wall position q as functions of z, the coordinate perpendicular to the thin-film material, and time. The effects of HBL's on domain wall motion are investigated by comparing results from wall oscillation experiments with those from the theory. In these experiments, a bias field pulse is used to make a step change in equilibrium position of either bubble or stripe domain walls, and the wall response is measured by using transient photography. During the initial response, the dynamic wall structure closely resembles the initial static structure. The wall accelerates to a relatively high velocity (≈20 m/sec), resulting in a short (≈22 nsec ) section of initial rapid motion. An HBL gradually forms near one of the film surfaces as a result of local dynamic properties, and moves along the wall surface toward the film center. The presence of this structure produces low-frequency, triangular-shaped oscillations in which the experimental wall velocity is nearly constant, v_s≈ 5-8 m/sec. If the HBL reaches the opposite surface, i.e., if the average internal angle reaches an integer multiple of π, the momentum stored in the HBL is lost, and the wall chirality is reversed. This results in abrupt transitions to overdamped motion and changes in wall chirality, which are observed as a function of bias pulse amplitude. The pulse amplitude at which the n^th punch- through occurs just as the wall reaches equilibrium is given within 0.2 0e by H_n = (2v_sH'/γ)^1/2 • (nπ)^1/2 + H_sv), where H' is the effective field gradient from the surrounding domains, and H_sv is a small (less than 0.03 0e), effective drag field. Observations of wall oscillation in the presence of in-plane fields parallel to the wall show that HBL formation is suppressed by fields greater than about 40 0e (≈2πM_s), resulting in the high-frequency, sinusoidal oscillations associated with a simple internal wall structure.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

INVESTIGATION OF BUBBLE DYNAMICS AND HEATING DURING FOCUSED ULTRASOUND INSONATION IN TISSUE-MIMICKING MATERIALS

The deposition of ultrasonic energy in tissue can cause tissue damage due to local heating. For pressures above a critical threshold, cavitation will occur in tissue and bubbles will be created. These oscillating bubbles can induce a much larger thermal energy deposition in the local region. Traditionally, clinicians and researchers have not exploited this bubble-enhanced heating since cavitation behavior is erratic and very difficult to control. The present work is an attempt to control and utilize this bubble-enhanced heating. First, by applying appropriate bubble dynamic models, limits on the asymptotic bubble size distribution are obtained for different driving pressures at 1 MHz. The size distributions are bounded by two thresholds: the bubble shape instability threshold and the rectified diffusion threshold. The growth rate of bubbles in this region is also given, and the resulting time evolution of the heating in a given insonation scenario is modeled. In addition, some experimental results have been obtained to investigate the bubble-enhanced heating in an agar and graphite based tissue- mimicking material. Heating as a function of dissolved gas concentrations in the tissue phantom is investigated. Bubble-based contrast agents are introduced to investigate the effect on the bubble-enhanced heating, and to control the initial bubble size distribution. The mechanisms of cavitation-related bubble heating are investigated, and a heating model is established using our understanding of the bubble dynamics. By fitting appropriate bubble densities in the ultrasound field, the peak temperature changes are simulated. The results for required bubble density are given. Finally, a simple bubbly liquid model is presented to estimate the shielding effects which may be important even for low void fraction during high intensity focused ultrasound (HIFU) treatment.

Are mortgage lenders guilty of the housing bubble? A UK perspective

Existing theoretical models of house prices and credit rely on continuous rationality of consumers, an assumption that has been frequently questioned in recent years. Meanwhile, empirical investigations of the relationship between prices and credit are often based on national-level data, which is then tested for structural breaks and asymmetric responses, usually with subsamples. Earlier author argues that local markets are structurally different from one another and so the coefficients of any estimated housing market model should vary from region to region. We investigate differences in the price–credit relationship for 12 regions of the UK. Markov-switching is introduced to capture asymmetric market behaviours and turning points. Results show that credit abundance had a large impact on house prices in Greater London and nearby regions alongside a strong positive feedback effect from past house price movements. This impact is even larger in Greater London and the South East of England when house prices are falling, which are the only instances where the credit effect is more prominent than the positive feedback effect. A strong positive feedback effect from past lending activity is also present in the loan dynamics. Furthermore, bubble probabilities extracted using a discrete Kalman filter neatly capture market turning points.