MECHANISM OF CRYSTALLIZATION IN AGITATED SOLUTIONS

A new concept of molecular cluster dynamics in turbulent flows is used to develop a physical model and mechanism of homogeneous and heterogeneous nucleation in solutions. The proposed physical model considers the motion of molecular clusters in centrifugal fields of free eddies induced in turbulent flows. The clusters are moved to an eddy boundary region where they are concentrated and aggregated to form critical nuclei.

It is shown theoretically that a critical fluid velocity exists, below which nuclei are not formed, and that its value depends upon the degree of solution supersaturation. The nucleation theory takes into account the hydrodynamic parameters of the process and permits calculation of the number of nuclei formed. Satisfactory agreement was obtained between calculated values and published experimental data.     

MECHANISM OF CRYSTALLIZATION IN AGITATED SOLUTIONS

A new concept of molecular cluster dynamics in turbulent flows is used to develop a physical model and mechanism of homogeneous and heterogeneous nucleation in solutions. The proposed physical model considers the motion of molecular clusters in centrifugal fields of free eddies induced in turbulent flows. The clusters are moved to an eddy boundary region where they are concentrated and aggregated to form critical nuclei.

It is shown theoretically that a critical fluid velocity exists, below which nuclei are not formed, and that its value depends upon the degree of solution supersaturation. The nucleation theory takes into account the hydrodynamic parameters of the process and permits calculation of the number of nuclei formed. Satisfactory agreement was obtained between calculated values and published experimental data.     

HYDRODYNAMICS OF LAMINAR LIQUID JETS. EXPERIMENTAL STUDY AND COMPARISON WITH TWO MODELS

Laminar jets of Newtonian liquids issuing from long vertical cylindrical nozzles and falling freely through stagnant air were studied experimentally for Reynolds numbers between 300 and 1000. Jet diameters were measured from still photographs, and radial distributions of axial velocity were obtained by laser Doppler anemometry. The effect of nozzle diameter, fluid viscosity and surface tension was investigated.

The experimental results were compared with numerical solutions of the Protean coordinate model developed by Duda and Vrentas. The boundary layer simplifications were confirmed to be valid only for the downstream region of the jet and for Reynolds numbers greater than 1000.

The experimental diameters were also compared with predictions from a form of the Bernoulli equation with a surface tension term. The asymptotic validity of the model was confirmed, provided that the dissipation term arising from fluid viscosity could be neglected.

Neither model correlated the jet formation region satisfactorily. For this region, an empirical correlation was developed which improves the diameter prediction and is complementary of either model.     

GAS MIXING IN SLUGGING AND TURBULENT FLUIDIZED BEDS

The gas phase mixing in a fluidized bed of glass beads (d_p = 0.362 mm) in the slugging and turbulent flow regimes has been studied in a 0.1 m-ID × 3.0 m high Plexiglas column.

The gas dispersion in the downstream of the bed has been described by a diffusion process with the axial and radial dispersion coefficients. The radial dispersion coefficient of the gas phase is nearly constant with the variation of gas velocity in the slugging flow regime, but it increases with an increase in gas velocity in the turbulent flow regime.

Appreciable backmixing of the gas phase is pronounced in the slugging flow regime whereas the lower gas backmixing is produced in the turbulent flow regime. The gas backmixing coefficient increases with an increase in gas velocity in the slugging flow regime, but it decreases slightly with an increase in gas velocity in the turbulent flow regime.

The radial mixing and backmixing coefficients of the gas in terms of Peclet numbers have been correlated with the relevant dimensionless parameters (U_g/U_mf, p_s/p_g, d_p/D_t).

The gas flow pattern in the bed has been well represented by a simplified model based on the two gas phases in the dilute and dense phases which are percolating through the bed in plug flow. The present model can predict the gas exchange coefficient between the phases, the fractions of the dilute phase, the interstitial gas in the dense phase, and the interstitial gas velocity in the bed.     

STEAM CONDENSATION IN A VERTICAL CORRUGATED DUCT

An experimental investigation was made on the condensation of water steam in a vertical corrugated duct. The data have been correlated as follows

Co=5.11Re^-0.431 150 ≤ Re ≤ 350 Co = 0.034311Re^-0.425350 ≤Re ≤l000

The vertical corrugated duct is constructed of two corrugated plates with corrugation inclination angles of β = 0 and β = 45° respectively (relative to the overall flow direction).

The condensation heat transfer coefficient in the corrugated duct is more than two times higher than that of bulk condensation on a vertical plate. A physical model was proposed to explain the heat transfer enhancement. Attention was also paid to the effect of exit steam velocity on the heat transfer during partial condensation. It was demonstrated that the heat transfer in the corrugated duct was strongly affected even at a low exit velocity, which was different from the case of bulk condensation on a vertical plate. Experimental apparatuses and the method for examining their reliability are described in detail.     

DROP DISPERSION IN SUSPENSION POLYMERIZATION

Four models, two based on laminar shear and two based on turbulent flow, are proposed to describe drop dispersion in non-coalescing systems. The models predict the largest surviving drop size _dmax as a function of geometry, speed and physical property variables.

Laboratory data including suspension polymerization runs support the boundary layer laminar shear model for drops larger than approximately 200 microns. Smaller drops support a turbulence model.

The boundary layer shear model was confirmed in scale-up suspension polymerization runs aimed at producing 1000 micron maximum bead sizes. Five approximately geometrically similar polymerizers were used, varying in size from 7.5 to 15000 liters.     

MASS TRANSFER THROUGH PERMEABLE WALLS: A NEW INTEGRAL EQUATION APPROACH FOR CYLINDRICAL TUBES WITH LAMINAR FLOW

Mass transfer through cylindrical semipermeable walls is analyzed. The solution is obtained in terms of integral equations. Despite the existence of a non-homogeneous boundary condition on the semipermeable wall, the solution thus obtained is particularly advantageous since the associated eigenvalue problem is independent of the Sherwood number. This parameter takes into account the main conductances at the tube wall.

The approach is applied to the case of mass transfer from the interior of a capillary tube with semipermeable walls to an external fluid. The flow in the tube is laminar, and the external flow is assumed turbulent.

The mathematical methodology employed provides a framework to develop numerical schemes of fast and sure convergence.     

LIQUID CIRCULATION PATTERNS AND THEIR EFFECT ON GAS HOLD-UP AND AXIAL MIXING IN BUBBLE COLUMNS

Liquid velocity profiles in a bubble column were measured with the aid of a hot-film anemometer. Two different types of profiles were detected that were characterized by the formation of a boundary layer at the column wall. Under slow flow conditions, the boundary layer is large and controls the liquid velocity profile, but in turbulent flow the boundary layer is very small and has no significant effect on the velocity profile. Microscopic and macroscopic balances were used to predict both the liquid velocity profile and the average liquid velocity.

The effect of the liquid velocity profile and average velocity on axial dispersion and gas-hold-up are analyzed and design procedures are recommended.     

AN ANALYSIS OF CELLULAR FLOW OBSERVED ON THE INTERFACE OF A PARTIALLY FORMED DROPLET†

In [16], during an experiment designed to model the internal circulation of a forming droplet, secondary surface flows were observed on the droplet interface.

After summarizing the experimental results of [16], we present one possible mechanism, based on the surface surfactant mass transport equation of Levich and the surface stress-strain boundary conditions at a free surface, that provides a good qualitative explanation of the origins and the nature of the secondary motion observed in [16]. The critical hypotheses in this mechanism are that the normal component of ihe vorticity at the free surface is determined primarily by the components of the velocity field tangential to the level lines of the surface surfactant density, near the maxima and minima of that density function and that the normal component of the fluid stress does not vanish at such points.

The consequent analysis of the mass transport equation in the interface shows that the resulting surface motion may be viewed as arising from a resonance phenomenon analogous to the forced vibrations of a spring at resonance.

The effects of adsorbtion-desorbtion and surface dilational viscosity may be incorporated in this mechanism. A method for the experimental measurement of surface dilational viscosities is proposed.     

APPLICATION OF TURBULENCE FUNDAMENTALS TO REACTOR MODELLING AND SCALEUP

Applications of the fundamentals of turbulent mixing become clear once those fundamentals are understood. The first article in this series presented those fundamentals, in order to show how to apply turbulent mixing fundamentals modelling and scaleup, this article covers the following topics:

1. reaction types and their interaction with mixing;

2. closure of the Reynolds equations for mixing and reactions;

3. application to complex geometries;

4. random coalescence-dispersion modelling;

5. application to complex chemistry.

The most difficult problem in applying our knowledge of turbulence to mixer modelling and scaleup is the choice of model complexity. The levels of model complexity available and how to apply them to various problems are presented following the introduction.     

A MONTE-CARLO APPROACH TO THE INTERPRETATION OF BUBBLE PROBE SIGNALS IN FREELY BUBBLING FLUID BEDS†

The interpretation of the signals generated by a double probe may be done by the three characteristic times method:

t₁ the time duration of the pulses;

t₂ the time shift between the signals of the two probes;

t₃ the time interval between two pulses on one channel.

Each of these times is largely dispersed and the corresponding histograms may be constructed.

The present work is a trial to go over from the time histograms to physical properties of the bubbles combining a Monte-Carlo simulation and a flexible simplex optimisation procedure.

As a result, the percentage of oblique bubbles cutting just one level, the bubble size distribution, the average velocity-size relation, and the individual dispersion around it may be defined.

The procedure is finally applied to experimental results obtained with a light probe in a fluid bed of glass beads.     

MULTIPHASE HYDRODYNAMICS AND SOLID-LIQUID MASS TRANSPORT IN AN EXTERNAL-LOOP AIRLIFT REACTOR—A COMPARATIVE STUDY

The solid-solid mass transfer performance of an external-loop airlift reactor was measured by dissolution of benzoic acid coated on nylon-6 particles, and the hydrodynamics of the gas-liquid-solid multiphase system in the airlift reactor were investigated. The solid-liquid system was designed to simulate the micro-carrier culture of animal cells, and some typical suspensions of immobilized enzyme particles.

The solid-liquid mass transfer coefficient remained constant below a superficial air velocity of 0.04 ms^-1 for the particles examined, but increased rapidly with further increase in gas velocity. Solids loading (0.3-3.5% w/w) did not affect the mass transfer coefficient in turbulent flow.

The mass transfer coefficient was correlated with energy dissipation rate in the airlift reactor. The mass transfer coefficient in stirred vessels, bubble columns, fluidized beds, and airlift reactors was compared.

Over an energy dissipation Reynolds number of 4-400, the solid-liquid mass transfer coefficient in the airlift device was comparable to that obtainable in fluidized beds. The performance of the airlift was distinctly superior to that of bubble columns and stirred tanks.     

PHASE HOLDUPS AND LIQUID-LIQUID EXTRACTION IN THREE PHASE FLUIDIZED BEDS

Effects of the continuous phase velocity (0.01-0.08 m/s(, the dispersed phase velocity (0.0-0.04 m/s) and particle size (1.0-3.0 mm) on the individual phase holdups and the mass transfer coefficient have been determined in two (liquid-liquid) and three (liquid-liquid-solid) phase fluidized beds.

In the beds, the dispersed phase holdup increased with dispersed phase velocity but it decreased with continuous phase velocity. Whereas the continuous phase holdup decreased with dispersed phase velocity but it increased with continuous phase velocity. The bed porosity increased with both the dispersed and continuous phase velocities in the beds of 1.7 and 3.0 mm particles. In addition, the continuous phase holdup decreased with the presence of solid particles in the bed, however, the dispersed phase holdup was not affected by the presence of the particles.

The overall mass transfer coefficients in the continuous and dispersed phases increased with increasing fluid velocities but it decreased with the bed height.

The continuous phase holdup and mass transfer coefficient data have been correlated with the operating variables and the dimensionless groups.     

NUMERICAL SIMULATION OF TURBULENT FLOW OF AGITATED LIQUID WITH PITCHED BLADE IMPELLER

The turbulent flow field in an agitated system with baffles was solved numerically using the standard k-e model, an algebraic Reynolds stress model (ASM) and a differential Reynolds stress model (RSM). The commercial software FLOW3D (CFDS, Harwell Laboratories, 1991) was used for this purpose. The aim of the study was to investigate the influence of the impeller boundary conditions and turbulence models to the agreement with experimentally obtained laser-Doppler anemometry data. The boundary conditions for the impeller discharge used in the numerical calculations were obtained as whole-cycle-ensemble averages from experimental LDA-measurements (Fort et al., 1992). Since measurements of the rate of dissipation of turbulent kinetic energy ( ε) was not available the dissipation rate per unit mass in the impeller discharge was estimated from the expression:

where k is the turbulent kinetic energy per unit mass and L the macroscale of turbulence in the pitched blade impeller discharge. The macroscale of turbulence (L) in the impeller boundary condition for e was varied in order to optimize the fit of theoretically obtained profiles of turbulent kinetic energy with experimental data. The constant A was fixed to 0.85 according to Wu and Patterson (1989). The optimal values of L for the different turbulence models were compared with the projected height of the impeller blade (h). All three components of the mean velocity were compared with experimental data for the optimal ratio of L/h for six radial cross-sections in the tank.

The mean velocity field obtained from simulations showed good agreement with experimental data for all models, with somewhat better agreement for the k — e model. An optimal value of the ratio L/h was found to be equal to 2.0 for the k — ε model and 1.3 for the ASM. However, no such optimal value for the RSM could be determined in this study.     

Drying Sawdust in a Pulsed Fluidized Bed

The objective of this work was the experimental and theoretical study of sawdust drying, in batch and continuous experiences, using a pulsed fluidized bed dryer.

In the batch experiences, a 2³ factorial design was used to determine the kinetics of drying, the critical moisture content, and the effective coefficients of both diffusivity and heat transfer, all of them as a function of the velocity and temperature of the air, the speed of turning of the slotted plate that generates the air pulses in the dryer, using sawdust with 65% moisture in each run.

In the continuous operation, a 2³ factorial design was used to study the effect of the solid flow and the velocity and temperature of the air on both the product moisture and the distribution of residence times. In order to determine these last ones, digital image processing was used, utilizing sawdust colored by a solution of methylene blue as tracer.

The statistically significant factors were the velocity and the temperature of the heating air, for both the continuous and batch operations. Although the speed of turn of the slotted plate was not significant, it was observed that the air pulses increased the movement of particles, facilitating its fluidization, especially at the beginning of drying.

The heat transfer coefficients were adjusted according to the equation Nu = 0.0014 Re_p^1.52, whose standard deviation of fit is 0.145.

The period of decreasing rate was adjusted to several diffusivity models, giving the best fit the simplified variable diffusivity model (SVDM). The curve of distribution of residence times was adjusted using the model of tanks in series, with values between 2.6 and 5 tanks.     

EFFECTS OF GEOMETRY ON HYDRODYNAMICS IN EXTERNAL-LOOP AIRLIFT REACTORS

Experimental investigations were carried out in model external-loop airlift reactors. Two reactors of laboratory scale (riser liquid height ranged between 1.16-1.56 m, riser diameter 0.03 m, A_D/A_R ratio between 0.111-1,000, total liquid volume V_T = (1.189-2.446).10^-3m³) and pilot-plant scale (riser liquid height of 4.4 and 4.7 m, respectively, riser diameter 0.200 m, A_D/A_R ratio of 0.1225 and 0.040 m, total liquid volume, V_T = (0.144-0.170) m³) were used.

The influences of reactor geometry characterized by some parameter as: A_D/A_R ratio, liquid height in riser and downcomer and liquid height in gas separator, together with the amount of introduced air, on the basic hydrodynamic design parameters: gas holdup and liquid circulation velocity were analysed.

The influence of gas sparger design on gas holdup and liquid velocity was found to be negligible.

The experimental liquid circulation velocity was correlated using a simplified form of the energy balance in airlift reactors, valid for external-loop airlift reactors with almost complete phase separation at the top.

An original dimensionless correlation for gas holdup prediction involving superficial velocities of gas and liquid, cross sectional areas, dispersion height, riser diameter, as well as Froude number, was obtained.     

Welding Mechanisms of Plastics: A Review

Strength of welded joints is a function of technological parameters of the production process. The type of function is dependent on the welding mechanism. Different mechanisms were found under various welding conditions. The processes included in the plastic welding mechanism are divided into two groups:

1) Processes which realize the joining of the parts.

2) Processes which create conditions for the first group to proceed. The first series of processes includes:

a) diffusion of macroradicals, molecular segments or molecules of the polymer which can be either in a solid, melted or dissolved state.

b) convective mass transfer.

c) recombination of macroradicals across the contact surface.

d) physical (surface) interaction.

e) any combination of processes described above.

The second group contains:

a) formation of the real contact surface.

b) formation of the macroradicals.

c) destruction and removal of inert layers which prevent real contact of active material.

Each process and the conditions of its proceeding are discussed individually.     

TAYLOR DISPERSION OF MULTICOMPONENT SOLUTES

Coupled transport of multicomponent solutes in globally continuous systems is considered in the framework of the Generalized Taylor dispersion theory. Coupling between transports of n different species at the local (or micro-) scale, is considered to result from first-order irreversible surface reactions occurring on the local space boundaries, or from the off-diagonal terms of the solute diffusivity matrices.

General expressions are obtained for the global effective (long-time) solute dispersion matrix cofficients: mean global scalar reactivity, velocity vector and dispersivity dyadic.

The effect of surface chemical reactions is to partition the matter between different solute constituents. This is manifested in a coupling of the global transport coefficients, which may be mathematically removed by a linear (canonic) transformation applied to the effective global transport equation. This type of coupling does not exist for inert solutes.

The second type of the global coupling is represented by the off-diagonal terms of the global velocity and dispersivity matrices. It exists for both reactive and inert solutes. This coupling stems from the convective dispersion process (dependence or the global velocity vector on the local space coordinate). Is shown to be irremovable from the global transport equation by any linear transformation via the solute partition matrix. In the canonic form of the global equation the irremovable coupling is manifested by the traceless parts of the global solute velocity matrix and the global solute dispersivity.

The solution scheme is illustrated by calculating the mean global diffusivity of a solute consisting of two components, transport of which is coupled at the microscale via the molecular diffusivity matrix. At the macroscale the coupling is shown to be represented by negative off-diagonal terms of the global diffusivity matrix,