Solution of a heterogeneous modeling of a horizontal-flow anaerobic immobilized biomass (HAIB) reactor by the sequencing method

A modeling study was completed to develop a methodology that combines the sequencing and finite difference methods for the simulation of a heterogeneous model of a tubular reactor applied in the treatment of wastewater. The system included a liquid phase (convection diffusion transport) and a solid phase (diffusion reaction) that was obtained by completing a mass balance in the reactor and in the particle, respectively. The model was solved using a pilot-scale horizontal-flow anaerobic immobilized biomass (HAIB) reactor to treat domestic sewage, with the concentration results compared with the experimental data. A comparison of the behavior of the liquid phase concentration profile and the experimental results indicated that both the numerical methods offer a good description of the behavior of the concentration along the reactor. The advantage of the sequencing method over the finite difference method is that it is easier to apply and requires less computational time to model the dynamic simulation of outlet response of HAIB.

The influence of immobilized biomass and particle size on the fluid dynamics of an anaerobic fluidized bed reactor

Fluid dynamic analysis is an important branch of several chemical engineering related areas, such as drying processes and chemical reactors. However, aspects concerning fluid dynamics in wastewater treatment bioreactors still require further investigation, as they highly influence process efficiency. Therefore, it is essential to evaluate the influence of biofilm on the reactor fluid dynamic behavior, through the analysis of a few important parameters, such as minimum fluidization velocity, bed expansion and porosity, and particle terminal velocity. The main objective of the present work was to investigate the fluid dynamics of an anaerobic fluidized bed reactor, having activated carbon particles as support media for biomass immobilization. Reactor performance was tested using synthetic residual water, which was prepared using the solution employed in BOD determination. The results showed that the presence of immobilized biomass increased particle density and altered the main fluid dynamic parameters investigated.

Macro-Fiber Composite Actuators for Flow Control of a Variable Camber Airfoil

This research employs solid-state actuators for delay of flow separation seen in airfoils at low Reynolds numbers. The flow control technique investigated here is aimed for a variable camber airfoil that employs two active surfaces and a single four-bar (box) mechanism as the internal structure. To reduce separation, periodic excitation to the flow around the leading edge of the airfoil is induced by a total of nine piezocomposite actuated clamped-free unimorph benders distributed in the spanwise direction. An electromechanical model is employed to design an actuator capable of high deformations at the desired frequency for lift improvement at post-stall angles. The optimum spanwise distribution of excitation for increasing lift coefficient is identified experimentally in the wind tunnel. A 3D (non-uniform) excitation distribution achieved higher lift enhancement in the post-stall region with lower power consumption when compared to the 2D (uniform) excitation distribution. A lift coefficient increase of 18.4% is achieved with the identified non-uniform excitation mode at the bender resonance frequency of 125 Hz, the flow velocity of 5 m/s and at the reduced frequency of 3.78. The maximum lift (Clmax) is increased 5.2% from the baseline. The total power consumption of the flow control technique is 639 mW(RMS).

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.

Piezoceramic Composite Actuators for Flow Control in Low Reynolds Number Airflow

The research presented here employs solid-state actuators for flow separation delay or for forced attachment of separated flow seen in airfoils at low Reynolds numbers. To reduce separation, periodic excitation to the flow around the leading edge of the airfoil is induced by Macro-Fiber Composite actuated clamped-free unimorph benders. An electromechanical model of the unimorph is briefly presented and parametric study is conducted to aid the design of a unimorph to output high deformation at a desired frequency. The optimum frequency and amplitude for lift improvement at post-stall angles are identified experimentally. Along with aerodynamic force and structural displacement measurements, helium bubble flow visualization is used to verify existing separated flow, and the attached flow induced by flow control. The lift enhancement induced by several flow control techniques is compared. A symmetric and non-uniform (3D) flow excitation results in the maximum lift enhancement at post-stall region at the lowest power consumption level. A maximum lift coefficient increase of 27.5% (in the post-stall region) is achieved at 125 Hz periodic excitation, with the 3D symmetric actuation mode at 5 m/s and the reduced frequency of 3.78. C(l,max) is increased 7.6% from the baseline.

Optimization of the flow in the catalytic converter of internal combustion engines by means of screens

The flow in the automotive catalytic converter is, in general, not uniform. This significantly affects cost, service life, and performance, in particular, during cold startup. The current paper reports on a device that provided a large improvement in flow uniformity. The device is to be placed in the converter inlet diffuser and is constructed out of ordinary screens. It is cheap and easy to install. Moreover, the device does not present most of the undesired effects, such as increase in pressure drop and time to light off, often observed in other devices developed for the same purpose.

Modeling and experimental study of nucleate boiling on a vertical array of horizontal plain tubes

An investigation of nucleate boiling on a vertical array of horizontal plain tubes is presented in this paper. Experiments were performed with refrigerant RI 23 at reduced pressures varying from 0.022 to 0.64, tube pitch to diameter ratios of 1.32, 1.53 and 2.00, and heat fluxes from 0.5 to 40 kW/m(2). Brass tubes with external diameters of 19.05 mm and average roughness of 0.12 mu m were used in the experiments. The effect of the tube spacing on the local heat transfer coefficient along the tube array was negligible within the present range of experimental conditions. For partial nucleate boiling, characterized by low heat fluxes, and low reduced pressures, the tube positioning shows a remarkable effect on the heat transfer coefficient. Based on these data, a general correlation for the prediction of the nucleate boiling heat transfer coefficient on a vertical array of horizontal tubes under flooded conditions was proposed. According to this correlation, the ratio between the heat transfer coefficients of a given tube and the lowest tube in the array depends only on the tube row number, the reduced pressure and the heat flux. By using the proposed correlation, most of the experimental heat transfer coefficients obtained in the present study were predicted within +/- 15%. The new correlation compares reasonably well with independent data from the literature. (C) 2008 Elsevier Inc. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part I - A two-phase flow pattern map and a flow pattern based phenomenological model for two-phase flow frictional pressure drops

An updated flow pattern map was developed for CO2 on the basis of the previous Cheng-Ribatski-Wojtan-Thome CO2 flow pattern map [1,2] to extend the flow pattern map to a wider range of conditions. A new annular flow to dryout transition (A-D) and a new dryout to mist flow transition (D-M) were proposed here. In addition, a bubbly flow region which generally occurs at high mass velocities and low vapor qualities was added to the updated flow pattern map. The updated flow pattern map is applicable to a much wider range of conditions: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to +25 degrees C (reduced pressures from 0.21 to 0.87). The updated flow pattern map was compared to independent experimental data of flow patterns for CO2 in the literature and it predicts the flow patterns well. Then, a database of CO2 two-phase flow pressure drop results from the literature was set up and the database was compared to the leading empirical pressure drop models: the correlations by Chisholm [3], Friedel [4], Gronnerud [5] and Muller-Steinhagen and Heck [6], a modified Chisholm correlation by Yoon et al. [7] and the flow pattern based model of Moreno Quiben and Thome [8-10]. None of these models was able to predict the CO2 pressure drop data well. Therefore, a new flow pattern based phenomenological model of two-phase flow frictional pressure drop for CO2 was developed by modifying the model of Moreno Quiben and Thome using the updated flow pattern map in this study and it predicts the CO2 pressure drop database quite well overall. (C) 2007 Elsevier Ltd. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part II - An updated general flow boiling heat transfer model based on flow patterns

Corresponding to the updated flow pattern map presented in Part I of this study, an updated general flow pattern based flow boiling heat transfer model was developed for CO2 using the Cheng-Ribatski-Wojtan-Thome [L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes, Int. J. Heat Mass Transfer 49 (2006) 4082-4094; L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, Erratum to: ""New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside tubes"" [Heat Mass Transfer 49 (21-22) (2006) 4082-4094], Int. J. Heat Mass Transfer 50 (2007) 391] flow boiling heat transfer model as the starting basis. The flow boiling heat transfer correlation in the dryout region was updated. In addition, a new mist flow heat transfer correlation for CO2 was developed based on the CO2 data and a heat transfer method for bubbly flow was proposed for completeness sake. The updated general flow boiling heat transfer model for CO2 covers all flow regimes and is applicable to a wider range of conditions for horizontal tubes: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to 25 degrees C (reduced pressures from 0.21 to 0.87). The updated general flow boiling heat transfer model was compared to a new experimental database which contains 1124 data points (790 more than that in the previous model [Cheng et al., 2006, 2007]) in this study. Good agreement between the predicted and experimental data was found in general with 71.4% of the entire database and 83.2% of the database without the dryout and mist flow data predicted within +/-30%. However, the predictions for the dryout and mist flow regions were less satisfactory due to the limited number of data points, the higher inaccuracy in such data, scatter in some data sets ranging up to 40%, significant discrepancies from one experimental study to another and the difficulties associated with predicting the inception and completion of dryout around the perimeter of the horizontal tubes. (C) 2007 Elsevier Ltd. All rights reserved.

Debris flow affecting the Cubatauo Oil Refinery, Brazil

On February 6, 1994, a large debris flow developed because of intense rains in a 800-m-high mountain range called Serra do Cubatao, the local name for the Serra do Mar, located along the coast of the state of Sao Paulo, Brazil. It affected the Presidente Bernardes Refinery, owned by Petrobras, in Cubatao. The damages amounted to about US $40 million because of the muck cleaning, repairs, and 3-week interruption of the operations. This prompted Petrobras to conduct studies, carried out by the authors, to develop protection works, which were done at a cost of approximately US $12 million. The paper describes the studies conducted on debris flow mechanics. A new criteria to define rainfall intensities that trigger debris flows is presented, as well as a correlation of slipped area with soil porosity and rain intensity. Also presented are (a) an actual grain size distribution of a deposited material, determined by laboratory and a large-scale field test, and (b) the size distribution of large boulders along the river bed. Based on theory, empirical experience and back-analysis of the events, the main parameters as the front velocity, the peak discharge and the volume of the transported sediments were determined in a rational basis for the design of the protection works. Finally, the paper describes the set of the protection works built, emphasizing their concept and function. They also included some low-cost innovative works.

A comparison of dynamic and quasi-static results for the flow towards a well

The present analysis takes into account the acceleration term in the differential equation of motion to obtain exact dynamic solutions concerning the groundwater flow towards a well in a confined aquifer. The results show that the error contained in the traditional quasi-static solution is very small in typical situations.

Volumetric reconstruction of the mean flow around circular cylinders fitted with strakes

The volumetric reconstruction technique presented in this paper employs a two-camera stereoscopic particle image velocimetry (SPIV) system in order to reconstruct the mean flow behind a fixed cylinder fitted with helical strakes, which are commonly used to suppress vortex-induced vibrations (VIV). The technique is based on the measurement of velocity fields at equivalent adjacent planes that results in pseudo volumetric fields. The main advantage over proper volumetric techniques is the avoidance of additional equipment and complexity. The averaged velocity fields behind the straked cylinders and the geometrical periodicity of the three-start configuration are used to further simplify the reconstruction process. Two straked cylindrical models with the same pitch (p = 10d) and two different heights (h = 0.1 and 0.2d) are tested. The reconstructed flow shows that the strakes introduce in the wake flow a well-defined wavelength of one-third of the pitch. Measurements of hydrodynamic forces, fluctuating velocity, vortex formation length, and vortex shedding frequency show the interdependence of the wake parameters. The vortex formation length is increased by the strakes, which is an important effect for the suppression of vortex-induced vibrations. The results presented complement previous investigations concerning the effectiveness of strakes as VIV suppressors and provide a basis of comparison to numerical simulations.

On the wake-induced vibration of tandem circular cylinders: the vortex interaction excitation mechanism

The mechanism of wake-induced vibrations (WIV) of a pair of cylinders in a tandem arrangement is investigated by experiments. A typical WIV response is characterized by a build-up of amplitude persisting to high reduced velocities; this is different from a typical vortex-induced vibration (VIV) response, which occurs in a limited resonance range. We suggest that WIV of the downstream cylinder is excited by the unsteady vortex-structure interactions between the body and the upstream wake. Coherent vortices interfering with the downstream cylinder induce fluctuations in the fluid force that are not synchronized with the motion. A favourable phase lag between the displacement and the fluid force guarantees that a positive energy transfer from the flow to the structure sustains the oscillations. If the unsteady vortices are removed from the wake of the upstream body then WIV will not be excited. An experiment performed in a steady shear flow turned out to be central to the understanding of the origin of the fluid forces acting on the downstream cylinder.

3D CFD Simulation of Vortex-induced Vibration of Cylinder

The 3D flow around a circular cylinder free to oscillate transversely to the free stream was simulated using Computational Fluid Dynamics (CFD) and the Spalart-Allmaras Detached Eddy Simulation (DES) turbulence model for a Reynolds number Re = 10(4). Simulations were carried out for a small mass-damping parameter m*zeta = 0.00858, where m* = 3.3 and zeta = 0.0026. We found good agreement between the numerical results and experimental data. The simulations predicted the high observed amplitudes of the upper branch of vortex-induced vibrations for low mass-damping parameters.

Numerical solution of the two-phase expansion of a metastable flashing liquid jet using the dispersion-controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two-dimensional axisymmetric domain, while in the two-phase region. In general, the phenomenon starts off when a high-pressure, high-temperature liquid jet emerges from a small nozzle or orifice expanding into a low-pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two-phase high-speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation-of-state for the substance. A false-transient model is implemented using the shock capturing scheme: dispersion-controlled dissipative (DCD), which was used to calculate the flow conditions as the steady-state condition is reached. Numerical results with computational code DCD-2D vI have been analyzed. Copyright (C) 2009 John Wiley & Sons, Ltd.