Academic literature on the topic 'Computational methods in fluid flow, heat and mass transfer (incl. computational fluid dynamics)'

Computational Fluid Dynamics (CFD) has numerous applications in the field of energy research, in modelling the basic physics of combustion, multiphase flow and heat transfer; and in the simulation of mechanical devices such as turbines, wind wave and tidal devices, and other devices for energy generation. With the constant increase in available computing power, the fidelity and accuracy of CFD simulations have constantly improved, and the technique is now an integral part of research and development. In the past few years, the development of multiscale methods has emerged as a topic of intensive research. The variable scales may be associated with scales of turbulence, or other physical processes which operate across a range of different scales, and often lead to spatial and temporal scales crossing the boundaries of continuum and molecular mechanics. In this paper, we present a short review of multiscale CFD frameworks with potential applications to energy problems.

Flow, heat, and mass transfer in fixed beds of catalyst particles are complex phenomena and, when combined with catalytic reactions, are multiscale in both time and space; therefore, advanced computational techniques are being applied to fixed bed modeling to an ever-greater extent. The fast-growing literature on the use of computational fluid dynamics (CFD) in fixed bed design reflects the rapid development of this subfield of reactor modeling. We identify recent trends and research directions in which successful methodology has been established, for example, in computer generation of packings of complex particles, and where more work is needed, for example, in the meshing of nonsphere packings and the simulation of industrial-size packed tubes. Development of fixed bed reactor models, by either using CFD directly or obtaining insight, closures, and parameters for engineering models from simulations, will increase confidence in using these methods for design along with, or instead of, expensive pilot-scale experiments.

This paper focuses on the Computational Fluid Dynamics (CFD) studies on one of the commonly used drying processes for different applications. First, a brief information about drying is given with determining important properties that effect drying characteristics. Next, basic principles of CFD modelling are explained while capabilities of computational processing are presented. A detailed literature survey about CFD studies in convective drying process is then conducted. Finally, some sound concluding remarks are listed. It may be concluded that the CFD is a powerful and flexible tool that can be adopted to many different physical situations including complex scenarios, results of CFD simulations represent good predictions for fluid flow, heat and mass transfer of various drying methods and those numerical studies can be used for validation and controlling of applicability of new drying systems.

Computational fluid dynamics was born principally in the aerospace field as a method for fluid flow and heat transfer research methods following experimental and analytical approaches. Along with progress in the cost performance of computers, computational fluid dynamics is now establishing itself as a tool to improve production processes and product quality in the steel, nonferrous metals, glass, plastics, and composite materials industries.Materials manufacturers use computational fluid dynamics for diverse purposes:1. Reduction in experimental conditions and costs;2. Detailed analysis of mechanisms with multifaceted information unobtainable through experimentation;3. Universal tool for scale-up; and4. Evaluation of novel processes.It can be readily imagined that accuracy, flexibility, and other requirements of computational fluid dynamics should vary with specific applications.Fluids generally observed in materials manufacturing processes are molten materials such as metal, glass, and plastics, and gases for stirring and refining. In the flow of such fluids, materials quality and process characteristics are governed by the following:1. Transport phenomena in the bulk region (where fluid flow is normally turbulent);2. Chemical reaction at interfaces;3. Transport phenomena in boundary layers near the interfaces; and4. Complex coupled phenomena (heat transfer, diffusion, chemical reaction, phase transformation like solidification, free surface, electromagnetic force, and bubble flow).

The heat transfer in annular heat exchanger with titanium oxide of 1.0 volume % concentration as the medium of heat exchanger is considered in this study. The heat transfer simulation of the flow is performed by using Computational Fluid Dynamics package, Ansys Fluent. The heat transfer coefficients of water to titanium oxide nanofluid flowing in a horizontal counter-flow heat exchanger under turbulent flow conditions are investigated. The results show that the convective heat transfer coefficient of the nanofluid is slightly higher than that of the base fluid by several percents. The heat transfer coefficient increases with the increase of the mass flow rate of hot water and also the nanofluid.

Stirred tanks are prevalent in various industries, including chemical, biochemical, and pharmaceutical industries. These reactors are suitable for ensuring efficient mass and heat transfer because adequate mixing can be achieved. Numerous studies have been conducted on small-scale stirred-tank reactors. However, upscaling such reactors is challenging because of the complex flow behavior inside the system, especially for the mixing of immiscible liquid–liquid systems. Thus, the objectives of this study were to examine the flow behavior and upscale an immiscible liquid–liquid stirred tank using CFD simulation by investigating a flat-bottomed stirred tank reactor, equipped with a six-blade Rushton turbine. The simulated results were in good agreement with those obtained experimentally. The scale of the reactor significantly affects the hydrodynamic behavior, and the uniformity of the radial distribution of the velocity decreases with increasing Reynolds number. Furthermore, the upscaling criteria were evaluated for geometric similarity and equal mixing times. The proposed scaling law reliably scaled up the immiscible liquid–liquid mixing in a stirred tank with a difference in the range of ±10%.

The Nusselt number for a Shell and tube Heat Exchanger with segmental baffles for different nanofluids, for different mass flow rate are discussed in the present paper. A shell and tube heat exchanger with 7 tubes and 4 segmental baffles modelling is done using SOLIDWORKS and simulation is done by the Computational Fluid Dynamic (CFD) software; ANSYS-FLUENT. By using Fluent, computational fluid dynamics software the heat transfer coefficient and various heat characteristics of Al2O3–H2O, TiO2–H2O and CuO–H2O for 1% volume of concentration nanofluids are estimated in the Shell and Tube Heat Exchanger considering the turbulent flow.

A detailed three-dimensional computational fluid dynamics (CFD) analysis on gas-particle flow and heat transfer inside a solid-particle solar receiver, which utilizes free-falling particles for direct absorption of concentrated solar radiation, is presented. The two-way coupled Euler-Lagrange method is implemented and includes the exchange of heat and momentum between the gas phase and solid particles. A two-band discrete ordinate method is included to investigate radiation heat transfer within the particle cloud and between the cloud and the internal surfaces of the receiver. The direct illumination energy source that results from incident solar radiation was predicted by a solar load model using a solar ray-tracing algorithm. Two kinds of solid-particle receivers, each having a different exit condition for the solid particles, are modeled to evaluate the thermal performance of the receiver. Parametric studies, where the particle size and mass flow rate are varied, are made to determine the optimal operating conditions. The results also include detailed information for the gas velocity, temperature, particle solid volume fraction, particle outlet temperature, and cavity efficiency.

In this numerical study, the heat transfer performance of shell-and-tube heat exchangers (STHXs) has been compared for two different tube arrangements. STHX having 21 and 24 tubes arranged in the inline and staggered grid has been considered for heat transfer analysis. Shell-and-tube heat exchanger with staggered grid arrangement has been observed to provide lesser thermal stratification as compared to the inline arrangement. Further, the study of variation in the mass flow rate of shell-side fluid having constant tube-side flow rate has been conducted for staggered grid structure STHX. The mass flow rate for the shell side has been varied from 0.1 kg/s to 0.5 kg/s, respectively, keeping the tube-side mass flow rate as constant at 0.25 kg/s. The influence of bulk mass-influx transfer rate on heat transfer efficiency, effectiveness, and pressure drop of shell-tube heat exchangers has been analyzed. CFD results were compared with analytical solutions, and it shows a good agreement between them. It has been observed that pressure drop is minimum for the flow rate of 0.1 kg/s, and outlet temperatures at the shell side and tube side have been predicted to be 40.94°C and 63.63°C, respectively.

Thesis (Ph. D.)--State University of New York at Binghamton, Thomas J. Watson School of Engineeering and Applied Science, Department of Mechanical Engineering, 2009.
Includes bibliographical references.

Electrostatic Precipitators (ESPs) are widely used to capture particulate matter from flue gas. In coal-based power stations, they are used for capturing fly ash before the flue gas is released to the environment. Coal-based power plants are still one of the major suppliers of energy because they are more reliable and have lower unit cost of power generation. Under the current environmental protection regulation controlled by the Environment Protection Agency (EPA), only the finer particles can be released to the environment. However, this is likely to change and coal-based power plants will then have to face stricter rules about permissible size limits for particulate matter (PM) discharged in flue gas, namely the particle size of PM 2.5 (micron) or less. It is therefore required that the capabilities of ESPs are enhanced so that they will be able to capture these finer particles. The main aim of the research is to investigate the micro-size particulate matter capture ability of existing ESPs and determine the operational parameter relationships to improve the collection efficiency of ESPs. In particular, this research focuses specifically on the flow phenomena of the finer fly ash particles inside the ESP model and how they are impacted by the changed geometries and varied electric fields. This involves studying the flow velocity and forces associated with the flow and the electric field and the relevant parameters affecting the dust collection and thus establishing and validating a relationship between the interaction of two phase flow and electric field to reveal the underlying physics for collecting finer particles. To achieve the aim, a laboratory scale ESP was constructed for undertaking various tests and measurements using a novel method. This method involved flow measurement in the ESP chamber using a pitot tube and a cobra tube, whilst employing different shaped baffles in the chamber, varying production of electrostatic field in the ESP model and testing its capturing capacity. This research investigated the influence of internal geometry of the ESP on the flow in the ESP chamber. Two different shaped baffles – semicircular and arrow shaped - were designed, fabricated and inserted in the ESP chamber to effect changes to the flue gas pathway to enhance collection efficiency and collection capability of submicron particles. The flow measurements and experimental results were compared and validated with the 2D ii simulation results. Results using baffles indicate that internal geometry of the ESP has an influence on collection efficiency and changing the internal shape produces swirling flow inside ESP, which, in turn, improves collection efficiency. In addition, baffles increase residence time, which allows capture of sub-micron particles. A high voltage transformer and associated electrode plates and rods were designed, constructed and fitted into the model ESP for measuring and investigating particle collection efficiency under various velocities and electric/voltage characteristics. Production of electric field in a lab model ESP of this type and its testing constitutes a novel approach as such work is not found in the public domain. The experimental results show that ESP collection efficiency is higher at high voltages and at low fly ash velocity and the collection efficiency rapidly decreases when voltage reduces. A mathematical model was developed and validated with the experimental measurements to confirm the collection efficiency. By analysing the various conditions and scenarios, an optimum operational condition within an operational range were developed and recommended for future ESP operation. By implementing a TR (Transformer–Rectifier) in different collection chambers, power consumption of the ESP can be reduced. The research also revealed new information on the particulate matter size distribution and the collection of submicron particles from flue gas of coal-fired power plants. Particle size distribution analysis was conducted using a Mastersizer and the morphology of the particles was analysed using a Scanning Electron Microscope (SEM). Size distribution analysis suggested that higher voltage and lower flue gas velocity will be more suitable to capture submicron particles. The morphology study indicated that smaller particles have a tendency to agglomerate with bigger particles. Overall, this thesis provides new knowledge about Electrostatic Precipitator operation with new geometries and under various electric field conditions at a laboratory scale, whilst achieving operational efficiency improvement and improving the capture of sub-micron particulate matter. The knowledge obtained from this research would be a good basis to operate industrial ESPs for future sustainable coal-fired power generation.

Particulate matter emission is one of the major air pollution problems of coal fired power plants. Fine particulates constitute a smaller fraction by weight of the total suspended particle matter in a typical particulate emission, but they are considered potentially hazardous to health because of the high probability of deposition in deeper parts of the respiratory tract. Electrostatic precipitators (ESP) are the most widely used devices that are capable of controlling particulate emission effectively from power plants and other process industries. Although the dust collection efficiency of the industrial precipitator is reported as about 99.5%, an anticipation of future stricter environmental protection agency (EPA) regulations have led the local power station seeking new technologies to achieve the new requirements at minimum cost and thus control their fine particulate emissions to a much greater degree than ever before.

This study aims to identify the options for controlling fine particle emission through improvement of the ESP performance efficiency. An ESP system consists of flow field, electrostatic field and particle dynamics. The performance of an ESP is significantly affected by its complex flow distribution arising as a result of its complex internal geometry, hence the aerodynamic characteristics of the flow inside an ESP always need considerable attention to improve the efficiency of an ESP. Therefore, a laboratory scale ESP model, geometrically similar to an industrial ESP, was designed and fabricated at the Thermodynamics Laboratory of CQUniversity, Australia to examine the flow behaviour inside the ESP. Particle size and shape morphology analyses were conducted to reveal the properties of the fly ash particles which were used for developing numerical models of the ESP.

Numerical simulations were carried out using Computational Fluid Dynamics (CFD) code FLUENT and comparisons were made with the experimental results. The ESP was modelled in two steps. Firstly, a novel 3D fluid (air) flow was modelled considering the detailed geometrical configuration inside the ESP. A novel boundary condition was applied at the inlet boundary of this model to overcome all previous assumptions on uniform velocity at the inlet boundary. Numerically predicted velocity profiles inside the ESP model are compared with the measured data obtained from the laboratory experiment. The model with a novel boundary condition predicted the flow distribution more accurately. In the second step, as the complete ESP system consists of an electric field and a particle phase in addition to the fluid flow field, a two dimensional ESP model was developed. The electrostatic force was applied to the flow equations using User Defined Functions (UDF). A discrete phase model was incorporated with this two dimensional model to study the effect of particle size, electric field and flue gas flow on the collection efficiency of particles inside the ESP. The simulated results revealed that the collection efficiency cannot be improved by the increased electric force only unless the flow velocity is optimized.

The CFD model was successfully applied to a prototype ESP at the power plant and used to recommend options for improving the efficiency of the ESP. The aerodynamic behaviour of the flow was improved by geometrical modifications in the existing 3D numerical model. In particular, the simulation was performed to improve and optimize the flow in order to achieve uniform flow and to increase particle collection inside the ESP. The particles injected in the improved flow condition were collected with higher efficiency after increasing the electrostatic force inside the 2D model. The approach adopted in this study to optimize flow and electrostatic field properties is a novel approach for improving the performance of an electrostatic precipitator.

The thermal efficiency of gas turbines increases with the temperature of the gas entering its turbine component. To enable high inlet temperatures, even those that far exceed the melting point of the turbine materials, the turbine must be cooled. One way is by internal cooling, where cooler air passes through U-ducts embedded inside turbine vanes and blades. Since the flow and heat transfer in these ducts are highly complicated, computational fluid dynamics (CFD) based on RANS have been used extensively to explore and assess design concepts. However, RANS have been found to be unreliable – giving accurate results for some designs but not for others. In this study, large-eddy simulations (LES) were performed for a U-duct with a trapezoidal cross section to assess four widely used RANS turbulence models: realizable k-ε (k-ε), shear-stress transport (SST), Reynolds stress model with linear pressure strain (RSM-LPS), and the seven-equation stress-omega full Reynolds stress model (RSM).

Αντικείμενο της διατριβής είναι η υπολογιστική και πειραματική διερεύνηση των φαι-νομένων μεταφοράς θερμότητας και μάζας σε πρότυπη εργαστηριακή μονάδα μηχανικής ξήρανσης δια συναγωγής, η οποία σχεδιάσθηκε, κατασκευάσθηκε και εξοπλίσθηκε με μετρητικό εξοπλισμό και ολοκληρωμένο σύστημα ελέγχου, συλλογής και επεξεργασίας των πειραματικών μετρήσεων. Στην εργαστηριακή μονάδα ξήρανσης παρέχεται η δυνατότητα μεταβολής και ελέγχου των βασικών παραμέτρων από τις οποίες επηρεάζεται η θερμική διεργασία της ξήρανσης, όπως η ταχύτητα, η θερμοκρασία και η υγρασία του αέρα ξήρανσης. Η μέτρηση της αποβολής της περιεχομένης υγρασίας στο υπό ξήρανση προϊόν πραγματοποιείται μέσω δυναμοκυψελών υψηλής ακρίβειας, ενώ η χωρική κατανομή της ταχύτητας ροής στην είσοδο του θαλάμου ξήρανσης κατά τη διεξαγωγή κάθε κύκλου πειραμάτων συνεχώς μετρείται με συστοιχία σωλήνων pitot και ενός συστήματος συγκροτούμενου από ηλεκτροβαλβίδες και μορφομετατροπέα πίεσης. Η χωρική κατανομή της θερμοκρασίας και της ταχύτητας στον θάλαμο ξήρανσης είναι δυνατή μέσω αισθητηρίων που προσαρμόζονται σε καρτεσιανό σύστημα μετακίνησης ελεγχόμενου από υπολογιστή το οποίο σχεδιάσθηκε, κατασκευάσθηκε και τοποθετήθηκε στην έξοδο του κατακόρυφου θαλάμου ξήρανσης αποτελώντας πλέον αναπόσπαστο τμήμα της εργαστηριακής μονάδας. Όλες οι μετρήσεις πραγματοποιήθηκαν στον κατακόρυφο θάλαμο ξήρανσης, ενώ παρέχεται η δυνατότητα διεξαγωγής μετρήσεων και στο θάλαμο ξήρανσης οριζόντιας διάταξης. Στα πλαίσια της διατριβής, έγινε συστηματική πειραματική διερεύνηση της ξήρανσης δια συναγωγής σε φέτες κυδωνιών και μελετήθηκε η επίδραση διαφόρων παραμέτρων που επηρεάζουν τη θερμική διεργασία της ξήρανσης σε αυτό το αγροτικό προϊόν, για θερμοκρασίες αέρα 40, 50 και 60οC και ταχύτητες αέρα 1, 2 και 3 m/s. Σκοπός των με-τρήσεων ήταν ο προσδιορισμός: (i) της επίδρασης της θερμοκρασίας και της ταχύτητας του αέρα στις καμπύλες ξήρανσης κυλινδρικών φετών κυδωνιού, (ii) της επίδρασης του πάχους των κυλινδρικών φετών του κυδωνιού στις καμπύλες ξήρανσης, (iii) της επίδρασης του προσανατολισμού των κυλινδρικών φετών κυδωνιού, ως προς τη διεύθυνση της προσπίπτουσας ροής, στις καμπύλες ξήρανσης, (iv) της προσαρμογής των καμπύλων ξήρανσης σε διάφορα απλά μοντέλα ξήρανσης λεπτού στρώματος, v) των ενεργών συντελεστών διάχυσης υγρασίας για την κάθε περίπτωση με την μέθοδο της κλίσης (slope method) και οι οποίοι συσχετίστηκαν με τη θερμοκρασία του αέρα ξήρανσης έτσι ώστε ο συντελεστής διάχυσης της υγρασίας να εκφρασθεί με την εξίσωση μορφής τύπου Arrhenius και vi) των διεπιφανειακών συντελεστών μεταφοράς θερμότητας και μάζας οι οποίοι στη εκφράζονται ως συνάρτηση των αδιάστατων αριθμών Nu, Re και Pr με τη μορφή Nu=aRebPr1/3. Η προσομοίωση του ρευστοθερμικού πεδίου στο θάλαμο ξήρανσης και ο υπολογισμός των διεπιφανειακών συντελεστών μεταφοράς θερμότητας και μάζας γύρω από την επι-φάνειας του προϊόντος πραγματοποιήθηκε με χρήση της εργαλείων της υπολογιστικής ρευστοδυναμικής (CFD). Έγιναν προσομοιώσεις CFD μόνιμης κατάστασης (steady-state), θεωρώντας τυρβώδη ροή ενώ ο θάλαμος ξήρανσης και η κυλινδρική φέτα του κυδωνιού εξιδανικεύθηκε ως μια δισδιάστατη αξονοσυμμετρική διαμόρφωση. Ως μοντέλο τύρβης χρησιμοποιήθηκε το μοντέλο SST (Shear Stress Transport) k-ω, ενώ για την προσέγγιση του οριακού στρώματος στα τοιχώματα του προϊόντος επιλέχθηκε το μοντέλο LRNM (Low Reynolds Number Model). Από την επίλυση των πεδίων ροής και θερμοκρασίας προσδιορίσθηκαν οι κατανομές των διεπιφανειακών συντελεστών στην προσήνεμη και στην υπήνεμη πλευρά της κυλινδρικής φέτας του κυδωνιού για όλες τις πειραματικές συνθήκες. Από τον υπολογισμό του μέσων σταθμισμένων τιμών του διεπιφανειακού συντελεστή μεταφοράς θερμότητας συνάγεται μια συσχέτιση των αδιάστατων αριθμών Nu, Re και Pr, στη μορφή Nu=aRebPr1/3, που ως εύρημα εμπλουτίζει την υφιστάμενη βιβλιογραφία. Στο τελικό στάδιο της διατριβής, αναπτύχθηκε και αποτιμήθηκε σε σύγκριση με τις πειραματικές μετρήσεις ένα μονοδιάστατο αριθμητικό μοντέλο μη-μόνιμης μεταφοράς θερμότητας και μάζας για την προσομοίωση των καμπυλών ξήρανσης σε κυλινδρικές φέτες κυδωνιών. Στο μοντέλο, η μεταφορά θερμότητας εντός του προϊόντος γίνεται με αγωγή ενώ η μεταφορά μάζας γίνεται με υγρή διάχυση, με την εξάτμιση του περιεχόμενου νερού στις φέτες του κυδωνιών να λαμβάνει χώρα από την προσήνεμη και την υπήνεμη επιφάνεια. Στο αριθμητικό μοντέλο, λαμβάνεται υπόψη η συρρίκνωση της κυλινδρικής φέτας του κυδωνιού, θεωρώντας ότι ο όγκος της προϊόντος μειώνεται κάθε φορά κατά τον όγκο του νερού που εξατμίζεται τις δύο επιφάνειες της φέτας. Στον αριθμητικό κώδικα, οι θερμοφυσικές ιδιότητες του κυδωνιού και του αέρα προσδιορίζονται από σχέσεις που συναντώνται στη βιβλιογραφία, ο ενεργός συντελεστής διάχυσης της υγρασίας εισάγεται ως αποτέλεσμα της επεξεργασίας των πειραματικών μετρήσεων, ενώ για τους συντελεστές μεταφοράς χρησιμοποιήθηκαν οι μέσες σταθμισμένες τιμές των διεπιφανειακών συντελεστών μεταφοράς θερμότητας και μάζας, ως αποτέλεσμα των CFD προσομοιώσεων και για περίπτωση μη-συζυγούς προσέγγισης (non-conjugated approach). Στοχεύοντας στην καλύτερη προσαρμογή των πειραματικών μετρήσεων και των υπολογιστικών αποτελεσμάτων, χρησιμοποιήθηκε η ανάλυση μη-γραμμικής παλινδρόμησης, με τους αλγόριθμους SQP (Sequential Quadratic Programming) και εσωτερικού σημείου (internal point), για τον προσδιορισμό των συντελεστών της εξίσωσης διάχυσης της υγρασίας, με μικρή όμως βελτίωση των υπολογιστικών αποτελεσμάτων, σε συνδυασμό με τη σημαντική αύξηση του χρόνου υπολογισμού. Συμπερασματικά, από τη συνολική αξιολόγηση των αποτελεσμάτων του αριθμητικού κώδικα αποδείχθηκε ότι το προτεινόμενο αριθμητικό μοντέλο που βασίζεται στη διάχυση είναι ικανό να περιγράψει αποτελεσματικά τη σύζευξη της μεταφοράς της θερμότητας και της μάζας όπως και να αποτυπώσει ικανοποιητικά τη χρονική εξέλιξη της περιεχόμενης υγρασίας και θερμοκρασίας εντός του προϊόντος, με την ελάχιστη χρήση πειραματικών μεταβλητών εισόδου ενώ έχει ελάχιστες υπολογιστικές απαιτήσεις. Για αυτούς τους λόγους μπορεί να θεωρηθεί κατάλληλο για την ανάλυση της διαδικασίας της ξήρανσης δια συναγωγής σε οποιοδήποτε οργανικό ή μη-οργανικό προϊόν.
The objective subject of this thesis is the computational and experimental investigation of heat and mass transfer phenomena in a new laboratory mechanical convection drying unit, which was designed, constructed and equipped with measuring equipment and an integrated control system of collection and processing of experimental measurements. In laboratory drying unit there is an option to change and control the main parameters of which affected the thermal drying process, such as speed, temperature and humidity of the drying air. Measurement of the removal of moisture content in the dried product is carried out through high-precision load cells, and the spatial distribution of the flow velocity at the entrance of the drying chamber during of each experiment, is continuously measured by pitot tube array and a system composed of solenoids and a pressure transducer. The spatial distribution of temperature and velocity in the drying chamber is possible by means of sensors fitted to a computer controlled cartesian motion system which is designed, constructed and placed at the outlet of the vertical drying chamber, constituting an integral part of the facility. All measurements were performed on the vertical drying chamber while it is possible to conduct measurements in a horizontal layout of the drying chamber. In this thesis became systematic experimental investigation of convective drying sliced quince and studied the effect of various parameters affecting the thermal drying process in this agricultural product, for air temperatures of 40, 50 and 60°C and air velocities 1, 2 and 3 m/s. The purpose of the measurements was to determine: (i) the effect of temperature and air velocity in drying curves of cylindrical quince slice, (ii) the effect of the thickness of the cylindrical slice of quince in drying curves, (iii) the effect of the orientation of the cylindrical quince slice, in the direction of incident flow, in the drying curves (iv) the adjusting of the drying curves in several simple thin layer drying models v) the effective moisture diffusivity coefficients for each case with the slope method which correlated with the temperature of the drying air so that the diffusion coefficient of moisture be expressed by Arrhenius type equation form and vi ) the interfacial heat and mass transfer coefficients which expressed as a function of dimensionless numbers Nu, Re and Pr in the form Nu = aRebPr1/3. The simulation of the flow and temperature fields in the drying chamber and the calcu-lation of the interfacial heat and mass transfer coefficients around the surface of the product were performed using the tools of Computational Fluid Dynamics (CFD). CFD simulations were steady state, considering turbulent flow while drying chamber and cy-lindrical slice of quince specialized as an axisymmetric two-dimensional configuration. As turbulence model was used the SST k-ω model while on the approximation of the boundary layer near the walls of the product the LRNM was chosen. By solving the flow and temperature fields determined distributions of interfacial heat and mass transfer coefficients in front and rear of the cylindrical slice of quince for all experimental conditions. The calculation of the weighted average prices of the interfacial heat transfer coefficient indicates a correlation between dimensionless numbers Nu, Re and Pr, in the form Nu = aRebPr1/3, which as finding enriches the existing literature. In the final stage of the thesis, developed and evaluated in comparison with the experi-mental measurements, a one-dimensional transient numerical model of heat and mass transfer to simulate drying curves in cylindrical slices of quince. The heat transfer inside the quince is considered to be by conduction while the moisture transfer is considered to be governed solely by liquid diffusion. Evaporation is considered to take place only from the windward and leeward surface of the quince slice. The numerical model takes into account the shrinkage of the cylindrical slice of quince, assuming that the cylindrical volume decreases each time as much as the volume of water that evaporates on both surfaces of the slice. The numerical code used the thermophysical properties of quince and air from the literature, the effective diffusion coefficient of moisture experimentally determined by the method of the slopes, while the transfer coefficients used the weighted average prices of interfacial heat and mass transfer coefficients derived from the simulations with CFD (non-conjugated approach). In order to achieve higher accuracy between experimental data and predictions, a non-linear regression analysis, using an Arrhenius type effective diffusion equation, was also performed. However, preliminary result, obtained using the SQP (Sequential Quadratic Programming) and Interior Point algorithms for the minimization of the Chi-square function (χ2) showed only small improvement of the calculated results with a significant increase of the computational cost. In conclusion, the overall assessment of the results of the numeric code shown that the proposed numerical model based on diffusion is able to effectively describe the coupling of heat transfer and mass, as to capture the time evolution of moisture content and temperature within the product, with minimum use of experimental input variables and minimum computational requirements. For these reasons it may be considered appropriate to analyze the convective drying process in any organic or non-organic product.

This chapter describes a study designed to evaluate the spectrum of the residence time of the water at different depths of a deep lake, and to examine the mechanisms governing the seasonal cycle of thermal stratification and destratification, with the ultimate aim of assessing the actual exchange time of the lake water. The study was performed on Lake Maggiore (depth 370m) using a multidimensional mathematical model and computer codes for the heat and mass transfer in very large natural water bodies. A 3D Eulerian time-dependent CFD (Computational Fluid Dynamics) code was applied under real conditions, taking into account the effects of the monthly mean values of the mass flow rates and temperatures of all the tributaries, mass flow rate of the Ticino effluent and meteorological, hydrological, and limnological parameters available from the rich data-base of the CNR-ISE (Pallanza). The velocity distributions from these simulations were used to compute the paths of a large number of massless markers with different initial positions and evaluate their residence times in the lake.

The U. S. Department of Energy (DOE) is supporting the development of a next generation nuclear plant (NGNP), which will be based on a very high temperature reactor (VHTR) design. The VHTR is a single-phase helium-cooled reactor wherein the helium will be heated initially to 750 °C and later to temperatures approaching 1000 °C. The high temperatures are desired to increase reactor efficiency and to provide a heat source for the manufacture of hydrogen and other applications. While computational fluid dynamics (CFD) has not been used in the past to design or license nuclear reactors in the U. S., it is expected that CFD will be used in the design and safety analysis of forthcoming designs. This is partly because of the maturity of CFD and partly because detailed information is desired of the flow and heat transfer inside the reactor to avoid hot spots and other conditions that might compromise reactor safety. Numerical computations of turbulent flow should be validated against experimental data for flow conditions that contain some or all of the physics expected in the thermal fluid machinery of interest. To this end, a scaled model of a narrow slice of the lower plenum of the prismatic VHTR was constructed and installed in the Idaho National Laboratory’s (INL) matched index of refraction (MIR) test facility and data were taken. The data were then studied and compared to CFD calculations to help determine their suitability for validation data. One of the main findings was that the inlet data, which were measured and controlled by calibrated mass flow rotameters and were also measured using detailed stereo particle image velocimetry (PIV) showed considerable discrepancies in mass flow rate between the two methods. The other finding was that a randomly unstable recirculation zone occurs in the flow. This instability has a very significant effect on the flow field in the vicinity of the inlet jets. Because its time scale is long and because it is apparently a random instability, it was deemed undesirable for a validation data set. It was predicted using CFD that by eliminating the first of the four jets, the recirculation zone could be stabilized. The present paper reports detailed results for the three-jet case with comparisons to the four-jet data inasmuch as three-jet data are still unavailable. Hence, the present simulations are true or blind predictions.

Development of geometry-independent computational method and educational codes for simulation of 2D flows around objects of complex geometry is presented. Referred as immersed boundary method, it introduces virtual forcing to governing equations to represent the effect of physical boundaries. The present method is based on a finite-volume approach on a staggered grid with a fractional-step method to solve Navier-Stokes equation and continuity equation. Both momentum and mass forcings are introduced on and inside the object to satisfy no-slip condition and mass conservation. Since Cartesian grid lines in general do not coincide with the immersed boundaries, several interpolation schemes are employed. Several examples are simulated using the method presented in this study and the results agree well with other results. Both user-friendly preprocessor with GUI and FORTRAN-based solver are open to the public for educational purposes.

Industrial energy systems involve many components where fluid flow, heat and mass transfer are the important transport processes. For design and development as well as investigation of innovative ideas, computational methods are of vital importance. Since the early pioneering works (by, e.g., the Spalding group) in the 1960s and 1970s, the development has been tremendous and nowadays CFD (computational fluid dynamics) is an established methodology. During the same time period computer capacities in addition have been brought extensively forward. Although still many topics need further development, applications to industrial problems are constantly increasing. In this paper examples of computational heat transfer, fluid flow and related transport phenmena applied to plate heat exchangers (PHEs), gas turbine heat transfer, and fuel cells are highlighted.

Increasingly, different concepts such as safety, hygiene and comfort interact in the characterization of the workplaces. Being comfortable during periods of low activity, seems to be a requirement for most people and, secondly, in other sectors where performance is critical, the priorities are different. In both cases, the intellectual and physical performance is strongly affected by the sensation of thermal comfort. Thus, various approaches can be applied to provide comfort, for example, in the design of buildings and the selection of appropriate clothing. Comfortable clothing is a complex and interdisciplinary concept, consisting of a balance of the sensorial, psychological and physiological aspects [1, 2]. In objective terms, the behavior of the human body depends on several factors: temperature, air velocity and humidity, production of metabolic heat and clothing insulation. All these factors determine the heat and mass transfer processes between the human body and the environment. The Computational Fluid Dynamics (CFD) numerical simulation has been a powerful tool in this investigation field. A combined simulation of a room and a thermal manikin has been developed in the FLUENT code. Using a manikin with real dimensions, divided into parts with different temperatures, seems important to give accurate fluid flow and moisture distributions.

The flow modeling approaches employed in Computational Fluid Dynamics (CFD) codes dedicated to the thermal analysis of electronic equipment are generally not specific for the analysis of forced airflows over populated Printed Circuit Boards. This limitation has been previously highlighted [1], with component junction temperature prediction errors of up to 35% reported. This study evaluates the predictive capability of candidate turbulence models more suited to the analysis of electronic component heat transfer. Significant improvements in component junction temperature prediction accuracy are obtained, relative to a standard high-Reynolds number k-e model, which are attributed to better prediction of both board leading edge heat transfer and component thermal interaction. Such improvements would enable parametric analysis of product thermal performance to be undertaken with greater confidence in the thermal design process, and the generation of more accurate temperature boundary conditions for use in Physics-of-Failure based reliability prediction methods. The case is made for vendors of CFD codes dedicated to the thermal analysis of electronics to consider the adoption of eddy viscosity turbulence models more suited to board-level analysis.

Interfacial dynamics in micron-sized geometries is getting more attention as technologies evolve towards small scales. In pharmaceutical and biochemical applications, fluid droplets moving in microfluidic channels are utilized as “micro-reactors” to precisely control the reaction rate and enhance the mixing efficiency of the reagents. In emulsion generations, low polydispersity of droplets is achieved by the “passive” fission, fusion and sorting of droplets in microfluidic channels. In nanolithography, a micron/nano sized fluid bridge forms between the substrate and the depositing tip. The ability to numerically simulate in a fast and accurate manner the behavior of a small volume of fluid in microfluidic systems is thus demanded. We utilize a 3D Spectral Boundary Element algorithm for interfacial dynamics in microfluidics that exploits all the benefits of the spectral methods (i.e. exponential convergence and numerical stability) but it is not affected by the disadvantage of the spectral methods used in volume discretization to create denser systems. We consider the deformation and migration of droplets in a square microfluidic channel. We investigate the behavior for a single droplet. The effects of droplet size, flow rate and viscosity ratio will be presented.