Growth of a Bose-Einstein condensate: a detailed comparison of theory and experiment

We extend the earlier model of condensate growth of Davis et at (Davis M J, Gardiner C W and Ballagh R J 2000 Phys. Rev. A 62 063608) to include the effect of gravity in a magnetic trap. We carry out calculations to model the experiment reported by Kohl et al (Kohl M, Davis M J, Gardiner C W, Hansch T and Esslinger T 2001 Preprint cond-mat/0106642) who study the formation of a rubidium Bose-Einstein condensate for a range of evaporative cooling parameters. We find that, in the regime where our model is valid, the theoretical curves agree with all the experimental data with no fitting parameters. However, for the slowest cooling of the gas the theoretical curve deviates significantly from the experimental curves. It is possible that this discrepancy may be related to the formation of a quasicondensate.

Spatial structure and habitat variation in a grasshopper hybrid zone

A hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae) in northern Spain has been analyzed for variation in morphology and ecology. These species are readily distinguished by the number of stridulatory pegs on the hind femur. Both sexes are fully winged and inhabit disturbed habitats throughout the study area. We develop a maximum-likelihood approach to fitting a two-dimensional cline to geographical variation in quantitative traits and for estimating associations of population mean with local habitat. This method reveals a cline in peg number approximately 30 km south of the Picos de Europa Mountains that shows substantial deviations in population mean compared with the expectations of simple tension zone models. The inclusion of variation in local vegetation in the model explains a significant proportion of the residual variation in peg number, indicating that habitat-genotype associations contribute to the observed spatial pattern. However, this association is weak, and a number of populations continue to show strong deviations in mean even after habitat is included in the final model. These outliers may be the result of long-distance colonization of sites distant from the cline center or may be due to a patchy pattern of initial contact during postglacial expansion. As well as contrasting with the smooth hybrid zones described for Chorthippus parallelus, this situation also contrasts with the mosaic hybrid zones observed in Gryllus crickets and in parts of the hybrid zone between Bombina toad species, where habitat-genotype associations account for substantial amounts of among-site variation.

Turbulent flow around a rotating stepped cylinder

Direct numerical simulation (DNS) of turbulent flow around a rotating cylinder with two backward-facing steps axisymmetrically mounted in the circumferential direction was performed and compared with DNS of plane backward-facing step flow (PBSF) of Le [J. Fluid Mech. 330, 349 (1997)]. The original motivation of this work stemmed from the efforts to design a simple device which can generate flows of high turbulence intensity at low cost for corrosion researchers. It turned out that the current flow shows flow structures quite similar to those of PBSF downstream of the step, even though configurations of the two flows are totally different from one another. The stepped cylinder appears to be a cost-effective tool in the generation of flow structures similar to those of PBSF. (C) 2002 American Institute of Physics.

Growth and compaction behaviour of copper concentrate granules in a rotating drum

In order to understand the growth and compaction behaviour of chalcopyrite (copper concentrate), batch granulation tests were carried out using a rotating drum. The granule growth exhibited induction-type behaviour, as defined by Iveson and Litster [AIChE J. 44 (1998) 15 10]. There were two consecutive stages during granulation: the induction stage, during which the granules are gradually being compacted and little or no growth occurs, and the rapid growth stage, which starts when the granules have become surface wet and are rapidly growing. In agreement with earlier findings. an increased amount of binder liquid shortened the induction time. The compaction behaviour was also investigated. A displaced volume method was adopted to determine the porosity of the granules. It was shown that this technique had a limitation as it was unable to detect the reduction of the volumes of the granule pores after the granules had become surface wet. Due to this, some of the measurements were not suited for fitting a three-parameter empirical model. Attempts were made to determine whether the rapid growth stage started with the pore saturation exceeding a certain critical value, but due to the scatter in the porosity measurements and the fact that some of the measurements could not be used, it was not possible to determine a critical pore saturation, However, the porosity measurements clearly demonstrated that the porosity of the granules decreased during the induction stage of an experiment and that when rapid growth occurred, the granules had a pore saturation was around 0.85. This value was slightly lower than unity, which is most likely due to trapped air bubbles. (C) 2002 Published by Elsevier Science B.V.

Surface alloying of AZ91D alloy by diffusion coating

A new technique of surface modification by diffusion coating for AZ91D alloy was developed. A 1.0-2.0-mm alloy layer, which has hardness four to five times higher than the substrate metal, was formed after the treatment. Consequent solution treatment and aging could further improve the hardness of the alloy layer. Microstructure and chemical composition were investigated using optical microscope and electron probe.

Vacancy solution theory for binary adsorption equilibria in heterogeneous carbon

A heterogeneous modified vacancy solution model of adsorption developed is evaluated. The new model considers the adsorption process through a mass-action law and is thermodynamically consistent, while maintaining the simplicity in calculation of multicomponent adsorption equilibria, as in the original vacancy solution theory. It incorporates the adsorbent heterogeneity through a pore-width-related potential energy, represented by Steele's 10-4-3 potential expression. The experimental data of various hydrocarbons, CO2 and SO2 on four different activated carbons - Ajax, Norit, Nuxit, and BPL - at multiple temperatures over a wide range of pressures were studied by the heterogeneous modified VST model to obtain the isotherm parameters and micropore-size distribution of carbons. The model successfully correlates the single-component adsorption equilibrium data for all compounds studied on various carbons. The fitting results for the vacancy occupancy parameter are consistent with the pressure change on different carbons, and the effect of pore heterogeneity is important in adsorption at elevated pressure. It predicts binary adsorption equilibria better than the IAST scheme, reflecting the significance of molecular size nonideality.

Modeling of phase equilibrium and vapor adsorption on carbon black based on a combination of a lattice theory and equation of state

A thermodynamic approach is developed in this paper to describe the behavior of a subcritical fluid in the neighborhood of vapor-liquid interface and close to a graphite surface. The fluid is modeled as a system of parallel molecular layers. The Helmholtz free energy of the fluid is expressed as the sum of the intrinsic Helmholtz free energies of separate layers and the potential energy of their mutual interactions calculated by the 10-4 potential. This Helmholtz free energy is described by an equation of state (such as the Bender or Peng-Robinson equation), which allows us a convenient means to obtain the intrinsic Helmholtz free energy of each molecular layer as a function of its two-dimensional density. All molecular layers of the bulk fluid are in mechanical equilibrium corresponding to the minimum of the total potential energy. In the case of adsorption the external potential exerted by the graphite layers is added to the free energy. The state of the interface zone between the liquid and the vapor phases or the state of the adsorbed phase is determined by the minimum of the grand potential. In the case of phase equilibrium the approach leads to the distribution of density and pressure over the transition zone. The interrelation between the collision diameter and the potential well depth was determined by the surface tension. It was shown that the distance between neighboring molecular layers substantially changes in the vapor-liquid transition zone and in the adsorbed phase with loading. The approach is considered in this paper for the case of adsorption of argon and nitrogen on carbon black. In both cases an excellent agreement with the experimental data was achieved without additional assumptions and fitting parameters, except for the fluid-solid potential well depth. The approach has far-reaching consequences and can be readily extended to the model of adsorption in slit pores of carbonaceous materials and to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Analysis of adsorption data of graphitized thermal carbon black with a DFT-lattice gas theory

In this paper we analyzed the adsorption of gases and vapors on graphitised thermal carbon black by using a modified DFT-lattice theory, in which we assume that the behavior of the first layer in the adsorption film is different from those of second and higher layers. The effects of various parameters on the topology of the adsorption isotherm were first investigated, and the model was then applied in the analysis of adsorption data of numerous substances on carbon black. We have found that the first layer in the adsorption film behaves differently from the second and higher layers in such a way that the adsorbate-adsorbate interaction energy in the first layer is less than that of second and higher layers, and the same is observed for the partition function. Furthermore, the adsorbate-adsorbate and adsorbate-adsorbent interaction energies obtained from the fitting are consistently lower than the corresponding values obtained from the viscosity data and calculated from the Lorentz-Berthelot rule, respectively.

Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

Zygosity Diagnosis in the Absence of Genotypic Data: An Approach Using Latent Class Analysis

For zygosity diagnosis in the absence of genotypic data, or in the recruitment phase of a twin study where only single twins from same-sex pairs are being screened, or to provide a test for sample duplication leading to the false identification of a dizygotic pair as monozygotic, the appropriate analysis of respondents' answers to questions about zygosity is critical. Using data from a young adult Australian twin cohort (N = 2094 complete pairs and 519 singleton twins from same-sex pairs with complete responses to all zygosity items), we show that application of latent class analysis (LCA), fitting a 2-class model, yields results that show good concordance with traditional methods of zygosity diagnosis, but with certain important advantages. These include the ability, in many cases, to assign zygosity with specified probability on the basis of responses of a single informant (advantageous when one zygosity type is being oversampled); and the ability to quantify the probability of misassignment of zygosity, allowing prioritization of cases for genotyping as well as identification of cases of probable laboratory error. Out of 242 twins (from 121 like-sex pairs) where genotypic data were available for zygosity confirmation, only a single case was identified of incorrect zygosity assignment by the latent class algorithm. Zygosity assignment for that single case was identified by the LCA as uncertain (probability of being a monozygotic twin only 76%), and the co-twin's responses clearly identified the pair as dizygotic (probability of being dizygotic 100%). In the absence of genotypic data, or as a safeguard against sample duplication, application of LCA for zygosity assignment or confirmation is strongly recommended.

Direction of Causation Modeling Between Cross-Sectional Measures of Parenting and Psychological Distress in Female Twins

Under certain conditions, cross-sectional analysis of cross-twin intertrait correlations can provide important information about the direction of causation (DOC) between two variables. A community-based sample of Australian female twins aged 18 to 45 years was mailed an extensive Health and Lifestyle Questionnaire (HLQ) that covered a wide range of personality and behavioral measures. Included were self-report measures of recent psychological distress and perceived childhood environment (PBI). Factor analysis of the PBI yielded three interpretable dimensions: Coldness, Overprotection, and Autonomy. Univariate analysis revealed that parental Overprotection and Autonomy were best explained by additive genetic, shared, and nonshared environmental effects (ACE), whereas the best-fitting model for PBI Coldness and the three measures of psychological distress (Depression, Phobic Anxiety, and Somatic Distress) included only additive genetic and nonshared environmental effects (AE). A common pathway model best explained the covariation between (1) the three PBI dimensions and (2) the three measures of psychological distress. DOC modeling between latent constructs of parenting and psychological distress revealed that a model which specified recollected parental behavior as the cause of psychological distress provided a better fit than a model which specified psychological distress as the cause of recollected parental behavior. Power analyses and limitations of the findings are discussed.

Conformational changes in SP-B as a function of surface pressure

X-ray reflectivity of bovine and sheep surfactant-associated protein B (SP-B) monolayers is used in conjunction with pressure-area isotherms and protein models to suggest that the protein undergoes changes in its tertiary structure at the air/water interface under the influence of surface pressure, indicating the likely importance of such changes to the phenomena of protein squeeze out as well as lipid exchange between the air-water interface and subphase structures. We describe an algorithm based on the well-established box- or layer-models that greatly assists the fitting of such unknown scattering-length density profiles, and which takes the available instrumental resolution into account. Scattering-length density profiles from neutron reflectivity of bovine SP-B monolayers on aqueous subphases are shown to be consistent with the exchange of a large number of labile protons as well as the inclusion of a significant amount of water, which is partly squeezed out of the protein monolayer at elevated surface pressures.

Prediction of multilayer adsorption and capillary condensation phenomena in cylindrical mesopores

MCM-41 periodic mesoporous silicates with a high degree of structural ordering are synthesized and used as model adsorbents to study the isotherm prediction of nitrogen adsorption. The nitrogen adsorption isotherm at 77 K for a macroporous silica is measured and used in high-resolution alpha(s)-plot comparative analysis to determine the external surface area, total surface area and primary mesopore volume of the MCM-41 materials. Adsorption equilibrium data of nitrogen on the different pore size MCM-41 samples (pore diameters from 2.40 to 4.92 nm) are also obtained. Based on the Broekhoff and de Boer' thermodynamic analysis, the nitrogen adsorption isotherms for the different pore size MCM-41 samples are interpreted using a novel strategy, in which the parameters of an empirical expression, used to represent the potential of interaction between the adsorbate and adsorbent, are obtained by fitting only the multilayer region prior to capillary condensation for C-16 MCM-41. Subsequently the entire isotherm, including the phase transition, is predicted for all the different pore size MCM-41 samples without any fitting. The results show that the prediction of multilayer adsorption and total adsorbed amount are in good agreement with the experimental isotherms. The predictions of the relative pressure corresponding to capillary equilibrium (coexistence) transition agree remarkably with experimental data on the adsorption branch even for hysteretic isotherms, confirming that this is the branch appropriate for pore size distribution analysis. The impact of pore radius on the adsorption film thickness and capillary coexistence pressure is also investigated, and found to agree with the experimental data. (C) 2003 Elsevier Inc. All rights reserved.

Estudo da corrosão dos aços AISI 1010 e AISI 316L em soluções com elevado teor de cloreto e em ácido naftênico

A mudança no mercado global do petróleo nos últimos anos, com o declínio das reservas de óleo leve, têm forçado a busca por novos campos petrolíferos em ambientes mais remotos, como nos campos localizados na camada pré-sal, e a exploração de óleos pesados que possuem elevado teor de ácidos naftênicos. Isso acarreta em grandes desafios para a previsão do desempenho de materiais frente às novas condições ambientais em que estão inseridos. No presente trabalho, o comportamento da corrosão do aço carbono AISI 1010 e do aço inoxidável AISI 316L foi estudado em soluções aquosas com elevado teor de cloreto e em solução de ácido naftênico ciclopentanóico a fim de ter melhor entendimento da ação dessas espécies no processo de corrosão e simular a corrosão pela água de produção na indústria petrolífera. Foram aplicadas as técnicas de potencial de circuito aberto, polarização potenciodinâmica, voltametria cíclica, espectroscopia de impedância eletroquímica, espectroscopia Raman, microscopia eletrônica de varredura e microscopia de força atômica, usadas, em cada caso, de acordo com a conveniência. O aumento da [Cl-] na faixa de 1,2–2,8 mol.L-1 não altera os processos catódicos e anódicos perto do Ecorr para os aços AISI 1010 e AISI 316L. Em condições de sobrepotenciais afastados do Ecorr, o aumento de [Cl-] aumenta os processos oxidativos de corrosão, o que é expresso pelas maiores densidades de corrente e carga anódica e aumento da perda de massa sofrida pelos eletrodos de ambos os aços. Portanto, os danos da corrosão são mais intensos quando se aumenta a [Cl-]. O aço AISI 1010 é ativo nas soluções de NaCl e a corrosão se propaga livremente de forma uniforme. Para o aço AISI 316L, uma ampla faixa de passividade pode ser vista nas soluções de NaCl; no Epit ocorre a ruptura do filme passivo e o crescimento de pites estáveis. Após 24 h de imersão em soluções de sulfato de sódio (branco) e de ácido naftênico ciclopentanóico ocorre crescimento de filme de óxido e as fases α-Fe2O3, Fe3O4 e δ-FeO(OH) foram identificadas nos espécimes de aço AISI 1010 e Fe3O4 foi identificado nos defeitos do filme prévio presente na superfície do aço AISI 316L. Os filmes formados em solução de ácido ciclopentanóico possuem menor resistência à polarização, maior rugosidade e maior taxa de corrosão quando comparado aos filmes crescidos na solução branco, para ambos os aços. A presença do ácido naftênico muda a forma como a reação de corrosão se procede e contribui para o aumento da corrosão. A corrosão naftênica foi mais pronunciada no aço carbono porque a presença dos elementos de liga no aço inox reduzem o número de sítios ativos ricos em Fe e tornam menos oportuna a ligação do Fe com o naftenato.

Estudo da corrosão e inibição dos aços carbono AISI 1010, inoxidável AISI 316 e duplex UNS S31803 em meio de solução de íons cloreto

O comportamento da corrosão e inibição à corrosão dos aços carbono AISI 1010, inox AISI 316 e duplex UNS S31803 foi estudado em meio de solução de íons cloreto à 3,0% (m/v), na ausência e presença do benzimidazol e imidazol como inibidores. A caracterização química e morfológica dos aços foi realizada por meio das técnicas de espectrometria de emissão ótica, difração de raios X (DRX), microscopia ótica, microscopia eletrônica de varredura (MEV) e energia dispersiva de raios X (EDX). As análises eletroquímicas foram realizadas através das técnicas de polarização potenciodinâmica e espectroscopia de impedância eletroquímica. As análises de DRX e de metalografia mostraram as fases presentes em cada aço, sendo o aço AISI 1010 composto pela fase ferrita, o aço AISI 316 pelas fases de FeNi e Cr e o aço UNS S31803 pelas fases austenita e ferrita. Além disso, a metalografia e as análises de MEV e EDX permitiram identificar regiões e certos elementos presentes nos aços que propiciam à ocorrência da corrosão, tais como inclusões. Os inibidores foram testados em diferentes concentrações (25 ppm, 50 ppm, 100 ppm, 500 ppm e 1000 ppm) para os três aços, através das curvas de polarização e impedância eletroquímica, e verificou-se que para todas as concentrações houve aumento da resistência à corrosão dos aços. Pelas curvas de polarização verificou-se que o benzimidazol proporcionou aos aços AISI 1010, AISI 316 e UNS S31803, eficiências de inibição de cerca de 51%, 71% e 75%, respectivamente. Enquanto que o imidazol apresentou eficiência de cerca de 73%, 95% e 86%, respectivamente. Os resultados de impedância eletroquímica mostraram que as eficiências de inibição do benzimidazol foram de aproximadamente 52%, 73% e 71%, respectivamente, para os aços AISI 1010, AISI 316 e UNS S31803. E por sua vez, o imidazol apresentou eficiências de aproximadamente 96% para os aços AISI 1010 e AISI 316 e 85% para o aço UNS S31803. O teste de perda de massa mostrou que para o aço AISI 1010 tanto o benzimidazol quanto e o imidazol inibiram a corrosão, sendo que reduziram a corrosão em cerca de 17% e 24%, respectivamente. Nas análises das curvas de polarização em estudos com a água do mar observou-se que os inibidores foram menos eficientes do que em meio de solução de cloreto. O benzimidazol obteve eficiências de cerca de 14%, 50% e 33%, respectivamente, para os aços AISI 1010, AISI 316 e UNS S31803. Enquanto que o imidazol apresentou eficiências de aproximadamente 21%, 59% e 34%, respectivamente. Em todas as análises eletroquímicas e análise de perda de massa, o imidazol se mostrou o melhor inibidor para os aços estudados.