AP05_4.vp 1 Introduction Most materials of sculpture works, paintings, plasters, stuccos, etc., have absorbed a moisture content that settles at an equilibrium level corresponding to the relative humidity and temperature of the surrounding air (Cassar, 1995). Wood is one of the most sensitive materials, furthermore with con- siderably different mechanical properties and extensibility in each of its three primary axes. Variations in moisture content therefore result in stresses leading to deformations or even to cracks opening in the exhibits (Camuffo, 1998). The internal environment standard in leading galleries and museums in- volves keeping constant both the interior temperature and the air humidity in the exhibition rooms, but due to impracti- cable costs, such air-conditioning cannot realistically be im- plemented as far as remote historical objects are concerned (Cassar, 1993). Owing to this, an invaluable part of European cultural heritage is more or less exposed to the damaging in- fluence of varying moisture caused by air humidity fluctua- tions (Kotterer, 2002). However, when it is not possible to keep constant both air temperature and humidity, the primary goal is to avoid dimensional extensions or shrinking of the exhibits arising from variations in air humidity and temperature. In principle, it is possible to accept modest temperature changes if they are compensated simultaneously by some adequate air humidity adjustments. In this way it is possible to protect the exhibits against the impact of moisture change and to ful- fil the crucial conservation requirement, particularly with respect to the usually anisotropic character of these phenom- ena. Apparently, this kind of compensation can meet the preservation demands only if both the temperature and the air humidity changes are smooth and slow enough for the internal atmosphere to remain in an almost steady state and the absorbed moisture in the exhibit materials to remain constant. 2 Air humidity and equilibrium moisture content In view of the preservation demands it is a hard task to de- fine unambiguously the proper parameters of the desirable environment. First of all, it is not the visitor s comfort, but rather the benefit of the exhibits that should be the priority in exhibition rooms (Camuffo et al., 2002). Recent research in preventive conservation has introduced the complex con- cept of the microclimate of a room, in which monuments or exhibits are placed (Camuffo, 1998). Among the numerous parameters of the microclimate, air humidity and tempera- ture are considered as primary and as the most important attributes of the microclimate (Cassar, 1993). Particularly in remote sites of deposits, where neither heating nor air han- dling devices are in operation, the humidity impact is the most dangerous exposure from the preservation point of view. The decisive role of moisture sorption impact is typical for most of the materials that artistic works are made of, i.e. wood, paper, parchment, leather, ivory, bone, paints, plas- ter, stucco or stones containing abundant clay minerals, etc. These materials absorb a specific amount of water and are rather sensitive to variations in it. The steady value of this amount, corresponding to the surrounding air humidity and temperature, is called the equilibrium moisture content (EMC), usually expressed as the ratio of the mass of water to the mass of anhydrous material. After the ambient temperature or humidity have changed their values, the absorbed moisture content changes accordingly (Massari, 1993). An increase in EMC increase is then followed by swelling of the material and, conversely a decrease results in contraction. Due to the non-isotropic character of these size changes, harmful defor- mations or destructive cracks appear as a result. The material extensions resulting from growing EMC are relatively high, and it is important to note that this extension is much higher than just the thermal extension of the dry wood. In fact the expansion phenomenon is somewhat more complex. For example, a rise in temperature induces thermal extension of dry wood alone, but consequently it results in a drop in relative humidity and therefore also in a drop in EMC, which brings about a material contraction, and vice-versa (Camuffo, 1998). In this way, thermal expansion and EMC contraction are of opposite character and the shrinkage is partially mitigated by the expansion. However, the dimension change due to relative humidity is largely dominant, since only the temperature expansion itself is more than ten times weaker than that of relative humidity (Kowalski, 2003). © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 5 Czech Technical University in Prague Acta Polytechnica Vol. 45 No. 4/2005 Model Based Control of Moisture Sorption in a Historical Interior P. Zítek, T. Vyhlídal This paper deals with a novel scheme for microclimate control in historical exhibition rooms, inhibiting moisture sorption phenomena that are inadmissible from the preventive conservation point of view. The impact of air humidity is the most significant harmful exposure for a great deal of the cultural heritage deposited in remote historical buildings. Leaving the interior temperature to run almost its spontaneous yearly cycle, the proposed non-linear model-based control protects exhibits from harmful variations in moisture content by compensating the temperature drifts with an adequate adjustment of the air humidity. Already implemented in a medieval interior since 1999, the proposed microclimate control has proved capable of permanently maintaining constant a desirable moisture content in organic or porous materials in the interior of a building. Keywords: equilibrium moisture sorption, humidity control, non-linear cascade control. 3 Equilibrium moisture content models For each of the considered materials the equilibrium moisture content settles at a level appropriate to the ambient air humidity and temperature. Although the EMC levels are different for various materials, the following properties are common to all of them � the EMC always increases with growing � and decreases with growing T, � the EMC value is much more sensitive to changes in air hu- midity than to variations in temperature. The relationship between EMC as u, and both air temper- ature T and relative humidity �, u T� �( , )� has been fitted by several formulae developed for various areas of application. Usually this relationship is plotted in the coordinates � and u, as so-called sorption isotherms, with temperature considered as a parameter. In particular the mathematical models by Day and Nelson and by Simpson (Ball et al., 2001) were found to fit well the experimental data for professionally evaluating EMC. The Day and Nelson model modified by Ball et al. (2001) is of the form u A T A T k T DN K B T DN K K DN K � � �� � � � � � log( ) ( ) , ( ) exp( ) ( ) 1 1 1 � 300 3002 � � � � � � � � � � � , ( ) exp( ) ,B T k T DN K K (1) where TK is the temperature in degrees Kelvin, � is the rela- tive humidity expressed as the mass ratio and k1, k2, �, � are model parameters. The Simpson model, in the version for wood moisture evaluation is as follows u W K K K K K K K K K K K K � � � � �1800 1 1 0 0 0 1 0 2 1 2 2 0 1 0 2 1 2 2 � � � � � � � � �� � � � , (2) where the parameters W, K0, K1, K2 are given as quadratic polynomial functions of temperature T. Using models (1) or (2), high accuracy can be achieved in fitting the experimental data. However, for the microclimate control idea proposed below these models are less suitable because their derivatives result in fairly complicated formulae. That is why we chose the following Henderson model with only three parameters u A T B T C � � � � � � � � � � ln( ) ( ) ( , ) 1 � �� , (3) to derive the following model based controller of moisture sorption. In this model the relative humidity is considered as the mixing ratio � � 0 1, , temperature T and parameter B are in °C, C is a positive dimensionless exponent, A is a pa- rameter in °C�1 and u is the dimensionless ratio of moisture mass content to the mass of dry material. As regards the function �(.) the mutual dependence between the variables � and T should be noted. If x is the absolute humidity in kg of water per one kg of dry air, its value, corresponding to a certain given � and T, results, e.g., from the Magnus relationship (Camuffo, 1998), as follows x aT b T� � � �3 795 10 103. ( )� , (4) where a � 7.5 and the additive temperature constant b � 237.3 °C. Only from x can the actual water vapour content be assessed for a given air volume. On the other hand the rel- ative humidity changes inversely to the temperature change if the absolute humidity x is maintained at a constant level. Then for two states �1, T1 and �2, T2 of humid air with the same x it holds � � 2 1 10 10 1 1 2 2 � � � aT b T aT b T ( ) ( ) . (5) Hence, if the humid air temperature has changed without changing the water vapour content x, its relative humidity moves in the opposite way. For example, a temperature drop under x � const. necessitates the relative humidity increase given by (5). This phenomenon can be observed in simulta- neous records of � and T, where most of the relative humidity fluctuations are usually due to temperature changes, while only a minor part is due to the actual change in water vapour content. Another consequence can be derived from the Magnus formula. If a considerable change in relative humidity from �0 to � has occurred under x = const., the corresponding temperature change is relatively small, given as follows T T b T a0 0 0� � � �(log log )� � , (6) where, unlike (3), the decimal logarithm is used. This means that, for example, a drop of relative humidity from 60 to 40 per cent, where x does not change, corresponds to a tem- perature increase of 6 °C only. Since the sorption isotherms describing the relationship u T� �( , )� are rather close to each other, the small temperature increments given by (6) re- sult in the property that the lines of x � const. are of a rather similar shape to the isotherms in the diagram of u T� �( , )� and are relatively close to them. A comparison of several iso- therms with the lines of x � const. can be seen in Fig. 2. Model (3) is marked out with advantageous analytical properties, particularly with regard to its derivatives. The model parameters A, B, C can be assessed to fit a set of data obtained either from experiments or from a more pre- cise EMC model (1) or (2). For a given set of data �i, Ti, ui, i � 1, 2, …, n the most fitting parameters A, B, C can be identi- fied by means of the least-square approach. For wood the parameter C is about 0.65, B is from 260 to 310 °C and the value of A is about 0.28 °C�1 if u is expressed in kg/kg. In this way either the experimental data or the more complicated formulae (1), (2) can be fitted by the relatively simple model (3) with a good coincidence over sufficiently broad intervals, i.e. for air temperature from 5 to 40 °C and for air humidity from 10 % to 90 % in the controlled interior. An example of the attainable conformity of the Simpson expression of a sorption isotherm with the Henderson model is shown in Fig. 1. On the other hand, it is obvious from the nature of function (3) that it cannot fit EMC for relative humidity values close to 100 %, since then ln( )1 � � ��� . Obviously, due to varying sorption properties, parameters A, B, C result in dif- ferent values for various materials. In order to fit the three parameters of the Henderson model to the more precise Simpson model, the iterative least square method was used. Applying this method to fit the points of the characteristics within the intervals T � 0 25, °C 6 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 4/2005 Czech Technical University in Prague and � � 10 90, %, which define the region of interest from a practical point of view, the following parameters result A � 6.962 10�5, B � 305.7 °C and C � 0.64. The accordance of these two models can be seen in the comparison of one of the isotherms in Fig. 1. 4 Equal-sorption principle of air humidity control Let us assume the well-being of the preserved exhibits as the priority in setting up the parameters of the exhibi- tion room microclimate. The following facts and claims are then to be taken into account in controlling the internal air parameters � The impact of air humidity changes is more significant than the impact of usual temperature fluctuations. � The primary harmful impact on the state of the exhibits is due to the moisture content absorbed in their porous or organic materials. � Changes in moisture content have the crucial harmful im- pact, bringing about an anisotropic swelling or shrinking that results in deformations or cracks in exhibits. � Furthermore, a higher level (> 12 %) of moisture content becomes favourable for microbiologically harmful organisms to attack the surface of the exhibit. Hence, it is not the air temperature and humidity, but the moisture content in the absorbent materials that is the deci- sive parameter to be kept as constant as possible in order to avoid any swelling or shrinking of the exhibits. However, this parameter can hardly be considered as a usual controlled vari- able, since its continual measurements are not available. First of all, the harmful impact is brought about by the anisotropic behaviour of the exhibit materials. On the other hand, it is important to realize that an exactly constant EMC can be kept only in one sort of exhibit due to the differences among the sorption isotherms for various materials. Moreover, the idea of compensating the spontaneous temperature fluctua- tions by humidity adjustments is feasible if, and only if, the variations of the interior temperature are smooth and slow enough, i.e. sufficiently close to a steady state. It should be emphasized that the moisture content is only slightly sensitive to temperature variations, but much more dependent on the air humidity. Let the following example be mentioned: A great temperature increase, say from 5 °C to 40 °C, causes approximately the same change in moisture content in wood as results from a relative humidity drop as small as about 6 %, (Camuffo, 1998). Although no direct measurement of EMC is available, the Henderson model allows us to express its value analytically from the air humidity and temperature measurement, pro- vided that the state of the environment is all the time close to a steady state. With the help of this model, a control action can be assessed protecting exhibits made of the selected material from harmful EMC changes. Let us assume that the tempera- ture change is slow and smooth, and let a simultaneous air humidity correction be provided in order to keep zero incre- ment of the moisture content. In other words, although both temperature T and relative humidity � of the air change, the moisture content in a specific material can be kept constant. If small increments of the state variables are considered, the compensating �� correction results from the derivative of model (3), i.e. from the condition (Zítek and Němeček, 2004) � � � � �u T T� � � � 0 (7) Using the Henderson model, the derivatives in this equa- tion can be obtained from (3) in an analytical form. For the air temperature, the derivative is always negative �u T A C B T C C � � �� �� � � � � � � ln( ) ( ) ( ) 1 01 (8) (0< � < 1) and, conversely, for the humidity the derivative is always positive � � � u B T A C A C C � � � � � � � � �� �� � � � � �1 1 1 0 1 ( ) ln( ) ( ) . (9) The opposite signs of the two derivatives correspond to the opposite character of temperature and relative humidity increments in case of unchanged absolute humidity x given as the ratio of water vapour mass to the mass of dry air. Following the requirement to keep the moisture content constant in the selected material, it is necessary to provide a compensating air humidity correction � � corresponding to the admissible temperature increment �T according to the following ratio of the derivatives © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 7 Czech Technical University in Prague Acta Polytechnica Vol. 45 No. 4/2005 � [%] Fig. 2: Comparison of the shapes of isotherms (solid) and the lines of x � const. (dashed) (using the Magnus and Hen- derson models) � [%] Fig. 1: Comparison of fitting the Simpson model by the Hen- derson model, with parameters A � 6.96 10�5, B � 305.7, C � 0.651, isotherms T � 15 °C (solid – Simpson, dashed – Henderson) � � �� � �� � � � u T u T K T TC( , ) , (10) where the compensation gain K TC( , )� , given by K T B TC ( , ) ln( )( ) � � � � � � � 1 1 (11) is a function of both the temperature and the air humidity. However, from the three parameters of the Henderson mod- el, K TC( , )� is dependent only on B. Obviously, due to the dominant role of B in the denominator value, the relative change of with respect to the temperature variations is fairly low. The increment-proportional character of relationship (11) does not guarantee that a constant desired EMC value u uD� can be maintained by this kind of control. In order to prevent the moisture content from drifting from uD, the in- tegration action has to be added into the controller and the following feed-forward-feedback control action results � � � �� � d d d d � � � ( ) ( ), ( ) ( ) ( ), ( ) t t K t T t T t t K u T t tC I D� � � � � , (12) where the feed-forward gain K TC( , )� , compensating the temperature trend, is given by (11), the unavailable EMC value is computed from the � and T measurements by means of model (3), and KI is an integration gain. This controller is considered as a master controller in a cascade scheme, see Fig. 3, where the slave is a humidity control with a variable desired � given by (12). With regard to the relative hu- midity to be kept in the controlled environment (within the interval � � 40 60, %), the weak variability of the feed- -forward gain can be neglected and a constant coefficient K TC( , )� � � 0.0013 °C �1 can be considered. As regards the integration gain KI, its value should be assigned low enough to prevent the control process from obtaining an oscillatory character. In any case, it is desirable to set KI to achieve a suffi- cient filtering effect in the day cycle fluctuations of the relative air humidity. In general, the dynamics of the interior air humidity can be considered subject to the model � � � �V x t t Q x t x t S x t x t q tI L E I x W W I M� � d d ( ) ( ) ( ) ( ) ( ) ( )� � � � � (13) where V is the interior volume, � air density, Q L the in-leak- age of the ambient air, �x the effective diffusion coefficient between the walls and the internal air, SW is the entire wall surface and qM the water input provided by the humidifier ( qM > 0) or dehumidifier (qM < 0). The values of absolute humidity of the external and internal air and the walls respec- tively, are xE, xI, xW, where xI is computed from the measured interior temperature T(t) and the interior relative humid- ity by means of the Magnus formula (4). The in-leakage flow of the outdoor air is usually referred to as being ap- proximately 25 to 30 % of the internal volume per hour ( ( ) ( . . ) )Q VL � � �0 25 0 3 1to h . The moisture diffusion from the walls is highly variable. According to our measurements, the moisture release from damp walls may represent as signifi- cant a water source for the interior as window leakage does. The last term qM(t) of (13) is the water transfer amount. This control actuating effect has to be powerful enough to counter- act even the highest demands of moisture transport. The time constant of the humidity process (13) in the interior is as follows � � � V Q KL D . (14) Since the on-off control is commonly used for the humidi- fying device (i.e. as the slave controller) this time constant is valid approximately also for the humidity-controlled inte- rior. Unlike the master controller output, equation (13) is expressed for absolute humidity x, but the interior output measurement is again provided for the relative humidity. Hence the two conversions � � x and x � � cancel each other in the overall plant model. The only conversion left is given by the Henderson model. After linearizing this model as in (7) the following characteristic equation of the control loop is obtained � s s K u s s K kI I H 2 2 0� � � � � � (15) and the non-oscillating character of the control process is obtained if the integration gain satisfies the condition K kI H� �( )4 1 . In considering the usual parameters for an interior about 100 m2 in size or larger, this involves adjusting the integration gain approximately as low as K I � �1 1h . 5 Application case study Within the framework of the research project several his- torical interiors in the Czech Republic are being monitored. One of them is St. Václav Chapel in Plasy Monastery in West Bohemia. As in the other sites the selected interior microcli- mate is monitored by dataloggers measuring the tempera- tures and relative humidities each 30 minutes. In the Plasy application, a model of the interior EMC pro- cess has been identified from the collected data, and an equal-sorption humidity controller has been designed. Corre- lating the interior and exterior absolute humidity (computed by means of the Magnus formula) and using the least-square- -based identification method, the time constant results at � 2.44 hours. The interval in Fig. 5 was selected for identifi- cation in order to obtain the day cycle fluctuations of the abso- lute humidity with the least influence of long-term transient phenomena. Due to this selection, wall moisture release, which has much slower dynamics than the variations caused by infiltration (the dynamics of which we want to model), has a minimum effect on the system dynamics. A comparison of the interior absolute humidity simulated by model (13) (the exterior absolute humidity is considered as the model input) and the real interior absolute humidity can be seen in Fig. 5. 8 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 4/2005 Czech Technical University in Prague Fig. 3: Scheme of microclimate control Even though there are considerable discrepancies between the modeled and measured data (the process itself is non- linear and its dynamics is much more complicated than the dynamics of model (13)), the model fits the character of the dynamics sufficiently. Using the model to transfer the measured exterior hu- midity into the interior humidity, and considering the inte- rior temperature to be the temperature measured on the same time interval as the exterior humidity, see Fig. 5, the simulation of the control in keeping EMC has been per- formed. The results of the simulation are shown in Fig. 8. The parameters of the master controller have been chosen KC � 0.0011 °C �1, KI � 0.007 hour �1. The maximum dehu- midifier performance needed for counteracting the worst moisture sources is approximately 500 g of water per hour. 6 Conclusions A novel microclimate control scheme has been proposed, for protecting preserved exhibits from changes in equilibrium moisture content in the materials they are made of. Although the EMC as an alternative controlled variable is unavailable for continual measurement, changes in EMC can be esti- © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 9 Czech Technical University in Prague Acta Polytechnica Vol. 45 No. 4/2005 Fig. 4: Temperature in St. Václav’s Chapel and ambient tempera- ture measurements Fig. 5: Relative humidity in St. Václav’s Chapel and ambient rela- tive humidity measurements Fig. 6: Process/model comparison, model time constant � 2.44 hours and static gain Q Q KL L D( ) .� � 014 (absolute humidity mean values: xI � 5.74 g/kg, xE � 6.45 g/kg) Fig. 7: The ambient relative humidity and interior temperature used as the model inputs in simulating control on EMC Fig. 8: Results of simulating the control on EMC in St. Václav’s Chapel. Upper: un – EMC if no control is used (computed from the measured data), uD – desired value of EMC, uc – EMC if control is used (resulting from simulation). Lower: the amount of water to be extracted mated by a non-linear model (3) from the humidity and tem- perature measurements. Compensating the natural interior temperature variations simultaneously by the appropriate adjustment of air humidity proves efficient in maintaining a reference EMC constant during the entire annual weather cycle. Basically, the non-linear model-based microclimate control is marked out with a feed-forward character since the changes in the thermodynamic air state variables are much faster than the sorption phenomena, and therefore they in- hibit the possible sorption process. The proposed moisture content control scheme has already been implemented and has been in successful operation in The Holy Cross Chapel at Karlštejn Castle in the Czech Republic (Zítek and Němeček, 2004) since 1999. Incited by the alarming deterioration of the collection of 29 precious medieval paintings by Master Theodoricus, the implementation was based on a thorough investigation of the interior microclimate. In heavy inertial structures, like those in this implementation, the tempera- ture course can be allowed to run almost its natural yearly cycle, without considerable heating influence, since the temperature does not drop below six or seven °C. However, in mansions and similar buildings which do not have such thick walls, the winter course of the heating-unaffected interior temperature would drop to zero, or even below freezing point during the winter season, and it is essential to maintain the temperature safely above zero. Hence the temperature cycle has to be corrected by a low-level heating during the winter season in buildings where such low temperatures occur. On the other hand, moisture control is able to safely protect the interior from condensation effects if the interior temperature does not drop below 5 °C. Although the proposed control scheme is based on the EMC model for a single material selected as decisive for preservation purposes, the differences between the sorption isotherms of various materials affect the control action only weakly, also with regard to the low level of humidity correc- tions needed. Humidity compensation therefore works well not only for the selected material but also for other materials. Although the sorption differences between materials are ex- pressed in three parameters A, B, C, only B (in °C) influences the main part of the control action (13). Only this property of the Henderson model turned out to be significant in helping to show that the control actions (13) computed for various materials would not be very distinct. The value of B is rela- tively high (more than 200 °C) and therefore the compensa- tion control parameter KC(T, �) is rather low and makes the humidity correction relatively small. In this way the sorption differences between various materials have an inconsiderable effect on the resulting control action (13). On the other hand, the following two factors are to be emphasized: � The equilibrium assumption. It is of primary importance to avoid any rapid air temperature and humidity changes in the interior. � Providing sufficient internal air circulation. The control of the internal air humidity will have little effect without ensuring sufficient air circulation, particularly near the surface of the walls. It should be emphasized that the most important advan- tage of equal-sorption control is not economy in power consumption but, above all, the careful character of its opera- tion from the preventive conservation point of view. Due to the lower temperature and controlled humidity during the winter season the exhibits do not suffer such exposure to high gradients of both temperature and humidity as they do in the case of usual air-conditioning, or in historical buildings where no technical care of the internal environment is provided. 7 Acknowledgement The presented research has been supported by the Grant Agency of the Czech Republic under Project No. 101-03-1365 and by the Ministry of Education of the Czech Republic under Project LN00B096. References [1] Camuffo, D.: Microclimate for Cultural Heritage. Elsevier Science Ltd., Amsterdam, London, 1998 [2] Camuffo, D., Bernardi, A., Sturaro, G., Valentino, A.: “The Microclimate inside the Pollaiolo and Botticelli Rooms in the Uffizi Gallery, Florence.“ Journal for Cul- tural Heritage, Vol. 3 (2002), p. 155–161. [3] Cassar, M.: Environmental Management. Guidelines for Museums and Galleries. Routledge, London 1995 [4] Cassar, M.: “A Pragmatic Approach to Environmen- tal Improvements in the Courtauld Institute Galleries in Somerset House.” ICOM Committee for Conservation, Vol. 2 (1993), p. 595–600. [5] Kotterer, M.: “Research Report of the Project EU 1383 PREVENT”, Museum Ostdeutsche Galerie Regensburg, 2002. [6] Kowalski, S. J.: Thermomechanics of Drying Processes. Springer, Berlin, 2003. [7] Massari, G., Massari I.: Damp Buildings. Old and New. ICCROM, Rome, 1993. [8] Ball, R. D., Simpson, I. G. and Pang, S.: “Measurement, Modelling and Prediction of E.M.C. in Pinus Radiata Heartwood and Sapwood.” Holz als Roh- und Werkstoff 59 (2001), Springer-Verlag 2001, p. 457–462. [9] Zítek, P., Němeček, M.: “Stabilizing Microclimate of His- torical Interior by Moisture Sorption Control.” IEEE Conference on Math. Modeling, Automation and Ro- botics 2004, Szczecin, Poland (to appear). Prof. Ing. Pavel Zítek, DrSc. e-mail: zitek@fsid.cvut.cz Ing. Tomáš Vyhlídal, Ph.D. e-mail: vyhlidal@fsid.cvut.cz Centre for Applied Cybernetics Institute of Instrumentation and Control Engineering Czech Technical University in Prague Faculty of Mechanical Engineering Technická 4 166 07 Praha 6, Czech Republic 10 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 4/2005 Czech Technical University in Prague