A semi-Lagrangian finite volume method for Newtonian contraction flows

A new finite volume method for solving the incompressible Navier--Stokes equations is presented. The main features of this method are the location of the velocity components and pressure on different staggered grids and a semi-Lagrangian method for the treatment of convection. An interpolation procedure based on area-weighting is used for the convection part of the computation. The method is applied to flow through a constricted channel, and results are obtained for Reynolds numbers, based on half the flow rate, up to 1000. The behavior of the vortex in the salient corner is investigated qualitatively and quantitatively, and excellent agreement is found with the numerical results of Dennis and Smith [Proc. Roy. Soc. London A, 372 (1980), pp. 393-414] and the asymptotic theory of Smith [J. Fluid Mech., 90 (1979), pp. 725-754].

Moment decay rates of solutions of stochastic differential equations

The objective of this paper is to investigate the p-ίh moment asymptotic stability decay rates for certain finite-dimensional Itό stochastic differential equations. Motivated by some practical examples, the point of our analysis is a special consideration of general decay speeds, which contain as a special case the usual exponential or polynomial type one, to meet various situations. Sufficient conditions for stochastic differential equations (with variable delays or not) are obtained to ensure their asymptotic properties. Several examples are studied to illustrate our theory.

Uniqueness criteria for continuous-time Markov chains with general transition structures

We derive necessary and sufficient conditions for the existence of bounded or summable solutions to systems of linear equations associated with Markov chains. This substantially extends a famous result of G. E. H. Reuter, which provides a convenient means of checking various uniqueness criteria for birth-death processes. Our result allows chains with much more general transition structures to be accommodated. One application is to give a new proof of an important result of M. F. Chen concerning upwardly skip-free processes. We then use our generalization of Reuter's lemma to prove new results for downwardly skip-free chains, such as the Markov branching process and several of its many generalizations. This permits us to establish uniqueness criteria for several models, including the general birth, death, and catastrophe process, extended branching processes, and asymptotic birth-death processes, the latter being neither upwardly skip-free nor downwardly skip-free.

Differences in performance among four indices used to evaluate diversity in planktonic ecosystems

The performance of four common estimators of diversity are investigated using calanoid copepod data from the Continuous Plankton Recorder (CPR) survey. The region of the North Atlantic and the North Sea was divided into squares of 400 nautical miles for each 2-month period. For each 144 possible cases, Pielou's pooled quadrat method was performed with the aims of determining asymptotic diversity and investigating the CPR sample-size dependence of diversity estimators. It is shown that the performance of diversity indices may greatly vary in space and time (at a seasonal scale). This dependence is more pronounced in higher diverse environments and when the sample size is small. Despite results showing that all estimators underestimate the `actual' diversity, comparison of sites remained reliable from a few pooled CPR samples. Using more than one CPR sample, the Gini coefficient appears to be a better diversity estimator than any other indices and spatial or temporal comparisons are highly satisfactory. In situations where comparative studies are needed but only one CPR sample is available, taxonomic richness was the preferred method of estimating diversity. Recommendations are proposed to maximise the efficiency of diversity estimations with the CPR data.

Towards an integrated forecasting system for pelagic fisheries

First results of a coupled modeling and forecasting system for the pelagic fisheries are being presented. The system consists currently of three mathematically fundamentally different model subsystems: POLCOMS-ERSEM providing the physical-biogeochemical environment implemented in the domain of the North-West European shelf and the SPAM model which describes sandeel stocks in the North Sea. The third component, the SLAM model, connects POLCOMS-ERSEM and SPAM by computing the physical-biological interaction. Our major experience by the coupling model subsystems is that well-defined and generic model interfaces are very important for a successful and extendable coupled model framework. The integrated approach, simulating ecosystem dynamics from physics to fish, allows for analysis of the pathways in the ecosystem to investigate the propagation of changes in the ocean climate and lower trophic levels to quantify the impacts on the higher trophic level, in this case the sandeel population, demonstrated here on the base of hindcast data. The coupled forecasting system is tested for some typical scientific questions appearing in spatial fish stock management and marine spatial planning, including determination of local and basin scale maximum sustainable yield, stock connectivity and source/sink structure. Our presented simulations indicate that sandeels stocks are currently exploited close to the maximum sustainable yield, but large uncertainty is associated with determining stock maximum sustainable yield due to stock eigen dynamics and climatic variability. Our statistical ensemble simulations indicates that the predictive horizon set by climate interannual variability is 2–6 yr, after which only an asymptotic probability distribution of stock properties, like biomass, are predictable.

A multicomponent model of phytoplankton size structure

Size-fractionated filtration (SFF) is a direct method for estimating pigment concentration in various size classes. It is also common practice to infer the size structure of phytoplankton communities from diagnostic pigments estimated by high-performance liquid chromatography (HPLC). In this paper, the three-component model of Brewin et al. (2010) was fitted to coincident data from HPLC and from SFF collected along Atlantic Meridional Transect cruises. The model accounted for the variability in each data set, but the fitted model parameters differed for the two data sets. Both HPLC and SFF data supported the conceptual framework of the three-component model, which assumes that the chlorophyll concentration in small cells increases to an asymptotic maximum, beyond which further increase in chlorophyll is achieved by the addition of larger celled phytoplankton. The three-component model was extended to a multicomponent model of size structure using observed relationships between model parameters and assuming that the asymptotic concentration that can be reached by cells increased linearly with increase in the upper bound on the cell size. The multicomponent model was verified using independent SFF data for a variety of size fractions and found to perform well (0.628 ≤ r ≤ 0.989) lending support for the underlying assumptions. An advantage of the multicomponent model over the three-component model is that, for the same number of parameters, it can be applied to any size range in a continuous fashion. The multicomponent model provides a useful tool for studying the distribution of phytoplankton size structure at large scales.

Proof of polar ejection from the close-binary core of the planetary nebula Abell 63

We present the first detailed kinematical analysis of the planetary nebula Abell 63, which is known to contain the eclipsing close-binary nucleus UU Sge. Abell 63 provides an important test case in investigating the role of close-binary central stars on the evolution of planetary nebulae. Longslit observations were obtained using the Manchester echelle spectrometer combined with the 2.1-m San Pedro Martir Telescope. The spectra reveal that the central bright rim of Abell 63 has a tube-like structure. A deep image shows collimated lobes extending from the nebula, which are shown to be high-velocity outflows. The kinematic ages of the nebular rim and the extended lobes are calculated to be 8400 +/- 500 and 12900 +/- 2800 yr, respectively, which suggests that the lobes were formed at an earlier stage than the nebular rim. This is consistent with expectations that disc-generated jets form immediately after the common envelope phase. A morphological-kinematical model of the central nebula is presented and the best-fitting model is found to have the same inclination as the orbital plane of the central binary system; this is the first proof that a close-binary system directly affects the shaping of its nebula. A Hubble-type flow is well-established in the morphological-kinematical modelling of the observed line profiles and imagery. Two possible formation models for the elongated lobes of Abell 63 are considered, (i) a low-density, pressure-driven jet excavates a cavity in the remnant asymptotic giant branch (AGB) envelope; (ii) high-density bullets form the lobes in a single ballistic ejection event.

Early-type stars in the Galactic halo from the Palomar-Green survey - III. Completion of a magnitude range limited sample

High-resolution (R approximate to 40 000) echelle spectroscopic observations of 13 high-latitude early-type stars are presented. These stars comprise the final part of a complete magnitude range limited sample based on low-resolution spectroscopy of targets drawn from the Palomar-Green survey. The magnitude range under consideration is 13 less than or equal to B-PG less than or equal to 14.6, corresponding to an approximate distance limit for main-sequence B-type objects of 5 less than or equal to d less than or equal to 40 kpc. Three stars are found to be apparently normal, young stars, based on their positions on the (T-eff, log g) diagram, normal abundance patterns and relatively large projected rotational velocities. A further star, PG 1209+263, was found to belong to the chemically peculiar (CP) silicon star class of objects. The remainder are evolved subluminous stars lying on post- horizontal branch (post-HB) tracks, with the exception of PG 2120+062, which appears to be in a post-asymptotic giant branch evolutionary stage. For the young stars in the sample, we have derived distance and age estimates through comparison of the atmospheric parameters with recent theoretical evolutionary models. We discuss formation scenarios by comparing times-of- flight and evolutionary time-scales. It is found that all stars could have formed in the Galactic disc and been ejected from there soon after their birth, with the exception of PG 1209+263. The adopted proper motion is found to be a crucial factor in the kinematical analysis. We also present some number densities for young B-type halo stars, which indicate that they are extremely scarce objects.

A search for circumstellar material around B-type stars in the Galactic halo

19 B-type stars, selected from the Palomar-Green Survey, have been observed at infrared wavelengths to search for possible infrared excesses, as part of an ongoing programme to investigate the nature of early-type stars at high Galactic latitudes. The resulting infrared fluxes, along with Stromgren photometry, are compared with theoretical flux profiles to determine whether any of the targets show evidence of circumstellar material, which may be indicative of post-main- sequence evolution. Eighteen of the targets have flux distributions in good agreement with theoretical predictions. However, one star, PG 2120 + 062, shows a small near-infrared excess, which may be due either to a cool companion of spectral type F5-F7, or to circumstellar material, indicating that it may be an evolved object such as a post-asymptotic giant branch star, in the transition region between the asymptotic giant branch and planetary nebula phase, with the infrared excess due to recent mass loss during giant branch evolution.

A detailed abundance analysis of the hot post-AGB star ZNG-1 in M10

We present a model-atmosphere analysis for the bright (V similar to 13) star ZNG-1, in the globular cluster M10. From high-resolution (R similar to 40 000) optical spectra we confirm ZNG-1 to be a post-asymptotic giant branch (post-AGB) star. The derived atmospheric parameters are T-eff = 26 500 +/- 1000 K and log g = 3.6 +/- 0.2 dex. A differential abundance analysis reveals a chemical composition typical of hot post-AGB objects, with ZNG-1 being generally metal poor, although helium is approximately solar. The most interesting feature is the large carbon underabundance of more than 1.3 dex. This carbon deficiency, along with an observed nitrogen enhancement relative to other elements, may suggest that ZNG-1 evolved off the AGB before the third dredge-up occurred. Also, iron depletions observed in other similar stars suggest that gas- dust fractionation in the AGB progenitor could be responsible for the observed composition of these objects. However, we need not invoke either scenario since the chemical composition of ZNG-1 is in good agreement with abundances found for a Population II star of the same metallicity.

Boron abundances in B-type stars: A test of rotational depletion during main-sequence evolution

Boron abundances have been derived for seven main-sequence B- type stars from Hubble Space Telescope STIS spectra around the B III lambda2066 line. In two stars, boron appears to be undepleted with respect to the presumed initial abundance. In one star, boron is detectable but is clearly depleted. In the other four stars, boron is undetectable, implying depletions of 1-2 dex. Three of these four stars are nitrogen enriched, but the fourth shows no enrichment of nitrogen. Only rotationally induced mixing predicts that boron depletions are unaccompanied by nitrogen enrichments. The inferred rate of boron depletion from our observations is in good agreement with these predictions. Other boron-depleted nitrogen-normal stars are identified from the literature. In addition, several boron- depleted nitrogen-rich stars are identified, and while all fall on the boron-nitrogen trend predicted by rotationally induced mixing, a majority have nitrogen enrichments that are not uniquely explained by rotation. The spectra have also been used to determine iron group (Cr, Mn, Fe, and Ni) abundances. The seven B-type stars have near-solar iron group abundances, as expected for young stars in the solar neighborhood. We have also analyzed the halo B-type star PG 0832 + 676. We find [Fe/H] = -0.88 +/- 0.10, and the absence of the B III line gives the upper limit [B/H] <-2.5. These and other published abundances are used to infer the star's evolutionary status as a post-asymptotic giant branch star.

Plane wave and Coulomb asymptotics

A simple plane wave solution of the Schrodinger-Helmholtz equation is a quantum eigenfunction obeying both energy and linear momentum correspondence principles. Inclusion of the outgoing wave with scattering amplitude f asymptotic development of the plane wave, we show that there is a problem with angular momentum when we consider forward scattering at the point of closest approach and at large impact parameter given semiclassically by (l + 1/2)/k where l is the azimuthal quantum number and may be large (J. Leech et al., Phys. Rev. Lett. 88. 257901 (2002)). The problem is resolved via non- uniform, non-standard analysis involving the Heaviside step function, unifying classical, semiclassical and quantum mechanics, and the treatment is extended to the case of pure Coulomb scattering.

Semiclassical treatment of Wannier's theory when the exponent diverges

We consider a non-standard application of the Wannier model. A physical example is the single ionization of a hydrogenic beryllium ion with a fully stripped beryllium ion, where the ratio of the charge of the third particle to the charges of the escaping particles is 1/4; we investigate the single ionization by an electron of an atom comprising an electron and a nucleus of charge 1/4. An infinite exponent is obtained suggesting that this process is not tractable within the Wannier model. A modified version of Crothers' uniform semiclassical wavefunction for the outgoing particles has been adopted, since the Wannier exponents and are infinite for an effective charge of Z = 1/4. We use Bessel functions to describe the Peterkop functions u and u and derive a new turning point ?. Since u is well behaved at infinity, there exists only the singularity in u at infinity, thus we employ a one- (rather than two-) dimensional change of dependent variable, ensuring that a uniform solution is obtained that avoids semiclassical breakdown on the Wannier ridge. The regularized final-state asymptotic wavefunction is employed, along with a continuum-distorted-wave approximation for the initial-state wavefunction to obtain total cross sections on an absolute scale. © 2006 IOP Publishing Ltd.

A chemical analysis of five hot stars towards the Galactic centre

High resolution echelle spectroscopy is presented for thirteen stars lying in the direction of the Galactic centre which, on the basis of photographic photometry and low dispersion spectroscopy, have been classified as early-B-type. Eight of these stars have large rotational velocities which preclude a detailed analysis. The five stars with moderate to low projected rotational velocities have been analysed using model atmosphere techniques to determine atmospheric parameters and chemical compositions. Two of these stars appear to be evolved blue horizontal branch objects on the basis of their chemical compositions and small projected rotational velocity. The evolutionary status of a third is ambiguous but it is probably a post-asymptotic-giant branch star. The remaining two objects are probably young massive stars and show enhanced abundances of N, C, Mg and Si, consistent with their formation in the inner part of the Galactic disk. However their O abundances are normal, confirming results found previously for other early- type stars, which would imply a flat abundance gradient for this element in the inner region of our Galaxy.