Use of computational fluid dynamics in the design of bubble column reactors

Investigations into the modelling techniques that depict the transport of discrete phases (gas bubbles or solid particles) and model biochemical reactions in a bubble column reactor are discussed here. The mixture model was used to calculate gas-liquid, solid-liquid and gasliquid-solid interactions. Multiphase flow is a difficult phenomenon to capture, particularly in bubble columns where the major driving force is caused by the injection of gas bubbles. The gas bubbles cause a large density difference to occur that results in transient multi-dimensional fluid motion. Standard design procedures do not account for the transient motion, due to the simplifying assumptions of steady plug flow. Computational fluid dynamics (CFD) can assist in expanding the understanding of complex flows in bubble columns by characterising the flow phenomena for many geometrical configurations. Therefore, CFD has a role in the education of chemical and biochemical engineers, providing the examples of flow phenomena that many engineers may not experience, even through experimentation. The performance of the mixture model was investigated for three domains (plane, rectangular and cylindrical) and three flow models (laminar, k-e turbulence and the Reynolds stresses). mThis investigation raised many questions about how gas-liquid interactions are captured numerically. To answer some of these questions the analogy between thermal convection in a cavity and gas-liquid flow in bubble columns was invoked. This involved modelling the buoyant motion of air in a narrow cavity for a number of turbulence schemes. The difference in density was caused by a temperature gradient that acted across the width of the cavity. Multiple vortices were obtained when the Reynolds stresses were utilised with the addition of a basic flow profile after each time step. To implement the three-phase models an alternative mixture model was developed and compared against a commercially available mixture model for three turbulence schemes. The scheme where just the Reynolds stresses model was employed, predicted the transient motion of the fluids quite well for both mixture models. Solid-liquid and then alternative formulations of gas-liquid-solid model were compared against one another. The alternative form of the mixture model was found to perform particularly well for both gas and solid phase transport when calculating two and three-phase flow. The improvement in the solutions obtained was a result of the inclusion of the Reynolds stresses model and differences in the mixture models employed. The differences between the alternative mixture models were found in the volume fraction equation (flux and deviatoric stress tensor terms) and the viscosity formulation for the mixture phase.

Otimização dos parâmetros de um modelo de bomba em sistemas de bombeio centrífugo submerso

One of several techniques applied to production processes oil is the artificial lift, using equipment in order to reduce the bottom hole pressure, providing a pressure differential, resulting in a flow increase. The choice of the artificial lift method depends on a detailed analysis of the some factors, such as initial costs of installation, maintenance, and the existing conditions in the producing field. The Electrical Submersible Pumping method (ESP) appears to be quite efficient when the objective is to produce high liquid flow rates in both onshore and offshore environments, in adverse conditions of temperature and in the presence of viscous fluids. By definition, ESP is a method of artificial lift in which a subsurface electric motor transforms electrical into mechanical energy to trigger a centrifugal pump of multiple stages, composed of a rotating impeller (rotor) and a stationary diffuser (stator). The pump converts the mechanical energy of the engine into kinetic energy in the form of velocity, which pushes the fluid to the surface. The objective of this work is to implement the optimization method of the flexible polyhedron, known as Modified Simplex Method (MSM) applied to the study of the influence of the modification of the input and output parameters of the centrifugal pump impeller in the channel of a system ESP. In the use of the optimization method by changing the angular parameters of the pump, the resultant data applied to the simulations allowed to obtain optimized values of the Head (lift height), lossless efficiency and the power with differentiated results.

Dysprosium transport in Nd-Fe-B pellets

The addition of heavy rare earth (RE) elements to Nd2Fe14B based magnets to form (Nd,Dy)2Fe14B is known to increase the coercivity and high temperature performance required for hybrid vehicle electric motors and other extreme temperature applications. Attempts to conserve heavy rare earth elements for high temperature (RE)2Fe14B based magnets have led to the development of a grain boundary diffusion process for bulk magnets. This process relies on transport of a heavy rare earth, such as Dy, into a bulk Nd2Fe14B magnet along pores, a low volume fraction of eutectic liquid along grain boundary grain triple junctions and grain boundaries. This enriches the grain surfaces in Dy through the thickness of the bulk magnet, leading to larger increases coercivity with a smaller Dy concentration than can be achieved with homogeneous alloys. Attempts to carry out the same process during sintering require significant control of Dy transport efficiency. The macroscopic transport of Dy in Nd2.7Fe14B1.4 based powder packs is studied using a 'layered' pellet, where Nd2.7Fe14B1.4powder is an interlayer and Dy source as a center layer. The sintering of this layered pellet provided evidence for very large effective diffusion lengths aided by Dy rich liquid flow through connected porosity. Approaches to controlling Dy transportation include decreasing the liquid phase transport capability of the powder pack by increasing the melting point of the Dy source and the decreasing amount of RE rich liquid in the powder packs. The solid-liquid reaction is studied in which melt spun Nd2.7Fe14B1.4 ribbons are PVD coated with Dy-Fe eutectic composition and then thermally treated. The resulting microstructure from the reaction between Dy-Fe eutectic coating and Nd2.7Fe14B1.4 ribbon is interpreted as support for a proposed dissolution/reprecipitation process between solid and liquid phases. The estimate the diffusion coefficient and the effective diffusion length of Dy sources in Nd2.7Fe14B1.4 layered pellets and melt spun ribbons were obtained from the calculation of Fick's second law combined with EDS results from the experiment. The results indicate that the effective diffusion coefficient of Dy in the layered pellets is higher than the diffusion in ribbons due to its higher porosity than ribbons.

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.

Liquid-liquid microextraction without phase separation in a multicommuted flow system for diltiazem determination in pharmaceuticals

Liquid-liquid microextraction without phase segmentation was implemented in a multicommuted flow system for determination of the anti-hypertensive diltiazem. The procedure was based on ion pair formation between the drug and the dye bromothymol blue at pH 3.5. The detection was performed without phase separation in a glass tube coupled to a fiber-optics spectrophotometer. The total volume of chloroform was reduced to 50 mu L in comparison with 10 mL consumed in batch. A linear response was observed between 9 and 120 mu mol L(-1), with a detection limit of 0.9 mu mol L(-1) (99.7% confidence level). The coefficient of variation (n = 10), sampling rate and extraction efficiency were estimated as 0.6%, 78 determinations per hour and 61%, respectively. About 30 mu g of bromothymol blue was consumed and the waste volume was 380 mu L per determination. The results for pharmaceutical samples agreed with those obtained by the reference procedure at the 95% confidence level. (C) 2011 Elsevier B.V. All rights reserved.

A critical evaluation of a flow-cell based on a liquid core waveguide for chemiluminescence measurements

Liquid-core waveguides (LCWs), devices that constrain the emitted radiation minimizing losses during the transport, are an alternative to maximize the amount of detected radiation in luminescence. In this work, the performance of a LCW flow-cell was critically evaluated for chemiluminescence measurements, by using as model the oxidation of luminol by hydrogen peroxide or hypochlorite. An analytical procedure for hypochlorite determination was also developed, with linear response in the range 0.2-3.8 mg/L (2.7-51 mu mol/L), a detection limit estimated as 8 mu g/L (0.64 mu mol/L) at the 99.7% confidence level and luminol consumption of 50 mu g/determination. The coefficients of variation were 3.3% and 1.6% for 0.4 and 1.9 mg/L CIO(-), respectively, with a sampling rate of 164 determinations/h. The procedure was applied to the analysis of Dakin`s solution samples, yielding results in agreement with those obtained by iodometric titration at the 95% confidence level. Copyright (c) 2008 John Wiley & Sons, Ltd.

An improved procedure for flow-based turbidimetric sulphate determination based on a liquid core waveguide and pulsed flows

An improved flow-based procedure is proposed for turbidimetric sulphate determination in waters. The flow system was designed with solenoid micro-pumps in order to improve mixing conditions and minimize reagent consumption as well as waste generation. Stable baselines were observed in view of the pulsed flow characteristic of the systems designed with solenoid micro-pumps, thus making the use of washing solutions unnecessary. The nucleation process was improved by stopping the flow prior to the measurement, thus avoiding the need of sulphate addition. When a 1-cm optical path flow cell was employed, linear response was achieved within 20-200 mg L(-1), described by the equation S = -0.0767 + 0.00438C (mg L(-1)), r = 0.999. The detection limit was estimated as 3 mg L(-1) at the 99.7% confidence level and the coefficient of variation was 2.4% (n = 20). The sampling rate was estimated as 33 determinations per hour. A long pathlength (100-cm) flow cell based on a liquid core waveguide was exploited to increase sensitivity in turbidimetry. Baseline drifts were avoided by a periodical washing step with EDTA in alkaline medium. Linear response was observed within 7-16 mg L(-1), described by the equation S = -0.865 + 0.132C (mg L(-1)), r = 0.999. The detection limit was estimated as 150 mu g L(-1) at the 99.7% confidence level and the coefficient of variation was 3.0% (n = 20). The sampling rate was estimated as 25 determinations per hour. The results obtained for freshwater and rain water samples were in agreement with those achieved by batch turbidimetry at the 95% confidence level. (C) 2008 Elsevier B.V All rights reserved.

A flow-based analytical procedure for salbutamol determination exploiting chemiluminescence in a liquid-core waveguide

Salbutamol is a bronchodilator whose use is restricted due to its anabolic effects. A flow-based procedure for salbutamol determination based on the inhibition of chemiluminescence of the luminol/hypochlorite system was developed. A flow cell constructed with a liquid-core waveguide was employed to constrain the emitted radiation, minimizing losses during transport to detector. Linear response was observed within 2.5 x 10(-6) and 1.0 x 10(-5) mol L-1 with a detection limit estimated as 1 x 10(-7) mol L-1 at the 99.7% confidence level. The coefficient of variation (n = 20), sampling rate, and luminol consumption per determination were estimated as 2.8%, 164 determinations h(-1), and 50 mu g, respectively. Results for pharmaceutical samples were in agreement with those obtained by reference procedures at the 95% confidence level.

Rheo-optical characterization of liquid crystalline acetoxypropylcellulose melt undergoing large shear flow and relaxation after flow cessation

The rheological and structural characteristics of acetoxypropylcellulose (APC) nematic melt are studied at shear rates ranging from 10 s(-1) to 1000 s(-1) which are relevant to extrusion based processes. APC shows a monotonic shear thinning behavior over the range of shear rates tested. The negative extrudate-swell shows a minimum when a critical shear rate (gamma) over dot(c) is reached. For shear rates smaller than (gamma) over dot(c), the flow-induced texture consists of two set of bands aligned parallel and normal to the flow direction. At shear rates larger than (gamma) over dot(c), the flow induced texture is reminiscent of a 2 fluids structure. Close to the shearing walls, domains elongated along the flow direction and stacked along the vorticity are imaged with POM, whereas SALS patterns indicate that the bulk of the sheared APC is made of elliptical domains oriented along the vorticity. No full nematic alignment is achieved at the largest shear rate tested. Below (gamma) over dot(c), the stress relaxation is described by a stretched exponential. Above (gamma) over dot(c), the stress relaxation is described by a fast and a slow process. The latter coincides with the growth of normal bands thicknesses, as the APC texture after flow cessation consists of two types of bands with parallel and normal orientations relative to the flow direction. Both bands thicknesses do not depend on the applied shear rate, in contrast to their orientation. (C) 2015 Elsevier Ltd. All rights reserved.

Extensional flow of nematic liquid crystal under electric field gradient

Systematic asymptotic methods are used to formulate a model for the extensional flow of a thin sheet of nematic liquid crystal. With no external body forces applied, the model is found to be equivalent to the so-called Trouton model for Newtonian sheets (and fi bers), albeit with a modi fied "Trouton ratio". However, with a symmetry-breaking electric field gradient applied, behavior deviates from the Newtonian case, and the sheet can undergo fi nite-time breakup if a suitable destabilizing field is applied. Some simple exact solutions are presented to illustrate the results in certain idealized limits, as well as sample numerical results to the full model equations.

High-performance liquid chromatography of the renal blood flow marker p-aminohippuric acid (PAH) and its metabolite N-acetyl PAH improves PAH clearance measurements.

PAH (N-(4-aminobenzoyl)glycin) clearance measurements have been used for 50 years in clinical research for the determination of renal plasma flow. The quantitation of PAH in plasma or urine is generally performed by colorimetric method after diazotation reaction but the measurements must be corrected for the unspecific residual response observed in blank plasma. We have developed a HPLC method to specifically determine PAH and its metabolite NAc-PAH using a gradient elution ion-pair reversed-phase chromatography with UV detection at 273 and 265 nm, respectively. The separations were performed at room temperature on a ChromCart (125 mmx4 mm I.D.) Nucleosil 100-5 microm C18AB cartridge column, using a gradient elution of MeOH-buffer pH 3.9 1:99-->15:85 over 15 min. The pH 3.9 buffered aqueous solution consisted in a mixture of 375 ml sodium citrate-citric acid solution (21.01 g citric acid and 8.0 g NaOH per liter), added up with 2.7 ml H3PO4 85%, 1.0 g of sodium heptanesulfonate and completed ad 1000 ml with ultrapure water. The N-acetyltransferase activity does not seem to notably affect PAH clearances, although NAc-PAH represents 10.2+/-2.7% of PAH excreted unchanged in 12 healthy subjects. The performance of the HPLC and the colorimetric method have been compared using urine and plasma samples collected from healthy volunteers. Good correlations (r=0.94 and 0.97, for plasma and urine, respectively) are found between the results obtained with both techniques. However, the colorimetric method gives higher concentrations of PAH in urine and lower concentrations in plasma than those determined by HPLC. Hence, both renal (ClR) and systemic (Cls) clearances are systematically higher (35.1 and 17.8%, respectively) with the colorimetric method. The fraction of PAH excreted by the kidney ClR/ClS calculated from HPLC data (n=143) is, as expected, always <1 (mean=0.73+/-0.11), whereas the colorimetric method gives a mean extraction ratio of 0.87+/-0.13 implying some unphysiological values (>1). In conclusion, HPLC not only enables the simultaneous quantitation of PAH and NAc-PAH, but may also provide more accurate and precise PAH clearance measurements.

Gas-liquid mass transfer in bubbly flow: Estimation of mass transfer, bubble size and reactor performance in various applications

Gas-liquid mass transfer is an important issue in the design and operation of many chemical unit operations. Despite its importance, the evaluation of gas-liquid mass transfer is not straightforward due to the complex nature of the phenomena involved. In this thesis gas-liquid mass transfer was evaluated in three different gas-liquid reactors in a traditional way by measuring the volumetric mass transfer coefficient (kLa). The studied reactors were a bubble column with a T-junction two-phase nozzle for gas dispersion, an industrial scale bubble column reactor for the oxidation of tetrahydroanthrahydroquinone and a concurrent downflow structured bed.The main drawback of this approach is that the obtained correlations give only the average volumetric mass transfer coefficient, which is dependent on average conditions. Moreover, the obtained correlations are valid only for the studied geometry and for the chemical system used in the measurements. In principle, a more fundamental approach is to estimate the interfacial area available for mass transfer from bubble size distributions obtained by solution of population balance equations. This approach has been used in this thesis by developing a population balance model for a bubble column together with phenomenological models for bubble breakage and coalescence. The parameters of the bubble breakage rate and coalescence rate models were estimated by comparing the measured and calculated bubble sizes. The coalescence models always have at least one experimental parameter. This is because the bubble coalescence depends on liquid composition in a way which is difficult to evaluate using known physical properties. The coalescence properties of some model solutions were evaluated by measuring the time that a bubble rests at the free liquid-gas interface before coalescing (the so-calledpersistence time or rest time). The measured persistence times range from 10 msup to 15 s depending on the solution. The coalescence was never found to be instantaneous. The bubble oscillates up and down at the interface at least a coupleof times before coalescence takes place. The measured persistence times were compared to coalescence times obtained by parameter fitting using measured bubble size distributions in a bubble column and a bubble column population balance model. For short persistence times, the persistence and coalescence times are in good agreement. For longer persistence times, however, the persistence times are at least an order of magnitude longer than the corresponding coalescence times from parameter fitting. This discrepancy may be attributed to the uncertainties concerning the estimation of energy dissipation rates, collision rates and mechanisms and contact times of the bubbles.

Identification of new transformation products during enzymatic treatment of tetracycline and erythromycin antibiotics at laboratory scale by an on-line turbulent flow liquid-chromatography coupled to a high resolution mass spectrometer LTQ-Orbitrap

This work describes the formation of transformation products (TPs) by the enzymatic degradation at laboratory scale of two highly consumed antibiotics: tetracycline (Tc) and erythromycin (ERY). The analysis of the samples was carried out by a fast and simple method based on the novel configuration of the on-line turbulent flow system coupled to a hybrid linear ion trap – high resolution mass spectrometer. The method was optimized and validated for the complete analysis of ERY, Tc and their transformation products within 10 min without any other sample manipulation. Furthermore, the applicability of the on-line procedure was evaluated for 25 additional antibiotics, covering a wide range of chemical classes in different environmental waters with satisfactory quality parameters. Degradation rates obtained for Tc by laccase enzyme and ERY by EreB esterase enzyme without the presence of mediators were ∼78% and ∼50%, respectively. Concerning the identification of TPs, three suspected compounds for Tc and five of ERY have been proposed. In the case of Tc, the tentative molecular formulas with errors mass within 2 ppm have been based on the hypothesis of dehydroxylation, (bi)demethylation and oxidation of the rings A and C as major reactions. In contrast, the major TP detected for ERY has been identified as the “dehydration ERY-A”, with the same molecular formula of its parent compound. In addition, the evaluation of the antibiotic activity of the samples along the enzymatic treatments showed a decrease around 100% in both cases

Modeling of Liquid-Solid Flow in Industrial Scale

In the present work, liquid-solid flow in industrial scale is modeled using the commercial software of Computational Fluid Dynamics (CFD) ANSYS Fluent 14.5. In literature, there are few studies on liquid-solid flow in industrial scale, but any information about the particular case with modified geometry cannot be found. The aim of this thesis is to describe the strengths and weaknesses of the multiphase models, when a large-scale application is studied within liquid-solid flow, including the boundary-layer characteristics. The results indicate that the selection of the most appropriate multiphase model depends on the flow regime. Thus, careful estimations of the flow regime are recommended to be done before modeling. The computational tool is developed for this purpose during this thesis. The homogeneous multiphase model is valid only for homogeneous suspension, the discrete phase model (DPM) is recommended for homogeneous and heterogeneous suspension where pipe Froude number is greater than 1.0, while the mixture and Eulerian models are able to predict also flow regimes, where pipe Froude number is smaller than 1.0 and particles tend to settle. With increasing material density ratio and decreasing pipe Froude number, the Eulerian model gives the most accurate results, because it does not include simplifications in Navier-Stokes equations like the other models. In addition, the results indicate that the potential location of erosion in the pipe depends on material density ratio. Possible sedimentation of particles can cause erosion and increase pressure drop as well. In the pipe bend, especially secondary flows, perpendicular to the main flow, affect the location of erosion.

A SAXS study of flow alignment of thermotropic liquid crystal mixtures

We report on the capillary flow behaviour of thermotropic liquid crystal mixtures containing 4-n-octyl-4'-cyanobiphenyl (8CB) and 4-n-pentyl-4'-cyanobiphenyl (5CB). The liquid crystal mixtures are studied in the Nematic (N) and Smectic (SA) phases at room temperature. Polarised optical microscopy (POM), rheology and simultaneous X-ray diffraction (XRD)/capillary flow experiments are performed to characterise the system. Polarised optical microscopy reveals a dramatic change in optical texture when the 5CB content is increased from 20 to 30% in the mixtures. X-ray diffraction results show that the system goes through a SA-N phase transition, such that the mixtures are smectic for 10-20% 5CB and nematic for 30-90% 5CB. Smectic mixtures flow with the layers aligned along the flow direction (mesogens perpendicular to flow) while nematic mixtures flow with the mesogens aligned in the flow direction. Simultaneous XRD/shear flow experiments show that the SA-N transition is independent of the flow rate in the range 1-6 ml min-1. The correlation length of the liquid crystal order decreases with increasing 5CB content. Rheology is used to prove that the correlation length behaviour is related to a reduction in the viscosity of the mixture.