Evaluation of a bounded high order upwind scheme for 3D incompressible free surface flow computations

fit the context of normalized variable formulation (NVF) of Leonard and total variation diminishing (TVD) constraints of Harten. this paper presents an extension of it previous work by the authors for solving unsteady incompressible flow problems. The main contributions of the paper are threefold. First, it presents the results of the development and implementation of a bounded high order upwind adaptative QUICKEST scheme in the 3D robust code (Freeflow), for the numerical solution of the full incompressible Navier-Stokes equations. Second, it reports numerical simulation results for 1D hock tube problem, 2D impinging jet and 2D/3D broken clam flows. Furthermore, these results are compared with existing analytical and experimental data. And third, it presents the application of the numerical method for solving 3D free surface flow problems. (C) 2007 IMACS. Published by Elsevier B.V. All rights reserved,

Galerkin finite element method for incompressible thermofluid flows framed within the bond graph theory

In this paper a bond graph methodology is used to model incompressible fluid flows with viscous and thermal effects. The distinctive characteristic of these flows is the role of pressure, which does not behave as a state variable but as a function that must act in such a way that the resulting velocity field has divergence zero. Velocity and entropy per unit volume are used as independent variables for a single-phase, single-component flow. Time-dependent nodal values and interpolation functions are introduced to represent the flow field, from which nodal vectors of velocity and entropy are defined as state variables. The system for momentum and continuity equations is coincident with the one obtained by using the Galerkin method for the weak formulation of the problem in finite elements. The integral incompressibility constraint is derived based on the integral conservation of mechanical energy. The weak formulation for thermal energy equation is modeled with true bond graph elements in terms of nodal vectors of temperature and entropy rates, resulting a Petrov-Galerkin method. The resulting bond graph shows the coupling between mechanical and thermal energy domains through the viscous dissipation term. All kind of boundary conditions are handled consistently and can be represented as generalized effort or flow sources. A procedure for causality assignment is derived for the resulting graph, satisfying the Second principle of Thermodynamics. (C) 2007 Elsevier B.V. All rights reserved.

Numerical investigation of the three-dimensional secondary instabilities in the time-developing compressible mixing layer

Mixing layers are present in very different types of physical situations such as atmospheric flows, aerodynamics and combustion. It is, therefore, a well researched subject, but there are aspects that require further studies. Here the instability of two-and three-dimensional perturbations in the compressible mixing layer was investigated by numerical simulations. In the numerical code, the derivatives were discretized using high-order compact finite-difference schemes. A stretching in the normal direction was implemented with both the objective of reducing the sound waves generated by the shear region and improving the resolution near the center. The compact schemes were modified to work with non-uniform grids. Numerical tests started with an analysis of the growth rate in the linear regime to verify the code implementation. Tests were also performed in the non-linear regime and it was possible to reproduce the vortex roll-up and pairing, both in two-and three-dimensional situations. Amplification rate analysis was also performed for the secondary instability of this flow. It was found that, for essentially incompressible flow, maximum growth rates occurred for a spanwise wavelength of approximately 2/3 of the streamwise spacing of the vortices. The result demonstrated the applicability of the theory developed by Pierrehumbet and Widnall. Compressibility effects were then considered and the maximum growth rates obtained for relatively high Mach numbers (typically under 0.8) were also presented.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

A geometric mass-preserving redistancing scheme for the level set function

In this paper we describe and evaluate a geometric mass-preserving redistancing procedure for the level set function on general structured grids. The proposed algorithm is adapted from a recent finite element-based method and preserves the mass by means of a localized mass correction. A salient feature of the scheme is the absence of adjustable parameters. The algorithm is tested in two and three spatial dimensions and compared with the widely used partial differential equation (PDE)-based redistancing method using structured Cartesian grids. Through the use of quantitative error measures of interest in level set methods, we show that the overall performance of the proposed geometric procedure is better than PDE-based reinitialization schemes, since it is more robust with comparable accuracy. We also show that the algorithm is well-suited for the highly stretched curvilinear grids used in CFD simulations. Copyright (C) 2010 John Wiley & Sons, Ltd.

Hybrid modeling of convective laminar flow in a permeable tube associated with the cross-flow process

The confined flows in tubes with permeable surfaces arc associated to tangential filtration processes (microfiltration or ultrafiltration). The complexity of the phenomena do not allow for the development of exact analytical solutions, however, approximate solutions are of great interest for the calculation of the transmembrane outflow and estimate of the concentration, polarization phenomenon. In the present work, the generalized integral transform technique (GITT) was employed in solving the laminar and permanent flow in permeable tubes of Newtonian and incompressible fluid. The mathematical formulation employed the parabolic differential equation of chemical species conservation (convective-diffusive equation). The velocity profiles for the entrance region flow, which are found in the connective terms of the equation, were assessed by solutions obtained from literature. The velocity at the permeable wall was considered uniform, with the concentration at the tube wall regarded as variable with an axial position. A computational methodology using global error control was applied to determine the concentration in the wall and concentration boundary layer thickness. The results obtained for the local transmembrane flux and the concentration boundary layer thickness were compared against others in literature. (C) 2007 Elsevier B.V. All rights reserved.

Gas flow analysis in a Kraft recovery boiler

Demands for optimal boiler performance and increased concerns in lowering emission have always been the driving force in the reevaluation and evolution of the Kraft boiler: specifically the air distribution strategies that are directly related to achieving increased residence time of flue gas combustion inside the furnace which in turn lowers atmosphere emission levels and enhances boiler operation. This paper presents the results of a study that analyzes the interaction of the different multilevel air injections have on flue gas flow patterns including various quaternary air supply arrangements. Additionally, this study assesses the performance of the CFD (Computational Fluid Dynamics) model against data available in literature. Simulations were performed considering isothermal and incompressible flows, and did not take into account thermal phenomena or chemical reactions. The numerical solutions generated proved to be coherently related to the data available in literature, and provided proof of the efficiency of tertiary level air injection, as well as revealed that quaternary air injection ports arranged in a symmetrical configuration is most suitable for optimal equipment operation. (C) 2010 Elsevier B.V. All rights reserved.

Multiphase flow in the vascular system of wood: From microscopic exploration to 3-D Lattice Boltzmann experiments

This paper provides insights into liquid free water dynamics in wood vessels based on Lattice Boltzmann experiments. The anatomy of real wood samples was reconstructed from systematic 3-D analyses of the vessel contours derived from successive microscopic images. This virtual vascular system was then used to supply fluid-solid boundary conditions to a two-phase Lattice Boltzmann scheme and investigate capillary invasion of this hydrophilic porous medium. Behavior of the liquid phase was strongly dependent on anatomical features, especially vessel bifurcations and reconnections. Various parameters were examined in numerical experiments with ideal vessel bifurcations, to clarify our interpretation of these features. (c) 2010 Elsevier Ltd. All rights reserved.

Uso do teste ML-Flow como auxiliar na classificação e tratamento da hanseníase

FUNDAMENTOS: O tratamento da hanseníase é definido pela classificação de pacientes em paucibacilares (PB) e multibacilares (MB). A OMS (Organização Mundial de Saúde) classifica os doentes de acordo com o número de lesões, mas Ridley-Jopling (R&J) utiliza também exames complementares, porém é de difícil utilização fora dos serviços de referência. Em 2003 foi desenvolvido um teste denominado ML-Flow, uma alternativa à sorologia por ELISA para auxiliar na classificação de pacientes em PB e MB e auxiliar na decisão terapêutica. OBJETIVOS: Observar a concordância entre o teste de ML-Flow e baciloscopia de linfa, exame já consagrado para detecção de MB. Analisar a utilidade do teste de ML-Flow em campo. MATERIAL E MÉTODOS: Estudo retrospectivo avaliando prontuário de 55 pacientes virgens de tratamento, diagnosticados como PB ou MB por R&J. Submetidos à baciloscopia e ao teste de ML-Flow. RESULTADOS: Nos MB, a baciloscopia foi positiva em 80% dos casos, o ML-flow foi positivo em 82,5%. Entre os PB, o ML-Flow foi positivo em 37,5% e a baciloscopia do esfregaço foi negativa em 100% dos casos. A concordância entre os resultados da baciloscopia do esfregaço e ML-Flow foi de 87,5%, kappa=0,59, p<0,001. CONCLUSÃO: Nenhum teste laboratorial é 100% sensível e específico para a correta classificação de todas as formas de hanseníase. O ML-Flow é um teste rápido, de fácil manuseio em campo, menos invasivo que a baciloscopia podendo ser útil para auxiliar na decisão terapêutica em locais de difícil acesso a serviços de referência.

Forced convection to laminar flow of liquid egg yolk in circular and annular ducts

The steady-state heat transfer in laminar flow of liquid egg yolk - an important pseudoplastic fluid food - in circular and concentric annular ducts was experimentally investigated. The average convection heat transfer coefficients, determined by measuring temperatures before and after heating sections with constant temperatures at the tube wall, were used to obtain simple new empirical expressions to estimate the Nusselt numbers for fully established flows at the thermal entrance of the considered geometries. The comparisons with existing correlations for Newtonian and non-Newtonian fluids resulted in excellent agreement. The main contribution of this work is to supply practical and easily applicable correlations, which are, especially for the case of annulus, rather scarce and extensively required in the design of heat transfer operations dealing with similar shear-thinning products. In addition, the experimental results may support existing theoretical analyses.

An Eulerian Immersed Boundary Method for flow simulations over stationary and moving rigid bodies

The fluid flow over bodies with complex geometry has been the subject of research of many scientists and widely explored experimentally and numerically. The present study proposes an Eulerian Immersed Boundary Method for flows simulations over stationary or moving rigid bodies. The proposed method allows the use of Cartesians Meshes. Here, two-dimensional simulations of fluid flow over stationary and oscillating circular cylinders were used for verification and validation. Four different cases were explored: the flow over a stationary cylinder, the flow over a cylinder oscillating in the flow direction, the flow over a cylinder oscillating in the normal flow direction, and a cylinder with angular oscillation. The time integration was carried out by a classical 4th order Runge-Kutta scheme, with a time step of the same order of distance between two consecutive points in x direction. High-order compact finite difference schemes were used to calculate spatial derivatives. The drag and lift coefficients, the lock-in phenomenon and vorticity contour plots were used for the verification and validation of the proposed method. The extension of the current method allowing the study of a body with different geometry and three-dimensional simulations is straightforward. The results obtained show a good agreement with both numerical and experimental results, encouraging the use of the proposed method.

Electrochemical degradation of the dye reactive orange 16 using electrochemical flow-cell

Electrochemical removals of color and organic load from solutions containing the dye reactive orange 16 (RO16) were performed in an electrochemical flow-cell, using a platinum working electrode. The influence of the process variables flow-rate, such as NaCl concentration, applied potential and solution pH, were studied. The best color removal achieved was 93% (λ = 493 nm) after 60 min at 2.2 V vs. RHE electrolysis, using 1.00 g L-1 NaCl as supporting electrolyte. The rises in the concentration of NaCl and applied potential increased the color removal rate. The best total organic carbon removal (57%) was obtained at 1.8 V, without the separating membrane, indicating that the ideal conditions for the color removal are not necessarily the same as those to remove the total organic carbon. The degradation efficiency decreased with the solution pH decrease.

Evaluation of a simple hyphenated system for flow injection solid-phase pre-concentration and capillary electrophoresis

In this work, the development and evaluation of a hyphenated flow injection-capillary electrophoresis system with on-line pre-concentration is described. Preliminary tests were performed to investigate the influence of flow rates over the analytical signals. Results revealed losses in terms of sensitivity of the FIA-CE system when compared to the conventional CE system. To overcome signal decrease and to make the system more efficient, a lower flow rate was set and an anionic resin column was added to the flow manifold in order to pre-concentrate the analyte. The pre-concentration FIA-CE system presented a sensitivity improvement of about 660% and there was only a small increase of 8% in total peak dispersion. These results have confirmed the great potential of the proposed system for many analytical tasks especially for low concentration samples.

Determination of paracetamol at a graphite-polyurethane composite electrode as an amperometric flow detector

A bare graphite-polyurethane composite was evaluated as an amperometric flow injection detector in the determination of paracetamol (APAP) in pharmaceutical formulations. A linear analytical curve was observed in the 5.00 x 10-5 to 5.00 x 10-3 mol L-1 range with a minimum detectable net concentration of 18.9 µmol L-1 and 180 determinations h-1, after optimization of parameters such as the detection potential, sample loop volume, and carrier solution flow rate. Interference of ascorbic acid was observed, however, it was possible overcome the interference, reaching results that agreed with HPLC within 95% confidence level. These results showed that the graphite-polyurethane composite can be used as an amperometric detector for flow analysis in the determination of APAP.

Development of a flow through photo-reactor to study degradation of organic compounds by Sequential Injection Analysis (SIA)

This work describes a photo-reactor to perform in line degradation of organic compounds by photo-Fenton reaction using Sequential Injection Analysis (SIA) system. A copper phthalocyanine-3,4',4²,4²¢-tetrasulfonic acid tetrasodium salt dye solution was used as a model compound for the phthalocyanine family, whose pigments have a large use in automotive coatings industry. Based on preliminary tests, 97% of color removal was obtained from a solution containing 20 µmol L-1 of this dye.