A continuously differentiable upwinding scheme for the simulation of fluid flow problems

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

Simplified architecture of a neurofuzzy controller applied on industrial systems for fluid flow control

This work presents a simplified architecture of a neurofuzzy controller for general purpose applications that tries to minimize the processing used in the several stages of hazy modeling of systems. The basic procedures of fuzzification and defuzzification are simplified to the maximum while the inference procedures are computed in a private way. The simplified architecture allows a fast and easy configuration of the neurofuzzy controller and the structuring rules that define the control actions is automatic. Th controller's Limits and performance are standardized and the control actions are previously calculated. For application, the industrial systems of fluid flow control will be considered.

Entropy variation in isothermal fluid flow considering real gas effects

The present paper concerns on the estimative of the pressure loss and entropy variation in an isothermal fluid flow, considering real gas effects. The 1D formulation is based on the isothermal compressibility module and on the thermal expansion coefficient in order to be applicable for both gas and liquid as pure substances. It is emphasized on the simple methodology description, which establishes a relationship between the formulation adopted for ideal gas and another considering real gas effects. A computational procedure has been developed, which can be used to determine the flow properties in duct with a variable area, where real gas behavior is significant. In order to obtain quantitative results, three virial coefficients for Helium equation of state are employed to determine the percentage difference in pressure and entropy obtained from different formulations. Results are presented graphically in the form of real gas correction factors, which can be applied to perfect gas calculations.

Development of a low-Reynolds-number k-ω model for FENE-P fluids

A low-Reynolds-number k-ω model for Newtonian fluids has been developed to predict drag reduction of viscoelastic fluids described by the FENE-P model. The model is an extension to viscoelastic fluids of the model for Newtonian fluids developed by Bredberg et al. (Int J Heat Fluid Flow 23:731-743, 2002). The performance of the model was assessed using results from direct numerical simulations for fully developed turbulent channel flow of FENE-P fluids. It should only be used for drag reductions of up to 50 % (low and intermediate drag reductions), because of the limiting assumption of turbulence isotropy leading to an under-prediction of k, but compares favourably with results from k-ε models in the literature based on turbulence isotropy. © 2012 Springer Science+Business Media Dordrecht.

Transient non-Darcy forced convective heat transfer from a flat plate embedded in a fluid-saturated porous medium

Transient non-Darcy forced convection on a flat plate embedded in a porous medium is investigated using the Forchheimer-extended Darcy law. A sudden uniform pressure gradient is applied along the flat plate, and at the same time, its wall temperature is suddenly raised to a high temperature. Both the momentum and energy equations are solved by retaining the unsteady terms. An exact velocity solution is obtained and substituted into the energy equation, which then is solved by means of a quasi-similarity transformation. The temperature field can be divided into the one-dimensional transient (downstream) region and the quasi-steady-state (upstream) region. Thus the transient local heat transfer coefficient can be described by connecting the quasi-steady-state solution and the one-dimensional transient solution. The non-Darcy porous inertia works to decrease the velocity level and the time required for reaching the steady-state velocity level. The porous-medium inertia delays covering of the plate by the steady-state thermal boundary layer. © 1990.

Experimental analysis of heat transfer in packed beds with air flow

Heat-transfer studies were carried out in a packed bed of glass beads, cooled by the wall, through which air percolated. Tube-to-particle diameter ratios (D/dp) ranged from 1.8 to 55, while the air mass flux ranged from 0.204 to 2.422 kg/m2·s. The outlet bed temperature (TL) was measured by a brass ring-shaped sensor and by aligned thermocouples. The resulting radial temperature profiles differed statistically. Angular temperature fluctuations were observed through measurements made at 72 angular positions. These fluctuations do not follow a normal distribution around the mean for low ratios D/dp. The presence of a restraining screen, as well as the increasing distance between the temperature measuring device and the bed surface, distorts TL. The radial temperature profile at the bed entrance (T0) was measured by a ring-shaped sensor, and T 0 showed to be a function of the radial position, the particle diameter, and the fluid flow rate.

Experimental analysis of heat recoveries with heat pipes operating in low inclination angle

Heat recovery devices are important in the optimization of thermal systems, since they can be used to reduce thermal losses to the environment. The use of heat pipes in these types of equipment can provide heat recoveries of higher efficiency, since both fluid flows are external and there are less contamination risks between the hot and cold fluids. The objective of this work is to study a heat recovery unit constructed with heat pipes and mainly, to analyze the influence of the inclination of the heat pipes on the performance of the equipment. For this analysis, a heat recovery unit was constructed which possesses 48 finned heat pipes in triangular geometry, the evaporator and condenser being of the same length. This unit was tested in an air-air system simulating a heat recovery process in which heat was supplied to the hot fluid by electrical resistances. The results have shown that there exists an inclination at which the system has a better performance, but for higher inclinations there is no significant increase of the efficiency of the system. This paper also presents the influence of inclination of heat pipes on effectiveness and NTU parameters which are important in heat exchanger design.

Influência da lubrirrefrigeração na qualidade superficial de metais retificados

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

Criteria for mixed grids in computational fluid dynamics

The finite volume method is used as a numerical method for solving the fluid flow equations. This method is appropriate to employ under structured and unstructured meshes. Mixed grids, combining both types of grids, are investigated. The coupling of different grids is done by overlapping strategy. The computational effort for the mixed grid is evaluated by the CPU-time, with different percentage of covering area of the unstructured mesh. The present scheme is tested for the driven cavity problem, where the incompressible fluid is integrated by calculating the velocity fields and computing the pressure field in each time step. Several schemes for unstructured grid are examined, and the compatibility condition is applied to check their consistency. A scheme to verify the compatibility condition for the unstructured grids is presented. (c) 2006 IMACS. Published by Elsevier B.V. All rights reserved.

Simulação numérica de trocadores de calor multitubulares equipados com fitas helicoidais externamente aos tubos

Classical shell-and-tube heat exchangers are usually equipped with segmental baffles. These baffles serve two basic functions: (a) they provide tube supports, thereby preventing or reducing mechanical problems, such as sagging or vibration; (b) they direct the fluid flow over the tubes so as to introduce a cross-flow component, thereby increasing the heat transfer. Segmented baffles have several sources of performance loss, some due to various leakage flows and others caused by stagnation zones. A new concept of longitudinal flow heat exchanger - based on placing twisted tapes along the tube bundle subchannels - was developed to mitigate drawbacks of other types of tubular heat exchangers. In this paper, a numerical model has been implemented in order to simulate the thermal-hydraulic feature of tubular heat exchangers equipped either with segmental baffles or with subchannel twisted tapes. The tube bundle has been described by means of an equivalent porous medium type model, allowing a macroscopic description of the shell-side flow. The basic equations - continuity, momentum and energy - have been solved by using the finite volume method. Typical numerical results have been compared with experimental data, reaching a very good agreement. A comparative analysis of different types of heat exchangers has been carried out, revealing the satisfactory thermal-hydraulic efficiency level of the twisted tapes heat exchangers.

Study on transient ice formation of laminar flow inside externally supercooled rectangular duct

The transient process of solidification of laminar liquid flow (water) submitted to super-cooling was investigated both theoretically and experimentally. In this study an alternative analytical formulation and numerical approach were adopted resulting in the unsteady model with temperature dependent thermophysical properties in the solid region. The proposed model is based upon the fundamental equations of energy balance in the solid and liquid regions as well as across the solidification front. The basic equations and the associated boundary and initial conditions were made dimensionless by using the Landau transformation to immobilize the moving front and render the problem to a fixed plane type problem. A laminar velocity profile is admitted in the liquid domain and the resulting equations were discretized using the finite difference approach. The numerical predictions obtained were compared with the available results based on other models and concepts such as Neumann analytical model, the apparent thermal capacity model due to Bonacina and the conventional fixed grid energy model due to Goodrich. To obtain further comparisons and more validation of the model and the numerical solution, an experimental rig was constructed and instrumented permitting very well controlled experimental measurements. The numerical predictions were compared with the experimental results and the agreement was found satisfactory.

Application of a bounded upwinding scheme to complex fluid dynamics problems

This paper is concerned with an overview of upwinding schemes, and further nonlinear applications of a recently introduced high resolution upwind differencing scheme, namely the ADBQUICKEST [V.G. Ferreira, F.A. Kurokawa, R.A.B. Queiroz, M.K. Kaibara, C.M. Oishi, J.A.Cuminato, A.F. Castelo, M.F. Tomé, S. McKee, assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems, International Journal for Numerical Methods in Fluids 60 (2009) 1-26]. The ADBQUICKEST scheme is a new TVD version of the QUICKEST [B.P. Leonard, A stable and accurate convective modeling procedure based on quadratic upstream interpolation, Computer Methods in Applied Mechanics and Engineering 19 (1979) 59-98] for solving nonlinear balance laws. The scheme is based on the concept of NV and TVD formalisms and satisfies a convective boundedness criterion. The accuracy of the scheme is compared with other popularly used convective upwinding schemes (see, for example, Roe (1985) [19], Van Leer (1974) [18] and Arora & Roe (1997) [17]) for solving nonlinear conservation laws (for example, Buckley-Leverett, shallow water and Euler equations). The ADBQUICKEST scheme is then used to solve six types of fluid flow problems of increasing complexity: namely, 2D aerosol filtration by fibrous filters; axisymmetric flow in a tubular membrane; 2D two-phase flow in a fluidized bed; 2D compressible Orszag-Tang MHD vortex; axisymmetric jet onto a flat surface at low Reynolds number and full 3D incompressible flows involving moving free surfaces. The numerical simulations indicate that this convective upwinding scheme is a good generic alternative for solving complex fluid dynamics problems. © 2012.

Análise numérica de condensadores do tipo arame-sobre-tubo usados em refrigeradores domésticos

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

Análise numérica de evaporadores do tipo Roll-Bond usados em refrigeradores domésticos

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

A Numerical Study of the FENE-CR Model Applied to a Jet Flow Problem

The FENE-CR model is investigated through a numerical algorithm to simulate the time-dependent moving free surface flow produced by a jet impinging on a flat surface. The objective is to demonstrate that by increasing the extensibility parameter L, the numerical solutions converge to the solutions obtained with the Oldroyd-B model. The governing equations are solved by an established free surface flow solver based on the finite difference and marker-and-cell methods. Numerical predictions of the extensional viscosity obtained with several values of the parameter L are presented. The results show that if the extensibility parameter L is sufficiently large then the extensional viscosities obtained with the FENE-CR model approximate the corresponding Oldroyd-B viscosity. Moreover, the flow from a jet impinging on a flat surface is simulated with various values of the extensibility parameter L and the fluid flow visualizations display convergence to the Oldroyd-B jet flow results.