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.

Statistical Optimum Design of Heat Exchangers

The optimal design of a heat exchanger system is based on given model parameters together with given standard ranges for machine design variables. The goals set for minimizing the Life Cycle Cost (LCC) function which represents the price of the saved energy, for maximizing the momentary heat recovery output with given constraints satisfied and taking into account the uncertainty in the models were successfully done. Nondominated Sorting Genetic Algorithm II (NSGA-II) for the design optimization of a system is presented and implemented inMatlab environment. Markov ChainMonte Carlo (MCMC) methods are also used to take into account the uncertainty in themodels. Results show that the price of saved energy can be optimized. A wet heat exchanger is found to be more efficient and beneficial than a dry heat exchanger even though its construction is expensive (160 EUR/m2) compared to the construction of a dry heat exchanger (50 EUR/m2). It has been found that the longer lifetime weights higher CAPEX and lower OPEX and vice versa, and the effect of the uncertainty in the models has been identified in a simplified case of minimizing the area of a dry heat exchanger.

Heat Transfer in Plate Heat Exchangers

This study illustrates the different types of plate heat exchangers that are commonly used in various domestic and industrial applications. The main purpose of this paper was to devise a methodology that is capable of calculating optimum number of plates in the design of a plate heat exchanger. To obtain the appropriate number of plates, typically several iterations must be made before a final acceptable design is completed, since plate amount depends on many factors such as, flow velocities, physical properties of the streams, flow channel geometry, allowable pressure drop, plate dimensions, and the gap between the plates. The methodology presented here can be used as a general guide for designing a plate heat exchanger. To investigate the effects of relevant parameters on the thermal-hydraulic design of a plate heat exchanger, several experiments were carried out for single-phase and counter flow arrangement with two brazed plate heat exchangers by varying the flow rates and the inlet temperatures of the fluid streams. The actual heat transfer coefficients obtained based on the experiment were nearly close to the calculated values and to improve the design, a correction factor was introduced. Besides, the effect of flow channel velocity on the pressure drop inside the unit is presented.

Offdesign operation of kraft recovery boiler

A model to solve heat and mass balances during the offdesign load calculations was created. These equations are complex and nonlinear. The main new ideas used in the created offdesign model of a kraft recovery boiler are the use of heat flows as torn iteration variables instead of the current practice of using the mass flows, vectorizing equation solving, thus speeding up the process, using non dimensional variables for solving the multiple heat transfer surface problem and using a new procedure for calculating pressure losses. Recovery boiler heat and mass balances are reduced to vector form. It is shown that these vectorized equations can be solved virtually without iteration. The iteration speed is enhanced by the use of the derived method of calculating multiple heat transfer surfaces simultaneously. To achieve this quick convergence the heat flows were used as the torn iteration parameters. A new method to handle pressure loss calculations with linearization was presented. This method enabled less time to be spent calculating pressure losses. The derived vector representation of the steam generator was used to calculate offdesign operation parameters for a 3000 tds/d example recovery boiler. The model was used to study recovery boiler part load operation and the effect of the black liquor dry solids increase on recovery boiler dimensioning. Heat flows to surface elements for part load calculations can be closely approximated with a previously defined exponent function. The exponential method can be used for the prediction of fouling in kraft recovery boilers. For similar furnaces the firing of 80 % dry solids liquor produces lower hearth heat release rate than the 65 % dry solids liquor if we fire at constant steam flow. The furnace outlet temperatures show that capacity increase with firing rate increase produces higher loadings than capacity increase with dry solids increase. The economizers, boiler banks and furnaces can be dimensioned smaller if we increase the black liquor dry solids content. The main problem with increased black liquor dry solids content is the decrease in the heat available to superheat. Whenever possible the furnace exit temperature should be increased by decreasing the furnace height. The increase in the furnace exit temperature is usually opposed because of fear of increased corrosion.

Heat transfer in superheaters of the bubbling fluidized bed boiler conversions

It is often reasonable to convert old boiler to bubbling fluidized bed boiler instead of building a new one. Converted boiler consists of old and new heat surfaces which must be fitted to operate together. Prediction of heat transfer in not so ideal conditions sets challenges for designers. Two converted boilers situated in Poland were studied on the grounds of acceptance tests and further studies. Calculation of boiler process was performed with boiler design program. Main interest was heat transfer in superheaters and factors affecting it. Theory for heat transfer is presented according to information found from literature. Results obtained from experimental studies and calculations have been compared. With correct definitions calculated parameters corresponded well to measured data at boiler maximum design load. However overload situations revealed to be difficult to model at least without considering changes in the combustion process which requires readjustments to the design program input values.

Removal of chromium (vi) in aqueous environments using cork and heat-treated cork samples from quercus cerris and quercus suber

Chromium (VI) removal and its reduction to chromium (III) from aqueous solution by untreated and heat-treated Quercus cerris and heat-treated Quercus suber black agglomerate cork granules was investigated. Initial screening studies revealed that among the sorbents tested, untreated Q. cerris and Q. suber black agglomerate are the most efficient in the removal of Cr(VI) ions and were selected for adsorption essays. Heat treatment adversely affected chromium adsorption and chromium (VI) reduction in Q. cerris cork. The highest metal uptake was found at pH 3.0 for Q. cerris and pH 2.0 for black agglomerate. The experimental data fitted the Langmuir model and the calculated qmax was 22.98 mg/g in black agglomerate and 21.69 mg/g in untreated Q. cerris cork. The FTIR results indicated that while in black agglomerate, lignin is the sole component responsible for Cr(VI) sorption, and in untreated Q. cerris cork, suberin and polysaccharides also play a significant role on the sorption. The SEM-EDX results imply that chromium has a homogenous distribution within both cork granules. Also, phloemic residues in Q. cerris granules showed higher chromium concentration. The results obtained in this study show that untreated Q. cerris and black agglomerate cork granules can be an effective and economical alternative to more costly materials for the treatment of liquid wastes containing chromium

A three dimensional study of heat transfer in a supercritical circulating fluidized bed combustor

From the boiler design point of view, it is imperative to know and understand the operation of the boiler. Since comprehensive measurement of a large furnace is impossible, the furnace can be modeled in order to study its behavior and phenomena. This requires the used model to be validated to correspond with the physical furnace behavior. In this thesis, a three dimensional furnace model is validated to match a bituminous coal utilizing, supercritical once-through circulating fluidized bed combustor based on measurement data. The validated model is used for analyzing the furnace heat transfer. Other heat transfer analysis methods are energy balance method based on tube surface temperature measurements and a method based on measured temperature difference between the tube crest and the fin. The latter method was developed in the thesis using Fluent-software. In the theory part, literature is reviewed and the fundamental aspects of circulating fluidized bed are discussed. These aspects are solid particle behavior in fluidization known as hydrodynamics, behavior of fuel and combustion and heat transfer. Fundamental aspects of modeling are also presented.

Statistical Analysis and Numerics of Heat Exchanger Models

The identifiability of the parameters of a heat exchanger model without phase change was studied in this Master’s thesis using synthetically made data. A fast, two-step Markov chain Monte Carlo method (MCMC) was tested with a couple of case studies and a heat exchanger model. The two-step MCMC-method worked well and decreased the computation time compared to the traditional MCMC-method. The effect of measurement accuracy of certain control variables to the identifiability of parameters was also studied. The accuracy used did not seem to have a remarkable effect to the identifiability of parameters. The use of the posterior distribution of parameters in different heat exchanger geometries was studied. It would be computationally most efficient to use the same posterior distribution among different geometries in the optimisation of heat exchanger networks. According to the results, this was possible in the case when the frontal surface areas were the same among different geometries. In the other cases the same posterior distribution can be used for optimisation too, but that will give a wider predictive distribution as a result. For condensing surface heat exchangers the numerical stability of the simulation model was studied. As a result, a stable algorithm was developed.

Development of Color Difference Equations Matching to Human Vision System

Research on color difference evaluation has been active in recent thirty years. Several color difference formulas were developed for industrial applications. The aims of this thesis are to develop the color density which is denoted by comb g and to propose the color density based chromaticity difference formulas. Color density is derived from the discrimination ellipse parameters and color positions in the xy , xyY and CIELAB color spaces, and the color based chromaticity difference formulas are compared with the line element formulas and CIE 2000 color difference formulas. As a result of the thesis, color density represents the perceived color difference accurately, and it could be used to characterize a color by the attribute of perceived color difference from this color.

Moisture sorption isotherms and isosteric heat determination in Chilean papaya (Vasconcellea pubescens)

The moisture sorption isotherms of Chilean papaya were determined at 5, 20, and 45 ºC, over a relative humidity range of 10-95%. The GAB, BET, Oswin, Halsey, Henderson, Smith, Caurie and Iglesias-Chirife models were applied to the sorption experimental data. The goodness of fit of the mathematical models was statistically evaluated by means of the determination coefficient, mean relative percentage deviation, sum square error, root-mean-square error, and chi-square values. The GAB, Oswin and Halsey models were found to be the most suitable for the description of the sorption data. The sorption heats calculated using the Clausius-Clapeyron equation were 57.35 and 59.98 kJ·mol-1, for adsorption and desorption isotherms, respectively.

A Theoretical investigation of a potential high energy density compound 3,6,7,8-tetranitro-3,6,7,8-tetraaza-tricyclo[3.1.1.1(2,4)]octane

The B3LYP/6-31G (d) density functional theory (DFT) method was used to study molecular geometry, electronic structure, infrared spectrum (IR) and thermodynamic properties. Heat of formation (HOF) and calculated density were estimated to evaluate detonation properties using Kamlet-Jacobs equations. Thermal stability of 3,6,7,8-tetranitro-3,6,7,8-tetraaza-tricyclo [3.1.1.1(2,4)]octane (TTTO) was investigated by calculating bond dissociation energy (BDE) at the unrestricted B3LYP/6-31G(d) level. Results showed the N-NO2 bond is a trigger bond during the thermolysis initiation process. The crystal structure obtained by molecular mechanics (MM) methods belongs to P2(1)/C space group, with cell parameters a = 8.239 Å, b = 8.079 Å, c = 16.860 Å, Z = 4 and r = 1.922 g cm-3. Both detonation velocity of 9.79 km s-1 and detonation pressure of 44.22 GPa performed similarly to CL-20. According to the quantitative standards of energetics and stability, TTTO essentially satisfies this requirement as a high energy density compound (HEDC).

Influence of Heat Treatment on Strength Properties of Finger-jointed Components

The purpose for the thesis was to study the thermo treatment of finger-jointed wood. The thesis concentrated on examining the tensile and bending strength of finger-jointed and thermo treated wood. The aim was to find out how different treatment temperature levels and adhesives influence the strength of wood that has been finger-jointed before heat treatment. Secondary objectives were to examine the influence of the treatment time at one temperature, determine differences in the strength between the joints in heartwood and sapwood, examine the visual appearance of the finger joints after the treatment and establish possibilities to reach a characteristic strength level corresponding to C14. Only minor differences in strength properties were measured between the finger-jointed wood treatments II and III. A greater difference was shown between these two treatment temperatures I, which lead to reduced strength. The average strength of joints glued with adhesive 2 was higher after treatments II and III compared to those glued with the adhesive 1. At the treatment temperature I, the adhesive 1 strength properties were at the same level compared to the adhesive 2 or better. There were not any significant differences.