Numerical analysis of refrigerant flow along non-adiabatic capillary tubes using a two-fluid model

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Limits of validity for a semiclassical mean-field two-fluid model for Bose-Einstein condensation thermodynamics

We reinvestigate the Bose-Einstein condensation (BEC) thermodynamics of a weakly interacting dilute Bose gas under the action of a trap using a semi-classical two-fluid mean-field model in order to find the domain of applicability of the model. Such a model is expected to break down once the condition of diluteness and weak interaction is violated. We find that this breakdown happens for values of coupling and density near the present experimental scenario of BEG. With the increase of the interaction coupling and density the model may lead to unphysical results for thermodynamic observables. (C) 2000 Published by Elsevier B.V. B.V, All rights reserved.

A two-fluid model for refrigerant flow through adiabatic capillary tubes

This work presents a numerical model to simulate refrigerant flow through capillary tubes, commonly used as expansion devices in refrigeration systems. The flow is divided in a single-phase region, where the refrigerant is in the subcooled liquid state, and a region of two-phase flow. The capillary tube is considered straight and horizontal. The flow is taken as one-dimensional and adiabatic. Steady-state condition is also assumed and the metastable flow phenomena are neglected. The two-fluid model, considering the hydrodynamic and thermal non-equilibrium between the liquid and vapor phases, is applied to the two-phase flow region. Comparisons are made with experimental measurements of the mass flow rate and pressure distribution along two capillary tubes working with refrigerant R-134a in different operating conditions. The results indicate that the present model provides a better estimation than the commonly employed homogeneous model. Some computational results referring to the quality, void fraction, velocities, and temperatures of each phase are presented and discussed.

Notes on the two-brane model with variable tension

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A two band model for superconductivity: probing interband pair formation

We propose a two band model for superconductivity. It turns out that the simplest nontrivial case considers solely interband scattering, and both bands can be modeled as symmetric (around the Fermi level) and flat, thus each band is completely characterized by its half-band width Wn (n=1,2). A useful dimensionless parameter is d, proportional to W2 - W1. The case delta = 0 retrieves the conventional BCS model. We probe the specific heat, the ratio gap over critical temperature, the thermodynamic critical field and tunneling conductance as functions of d and temperature (from zero to Tc). We compare our results with experimental results for MgB2 and good quantitative agreement is obtained, indicating the relevance of interband coupling. Work in progress also considers the inclusion of band hybridization and general interband as well as intra-band scattering mechanisms.

Charge-orbital ordering and phase separation in the two-orbital model for manganites: Roles of Jahn-Teller phononic and Coulombic interactions

The main properties of realistic models for manganites are studied using analytic mean-field approximations and computational numerical methods, focusing on the two-orbital model with electrons interacting through Jahn-Teller (JT) phonons and/or Coulombic repulsions. Analyzing the model including both interactions by the combination of the mean-field approximation and the exact diagonalization method, it is argued that the spin-charge-orbital structure in the insulating phase of the purely JT-phononic model with a large Hund couphng J(H) is not qualitatively changed by the inclusion of the Coulomb interactions. As an important application of the present mean-held approximation, the CE-type antiferromagnetic state, the charge-stacked structure along the z axis, and (3x(2) - r(2))/(3y(2) - r(2))-type orbital ordering are successfully reproduced based on the JT-phononic model with large JH for the half-doped manganite, in agreement with recent Monte Carlo simulation results. Topological arguments and the relevance of the Heisenberg exchange among localized t(2g) spins explains why the inclusion of the nearest-neighbor Coulomb interaction does not destroy the charge stacking structure. It is also verified that the phase-separation tendency is observed both in purely JT-phononic (large JH) and purely Coulombic models in the vicinity of the hole undoped region, as long as realistic hopping matrices are used. This highlights the qualitative similarities of both approaches and the relevance of mixed-phase tendencies in the context of manganites. In addition, the rich and complex phase diagram of the two-orbital Coulombic model in one dimension is presented. Our results provide robust evidence that Coulombic and JT-phononic approaches to manganites are not qualitatively different ways to carry out theoretical calculations, but they share a variety of common features.

Análise do escoamento de HFC-134a em tubos capilares usando o modelo de dois fluidos

This work presents a numerical model to simulate refrigerant flow through capillary tubes, commonly used as expansion devices in refrigeration systems. The capillary tube is considered straight and horizontal. The flow is taken as one-dimensional and adiabatic. Steady state and thermodynamic equilibrium conditions are assumed. The two-fluid model, involving four conservation equations and considering the hidrodynamic nonequilibrium between the liquid and vapor phases is applied to the flow region. The pressure profiles and the mass flow rates given by the model are compared with experimental data.

Efeitos numéricos na simulação de escoamentos gás-sólido em leito fuidizado borbulhante utilizando a teoria cinética dos escoamentos granulares

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Numerical simulation of drop impact and jet buckling problems using the eXtended Pom-Pom model

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A Reynolds stress model for turbulent flows of viscoelastic fluids

A second-order closure is developed for predicting turbulent flows of viscoelastic fluids described by a modified generalised Newtonian fluid model incorporating a nonlinear viscosity that depends on a strain-hardening Trouton ratio as a means to handle some of the effects of viscoelasticity upon turbulent flows. Its performance is assessed by comparing its predictions for fully developed turbulent pipe flow with experimental data for four different dilute polymeric solutions and also with two sets of direct numerical simulation data for fluids theoretically described by the finitely extensible nonlinear elastic - Peterlin model. The model is based on a Newtonian Reynolds stress closure to predict Newtonian fluid flows, which incorporates low Reynolds number damping functions to properly deal with wall effects and to provide the capability to handle fluid viscoelasticity more effectively. This new turbulence model was able to capture well the drag reduction of various viscoelastic fluids over a wide range of Reynolds numbers and performed better than previously developed models for the same type of constitutive equation, even if the streamwise and wall-normal turbulence intensities were underpredicted.

A model to predict R134a refrigerant leakage through the radial clearance of rolling piston compressors

In this paper a non-isothermal two-phase model for oil-R134a refrigerant mixture flow is presented to predict the R134a leakage through the radial clearance of rolling piston compressors. The flow is divided in a liquid single-phase region and in a two-phase region, in which the homogeneous model is used to simulate the flow. The refrigerant leakage is determined using the mixture mass flow rate and the refrigerant mass fraction variation along the flow. The results are obtained for inlet pressures varying from 200 to 700 kPa, inlet temperatures ranging from 40 to 60 degrees C, and minimal clearances between 10 and 60 mu m. The results are firstly compared to existing isothermal model data, showing that there is a significant difference between the leakage flow rates predicted by isothermal and non-isothermal models. Finally, a useful general equation for compressor designers is proposed to calculate the refrigerant leakage for a large range of operation conditions. (C) 2012 Elsevier Ltd and IIR. All rights reserved.

Bose-Einstein condensation thermodynamics of a trapped gas with attractive interaction

We study the Bose-Einstein condensation of an interacting gas with attractive interaction confined in a harmonic trap using a semiclassical two-fluid mean-field model. The condensed state is described by the converged numerical solution of the Gross-Pitaevskii equation. By solving the system of coupled equations of this model iteratively we obtain the converged results for the temperature dependencies of the condensate fraction, chemical potential, and internal energy for the Bose-Einstein condensate of Li-7 atoms. (C) 2000 Elsevier B.V. B.V. All rights reserved.

Noncommutative fluid dynamics in the Snyder space-time

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The effect of hydrothermal alteration on the Sr and Nd isotopic signatures of the Barra do Itapirapua Carbonatite, Southern Brazil

The Cretaceous Barra do Itapirapua carbonatite in southern Brazil experienced extensive postmagmatic hydrothermal alteration. In this article, Sr and Nd isotope ratios of coexisting samples of hydrothermally overprinted and of preserved, nonoverprinted carbonatite are presented. Hydrothermal alteration caused strong REE enrichment, leading to the formation of minerals of the bastnaesite group. In the overprinted carbonatite, Nd contents reach 4000 ppm, two orders of magnitude higher than in the fresh carbonatite, but epsilon(Nd) varies only within a range of 3.4 units. In contrast, Sr was leached from the carbonatite during the postmagmatic alteration; hence values of around 10,000 ppm in the fresh carbonatite drop to about 1000 ppm in the overprinted samples. Leaching is accompanied by a variation of Sr isotopic composition toward more radiogenic values, resulting in an increase of 15 units in epsilon(Sr). Variation of Sr isotopic composition is related to postmagmatic alteration and is decoupled from the variation of Nd isotopic composition, ruling out heterogeneities in the mantle source as the main cause of isotopic variability in the data set. Furthermore, this cannot be explained by bulk crustal contamination. A two-step model is proposed in which (1) a REE-rich, carbonatite-derived hydrothermal fluid overprinted the pristine carbonatite, causing REE-enrichment with a relative small change of isotopic composition; and (2) crust-derived hydrothermal fluids percolated the cooling carbonatite, leaching the original Sr from the carbonatite and introducing a more radiogenic Sr isotopic signature. The amounts of carbonatite-derived Nd with primitive, carbonatite-like Nd isotope ratios introduced during the first stage of hydrothermal alteration are high enough to buffer the effect of crust-derived Nd on the Nd isotopic composition of the overprinted carbonatite.

Heat transfer in fixed bed: a model non-linearity approach

Studies have been carried out on the heat transfer in a packed bed of glass beads percolated by air at moderate flow rates. Rigorous statistic analysis of the experimental data was carried out and the traditional two parameter model was used to represent them. The parameters estimated were the effective radial thermal conductivity, k, and the wall coefficient, h, through the least squares method. The results were evaluated as to the boundary bed inlet temperature, T-o, number of terms of the solution series and number of experimental points used in the estimate. Results indicated that a small difference in T-o was sufficient to promote great modifications in the estimated parameters and in the statistical properties of the model. The use of replicas at points of high parametric information of the model improved the results, although analysis of the residuals has resulted in the rejection of this alternative. In order to evaluate cion-linearity of the model, Bates and Watts (1988) curvature measurements and the Box (1971) biases of the coefficients were calculated. The intrinsic curvatures of the model (IN) tend to be concentrated at low bed heights and those due to parameter effects (PE) are spread all over the bed. The Box biases indicated both parameters as responsible for the curvatures PE, h being somewhat more problematic. (C) 2000 Elsevier B.V. Ltd. All rights reserved.