CFD simulation of complex phenomena containing suspensions and flow throughporous media

There is an increasing reliance on computers to solve complex engineering problems. This is because computers, in addition to supporting the development and implementation of adequate and clear models, can especially minimize the financial support required. The ability of computers to perform complex calculations at high speed has enabled the creation of highly complex systems to model real-world phenomena. The complexity of the fluid dynamics problem makes it difficult or impossible to solve equations of an object in a flow exactly. Approximate solutions can be obtained by construction and measurement of prototypes placed in a flow, or by use of a numerical simulation. Since usage of prototypes can be prohibitively time-consuming and expensive, many have turned to simulations to provide insight during the engineering process. In this case the simulation setup and parameters can be altered much more easily than one could with a real-world experiment. The objective of this research work is to develop numerical models for different suspensions (fiber suspensions, blood flow through microvessels and branching geometries, and magnetic fluids), and also fluid flow through porous media. The models will have merit as a scientific tool and will also have practical application in industries. Most of the numerical simulations were done by the commercial software, Fluent, and user defined functions were added to apply a multiscale method and magnetic field. The results from simulation of fiber suspension can elucidate the physics behind the break up of a fiber floc, opening the possibility for developing a meaningful numerical model of the fiber flow. The simulation of blood movement from an arteriole through a venule via a capillary showed that the model based on VOF can successfully predict the deformation and flow of RBCs in an arteriole. Furthermore, the result corresponds to the experimental observation illustrates that the RBC is deformed during the movement. The concluding remarks presented, provide a correct methodology and a mathematical and numerical framework for the simulation of blood flows in branching. Analysis of ferrofluids simulations indicate that the magnetic Soret effect can be even higher than the conventional one and its strength depends on the strength of magnetic field, confirmed experimentally by Völker and Odenbach. It was also shown that when a magnetic field is perpendicular to the temperature gradient, there will be additional increase in the heat transfer compared to the cases where the magnetic field is parallel to the temperature gradient. In addition, the statistical evaluation (Taguchi technique) on magnetic fluids showed that the temperature and initial concentration of the magnetic phase exert the maximum and minimum contribution to the thermodiffusion, respectively. In the simulation of flow through porous media, dimensionless pressure drop was studied at different Reynolds numbers, based on pore permeability and interstitial fluid velocity. The obtained results agreed well with the correlation of Macdonald et al. (1979) for the range of actual flow Reynolds studied. Furthermore, calculated results for the dispersion coefficients in the cylinder geometry were found to be in agreement with those of Seymour and Callaghan.

The influence of crystallization conditions on the filtration characteristics of sulphathiazole suspensions

Cooling crystallization is one of the most important purification and separation techniques in the chemical and pharmaceutical industry. The product of the cooling crystallization process is always a suspension that contains both the mother liquor and the product crystals, and therefore the first process step following crystallization is usually solid-liquid separation. The properties of the produced crystals, such as their size and shape, can be affected by modifying the conditions during the crystallization process. The filtration characteristics of solid/liquid suspensions, on the other hand, are strongly influenced by the particle properties, as well as the properties of the liquid phase. It is thus obvious that the effect of the changes made to the crystallization parameters can also be seen in the course of the filtration process. Although the relationship between crystallization and filtration is widely recognized, the number of publications where these unit operations have been considered in the same context seems to be surprisingly small. This thesis explores the influence of different crystallization parameters in an unseeded batch cooling crystallization process on the external appearance of the product crystals and on the pressure filtration characteristics of the obtained product suspensions. Crystallization experiments are performed by crystallizing sulphathiazole (C9H9N3O2S2), which is a wellknown antibiotic agent, from different mixtures of water and n-propanol in an unseeded batch crystallizer. The different crystallization parameters that are studied are the composition of the solvent, the cooling rate during the crystallization experiments carried out by using a constant cooling rate throughout the whole batch, the cooling profile, as well as the mixing intensity during the batch. The obtained crystals are characterized by using an automated image analyzer and the crystals are separated from the solvent through constant pressure batch filtration experiments. Separation characteristics of the suspensions are described by means of average specific cake resistance and average filter cake porosity, and the compressibilities of the cakes are also determined. The results show that fairly large differences can be observed between the size and shape of the crystals, and it is also shown experimentally that the changes in the crystal size and shape have a direct impact on the pressure filtration characteristics of the crystal suspensions. The experimental results are utilized to create a procedure that can be used for estimating the filtration characteristics of solid-liquid suspensions according to the particle size and shape data obtained by image analysis. Multilinear partial least squares regression (N-PLS) models are created between the filtration parameters and the particle size and shape data, and the results presented in this thesis show that relatively obvious correlations can be detected with the obtained models.

Separation of Gas Bubbles from Circulation Waters and Pulp Suspensions in Open Channel Flow

The objective of this thesis was to study the removal of gases from paper mill circulation waters experimentally and to provide data for CFD modeling. Flow and bubble size measurements were carried out in a laboratory scale open gas separation channel. Particle Image Velocimetry (PIV) technique was used to measure the gas and liquid flow fields, while bubble size measurements were conducted using digital imaging technique with back light illumination. Samples of paper machine waters as well as a model solution were used for the experiments. The PIV results show that the gas bubbles near the feed position have the tendency to escape from the circulation channel at a faster rate than those bubbles which are further away from the feed position. This was due to an increased rate of bubble coalescence as a result of the relatively larger bubbles near the feed position. Moreover, a close similarity between the measured slip velocities of the paper mill waters and that of literature values was obtained. It was found that due to dilution of paper mill waters, the observed average bubble size was considerably large as compared to the average bubble sizes in real industrial pulp suspension and circulation waters. Among the studied solutions, the model solution has the highest average drag coefficient value due to its relatively high viscosity. The results were compared to a 2D steady sate CFD simulation model. A standard Euler-Euler k-ε turbulence model was used in the simulations. The channel free surface was modeled as a degassing boundary. From the drag models used in the simulations, the Grace drag model gave velocity fields closest to the experimental values. In general, the results obtained from experiments and CFD simulations are in good qualitative agreement.

Pressure filtration characteristics of enzymatically hydrolyzed biomass suspensions

Enzymatic hydrolysis of lignocellulosic polymers is likely to become one of the key technologies enabling industrial production of liquid biofuels and chemicals from lignocellulosic biomass. Certain types of enzymes are able to hydrolyze cellulose and hemicellulose polymers to shorter units and finally to sugar monomers. These monomeric sugars are environmentally acceptable carbon sources for the production of liquid biofuels, such as bioethanol, and other chemicals, such as organic acids. Liquid biofuels in particular have been shown to contribute to the reduction of net emissions of greenhouse gases. The solid residue of enzymatic hydrolysis is composed mainly of lignin and partially degraded fibers, while the liquid phase contains the produced sugars. It is usually necessary to separate these two phases at some point after the hydrolysis stage. Pressure filtration is an efficient technique for this separation. Solid-liquid separation of biomass suspensions is difficult, because biomass solids are able to retain high amounts of water, which cannot be readily liberated by mechanical separation techniques. Most importantly, the filter cakes formed from biomaterials are compressible, which ultimately means that the separation may not be much improved by increasing the filtration pressure. The use of filter aids can therefore facilitate the filtration significantly. On the other hand, the upstream process conditions have a major influence on the filtration process. This thesis investigates how enzymatic hydrolysis and related process conditions affect the filtration properties of a cardboard suspension. The experimental work consists of pressure filtration and characterization of hydrolysates. The study provides novel information about both issues, as the relationship between enzymatic hydrolysis conditions and subsequent filtration properties has so far not been considered in academic studies. The results of the work reveal that the final degree of hydrolysis is an important factor in the filtration stage. High hydrolysis yield generally increases the average specific cake resistance. Mixing during the hydrolysis stage resulted in undefined changes in the physical properties of the solid residue, causing a high filtration resistance when the mixing intensity was high. Theoretical processing of the mixing data led to an interesting observation: the average specific cake resistance was observed to be linearly proportional to the mixer shear stress. Another finding worth attention is that the size distributions of the solids did not change very dramatically during enzymatic hydrolysis. There was an observable size reduction during the first couple of hours, but after that the size reduction was minimal. Similarly, the size distribution of the suspended solids remained almost constant when the hydrolyzed suspension was subjected to intensive mixing. It was also found that the average specific cake resistance was successfully reduced by the use of filter aids. This reduction depended on the method of how the filter aids were applied. In order to obtain high filtration capacity, it is recommended to use the body feed mode, i.e. to mix the filter aid with the slurry prior to filtration. Regarding the quality of the filtrate, precoat filtration was observed to produce a clear filtrate with negligible suspended solids content, while the body feed filtrates were turbid, irrespective of which type of filter aid was used.