it is well known that transformers are regarded as indispensable and crucial components in power systems . due to market globalization , and in some cases to accommodate particular specification requests this competitive environment mandates the adoption of design strategies yielding better performance at lower costs . in the past , adly and abd - el - hafiz   demonstrated that feed - forward neural networks may be utilized to predict design details of power transformers after being trained using dimensional and winding details of a set of actual transformers . alternatively , finite element analysis ( fea ) coupled to an educated trial and error approach was introduced . furthermore , a computer - aided trial search looping algorithm aiming at minimizing transformer design cost has been demonstrated . other approaches coupling fea to a knowledge - based design optimization strategy and genetic algorithms were presented [ 57 ] . hernndez and arjona   proposed another approach that couples classical design equations to an intelligent design search algorithm . a quick review of these methodologies reveals that a wide span of design strategies could be utilized to achieve an optimum power transformer design . for instance , analytical formulations may be utilized for the quick estimation of transformer dimensions and design details . methodologies based upon more accurate fea computations offer precise estimation of transformer performance measures , provided that design specifics are suggested a priori . other methodologies , on the other hand , may utilize a hybrid strategy or even non - traditional heuristic and/or evolutionary computation approaches . several techniques have addressed transformer design problems using single - objective particle swarm optimization ( pso ) . demonstrated that the two objectives of minimizing power loss and leakage inductance were combined into one objective function using weighted aggregation . single - objective evolutionary optimization was , then applied using a hybrid algorithm of pso and differential evolution . jalilvand and bagheri   also utilized single - objective pso in the optimal design of protective current transformers . the objectives of making current measurements more accurate and designing more efficient current transformers in terms of both size and cost were formulated as an optimization problem to be solved by pso . on the other hand , du et al .   focused on improving the standard single - objective pso algorithm and utilizing the improved version in the optimal design of rectifier transformers . the reduction of a multi - objective optimization problem to a single - objective problem is usually performed by constructing a weighted sum of the original objective functions . while such methods are easy to implement and use , it is difficult to determine the appropriate weight coefficients when enough information about the problem is not available . another drawback of such approaches is that several runs of the algorithm are needed in order to obtain a set of optimal compromise solutions to choose from . furthermore , some optimal solutions can not be obtained , in some cases , regardless of the weight combinations used . hence , multi - objective pso becomes useful as it enables finding several optimal compromise solutions in a single run of the algorithm instead of having to perform a series of separate runs as in the case of classical optimization methods . the purpose of this paper is to present a power transformer design methodology using multi - objective evolutionary optimization . using this methodology , which is tailored to be target performance design - oriented estimated design parameters and details using the proposed methodology may also be considered for further refinement by other fea approaches . it should be stated that while the proposed methodology is analytical in nature , some parameter range settings have utilized previously reported power transformer field computation results . details of the proposed methodology as well as sample design results are reported in the following sections . in addition to the mandated primary line voltage vl1 , secondary line voltage vl2 and supply frequency f , a three - phase power transformer design is usually optimized to meet volt - ampere rating s , total copper losses pcu , no - load losses pnl and equivalent reactance x requirements . in other words , a performance - oriented design problem reduces to the proper selection of windings and dimensional details that would lead to a set of targeted performance figures . expressions linking the above - mentioned performance figures to the windings and dimensional details of a three - phase power transformer may be deduced in a systematic way as given below ( please refer , for instance , to [ 1517]).(1)vl13=4.44fbkfkc4d2n1,where b is the core maximum flux density ( magnetic loading ) , kf is the laminations stack factor , kc is the gross area to maximum circular area ratio , d is the core bounding diameter and n1 is the primary winding number of turns . it is also known that the window space factor of a three - phase transformer sw may be expressed as:(2)sw=2n1ac1 + 2n2ac2hwww , where n2 is the secondary winding number of turns , hw is the window height , ww is the window width , while ac1 and ac2 represent the primary and secondary winding cross sectional areas , respectively . denoting the window height to width ratio by kw and assuming a common current density ( electric loading ) j in both windings while n1iph1  =  n2iph2 ( where iph1 and iph2 are the primary and secondary phase currents ) , expression ( 2 ) may be rewritten in the form:(3)sw=4n1ac1kwhw2 . it should be pointed out here that , usually , current densities in low and high voltage windings are not identical due to standard wire size availability and/or other design factor constraints . nevertheless , the assumed current density j may be regarded as an average figure for both windings . from expressions ( 1 ) and ( 3 ) , the volt - ampere rating of a three - phase transformer may thus be expressed as:(4)s=34.44fbkfkc4d2n1jswhw24kwn1=3.33kfkcswf4kwjbd2hw2 . namely , these three components are the ohmic winding losses pcu - ohmic , the eddy current losses in the windings pcu - eddy and the copper terminals connection losses pcu - con . while designing a transformer to meet pre - mandated specification , maintaining the total copper losses below the threshold values becomes a must . in order to achieve this goal alternatively , appropriate computational safety factors may be applied to fast analytical design methodologies . while pcu - con    0.05pcu - ohmic , eddy current losses in transformer windings are dependent on the window height to width ratio kw . as previously reported by saleh et al . , electromagnetic field computation results suggest that , taking 2    kw    2.5 , winding eddy current losses may be estimated as pcu - eddy    0.15pcu - ohmic . from the previous common electric loading assumptions and for the aforementioned kw range , total copper losses may thus be given by the expression:(5)pcu1.23iph12cun1lmt1ac1 + 3iph22cun2lmt2ac2=7.2j2cun1ac1lmt=1.2j2cuvolcu , where cu is the specific resistivity of copper , volcu is the overall copper volume , while lmt1 , lmt2 and lmt represent the average turn length of the primary , secondary and both windings , respectively . no - load losses pnl , on the other hand , may also be regarded as a superposition of two components . more specifically , these two components are the core losses pfe and stray losses pstray . once more , it should be stated that accurate estimation of the stray losses requires massive computational resources that involve complex models of coupled magnetic , thermal and mechanical variables ( please refer , for instance , to [ 1921 ] ) . by adopting standard fabrication methodologies [ 1517 ] , stray losses may be estimated in accordance with the inequality pstray    0.3pfe . consequently , an upper limit for the no - load losses may be expressed in the form:(6)pnl1.3pfe=1.3wfe(b , f)fevolfe1.3wfe(b , f)fekfkc4d2(3hw+4ww+6d)1.3wfe(b , f)fekfkc4d23kw+4kwhw+6d , where fe is the steel lamination density and wfe(b ,  f ) is the specific core losses as a function of flux density and frequency that may be deduced by referring to the core lamination specifications data sheet . please note that explicit function formulation for wfe(b ,  f ) is either given in manufacturers specification sheets or simply inferred by fitting reported curves . referring to , the equivalent transformer reactance may be computed from:(7)x=2fon12lmtlwha+b1+b23,where o is the permeability of free space , lwh is the windings height , a is the spacing between the low and high voltage windings , while b1 and b2 represent the gross primary and secondary winding thicknesses , respectively . following the assumption of identical current densities for both windings and assuming similar winding heights , winding thicknesses may be assumed equal such that b1  =  b2  =  b ( please refer to fig . 1 ) . usually , the spacing between high voltage ( outer ) windings is double the distance between a low voltage ( inner ) winding and its corresponding high voltage winding . in other words , 4(a  +  b ) . from practical industrial considerations a    b/4 . in this case , the winding thickness may be correlated to the window dimensions according to:(8)bww5hw5kw . denoting the winding height lwh to the window height hw ratio by kh and following the previously stated practical assumptions as well as ( 8) , expression ( 7 ) may be rewritten in the form:(9)x=2fon12(d+2b+a)khhwa+2b3=112fon1230khkwd+9hw20kw . from ( 3 ) , n1 may be expressed in the form:(10)n1=swhw24kwac1=swhw2j4kwiph1 . substituting ( 10 ) into ( 9 ) , we obtain:(11)x=112fosw2480khkw3iph12j2hw4d+9hw20kw . following the same window configuration assumptions , expression ( 5 ) ( 4 ) , ( 6 ) , ( 11 ) and ( 12 ) , it is clear that the target performance oriented three - phase transformer design problem may be reduced to the proper selection of four unknowns . namely , those unknowns are the current density j , the maximum core magnetic flux density b , the transformer core diameter d , and the window height hw . dividing ( 11 ) by ( 12 ) , we get:(13)xpcu=11fosw864cukhkw2iph12hw2 . consequently , hw may be deduced from the expression:(14)hw=864cukhkw2iph1211foswxpcu.after obtaining hw value , remaining unknowns may be deduced by solving ( 4 ) , ( 6 ) and ( 12 ) . given the highly nonlinear nature of the equations under consideration , multi - objective optimization is needed to achieve a minimum cost design subject to the range restrictions for unknowns j , b and d. details of the multi - objective problem formulation are given in the following section . in engineering design problems , computational models are often used to describe the complex behaviors of physical systems and optimal solutions are sought with respect to some performance criteria . hence , multi - objective optimization becomes useful in obtaining a set of optimal compromise solutions ( pareto - optimal front ) so that the designer can select the best choice . the basic concepts of multi - objective optimization are introduced using a d - dimensional search space , srd , and k objective functions defined over s as given by bui and alam : ( 15)f(x)=[f1(x),f2(x),  ,fk(x)],subject to m inequality constraints:(16)gi(x)0,i=1,  ,m . the objective functions fi(x ) may be conflicting with each other , thereby preventing the detection of a single global minimum at the same point in  . consequently , optimality of a solution in multi - objective problems is defined differently . a vector v  =  ( v1 ,  v2 ,   , vk ) is said to dominate a vector u  =  ( u1 ,  u2 ,   ,  uk ) for a multi - objective minimization problem if and only if vi    ui for all i  =  1 ,  2 ,  ... ,  k and vi  <  ui for at least one component , where k is the dimension of the objective space . a solution u    u , where u is the universe , is said to be pareto optimal if and only if there exists no other solution v    u , such that u is dominated by v. such solutions , u , are called non - dominated solutions . the set of all such non - dominated solutions constitutes the pareto - optimal front . for the transformer design approach under consideration , s , pcu , pnl and x multi - objective optimization is then utilized to determine the other leading design parameters ; j , b and d. hence , x  =  ( j ,  b ,  d ) . within the current implementation , expressions ( 15 ) and ( 16 ) are formulated as:(17)f(x)=scmp(x)-s24,(pcucmp(x)-pcu)4,pnlcmp(x)-pnl44,subject to the following x inequality constraints:(18)jlbldljbdjubudu , where s , pcucmp and pnlcmp are the computed volt - ampere rating , the computed total copper losses and the computed no - load losses , respectively . it should be pointed out that the inequality ranges given in ( 18 ) should be in accordance with the typical lower and upper limits of j , b and d for power transformers in the range under consideration . inspired by the behavior of bird flocks or insect swarms , kennedy and eberhart first proposed pso in 1995 . pso is a population based heuristic , where the population of the potential solutions is called a swarm and each individual solution within the swarm is called a particle . considering a d - dimensional search space , an ith particle is associated with a position in the search space xi  =  ( xi,1 ,   ,  xi , d ) , a velocity vi  =  ( vi,1 ,   ,  vi , d ) and an individual experience vector pbi  =  ( pbi,1 ,   ,  pbi , d ) storing the position corresponding to the particle s personal best performance . experience of the whole swarm is captured in the vector gb  =  ( gb1 ,   ,  gbd ) , which corresponds to the position of the global best performance in the swarm . the movement of a particle toward the optimum solution is governed by updating its velocity and position according to eqs . 2 , respectively . while the parameter w is the inertia weight , parameters c1 and c2 are acceleration coefficients . parameters r1 and r2 are random numbers , generated uniformly in the interval [ 0 ,  1 ] and are responsible for providing randomness to the flight of the swarm . the second term in eq . 2 is the cognition term , which takes into account only the particle s individual experience . the third term in eq . the values of c1 and c2 allow the particle to tune the cognition and social terms , respectively . a larger value of c1 allows exploration , while a larger value of c2 encourages exploitation . single objective pso has been successfully utilized in many engineering applications such as the optimization of devices and systems   and field computation in nonlinear magnetic media . in order to handle multi - objective optimization , several approaches adapt single objective pso using the pareto dominance concept to determine the best positions that will guide the swarm during search . additional criteria are also imposed to take into consideration further issues such as swarm diversity and pareto front spread . in this paper , the time variant multi - objective particle swarm optimization ( tv - mopso ) algorithm is utilized . to achieve good balance between exploration and exploitation of the search space , tv - mopso is adaptive in nature with respect to its inertia weight and acceleration coefficients . a mutation operator is incorporated to resolve the problem of premature convergence to the local pareto - optimal front that is often observed in multi - objective pso . an archive is also maintained to store the non - dominated solutions found during execution . the first step generates an initial swarm swrmo of size ms with zero velocities and random values for the coordinates from the respective domains of each dimension . an archive of maximum size ma is initialized to contain the non - dominated solutions from swrmo . the second step represents the main iteration cycle in which the swarm is updated , the archive is updated and the swarm is mutated at each iteration t. the swarm swrmt is updated , in step ( 2.1 ) of the algorithm , by updating the velocity and coordinates of each particle using eqs . 2 , respectively . to update the velocity , the global best solution is obtained from the archive using a diversity consideration . the method for computing diversity of the solutions is based on a nearest neighbor concept . the present solution is compared with the personal best solution , and replaces the latter only if it dominates that solution . moreover , these parameters include an inertia coefficient , w , a local acceleration coefficient , c1 , and a global acceleration coefficient , c2 . the inertia coefficient w is decreased linearly with each iteration from an initial value wi to a final value wf . the value of w at iteration number t is calculated as:(19)w=(wi - wf)(m - t)m+wf , where m is the maximum number of iterations . to compromise between exploration and exploitation of the search space , the cognitive acceleration coefficient c1 and the social acceleration coefficient c2 are varied linearly with each iteration as given by ( 20 ) and ( 21 ) , respectively . while c1 decreases from the initial value c1i to the final value c1f , c2 increases from c2i to c2f.(20)c1=(c1f - c1i)tm+c1i,(21)c2=(c2f - c2i)tm+c2i . the archive at is updated , in step ( 2.2 ) of the algorithm , by including the non - dominated solutions from the combined population of the swarm and the archive . if the size of the archive exceeds the maximum limit ( ma ) , it is truncated using the diversity consideration . to explore the search space to a greater extent , while obtaining better diversity , a mutation operator is used in step ( 2.3 ) of the algorithm shown in fig . mutation is performed with probability inversely proportional to the chromosome length d. given a particle p , a randomly chosen coordinate ( i.e. , variable ) of the particle , pk , is mutated as follows:(22)pk=pk+(t , pku - pk)if flip=0pk-(t , pk - pkl)if flip=1,where flip , pkl and pku denote the random event of returning 0 or 1 , the lower and the upper limits of pk . the function  is defined by:(23)(t , x)=x1-r(1-tm)q , where r is a random number in the range [ 0 ,  1 ] , m is the maximum number of iterations and t is the iteration number . the parameter q determines the mutation s dependence level on the iteration number . after executing the specified number of iterations this archive contains the final non - dominated front ( i.e. , pareto optimal front ) . to serve the testing and estimation purposes , the proposed design approach has been implemented in digital form . the methodology has been particularly utilized to design 2550  mva , 66  kv/11  kv , dyn11 , 50  hz power transformers subject to a variety of design performance constraints . as per practical transformer stacking and assembly measures for the mva range under consideration , it was decided to set throughout the computations sw  =  0.2 , kh  =  0.9 , kf  =  0.95 , and consider 11-step cores ( leading to kc  =  0.958 ) [ 1517 ] . common values for o , cu and fe were taken as 4    10 h / m , 2.1    10    m and 7.65    10  kg / m , respectively . using armco steel tran - cor - h0 carlite-3 core laminations , an expression for wfe(b ,  50  hz ) was inferred from data offered by the manufacturer ( please refer to ) . the typical values of jl , ju , bl , bu , dl and du for power transformers in the range under consideration are set as 1.1    10 , 3.2    10 , 1.0 , 1.8 , 0.1 , 0.7 , respectively . in order to test the proposed performance - oriented design methodology , two transformers ( rated 25  mva and 40  mva ) of known manufacturer design details and measured performances are considered . it should be stated here that a considerable number of units of these particular transformer designs , which obviously passed all standard routine tests , has been acquired and installed in several national and regional grid sub - stations . as previously discussed , measured performance figures of actual transformers are taken as the target design requirements for the design methodology . the tv - mopso algorithm is executed using a swarm of ms  =  50 particles , a maximum archive size of ma  =  200 and for m  =  1000 iterations . the parameters used in the reported results are wi  =  0.7 , wf  =  0.4 , c1i  =  2.5 , c1f  =  0.5 , c2i  =  0.5 , c2f  =  2.5 and q  =  5 . 3 . out of a set of design parameters inferred by the tv - mopso implementation , the design corresponding to a minimum iron core volume is taken as the optimum choice . comparisons between design parameters of the actual transformers and those proposed by the design methodology under consideration are given in tables 1 and 2 . variations between actual and computed performance ( as well as cost ) figures are also given in the same tables . with the exception of the suggested current density j , it is clear that the proposed approach leads to good qualitative and quantitative performance - oriented design results . moreover , the proposed higher j value by the suggested methodology may be regarded as a possible cost minimization option as indicated in tables 1 and 2 by the possible reduction in the transformers copper volume . using the proposed approach , computations are also carried out to investigate the design parameters deviation from those of the 40  mva test transformer as a result of changing volt ampere rating and power loss requirements . in the first computation set , design parameters corresponding to 30  mva and 50  mva transformer ratings having the same specifics , per - unit reactances and percentage total copper and no - load losses ( i.e. , efficiencies ) were computed . 4 , almost all design parameter values increase as the transformer rating is increased while maintaining the same percentage total copper and no - load losses . it should be mentioned here that hw variation is minimal in these cases since transformer voltages are assumed unchanged . this suggests that any rating variation will similarly affect the phase currents squared and total copper loss values , thus minimally affecting hw as indicated by expression ( 14 ) . in the second computation set , design parameters corresponding to 40  mva transformer ratings having the same specifics , per - unit reactances but while varying the percentage total copper and no - load losses ( i.e. , efficiencies ) are computed . as shown in fig . 5 , a smaller loss restriction is achieved by a larger size transformer with reduced current and flux density values . this is a particularly encountered design trade - off between capital and running costs for a power transformer . in this paper , a performance - oriented power transformer design methodology using multi - objective evolutionary optimization has been introduced in detail . experimental testing as well as other presented computational results clearly demonstrates the qualitative and quantitative accuracy of the methodology . one advantage of using multi - objective evolutionary optimization is that it deals simultaneously with a set of possible solutions ( i.e. , a population ) . this enables finding several members of the pareto front in a single run of the algorithm instead of having to perform a series of separate runs as in the case of classical optimization methods . these options can be extremely useful to minimize overall production costs in view of the changing global prices of different transformer components , especially copper and steel laminations . the proposed methodology may be easily utilized to obtain a quick first guess design details for more sophisticated design approaches such as those utilizing fea packages . moreover , in the presence of detailed design strategies , the proposed methodology may be easily improved to relax some assumptions by including those strategies . future work is planned to enhance the accuracy of the proposed methodology as well as to extend its applicability to cover more detailed transformer design aspects .
transformers are regarded as crucial components in power systems . due to market globalization 
 , power transformer manufacturers are facing an increasingly competitive environment that mandates the adoption of design strategies yielding better performance at lower costs . in this paper , a power transformer design methodology using multi - objective evolutionary optimization is proposed . using this methodology , which is tailored to be target performance design - oriented 
 , quick rough estimation of transformer design specifics may be inferred . 
 testing of the suggested approach revealed significant qualitative and quantitative match with measured design and performance values . 
 details of the proposed methodology as well as sample design results are reported in the paper .