Estudos de técnicas de reinicialização em métodos level-set e suas aplicações em escoamentos bifásicos

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

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)

Numerical modeling of flow through an industrial burner orifice

This paper presents numerical modeling of a turbulent natural gas flow through a non-premixed industrial burner of a slab reheating furnace. The furnace is equipped with diffusion side swirl burners capable of utilizing natural gas or coke oven gas alternatively through the same nozzles. The study is focused on one of the burners of the preheating zone. Computational Fluid Dynamics simulation has been used to predict the burner orifice turbulent flow. Flow rate and pressure at burner upstream were validated by experimental measurements. The outcomes of the numerical modeling are analyzed for the different turbulence models in terms of pressure drop, velocity profiles, and orifice discharge coefficient. The standard, RNG, and Realizable k-epsilon models and Reynolds Stress Model (RSM) have been used. The main purpose of the numerical investigation is to determine the turbulence model that more consistently reproduces the experimental results of the flow through an industrial non-premixed burner orifice. The comparisons between simulations indicate that all the models tested satisfactorily and represent the experimental conditions. However, the Realizable k-epsilon model seems to be the most appropriate turbulence model, since it provides results that are quite similar to the RSM and RNG k-epsilon models, requiring only slightly more computational power than the standard k-epsilon model. (C) 2014 Elsevier Ltd. 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.

Parameterized fast decoupled power flow methods for obtaining the maximum loading point of power systems. Part I. Mathematical modeling

The conventional Newton and fast decoupled power flow (FDPF) methods have been considered inadequate to obtain the maximum loading point of power systems due to ill-conditioning problems at and near this critical point. It is well known that the PV and Q-theta decoupling assumptions of the fast decoupled power flow formulation no longer hold in the vicinity of the critical point. Moreover, the Jacobian matrix of the Newton method becomes singular at this point. However, the maximum loading point can be efficiently computed through parameterization techniques of continuation methods. In this paper it is shown that by using either theta or V as a parameter, the new fast decoupled power flow versions (XB and BX) become adequate for the computation of the maximum loading point only with a few small modifications. The possible use of reactive power injection in a selected PV bus (Q(PV)) as continuation parameter (mu) for the computation of the maximum loading point is also shown. A trivial secant predictor, the modified zero-order polynomial which uses the current solution and a fixed increment in the parameter (V, theta, or mu) as an estimate for the next solution, is used in predictor step. These new versions are compared to each other with the purpose of pointing out their features, as well as the influence of reactive power and transformer tap limits. The results obtained with the new approach for the IEEE test systems (14, 30, 57 and 118 buses) are presented and discussed in the companion paper. The results show that the characteristics of the conventional method are enhanced and the region of convergence around the singular solution is enlarged. In addition, it is shown that parameters can be switched during the tracing process in order to efficiently determine all the PV curve points with few iterations. (C) 2003 Elsevier B.V. All rights reserved.

Mathematical modeling of the ester oil-refrigerant R134A mixture two-phase flow with foam formation through a small diameter tube

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

An exact equation for the free surface of a fluid in a porous medium

We study the problem of the evolution of the free surface of a fluid in a saturated porous medium, bounded from below by a. at impermeable bottom, and described by the Laplace equation with moving-boundary conditions. By making use of a convenient conformal transformation, we show that the solution to this problem is equivalent to the solution of the Laplace equation on a fixed domain, with new variable coefficients, the boundary conditions. We use a kernel of the Laplace equation which allows us to write the Dirichlet-to-Neumann operator, and in this way we are able to find an exact differential-integral equation for the evolution of the free surface in one space dimension. Although not amenable to direct analytical solutions, this equation turns out to allow an easy numerical implementation. We give an explicit illustrative case at the end of the article.

Supercritical fluid extraction from fennel (Foeniculum vulgare): global yield, composition and kinetic data

Supercritical fluid extraction (SFE) from solids has proven to be technically feasible for almost any system; nonetheless, its economical viability has been proven for a restricted number of systems. A common practice is to compare the cost of manufacturing of vegetable extracts by a variety of techniques without deeply considering the huge differences in composition and functional properties among the various types of extracts obtained; under this circumstance, the cost of manufacturing do not favor SFE. Additionally, the influence of external parameters such as the agronomic conditions and the SFE system geometry are not considered. In the present work, these factors were studied for the system fennel seeds + CO2. The effects of the harvesting season and the degree of maturation on the global yields for the system fennel seeds + CO2 were analyzed at 300 bar and 40 degrees C. The effects of the pressure on the global yields were determined for the temperatures of 30 and 40 degrees C. Kinetics experiments were done for various ratios of bed height to bed diameter. Fennel extracts were also obtained by hydrodistillation and low-pressure solvent extraction. The chemical composition of the fennel extracts were determined by gas chromatography. The SFE maximum global yield (12.5%, dry basis) was obtained with dry harvested fennel seeds. Anethole and fenchone were the major constituents of the extract; the following fat acids palmitic (C16H32O2), palmitoleic stearic (C18H36O2), oleic (C18H34O2), linoleic (C18H32O2) and linolenic (C18H30O2) were also detected in the extracts. A relation between amounts of feed and solvent, bed height and diameter, and solvent flow rate was proposed. The models of Sovova, Goto et al. and Tan and Lion were capable of describing the mass transfer kinetics. (c) 2005 Elsevier B.V. All rights reserved.

Study of alternative schemes for the parameterization step of the continuation power flow method based on physical parameters, part I: Mathematical modeling

The conventional power flow method is considered to be inadequate to obtain the maximum loading point because of the singularity of Jacobian matrix. Continuation methods are efficient tools for solving this kind of problem since different parameterization schemes can be used to avoid such ill-conditioning problems. This paper presents the details of new schemes for the parameterization step of the continuation power flow method. The new parameterization options are based on physical parameters, namely, the total power losses (real and reactive), the power at the slack bus (real or reactive), the reactive power at generation buses, and transmission line power losses (real and reactive). The simulation results obtained with the new approach for the IEEE test systems (14, 30, 57, and 118 buses) are presented and discussed in the companion paper. The results show that the characteristics of the conventional method are not only preserved but also improved.

Distribution Transformer Modeling for Application in Three-Phase Power Flow Algorithm

In this paper, it is presented a methodology for three-phase distribution transformer modeling, considering several types of transformer configuration, to be used in algorithms of power flow in three-phase radial distribution networks. The paper provides a detailed discussion about the models and the results from an implementation of the power flow algorithm. The results, taken from three different networks, are presented for several transformer configurations and for voltage regulators as well.

Simulations of incompressible fluid flows by a least squares finite element method

In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using velocity-pressure-vorticity and velocity-pressure-stress formulations, named u-p-ω) and u-p-τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is convenient for implementation by LSFEM. The main purposes of this work are the numerical computation of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier-Stokes equations in u-p-ω and u-p-τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate the numerical code and the preliminary results are presented and compared with available results from the literature. Copyright © 2005 by ABCM.

Voltage regulator modeling for the three-phase power flow in distribution networks

In this paper a three-phase power flow for electrical distribution systems considering different models of voltage regulators is presented. A voltage regulator (VR) is an equipment that maintains the voltage level in a predefined value in a distribution line in spite of the load variations within its nominal power. Three different types of connections are analyzed: 1) wye-connected regulators, 2) open delta-connected regulators and 3) closed delta-connected regulators. To calculate the power flow, the three-phase backward/forward sweep algorithm is used. The methodology is tested on the IEEE 34 bus distribution system. ©2008 IEEE.

Energetic, ecologic and fluid-dynamic analysis of a fluidized bed gasifier operating with sugar cane bagasse

This work aims to study the thermodynamic, ecological and fluid-dynamic aspects of a circulating fluidized bed gasifier using sugar cane bagasse as biomass, in order to estimate a model of its normal operation. In the initial stage was analysed the composition of biomass selected (sugar cane bagasse) and its lower heating value (LHV) was calculated. The energy balance of the gasifier was done, being the volumetric flow of air, synthesis gas and biomass estimated. Also the power produced by this gasifier was theoretically estimated. Then the circulating fluidized bed gasifier was designed for operation with approximately 100 kg/h of processed biomass. Cross-sectional area of the reactor, feeder size, diameter of the exit zone of the gases and minimum height of the expanded bed were selected. Some bed gasifier hydrodynamic factors were also studied. The minimum fluidization velocity, fluidization terminal velocity, and average fluidizing velocity were calculated, in order to understand the fluid-dynamic behaviour of gasification of this fuel. It was obtained a theoretical model that can support a possible prototype of circulating fluidized bed gasifier biomass. Finally, there were studied the ecological aspects of the gasifier, through an overall methodology. Ecological efficiencies were estimated for two scenarios: first considering the carbon cycle and thereafter disregarding the carbon cycle. In both cases, it can be proved the ecological viability of the project. © 2013 Elsevier Ltd. All rights reserved.

Enhancement of optical absorption by modulation of the oxygen flow of TiO2 films deposited by reactive sputtering

Oxygen-deficient TiO2 films with enhanced visible and near-infrared optical absorption have been deposited by reactive sputtering using a planar diode radio frequency magnetron configuration. It is observed that the increase in the absorption coefficient is more effective when the O-2 gas supply is periodically interrupted rather than by a decrease of the partial O-2 gas pressure in the deposition plasma. The optical absorption coefficient at 1.5 eV increases from about 1 x 10(2) cm(-1) to more than 4 x 10(3) cm(-1) as a result of the gas flow discontinuity. A red-shift of similar to 0.24 eV in the optical absorption edge is also observed. High resolution transmission electron microscopy with composition analysis shows that the films present a dense columnar morphology, with estimated mean column width of 40nm. Moreover, the interruptions of the O-2 gas flow do not produce detectable variations in the film composition along its growing direction. X-ray diffraction and micro-Raman experiments indicate the presence of the TiO2 anatase, rutile, and brookite phases. The anatase phase is dominant, with a slight increment of the rutile and brookite phases in films deposited under discontinued O-2 gas flow. The increase of optical absorption in the visible and near-infrared regions has been attributed to a high density of defects in the TiO2 films, which is consistent with density functional theory calculations that place oxygen-related vacancy states in the upper third of the optical bandgap. The electronic structure calculation results, along with the adopted deposition method and experimental data, have been used to propose a mechanism to explain the formation of the observed oxygen-related defects in TiO2 thin films. The observed increase in sub-bandgap absorption and the modeling of the corresponding changes in the electronic structure are potentially useful concerning the optimization of efficiency of the photocatalytic activity and the magnetic doping of TiO2 films. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4724334]

Heat transfer in laminar flow of non-Newtonian fluids in ducts of elliptical section

Laminar forced convection inside tubes of various cross-section shapes is of interest in the design of a low Reynolds number heat exchanger apparatus. Heat transfer to thermally developing, hydrodynamically developed forced convection inside tubes of simple geometries such as a circular tube, parallel plate, or annular duct has been well studied in the literature and documented in various books, but for elliptical duct there are not much work done. The main assumption used in this work is a laminar flow of a power flow inside elliptical tube, under a boundary condition of first kind with constant physical properties and negligible axial heat diffusion (high Peclet number). To solve the thermally developing problem, we use the generalized integral transform technique (GITT), also known as Sturm-Liouville transform. Actually, such an integral transform is a generalization of the finite Fourier transform where the sine and cosine functions are replaced by more general sets of orthogonal functions. The axes are algebraically transformed from the Cartesian coordinate system to the elliptical coordinate system in order to avoid the irregular shape of the elliptical duct wall. The GITT is then applied to transform and solve the problem and to obtain the once unknown temperature field. Afterward, it is possible to compute and present the quantities of practical interest, such as the bulk fluid temperature, the local Nusselt number and the average Nusselt number for various cross-section aspect ratios. (C) 2006 Elsevier. SAS. All rights reserved.