Exoplanets or Dynamic Atmospheres? The Radial Velocity and Line Shape Variations of 51 Pegasi and tau Bootis

The stars 51 Pegasi and tau Bootis show radial velocity variations that have been interpreted as resulting from companions with roughly Jovian mass and orbital periods of a few days. Gray and Gray & Hatzes reported that the radial velocity signal of 51 Peg is synchronous with variations in the shape of the line lambda 6253 Fe I; thus, they argue that the velocity signal arises not from a companion of planetary mass but from dynamic processes in the atmosphere of the star, possibly nonradial pulsations. Here we seek confirming evidence for line shape or strength variations in both 51 Peg and tau Boo, using R = 50,000 observations taken with the Advanced Fiber Optic Echelle. Because of our relatively low spectral resolution, we compare our observations with Gray's line bisector data by fitting observed line profiles to an expansion in terms of orthogonal (Hermite) functions. To obtain an accurate comparison, we model the emergent line profiles from rotating and pulsating stars, taking the instrumental point-spread function into account. We describe this modeling process in detail. We find no evidence for line profile or strength variations at the radial velocity period in either 51 Peg or in tau Boo. For 51 Peg, our upper limit for line shape variations with 4.23 day periodicity is small enough to exclude with 10 sigma confidence the bisector curvature signal reported by Gray & Hatzes; the bisector span and relative line depth signals reported by Gray are also not seen, but in this case with marginal (2 sigma ) confidence. We cannot, however, exclude pulsations as the source of 51 Peg's radial velocity variation because our models imply that line shape variations associated with pulsations should be much smaller than those computed by Gray & Hatzes; these smaller signals are below the detection limits both for Gray & Hatzes's data and for our own. tau Boo's large radial velocity amplitude and v sin i make it easier to test for pulsations in this star. Again we find no evidence for periodic line shape changes, at a level that rules out pulsations as the source of the radial velocity variability. We conclude that the planet hypothesis remains the most likely explanation for the existing data.

Liquid-liquid slug flow: Hydrodynamics and pressure drop

In this paper, the hydrodynamics and the pressure drop of liquid-liquid slug flow in round microcapillaries are presented. Two liquid-liquid flow systems are considered, viz. water-toluene and ethylene glycol/water-toluene. The slug lengths of the alternating continuous and dispersed phases were measured as a function of the slug velocity (0.03-0.5 m/s), the organic-to-aqueous flow ratio (0.1-4.0), and the microcapillary internal diameter (248 and 498 mu m). The pressure drop is modeled as the sum of two contributions: the frictional and the interface pressure drop. Two models are presented, viz, the stagnant film model and the moving film model. Both models account for the presence of a thin liquid film between the dispersed phase slug and the capillary wall. It is found that the film velocity is of negligible influence on the pressure drop. Therefore, the stagnant film model is adequate to accurately predict the liquid-liquid slug flow pressure drop. The influence of inertia and the consequent change of the slug cap curvature are accounted for by modifying Bretherton's curvature parameter in the interface pressure drop equation. The stagnant film model is in good agreement with experimental data with a mean relative error of less than 7%.

Implementation of Menter's Transition Model on an Isolated Natural Laminar Flow Nacelle

This study evaluates the implementation of Menter's gamma-Re-theta Transition Model within the CFX12 solver for turbulent transition prediction on a natural laminar flow nacelle. Some challenges associated with this type of modeling have been identified. The computational fluid dynamics transitional flow simulation results are presented for a series of cruise cases with freestream Mach numbers ranging from 0.8 to 0.88, angles of attack from 2 to 0 degrees, and mass flow ratios from 0.60 to 0.75. These were validated with a series of wind-tunnel tests on the nacelle by comparing the predicted and experimental surface pressure distributions and transition locations. A selection of the validation cases are presented in this paper. In all cases, computational fluid dynamics simulations agreed reasonably well with the experiments. The results indicate that Menter's gamma-Re-theta Transition Model is capable of predicting laminar boundary-layer transition to turbulence on a nacelle. Nonetheless, some limitations exist in both the Menter's gamma-Re-theta Transition Model and in the implementation of the computational fluid dynamics model. The implementation of a more comprehensive experimental correlation in Menter's gamma-Re-theta Transition Model, preferably the ones from nacelle experiments, including the effects of compressibility and streamline curvature, is necessary for an accurate transitional flow simulation on a nacelle. In addition, improvements to the computational fluid dynamics model are also suggested, including the consideration of varying distributed surface roughness and an appropriate empirical correction derived from nacelle experimental transition location data.

Response characteristics of optical sensors for oxygen: models based on a distribution in tau(o) or kappa(q)

The features of two popular models used to describe the observed response characteristics of typical oxygen optical sensors based on luminescence quenching are examined critically. The models are the 'two-site' and 'Gaussian distribution in natural lifetime, tau(o),' models. These models are used to characterise the response features of typical optical oxygen sensors; features which include: downward curving Stern-Volmer plots and increasingly non-first order luminescence decay kinetics with increasing partial pressures of oxygen, pO(2). Neither model appears able to unite these latter features, let alone the observed disparate array of response features exhibited by the myriad optical oxygen sensors reported in the literature, and still maintain any level of physical plausibility. A model based on a Gaussian distribution in quenching rate constant, k(q), is developed and, although flawed by a limited breadth in distribution, rho, does produce Stern-Volmer plots which would cover the range in curvature seen with real optical oxygen sensors. A new 'log-Gaussian distribution in tau(o) or k(q)' model is introduced which has the advantage over a Gaussian distribution model of placing no limitation on the value of rho. Work on a 'log-Gaussian distribution in tau(o)' model reveals that the Stern-Volmer quenching plots would show little degree in curvature, even at large rho values and the luminescence decays would become increasingly first order with increasing pO(2). In fact, with real optical oxygen sensors, the opposite is observed and thus the model appears of little value. In contrast, a 'log-Gaussian distribution in k(o)' model does produce the trends observed with real optical oxygen sensors; although it is technically restricted in use to those in which the kinetics of luminescence decay are good first order in the absence of oxygen. The latter model gives a good fit to the major response features of sensors which show the latter feature, most notably the [Ru(dpp)(3)(2+)(Ph4B-)(2)] in cellulose optical oxygen sensors. The scope of a log-Gaussian model for further expansion and, therefore, application to optical oxygen sensors, by combining both a log-Gaussian distribution in k(o) with one in tau(o) is briefly discussed.

A model for radiation-induced bystander effects, with allowance for spatial position and the effects of cell turnover

Bystander effects, whereby cells that are not directly exposed to ionizing radiation exhibit adverse biological effects, have been observed in a number of experimental systems. A novel stochastic model of the radiation-induced bystander effect is developed that takes account of spatial location, cell killing and repopulation. The ionizing radiation dose- and time-responses of this model are explored, and it is shown to exhibit pronounced downward curvature in the high dose-rate region, similar to that observed in many experimental systems, reviewed in the paper. It is also shown to predict the augmentation of effect after fractionated delivery of dose that has been observed in certain experimental systems. It is shown that the generally intractable solution of the full stochastic system can be considerably simplified by assumption of pairwise conditional dependence that varies exponentially over time. (C) 2004 Elsevier Ltd. All rights reserved.

The geometry of unstented and stented pig common carotid artery bypass grafts

The long-term success of arterial bypass grafting with autologous saphenous veins is limited by neointimal hyperplasia (NIH), which seemingly develops preferentially at sites where hydrodynamic wall shear is low. Placement of a loose-fitting, porous stent around end-to-end, or end-to-side, autologous saphenous vein grafts on the porcine common carotid artery has been found significantly to reduce NIH, but the mechanism is unclear. In a preliminary study, we implanted autologous saphenous vein grafts bilaterally on the common carotid arteries of pigs, placing a stent around one graft and leaving the contralateral graft unstented. At sacrifice 1 month post implantation, the grafts were pressure fixed in situ and resin casts were made. Unstented graft geometry was highly irregular, with non-uniform dilatation, substantial axial lengthening, curvature, kinking, and possible long-pitch helical distortion. In contrast, stented grafts showed no major dilatation, lengthening or curvature, but there was commonly fine corrugation, occasional slight kinking or narrowing of segments, and possible long-pitch helical distortion. Axial growth of grafts against effectively tethered anastomoses could account for these changes. CFD studies are planned, using 3D MR reconstructions, on the effects of graft geometry on the flow. Abnormality of the flow could favour the development of vascular pathology, including NIH.

The behaviour of damage tolerant hat-stiffened composite panels loaded in uniaxial compression

Damage tolerant hat-stiffened thin-skinned composite panels with and without a centrally located circular cutout, under uniaxial compression loading, were investigated experimentally and analytically. These panels incorporated a highly postbuckling design characterised by two integral stiffeners separated by a large skin bay with a high width to skin-thickness ratio. In both configurations, the skin initially buckled into three half-wavelengths and underwent two mode-shape changes; the first a gradual mode change characterised by a central deformation with double curvature and the second a dynamic snap to five half-wavelengths. The use of standard path-following non-linear finite element analysis did not consistently capture the dynamic mode change and an approximate solution for the prediction of mode-changes using a Marguerre-type Rayleigh-Ritz energy method is presented. Shortcomings with both methods of analysis are discussed and improvements suggested. The panels failed catastrophically and their strength was limited by the local buckling strength of the hat stiffeners. (C) 2001 Elsevier Science Ltd. All rights reserved.

Whole wafer magnetostriction metrology for magnetic films and multilayers

The requirements for metrology of magnetostriction in complex multilayers and on whole wafers present challenges. An elegant technique based on radius of curvature deformation of whole wafers in a commercial metrology tool is described. The method is based on the Villari effect through application of strain to a film by introducing a radius of curvature. Strain can be applied tensilely and compressively depending on the material. The design, while implemented on 3'' wafers, is scalable. The approach removes effects arising from any shape anisotropy that occurs with smaller samples, which can lead to a change in magnetic response. From the change in the magnetic anisotropy as a function of the radius, saturation magnetostriction ?s can be determined. Dependence on film composition and film thickness was studied to validate the radius of curvature approach with other techniques. ?s decreases from positive values to negative values through an increase in Ni concentration around the permalloy composition, and ?s also increases with a decrease in film thickness, in full agreement with previous reports. We extend the technique by demonstrating the technique applied to a multi-layered structure. These results verify the validity of the method and are an important step to facilitate further work in understanding how manipulation of multilayered films can offer tailored magnetostriction.

Tip-enhanced Raman microscopy: practicalities and limitations

The feasibility of apertureless scanning near-field Raman microscopy, exploiting the local enhancement in Raman scattering in the vicinity of a silver or gold tip, was investigated. Using the finite difference time domain method we calculated the enhancement of electric field strength, and hence Raman scattering, achieved through the resonant excitation of local modes in the tip. By modelling the frequency-dependent dielectric response of the metal tip we were able to highlight the resonant nature of the tip-enhancement and determine the excitation wavelength required for the strongest electric field enhancement, and hence Raman scattering intensity, which occurs for the excitation of modes localized at the tip apex. It is demonstrated that a peak Raman enhancement of 10(7)-fold should be achievable with <5 nm spatial resolution. We show that surface-enhanced Raman scattering from carbon contamination on a silver or gold tip can be significant. However, we find for a tip of radius of curvature 20 nm that the Raman enhancement should decay totally within 20 nm from the tip. Hence withdrawal of the tip by this distance should lead to the disappearance of the tip-enhanced signal, leaving only that from carbon contamination on the tip itself and the intrinsic signal from the sample. Copyright (C) 2003 John Wiley Sons, Ltd.

Inter-observer agreement in evaluating the anterior chamber angle by ultrasound biomicroscopy

Purpose. To determine the agreement between observers in estimating the configuration of the human anterior chamber angle (ACA) using ultrasound biomicroscopy (UBM). Methods. Two masked clinicians used UBM images to estimate, in 41 eyes, the configuration of the ACA, especially (1) the position of contact between the peripheral iris and the inside of the eye wall, (2) the angularity of the approach to the ACA, and (3) the curvature of the peripheral iris. Agreement between observers was evaluated by the kappa statistic. Results. Inter-observer agreement in assessing the iris insertion (kappa = 0.79), angular width (Kappa = 0.95), and the peripheral iris curvature (kappa = 0.84) was high. Conclusions. The agreement between observers in evaluating the anterior chamber angle configuration by UBM was excellent.

Intraobserver and interobserver agreement in evaluating the anterior chamber angle configuration by ultrasound biomicroscopy

Purpose: To determine the intra- and interobserver agreement in assessing the configuration of the human anterior chamber angle using ultrasound biomicroscopy (UBM). Methods: Two masked clinicians used ubm images to estimate, in 41 eyes, (a) the position of contact between the peripheral iris and the inside of the eye wall, (b) the angular size of the anterior chamber angle (ACA), and (c) the curvature of the peripheral iris. Both observers, masked to the previous results, examined the same images in a second session. Agreement was evaluated using the unweighted ? statistic. Results: Intraobserver agreement in assessing the iris insertion, angular width, and the iris curvature was high (range of ? values, 0.83-0.92). Interobserver agreement in evaluating the level of iris insertion (? = 0.79), the angular width (? = 0.95), and the iris curvature (? = 0.84) was also high. Conclusion: The agreement within the same observer and between observers in evaluating the ACA configuration by UBM was excellent.

LINEAR ACCELERATION EMISSION IN PULSAR MAGNETOSPHERES

Linear acceleration emission occurs when a charged particle is accelerated parallel to its velocity. We evaluate the spectral and angular distribution of this radiation for several special cases, including constant acceleration (hyperbolic motion) of finite duration. Based on these results, we find the following general properties of the emission from an electron in a linear accelerator that can be characterized by an electric field E acting over a distance L: (1) the spectrum extends to a cutoff frequency (h) over bar omega(c)/mc(2) approximate to L(E/E(Schw))(2)/(lambda) over bar (C), where E(Schw) = 1.3 x 10(18) V m(-1) is the Schwinger critical field and (lambda) over bar (C) = (h) over bar /mc = 3.86 x 10(-13) m is the Compton wavelength of the electron, (2) the total energy emitted by a particle traversing the accelerator is 4/3 alpha(f)(h) over bar omega(c) in accordance with the standard Larmor formula where alpha(f) is the fine-structure constant, and (3) the low frequency spectrum is flat for hyperbolic trajectories, but in general depends on the details of the accelerator. We also show that linear acceleration emission complements curvature radiation in the strongly magnetized pair formation regions in pulsar magnetospheres. It dominates when the length L of the accelerator is less than the formation length rho/gamma of curvature photons, where rho is the radius of curvature of the magnetic field lines and gamma the Lorentz factor of the emitting particle. In standard static models of pair creating regions linear acceleration emission is negligible, but it is important in more realistic dynamical models in which the accelerating field fluctuates on a short length scale.

Implementation of a non-orthogonal model for the finite element simulation of textile composite draping

In the pursuit of producing high quality, low-cost composite aircraft structures, out-of-autoclave manufacturing processes for textile reinforcements are being simulated with increasing accuracy. This paper focuses on the continuum-based, finite element modelling of textile composites as they deform during the draping process. A non-orthogonal constitutive model tracks yarn orientations within a material subroutine developed for Abaqus/Explicit, resulting in the realistic determination of fabric shearing and material draw-in. Supplementary material characterisation was experimentally performed in order to define the tensile and non-linear shear behaviour accurately. The validity of the finite element model has been studied through comparison with similar research in the field and the experimental lay-up of carbon fibre textile reinforcement over a tool with double curvature geometry, showing good agreement.

Barrier characteristics of PtSi/p-Si schottky diodes as determined from I-V-T measurements

The current-voltage-temperature characteristics of PtSi/p-Si Schottky barrier diodes were measured in the temperature range 60-115 K. Deviation of the ideality factor from unity below 80 K may be modelled using the so-called T-0 parameter with T-0 = 18 K. It is also shown that the curvature in the Richardson plots may be remedied by using the flatband rather than the zero-bias saturation current density. Physically, the departure from ideality is interpreted in terms of an inhomogeneous Schottky contact. Here we determine a mean barrier height at T = 0 K, phi(b)(-0) = 223 mV, with an (assumed) Gaussian distribution of standard deviation sigma(phi) = 12.5 mV. These data are correlated with the zero-bias barrier height, phi(j)(0) = 192 mV (at T = 90 K), the photoresponse barrier height, phi(ph) = 205 mV, and the flatband barrier height, phi(fb) = 214 mV. Finally, the temperature coefficient of the flatband barrier was found to be -0.121 mV K-1, which is approximately equal to 1/2(dE(g)(i)/dT), thus suggesting that the Fermi level at the interface is pinned to the middle of the band gap.

Kinetics of levoglucosan and formaldehyde formation during cellulose pyrolysis process

The mechanisms and kinetics studies of the formation of levoglucosan and formaldehyde from anhydroglucose radical have been carried out theoretically in this paper. The geometries and frequencies of all the stationary points are calculated at the B3LYP/6-31+G(D,P) level based on quantum mechanics, Six elementary reactions are found, and three global reactions are involved. The variational transition-state rate constants for the elementary reactions are calculated within 450-1500 K. The global rate constants for every pathway are evaluated from the sum of the individual elementary reaction rate constants. The first-order Arrhenius expressions for these six elementary reactions and the three pathways are suggested. By comparing with the experimental data, computational methods without tunneling correction give good description for Path1 (the formation of levoglucosan); while methods with tunneling correction (zero-curvature tunneling and small-curvature tunneling correction) give good results for Path2 (the first possibility for the formation of formaldehyde), all the test methods give similar results for Path3 (the second possibility for the formation of formaldehyde), all the modeling results for Path3 are in good agreement with the experimental data, verifying that it is the most possible way for the formation of formaldehyde during cellulose pyrolysis. © 2012 Elsevier Ltd. All rights reserved.