Dynamic Model of a Small Once-Through Boiler

In this study, a model for the unsteady dynamic behaviour of a once-through counter flow boiler that uses an organic working fluid is presented. The boiler is a compact waste-heat boiler without a furnace and it has a preheater, a vaporiser and a superheater. The relative lengths of the boiler parts vary with the operating conditions since they are all parts of a single tube. The present research is a part of a study on the unsteady dynamics of an organic Rankine cycle power plant and it will be a part of a dynamic process model. The boiler model is presented using a selected example case that uses toluene as the process fluid and flue gas from natural gas combustion as the heat source. The dynamic behaviour of the boiler means transition from the steady initial state towards another steady state that corresponds to the changed process conditions. The solution method chosen was to find such a pressure of the process fluid that the mass of the process fluid in the boiler equals the mass calculated using the mass flows into and out of the boiler during a time step, using the finite difference method. A special method of fast calculation of the thermal properties has been used, because most of the calculation time is spent in calculating the fluid properties. The boiler was divided into elements. The values of the thermodynamic properties and mass flows were calculated in the nodes that connect the elements. Dynamic behaviour was limited to the process fluid and tube wall, and the heat source was regarded as to be steady. The elements that connect the preheater to thevaporiser and the vaporiser to the superheater were treated in a special way that takes into account a flexible change from one part to the other. The model consists of the calculation of the steady state initial distribution of the variables in the nodes, and the calculation of these nodal values in a dynamic state. The initial state of the boiler was received from a steady process model that isnot a part of the boiler model. The known boundary values that may vary during the dynamic calculation were the inlet temperature and mass flow rates of both the heat source and the process fluid. A brief examination of the oscillation around a steady state, the so-called Ledinegg instability, was done. This examination showed that the pressure drop in the boiler is a third degree polynomial of the mass flow rate, and the stability criterion is a second degree polynomial of the enthalpy change in the preheater. The numerical examination showed that oscillations did not exist in the example case. The dynamic boiler model was analysed for linear and step changes of the entering fluid temperatures and flow rates.The problem for verifying the correctness of the achieved results was that there was no possibility o compare them with measurements. This is why the only way was to determine whether the obtained results were intuitively reasonable and the results changed logically when the boundary conditions were changed. The numerical stability was checked in a test run in which there was no change in input values. The differences compared with the initial values were so small that the effects of numerical oscillations were negligible. The heat source side tests showed that the model gives results that are logical in the directions of the changes, and the order of magnitude of the timescale of changes is also as expected. The results of the tests on the process fluid side showed that the model gives reasonable results both on the temperature changes that cause small alterations in the process state and on mass flow rate changes causing very great alterations. The test runs showed that the dynamic model has no problems in calculating cases in which temperature of the entering heat source suddenly goes below that of the tube wall or the process fluid.

Theoretical and numerical study of thermomagnetic convection in magnetic fluids

In this thesis, the magnetic field control of convection instabilities and heat and mass transfer processesin magnetic fluids have been investigated by numerical simulations and theoretical considerations. Simulation models based on finite element and finite volume methods have been developed. In addition to standard conservation equations, themagnetic field inside the simulation domain is calculated from Maxwell equations and the necessary terms to take into account for the magnetic body force and magnetic dissipation have been added to the equations governing the fluid motion.Numerical simulations of magnetic fluid convection near the threshold supportedexperimental observations qualitatively. Near the onset of convection the competitive action of thermal and concentration density gradients leads to mostly spatiotemporally chaotic convection with oscillatory and travelling wave regimes, previously observed in binary mixtures and nematic liquid crystals. In many applications of magnetic fluids, the heat and mass transfer processes including the effects of external magnetic fields are of great importance. In addition to magnetic fluids, the concepts and the simulation models used in this study may be applied also to the studies of convective instabilities in ordinary fluids as well as in other binary mixtures and complex fluids.

O experimento de Clement-Desormes no século XXI

The Clement-Desormes experiment is reviewed. By reason of a finite difference between the pressure within the system and its surroundings, Bertrand and McDonald have criticized the usual consideration of the adiabatic expansion as reversible. Garland, Nibler and Shoemaker oppose, defining regions through virtual boundaries where the surroundings do not operate. For Holden, the use of virtual boundaries is expendable. Experiments cannot support a hypothesis testing due to experiment's intrinsic uncertainty. The role of polytropy in uncertainty is discussed. Both thermodynamic definitions and kinetic model depict the real processes as irreversible phenomena and the reversible ones as a limiting hypothetical case.

Kestomagneettigeneraattorin lämmönsiirron mallinnus tuulivoimakäytössä

Pyörivien sähkökoneiden suunnittelussa terminen suunnittelu on yhtä tärkeää kuin sähköinen ja mekaaninen suunnittelukin. Tässä diplomityössä tarkoituksena on kehittää ilmajäähdytteisten kestomagneettigeneraattorien laskentaan soveltuva lämmönsiirtymismalli, jolla staattorin lämpötilajakauma voitaisiin selvittää. Kehitetty lämmönsiirtymismalli perustuu kolmiulotteiseen äärellisen erotuksen (finite difference) menetelmään. Malli ottaa huomioon lämmönjohtumisen staattorin aktiiviosissa ja konvektion jäähdytysilmavirtaan. Mallissa on myös yksinkertainen painehäviölaskenta jäähdytysjärjestelmän komponenttien mitoittamista varten. Laskentamallilla lasketaan esimerkkitapauksena 4,3 MW:n kestomagneettigeneraattorin jäähdytystä eri toimintapisteissä. Tuloksia verrataan CFD-mallinnuksen antamiin tuloksiin sekä kokeellisten mittausten antamiin tuloksiin.

Insulation System in an Integrated Motor Compressor

A high-speed and high-voltage solid-rotor induction machine provides beneficial features for natural gas compressor technology. The mechanical robustness of the machine enables its use in an integrated motor-compressor. The technology uses a centrifugal compressor, which is mounted on the same shaft with the high-speed electrical machine driving it. No gearbox is needed as the speed is determined by the frequency converter. The cooling is provided by the process gas, which flows through the motor and is capable of transferring the heat away from the motor. The technology has been used in the compressors in the natural gas supply chain in the central Europe. New areas of application include natural gas compressors working at the wellheads of the subsea gas reservoir. A key challenge for the design of such a motor is the resistance of the stator insulation to the raw natural gas from the well. The gas contains water and heavy hydrocarbon compounds and it is far harsher than the sales gas in the natural gas supply network. The objective of this doctoral thesis is to discuss the resistance of the insulation to the raw natural gas and the phenomena degrading the insulation. The presence of partial discharges is analyzed in this doctoral dissertation. The breakdown voltage of the gas is measured as a function of pressure and gap distance. The partial discharge activity is measured on small samples representing the windings of the machine. The electrical field behavior is also modeled by finite element methods. Based on the measurements it has been concluded that the discharges are expected to disappear at gas pressures above 4 – 5 bar. The disappearance of discharges is caused by the breakdown strength of the gas, which increases as the pressure increases. Based on the finite element analysis, the physical length of a discharge seen in the PD measurements at atmospheric pressure was approximated to be 40 – 120 m. The chemical aging of the insulation when exposed to raw natural gas is discussed based on a vast set of experimental tests with the gas mixture representing the real gas mixture at the wellhead. The mixture was created by mixing dry hydrocarbon gas, heavy hydrocarbon compounds, monoethylene glycol, and water. The mixture was chosen to be more aggressive by increasing the amount of liquid substances. Furthermore, the temperature and pressure were increased, which resulted in accelerated test conditions. The time required to detect severe degradation was thus decreased. The test program included a comparison of materials, an analysis of the e ects of di erent compounds in the gas mixture, namely water and heavy hydrocarbons, on the aging, an analysis of the e ects of temperature and exposure duration, and also an analysis on the e ect of sudden pressure changes on the degradation of the insulating materials. It was found in the tests that an insulation consisting of mica, glass, and epoxy resin can tolerate the raw natural gas, but it experiences some degradation. The key material in the composite insulation is the resin, which largely defines the performance of the insulation system. The degradation of the insulation is mostly determined by the amount of gas mixture di used into it. The di usion was seen to follow Fick’s second law, but the coe cients were not accurately defined. The di usion was not sensitive to temperature, but it was dependent upon the thermodynamic state of the gas mixture, in other words, the amounts of liquid components in the gas. The weight increase observed was mostly related to heavy hydrocarbon compounds, which act as plasticizers in the epoxy resin. The di usion of these compounds is determined by the crosslink density of the resin. Water causes slight changes in the chemical structure, but these changes do not significantly contribute to the aging phenomena. Sudden changes in pressure can lead to severe damages in the insulation, because the motion of the di used gas is able to create internal cracks in the insulation. Therefore, the di usion only reduces the mechanical strength of the insulation, but the ultimate breakdown can potentially be caused by a sudden drop in the pressure of the process gas.

On the development of an unstructured grid solver for inert and reactive high speed flow simulations

An unstructured grid Euler solver for reactive compressible flow applications is presented. The method is implemented in a cell centered, finite volume context for unstructured triangular grids. Three different schemes for spatial discretization are implemented and analyzed. Time march is implemented in a time-split fashion with independent integrators for the flow and chemistry equations. The capability implemented is tested for inert flows in a hypersonic inlet and for inert and reactive supersonic flows over a 2-D wedge. The results of the different schemes are compared with each other and with independent calculations using a structured grid code. The strengths and the possible weaknesses of the proposed methods are discussed.

Mechanical design and analysis of modular active magnetic bearing test rig

Tässä työssä on tutkittu modulaarisen aktiivimagneettilaakeroidun koelaitteen mekaanista suunnittelua ja analysointia. Suurnopeusroottorin suunnittelun teoria on esitelty. Lisäksi monia analyyttisiä mallinnusmenetelmiä mekaanisten kuormitusten mallintamiseksi on esitelty. Koska kyseessä on suurnopeussähkökone, roottoridynamiikka ja sen soveltuvuus suunnittelussa on esitelty. Magneettilaakerien rakenteeseen ja toimintaan on tutustuttu osana tätä työtä. Kirjallisuuskatsaus nykyisistä koelaitteista esimerkiksi komponenttien ominaisuuksien tunnistamiseen ja roottoridynamiikan tutkimuksiin on esitelty. Työn rajauksena on konseptisuunnittelu muunneltavalle magneettilaakeroidulle (AMB) koelaitteelle ja suunnitteluprosessin dokumentointi. Muunneltavuuteen päädyttiin, koska se mahdollistaa erilaisten komponenttiasetteluiden testaamisen erilaisille magneettilaakerikokoonpanoille ja roottoreille. Pääpaino tässä työssä on suurnopeus induktiokoneen roottorin suunnittelussa ja mallintamisessa. Modulaaristen toimilaitteiden kuten magneettilaakerien ja induktiosähkömoottorin rakenne on esitelty ja modulaarisen rakenteen käytettävyyden hyödyistä koelaitekäytössä on dokumentoitu. Analyyttisiä ja elementtimenetelmään perustuvia tutkimusmenetelmiä on käytetty tutkittaessa suunniteltua suurnopeusroottoria. Suunnittelun ja analysoinnin tulokset on esitelty ja verrattu keskenään eri mallinnusmenetelmien välillä. Lisäksi johtopäätökset sähkömagneettisten osien liittämisen monimutkaisuudesta ja vaatimuksista roottoriin ja toimilaitteisiin sekä mekaanisten että sähkömagneettisten ominaisuuksien optimoimiseksi on dokumentoitu.

THE EFFECTIVE DIFFUSIVITY OF CLOVE(Eugenia caryophyllus) ESSENTIAL OIL IN PRESSURIZED CO2

The effective diffusivity of clove essential oil in subcritical liquid CO2 was estimated. The experimental apparatus employed was a fixed-bed extractor. The fixed bed was formed with grounded (mesh -32 + 65) and compacted clove buds which were considered a solid element. The effective diffusion coefficient was evaluated by fitting the experimental concentration profile to the unsteady state mass balance equation for unidirectional diffusion in a finite solid medium. The diffusion coefficient was related to the concentration of oil in the solid by an exponential function. The estimated values of the effective diffusion coefficient varied from 3.64 to 5.22x10-10 m2/s. The average relative errors were lower than 3.1%.

Etudes ab initio des effets de la température sur le spectre optique des semi-conducteurs

Thèse réalisée en cotutelle avec l'Université Catholique de Louvain (Belgique)

Modélisation mathématique de la propagation de la malaria.

Un modèle mathématique de la propagation de la malaria en temps discret est élaboré en vue de déterminer l'influence qu'un déplacement des populations des zones rurales vers les zones urbaines aurait sur la persistance ou la diminution de l'incidence de la malaria. Ce modèle, sous la forme d'un système de quatorze équations aux différences finies, est ensuite comparé à un modèle analogue mais en temps continu, qui prend la forme d'équations différentielles ordinaires. Une étude comparative avec la littérature récente permet de déterminer les forces et les faiblesses de notre modèle.

Invariant discretizations of partial differential equations

Un algorithme permettant de discrétiser les équations aux dérivées partielles (EDP) tout en préservant leurs symétries de Lie est élaboré. Ceci est rendu possible grâce à l'utilisation de dérivées partielles discrètes se transformant comme les dérivées partielles continues sous l'action de groupes de Lie locaux. Dans les applications, beaucoup d'EDP sont invariantes sous l'action de transformations ponctuelles de Lie de dimension infinie qui font partie de ce que l'on désigne comme des pseudo-groupes de Lie. Afin d'étendre la méthode de discrétisation préservant les symétries à ces équations, une discrétisation des pseudo-groupes est proposée. Cette discrétisation a pour effet de transformer les symétries ponctuelles en symétries généralisées dans l'espace discret. Des schémas invariants sont ensuite créés pour un certain nombre d'EDP. Dans tous les cas, des tests numériques montrent que les schémas invariants approximent mieux leur équivalent continu que les différences finies standard.

FDTD analysis of a symmetric T-strip fed wideband rectangular microstrip antenna

A theoretical analysis of a symmetric T-shaped rnicrostripfed rectangular microstrip antenna using the finite-difference titnedoniain (FDTD) method is presented in this paper. The resonant frequency, return loss, impedance bandwidth, and radiation patterns are predicted and are in good agreement with the measured results

INVESTIGATIONS ON SLOT-LOADED SQUARE MICROSTRIP PATCH ANTENNAS FOR COMPACT DUAL FREQUENCY DUAL POLARIZED OPERATIONS

The thesis is the outcome of the experimental and theoretical investigations carried out on a novel slotted microstrip antenna.The antenna excites two resonance frequencies and provides orthogonal polarization. The radiation characteristics of the antenna are studied in detail. The antenna design is optimized using IE3D electromagnetic simulation tool. The frequency-Difference Time-Domain (FDTD) method is employed for the analysis of the antenna.The antenna can be used for personal and satellite communication applications.

DEVELOPMENT AND ANLAYSIS OF BROADBAND L·STRIP FED MICROSTRIP ANTENNAS

The thesis is the outcome of the experimental and theoretical Investigations on novel feeding techniques for bandwidth enhancement of microstrip patches. The new feeding techniques provide bandwidth enhancement without deteriorating the radiation characteristics of the antenna. The antenna is analysed using finite Difference Time Domain (FDTD) method. The predicated results are compared with the experimental results and excellent agreement is observed. The results are also verified using IE3D simulation software. The antenna is suitable for personal and broadband communications.

Development and Analysis of Microstrip Antennas For Dual Band Microwave Communication

The thesis explores the outcome of the exhaustive theoretical and experimental investigations performed on Octagonal Microstrip Antenna configurations. Development of the MATLAB TM backed 3D-Conformal Finite Difference Time Domain (CFDTD)Modeller for the numerical computation of the radiation characteristics of the antenna is the theme of the work. The predicted results are verified experimentally and by IE3D TM simulation. The influence of the patch dimensions,feed configurations,feed dimensions and feed positions upon the radiation performance of the antenna is studied in detail. Octagonal Microstrip Antenna configurations suitable for Mobile-Bluetooth application is dealt in detail. A simple design formula for the regular Octagonal geometry is also presented. A compact planar multi band antenna for GPS/DCS/2.4/5.8GHz WLAN application is included as appendix A. Planar near field measurement technique is explained in appendix B.