Design of axial-flux permanent-magnet low-speed machines and performance comparison between radial-flux and axial-flux machines

This thesis presents an alternative approach to the analytical design of surface-mounted axialflux permanent-magnet machines. Emphasis has been placed on the design of axial-flux machines with a one-rotor-two-stators configuration. The design model developed in this study incorporates facilities to include both the electromagnetic design and thermal design of the machine as well as to take into consideration the complexity of the permanent-magnet shapes, which is a typical requirement for the design of high-performance permanent-magnet motors. A prototype machine with rated 5 kW output power at 300 min-1 rotation speed has been designed and constructed for the purposesof ascertaining the results obtained from the analytical design model. A comparative study of low-speed axial-flux and low-speed radial-flux permanent-magnet machines is presented. The comparative study concentrates on 55 kW machines with rotation speeds 150 min-1, 300 min-1 and 600 min-1 and is based on calculated designs. A novel comparison method is introduced. The method takes into account the mechanical constraints of the machine and enables comparison of the designed machines, with respect to the volume, efficiency and cost aspects of each machine. It is shown that an axial-flux permanent-magnet machine with one-rotor-two-stators configuration has generally a weaker efficiency than a radial-flux permanent-magnet machine if for all designs the same electric loading, air-gap flux density and current density have been applied. On the other hand, axial-flux machines are usually smaller in volume, especially when compared to radial-flux machines for which the length ratio (axial length of stator stack vs. air-gap diameter)is below 0.5. The comparison results show also that radial-flux machines with alow number of pole pairs, p < 4, outperform the corresponding axial-flux machines.

Design Aspects of Megawatt-Range Direct-Driven Permanent MagnetWind Generators

Direct-driven permanent magnet synchronous generator is one of the most promising topologies for megawatt-range wind power applications. The rotational speed of the direct-driven generator is very low compared with the traditional electrical machines. The low rotational speed requires high torque to produce megawatt-range power. The special features of the direct-driven generators caused by the low speed and high torque are discussed in this doctoral thesis. Low speed and high torque set high demands on the torque quality. The cogging torque and the load torque ripple must be as low as possible to prevent mechanical failures. In this doctoral thesis, various methods to improve the torque quality are compared with each other. The rotor surface shaping, magnet skew, magnet shaping, and the asymmetrical placement of magnets and stator slots are studied not only by means of torque quality, but also the effects on the electromagnetic performance and manufacturability of the machine are discussed. The heat transfer of the direct-driven generator must be designed to handle the copper losses of the stator winding carrying high current density and to keep the temperature of the magnets low enough. The cooling system of the direct-driven generator applying the doubly radial air cooling with numerous radial cooling ducts was modeled with a lumped-parameter-based thermal network. The performance of the cooling system was discussed during the steady and transient states. The effect of the number and width of radial cooling ducts was explored. The large number of radial cooling ducts drastically increases the impact of the stack end area effects, because the stator stack consists of numerous substacks. The effects of the radial cooling ducts on the effective axial length of the machine were studied by analyzing the crosssection of the machine in the axial direction. The method to compensate the magnet end area leakage was considered. The effect of the cooling ducts and the stack end area effects on the no-load voltages and inductances of the machine were explored by using numerical analysis tools based on the three-dimensional finite element method. The electrical efficiency of the permanent magnet machine with different control methods was estimated analytically over the whole speed and torque range. The electrical efficiencies achieved with the most common control methods were compared with each other. The stator voltage increase caused by the armature reaction was analyzed. The effect of inductance saturation as a function of load current was implemented to the analytical efficiency calculation.

Painatuksen UV-stabiilisuuden optimointi fenolihartsatussa paperissa

T��ss�� diplomity��ss�� tutkittiin painetun paperin ja siit�� fenoliformaldehydihartsilla impregnoimalla valmistetun pinnoituskalvon UV-stabiilisuuden parantamis-mahdollisuuksia. Ty��n kirjallisuusosassa k��sitell����n painetun pinnoituskalvon valmistusprosessia ja painatuksen UV-valonkestoon vaikuttavia tekij��it��. Painov��rin pigmentti, sen m����r�� ja k��sittely, painov��rin sideaine sek�� fenoliformaldehydihartsi ja sen lis��aineet vaikuttavat pinnoitetun betonoimisvanerin s����nkesto-ominaisuuksiin. Erilaisilla ep��orgaanisilla valkoisilla pigmenteill�� ja kidemuodoilla on erilainen UV-valonkesto ja taitekerroin. P����llyst��m��ll�� titaanidioksidi esimerkiksi alumiini- tai zirkoniumoksideilla sen UV-valonkestoa voidaan parantaa merkitt��v��sti. UV-hajoaminen voidaan havaita painetun pinnoitteen liituuntumisena. Liituuntumista voidaan pit���� veden ja hapen v��lisen�� reaktiona, jota titaanidioksidi ja UV-s��teily katalysoivat. Sen takia my��s muiden valkoisten ep��orgaanisten pigmenttien ominaisuuksia ja k��ytt���� selvitettiin. Kokeissa k��ytettiin yhdeks���� eri painov��ri��, kahta eri paksuista paperia ja kahta eri tyyppist�� hartsia. Painov��rej�� ohennettiin vedell�� ja paperin painopuolta vaihdeltiin. Kaikissa painatuksissa k��ytettiin kolmea eri rasterointiasteen laattaa, jolloin painov��rin m����r���� paperissa saatiin v��hennetty��. Painetuista papereista mitattiin densiteetti, v��rim����r��, pisara-absorptio vedell�� ja kontaktikulma hartsilla. My��s painov��rin tunkeumaa selvitettiin paperin poikkileikeist�� tehtyjen SEM-kuvien avulla. Painetut paperit impregnoitiin fenoliformaldehydihartsilla kalvoksi. Pinnoituskalvot puristettiin vanerin pinnalle laboratoriopuristimella. Koekappaleet altistettiin UV-valolle, sateelle ja pakkaselle s����kaapissa 400 h ajan, mik�� vastaa noin 1,5 vuotta ulkona Suomen oloissa. Kappaleista mitattiin kiilto, v��rinmuutos ja liituuntuminen. Pinnoitteen liituuntumista tapahtui v��hiten niiss�� koepisteiss��, joissa painatus oli tehty 30 % rasteroidulla laatallaS����nkest��v�� TiO2 osoittautui hyv��ksi, mutta my��s ZnO-pigmentill�� saatiin hyvi�� tuloksia. ZnO-koepisteess�� liituuntumisreaktio ei ole niin voimakkaasti katalysoitu kuin TiO2-koepisteiss��. Paksun paperin painatuspuolella n��ytti olevan merkityst�� s����nkestoon. Huopapuolelle painettuna pinnoitteen liituuntuminen oli v��h��isemp����

Conventional and nonconventional Kramers-Kronig analysis in optical spectroscopy

This Master���s Thesis is dedicated to the investigation and testing conventional and nonconventional Kramers-Kronig relations on simulated and experimentally measured spectra. It is done for both linear and nonlinear optical spectral data. Big part of attention is paid to the new method of obtaining complex refractive index from a transmittance spectrum without direct information of the sample thickness. The latter method is coupled with terahertz tome-domain spectroscopy and Kramers-Kronig analysis applied for testing the validity of complex refractive index. In this research precision of data inversion is evaluated by root-mean square error. Testing of methods is made over different spectral range and implementation of this methods in future is considered.

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.

Effects of interband and intraband impurity scattering on coherence length in the pnictide superconductors

After discovery of cuprates, a search for new high temperature superconducting families began and it led to the discovery of layered pnictide compounds with critical temperatures limited up to ���56 K. Pnictides consist elements from Group V of Periodic Table (nitrogen, phosphorus, arsenic, antimony and bismuth). In this work coherence length ���h in mixed state of pnictide superconductors is calculated numerically. In calculation is taken into account interband and intraband impurity scattering in framework of quasiclassical Eilenberger theory for s�� pairing symmetry. Differences between Ginzburg-Landau and Eilenberger theories is shown and the comparison with existing models is done.

Effects of interband and intraband impurity scattering on coherence length in the pnictide superconductors

After discovery of cuprates, a search for new high temperature superconducting families began and it led to the discovery of layered pnictide compounds with critical temperatures limited up to ~56 K. Pnictides consist elements from Group V of Periodic Table (nitrogen, phosphorus, arsenic, antimony and bismuth). In this work coherence length ��h in mixed state of pnictide superconductors is calculated numerically. In calculation is taken into account interband and intraband impurity scattering in framework of quasiclassical Eilenberger theory for s�� pairing symmetry. Differences between Ginzburg-Landau and Eilenberger theories is shown and the comparison with existing models is done.

Performance Characteristics of an Axial-Flux Solid-Rotor-Core Induction Motor

The integration of electric motors and industrial appliances such as pumps, fans, and compressors is rapidly increasing. For instance, the integration of an electric motor and a centrifugal pump provides cost savings and improved performance characteristics. Material cost savings are achieved when an electric motor is integrated into the shaft of a centrifugal pump, and the motor utilizes the bearings of the pump. This arrangement leads to a smaller configuration that occupies less floor space. The performance characteristics of a pump drive can be improved by using the variable-speed technology. This enables the full speed control of the drive and the absence of a mechanical gearbox and couplers. When using rotational speeds higher than those that can be directly achieved by the network frequency the structure of the rotor has to be mechanically durable. In this thesis the performance characteristics of an axial-flux solid-rotor-core induction motor are determined. The motor studied is a one-rotor-one-stator axial-flux induction motor, and thus, there is only one air-gap between the rotor and the stator. The motor was designed for higher rotational speeds, and therefore a good mechanical strength of the solid-rotor-core rotor is required to withstand the mechanical stresses. The construction of the rotor and the high rotational speeds together produce a feature, which is not typical of traditional induction motors: the dominating loss component of the motor is the rotor eddy current loss. In the case of a typical industrial induction motor instead the dominating loss component is the stator copper loss. In this thesis, several methods to decrease the rotor eddy current losses in the case of axial-flux induction motors are presented. A prototype motor with 45 kW output power at 6000 min-1 was designed and constructed for ascertaining the results obtained from the numerical FEM calculations. In general, this thesis concentrates on the methods for improving the electromagnetic properties of an axial-flux solid-rotor-core induction motor and examines the methods for decreasing the harmonic eddy currents of the rotor. The target is to improve the efficiency of the motor and to reach the efficiency standard of the present-day industrial induction motors equipped with laminated rotors.

CO2 Laser Cutting of Particleboard, HDF and MDF

The purpose of this thesis is to reveal how the laser cutting parameters influence lasercutting of particleboard, HDF and MDF. The literature review introduces the basic principle of CO2 laser, CO2 laser equipment and its usage in cutting of wood-based materials. The experimental part focuses on the discussion and analysis ofthe test data and attempts to draw conclusions on the influence of various parameters, including laser power, focal length of the lens and cutting gas, on the cutting speed and kerf quality. The tested materials include various thicknesses of particleboard, HDF and MDF samples. A TRUMPF TLF2700 HQ laser equipment was used for the experiments. To obtain valid data, the test samples must be completely cut through without any bonding of wood fibre. The maximum cutting speed is linear dependent on the laser power in thecondition that the other parameters are constant. For each thickness of a specific material type, there is a minimum laser power for cutting. Normally, the topand bottom kerf widths increase with the enhancement of laser power. There may be a critical laser power which can generate the minimum cross-sectional kerf width. Lens of larger focal length may achieve higher cutting speed. As the focal length becomes larger, the top kerf width tends to increase while the bottom andcross-sectional kerf widths to the opposite. Of all cutting gases, oxygen can help achieve higher cutting speed. The gas pressure of nitrogen does not seem to have strong influence on the cutting result. Generally, 2 bar air is more preferable for higher cutting speed. For particleboard and MDF samples of larger thickness than 12 mm, 2 bar argon can be used to reach remarkably higher cutting speed than the 5 bar. Generally, the 190.5 mm lens can produce smallest total kerf width. The kerf sides of thicker samples are darker than the thinner ones. The sample darkness tends to be lower as laser power increased. 63.5 mm lens seemed tocause more darkness than other lens. 5 bar cutting gases can produce less dark side kerfs than 2 bar ones. Oxygen normally causes darker kerfs than other gases. No distinct differences were found between nitrogen and argon.

Concentrated Winding Multiphase Permanent Magnet Machine Design and Electromagnetic Properties ��� Case Axial Flux Machine

Concentrated winding permanent magnet machines and their electromagnetic properties are studied in this doctoral thesis. The thesis includes a number of main tasks related to the application of permanent magnets in concentrated winding open slot machines. Suitable analytical methods are required for the first design calculations of a new machine. Concentrated winding machines differ from conventional integral slot winding machines in such a way that adapted analytical calculation methods are needed. A simple analytical model for calculating the concentrated winding axial flux machines is provided. The next three main design tasks are discussed in more detail in the thesis. The magnetic length of the rotor surface magnet machines is studied, and it is shown that the traditional methods have to be modified also in this respect. An important topic in this study has been to evaluate and minimize the rotor permanent magnet Joule losses by using segmented magnets in the calculations and experiments. Determination of the magnetizing and leakage inductances for a concentrated winding machine and the torque production capability of concentrated winding machines with different pole pair numbers are studied, and the results are compared with the corresponding properties of integral slot winding machines. The thesis introduces a new practical permanent magnet motor type for industrial use. The special features of the machine are based on the option of using concentrated winding open slot constructions of permanent magnet synchronous machines in the normal speed ranges of industrial motors, for instance up to 3000 min-1, without excessive rotor losses. By applying the analytical equations and methods introduced in the thesis, a 37 kW 2400 min-1 12-slot 10-pole axial flux machine with rotor-surfacemounted magnets is designed. The performance of the designed motor is determined by experimental measurements and finite element calculations.

Fatigue analysis of brackets with and without scallop for bridge and other deck structure using effective notch stress.

Due to functional requirement of a structural detail brackets with and without scallop are frequently used in bridges, decks, ships and offshore structure. Scallops are designed to serve as passage way for fluids, to reduce weld length and plate distortions. Moreover, scallops are used to avoid intersection of two or more welds for the fact that there is the presence of inventible inherent initial crack except for full penetrated weld and the formation of multi-axial stress state at the weld intersection. Welding all around the scallop corner increase the possibility of brittle fracture even for the case the bracket is not loaded by primary load. Avoiding of scallop will establish an initial crack in the corner if bracket is welded by fillet welds. If the two weld run pass had crossed, this would have given a 3D residual stress situation. Therefore the presences and absence of scallop necessitates the 3D FEA fatigue resistance of both types of brackets using effective notch stress approach ( ). FEMAP 10.1 with NX NASTRAN was used for the 3D FEA. The first and main objective of this research was to investigate and compare the fatigue resistance of brackets with and without scallop. The secondary goal was the fatigue design of scallops in case they cannot be avoided for some reason. The fatigue resistance for both types of brackets was determined based on approach using 1 mm fictitiously rounded radius based on IIW recommendation. Identical geometrical, boundary and loading conditions were used for the determination and comparison of fatigue resistance of both types of brackets using linear 3D FEA. Moreover the size effect of bracket length was also studied using 2D SHELL element FEA. In the case of brackets with scallop the flange plate weld toe at the corner of the scallop was found to exhibit the highest and made the flange plate weld toe critical for fatigue failure. Whereas weld root and weld toe at the weld intersections were the highly stressed location for brackets without scallop. Thus weld toe for brackets with scallop, and weld root and weld toe for brackets without scallop were found to be the critical area for fatigue failure. Employing identical parameters on both types of brackets, brackets without scallop had the highest except for full penetrated weld. Furthermore the fatigue resistance of brackets without scallop was highly affected by the lack of weld penetration length and it was found out that decreased as the weld penetration was increased. Despite the fact that the very presence of scallop reduces the stiffness and also same time induce stress concentration, based on the 3D FEA it is worth concluding that using scallop provided better fatigue resistance when both types of brackets were fillet welded. However brackets without scallop had the highest fatigue resistance when full penetration weld was used. This thesis also showed that weld toe for brackets with scallop was the only highly stressed area unlike brackets without scallop in which both weld toe and weld root were the critical locations for fatigue failure when different types of boundary conditions were used. Weld throat thickness, plate thickness, scallop radius, lack of weld penetration length, boundary condition and weld quality affected the fatigue resistance of both types of brackets. And as a result, bracket design procedure, especially welding quality and post weld treatment techniques significantly affect the fatigue resistance of both type of brackets.

Boundary identification in the domain of a parabolic partial differential equation

Bakgrunden och inspirationen till f��religgande studie ��r tidigare forskning i till��mpningar p�� randidentifiering i metallindustrin. Effektiv randidentifiering m��jligg��r mindre s��kerhetsmarginaler och l��ngre serviceintervall f��r apparaturen i industriella h��gtemperaturprocesser, utan ��kad risk f��r materielhaverier. I idealfallet vore en metod f��r randidentifiering baserad p�� uppf��ljning av n��gon indirekt variabel som kan m��tas rutinm��ssigt eller till en ringa kostnad. En dylik variabel f��r sm��ltugnar ��r temperaturen i olika positioner i v��ggen. Denna kan utnyttjas som insignal till en randidentifieringsmetod f��r att ��vervaka ugnens v��ggtjocklek. Vi ger en bakgrund och motivering till valet av den geometriskt endimensionella dynamiska modellen f��r randidentifiering, som diskuteras i arbetets senare del, framom en flerdimensionell geometrisk beskrivning. I de aktuella industriella till��mpningarna ��r dynamiken samt f��rdelarna med en enkel modellstruktur viktigare ��n exakt geometrisk beskrivning. L��sningsmetoder f��r den s.k. sidledes v��rmeledningsekvationen har m��nga saker gemensamt med randidentifiering. D��rf��r studerar vi egenskaper hos l��sningarna till denna ekvation, inverkan av m��tfel och n��got som brukar kallas f��rorening av m��tbrus, regularisering och allm��nnare f��ljder av icke-v��lst��lldheten hos sidledes v��rmeledningsekvationen. Vi studerar en upps��ttning av tre olika metoder f��r randidentifiering, av vilka de tv�� f��rsta ��r utvecklade fr��n en strikt matematisk och den tredje fr��n en mera till��mpad utg��ngspunkt. Metoderna har olika egenskaper med specifika f��rdelar och nackdelar. De rent matematiskt baserade metoderna karakteriseras av god noggrannhet och l��g numerisk kostnad, dock till priset av l��g flexibilitet i formuleringen av den modellbeskrivande partiella differentialekvationen. Den tredje, mera till��mpade, metoden k��nnetecknas av en s��mre noggrannhet f��rorsakad av en h��gre grad av icke-v��lst��lldhet hos den mera flexibla modellen. F��r denna gjordes ��ven en ansats till feluppskattning, som senare kunde observeras ��verensst��mma med praktiska ber��kningar med metoden. Studien kan anses vara en god startpunkt och matematisk bas f��r utveckling av industriella till��mpningar av randidentifiering, speciellt mot hantering av olinj��ra och diskontinuerliga materialegenskaper och pl��tsliga f��r��ndringar orsakade av ���nedfallande��� v��ggmaterial. Med de behandlade metoderna f��refaller det m��jligt att uppn�� en robust, snabb och tillr��ckligt noggrann metod av begr��nsad komplexitet f��r randidentifiering.

Heat-transfer improvements in an axial-flux permanent-magnet synchronous machine

Axial-flux machines tend to have cooling difficulties since it is difficult to arrange continuous heat path between the stator stack and the frame. One important reason for this is that no shrink fitting of the stator is possible in an axial-flux machine. Using of liquid-cooled end shields does not alone solve this issue. Cooling of the rotor and the end windings may also be difficult at least in case of two-stator-single-rotor construction where air circulation in the rotor and in the end-winding areas may be difficult to arrange. If the rotor has significant losses air circulation via the rotor and behind the stator yokes should be arranged which, again, weakens the stator cooling. In this paper we study a novel way of using copper bars as extra heat transfer paths between the stator teeth and liquid cooling pools in the end shields. After this the end windings still suffer of low thermal conductivity and means for improving this by high-heat-conductance material was also studied. The design principle of each cooling system is presented in details. Thermal models based on Computational Fluid Dynamics (CFD) are used to analyse the temperature distribution in the machine. Measurement results are provided from different versions of the machine. The results show that significant improvements in the cooling can be gained by these steps.

Natural freeze crystallization : purification of aqueous Na2SO4 solutions with cold air

The aim of this thesis was to examine efficiency of freeze crystallization and eutectic freeze crystallization in purification of wastewater by imitating natural freezing. In addition, a mathematic model based on heat transfer to determine ice thickness and ice growth rate was examined. Also, the amount of sodium sulfate crystallized at the eutectic point was under investigation. In literature part, advantages and applications of the freeze crystallization are discussed, and possibility to apply it in Northern hemisphere winter weather conditions is under study. Furthermore, main sources of sodium sulfate from Finnish industries are described. The experiments were carried out in modified chest freezer, where a fan was placed in order to obtain laminar air flow inside. Picolog PT-104 data logger was used to monitor temperature changes in the salt-water solution, and constant temperature was maintained in the crystallizer with Lauda RP 850 thermostat. The impurity of formed ice layer was determined by weighing ice samples after experiment and again after 24 hours drying to full dryness in oven. Volume of salt-water solution was also measured after experiment. The highest purity of formed ice layer was obtained with small temperature difference and with long freezing time. On the other hand, the amount of crystallized sodium sulfate was its greatest with long freezing time and higher temperature difference. The results obtained by the mathematic model and empirical results did not differ significantly in most of the experiments. However, the difference increased when salt-water mixture reached its eutectic point, leading to simultaneous ice and salt crystallization. Eutectic point was reached only with the highest salt concentration with one exception. In these cases, calculated values were in many cases greater than the experimental ones. In winter weather conditions freeze crystallization is cost-effective wastewater treatment method and rather simple. Nonetheless, the efficiency and separation rate are strongly depended on ambient temperature and its changes