A generalized alternating-direction implicit scheme for incompressible magnetohydrodynamic viscous flows at low magnetic Reynolds number

This paper presents numerical simulations of incompressible fluid flows in the presence of a magnetic field at low magnetic Reynolds number. The equations governing the flow are the Navier-Stokes equations of fluid motion coupled with Maxwell's equations of electromagnetics. The study of fluid flows under the influence of a magnetic field and with no free electric charges or electric fields is known as magnetohydrodynamics. The magnetohydrodynamics approximation is considered for the formulation of the non-dimensional problem and for the characterization of similarity parameters. A finite-difference technique is used to discretize the equations. In particular, an extension of the generalized Peaceman and Rachford alternating-direction implicit (ADI) scheme for simulating two-dimensional fluid flows is presented. The discretized conservation equations are solved in stream function-vorticity formulation. We compare the ADI and generalized ADI schemes, and show that the latter is more efficient in simulating low Reynolds number and magnetic Reynolds number problems. Numerical results demonstrating the applicability of this technique are also presented. The simulation of incompressible magneto hydrodynamic fluid flows is illustrated by numerical solution for two-dimensional cases. (c) 2007 Elsevier B.V. All rights reserved.

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.

Transition to turbulence in the Reynolds' experiment

In this study we simulate numerically the Reynolds' experiment for the transition from laminar to turbulent flow in a pipe. We present a discussion of the results from a dynamical systems perspective when a control parameter, the Reynolds number, is increased. The Landau scenario, where the transition is described by the excitation of infinite oscillatory modes within the fluid, is not observed. Instead what happens is best explained by the Ruelle-Takens scenario in terms of strange attractors. The Lyapunov exponent and fractal dimension for the attractor are calculated together with a measure of complex behaviour called the Lempel-Ziv complexity. (C) 2001 Elsevier B.V. B.V. All rights reserved.

NONLINEAR DIFFUSION PROCESS IN A BENARD SYSTEM AT THE CRITICAL-POINT FOR THE ONSET OF CONVECTION

The evolution equation governing surface perturbations of a shallow fluid heated from below at the critical Rayleigh number for the onset of convective motion, and with boundary conditions leading to zero critical wave number, is obtained. A solution for negative or cooling perturbations is explicitly exhibited, which shows that the system presents sharp propagating fronts.

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.

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

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

An investigation on parallel, divergent and convergent acetylene dual jet diffusion flames

Results are presented and discussed of an experimental investigation on acetylene turbulent dual jet diffusion flames. The study includes parameters of flames in parallel, divergent and convergent configurations. Tests with two parallel jets with addition of helium in the fuel stream were also performed and analysed. The variation of overall flame length and of other name physical characteristics, such as width, volume and conditions for lifting, are presented as functions of burner tip Reynolds number, jet distance from each other and inclination angle. The effects of diluent concentration in the fuel gas stream are presented for single and two parallel jets. (C) 1999 Elsevier B.V. Ltd.

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.

A study of a numerical solution of the steady two dimensions Navier-Stokes equations in a constricted channel problem by a compact fourth order method

We present a numerical solution for the steady 2D Navier-Stokes equations using a fourth order compact-type method. The geometry of the problem is a constricted symmetric channel, where the boundary can be varied, via a parameter, from a smooth constriction to one possessing a very sharp but smooth corner allowing us to analyse the behaviour of the errors when the solution is smooth or near singular. The set of non-linear equations is solved by the Newton method. Results have been obtained for Reynolds number up to 500. Estimates of the errors incurred have shown that the results are accurate and better than those of the corresponding second order method. (C) 2002 Elsevier B.V. All rights reserved.

Non-newtonian flow and pressure drop of pineapple juice in a plate heat exchanger

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

Friction losses in valves and fittings for liquid food products

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

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.

Friction losses in valves and fittings for power-law fluids

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

Demulsification of water/sunflower oil emulsions by a tangential filtration process using chemically impregnated ceramic tubes

A tangential filtration process was implemented in this study using porous ceramic tubes made of alpha-alumina produced by the slip-casting technique. These tubes were sintered at 1450 degrees C and characterized by mercury intrusion porosimetry, which revealed a mean pore size of 0.5 mu m. The tubes were chemically impregnated with a zirconium citrate solution, after which they were calcined and heat treated at temperatures of up to 600 and 900 degrees C to eliminate volatile organic compounds and transform the zirconium citrate into zirconium oxide impregnated in the alumina in the form of nanoparticle agglomerates. The microporous pipes were tested on a microfiltration hydraulic system to analyze their performance in the demulsification of sunflower oil and water mixtures. The fluid-dynamic parameters of Reynolds number and transmembrane pressure were varied in the process. The volume of permeate was analyzed by measuring the Total Organic Carbon concentration (TOC), which indicated 99% of oil phase retention. The emulsified mixture was characterized by optical microscopy, while the morphology and composition of the impregnated microporous tubes were analyzed by scanning electron microscopy (SEM). Quantification of the TOC values for the tube impregnated once at 600 degrees C showed the best demulsification performance, with the concentration on permeate smaller than 10 mg/L. The impregnated tube sintered once at 900 degrees C presented low carbon concentration (smaller than 20 mg/L), has the advantage of presenting the greatest trans-membrane flux in relation to the other microporous tube. (c) 2006 Elsevier B.V. All rights reserved.

Laminar flow of soursop juice through concentric annuli: Friction factors and rheology

Laminar axial flow of a pseudoplastic fluid food (soursop, juice) in annular ducts has been experimentally investigated. In the first part of the manuscript, the rheological behavior of soursop juice, being essential for the annular flow analysis, was completely determined from 9.3 to 49.4 degrees Brix and temperatures from 0.4 degrees C to 68.8 degrees C, using a rotational rheometer equipped with coaxial cylinders. In order to test the adequacy of the rheology results, pressure loss data in the laminar pipe flow were collected and then experimental and theoretical friction factors were compared, showing excellent agreement, which indicated the reliability of the Power-Law model for describing the soursop juices. In the second part, pressure loss in annular regions was measured and used to estimate friction factors, which were then compared to those resulted from analytical and semi-analytical equations. The principal contributions of this article are to provide a review on the determination of friction factors-Reynolds number of pseudoplastic fluids in annuli, and also supply extensive new experimental data on the rheological properties and pressure loss of an important shear-thinning fluid food, which is of particular interest for the food engineering process design. (c) 2006 Elsevier Ltd. All rights reserved.