Simulation results and applications of an advection bounded scheme to practical flows

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

A new approach for thermal performance calculation of cross-flow heat exchangers

A new numerical methodology for thermal performance calculation in cross-flow heat exchangers is developed. Effectiveness-number of transfer units (epsilon-NTU) data for several standard and complex flow arrangements are obtained using this methodology. The results are validated through comparison with analytical solutions for one-pass cross-flow heat exchangers with one to four rows and with approximate series solution for an unmixed-unmixed heat exchanger, obtaining in all cases very small errors. New effectiveness data for some complex configurations are provided. (c) 2005 Elsevier Ltd. All rights reserved.

Numerical simulation of a radial diffuser turbulent airflow

In the present work are presented results from numerical simulations performed with the ANSYS-CFX (R) code. We have studied a radial diffuser flow case, which is the main academic problem used to study the flow behavior on flat plate valves. The radial flow inside the diffuser has important behavior such as the turbulence decay downstream and recirculation regions inside the valve flow channel due to boundary layer detachment. These flow structures are present in compressor reed valve configurations, influencing to a greater extent the compressor efficiency. The main target of the present paper was finding the simulation set-up (computational domain, boundary conditions and turbulence model) that better fits with experimental data published by Tabatabai and Pollard. The local flow turbulence and velocity profiles were investigated using four different turbulence models, two different boundary conditions set-up, two different computational domains and three different flow conditions (Re-in - Reynolds number at the diffuser inlet). We used the Reynolds stress (BSL); the k-epsilon; the RNG k-epsilon; and the shear stress transport (SST) k-omega turbulence models. The performed analysis and comparison of the computational results with experimental data show that the choice of the turbulence model, as well as the choice of the other computational conditions, plays an important role in the results physical quality and accuracy. (c) 2007 Elsevier B.V. All rights reserved.

Simulation and economic analysis of bovine sex selection

A simulation model implemented in the programming software Delphi XE® was applied to evaluate sex selection in bovine. The hypothesis under investigation was that a dynamic model with stochastic and deterministic elements could detect the sexed semen technique to minimize pregnancy cost and to determine the adequate number of recipients required for in vivo (ET) and in vitro embryo production (IVP) in the proposed scenarios. Sex selection was compared through semen sexed using flow cytometry (C1) and density gradient centrifugation techniques (C2) in ET and IVP. Sensibility analyses were used to identify the adequate number of recipients for each scenario. This number was reinserted into the model to determine the biological and financial values that maximized ET and IVP using sexed semen (C1M and C2M). New scenarios showed that the density gradient technique minimized pregnancy cost based on the proposed scenarios. In addition, the adequate number of recipients (ET - C1M - 115 and C2M - 105)/(IVP - C1M - 145 and C2M - 140) per donor used was determined to minimize the pregnancy cost in all scenarios.

Continuation load flow using automatically determined branch megawatt losses as parameters

The conventional Newton's method is considered to be inadequate for the computation of the maximum loading point (MLP) of power systems since: (i) it encounters difficulties in the vicinity of the MLP: and (ii) the load flow Jacobian matrix becomes singular at the MLP. It is well known that continuation methods are powerful and useful tools that are able to trace the solution PV curve without experiencing such diffculties. However, continuation methods require a parameterisation so that a modified, well conditioned set of load flow equations is obtained. In particular, the Jacobian matrix associated with this modified set of equations should not be singular at the MLP. The authors propose that the actual power losses in transmission branches (lines and transformers) are used to parameterise the approach. Specific procedures for the automatic determination of the most appropriate parameter (branch) are proposed. Such procedures include the utilisation of fast voltage-stability indices. Simulation results are presented to show that the proposed method is able to trace the whole solution PV curve very efficiently.

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.

The use of simulation for forecasting the number of docks necessary in an expedition department

The simulation is a very powerful tool to develop more efficient systems, hence it is been widely used with the goal of productivity improvement. Its results, if compared with other methods, are not always optimum; however, if the experiment is rightly elaborated, its results will represent the real situation, enabling its use with a good level of reliability. This work used the simulation (through the ProModel (R) software) in order to study, understand, model and improve the expenditure system of an enterprise, with a premise of keeping the production-delivery flow considering quick, controlled and reliable conditions.

The STATFLUX code: a statistical method for calculation of flow and set of parameters, based on the Multiple-Compartment Biokinetical Model

The code STATFLUX, implementing a new and simple statistical procedure for the calculation of transfer coefficients in radionuclide transport to animals and plants, is proposed. The method is based on the general multiple-compartment model, which uses a system of linear equations involving geometrical volume considerations. Flow parameters were estimated by employing two different least-squares procedures: Derivative and Gauss-Marquardt methods, with the available experimental data of radionuclide concentrations as the input functions of time. The solution of the inverse problem, which relates a given set of flow parameter with the time evolution of concentration functions, is achieved via a Monte Carlo Simulation procedure.Program summaryTitle of program: STATFLUXCatalogue identifier: ADYS_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADYS_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputer for which the program is designed and others on which it has been tested: Micro-computer with Intel Pentium III, 3.0 GHzInstallation: Laboratory of Linear Accelerator, Department of Experimental Physics, University of São Paulo, BrazilOperating system: Windows 2000 and Windows XPProgramming language used: Fortran-77 as implemented in Microsoft Fortran 4.0. NOTE: Microsoft Fortran includes non-standard features which are used in this program. Standard Fortran compilers such as, g77, f77, ifort and NAG95, are not able to compile the code and therefore it has not been possible for the CPC Program Library to test the program.Memory, required to execute with typical data: 8 Mbytes of RAM memory and 100 MB of Hard disk memoryNo. of bits in a word: 16No. of lines in distributed program, including test data, etc.: 6912No. of bytes in distributed Program, including test data, etc.: 229 541Distribution format: tar.gzNature of the physical problem: the investigation of transport mechanisms for radioactive substances, through environmental pathways, is very important for radiological protection of populations. One such pathway, associated with the food chain, is the grass-animal-man sequence. The distribution of trace elements in humans and laboratory animals has been intensively studied over the past 60 years [R.C. Pendlenton, C.W. Mays, R.D. Lloyd, A.L. Brooks, Differential accumulation of iodine-131 from local fallout in people and milk, Health Phys. 9 (1963) 1253-1262]. In addition, investigations on the incidence of cancer in humans, and a possible causal relationship to radioactive fallout, have been undertaken [E.S. Weiss, M.L. Rallison, W.T. London, W.T. Carlyle Thompson, Thyroid nodularity in southwestern Utah school children exposed to fallout radiation, Amer. J. Public Health 61 (1971) 241-249; M.L. Rallison, B.M. Dobyns, F.R. Keating, J.E. Rall, F.H. Tyler, Thyroid diseases in children, Amer. J. Med. 56 (1974) 457-463; J.L. Lyon, M.R. Klauber, J.W. Gardner, K.S. Udall, Childhood leukemia associated with fallout from nuclear testing, N. Engl. J. Med. 300 (1979) 397-402]. From the pathways of entry of radionuclides in the human (or animal) body, ingestion is the most important because it is closely related to life-long alimentary (or dietary) habits. Those radionuclides which are able to enter the living cells by either metabolic or other processes give rise to localized doses which can be very high. The evaluation of these internally localized doses is of paramount importance for the assessment of radiobiological risks and radiological protection. The time behavior of trace concentration in organs is the principal input for prediction of internal doses after acute or chronic exposure. The General Multiple-Compartment Model (GMCM) is the powerful and more accepted method for biokinetical studies, which allows the calculation of concentration of trace elements in organs as a function of time, when the flow parameters of the model are known. However, few biokinetics data exist in the literature, and the determination of flow and transfer parameters by statistical fitting for each system is an open problem.Restriction on the complexity of the problem: This version of the code works with the constant volume approximation, which is valid for many situations where the biological half-live of a trace is lower than the volume rise time. Another restriction is related to the central flux model. The model considered in the code assumes that exist one central compartment (e.g., blood), that connect the flow with all compartments, and the flow between other compartments is not included.Typical running time: Depends on the choice for calculations. Using the Derivative Method the time is very short (a few minutes) for any number of compartments considered. When the Gauss-Marquardt iterative method is used the calculation time can be approximately 5-6 hours when similar to 15 compartments are considered. (C) 2006 Elsevier B.V. All rights reserved.

Modeling and simulation of cephalosporin C production in a fed-batch tower-type bioreactor

Immobilized cell utilization in tower-type bioreactor is one of the main alternatives being studied to improve the industrial bioprocess. Other alternatives for the production of beta -lactam antibiotics, such as a cephalosporin C fed-batch p recess in an aerated stirred-tank bioreactor with free cells of Cepha-losporium acremonium or a tower-type bioreactor with immobilized cells of this fungus, have proven to be more efficient than the batch profess. In the fed-batch process, it is possible to minimize the catabolite repression exerted by the rapidly utilization of carbon sources (such as glucose) in the synthesis of antibiotics by utilizing a suitable flow rate of supplementary medium. In this study, several runs for cephalosporin C production, each lasting 200 h, were conducted in a fed-batch tower-type bioreactor using different hydrolyzed sucrose concentrations, For this study's model, modifications were introduced to take into account the influence of supplementary medium flow rate. The balance equations considered the effect of oxygen limitation inside the bioparticles. In the Monod-type rate equations, eel concentrations, substrate concentrations, and dissolved oxygen were included as reactants affecting the bioreaction rate. The set of differential equations was solved by the numerical method, and the values of the parameters were estimated by the classic nonlinear regression method following Marquardt's procedure with a 95% confidence interval. The simulation results showed that the proposed model fit well with the experimental data,and based on the experimental data and the mathematical model an optimal mass flow rate to maximize the bioprocess productivity could be proposed.

Simulation of free surface flows in a distributed memory environment

A parallel technique, for a distributed memory machine, based on domain decomposition for solving the Navier-Stokes equations in cartesian and cylindrical coordinates in two dimensions with free surfaces is described. It is based on the code by Tome and McKee (J. Comp. Phys. 110 (1994) 171-186) and Tome (Ph.D. Thesis, University of Strathclyde, Glasgow, 1993) which in turn is based on the SMAC method by Amsden and Harlow (Report LA-4370, Los Alamos Scientific Laboratory, 1971), which solves the Navier-Stokes equations in three steps: the momentum and Poisson equations and particle movement, These equations are discretized by explicit and 5-point finite differences. The parallelization is performed by splitting the computation domain into vertical panels and assigning each of these panels to a processor. All the computation can then be performed using nearest neighbour communication. Test runs comparing the performance of the parallel with the serial code, and a discussion of the load balancing question are presented. PVM is used for communication between processes. (C) 1999 Elsevier B.V. B.V. All rights reserved.

Miniaturized structures for removal or selection of particles from liquid flow

The aim of this work was the development of miniaturized structures useful for retention and/or selection of particles and viscous substances from a liquid flow. The proposed low costs structures are similar to macroscopic wastewater treatment systems, named baffles, and allow disassemble. They were simulated using FEMLAB 3.2b package and manufactured in acrylic with conventional tools. Tests for retention or selection of particles in water or air and viscous fluids in water were carried out. Either in air or water particles with 50 mu m diameter will be retained but not with 13 mu m diameter. In aqueous flow, it is also possible the retention of viscous samples, such as silicone 350 cSt. The simulated results showed good agreement with experimental measurements. These miniaturized structures can be useful in sample pretreatment for chemical analysis and microorganism manipulation. (C) 2007 Elsevier B.V. All rights reserved.

Cluster identification and characterization in the riser of a circulating fluidized bed from numerical simulation results

A methodology of identification and characterization of coherent structures mostly known as clusters is applied to hydrodynamic results of numerical simulation generated for the riser of a circulating fluidized bed. The numerical simulation is performed using the MICEFLOW code, which includes the two-fluids IIT's hydrodynamic model B. The methodology for cluster characterization that is used is based in the determination of four characteristics, related to average life time, average volumetric fraction of solid, existing time fraction and frequency of occurrence. The identification of clusters is performed by applying a criterion related to the time average value of the volumetric solid fraction. A qualitative rather than quantitative analysis is performed mainly owing to the unavailability of operational data used in the considered experiments. Concerning qualitative analysis, the simulation results are in good agreement with literature. Some quantitative comparisons between predictions and experiment were also presented to emphasize the capability of the modeling procedure regarding the analysis of macroscopic scale coherent structures. (c) 2007 Elsevier B.V. All rights reserved.

Flow evolution of a metallic slurry during thixoextrusion of an automotive valve

The aim of this work is to characterize a metallic slurry (Al-4.5%Cu) flow during thixoforming of an automotive valve. The necessary globular structure was obtained by first inoculating the alloy with TIBAL (5%Ti, 1%B, Al - rest) at 750.0°C, and then reheating to a state between liquidus and solidus prior to thixoforming. Two metallic slurries, with a solid phase of approximately 86.1 and 78.2 percent, were used to generate different experimental flow patterns during the thixoforming process. The flow of the material into the die was observed for total, and partial displacement (2.7, 5.4, 7.5mm) of the punch. The first displacement shows formation of the valve rod. The patterns at each step of displacement of the punch were preserved by quenching in water, thus revealing the profile of the die fill and microstructural evolution. Degeneration of the globular phase was observed along the piece thixoextruded. Thixoextrusion forces versus time curves were generated for partial and full displacement of the punch. Porosity was visible along the billet prior to thixoforming. However, some areas show that the porosity gradually decreased to zero as the thixoextrusion pressure increased. Turbulent, transient and laminar flow are analyzed in this work.

Reliability assessment of energy limited systems using sequential montecarlo simulation

This paper explains why the reliability assessment of energy limited systems requires more detailed models for primary generating resources availability, internal and external generating dispatch and customer demand than the ones commonly used for large power systems and presents a methodology based on the full sequential Montecarlo simulation technique with AC power flow for their long term reliability assessment which can properly include these detailed models. By means of a real example, it is shown how the simplified modeling traditionally used for large power systems leads to pessimistic predictions if it is applied to an energy limited system and also that it cannot predict all the load point adequacy problems. © 2006 IEEE.

Development of a micro heat exchanger made with ceramic multi-layers (LTCC) and its setup to gas flow measurements

A green ceramic tape micro heat exchanger was developed using LTCC technology. The device was designed by using a CAD software and 2D and 3D simulations using a CFD package (COMSOL Multiphysics) to evaluate the fluid behavior in the microchannels. The micro heat exchanger is composed of five thermal exchange plates in cross flow arrangement and two connecting plates; heat exchanger dimensions are 26 × 26 × 6 mm3. Preliminary tests were carried out to characterize the device both in atmospheric pressure and in vacuum. The same techniques used in vacuum technology were applied to check the rotameters and to prevent device leakages. Thermal performance of the micro heat exchanger was experimentally tested. © 2009 Elsevier B.V. All rights reserved.