hungarian journal of industry and chemistry vol. 49(2) pp. 1–1 (2021) hjic.mk.uni-pannon.hu doi: 10.33927/hjic-2021-12 editorial preface to the special issue dedicated to the manufacturing_2021 conference istván gyurika1 1research centre for engineering sciences, university of pannonia, egyetem u. 10, veszprém, 8200, hungary the association of mechanical engineers in hungary, “gte”, has a long tradition in organizing a nation-wide, scientific-technical conference devoted to manufacturing. during the 29-30th october 2021, the university of pannonia offered to host the 24th gte manufacturing conference, where almost 50 presentations were delivered by engineering students with phd or post-doc level. participants from industry had also a good opportunity to present their new products and services, thus a fruitful collaboration platform had taken place at the premises of the university of pannonia. the highly appreciated hungarian journal of industry and chemistry is an excellent medium to convey the novelties and new scientific-technical achievements demonstrated during the conference. the phd students are also grateful to hjic for the chances to generate high-level publications at international level. the presentations were mainly offered in english language, giving additional opportunities for phd students to practice oral presentations of their scientific work and achieved results. all of the topics address one or more of the offered areas of manufacturing, including prodcorrespondence: gyurika@almos.uni-pannon.hu ucts, production, logistics, sensor-, monitoring-, roboticsthematic areas, as well as economic and human aspects, e.g., operator cooperation with robots. the review committee generated detailed feedback to the authors to show higher level of articulation of their ideas and results, and to maintain the scientific high level of the hjic journal. the selected papers are included in this special issue of hjic. the gte association will continue to offer its members opportunities to share up-to-date technical and scientific information and build a reliable social-technical network in the new ecosystem needed for the transition from indutry 4.0 to industry 5.0. we are looking ahead and inviting for the quad-centenarian, 25th jubilee manufacturing conference. istván gyurika chairman of gte’s division of manufacturing systems on behalf of the manufacturing: 2021 organising committee. https://doi.org/10.33927/hjic-2021-12 mailto:gyurika@almos.uni-pannon.hu hungarian journal of industrial chemistry veszprem vol. 30. pp. 281-284 (2002) stochastic model of first order chemical reactions taking place in particles t. buckle, a. bezegh1 and c. mihalyk62 (research institute of chemical and process engineering, university of kaposvar, egyetem u. 2., veszprem, h-8200, hungary 1department of environmental economics and technology, budapest university of economic sciences and public administration, fovam ter 8., budapest:, h-1093, hungary 2department of mathematics and computing, university ofveszprem, veszprem, p.o. box 158, h-8201, hungary) received: october 17, 2002 experimental observation had been made in the 1950s, that the application of the fluidisation method increases the formation rate and homogeneity of the iron oxide red product. results of the interpretation of this observation by the recently developed mathematical models had been presented. the interpretation is based on the stochastic treatment of the virtual rate constant due to the random temperature differences of the particles at heating processes. the stochastic model is derived from the first order kinetic equations assuming uniform distribution for the rate constant and includes the momentum equations. discussion of the model is described for different sets of parameters. keywords: stochastic model, particle dispersion processes, iron oxide red production, conversion rate, product homogeneity, momentum equation, particulate system introduction during the years of the '50s, an intensive research was lead by prof. k. polinszky, to study the properties of inorganic pigments and their production methods. this had been summarised elsewhere [1]. it is a crucial question what result can be obtained, if we treat the virtual rate constant as a stochastic parameter, due to the random temperature differences of the particles in thermal processes? the models of the first order processes are presented in the following: a group of inorganic pigments are produced by heat treatment of solid particles. for example: minium, from lead oxide, chromium oxide green, from potassium dichromate, iron oxide red from iron oxide black and yellow and from ferrous sulphate heptahydrate. at that time one of the authors, t. blickle, dealt with mathematical modelling and the application of fluidisation, respectively. as it was experienced, in a fluidised bed the conversion rate was higher than that in a rotating kiln, or in a steady layer. it was impossible to explain these phenomena by the available mathematical models. the mathematical models of particle dispersion processes, including stochastic processes as well, had been studied in the recent and past decades [2-7]. here we present the results of the interpretation of observations, gained by the recently developed mathematical model. stochastic model the kinetic equations: de~ * 1-=-k.c. dt j j (l) where c ~ is the concentration of the initial component, j ci is the concentration of the final component in the j4.11 particle, a. is the stochiometric factor. the virtual rate constant, k. is a stochastic variable with the probability density function, 8( ld and its ki values are assigned randomly to the j-dt particle by drawing according to this density function, i.e. ~ is a realisation of k. in the moment of the start, the 282 concentration frequency density function of the particles in the system is given by: n[c(o), c * (0),0] (3) from eqs.(l) and (2), the transition functions can be obtained: (4) (5) density function n[c(o),c*(o),o] and eqs.(4) and (5) yield the following model of the integral transformation equation: k [c(t)ac * (t)(ekt -1),] n~(t),c * (t),t ]= jn c * (t)ekt ,0 . (6) k., • ekt b(k)dk the momentum equation is: 1 1 m 1 (t) = j j c1 (t)n[c(t),c * (t),t]dc(t)dc * (t) (7) 00 using the eq.(6), one can get: 1 1 k 1 ) [c(t)ac * (t)(ekt -1),] 11 c (t n · m 1 (t) = j j j c* (t)ekt ,0 (8) ook .. kt * · e b(k )dkdc(t)dc (t) by applying eqs.(4) and (5), and after transformation: 1 1 ku[c(o) + c * (o)il(le-kt)]1 • m (t) = if i (9) 1 0 0 t .. · n[c(o).c * (o),ojb(k)dkdc(o)dc * (0) assuming that the distribution of k values are uniform in the range of km ku: kk'" +ku ai. k l1k ,. ak 2 _,oj{.= .y-k,n, ak = .jli, a;;= k values k.t and ku ar~ the lower and upper limits of a uniform distribution, k is the average value of the uniformly distributed virtual rate constant, l:jc is the width and ty~; is the standard deviation of its range. while a 1 is the coefficient of variation (relative standard deviation} of the uniformly distributed virtual rate constant. until now it was assumed, that the values of ki,, which had been drawn at t = 0 stay unchanged during the entire process. however, in many cases there are such effects in the system which make it necessary to draw again periodically, after each e. here, when the elapsed time, r between two drawings is 0 < r 50 and t time is: t=l·e>+-r (12) where lis the number of drawing. if c(o) = 0 and -r * 0, then: 1e-k't 2 . fe .juak [ sinh fe .jl2a k ] 1 c(t) = ilc * (0) 2 (13) ·[sinhf.r~l k jiiak 7:-2 if c(o) = 0 and -c = 0, then: [ [ sinh fe ..jyia k ] 1 ] c(t) = ilc * (0) 1e -kt ..fij? (13a) ke 12o-" 2 if ti(o) = u*2(0) = 0 and -c * 0, then: if rl(o) = u*2(0) = 0 and 7: = 0, then: o.s j.() 1.5 2.0 25 3.0 3.5 4.0 45 5.0 kt fig. i plot of conversion against time, at different relative ranges of the virtual rate constant ( sinh 'kejuak 'jt kejlia-k -2 'j/ cjz (t) = k c * (o)e-zkt [ . .j12 " ]~ smhkg--a 2 k (14a) because limx-to sinh x = 1 if r and 0 converges to 0, x while t is constant: (15) (16) discussion of the stochastic model eqs.(13a) and (14a) are used. r = 0; l = 1; c(o) = 0 (in this case t = 8 ). while discussing the function c(t) +kt let ak be 0; 0.25; ..1-c* (0) 0.5. (see fig.l.) t = 0; l = 1, c(o) = 0, tt\0) = 0'*2(0) = 0. here we discuss cr(t) = ~(t) +kt relationship, if 8k = 0.25; 0.5. c(t) (fig.2.) i rr(ij f(ll ojj 1.0 2.0 3jj 4.0 283 i i i i i ! kt s.u i fig.2 plot of relative concentration variance against time, at different relative ranges of the virtual rate constant i 8,0 ' i i i 7,0 i 6/.j i ~) l i.e • (0.) ~ 4.'j i i 2/! fig.3 plot of conversion against the relative range of the virtual rate constant, at different times here the c(t) is studied at different (0.8, 0.9, 0.95) ac * (0) values of the necessary timet (if r= 0, l = 1, ~(0) = 0 j as a function of 8k. (see fig.3.) we studied c and a 2 (t) at different l values. let r = 0 ; and for example kt = 1, &a = 2 . eq.( 13aj is then: -; . -1 • 1 [ (l c(t) = ac (0) 1-e (smh/) j (17) and eq.( 14a): 284 i 0.7 .-·--·-·--'"··~><······-· 0.65 1------1"----+----t-----;j 1 0.6v i ~ o.ss 1------1"----+----t-----; os '----....-l--~--'----...,.1---_...; 1 fig.4 plot of conversion against the number of drawings fig.4 shows the c(t) +l relationship according llc* (0) to eq.( 17), while the [ ) for 15 seconds and then were weighed by an analytical balance. from the values of weight losses, the gravimetric fi~re k (g m"2 h'1) and the penetration p (mm m2 year· ) were calculated. the metallic surfaces were visualised by electron microscopy on a tesla b300 microscope. for the values 4%, 6% and 8% of water concentration, the catodic and anodic polarisation curves were plotted on a t acussel s8r potentiometer with input impedance of 1012 0. from the shape of the polarisation curves the kinetic parameters (fist> scor and icor) of the corrosion process were calculated. the corrosion final compounds, for every value of water concentration, were insoluble in the system. they were analysed by x-ray diffraction on a hzg 4c karl zeis yena diffractometer using co (kj radiation, by lr spectroscopy on a specord m82. the chemical composition {c, h, 0 and fe) of the final compounds was also determined. results and discussion the ph-value of the corrosion solution, at the concentration of 8% add formic is 1.42 in an anhydrous medium. when hzo is added to the solution, the phvalue decreases to 1.37 for s% h20, fact indicating that the acid dissociation takes place in the presence of the two organic solvents, and water has no major contribution to this process. the values of the weight losses, converted into gravimetric figure k and penetration figure p, are presented in table 2. 19 table 3 the values of the corrosion parameters in the system methanol10% ethyleneglycol8% formic acidwater parameters of corrosion 4 ol37 %h20 6 8 4 ol50 ol60 %hz0 %hz0 6 8 4 6 8 e,.,mv ~>con mv icrm f!a/cm 2 -520 -493 -590 -475 -467 -542 -490 -448 -527 -535 -510 -600 -490 -495 -555 -510 -480 -550 25.36 20.36 18.87 23.86 18.12 16.35 20.12 17.25 12.78 ·1 em fig.3 the ir spectrum of the corrosion compounds resulted from ol 37 in the system methanol10% ethyleneglycol8% formic acid 4% water from these values, we can conclude that the corrosion rate is high for water concentration, higher than 4%, a fact indicating that this is the minimal percent of water necessary in passive layer formation. the metallic surface, analysed by electron microscopy (xl200), presented a pitting corrosion for ol 37 and ol 50, and for ol 60, a generalised corrosive process with a ~discontinuous compound layer (figure 1). the polarisation curves were plotted for the values of 4%, 6% and 8% water concentration. the curves have the same variation aspect for all concentrations, not presenting passive domains in the studied potential field. the polarisation curves for 4% water concentration are presented in figure 2. the values of the corrosion parameters calculated from these curves are presented in table 3. the density of the corrosion current decreases with the values of water concentration of 4% and the most stable steel in all water concentrations is ol 60. the corrosion products are insoluble in the system at au water concentrations and do not present x-ray diffraction spectra. the ir-spectra are the same for all the systems, meaning that all the compounds presented the same absorption wave numbers. the ir absorption spectra of the corrosion compounds resulted from ol 37 in the medium containing 4% h20 is presented in figure3. the existence of two peaks, one at 276 cm-1 and the other at 300 cm·i, characteristic of the fe-0 covalent and coordinative bonds are observed. also the displacement of the absorption band, characteristic of the carboxyl group from the formic acid, from the value 1724 cm·1 to the value 1656 cm·l, makes us to resume a stronger bonding of the acid [12-13]. also, the presence of the absorption band of 1080 cm1, characteristic of table 4 the chemical composition of the corrosion compounds in the system methanol 10% ethyleneglycol 8% formic acid4% water % ol37 ol50 ol60 c 17.87 18.05 18.04 h 3.56 4.07 3.78 0 45.06 44.76 44.83 fe 33.51 33.12 33.35 fig.4 the structure of the corrosion compounds fe(iii) polyformiate the group ho"from the alcohols and the peak characteristic of water are displaced at 3340 cm· 1 respectively 2900 cm-1• the rest of the peaks correspond to the vibration and rotation movements of the bonds c-c, c-h, a.s.o. [12-13]. chemical analysis of corrosion the compounds indicates that all the products present the same composition. in table 4, the chemical composition of the compounds obtained in the system with 4% h20 is presented. this chemical composition is very much alike with that of a compound in which the ratio fe/formiat radical is 1/2. from these data, for the corrosion of iron in the system methanol-ethylene glycol-formic acid-water, the following mechanism is proposed: *02 + h20 + 2e· --+ 2hofe2+ + 2ho. --+ fe(oh)2 nfe(ohh + 2nhcooh --+ [fe(hcoo)i!n + nh20 the structure of the polymeric compound of fe(il) with formic acid is presented in figure 4. the hexacoordination of iron is accomplished by water molecules and .-hogroups as bridge linkage. this structure is justified also by the tittle excess in the percent of oxygen and hydrogen from the elemental analysis of the corrosion compound. 20 • • • conclusion all three types of steel present a value for water concentration for which the corrosion rate has a maximum, and that is 4% h 20. at the value 4% h20, the molar ratio water/formic acid is 1.28/1, smaller than the ratio 4/1 that is specified in literature [14-15] as being the minimum ratio in which water can participate in oxyhydroxylic layer formation. this fact can be justified by the dissociation of the formic acid and through the organic solvent molecules, the molecules of water are ''free" to participate in the passive layer formation. the steel with the highest stability is ol 60, which has the highest carbon content. this high value of the stability may be connected to the existence in the steel's structure of a solid solution of ironcarbon, formed where the carbon concentration is over 0.3% [16]. the corrosion mechanism is a complex one, involving dissolved oxygen. the corrosion compounds being insoluble in the organic system are submitted to a polymerisation process. these products do not form crystalline compounds. references 1. lorentz w. j. and geona d.: electrochim. acta, 1975, 20, 273-278. 2. lorentz w. j. and helman h.: corrosion, 1979, 27, 101-108. 3. larbeer p. and lorentz w. j.: corrosion science, 1980, 20, 405-410. 4. larbeer p. and lorentz w. j.: electrochim. acta, 1980, 25, 375-401. 5. sato n. and cohen m.: j. electrochim. soc. 1964, 111,512-518. 6. heusler k. e.: z. electrochem, 1958, 62, 582-587. 7. kama to 0, saito h, shibata t : corrosion science, 1984, 24, 807-814. 8. lorentz w. j. and heusler k. e.: anodic dissolution of iron group metals in corrosion mechanism, ed. f. mansfeld, dekker, new york, 1987. 9. sutiman d.: bul. i.p.i., 1999, xlv 3-4, 19-23. 10. sutiman d. and cretescu i.: rev. de chim. 1999, 10, 50, 766-770. 11. sutiman d., georgescu 0. and cioroianu t.: hung. j. ind. chern., 1999,27,107-110. 12. avram m.: infrared spectroscopy, ed. tehnica, bucuresti, 1960 (in romanian). 13. the sadler of handbook of infrared spectra, sadler hayden, london, 1978. 14. stypula b.: corr. resit. alloy, 1965, l, 252-257. 15. banas j. and stypula b: metalurgy and foundry engineering, 1995, 6, 2, 112-119. 16. cartis g.: thermal trataments, ed. facia, timisoara, 1982 (in romanian). page 19 page 20 page 21 page 22 hungarian journal of industrial chemistry veszprem vol. 29. pp. 3538 (2001) a novel method for isolation of lipoxygenase-1 isoenzyme from soybean meal a. sz. nemeth, b. szaja.ni1, j. sz. marcz¥ and m. l. s1mon2 (research institute of chemical and process engineering, university of kaposvar p.o. box 125, veszprem, h-8201, hungary 1 ' sato ltd, szel street 4, budapest, hungary 2 department of biochemistry, attila j6zsefuniversity, p.o. box 533, szeged, h-6701, hungary) received: january 30, 2001 a novel method has been elaborated for the isolation of lipoxygenase-1 (lox-1) isoenzyme from soybean meal. the essential step of the procedure is the selective heat treatment of soybean meal extract at 70 "c in a medium of 0.05 ionic strength and ph 5.2 for 5 min. the heat treated extract practically free from lipoxygenase-2 (lox-2) and lipoxygenase-3 (lox-3) isoenzymes was pnrified by deae-cellulose ion exchange chromatography. the resulted lipoxygenase-1 isoenzyme was homogeneous as judged by page. a tenfold purification and a relatively high yield (39 %) was achieved. keywords: soybean; glycine max; lipoxygenase isoenzymes; selective heat treatment introduction soybean seeds contain an iron"containing dioxygenase namely lipoxygenase (linoleate : oxygen oxidoreductase, ec 1.13.11.12) catalyzing the hydroperoxidation of polyunsaturated fatty acids and esters containing a cis,cis-1,4-penta:diene system [1]. the hydroperoxides formed are cleaved by hydroperoxide lyase, resulting in c6 aldehydes in plants. there are at least three lipoxygenase isoenzymes in soybean seeds, named lipoxygenase-1 (lox-1), lipoxygenase-2 (lox-2) and lipoxygenase-3 (lox-3), which exhibit differences in substrate and product specificity, optimum ph for catalytic activity, isoelectric point and thermal stability [2-4]. the levels of the lipoxygenase isoenzymes are influenced to different degrees by both cultivar and climatic effects. the weather conditions can exert a greater effect than the cultivar effect £5]. with linoleic acid as substrate, linoleic acid 13hydroperoxide and 9-hydroperoxide can be formed in the reaction catalyzed by iipoxygenase. the ratio of isomers derived from the reaction depends on the isoenzyme(s) present and the reaction conditions {6]. soybean lipoxygenase can be used to produce flavour compounds, but in some cases the presence of lox-2 and lox-3 isoenzymes is undesirable because of side reactions. a large number of methods have been suggested for separation of the isoenzymes, e.g. ion exchange chromatography £7-9], chromatofocusing [loll], electrophoresis [12] and isoelectric focusing [7, 10]. we have found [13] that the selective heat treatment of the soybean meal extract makes the elimination of lox-2 and lox-3 isoenzymes possible. based on this finding, a novel method for the isolation of lox-1 isoenzyme has been elaborated. materials and methods reagents linoleic acid, tween 20, deae-cellulose and soybean lipoxygenase (type i-b) were purchased from sigma chemical co. (st. louis, mo, usa). acrylamide, n,n,n' ,n' -tetramethylethylenediamine and coomassie brilliant blue r-250 were supplied by reanal factory of laboratory chemicals (budapest, hungary) and n,n' • methylene-bis(acrylamide) was obtained from fluka a.g. (buchs, switzerland}. the other chemicals were of analytical grade. the sodium linoleate substrate was prepared according to axelrod et al. f14}. 36 apparatus a biochrom 4060 spectrophotometer was used for enzyme activity measurements and a 2117 multiphor ii electrophoresis system for electrophoresis (pharmacia lkb, uppsala, sweden). lipoxygenase assay the method of axelrod et al. [14] as modified by marczy et al. [5] was applied, except that 0.2 m sodium borate buffer was used for lox-1. one unit of activity was defined as the amount of enzyme producing 1 j.llllol of linoleic acid hydroperoxide per min at room temperature at ph 9.0 (lox-1) or 6.8 (lox-2 and lox-3). protein content protein contents were determined according to lowry et al. [15], with bovine serum albumin as standard or by absorbance at 280 nm. electrophoresis polyacrylamide gel electrophoresis (page) was performed in a continuous buffer system of 0.1 m sodium phosphate buffer (ph 7 .5). the total monomer concentration [acrylarnide + n,n'-methylenebis(acrylarnide)] in the gel was 5 %, the concentration of n,n' -methylene-bis(acrylarnide) being 2.7% of the total monomer concentration. samples of 15 j.tg protein in 0.01 m sodium phosphate buffer (ph 7 .5) were applied onto a gel slab {12 em x 4 em). electrophoresis was carried out for 2 h at 200 ma at 5 °c. the gel was stained for protein with 0.25 % of coomassie brilliant blue r-250 dissolved in methanol:water:acetic acid (23:23:4). experimental preparation of crude extracts soybean meal was prepared :&om a chandor cultivar harvested in 1993. kindly provided by the cereal research institute (szeged, hungary). ten grams of soybean meal was defatted with 20 cm3 portions of petroleum ether at 0 <>c. defatted meal was extracted with 5 volumes (v/w) of 0.05 m acetic acid by mechanical stirring for 1 h at 4 oc. the suspension was then filtered through cheesecloth and centrifuged at 10,000 g for 20 min. the resulting supernatant was used in the further experiments. selective heat treatment the crude extract was heated to 70 oc in a water bath under continuous stirring and was kept at this temperature for 5 min [13]. immediately after heat treatment the suspension was cooled to room temperature in an ice bath and then was centrifuged at 10,000 g for 20 min at room temperature. ion exchange chromatography anion exchange chromatography of heat treated soybean meal extract was performed according to axelrod et al. [14] and weber et al. [16] with modifications on deae-cellulose. 14.5 cm3 heat treated soybean meal extract was dialyzed against 3 x 500 cm3 0.01 m sodium phosphate buffer (ph 7.0). the enzymatically inactive precipitate was discarded and the supernatant was loaded on a deae-cellulose column (1.85 x 21 em) equilibrated with 0.01 m sodium phosphate buffer (ph 7.0). elution (22.5 cm31h) was performed with a linear gradient formed from equal volumes of 0.01 m and 0.22 m sodium phosphate buffer (ph 7 .0). fractions of 6 cm3 were collected and the protein contents were determined by absorbance at 280 nm. activities were assayed at ph 9.0 (lox-1) and 6.8 (lox-2 + lox-3). the enzymatically active fractions of lox-1 were pooled and concentrated by ultrafiltration. results and discussion· data from a typical lox-1 isoenzyme isolation and purification experiment are summarized in table i. the data show that the selective heat treatment of soybean meal extract at 70 oc is a crucial step in the isolation of lox-1. though the purification of lox-i was less than twofold in respect of specific activity as an effect of heat treatment, the amounts of lox-2 and lox-3 isoenzymes were markedly reduced by this simple method. the anion exchange chromatography makes the removal of further protein impurities possible (fig.l}. tenfold purification and a relatively high yield (39 %) was achieved by the two purification steps. the purified lox-l isoenzyme was electrophoretically homogeneus as judged by page (fig.2). 37 table 1 isolation oflipoxygenase-1 isoenzyme from soybean meal a volume totallox-1 lox-1 total lox-1 purification step (cm3) activity yield protein specific activity activity purification (units) (%) (mg) (units/mg erotein) ratio* extraction 17.0 582.6 100.0 317.6 1.83 5.4 1.0 heat treatment 14.5 392.5 67.4 130.8 3.00 20.1 1.6 deae-cellulose 30.0 227.3 39.0 chromatography 11.1 20.48 20.8 11.2 •rn the experiment, 5 g soybean meal was processed *lox-1 activity i (lox-2+lox-3 activity) """ 0.8 12 0.8 ] /"".. 1.0 ~ s -w ! 0.6 _§ l:l "" 0.6 0.8 0 :t:: ~:g <7 § 8 .::: ~ '-' 8 !:l v ~ i !':; 0.6 0 0.4 x 6 0.4 0 !:l ~ 0 .. 0 ....:) ~ 0 ...... .... 0 0.4 ~ ....:) ~ 4 g-i ...... 0.2 ::; 0 ·:~ 0.2 ..0 0 0 ~1 )· * ,;.. .. """' 1 to ' ; ' t lj ,.,.....,. * '-' -+ '-" fig. 2 polyacrylamide gel electrophoretic pattern of purified lox-1 isoenzyme. electrophoresis was carried out in 5 % polyacrylamide gel slab (12 x 4 em) for 2 h with 200 rna at 5 oc. samples 15 vg protein each. a: sigma type i-b iipoxygenase; b: lipoxygenase-lisoenzyme purified with the new method conclusion we have found earlier {17, 181 that the selective heat treatment of aqueous extracts of animal organs as kidney and pancreas could be advantageous in enzyme isolations. it was supposed that this simple method would be useful in the isolation of plant enzymes too. the difference in the heat stability of soybean lipoxygenase isoenzymes {19] and the practical demand on lox-1 isoenzyme make the application of heat treatment obvious. the heat treated soybean mea! extract is practically free from lox-2 and lox-3 isoenzymes. it can be used immediately for the production of linoleic acid 13-hydroperoxide, or it can be stored frozen without deterioration for a long time. in comparison with the known methods [7-12] the proposed novel procedure is more advantageous because it does not require special equipments and chemicals. therefore it is not expensive and is suitable for industrial purposes promoting the biocatalytic production of flavour compounds. references 1. hayashi 0., nozaki m. and abboti m. t.: the enzymes, 3rd ed.. ed.: boyer p.d., acad. press, newyork,pp.ll9-l89,1975 2. christopher. j. p., plstorius e. i{. and axelrod b.: biochim. biophys. acta, 1970, 198, 12-19 3. grosch w. and laskawy g.: j. agtic. food chern. 1975,23,791-794 4. zhuang h., hu..debrand d. f., andersen r. a and hamilton-kemp t. it: j. agric. food chern. 1991.39,1357·1364 5. marcz'l:' simon m. l., m6zs!k l. and szajanib.: j. agric. foodchem.j995,43, 313·315 3-8 6. shibata d., s1eczko j., dixon j. e., hermondson m., y azdanparast r. and axelrod b.: j. bioi. chem.1987, 262, 10080-10085 7. verhue w. m. and francke a.: biochim. biophys. acta 1972, 284,43-53 8. hausknecht e. c. and funk m. 0.: phytochemistry 1984,23, 1535-1539 9. shkarina t. n., kuhn h. and schewe t.: lipids 1992,27, 690-693 10. funk m. 0., whitney m. a., hauscknecht e. c. and o'brien e. m.: anal biochem.1985 ,146, 246251 11. fbiters m. c., veldink g. a. and vljegenthart j. f. g.: biochim. biophys. acta 1987,870,367-371 12. christopher j.p., pistorius e. k. and axelrod b.: biochim. biophys. acta, 1972,284, 54-621 13. nemeth a. sz., szaiani b., marczy j.sz. and simon m. l.: biotech. tech., 1998, 12(5), 389-392 14. axelrod b., cheesbrough t. m. and laakso s.: methods enzymol. 1981, 71, 441-451 15. lowry 0. h., rosebrough n. r., parr a. l. and randallr. j.: j. biol. chem.1951, 193,265-275 16. weber f., laskawy g. and grosch w.: z. lebensm. unters.-forsch.1974, 155, 142-150 17. szaiani b.: acta biochim. biophys. acad. sci. hung. 1980, 15, 287-293 18. dala e., kiss a., sum p. and szaiani b.: acta biochim.biophys. hung. 1986, 21, 327-333 19. 19. matoba t., hidaka h., narita h., kitamura k., kaizuma n. and kito m.: j. agric. food chern. 1985, 33, 852-855 page 36 page 37 page 38 page 39 hungarian journal of industrial chemistry veszprem vol. 29. pp. 61-65 (2001) toward global optimization of a difficult optimal control problem r.luus (department of chemical engineering, university of toronto, toronto, on m5s 3e5, canada) received: july 5, 2001 the optimal control problem introduced in 1980 by yeo [1] exhibits numerous local optima which have not been reported in the literature. when 50 time stages of equal length are used, there are only two local optima. when so time stages are taken, there are over 50 local optima. it was found that there are at least 15 local optima which are better than the one reported in the literature as the global optimum. with the use of 100 time stages of equal length, the number of local optima is increased substantially to over 150. with 100 time stages of equal length, the global optimum solution found with iterative dynamic progranuning (idp) gives i= 0.119044. this value of the performance index is 0.2% better than the value that had been accepted in the literature as the global optimum. idp provides a good means of getting the global optimum solution for this problem, but numerous runs are necessary due to the presence of a large number of local optima. keywords: optimal singular control, iterative dynamic progranuning, idp, global optimization, sensitivity introduction yeo's optimal control problem establishing the optimal control policy for nonlinear control problems is usually quite difficult, especially when the system exhibits low sensitivity of the performance index to a. change in the control policy. when the number of time stages is small, then the optimal control policy can usually be readily obtained either by iterative dynamic programming [2,3) or with the lj optimization procedure [4}. for many problems a small number of flexible stage lengths, where the stage lengths are incorporated into optimization, can be used to give very accurate results by enabling the optimal switching times to be accurately determined. however, there are some problems where the structure of the optimal control policy changes quite substantially when the number of time stages is increased. this is especially true with the system introduced by yeo [1]. the goal of this paper is to illustrate the interesting features presented by this system and to show how the optimal control policy can be established in a reliable way. one of the earliest examples tried with iterative dynamic programming [5,6], is the optimal control problem used by yeo fl] to illustrate the use of quasilinearization to solve nonlinear singular control problems. the system is described by dxl -=xz dt dx _3_::::u dt dxs =1 dt {2) {3) (5) 62 with x(o) = [ 0 -1 .j5 0 0 f, and final time t.f = 1. the constraints on the control are -4::::; u $10 (6) it is required to find the control policy that minimizes the performance index (7) this system was used for global optimization by rosen and luus [7] with the introduction of a line search into sequential quadratic programming (sqp). with the use of p = 80 time stages of equal length they obtained i= 0.119258, and concluded by extrapolation that the minimum value for the performance index for the continuous case (infinite number of stages) should be i= 0.11923. this system was recently investigated by esposito and floudas [8], with p = 10 time stages, who found two local optima as was found by luus [3]. the differential equations are well behaved, but to get six figure accuracy for the performance index, it was necessary to use double precision and to use a tight local error tolerance of 10·8 in the subroutine dverk [9]. the optimal value of the performance index depends very much on the number of time stages that are used. the goal here is to illustrate the nature of the optimal solution when the. number of time stages of equal length is increased. table i affect of the number of time stages p on the value of the performance index after 50 passes, each consisting of 20 iterations number of time staaes p 10 20 40 50 60 80 100 performance index i 0.120114 0.119438 0.119292 0.119277 0.119272 0.119134 0.119085 results and discussion cpu time seconds 5.9 9.0 28.9 47.1 61.1 110.5 180.7 in using iterative dynamic programming (idp), we used the multi~pass approach as outlined by luus [3, p. 106], using 20 iterations per pass. the region contraction factor was chosen as r= 0.95 and the region restoration factor 17 = 0.85. a single grid point was used at each time stage, and 3 candidates for control, chosen as o. ±r were used at each grid point. the initial value for control was chosen as 7.0 and the initial region size r = 4 was used. the computational results carried out on a pentiumiy400 personal computer in double precision after 50 passes for different number of time stages p are given in table i. it is noted that idp yields the optimum substantially faster than the global optimization procedure of esposito and floudas [8]. even with p = 100 the computation time is almost negligible. the use of p = 60 provides a very small refinement to the control policy obtained with p = 40 and p == 50 for which the optimal control policy is shown in fig.]. e "e 0 (.) 0.0 0.4 time,t 0.8 fig.l optimal control policy for p = 50 giving i = 0.119277 0.0 m time,t 0.8 fig.2 optimal control policy for p = 80, giving i ::: 0.1191325 ~ "c .e 8 c: til e .g c]) a. 0.11912 l---1---l..----jl--'--l----jl-_[__l__j 72 76 position of dip, stage number fig.3 local optima asa function of the location of the dip in the control pro.le for p = 80 e 'e 8 0.0 0.4 time, t 0.8 fig.4 control pro.le for a local optimum for p = 80, giving i= 0.1191603 since there is a significant improvement in the performance index for p = 80 stages, the run was repeated, using 40 iterations per pass and allowing 100 passes. the performance index was marginally decreased to i= 0.119133. as is shown in fig.2, the use of p = 80 yields a totally different optimal control policy. by using a greater number of time stages, the ~hotter stage lengths allow sudden changes to be made ill the control policy. these abrupt changes in the control policy here improve the performance index quite significantly. ::::. ~ 0 (.) -4 0.4 time,t 63 0.8 fig.5 control pro.le for local optimum for p = 100 giving i= 0.119250 e 'e 0 (.) -4 0.4 time,t 0.8 fig.6 control pro.le for local optimum for p = 100, giving i = 0.119092 there is a dip in the control policy to -4 at stage 59. and in the time interval (0.8, 1.0) there is a dip to -4 at stage 75, with a significant improvement in the resulting performance index over the value of 0.119258, where there are no dips to the minimum value of control as is shown in fig.3, the position of the latter dip affects the performance index. with p = 80 we can also get local optima with two dips in the time interval (0.8, 1.0}, as is shown in fig a, where the dips at stages 70 and 76 give a performance index i= 0.119160. 64 0.1190580 85 90 95 100 position of dip, stage number fig.7 local optima asa function of the location of a single dip in the control pro.le with p = 100 0.0 0.4 time,t 0.8 fig.8 control pro.le for a local optimum with p = 100. giving l = 0.119045 the gk1bal optimum for p = 80 with 1 = o.h91325 occurs with a single dip at stage 75. with the use of p = 100 time stages of equal length the number of local optima is quite large. the local optimum control policy in which there are no dips, yielding i= 0.119250 is shown in fig.s. the best local optimum with a single dip at stage 94, giving i = o.hooij2 is shown in fig.6. as is shown in fig.7. the position of the single dip affects the value of the performance index. as cae be seen 14 local optima with a sin,e!je dip give a better value for the performam:e index than the local tlptimum wiu!out any dips. 0.0 0.4 time, t fig.9 optimal control policy for p = 100 giving i= 0.1190437 0.0 0.4 time,t fig.jo control policy for p = 150, giving i = o.h8977 what makes this problem challenging is the existence of double dips in the time interval (0.8, 1.0), many of which give better values for the performance index. in fact, the best local optimum with dips at stages 88 and 95, giving 1 = 0.119045, is shown in fig.8. further improvement to 1 = 0.1190437 is obtained with the control policy shown in fig.9. where there are additional 3 dips in the time interval (0.2, 0.8). there are over 150 local optima for this system with p = 100 stages of equal length. here. with ldp only the most likely candidates for global optimum were scanned to yield the global optimum of l = 0.119044. 0.0 0.4 time, t 0.8 fig.ll state trajectories corresponding to the 65 this example illustrates the difficulty in determining with absolute certainty the global optimum of some nonlinear systems. idp proved to be an effective tool in the establishment of local optima. through the interpretation of the results obtained for the local optima, the global optimum for 100 time stages was readily obtained. however, the complexity of the problem is realized, so the results obtained here provide only a step toward obtaining the global optimum to the problem for which we could conceptually use an infinite number of time stages. acknowledgement financial support from the natural sciences and engineering research council of canada is gratefully appreciated. symbols control policy in x figure 10 i performance index to be minimized number of time stages of equal length the question of what happens when the number of time stages is further increased arises. it is obvious that the number of local optima will increase. in fact, a run was performed with p = 150, yielding i = 0.118977$, with the control policy shown in fig.jo. this result was obtained with the conditions used for table 1, so this result is expected to be close to the global optimum. it is interesting to note that now there are three dips in the time interval (0.8, 1.0). the corresponding trajectories for the first three state variables is given in fig.ll. there is no doubt that further noticeable improvement of the performance index can be obtained by using even a larger number of time stages, such as 200. alternatively, one could use time stages of varying length to refine the results obtained here. research in this area is continuing. conclusions by increasing the number of time stages to solve the optimal control problem presented by yeo [1], some very interesting features about this system were noted. with the use of 80 time stages numerous local optima, that had been observed for the frrst time, were found. a large number of these local optima were better than the best local optimum that was obtained by sqp. the use of 100 time stages yielded even a larger number of local optima. by scanning numerous local optima, the global optimum with 100 stages of equal length was determined as i= 0.119044. this value is 0.2% better than the best value reported in the literature. with the use of 150 stages, the performance index was further improved to 0.118977, showing that further increases in the number of time stages, or the use of stages of varying length should yield even better values for the performance index.. p r region size t time f_t final time u control variable x1 ith state variable x state vector greek letters y region contraction factor by which the 1 egion size is reduced after every iteration 1j restoration factor by which the region size is restored after every pass abbreviations idp iterative dynamic programming sqp sequential quadratic programming references 1. yeob.p.: int. j. control, 1980, 32, 723~730. 2. luus r.: control and intelligent systems, 1998, 26, 1~8. 3. luus r.: iterative dynamic programming, chapman & hall crc, london, uk, 2000. 4. luus r. and hennessy d.: ind. eng. chern. res., 1999,38, 1948·1955. 5. luus r.: hung. j, ind. chem., 1989, 17, 523~543. 6. luus r.: int. j. control1990, 19, 995~ 1013. 7. rosen 0. and luus r.: j. optimization theory & applk., 1992, 73, 547--562. 8. esposito w.r. 'lnd floudas c.a.: j. global optimization, 2()()(), 17.97-126. 9. hull t.e., enright w .d. and jackson k.r.: u5er guide to dverk •• a subroutine for solving nonstiff ode's, report 100. 1976, department of computer science, university tjff(.'fnnto, canada page 60 page 61 page 62 page 63 page 64 hungarian journal of industry and chemistry vol. 50 pp. 1–5 (2022) hjic.mk.uni-pannon.hu doi: 10.33927/hjic-2022-01 antioxidant effect of humic substances from hungarian leonardite attila csicsor*1,2 and etelka tombácz1,3 1doctoral school of environmental sciences, university of szeged, rerrich béla tér 1, szeged, 6720, hungary 2hymato products kft., kossuth lajos u. 33, szentkirályszabadja, 8225, hungary 3soós ernő water technology research and development center, university of pannonia, zrínyi miklós u. 18, nagykanizsa, 8800, hungary humic substances are natural substances that are continuously formed from the decay of plant residues. these materials have a very diverse range of properties, making them versatile. according to many new studies, these humic substances also exhibit antioxidant propensities. the aim of this paper was to shed light on whether humic substances really have antioxidant properties. keywords: humic substances, humic acid, fulvic acid, himatomelanic acid, antioxidant, leonardite 1. introduction the majority of the population of hungary suffers from one of the diseases associated with free radicals, e.g. diseases of civilization, obesity, cardiovascular disease, malignant neoplasms, etc. although it is well-known that our consumption of vegetables and fruit lags far behind those of the other member states of the european union, it has been scientifically proven that essential fresh vitamins and provitamins as well as minerals are fundamental not only for the smooth functioning of the body but also for their antioxidant components. furthermore, they can play a key role in terms of disease prevention. berries have become a major player in a number of related research papers as they are also essential in the prevention as well as aftercare of cancer and cardiovascular disease. what other options are available and what other substances may still have high antioxidant contents? based on the aforementioned points, it is also important and justified to study the antioxidant properties of new compounds such as humic substances and determine which of their components play a crucial role in the development of their antioxidant effect, both qualitatively and quantitatively. the aim of this research is to compare the antioxidant properties of humic acids with already wellknown natural antioxidants. firstly, different humic acid recieved: 26 august 2021; revised: 16 september 2021; accepted: 20 september 2021 *correspondence: csicsor.attila@gmail.com fractions were prepared. the expected results should help to more accurately interpret the complex behavior of humic acids, thereby helping to expand their range of applications. therefore, if the results are convincing, this could open up new fields with regard to the application of humic acids in the food and cosmetics industries [1]. 2. humic substances humic substances are a group of naturally occurring macromolecules found throughout nature, that is, in soil, air, water, carbon deposits and peat. chemically illdefined humic substances are natural organic colloids comprised of decomposition products of plant-derived biomass as the result of a process called humification. their conversion into stabilized humic material is one of the most complex and least understood biogeochemical processes of the carbon cycle. most natural humic acids (ha) are found in older peat, lignite and juvenile lignite. in hungary, the ha content of the so-called leonardite, a famous source of ha, is close to 70% in its natural form. (a near-surface deposit of leonardite is located in dudar, where this geological formation can be mined). in the humification process, dead plant matter that enters the soil is broken down by enzymes found in soil bacteria and fungi, so simple compounds like sugar and ammonia are formed from carbohydrates, fats, proteins and lignin, which on the one hand serve as a source of food for the soil microbes but on the other hand as a source in the formation of humic substances. as a result of biotic and abiotic (condensation and polymerization) processes, https://doi.org/10.33927/hjic-2022-01 mailto:csicsor.attila@gmail.com 2 csicsor and tombácz the decomposition products form high-molecular-weight humus compounds, the presence of which is characteristic of the soil [2]. on the basis of differences in solubility, humic substances can be divided into several groups such as the main fractions humin, ha, and fulvic acids (fa), as well as the alcohol soluble himatomelanic acids (hy). they are a mixture of similarly behaving, yellow-brown-black, acidic, high-molecular, natural organic substances that are operationally defined. humic substances themselves contain too many kinds of molecules that can be separated by changing the solubility conditions. even though these molecules behave similarly, each fraction has different properties and their molecular structure is not uniform [2]. the nomenclature of fractions that can be separated according to their solubility is suggested in the guidelines of the international humic substances society (ihss): • humin a black fraction of humic substances that is insoluble in water at any ph. • humic acids the fraction of humic substances that is insoluble in water under acidic conditions (ph < 2) but soluble in water at higher ph values. they can be extracted by various alkaline solutions before being precipitated by a strong acid and are dark brown or black in color. • fulvic acids the fraction of humic substances that is soluble in water at all ph values. fulvic acids are light yellow or yellowish brown in color. • himatomelanic acids the fraction of humic substances that is soluble in an alcohol. 3. antioxidants the most important physiological role of antioxidants is to neutralize the free radicals that are continuously formed in the szent-györgyi-krebs cycle and counteract the free radicals with different oxidizing forces that enter the body. an antioxidant is a substance that inhibits oxidation, more broadly speaking, retards or hinders oxidation. they are chemically reducing agents, i.e. electron donors. these materials are usually organic compounds but include metals and organometallic complexes. many different types of antioxidants exist, which usually work together and do not neutralize free radicals on their own. although our body itself is able to produce some antioxidants, sometimes it is necessary for us to ingest these substances from external sources. together, these antioxidants already form a very strong line of defence in our body, e.g. vitamin c, vitamin a, flavonoids, glutathione, resveratrol, unsaturated fatty acids, etc. [3, 4] 4. methods for measuring the antioxidant capacity the antioxidant capacity is the combined free radical scavenging effect of all antioxidant compounds in the system studied. since the need for its accurate numerical determination is growing, a number of analytical procedures and measurement systems have been developed. given that methodologies are constantly being modified and refined, nowadays, the number of applied methods exceeds one hundred. in the literature, the majority of studies use several methods to determine the antioxidant capacity [5]. the most common methods of measuring the antioxidant capacity can be divided into two main groups: electron transition-based ones (ferric iron reducing antioxidant power, total polyphenol content, copper ion reducing antioxidant capacity, trolox equivalent antioxidant capacity, and 2,2-diphenyl-1-picrylhydrazyl (dpph)) as well as those based on hydrogen atom transfer (oxygen radical absorptance capacity, total peroxyl radical-trapping potential, chemiluminescence-based methods, and photochemiluminescence measurements just to mention the most commonly used methods). although these two types of measurements determine the antioxidant capacity, the results obtained do not necessarily have to be correlated with each other, since the reducing power of a sample is not necessarily related to its ability to scavenge for the reaction of test compounds. electron transfer reactions involve colour changes, from which the antioxidant capacity can be deduced. the essence of these methods is to create a free radical as the result of a reaction. to this free radical, the antioxidant is added at various dilutions leading to a colour change, which is monitored by a spectrophotometer and then the antioxidant capacity of the test substance is calculated from the results obtained. methods involving hydrogen atom transfer are based on the kinetics of the reaction. tests measure how effective a sample is against a given free radical, namely its free radical scavenging capacity [5]. 5. antioxidant effect of humic substances the question may arise as to why humic substances would exhibit an antioxidant effect. humic acids are chemically very complex mixtures of composite molecules, that is, natural polymers formed during the varying degrees of polymerization of basic building blocks. according to their chemical structure, they are polyhydroxy carboxylic acids with quinone and semiquinone groups. in some respects, they are similar to flavonoids and phenols, in which the so-called flavone skeleton is polysubstituted by hydroxy groups. however, they also have a quinoid structure that is known to be responsible for antioxidant properties. these properties of ha have already been demonstrated in a number of scientific publications by both classical analytical methods (redox titrations) and instrumental analytical measurements (electron spin resonance) [6]. the general structure of ha and the formulae of some well-known antioxidants are presented in fig. 1 from hungarian journal of industry and chemistry antioxidant effect of humic substances 3 figure 1: model structure of ha and identical moieties of known antioxidants [2] their structure, it is clear that ha are comprised of a number of groups such as already well-known antioxidants, so will exhibit exceptionally high antioxidant capacities [7]. 6. results 6.1 measurement of total phenolic content all measurements of total phenolic content were made according to the method developed by shetty et al. [8]. during these measurements, only high-quality (a.r.) chemicals were used. the fractions of humic substances were extracted by ourselves from samples of leonardite retrieved from dudar. rather than indirectly measuring the antioxidant property of the sample, the method measures its total phenolic content, from which its antioxidant capacity can be deduced. the method consisted of diluting 1 ml of the sample in 5 ml of distilled water and 1 ml of 95% ethanol in a test tube before adding 0.5 ml of 50% folin-ciocalteu (fc) reagent (half of distilled water in half fc reagents) to each sample. after being stirred for 5 minutes, 1 ml of 5% na2co3 was added to the reaction mixture and allowed to stand for 1 hour before the absorbance values of the samples were measured at 725 nm. these absorbance values were then converted into µg gallic acid equivalents in order to compare the measured values with each other and with other data from the literature. whilst measuring the total phenolic content of the samples, a series of dilutions was made from our samples of known mass, for which the total phenolic content was measured. a series of dilutions had to be produced in order to measure an absorbance value of approximately one for each sample. this was necessary to be able to compare the samples because the three different fractions yielded one absorbance value at different concentrations once the reaction was complete. furthermore, since the samples are coloured (brownish dark brown), the measured absorbance values had to be calibrated in light of the background of the samples. therefore, the samples were diluted to the concentration present in the reaction volume. after measuring the series of samples, their total phenolic content was calculated from the calibration curve, which can be seen in fig. 2, and the results obtained are shown in table 1. the results show that different initial concentrations of the various samples are required to achieve a similar absorbance. in the humic acid samples, since the majority of the phenol groups are in this fraction, the lowest initial concentration is required to achieve an absorbance of approximately one. himatomelanic acid has an average number of phenol groups while fa is comprised of the figure 2: calibration curve of gallic acid 50 pp. 1–5 (2022) 4 csicsor and tombácz table 1: the measured and corrected absorbance of the samples and the gallic acid equivalent (gae) values sample corrected concentration [mg/ml] abs measured abs corrected gae [mg/ml] from the calibration curve gae [mg/g] hy 0.59 1.027 0.953 0.103 174.5 fa 1.18 0.999 0.948 0.102 86.8 ha 0.35 1.024 1.024 0.11 310.3 ihss fa 0.684 1 0.107 156.8 ihss ha 0.212 1 0.107 506.1 table 2: the ic50 readings of the different samples name of the sample original concentration ic50 [µg/ml] hy 1.5 mg/ml 200 fa 2 mg/ml 300 ha 3.25 mg/ml 460 fewest. the results suggest that the studied humic fractions are likely to exhibit an antioxidant effect. 6.2 dpph method measurement of the antioxidant capacity based on stable dpph radical scavenging is one of the first methods. the reaction proceeds as follows: the dark purple radical formed in the reaction mixture loses its colour when it reacts with antioxidants. this method is widely used because the radicalforming molecule dpph is commercially available, stable as well as not particularly reactive nor aggressive, which is beneficial in the reactions that take place as taking measurements is simple. on the other hand, a stable radical that is not found in the living organism is used instead of a radical formed during the normal metabolism in the cell. using this method, it is not possible to estimate how effective the sample is as an antioxidant with regard to biological radicals [8]. during the measurements, 6 ml of a dpph working solution (0.1 mg/ml) was added to one ml of the sample, vortexed and stored in the dark for 30 minutes until the colour reaction began. after 30 minutes, the absorbances were measured at 517 nm. from the measured values that had previously been corrected, the percentage of inhibition (inhibition%, inhib.%) and the value corresponding to 50% inhibition (ic50 value) were calculated. while the measurements were taken, in the same way as when the total phenolic content was measured, a series of dilutions was made from a sample of known concentration. from the measured absorbances, after background correction the percentages of inhibition of the samples were calculated and plotted as a function of the concentration in order to determine the ic50 values of the samples (fig. 3). from the data, it can be concluded that the materials prepared and tested exhibit antioxidant properties because they inhibit the decomposition of the dpph radical. figure 3: percentage of inhibition of himatomelanic acid (top), humic acid (middle), and fa (bottom) as functions of the concentration to determine the ic50 values their ic50 values (table 2) are promising and comparable with other data in the literature [9]. 7. conclusions although the results revealed that the samples made from domestic raw materials exhibit antioxidant properties, further studies are needed to gain a complete picture of the antioxidant capacity of these substances. furthermore, in the future, it is our intention to confirm the results collected so far with in-vitro and in-vivo experiments. hungarian journal of industry and chemistry antioxidant effect of humic substances 5 references [1] balogh, e.; hegedus, a.; stefanovits-banyai, e.: application of and correlation among antioxidant and antiradical assays for characterizing antioxidant capacity of berries sci. hortic., 2010, 125(3), 332– 336 doi: 10.1016/j.scienta.2010.04.015 [2] stevenson, f.j.: humus chemistry: genesis, composition, reactions (wiley & sons, new york, usa) 1994, pp. 188–210 isbn: 978-0-471-59474-1 [3] nccih: antioxidants: in depth, retrieved 20 june 2018 https://www.nccih.nih.gov/health/antioxidants-in-depth [4] škvára, p.; kezmanová, j.; mišl’anová, c.; vojs staňová, a.: analysis and identification of polyphenolic compounds in green foods using a combination of hplc-esi-it-tof-ms/ms hung. j. ind. chem., 2018 46(1), 35–38 doi: https://doi.org/10.1515/hjic-2018-0008 [5] huang, d.; ou, b.; prior, r.l.: the chemistry behind antioxidant capacity assays j. agric. food chem., 2015, 53(6), 1841–1856 doi: 10.1021/jf030723c [6] csicsor, j.; gerse, j.; tikos, a.: the biostimulant effect of different humic substance fractions on seed germination, in: humic substances in the global environment and implications on human health, n. senesi and t. m. miano (eds.) pp. 557–562; proc. 6th intern. meeting of the intern. humic substances soc., monopoli, bari, italy, sept. 20–25, 1994 [7] ponomarenko, e.v.; anderson, d.w.: importance of charred organic matter in black chernozem soils of saskatchewan can. j. soil sci., 2001, 81(3), 285– 297 doi: 10.4141/s00-075 [8] vattem, d.a.; shetty, k.: ellagic acid production and phenolic antioxidant activity in cranberry pomace (vaccinium macrocarpon) mediated by lentinus edodes using a solid-state system process biochem., 2003, 39(3), 367–379 doi: 10.1016/s00329592(03)00089-x [9] ruiz-torralba, a.; guerra-hernández, e.j.; garcíavillanova, b.: antioxidant capacity, polyphenol content and contribution to dietary intake of 52 fruits sold in spain. cyta journal of food, 2018, 16(1), 1131-1138 doi: 10.1080/19476337.2018.1517828 50 pp. 1–5 (2022) https://doi.org/10.1016/j.scienta.2010.04.015 https://www.nccih.nih.gov/health/antioxidants-in-depth https://doi.org/https://doi.org/10.1515/hjic-2018-0008 https://doi.org/https://doi.org/10.1515/hjic-2018-0008 https://doi.org/10.1021/jf030723c https://doi.org/10.1021/jf030723c https://doi.org/10.4141/s00-075 https://doi.org/10.1016/s0032-9592(03)00089-x https://doi.org/10.1016/s0032-9592(03)00089-x https://doi.org/10.1080/19476337.2018.1517828 introduction humic substances antioxidants methods for measuring the antioxidant capacity antioxidant effect of humic substances results measurement of total phenolic content dpph method conclusions hungarian journal of industry and chemistry vol. 49(2) pp. 65–70 (2021) hjic.mk.uni-pannon.hu doi: 10.33927/hjic-2021-24 using finite element analysis in the 3d printing of metals hussein alzyod*1 and péter ficzere1 1department of railway vehicles and vehicle system analysis, budapest university of technology and economics, műegyetem rkp. 3, budapest, 1111, hungary over the last few years, additive manufacturing, or as it is sometimes known, 3d printing, has become a significant research field for researchers worldwide. the necessity to increase the strength of materials and minimize their weight in the automotive and aviation industries has urged engineers and scientists to conduct more investigations and identify manufacturing methods to replace the classical ones. additive manufacturing involves building a geometry layer by layer from a wide range of materials, which helps to utilize materials efficiently while minimizing the amount of waste produced as well as build complex, large geometries and light-weight components. furthermore, it minimizes fabrication and processing times. in this paper, three different alloys were printed (tial6v4, alsi10mg and 316l) using msc simufact software to investigate the effect of changing machines on the effective stress and surface deviation. furthermore, thermal analysis as well as mechanical, thermal and thermomechanical calibrations were carried out to determine a parameter set consisting of the laser power, inherent strains, fraction of exposure energy and volumetric expansion factor. keywords: additive manufacturing, 3d printing, simufact additive. 1. introduction over the last three decades, research into printing technology/additive manufacturing (am) has progressed from fast prototyping to industry 4.0 [1, 2]. the third industrial revolution took place in 1987 with the commercial exploitation of the first stereolithography (sla) machine, which was developed by 3d systems and named sla-1. after many editions of sla machines in the early 90s, new am principle technologies were launched, namely solid ground curing, fused deposition modeling and laminated object manufacturing [3, 4]. this development was followed by many years of continuous improvement in am technology, from resin to metal powders and from non-functional molding applications to the fabrication of medical implants [5]. am methods have been characterized in the literature based on a variety of parameters, including direct or indirect process technology, the state of the raw materials, and the materials used. the state of the raw materials is the most commonly used as illustrated in fig. 1. am is a revolutionary manufacturing technique proliferating as major industries transition from conventional to advanced production methods. threedimensional printing is used in am to convert engineering design files into fully functional and durable structures. after the components in the first layer have been bonded together using glue or heat, the second layer *correspondence: hussein.alzyod@edu.bme.hu is constructed and the bonding procedure repeated, enabling the construction of previously unimaginable geometries. am creates samples of products quickly, which is crucial since it reduces the traditional trial-and-error process, allowing novel ideas to enter the market more quickly. it may also be used to rapidly create bespoke metal items to replace old or broken industrial parts. 1.1 part-building technology metal am can be classified into four main categories as shown in fig. 2: directed energy deposition, powder bed fusion, sheet lamination and binder jetting. table 1 shows the thickness of layers in each process [6]. 1.2 fields of application metal 3d printing is the fastest-growing sector. metal am is increasingly being utilized to create final products table 1: the thickness of layers in the metal am process [6]. processing the thickness of layers (mm) directed energy deposition 0.089 − 0.203 powder bed fusion 0.08 − 0.15 sheet lamination 0.1 − 0.19 binder jetting 0.089 − 0.203 https://doi.org/10.33927/hjic-2021-24 mailto:hussein.alzyod@edu.bme.hu 66 alzyod and ficzere figure 1: types of additive manufacturing processes based on the state of raw materials [1]. figure 2: classification of metal additive manufacturing technologies [5]. in many sectors like in the aviation, automotive, medical implants, material handling robotics as well as lifestyle and sport industries. the aviation industry companies in this sector are among the most well-known early investors of am technology since they were among the first to upgrade from modest research projects to large-scale production runs. norsk titanium, airbus, alcoa and ge amongst others have begun serial production or are planning to invest in significant manufacturing operations within am [7]. the automotive industry for many years, the automotive industry has utilized 3d printing to rapidly prototype single test components or entire visualization models. leading companies are experimenting with am. small series or customized production may be the first sectors to adopt this technology [7]. 1.3 modeling approach in simufact many commercial companies, e.g., autodesk, general electric (via geonx), msc software, ansys, etc., provide am simulation, which enables designers to understand where the critical zones in a particular geometry are located and produce approximate first findings within a few hours. with this information, designers can iteratively develop their concept by taking all aspects into consideration using simulation software rather than wasting time and money on the actual manufacture of a possibly unsuitable design. simufact is a piece of software hungarian journal of industry and chemistry brownian dynamics simulation of chain formation 67 figure 3: additive manufacturing workflow [own edition]. that uses finite element analysis to simulate metal am by focusing on powder bed fusion processes. the difficulties with metal powder bed fusion am can be addressed and analyzed by simulations [8]. as is shown in fig. 3, the am workflow to date has been one of the best traditional manufacturing technologies. the revision of am goes through those processes to select the best predicted design parameters, build orientation, removal cutting, supporting structure, etc. [9]. three-dimensional objects are constructed from several thin layers, each of which with a pattern defined by a computer-aided design (cad) that is exposed using a scanning laser for stereolithography or a pattern in the absence of a mask. manufacturing information is obtained from cad or other design software before being sliced into individual layers using the standard tessellation language [10]. sometimes, the building fails due to many reasons, e.g. supporting structures or a high degree of distortion. once the best predicted design parameters have been chosen, construction commences and an inspection is made. if any defects are observed, the component is subjected to a process of trial and error to achieve the desired part. alternatively, as illustrated in fig. 4, a virtual model can be built to minimize the number of required physical builds figure 4: am simulation workflow [own edition]. in the am simulation workflow. in this virtual model, different approaches can be applied, e.g. various parameters, changing supporting structures, building orientation, and keeping treatment cutting and support removal to achieve the desired product free from any adverse effects of the process the first time it is constructed. since a feedback loop is included in the simulation, the results of the process simulation can impact topology optimization. the significant benefits of this are that the formation of some of these parts may cost thousands of dollars and the trialand-error approach is quite undesirable. 2. design methodology and simulation given that am depends on many factors to determine the required parameter set to produce the desired part, e.g., laser power, fraction of exposure energy and volumetric expansion factor, experimental tests should be conducted to measure these parameters. 2.1 sample selection as is shown in fig. 5 of the simulation, an aircraft component was chosen because of the complexity of its geometry and the importance of am in the aerospace industry, 49(2) pp. 65–70 (2021) 68 alzyod and ficzere figure 5: an aircraft component [9]. e.g., in increasing asset uptime, reducing costs, manufacturing lighter parts, enhancing durability and improving customer satisfaction. 2.2 material selection three alloys (tial6v4, alsi10mg and 316l) were chosen with different properties to construct the sample for the purpose of investigating the effect of each material on the surface (shape) deviation with respect to the original sample. 2.3 machine selection different types of 3d printing machines are available using this piece of software, e.g. two kinds of electro optical systems (eos) have been used, namely m280 and m400, the size and power of both differ. 2.4 thermal analysis the laser power has to be adjusted, moreover, the printing temperature has to be greater than the melting point of each alloy but less than its boiling point. 2.5 mechanical calibration the purely mechanical, macro-scale analysis approach of simufact additive requires the input of the so-called inherent strain values before the simulation. as is shown in fig. 6, the inherent strain values are produced during the construction process as a result of plastic strain, thermal strain and phase transformations as shown in �inherent = �th + �pl + �ph. (1) moreover, they depend on material properties and process parameters. these strains can be calibrated from experiments, e.g. by measuring the degree of distortion after cutting a printed cantilever beam and running simulations to match the experimental distortion values that should be less than 3%, or estimated based on the process. this paper used the empirical calibration by printing two cantilevers on the xand y -axis as is shown in fig. 7 that are 72 mm figure 6: inherent strain [9]. figure 7: printing two cantilevers [9]. in length, 12 mm in width, and 9 mm high. afterwards, the cantilevers are cut in the middle of the teeth, so a particular part of them remains on the base plate. finally, the maximum observed distortion in the z-axis was measured as shown in fig. 8. three values of strains are inputted into the build stage dialogue (�xx,�yy,�zz). these values are dimensionless and defined in three directions, x, y , and z. default values will be shown in the dialogue [9]: �xx = −008, �yy = −0.003, �zz = −0.030. (2) the mechanical calibration was carried out for the three alloys, namely tial6v4, alsi10mg and 316l, to calculate the inherent strains. 2.6 thermomechanical calibration the thermomechanical analysis requires that the fraction of the exposure energy and the volumetric expansion factor have been calculated before the simulation is started. figure 8: measuring the distortion in the z-axis [9]. hungarian journal of industry and chemistry brownian dynamics simulation of chain formation 69 figure 9: thermal calibration since these values depend on material properties and machine parameters, they need to be carefully calibrated, which should be done experimentally. 2.7 fractions of exposure energy the fraction of exposure energy is a term that is responsible for the peak temperature of the process. during the exposure time, the energy predominantly (but not only) melts the powder. the rest of the powder reheats the solid material after the exposure time. for the purpose of thermomechanical calibration, it is recommended to carry out the thermal calibration first to determine the optimized exposure energy and then conduct the thermomechanical calibration with the optimized fraction of the exposure energy to calibrate the volumetric expansion factor. the calibration was carried out using a specimen of a cantilever. the peak temperature has to be measured during the process to determine the maximum degree of distortion in the z-axis after cutting, as is shown in fig. 8. the measured point was chosen on the upper surface because the temperature is the highest there and the average will be recorded between the two surfaces. thermal calibration of the three alloys, that is, tial6v4, alsi10mg and 316l, was carried out to determine the fraction of exposure energy. the thermal calibration of tial6v4 is shown in fig. 9. 2.8 volumetric expansion factor the volumetric expansion factor accounts for the thermal expansion and shrinkage effects caused by the modelling approach, that is, by combining multiple powder layers and repeated heating of the heat source in one element layer. the volumetric expansion factor can either be set to isotropic in all directions or anisotropic for each individual order and scales the thermal expansion factor of the material. volumetric expansion is the change in the volume of a given mass with temperature. the thermal expansion factor is a measure of the change in volume in three dimensions as is shown in eq. ??: v = 3αvo∆t (3) this equation is usually written as: v = βvo∆t, (4) where β = 3α denotes the volumetric expansion factor, vo represents the original volume and ∆t stands for the change in temperature. 3. results different machines and materials were used to print an aircraft component. as illustrated in table 2, changing the machine does not have any effect on the design or features of the sample. the effective stress in msc simufact was implemented once the printed aircraft component had been compared to calculations regarding residual stresses using stoney’s equation. the effective stress of both the m280 and m400 machines is identical. for the alloys tial6v4, alsi10mg and 316l, the stresses were 1260 mpa, 370 mpa and 680 mpa, respectively. when the m280 machine was used, the surface deviation, representing the difference between the designed and printed samples, of all three alloys did not change. however, when the m400 machine was used, the surface deviation in each alloy varied. the laser power of the machine, which is responsible for fusing the powder, is one of the most crucial elements to modify. as a result, the printing temperature of a material, which is roughly twice as high as its melting point, must be higher than its melting point but lower than its boiling point. for each alloy, the changes in laser power as well as the measured temperature are shown in table 3. for tial6v4, 200 w is a suitable power, while the optimal laser power for alsi10mg and 316l is 100 w. the mechanical calibration was carried out to identify the inherent stresses. for measuring the distortion in 49(2) pp. 65–70 (2021) 70 alzyod and ficzere table 2: the effects of the machines and alloys. material eos effective surface used machine stress [mpa] deviation [mm] max min tial6v4 m280 1262.03 0.17 −0.16 m400 1288.13 0.88 −0.93 alsi10mg m280 370.01 0.13 −0.12 m400 370.02 0.36 0.26 316l m280 683.19 0.17 −0.16 m400 680.65 0.16 −0.15 table 3: adjusting the laser power. alloy used melting laser printing point (c°) power (w) temperature (c°) tial6v4 1600 200 2650 alsi10mg 660 200 4220 150 3098 100 2084 316l 1400 200 4454 150 3500 100 2463 the beam, two points on the cantilever were chosen at heights of 3.09 and 2.01 mm. the simulation repeatedly printed the specimen until the permissible distortion error became less than 3%. the sample was printed seven times to achieve an acceptable degree of distortion. table 4 shows the inherent strains of the three alloys. since the volumetric expansion factor is a crucial parameter in 3d printing, the simulation carried out the thermomechanical calibration to determine the volumetric expansion factor, as is shown in table 5. 4. conclusion the simulation tool simufact was used to print a component of an airplane out of different alloys using am technology to determine the specified parameter and carry out virtual printing. as indicated in table 3, since the laser power has a massive effect on the printing temperature, it should be regulated carefully. furthermore, as indicated in tables 4 and 5, in order to determine the inherent strains and volumetric expansion factor, the component must be printed multiple times. references [1] hernandez korner, m.e.; lambán, m.p.; albajez, j.a.; santolaria, j.; ng corrales, l.d.c.; royo, j.: systematic literature review: integration of additive manufacturing and industry 4.0. metals, 2020. 10(8),1061 doi: 10.3390/met10081061 [2] alex, s.; árpád, t.: xbrl utilization as an automated industry analysis. hung. j. ind. chem., 2020, 48(1), 131–138 doi: 10.33927/hjic-2020-19 table 4: inherent strains of tial6v4, alsi10mg and 316l. metal �xx �yy �zz ns tial6v4 -0.00488099 -0.0022865 -0.3 7 alsi10mg -0.00637062 -0.0016680 -0.3 18 316l -0.0198353 -0.0025170 -0.3 3 ns: number of samples printed table 5: volumetric expansion factors for tial6v4, alsi10mg and 316l. metal volumetric expansion factor number of steps tial6v4 0.821869 18 alsi10mg 0.60999 2 316l 0.307238 24 [3] attaran, m.: the rise of 3-d printing: the advantages of additive manufacturing over traditional manufacturing. bus. horiz., 2017, 60(5), 677–688 doi: 10.1016/j.bushor.2017.05.011 [4] gábora, a.; varga, t. a.; kozma, i.; beke, s.; mankovits, t.: 3d geometrical reconstruction of closed-cell aluminum foam using ct images. in incell 2019: book of abstracts of the international conference on multifunctional cellular materials (ua editora universidade de aveiro.), page 74 isbn: 978-972-789-611-0 https://ria.ua.pt [5] li, j.z.; alkahari, m.r.; rosli, n.a.b.; hasan, r.; sudin, m.n.; ramli, f.r.: review of wire arc additive manufacturing for 3d metal printing, int. j. autom. technol., 2019, 13(3), 346–353 doi: 10.20965/ijat.2019.p0346 [6] loughborough, u.: the 7 categories of additive manufacturing, 2019, https://www.lboro.ac.uk [7] software solutions, s.s.s.a.s.: additive manufacturing technology, 2017, https://www.simufact.com [8] mankovits, t.: numerical analysis of unit cell models for orthopedic applications, iop conf. ser.: mater. sci. eng., 393, 012019 doi: 10.1088/1757899x/393/1/012019 [9] software solutions, s.s.s.a.s.: english webinar on the entire portfolio of msc software for additive manufacturing, 2017, https://www.simufact.com [10] rasiya, g.; shukla, a.; saran, k.: additive manufacturing-a review, mater. today: proc., 2021, 47(1), 6896–6901 doi: 10.1016/j.matpr.2021.05.181 [11] schelhorn, l.; gosch, m.; debeugny, l.; schroter, p.; schwarz, w.; soller, s.: optimal design and process simulation for additive manufacturing, in 8th european conference for aeronautics and space sciences, eucass2019-fp0354, madrid, spain, 1–15 doi: 10.13009/eucass2019-354 hungarian journal of industry and chemistry https://doi.org/10.3390/met10081061 https://doi.org/10.33927/hjic-2020-19 https://doi.org/10.1016/j.bushor.2017.05.011 https://ria.ua.pt/bitstream/10773/26710/1/abstractbook-incell2019.pdf https://doi.org/10.20965/ijat.2019.p0346 https://doi.org/10.20965/ijat.2019.p0346 https://www.lboro.ac.uk/research/amrg/about/the7categoriesofadditivemanufacturing/ https://www.simufact.com/additive-manufacturing.html#challenges-in-metal-am https://doi.org/10.1088/1757-899x/393/1/012019 https://doi.org/10.1088/1757-899x/393/1/012019 https://www.simufact.com/simufact-additive.html https://doi.org/10.1016/j.matpr.2021.05.181 https://doi.org/10.1016/j.matpr.2021.05.181 https://doi.org/10.13009/eucass2019-354 introduction part-building technology fields of application modeling approach in simufact design methodology and simulation sample selection material selection machine selection thermal analysis mechanical calibration thermomechanical calibration fractions of exposure energy volumetric expansion factor results conclusion hungarian journal of industry and chemistry vol. 50(2) pp. 17-21 (2022) hjic.mk.uni-pannon.hu doi: 10.33927/hjic-2022-13 investigations into the usage of the mineral alginite fermented with lactobacillus paracasei for cosmetic purposes pál tóth1 and áron németh1* 1 department of applied biotechnology and food science, budapest university of technology and economics, műegyetem rkp. 3, budapest, 1111, hungary a remarkable interplay between the skin and the fermentation of lactic acid bacteria (lab) occurs. the lactate and amino acids in the supernatant of this bacteria help to hydrate the skin. the fermentation broth of lactic acid bacteria, generally referred as "lactic acid bacteria ferment" according to international nomenclature of cosmetic ingredients (inci), has been used to make a variety of cosmetic components. the goal of this study was to evaluate new approaches to assess ferment filtrates (also known as cell-free fermentation broths) that could be used in cosmeceuticals. despite years of research on the production of lactic acid, aesthetic evaluations have not yet been performed. first, the lactobacillus paracasei strain ncaim b.01525 was employed in our research. furthermore, a fermentation broth was produced containing the hungaricum mineral alginite and the impact of hydration on human skin tested. the influence of alginite on the fermentation of lab was also examined. according to the results of the trials, although alginite might double and triple biomass and specific growth rates, respectively, it cannot facilitate hydration of the skin. these results might contribute to the development of more widely accessible, environmentally-friendly cosmetic components in the future. keywords: human skin, dermatoscope, moisturizing effect 1. introduction the use of cosmetics dates back 7,000 years to ancient egypt and were later used in the roman empire, for example, according to some legends, cleopatra bathed in goats' milk to preserve the youthfulness of her skin. later, milk baths became popular in the english royal court, for instance, catherine parr, the last queen of the house of tudor, and later elizabeth i, queen of england, regularly used a milk bath to preserve their beauty. in today's modern world, using natural ingredients in cosmetics has come to the fore again. although our laboratory has been dealing with lactic acid producers for many years, to date, no measurements have been made in terms of cosmetics. the structure of the skin and the composition of the stratum corneum (sc) are shown in fig. 1, including the natural moisturizing factor (nmf). the skin consists of two main layers, namely the dermis and epidermis. the epidermis is further divided into two main layers, that is, the viable epidermis and the stratum corneum (sc), which is the outermost layer mainly consisting of dead cells. the essential function of the sc is to act as a barrier, preventing dehydration caused by water loss from the body. the sc contains 30% nmf, which received: 2 sept 2022; revised: 15 oct 2022; accepted: 15 oct 2022 *correspondence: naron@f-labor.mkt.bme.hu consists of 40% amino acids such as serine, glycine and alanine as well as 12% lactate, which is capable of retaining water in the stratum corneum [2]. on the other hand, lactic acid bacteria (lab) are generally defined as gram-positive, non-spore-forming, catalase-negative, aerotolerant, acid-tolerant, nutrientdemanding and strictly fermentative organisms that lack cytochromes as well as produce lactic acid as the main end product of carbohydrate metabolism [3]. lab contains cell wall-bound proteinase that initiates the transformation of extracellular proteins into oligopeptides [4]. the protease activity as well as figure 1. composition of the skin and its natural moisturizing factor (nmf), where pca is pyrrolidone carboxylic acid [1] https://doi.org/10.33927/hjic-2022-13 mailto:naron@f-labor.mkt.bme.hu tóth and németh hungarian journal of industry and chemistry 18 catabolic production of proteins and peptides enable lab to produce amino acids in the fermentation broth [5]. since many amino acids are contained in various lab-fermented foods such as cheese [6], sausage [7] and even japanese sake [8], a lab-fermented ingredient contains large amounts of lactic acid and amino acids, which together form the nmf. therefore, these ingredients exert a hydration effect when they are applied on the skin. furthermore, a particular combination of a substrate and a strain of lab could provide another beneficial functional impact on the skin. our goal was to investigate how the mineral alginite affects fermentation and moisturization of the skin. another aim of this study was to screen for new methods to evaluate the ferment filtrate, i.e. cell-free broth, in terms of usage in cosmeceuticals. 2. experimental 2.1. fermentation the planned experiment concerning the production of lactobacillus ferment filtrate was performed using a strain of lactobacillus paracasei from the national collection of agricultural and industrial microorganisms (ncaim, hungary), namely ncaim lactobacillus paracasei b.01525. the fermentations were conducted in 250 ml shake flasks and 10 ml bactrac impedance tubes (sy-lab, austria fig. 2) at 37°c. only the shake flasks were shaken at 150 rpm. the fermentations could be monitored online using a sy-lab bactrac 4100 microbiological analyzer due to changes in the impedance of the medium (m%) and on the surface (e%) of the electrodes by following the same method previously reported by áron németh [9]. the measurements were performed using the bactrac analyzer and the results displayed by bacmonitor y 1.39er software. since the bactrac program only shows curves and does not give the corresponding points, these points had to be extracted in a different way. the curves were saved as qrp files, which were converted into jpeg format using smartqrp software. the jpeg files had already been recognised by digitizeit, which made it easy to determine the points associated with the curves. the points were copied into a microsoft excel spreadsheet and curves were fitted to the points using sigmaplot 2001 software for windows version 7.0. the generalized logistic equation by verhulst-pearl was used [10]. the nutrient solution was the following marierogosa-sharpe medium: peptone (10.0 g/l), yeast extract (5.0 g/l), beef extract (10.0 g/l), glucose (20.0 g/l), kh2po4 (2.0 g/l), sodium acetate (5.0 g/l), magnesium sulfate (0.2 g/l), manganese sulfate (0.05 g/l), tween 80 (1.08 g/l) and ammonium citrate (2.0 g/l). the fermentations with alginite were supplemented with 10g/l of alginite. during the fermentation, the changes in glucose and lactic acid concentrations were monitored by a waters breeze hplc system consisting of a bio-rad aminex hpx-87h column, waters 717 plus autosampler and waters 2414 ri detector. after appropriate dilution steps, the samples were mixed with a 0.2 µm-pore-size mixed ester syringe filter (vialab magyarország kft.). the dry matter was measured as follows: the contents of the bactrac tubes, that is, 10 ml of fermentation broth, were loaded into falcon conical centrifuge tubes and centrifuged at 6000 rpm for 10 mins. using a hermle z200a centrifuge. the supernatant was then decanted before the cells were suspended in distilled water and centrifuged once more. the supernatant was decanted again and the biomass poured into the crystallization cup with 2-3 ml of distilled water before being dried relatively quickly using our sartorius ma35 moisture analyzer. 2.2. determination of moisture content the short-term/immediate hydration effect of the ferment filtrates was determined by taking triplicate measurements using a dermatoscope (fig. 3). on the forearm of the subject, three 1 cm2 areas were marked out onto which 20 µl of fermented juice was pipetted. after 5 mins., these areas were wiped with a dry hand towel and then the hydration on that part of the subject’s skin was measured at given intervals with the installed corneometer (capacitance) sensor. in order to establish a basis for comparison, the hydration of the subject’s skin figure 2. bactrac and its vials with electrodes figure 3. measuring with a dermatoscope mineral alginite fermented with lactobacillus paracasei 50(2) pp. 17-21 (2022) 19 was recorded before measurements were taken and the values displayed here corrected according to this value. 3. results and discussion 3.1. bactrac results the effect of alginite on the kinetics of lab fermentation was also investigated. the fermentations were monitored by bactrac online by replicating both setups three times, that is, in the presence and and absence of alginite, supplemented with a noninoculated (blank) reference. the reference bactrac tube contained the medium mrs and 10 g/l of noninoculated alginite. the replicates were compared by excluding outliers from further processing and their curves are presented in fig. 4. the resulting curves clearly show the impedance of lactobacillus paracasei fermentations in the presence of alginite during the declining growth phase is twice as high as that of nonalginite fermentations. furthermore, the specific growth rate was also determined to quantitatively compare the effect of alginite on fermentation kinetics as is presented in fig. 5. the specific growth rates also differ significantly between the normal and alginite fermentations, namely 0.6 and 1.9 l/h, respectively (fig. 5). therefore, the specific growth rate of fermentation with alginite is three times higher than without. the data obtained from dry matter measurements also indicate that alginite fermentations achieved higher yields (fig. 6). 3.2. skin hydration since a lactic acid-producing strain was investigated, firstly the moisturizing effect of solutions of different lactic acid concentrations, namely 5, 10, 15 and 20 g/l of lactic acid in distilled water, was determined by a dermatoscope. how skin hydration changes over time once the given concentration of the lactic acid solution had been reduced is represented in fig. 7. it is clearly visible that the level of hydration rapidly decreases over time but eventually stabilizes after ca. 35 mins. since the initial phase is difficult to quantify, the steady-state values after 30, 35, 40 and 45 minutes were used before being plotted against concentration as presented in fig. 8. the trend is clear, that is, as the concentration of lactic acid increases, so does the moisturizing effect. the effect of the 5g/l lactic acid solution is so minimal that it has a relatively mild dehydration effect on the skin. a trendline was fitted to the points (r^2>0.99) and from the resulting equation, the hydration effect of the fermentation broth was predicted based on its lactic acid content. while our technique for measuring the moisturizing effect is first described and applied here, the lactic acid sting test (last), which also measures capacitance figure 4. lactobacillus paracasei fermentations in bactrac tubes figure 5. specific growth rates of l. paracasei fermentations figure 6. concentration of dry matter in the l. paracasei fermentations figure 7. the moisturizing effect of the lactic acid (la) solutions tóth and németh hungarian journal of industry and chemistry 20 using a corneometer (similarly to our dermatoscope’s corneometer) is highly comparable to it [11]-[13]. the only two differences are that the cited authors applied the ferment filtrates to the facial skin, to nasolabial folds to be exact, where the hydration effect was 6 times higher in contrast to on the forearm in our study, and that they only reported the measured level of hydration after 10 mins., which according to our time profile varies rapidly over time suggesting it may be more reliable. overall, a correlation with regard to the lactic acid solutions is observed between our results and those from the cited studies. the effect on skin hydration of the lactobacillus paracasei fermented broths was also measured and plotted against time (fig. 9). a comparison concerning the steady-state values of the fermentation samples in the presence and absence of alginite is presented in fig. 10. according to the results, the samples containing alginite were less hydrating than the cell-free fermentation broth. this may be due to the higher concentrations of lactic acid produced in the alginite-free fermentation broth, namely 23.8 g/l and 28.8 g/l, respectively. by applying the equation in fig. 5, hydration values of 15.2 and 11.7% were predicted for the alginite-free and alginitecontaining fermentation broths, respectively, based on these lactic acid concentrations. although the measured results in both cases are slightly lower, the difference is insignificant. 4. conclusion the current study investigated the effects of the mineral alginite on the fermentation of lactobacillus paracasei and its moisturizing effects on the skin. this study also aimed to test novel approaches for evaluating ferment filtrates (also known as cell-free fermentation broths) for use in cosmeceuticals. according to studies on the process, alginite increases the specific growth rate and dry matter content during fermentation. while the latter doubled, from 1.6 to 3.3 g/l to be exact, the former tripled from 0.60 to 1.91 l/h. based on these encouraging findings, it was predicted that the hydration effect of the filtrate in the presence of alginite would outweigh that without. however, dermatoscopic measurements provided evidence to the contrary. the measured hydration levels of the alginite-free and alginitecontaining fermentation broths were 13.9 and 8.1%, respectively. our future studies will focus on the effects of fermented broths containing alginite on additional aesthetic aspects such as antioxidants, skin whitening or the inhibition of the enzyme hyaluronidase. references [1] izawa, n.; sone, t.: cosmetic ingredients fermented by lactic acid bacteria, in: microbial production: from genome design to cell engineering, (springer, tokyo) 2014, pp. 233–242, doi: 10.1007/978-4-431-54607-8_20 [2] spier, h.w.; pascher, g.: zur analytischen und funktionellen physiologie der hautoberfl äche (analytical and functional physiology of the skin surface), hautarzt, 1956, 7(2), 55–60, pmid: 13318589 [3] axelsson, l.: lactic acid bacteria: classification and physiology, in: salminen, s.; von wright, a. (eds) lactic acid bacteria: microbiological and functional aspects, (marcel dekker, inc., new york) 1998, pp. 1–71, isbn: 0-8247-5332-1 [4] savijoki, k.; ingmer, h.; varmanen, p.: proteolytic systems of lactic acid bacteria, appl. microbiol. biotechnol., 2006, 71, 394–406, doi: 10.1007/s00253006-0427-1 figure 10. lactobacillus paracasei ferment filtrates moisturizing comparison in steady state figure 8. steady-state moisturizing effects of la solutions figure 9. a comparison of the moisturizing effect of lactobacillus paracasei ferment filtrates https://doi.org/10.1007/978-4-431-54607-8_20 https://doi.org/10.1007/978-4-431-54607-8_20 https://doi.org/10.1007/s00253-006-0427-1 https://doi.org/10.1007/s00253-006-0427-1 mineral alginite fermented with lactobacillus paracasei 50(2) pp. 17-21 (2022) 21 [5] lee, k.; lee, j.; kim, y.h.; moon, s.h.; park, y.h.: unique properties of four lactobacilli in amino acid production and symbiotic mixed culture for lactic acid biosynthesis, curr. microbiol., 2001, 43, 383– 390, doi: 10.1007/s002840010324 [6] visser, s.: proteolytic enzymes and their relation to cheese ripening and flavor: an overview, j. dairy science, 1993, 76, 329–350, doi: 10.3168/jds.s00220302(93)77354-3 [7] hierro, e.; de la hoz, l.; ordóñez, j.a.: contribution of the microbial and meat endogenous enzymes to the free amino acid and amine contents of dry fermented sausages, j. agric. food chem., 1999, 47, 1156–1161, doi: 10.1021/jf980834p [8] iwano, k.; takahashi, k.; ito, t.; nakazawa, n.: search for amino acids affecting the taste of japanese sake, j. brew. soc. jpn., 2004, 99(9), 659– 664, doi: 10.6013/jbrewsocjapan1988.99.659 [9] németh, á.: investigations into succinic acid fermentation, hung. j. ind. chem., 2019, 47(2), 1– 4, doi: 10.33927/hjic-2019-13 [10] edwardsa, w.h.; wilke, c.r.: mathematical representation of batch culture data, biotechnol. bioeng., 1968, 10, 205–232, doi: 10.1002/bit.260100208 [11] pan, y.; ma, x.; song, y.; zhao, j.; yan, s.: questionnaire and lactic acid sting test play different role on the assessment of sensitive skin: a cross-sectional study, clin. cosmet. investig. dermatology, 2021, 14, 1215–1225, doi: 10.2147/ccid.s325166 [12] an, s.; lee, e.; kim, s.; nam, g.; lee, h.; moon, s.; chang, i.: comparison and correlation between stinging responses to lactic acid and bioengineering parameters, contact derm., 2007, 57(3), 158–162, doi: 10.1111/j.1600-0536.2007.01182.x [13] wu, y.; wang, x.; zhou, y.; tan, y.; chen, d.; chen, y.; ye, m.: correlation between stinging, tewl and capacitance, skin res. technol., 2003, 9(2), 90–93, doi: 10.1034/j.1600-0846.2003.00026.x https://doi.org/10.1007/s002840010324 https://doi.org/10.3168/jds.s0022-0302(93)77354-3 https://doi.org/10.3168/jds.s0022-0302(93)77354-3 https://doi.org/10.1021/jf980834p https://doi.org/10.6013/jbrewsocjapan1988.99.659 https://doi.org/10.33927/hjic-2019-13 https://doi.org/10.1002/bit.260100208 https://doi.org/10.2147/ccid.s325166 https://doi.org/10.2147/ccid.s325166 https://doi.org/10.1111/j.1600-0536.2007.01182.x https://doi.org/10.1034/j.1600-0846.2003.00026.x hungarian journal of industrial chemistry veszprem . vol. 30. pp. 211-214 (2002) mathematical model of variable volume diafiltration m. n. tekic, z. zavarg6, d. krstic and m. djuric (faculty of technology, university ofnovi sad, bul. caralazara 1, 21000 novi sad, yugoslavia) received: april11, 2002; revised: august 13, 2002 the use of ultrafiltration as a separation technique is frequently accompanied by a diafiltration step to remove microsolutes. there are different ways to combine ultrafiltration and diafiltration to obtain the desired final concentration of components. variable volume diafiltration is a continuous process in which water is continuously added at a rate jess than the permeate flow rate so that the concentration of macrosolute continuously increases, and the final concentrations of macrosolute and microsolute can be reached simultaneously. in this paper a mathematical model of variable volume diafiltration is proposed. the model includes both initial volume and concentration, rejection coefficients, processing time as well as water volume added during diafiltration. the developed model could be served for design and optimization purposes. keywords: modeling, ultrafiltration, diafiltration introduction over two decades ultrafiltration has been recognized as a separation technique to fractionate species according to their size. this technique has several advantages over other separation techniques such as distillation, evaporation solvent extraction, chromatography, etc. ultrafiltration is widely used as a separation technique in many different fields such as chemical and biochemical engineering, food processing and pharmaceutical industry. in order to achieve macrosolute-microsolute separation different ultrafiltration processes can be used. the concentration of macrosolute by ultrafiltration is frequently accompanied by a diafiltration step to remove microsolute. it consists of a continuous or discontinuous addition of a pure solvent to retentate feed and can be applied to both batch and continuous ultrafiltration processes. both, batch and continuous diafilltration can be performed in different ways. batch diafiltration can be carried out depending on the way the pure solution is added as a continuous or discontinuous batch diafiltration. on the other hand, continuous diafiltration can be performed as a co-current or counter-current [1] process. an interesting process in which a pure solvent is continuously added at a rate less than permeate flow rate so that concentration of macrosolute and the removal of impurities occurs simultaneously. wa:c. proposed by jeffrin and charrier [2]. thi~ way a smgle diafiltration process with a continuou~ly decreasing volume makes possible that the final concentration!'> of components are reached simultaneously. the aim of this work is to develop a process model of a variable volume diafiltration. in the presemed model the initial solution concentration, initial volum\.,. retention, processing time as well as the membrane area are included. this model differs from the approach given by jaffrin and charrier [2] where the rejection coefficient was not included. thus, on the basis of incomplete macrosolute rejection and the constant flux assumption an analytical solution is obtained. model let us consider the diafiltration process schematically presented in fig. i. the equation of continuity is given by: dv -=qo-qf dt (1) here vis the volume of the solution in the tank and ql> and q, are diafiltration water flow rate and filtrate flow rate respectively. the mass balance of the macrosolute can be written as 212 d(vc) = qfc(lr) (2) dt table 1 values of model parameters and coefficients parameter unit value a ? 1 mj m3/m2s 2.5·10-5 vo m3 0.2 r 1 0.9-1 r 1 0-0.1 c/co 1 2-10 cclc1 2-10 pconst variable pvar 1 variable a 1 variable where c is concentration, and r is the macrosolute rejection coefficient. in order to obtain concentrated solution with the final volume \'j the flow rate of diafiltration water must be a constant smaller fraction of qf. qd =aqf a<1 after rearranging eq. (2) becomes 0 ~ ~ j d(lnc) = -jqf(~-r) dtj d(lnv) ~ 0 ~ from the third term of eq. ( 4) it follows v = v0 qf(l-a)t (3) (4) (5) and the first term on. the right hand side of eq.(4) becomes 's' qfc8 hydrocarbons) compounds. by trapping the compounds on a solid sorbent, this latter limitation is largely overcome. the concentrations of volatile organic hydrocarbons in a selected sample are given in table 3. different classes of organic compounds were found: alkanes, alkenes, aromatic hydrocarbons and terpenes. the precision of the measurements was ±5 to 20 %. conclusion different sampling techniques are available for the sampling of volatile organic compounds in air. this fact, in itself, indicates that each of them has its field of application. we have shown that each has advantages and limitations but they are not really comparable. the selection of a sampling technique is primarily governed by the objective of the analysis. if long-term timeweighted average concentrations of vocs are needed for monitoring, diffusive sampling is the only feasible method. sampling aboard an airplane requires instantaneous and more sensitive methods, such as grab sampling. active sampling is probably the most common of all air sampling techniques, 11ow being challenged by diffusive sampling on one side and automated semi-continuous sampling on another. one should remember, however that sampling is only a part of a system on which the quality of the analytical data relies. references 1. johansson i.: atmos. environ., 1978, 12, 1371 2. mm a.l., ghosh t.k., loyalka s.k. and warder r.c. jr.: indoor air quality and control, prentice-hall, englewood cliffs, nj, 1993 · 3. hough a.m and derwent r.g.: atmos. environ., 1987, 21, 2015-2033 4. lopez a., barthomeuf m.o. and huertas m.l.: atmos. environ., 1989,23, 1456-1478 5. mackenzffi a.r., harrison r.m., colbeck i. and hewrrrc.n.: atmos. environ., 1991, 25a, 351-359 6. warneck p.: chemistry of the natural atmosphere, new york, academic press, pp. 223-267, 1988 7. wentf.w.: air pollution, sci. a., 1955, 192,63-72 8. lamontage r.a., swinnerton j. and linnenboom v.j.: tehus, 1974, 26, 7h7 9. eicbmanr., ke'rseridis g., schebeskeg., jaenickbr., hahn j., warneck p. and junge c.: atmos. environ., 1980, 14, 695-703 10. gelencser a., siszler k. and hlavay j.: . environmental science & technology, 1997, 31,2869 page 50 page 51 page 52 page 53 hungarian journal of indus1rial chemistry veszprem vol. 30. pp. 299 ~ 303 (2002) hydrotreating of full range fcc gasoline j. hancsok, s. magyar and a. lengyel1 (department of hydrocarbon and coal processing, university of veszprem, h-8201 veszprem, p.o. box 158, hungary 1molhungarian oil and gas co., h-2443 szazhalombatta, p.o. box 1, hungary) received: november 20, 2002 sulphur content of engine gasoline must be reduced below 50 ppm in the european union from 2005, and gasoline containing 10 ppm sulphur will have tax allowance [1 ,2]. fcc gasoline is one of the blend stocks being applied in largest amount (about 20-50%). the sulphur content of this is significant (about 50-2000 ppm), therefore 50-95% of the sul~hur species of gasoline originates from this stream. selective hydrotreating of fcc gasoline may be a fav~urable techmque among the numerous new desulphurising methods. achievements of a research work, made for hydrotreating a £?11 r~ge of fcc gasoline, are presented in this paper. the authors were able to find combinations of process parameters bemg smtable to produce gasoline blend stock of 11 ppm sulphur content with only 2 units loss of octane number. keywords: fcc gasoline, desulphurization, olefin saturation, octane loss introduction further reduction of the automotive emission can be achieved effectively with complex development of fuels, engine construction, lubricants and other parts of vehicles (catalytic converter, tyre etc.). currently the sulphur specifications have dominant importance from the point of view of engine gasoline, because combustion products of sulphur species beside air pollutionare poison the vehicle catalysts. thus further significant reduction of sulphur content can be expected (table i) [3]. three main long-run methods are offered for reducing sulphur content of gasoline, each of which results in lower sulphur content of fcc gasoline: reduction by hydrofining of fcc feed, application of new catalysts and catalyst additives in the fcc unit and desulphurisation of fcc gasoline [4, 5]. these strategies can be used either in themselves or in combination. nevertheless, desulphurisation of fcc gasoline cannot generally be avoided to meet ultra low sulphur specifications of engine gasoline. the loss of octane number can be very significant (up to 10-15 units) applying conventional hydrotreating of fcc gasoline. accordingly, this process is not economical from two aspects: partly due to considerable loss of octane number, partly because significant hydrogen consumption. a number of desulphurising processes for fcc gasoline have ~n develope4 which can economically be used to produce low sulphur fcc gasoline with acceptable loss of octane number [6-1_7]. the new desulphurising processes are wldel y diversified in respect of their principle and technical configuration (selective hydrotreating, adsorption. extractive distillation, caustic extraction etc.). options for desulphurisation of fcc gasoline are summarized in table 2 [18]. some of these processes treat full range fcc gasoline, but others accomplish desulphurisation with only a portion of fcc gasoline. it is extremely important in the latter processes tha~ the colun:n for the distillation of gasoline has to be optunally destgned and the cut point well selected [19]. . fig.} illustrates the major optional pathways tor post-treating of fcc gasoline: the main features of the processes indicated on th1s figure . were . a~eady presented in table 2. in some cas~s th~ hght fracuon_ of fcc gasoline is sent to an ethenficatto~ or alkylatl?n unit for boosting the octane number, while the ~ea~1er fraction is hydrotreated. this kind of combmatmn reduces the overall octane loss of tx>shreating. these processes were not indicated on fig.l. during the research, the possibility of desulphurisation of a full. range fc~ gasoline on pt,pdlzeolite has been investigated. the rum of the study was to examine the effect of process parame~ers (primarily temperature and liquid hourly space veloc1ty} on the yield and quality of liquid product and to determine the advantageous process parameters. 300 region, country, state usa california eu japan process naphtha hydrotreating (nht) nht+octane increase selective nht selective nht + octane increase adsorption extractive distillation oxidation alkylation bio processing table 1 actual and planned motor gasoline specifications maximum sulphur content ppm maximum olefin content v n % (actual) (planned) (actual) (planned) 500 30 (2006) 25 no change 30 15 (end of2002) 4 no change 150 50 (2005) 18 still not decided 100 10 (2008?) not specified still not decided table 2 options for the desulphurisation of fcc gasoline key feature industrial application conventional yes zeolite+ isomerisation yes rt-225 yes dual catalyst yes catalytic distillation yes combination yes zn adsorbent yes alumina adsorbent pilot selective solvent sys. yes peroxyacid pilot ultrasou,nd pilot solid acid pilot bio catalysis no h2-consumption octane loss high high high low medium low medium low medium low medium low low low none none none low none low low low low low low low isal i octgaln i prlme-g+i scan fining conventional hydrotratlng i !sal i octgaln i prlm&-g+i name of process various octgain, isal scan fining prime-g+ cd hydro/ cdhds scanfining ii szorb irvad gt-desulf . . ced sulphco oats fig, j major optional pathways for the desu1phurisation of fcc gasoline licensors a number of firms exxonmobil, uop exxonmobil ifp cd tech exxonmobil philips alcoa gtc petro star bechtel bp enchira ~ydrogen, ·-t>~~ m1rogen 1 2 3 31 ~liquid product 301 fig.2 simplified drawing of the test apparatus. notations: 1, 6, 11, 13, 14, 18, 20, 22, 30, 34, 36, 37, 38: closing valves; 2, 8, 31, 39: control valves; 3, 7, 9, 15: manometers, 4: oxygen converter; 5: dryer; 10, 32: gas filter; 12: gas flow meter/controller; 16.23: back valve; 17, 19: liquid feeds burettes; 21: liquid pump; 24: pre-heater; 25: reactor, 26: sampling valve, 27,29: cooler, 28: separator; 33: pressure recorder; 35: pressure controller; 40: wet gas flow meter experimental apparatus desulphurisation of fcc gasoline has been carried out in a high-pressure reactor system (fig.2) at the department of hydrocarbon and coal processing, university of veszprem. this consists of a tubular reactor of 100 cm3 efficient volume and is free of back mixing. it contains the same equipments and devices applied in the reactor system of desulphurising plants (pumps, separators, heat exchangers, as well as regulators of temperature, pressure and gas flow). catalysts the hydrodesulphurising experiments were carried out on pt,pd/zeolite catalyst, applying 80 cm3 of it. feedstock as feedstock of the desulphurising experiments a full boiling range (data of simulated distillation: 6-228 °c) fcc gasoline were used. the major quality features are summarised in table 3. methods compositions of feedstock and liquid products were analysed by gas chromatography (chromocf ane} and the quality characteristics were calculated by a software from these compositions. composition of gas products was determined, by gas chromatography (astm d 5134-90). sulphur content was measured by pyrofluorescence method (astm d 5453). the experiments were carried out on catalyst of steady-state activity. by continuous operation. results and discussion process parameters of the experiments (table 4) were selected and based on literature data and on earlier results of the department. from the results of the experiments it can be stated that crack reactions resulting in lighter hydrocarbons have not proceeded in the investigated temperature range, because the yield of liquid products was high (>99.5 %) at every combination of process parameters. the degree of desulphurisation of fcc gasoline depended on the process conditions. sulphur content of the products (fig.3) became lower with increasing temperature and decreasing lhsv. the highest level of desulphurisation (80 %) was reached at 280 oc and irsv = 1.0 (fig.4). under these conditions the product contained 11 ppm sulphur. 302 table 3 main properties of the feedstock density (15,6°c), g/cm3 0.7423 sulphur, ppm 63 nitrogen, ppm 13 research octane number 93.4 motor octane number 81.7 (ron+ mon)/2 87.6 composition. % n-paraffins 4.0 i-paraffins 31.8 ole fins 24.9 aromatics 31.7 naphthenes 7.6 i 35 ~30+---~~~~~~~----------------~ .:r !~+---------~~~~~~~----------~ ~ 20 r-------------=~~...2~~::------~ ~ 15 5+---~--------~--~----~----~--~--~ 220 230 240 250 270 280 tih!lperatuno, •c fig.3 sulphur content of products as function of temperature 100 ~ 90 .: ao ~ ~ 70 i 00 l 50 40 :t: 30 20 220 230 fig.4 hydrodesulphurisation as function oftemperature taking into account the composition of the feedstock and products it can be stated that the olefin content of each product decreased in proportion to the feedstock. the degree of saturation of olefins as function of temperature is shown on fig.5. higher temperature and tower lhsv resulted in higher olefin saturation. the highest olefin saturation (approx. 50%) occwted when desulphurisation was the lowest. ~iainly paraffins have formed from the olefins. but in a less degree also aromatics and naphthenes. evaluating the change of the concentration of paraffins it was stated that more n~paraffins than i~paraffins were formed from the olefins. every product had lower iso/normal paraffin ratio than the feedstock (8.0). this ratio is presented on fig.6 as function of temperature. the ratio was lower at higher temperatures and lower lhsv. this can be attributed to thermodynamic reasons. because higher temperature binders isomerisation. table 4 applied process parameters parameter reaction temperature, oc reaction pressure, bar liquid hourly space velocity, h" 1 hzfhydrocarbon ratio, m3/m 3 60 50 ;.!: g 40 ~ 30 = c i 20 '6 10 0 220 230 240 250 260 270 property 230-280 30 1,0-3,0 300 280 290 temperature, •c fig.5 olefin saturation as function of temperature 7,0 .2 6,0 'e ~ 5,0 ; ~ .!! 4,0 3,0 220 230 240 250 260 270 280 290 temperatura, •c fig.6 iso/n-paraffin ratio as function of temperature 2,5 ~ z 0 2,0 ~ 0 1,5 !!:. ~ 1,0 . ~ 0,5 8 0,0 220 230 240 250 260 270 280 200 temperatuta, ·c fig.7 loss of octane number as function of temperature the outcome of the mentioned chemical changes was the lower octane number of the products. fig. 7 illustrates the loss of octane number as function of temperature. the largest loss of octane number (2 units in [ron+mon]/2: the average of ron and mon) was observed at 280 °c, lhsv = 1.0 h1• sensibility of every product became lower due to saturation of olefins. by desulphurisation to the same degree, the lowest loss of octane number could be reached with the largest lhsv (3.0 n·1). removal of the light fraction of fcc gasoline could result in significant reduction of the loss of octane number. or lower sulphur content could be reached with the same loss of octane number. conclusions from the results of the investigation carried out on pt,pd/zeolit catalyst with the use of full range fcc gasoline, having 63 ppm sulphur content it can be stated that at advantageous process conditions (280°c; jo bar; lhsv=l,o; h2/hc = 300) fcc gasoline of 11 ppm sulphur content can be produced with high yield and only 2 units decrease of octane number ([ron+mon]/2: average of ron and mon). about 75 %of the olefins are in the light fraction of the feedstock (below 70 °c). this light fraction may contain very little sulphur, because it was passed through a merox unit in the refinery, which extracts the mercaptans from the light fraction. below 70 oc thiophenes are not present. this means that if we would cut the feedstock at 70 oc and we would only hydrotreat the heavier fraction, significant octane loss reduction could be reached, but only about 2 or 3 ppm of sulphur wou~d bypass the desulphurisation with the light fractron. furthermore, we can raise the temperature of the reactor, and we can reach higher level of desulphurisation without facing further significant octane loss. however, we have to confirm this with further experiments, and this is the aim of our next research work. nevertheless, our results confirm the opinion that hydrodesulphurisation of fcc gasoline can mainly be accmr~plished effectively and economically by expedient refinement of light and heavy fractions gained by fractionation. acronyms fcc fluid catalytic cracking hcn heavy cracked naphtha lcn light cracked naphtha lhsv liquid hourly space velocity mcn medium cracked naphtha mon motor octane number ron research octane number references 1. dixon-decleve s.: world refining, 2001, 12(9), 8 303 2. anon.: oil gas european magazine 2001 27(1) 42-43 ' , ' 3. sweed n. h.: petroleum technology quarterly, autumn, 2001, 6(3), 45-51 4. reidt. a., brevord e. and laan m. n. t.: the challenge of meeting future gasoline specifications: pre-treating vs. post-treating options around the fccu, european catalyst technology conference, antwerp, 2001 5. bavaro v.: world refining, 2000, 10(2), 30-37 6. maple r. e.: hydrocarbon engineering, 2000. 5~,%~2 . 7. burneit p. a., huff g. a., pradhan v. r., glaseit j. a. and hurst p.: bp low gasoline technology oats™, ertc 5th annual meeting, rome,2000 8. gentry j., khanmamedov t., and wytcherley w.: hydrocarbon engineering, 2002, 7(2), 43-44 9. d~buissc~rt q., nocca j.l. and cariou j.p.: pnme-g+ : the key to fcc gasoline desulfurization, proceedings of the interfaces '2002 conference, budapest, 2002 10. upson l. l. and schnaith m. w.: petroleum and coal, 2001, 40(3), 139-146 11. shih s. s., owens p. j., palit s. and tryjanowski d. a.: mobil's octgain™ process: fcc gasoline desulfurization reaches a new performance level, npra 1999 annual meeting, san antonio, texas, 1999 12. stun1z g. f. and plantenga f. l.: new technologies to meet the low sulfur fuel challange. 17th world petroleum congress. block 2:excelling in refining and delivering qualily petrochemicals, rio de janerio, 2002 13. anon.: world refining, 2001, 12(8), 23 14. rock k. l.: cdhydro/cdhds for ultra low gasoline sulfur, ectc 2002, amsterdam. 2002 15. irvine r. l. and varraveto d. m.: petroleum technology quarterly, summer 1999, 37-44 16. gislason j.: hydrocarbon engineering, 2002, 7(2 l 39-42 17. turk b., gupta r. and arena b.. a new continuous catalytic process for desulfurization of syngas and hydrocarbons, 2002 npra annual meeting, san antonio, texas, 2002 18. o'connor p. and mayo s.: division of fuel chemistry preprints, 2001, 46(2), 381-386 19. goldens. w., hanson d., w. and fulton, s. a.: hydrocarbon processing, 2002, 81(2). 67-72 page 300 page 301 page 302 page 303 page 304 hungarian journal of industrial chemistry veszprem vol. 29. pp. 129134 (2001) identification of nonlinear systems using gaussian mixture of local models j. abonyi, t. chovan and f. szeifbrt (department of process engineering, university ofveszprem, p.o. box 158, h-8201, hungary) received: october 8, 2001 identification of operating regime based models of nonlinear dynamic systems is addressed. the operating regimes and the parameters of the local linear models are identified directly and simultaneously based on the expectation maximization (em) identification of gaussian mixture model (gmm). the proposed technique is demonstrated by means of the identification of a neutralization reaction in a continuously stirred tank reactor. keywords: operating regime based model, expectation maximization, takagi-sugeno fuzzy model, nonlinear system, neutralization reaction introduction the problem of a successful model based control application arises from difficulties in system modeling [1, 2]. this difficulty stems from lack of knowledge or understanding of the process to be controlled [3]. while it may not be possible to find process information that is universally applicable, .it would certainly be worthwhile to examine what types of process-knowledge would be most relevant for specific operating points of the process. this type of local understanding, in fact, will be a key to identifying reliable local models with a limited amount of data. the model that has a range of validity less than the operating regime of the process is called local model, as opposed to a global model that is valid in the full range of operation. global modeling is a complicated task because of the need to describe the interactions between a large number of phenomena that appear globally. local modeling, on the other hand, may be considerably simpler, because locally there may be a smaller number of phenomena that are relevant, and their interactions are simpler [4]. the modeling framework that is based on combining a number of local models. where each local model has a predefined operating region in which the local model is valid is called operating regime based model [5], where the local models are combined into a global model using an interpolation technique as it is illustrated in fig. i. the main advantage of this framework is its transparency. both the concept of operating regimes and the model structure are easy to understand. this is important, since the model structure can be interpreted in terms of operating regimes, but also quantitatively in terms of individual local models. the operating regime of the local models can be also represented by fuzzy sets [6]. this representation is appealing, since many systems change behaviors smoothly as a function of the operating point, and the soft transition between the regimes introduced by the fuzzy set representation captures this feature in an elegant fashion. fuzzy modeling and identification proved to be effective tools for the approximation of uncertain nonlinear systems because of the ability to combine expert knowledge and measured data. fuzzy models use if-then rules to describe the process through a collection of locally valid relationships. the antecedents (if-parts) of the rules divide the input space into several fuzzy subspaces, while the consequents (then-parts) describe the local behavior of the system in these fuzzy subspaces [7]. in this paper the local models are linear. the contribution of this paper is two-fold: • a new method for the identification of operating regime based models is proposed based on em identification of gaussian mixtures model. • method to transform the obtained mt1del imo takagi-sugeno fuzzy model is presented. the paper is organized as follows. section 2 pre~ents the structure of the operating regime based model along with the methods for its transformation into a fuzzy model. in section 3, the identification algorithm of the 130 y(k) where the l/j;(xk) function describes the operating regime of the i -th local linear model defined by the ei =[a; ,bi r parameter vector. the operating regime of the local models can also be represented by fuzzy sets [6]. hence, the entire global model eq.(3) can be conveniently represented by takagi-sugeno fuzzy rules: r;: if xk is a;(xk) then yk =a; xk +hi' i = l, ... ,c (4) u(k) where a; (xk) represents the multi variable membership function that describes the fuzzy set a; while a; and fig. i example for an operating regime based model. the operating region defined by the current input u(k) and output y(k) of the system is decomposed into four regimes. · model is proposed. an application example the identification of a ph process is given in section 3. conclusions are given in section 4. operating regime based modeling of dynamical system nonlinear dynamic systems are often represented in the nonlinear autoregressive with exogenous input (narx) model form, which establishes a nonlinear relationship between the past inputs and outputs and the predicted output: y(k +1) = f(y(k) •.•. ,y(k -ny),u(~-nd), .•. ,u(k -n")) (1) here, n.,. and n 11 denote the maximum lags considered for the output, and input terms, respectively, nd