Paperin pintaenergian vaikutus painovärin tarpeeseen

Työn tavoitteena oli selvittää pintaenergian vaikutus painovärin tarpeeseen IGT-laboratoriopainokoneella. Lisäksi selvitettiin onko karheudella vaikutusta kontaktikulmamittaukseen. Työ koostui kirjallisuusosasta ja kokeellisesta osasta. Kirjallisuusosassa tarkasteltiin aihepiiriä koskevia teorioita ja aiempia tutkimuksia. Kokeellisessa osassa valmistettiin laboratoriomittakaavassa pintaenergialtaan ja karheudeltaan erilaisia papereita. Pintaenergiaa muuteltiin kolmen erilaisen massaliiman ja kahden erilaisen pintaliiman avulla. Paperin karheutta säädettiin kalanteroinnilla. Papereista mitattiin pintaominaisuuksia, nesteen ja paperin välinen kontaktikulma, sekä suoritettiin painatuskokeita. Pintaenergia laskettiin kirjallisuudesta löytyneiden teorioiden perusteella ja näiden pintaenergian laskemismenetelmien käyttökelpoisuutta pintaenergian määrittämiseen selvitettiin. Lisäksi tutkittiin pintaenergian vaikutusta painovärin tarpeeseen. Tämä työ selvästi osoittaa, että yhden, kahden ja kolmen komponentin menetelmällä ja nesteitä vaihtelemalla saadaan erilaisia pintaenergian tuloksia. Pintaenergian määrittämiseen vaikuttaa myös kontaktikulmamittausten hankaluus ja kontaktikulman suuri vaihtelu, joka vaikuttaa koko mittaussarjan luotettavuustarkasteluun. Lisäksi paperin karheudella ja huokoisuudella oli merkitystä kontaktikulmamittaukseen.

Modelling variation of needle density of Scots pine at high latitudes

[Abstract]

Effect of compression wood on surface roughness and surface absorption of medium density fiberboard

Abstract

Torque vibration model of axial-flux surface-mounted permanent magnet synchronous machine

In order that the radius and thus ununiform structure of the teeth and otherelectrical and magnetic parts of the machine may be taken into consideration the calculation of an axial flux permanent magnet machine is, conventionally, doneby means of 3D FEM-methods. This calculation procedure, however, requires a lotof time and computer recourses. This study proves that also analytical methods can be applied to perform the calculation successfully. The procedure of the analytical calculation can be summarized into following steps: first the magnet is divided into slices, which makes the calculation for each section individually, and then the parts are submitted to calculation of the final results. It is obvious that using this method can save a lot of designing and calculating time. Thecalculation program is designed to model the magnetic and electrical circuits of surface mounted axial flux permanent magnet synchronous machines in such a way, that it takes into account possible magnetic saturation of the iron parts. Theresult of the calculation is the torque of the motor including the vibrations. The motor geometry and the materials and either the torque or pole angle are defined and the motor can be fed with an arbitrary shape and amplitude of three-phase currents. There are no limits for the size and number of the pole pairs nor for many other factors. The calculation steps and the number of different sections of the magnet are selectable, but the calculation time is strongly depending on this. The results are compared to the measurements of real prototypes. The permanent magnet creates part of the flux in the magnetic circuit. The form and amplitude of the flux density in the air-gap depends on the geometry and material of the magnetic circuit, on the length of the air-gap and remanence flux density of the magnet. Slotting is taken into account by using the Carter factor in the slot opening area. The calculation is simple and fast if the shape of the magnetis a square and has no skew in relation to the stator slots. With a more complicated magnet shape the calculation has to be done in several sections. It is clear that according to the increasing number of sections also the result will become more accurate. In a radial flux motor all sections of the magnets create force with a same radius. In the case of an axial flux motor, each radial section creates force with a different radius and the torque is the sum of these. The magnetic circuit of the motor, consisting of the stator iron, rotor iron, air-gap, magnet and the slot, is modelled with a reluctance net, which considers the saturation of the iron. This means, that several iterations, in which the permeability is updated, has to be done in order to get final results. The motor torque is calculated using the instantaneous linkage flux and stator currents. Flux linkage is called the part of the flux that is created by the permanent magnets and the stator currents passing through the coils in stator teeth. The angle between this flux and the phase currents define the torque created by the magnetic circuit. Due to the winding structure of the stator and in order to limit the leakage flux the slot openings of the stator are normally not made of ferromagnetic material even though, in some cases, semimagnetic slot wedges are used. In the slot opening faces the flux enters the iron almost normally (tangentially with respect to the rotor flux) creating tangential forces in the rotor. This phenomenon iscalled cogging. The flux in the slot opening area on the different sides of theopening and in the different slot openings is not equal and so these forces do not compensate each other. In the calculation it is assumed that the flux entering the left side of the opening is the component left from the geometrical centre of the slot. This torque component together with the torque component calculated using the Lorenz force make the total torque of the motor. It is easy to assume that when all the magnet edges, where the derivative component of the magnet flux density is at its highest, enter the slot openings at the same time, this will have as a result a considerable cogging torque. To reduce the cogging torquethe magnet edges can be shaped so that they are not parallel to the stator slots, which is the common way to solve the problem. In doing so, the edge may be spread along the whole slot pitch and thus also the high derivative component willbe spread to occur equally along the rotation. Besides forming the magnets theymay also be placed somewhat asymmetric on the rotor surface. The asymmetric distribution can be made in many different ways. All the magnets may have a different deflection of the symmetrical centre point or they can be for example shiftedin pairs. There are some factors that limit the deflection. The first is that the magnets cannot overlap. The magnet shape and the relative width compared to the pole define the deflection in this case. The other factor is that a shifting of the poles limits the maximum torque of the motor. If the edges of adjacent magnets are very close to each other the leakage flux from one pole to the other increases reducing thus the air-gap magnetization. The asymmetric model needs some assumptions and simplifications in order to limit the size of the model and calculation time. The reluctance net is made for symmetric distribution. If the magnets are distributed asymmetrically the flux in the different pole pairs will not be exactly the same. Therefore, the assumption that the flux flows from the edges of the model to the next pole pairs, in the calculation model from one edgeto the other, is not correct. If it were wished for that this fact should be considered in multi-pole pair machines, this would mean that all the poles, in other words the whole machine, should be modelled in reluctance net. The error resulting from this wrong assumption is, nevertheless, irrelevant.

Reconstruction and analysis of surface variation using photometric stereo

In many industrial applications, accurate and fast surface reconstruction is essential for quality control. Variation in surface finishing parameters, such as surface roughness, can reflect defects in a manufacturing process, non-optimal product operational efficiency, and reduced life expectancy of the product. This thesis considers reconstruction and analysis of high-frequency variation, that is roughness, on planar surfaces. Standard roughness measures in industry are calculated from surface topography. A fast and non-contact method to obtain surface topography is to apply photometric stereo in the estimation of surface gradients and to reconstruct the surface by integrating the gradient fields. Alternatively, visual methods, such as statistical measures, fractal dimension and distance transforms, can be used to characterize surface roughness directly from gray-scale images. In this thesis, the accuracy of distance transforms, statistical measures, and fractal dimension are evaluated in the estimation of surface roughness from gray-scale images and topographies. The results are contrasted to standard industry roughness measures. In distance transforms, the key idea is that distance values calculated along a highly varying surface are greater than distances calculated along a smoother surface. Statistical measures and fractal dimension are common surface roughness measures. In the experiments, skewness and variance of brightness distribution, fractal dimension, and distance transforms exhibited strong linear correlations to standard industry roughness measures. One of the key strengths of photometric stereo method is the acquisition of higher frequency variation of surfaces. In this thesis, the reconstruction of planar high-frequency varying surfaces is studied in the presence of imaging noise and blur. Two Wiener filterbased methods are proposed of which one is optimal in the sense of surface power spectral density given the spectral properties of the imaging noise and blur. Experiments show that the proposed methods preserve the inherent high-frequency variation in the reconstructed surfaces, whereas traditional reconstruction methods typically handle incorrect measurements by smoothing, which dampens the high-frequency variation.

Influence of plasma modification on surface properties and offset printability of coated paper

The properties of the paper surface play a crucial role in ensuring suitable quality and runnability in various converting and finishing operations, such as printing. Plasma surface modification makes it possible to modify the surface chemistry of paper without altering the bulk material properties. This also makes it possible to investigate the role of the surface chemistry alone on printability without influencing the porous structure of the pigment-coated paper. Since the porous structure of a pigment coating controls both ink setting and optical properties, surface chemical changes created by a plasma modification have a potential to decouple these two effects and to permit a better optimization of them both. The aim of this work was to understand the effects of plasma surface modification on paper properties, and how it influences printability in the sheet-fed offset process. The objective was to broaden the fundamental understanding of the role of surface chemistry on offset printing. The effects of changing the hydrophilicity/ hydrophobicity and the surface chemical composition by plasma activation and plasma coatings on the properties of coated paper and on ink-paper interactions as well as on sheet-fed offset print quality were investigated. In addition, the durability of the plasma surface modification was studied. Nowadays, a typical sheet-fed offset press also contains units for surface finishing, for example UVvarnishing. The role of the surface chemistry on the UV-varnish absorption into highly permeable and porous pigment-coated paper was also investigated. With plasma activation it was possible to increase the surface energy and hydrophilicity of paper. Both polar and dispersion interactions were found to increase, although the change was greater in the polar interactions due to induced oxygen molecular groups. The results indicated that plasma activation takes place particularly in high molecular weight components such as the dispersion chemicals used to stabilize the pigment and latex particles. Surface composition, such as pigment and binder type, was found to influence the response to the plasma activation. The general trend was that pilot-scale treatment modified the surface chemistry without altering the physical coating structure, whereas excessive laboratory-scale treatment increased the surface roughness and reduced the surface strength, which led to micro-picking in printing. It was shown that pilot-scale plasma activation in combination with appropriate ink oils makes it possible to adjust the ink-setting rate. The ink-setting rate decreased with linseed-oil-based inks, probably due to increased acid-base interactions between the polar groups in the oil and the plasma-treated paper surface. With mineral-oil-based inks, the ink setting accelerated due to plasma activation. Hydrophobic plasma coatings were able to reduce or even prevent the absorption of dampening water into pigmentcoated paper, even when the dampening water was applied under the influence of nip pressure. A uniform hydrophobic plasma coating with sufficient chemical affinity with ink gave an improved print quality in terms of higher print density and lower print mottle. It was also shown that a fluorocarbon plasma coating reduced the free wetting of the UV-varnish into the highly permeable and porous pigment coating. However, when the UV-varnish was applied under the influence of nip pressure, which leads to forced wetting, the role of the surface chemical composition seems to be much less. A decay in surface energy and wettability occurred during the first weeks of storage after plasma activation, after which it leveled off. However, the oxygen/carbon elemental ratio did not decrease as a function of time, indicating that ageing could be caused by a re-orientation of polar groups or by a contamination of the surface. The plasma coatings appeared to be more stable when the hydrophobicity was higher, probably due to fewer interactions with oxygen and water vapor in the air.

The role of van der Waals interactions in the case of H2 adsorption on ruthenium (0001) surface

Computational material science with the Density Functional Theory (DFT) has recently gained a method for describing, for the first time the non local bonding i.e., van der Waals (vdW) bonding. The newly proposed van der Waals-Density Functional (vdW-DF) is employed here to address the role of non local interactions in the case of H2 adsorption on Ru(0001) surface. The later vdW-DF2 implementation with the DFT code VASP (Vienna Ab-initio Simulation Package) is used in this study. The motivation for studying H2 adsorption on ruthenium surface arose from the interest to hydrogenation processes. Potential energy surface (PES) plots are created for adsorption sites top, bridge, fcc and hcp, employing the vdW-DF2 functional. The vdW-DF yields 0.1 eV - 0.2 eV higher barriers for the dissociation of the H2 molecule; the vdW-DF seems to bind the H2 molecule more tightly together. Furthermore, at the top site, which is found to be the most reactive, the vdW functional suggests no entrance barrier or in any case smaller than 0.05 eV, whereas the corresponding calculation without the vdW-DF does. Ruthenium and H2 are found to have the opposite behaviors with the vdW-DF; Ru lattice constants are overestimated while H2 bond length is shorter. Also evaluation of the CPU time demand of the vdW-DF2 is done from the PES data. From top to fcc sites the vdW-DF computational time demand is larger by 4.77 % to 20.09 %, while at the hcp site it is slightly smaller. Also the behavior of a few exchange correlation functionals is investigated along addressing the role of vdW-DF. Behavior of the different functionals is not consistent between the Ru lattice constants and H2 bond lengths. It is thus difficult to determine the quality of a particular exchange correlation functional by comparing equilibrium separations of the different elements. By comparing PESs it would be computationally highly consuming.

Surface transformation hardening of carbon steel with high power fiber laser

This study investigated the surface hardening of steels via experimental tests using a multi-kilowatt fiber laser as the laser source. The influence of laser power and laser power density on the hardening effect was investigated. The microhardness analysis of various laser hardened steels was done. A thermodynamic model was developed to evaluate the thermal process of the surface treatment of a wide thin steel plate with a Gaussian laser beam. The effect of laser linear oscillation hardening (LLOS) of steel was examined. An as-rolled ferritic-pearlitic steel and a tempered martensitic steel with 0.37 wt% C content were hardened under various laser power levels and laser power densities. The optimum power density that produced the maximum hardness was found to be dependent on the laser power. The effect of laser power density on the produced hardness was revealed. The surface hardness, hardened depth and required laser power density were compared between the samples. Fiber laser was briefly compared with high power diode laser in hardening medium-carbon steel. Microhardness (HV0.01) test was done on seven different laser hardened steels, including rolled steel, quenched and tempered steel, soft annealed alloyed steel and conventionally through-hardened steel consisting of different carbon and alloy contents. The surface hardness and hardened depth were compared among the samples. The effect of grain size on surface hardness of ferritic-pearlitic steel and pearlitic-cementite steel was evaluated. In-grain indentation was done to measure the hardness of pearlitic and cementite structures. The macrohardness of the base material was found to be related to the microhardness of the softer phase structure. The measured microhardness values were compared with the conventional macrohardness (HV5) results. A thermodynamic model was developed to calculate the temperature cycle, Ac1 and Ac3 boundaries, homogenization time and cooling rate. The equations were numerically solved with an error of less than 10-8. The temperature distributions for various thicknesses were compared under different laser traverse speed. The lag of the was verified by experiments done on six different steels. The calculated thermal cycle and hardened depth were compared with measured data. Correction coefficients were applied to the model for AISI 4340 steel. AISI 4340 steel was hardened by laser linear oscillation hardening (LLOS). Equations were derived to calculate the overlapped width of adjacent tracks and the number of overlapped scans in the center of the scanned track. The effect of oscillation frequency on the hardened depth was investigated by microscopic evaluation and hardness measurement. The homogeneity of hardness and hardened depth with different processing parameters were investigated. The hardness profiles were compared with the results obtained with conventional single-track hardening. LLOS was proved to be well suitable for surface hardening in a relatively large rectangular area with considerable depth of hardening. Compared with conventional single-track scanning, LLOS produced notably smaller hardened depths while at 40 and 100 Hz LLOS resulted in higher hardness within a depth of about 0.6 mm.

Heat transfer inside internal combustion engine: modelling and comparison with experimental data

The paper is devoted to study specific aspects of heat transfer in the combustion chamber of compression ignited reciprocating internal combustion engines and possibility to directly measure the heat flux by means of Gradient Heat Flux Sensors (GHFS). A one – dimensional single zone model proposed by Kyung Tae Yun et al. and implemented with the aid of Matlab, was used to obtain approximate picture of heat flux behavior in the combustion chamber with relation to the crank angle. The model’s numerical output was compared to the experimental results. The experiment was accomplished by A. Mityakov at four stroke diesel engine Indenor XL4D. Local heat fluxes on the surface of cylinder head were measured with fast – response, high – sensitive GHFS. The comparison of numerical data with experimental results has revealed a small deviation in obtained heat flux values throughout the cycle and different behavior of heat flux curve after Top Dead Center.

Application and development of numerical methods for the modelling of innovative gas cooled fission reactors

Innovative gas cooled reactors, such as the pebble bed reactor (PBR) and the gas cooled fast reactor (GFR) offer higher efficiency and new application areas for nuclear energy. Numerical methods were applied and developed to analyse the specific features of these reactor types with fully three dimensional calculation models. In the first part of this thesis, discrete element method (DEM) was used for a physically realistic modelling of the packing of fuel pebbles in PBR geometries and methods were developed for utilising the DEM results in subsequent reactor physics and thermal-hydraulics calculations. In the second part, the flow and heat transfer for a single gas cooled fuel rod of a GFR were investigated with computational fluid dynamics (CFD) methods. An in-house DEM implementation was validated and used for packing simulations, in which the effect of several parameters on the resulting average packing density was investigated. The restitution coefficient was found out to have the most significant effect. The results can be utilised in further work to obtain a pebble bed with a specific packing density. The packing structures of selected pebble beds were also analysed in detail and local variations in the packing density were observed, which should be taken into account especially in the reactor core thermal-hydraulic analyses. Two open source DEM codes were used to produce stochastic pebble bed configurations to add realism and improve the accuracy of criticality calculations performed with the Monte Carlo reactor physics code Serpent. Russian ASTRA criticality experiments were calculated. Pebble beds corresponding to the experimental specifications within measurement uncertainties were produced in DEM simulations and successfully exported into the subsequent reactor physics analysis. With the developed approach, two typical issues in Monte Carlo reactor physics calculations of pebble bed geometries were avoided. A novel method was developed and implemented as a MATLAB code to calculate porosities in the cells of a CFD calculation mesh constructed over a pebble bed obtained from DEM simulations. The code was further developed to distribute power and temperature data accurately between discrete based reactor physics and continuum based thermal-hydraulics models to enable coupled reactor core calculations. The developed method was also found useful for analysing sphere packings in general. CFD calculations were performed to investigate the pressure losses and heat transfer in three dimensional air cooled smooth and rib roughened rod geometries, housed inside a hexagonal flow channel representing a sub-channel of a single fuel rod of a GFR. The CFD geometry represented the test section of the L-STAR experimental facility at Karlsruhe Institute of Technology and the calculation results were compared to the corresponding experimental results. Knowledge was gained of the adequacy of various turbulence models and of the modelling requirements and issues related to the specific application. The obtained pressure loss results were in a relatively good agreement with the experimental data. Heat transfer in the smooth rod geometry was somewhat under predicted, which can partly be explained by unaccounted heat losses and uncertainties. In the rib roughened geometry heat transfer was severely under predicted by the used realisable k − epsilon turbulence model. An additional calculation with a v2 − f turbulence model showed significant improvement in the heat transfer results, which is most likely due to the better performance of the model in separated flow problems. Further investigations are suggested before using CFD to make conclusions of the heat transfer performance of rib roughened GFR fuel rod geometries. It is suggested that the viewpoints of numerical modelling are included in the planning of experiments to ease the challenging model construction and simulations and to avoid introducing additional sources of uncertainties. To facilitate the use of advanced calculation approaches, multi-physical aspects in experiments should also be considered and documented in a reasonable detail.

Analytic evaluation of three-phase short circuit demagnetization and hysteresis loss risk in rotor-surface-magnet Permanent-Magnet Synchronous Machines

Permanent magnet synchronous machines (PMSM) have become widely used in applications because of high efficiency compared to synchronous machines with exciting winding or to induction motors. This feature of PMSM is achieved through the using the permanent magnets (PM) as the main excitation source. The magnetic properties of the PM have significant influence on all the PMSM characteristics. Recent observations of the PM material properties when used in rotating machines revealed that in all PMSMs the magnets do not necessarily operate in the second quadrant of the demagnetization curve which makes the magnets prone to hysteresis losses. Moreover, still no good analytical approach has not been derived for the magnetic flux density distribution along the PM during the different short circuits faults. The main task of this thesis is to derive simple analytical tool which can predict magnetic flux density distribution along the rotor-surface mounted PM in two cases: during normal operating mode and in the worst moment of time from the PM’s point of view of the three phase symmetrical short circuit. The surface mounted PMSMs were selected because of their prevalence and relatively simple construction. The proposed model is based on the combination of two theories: the theory of the magnetic circuit and space vector theory. The comparison of the results in case of the normal operating mode obtained from finite element software with the results calculated with the proposed model shows good accuracy of model in the parts of the PM which are most of all prone to hysteresis losses. The comparison of the results for three phase symmetrical short circuit revealed significant inaccuracy of the proposed model compared with results from finite element software. The analysis of the inaccuracy reasons was provided. The impact on the model of the Carter factor theory and assumption that air have permeability of the PM were analyzed. The propositions for the further model development are presented.

Electromagnetic and thermal design of a multilevel converter with high power density and reliability

Electric energy demand has been growing constantly as the global population increases. To avoid electric energy shortage, renewable energy sources and energy conservation are emphasized all over the world. The role of power electronics in energy saving and development of renewable energy systems is significant. Power electronics is applied in wind, solar, fuel cell, and micro turbine energy systems for the energy conversion and control. The use of power electronics introduces an energy saving potential in such applications as motors, lighting, home appliances, and consumer electronics. Despite the advantages of power converters, their penetration into the market requires that they have a set of characteristics such as high reliability and power density, cost effectiveness, and low weight, which are dictated by the emerging applications. In association with the increasing requirements, the design of the power converter is becoming more complicated, and thus, a multidisciplinary approach to the modelling of the converter is required. In this doctoral dissertation, methods and models are developed for the design of a multilevel power converter and the analysis of the related electromagnetic, thermal, and reliability issues. The focus is on the design of the main circuit. The electromagnetic model of the laminated busbar system and the IGBT modules is established with the aim of minimizing the stray inductance of the commutation loops that degrade the converter power capability. The circular busbar system is proposed to achieve equal current sharing among parallel-connected devices and implemented in the non-destructive test set-up. In addition to the electromagnetic model, a thermal model of the laminated busbar system is developed based on a lumped parameter thermal model. The temperature and temperature-dependent power losses of the busbars are estimated by the proposed algorithm. The Joule losses produced by non-sinusoidal currents flowing through the busbars in the converter are estimated taking into account the skin and proximity effects, which have a strong influence on the AC resistance of the busbars. The lifetime estimation algorithm was implemented to investigate the influence of the cooling solution on the reliability of the IGBT modules. As efficient cooling solutions have a low thermal inertia, they cause excessive temperature cycling of the IGBTs. Thus, a reliability analysis is required when selecting the cooling solutions for a particular application. The control of the cooling solution based on the use of a heat flux sensor is proposed to reduce the amplitude of the temperature cycles. The developed methods and models are verified experimentally by a laboratory prototype.

The effect of cationic-anionic polyelectrolyte multilayer surface treatment on inkjet ink spreading and print quality

The focus of the work reported in this thesis was to study and to clarify the effect of polyelectrolyte multilayer surface treatment on inkjet ink spreading, absorption and print quality. Surface sizing with a size press, film press with a pilot scale coater, and spray coating, have been used to surface treat uncoated wood-free, experimental wood-free and pigmentcoated substrates. The role of the deposited cationic (polydiallydimethylammonium chloride, PDADMAC) and anionic (sodium carboxymethyl cellulose, NaCMC) polyelectrolyte layers with and without nanosilica, on liquid absorption and spreading was studied in terms of their interaction with water-based pigmented and dye-based inkjet inks. Contact angle measurements were made in attempt to explain the ink spreading and wetting behavior on the substrate. First, it was noticed that multilayer surface treatment decreased the contact angle of water, giving a hydrophilic character to the surface. The results showed that the number of cationic-anionic polyelectrolyte layers or the order of deposition of the polyelectrolytes had a significant effect on the print quality. This was seen for example as a higher print density on layers with a cationic polyelectrolyte in the outermost layer. The number of layers had an influence on the print quality; the print density increased with increasing number of layers, although the increase was strongly dependent on ink formulation and chemistry. The use of nanosilica clearly affected the rate of absorption of polar liquids, which also was seen as a higher density of the black dye-based print. Slightly unexpected, the use of nanosilica increased the tendency for lateral spreading of both the pigmented and dye-based inks. It was shown that the wetting behavior and wicking of the inks on the polyelectrolyte coatings was strongly affected by the hydrophobicity of the substrate, as well as by the composition or structure of the polyelectrolyte layers. Coating only with a cationic polyelectrolyte was not sufficient to improve dye fixation, but it was demonstrated that a cationic-anionic-complex structure led to good water fastness. A threelayered structure gave the same water fastness values as a five-layered structure. Interestingly, the water fastness values were strongly dependent not only on the formed cation-anion polyelectrolyte complexes but also on the tendency of the coating to dissolve during immersion in water. Results showed that by optimizing the chemistry of the layers, the ink-substrate interaction can be optimized.

Thermal modelling of commercial lithium-ion batteries

The accelerating adoption of electrical technologies in vehicles over the recent years has led to an increase in the research on electrochemical energy storage systems, which are among the key elements in these technologies. The application of electrochemical energy storage systems for instance in hybrid electrical vehicles (HEVs) or hybrid mobile working machines allows tolerating high power peaks, leading to an opportunity to downsize the internal combustion engine and reduce fuel consumption, and therefore, CO2 and other emissions. Further, the application of electrochemical energy storage systems provides an option of kinetic and potential energy recuperation. Presently, the lithium-ion (Li-ion) battery is considered the most suitable electrochemical energy storage type in HEVs and hybrid mobile working machines. However, the intensive operating cycle produces high heat losses in the Li-ion battery, which increase its operating temperature. The Li-ion battery operation at high temperatures accelerates the ageing of the battery, and in the worst case, may lead to a thermal runaway and fire. Therefore, an appropriate Li-ion battery cooling system should be provided for the temperature control in applications such as HEVs and mobile working machines. In this doctoral dissertation, methods are presented to set up a thermal model of a single Li-ion cell and a more complex battery module, which can be used if full information about the battery chemistry is not available. In addition, a non-destructive method is developed for the cell thermal characterization, which allows to measure the thermal parameters at different states of charge and in different points of cell surface. The proposed models and the cell thermal characterization method have been verified by experimental measurements. The minimization of high thermal non-uniformity, which was detected in the pouch cell during its operation with a high C-rate current, was analysed by applying a simplified pouch cell 3D thermal model. In the analysis, heat pipes were incorporated into the pouch cell cooling system, and an optimization algorithm was generated for the estimation of the optimalplacement of heat pipes in the pouch cell cooling system. An analysis of the application of heat pipes to the pouch cell cooling system shows that heat pipes significantly decrease the temperature non-uniformity on the cell surface, and therefore, heat pipes were recommended for the enhancement of the pouch cell cooling system.