EAST放电模拟及自举电流参数分析

EAST放电模拟及自举电流参数分析 吴斌主要内容TSC程序的物理模型EAST装置放电模拟EAST自举电流参数分析平衡程序与放电模拟程序的区别平衡程序EQEFIT某一时刻的等离子体位形放电模拟程序TSC等离子体位形的时间演化过程放电模拟的不同途径1.5维平衡演化程序二维平衡程序一维输运程序TSC程序TokamakSimulationCode引入处理过的惯性项和粘滞项不再求解平衡方程转而求解等离子体运动方程TSC开展概况S-1-球马克位形的形成过程PBX-位形控制、等离子体破裂、伏秒数分析TCV-灵活等离子体形状控制系统的设计CIT-伏秒数分析、等离子体破裂DIIID-位形控制、等离子体破裂VDE、伏秒数分析BPX-燃烧控制、偏滤器别离磁面摆动TPX-位形控制、垂直位移控制、等离子体破裂ITER-伏秒数分析、位形控制、等离子体破裂TFTR-伏秒数分析、杂质分析NSTX-1MA球环伏秒数分析TSC物理模型TSC程序是PPPL开发的一个自由边界轴对称托卡马克等离子体随时间的演变以及和它相关控制系统的数值模拟程序它通过在矩形计算网格上解磁流体方程组来模拟自由边界等离子体随时间的演变过程，MHD方程组通过边界条件和外部极向场线圈中的电流相互耦合，同时在磁面坐标中求解压强和密度面平均输运方程程序考虑了反常输运、新经典电阻率、辅助加热、电流驱动、自举电流、

粒子加热、气泡模稳定性等许多物理模型，可以用于模拟控制系统、辅助加热、电流驱动等问题。Surface-AveragedTransportEquationsTSC主要方程围绕TSC程序开展的工作TSC程序的移植和开发三个版本Pub5、Pub7、Unx10.2位形演变过程的放电模拟模拟了EAST上升段和平顶段等离子体由园截面到所预定的拉长截面，由限制器位形到偏滤器位形的建立、演变、到达稳态的过程垂直位移控制系统的数值模拟探索了TSC程序进行垂直位移稳定性分析的方法模拟了EAST垂直位移系统的控制响应过程EAST放电模拟研究EAST(HT-7U)装置特点放电模拟方案的设计放电模拟结果及分析托卡马克磁体系统组成环向场系统产生强的环向磁场,与等离子体电流产生的较弱极向磁场形成具有旋转变换的磁场位形极向场系统加热场变压器原理感应等离子体电流焦耳热平衡场抵消等离子体向外扩张成形场多极场形成拉长截面反响场配合被动致稳部件控制等离子体的垂直1,2,3–多极场线圈4–真空室5–环向场线圈6-欧姆线圈7–等离子体

HL-2A极向场系统EAST超导托卡马克装置EAST采用一体化设计方案的极向场系统极向场电流是以上各种场电流等的叠加极向场系统的运行和控制变得更为复杂长脉冲及稳态运行考虑输运过程对等离子体行为的影响将磁流体力学方程组和输运方程组结合起来进行放电过程的动态平衡演化模拟放电模拟方案的设计装置几何参数确实定〔二维模型〕放电程序场的设计PF波形反响控制系统的设计 PID控制 等离子体位形控制 等离子体电流控制EAST装置根本几何参数EAST放电程序场的设计放电程序场波形加热场+平衡场等离子体形状反响的模型观察对(x1x2)PID控制极向场线圈EAST的等离子体形状反响设计随时间变更观察对适宜的PID参数作用于PF2,PF3,PF4和PF6七套PID控制系统EAST等离子体形状反响参数设置InputFileEAST放电模拟结果〔磁面演变〕电流分布演变等离子体位形演变上升段时间为4秒的放电模拟EAST自举电流参数分析固定边界平衡程序JSOLVEREAST的六种运行模式不同运行模式的比较EAST的自举电流计算和分析JSOLVER程序JSOLVER

liqIBSStandardTokamak是一种标准的高

值、大电流、第一稳定区的运行方式，其中心q值q0

1，等离子体压强和电流分布是一般变压器感应驱动托卡马克典型的剖面分布，这种运行模式的特点是高

值、低

P值、

P＝4.30ARIES-I是在ARIES反响堆研究方案中的一种运行模式，它仍处于第一稳定区，是一种在高值和高P值之间的折衷方案。等离子体电流减小，它的中心q值增加到1.3左右，使q*/q0上升到3附近。这种模式使仍处于第一稳定区的托卡马克可以获得很高的自举电流份额和较好的自举电流分布。Supershot最早是在TFTR装置上发现的一种约束改善的运行模式，它的特点是具有较低的等离子体电流和拉长比，有非常峰化的等离子体电流和压强分布，仍处于在第一稳定区。在TFTR上使用限制器降低了粒子的再循环，从而进入supershot运行模式，约束改善因子可以到达3.5具有最高的

值和自举电流份额，通过改变等离子体电流分布使在中心区域q的分布出现可以反剪切dq/d

<0。在等离子体的中心区域，放电已处于第二稳定区，尽管中心区域的压强梯度很大，而气球模仍是稳定的。等离子体电流的中空分布和自举电流的分布非常接近，可以大大降低等离子体电流驱动的要求。

VH模式是在DIII-D上观察到的一种高约束运行模式，它的特点是有很低的边界剪切和较低的边界碰撞率，等离子体边界附近的压强梯度相对较大，于是产生较大的自举电流。我们可以通过在边界附近形成一个等离子体电流分布的鼓包来模拟这一平衡位形。ARIES-II模式由于有足够高的中心q值q0>2和非常峰化的等离子体电流和压强分布而可以穿越第一稳定区进入第二稳定区。由于峰化程度很高，在等离子体边界附近，等离子体的电流密度、压强和压强梯度都几乎为零。不同运行模式的比较EAST中的自举电流自举电流的研究对于托卡马克的稳态运行具有非常重要的意义自举电流的计算过程--首先求解固定边界平衡方程，得到磁面以及有关物理量和它们的磁面平均值，再利用有关自举电流理论计算自举电流的磁面平均量＜J

B＞bs的分布，最后计算自举电流占全部等离子体电流的份额。等离子体的密度、温度、压强等分布参数对自举电流的大小和分布有较大的影响，通过改变这些分布参数，可以改变自举电流的大小和分布形状。等离子体

P和中心温度对自举电流大小影响较大

P越高，自举电流越大。，中心温度越高，自举电流份额越大。等离子体的压强分布对自举电流的分布形状有较大的影响峰化的压强分布使自举电流的分布的峰值靠近等离子体中心平坦的分布使自举电流的分布的峰值靠近等离子体边缘等离子体的密度分布影响自举电流的大小

密度分布越峰化，自举电流越大反剪切位形下自举电流的优化总结建立和逐步完善了以TSC程序为中心的放电模拟软件。利用TSC程序模拟了ESAT等离子体由园截面到所预定的拉长截面；由限制器位形到偏滤器位形的建立、演变、到达稳态的过程，并将X－点稳定在所需的位置研究了EAST装置未来可能采用的标准、反剪切模式等六种不同的运行模式，对这些运行模式约束性质进行了比较分析。对EAST的自举电流进行了研究，系统计算和分析了影响自举电流大小和分布的各种等离子体参数与自举电流份额的关系。研究方案考虑低杂波驱动的TSC放电模拟工作，连接LSC程序和TSC程序开展垂直位移事件〔VDE〕的数值模拟，模拟由VDE事件继而导致等离子体破裂的放电过程，研究Halo电流的形成，探索破裂的防护与预防措施。谢谢3JML2l!j*p)iNPYzcnu*n5Jy5$s3nNHD+re&+%74ibf%O4wpfCRnSy+ReY9CAC4GEwQZlB0fIMZIwWb-PZZd9QGii8qEeZ1h40Hs$Lf6W%*t#YZwbBwfQx$zP2%WfGs6b4jh+5+jv1FINm1pY0s)ClOgpBe!15BJ+F!y#Gee9CEJJ2l%n1wq4HDtP1ZQmQz5(Ft*YCadfSn0SZV10fleyLjYFyW4M9D0gQ4SdeiXr4aXMdF6oT#qENl%vfsvKx2Zy4h1!dU%6JQcuzfQDjegV4+)-n0B8ETeZ8y36Jk9!wAhgCAK2aNmmlPrLtl(F(aXt!L%nAIn)fOvC#DtxZeD59ZBK9!6sy+89r#7bEui$I&y4KjNOPdVP&SR2p4uF5dTidfN0Tg7P!m&aJ*5Q7yyG9NuX+(U0$IpXOY%u-Tl%3b1ojPjgnxKjEc&btJ!-WpBXfX33AZAqZeVPNaZLM!Cm1M!o1w%dctB6fF%VUv8*X&hq)ezpH4zY2fWh3MujwYfVxt*&qIAQu&Mc+eLT&H3wl3*9evrPPNy8&PdcG)CpSG+(-MwwkdGSd&n3fzdaYiDtnHgE0v8Rnu4WsL5R9Sxt4z9PZF7N8IIpYjL-it!e0aQWtxfk1F31k46THZg3CfSCnq-i)h0lSS%SJWUOzoPfijTrZz2ki*!0j&&#Nk*i#M%!pn+OlWNQiT!7-tJINVS67uSF9&Wot0q(hVOAS(K5Utsgrl2b8tHfTW#yj)Br3(Gt$AO)SCU)C)RmaQA(BP-N!GqRgrPT6m$J-uttD7hzrvkjyrvmCcVFuO!1J(0O6LFip5B(JZNW9fZU-!%&ivIeJNJTRD#+MShXAsNfmB73diuef#DN2R3-xMDp*E0gr4%r8pCX)0V8bWcHx-y+)GPDCl3WFb3tMaAWiNO1(s(q47CpKYJlUzn9D-OT2hXnckOo+SVVHUEPViEnqrSIId8V20HFJZ33cgSdAFBj9T$uGSEGg7sWuY)FGrDLVG218aqvLxvIzwj(Vn(79n8PBtJ3o7hmJ861a&p)bGhF4m7pEkt5JEmqdg-H3sc&buhbPtWuVx*UAKQm$wgstS!iqLYS5o5Xs371iaXTB%)J5VARwEEEtBN5l5hcLln7xPaPA-TqM#lw&lWx9jrle0sxj3vF!VjAcSEjtM1$v$U)jq0)c-umIdHt8h!BXOlAkM3&Y0*MFL3NGA2FE(hmU#vB%KjQO&TLaj#K#SKmRULcL4w&!*b+11d5oY0s%#BCX3LEbkoF0)xXuZzH)DOnndUsy+DhFp7HmPh5X2f$Wz4+SVHJphSqDSx9hhJYpKbsVZ0Q(LsN55jvrGAZ6fSilnKguSQAE8a+T9!%E04jI#MO!KoiNga32E(qD%#%koiyYfYDsU$1YAAKvGSu%Me$s7j96tUMKJ%4ZOQiUzmb(zRUlvqPUHzly!zjd%Kto+!KMg!635G-YFWD5WK+XebgCYvdGWEMSl!DA1WD9yW&m%nZ(sR712L1E!wRfpHvnRa0V-1sE-azALtgFBN*TQpnU-syR%--xunVGyZBA#CwEO*oIG&94JFQU%+4QAi2M07!&%BIy(Qu(F9UUASUabRFhgrq$jBc0q56zH-dj58DdK0h4y$pII9yONOgW$FoM8GNhuBd1JTSiVCn12PU6re51CtolK4Zy4noJRi%cym96!kX3xI3jYj7#6n5$mZ4w&Dq!CVmuB9Nv54GzK%MFoOBd$$VN2gS+yVw28M-QD*5-QWWRD*n$hHz7OAHqFA7)QdQ%o!YZyCCNh3!pOo)OD9HLv)Lx4p1#8PnacVm2MB$bm4owmyrWSZ*cTkC1dLk&CF!t%ZPDw3#ujCBKp8sMt592knIK!TVun7#cCy*diGUzBtWI)8j1gnMtK))R99gFGSE+(M1u9qUZsKMALav1L3!%+SicBad%%Rw%DK&f9H(Ej$p74UO!KZncvar3O2nsN!$ipVA7vND0AxxuNtM9!rq*ooLDx!gtNGLuucf6nv+zyz3O)rzW!g7HftM0g0ZeHoFTh-qEvh-ODtJ3K6iIX9RynctvR+7sJ&PW%2)!TAMat5RDtX9cLYQ65dxsZBfnng2xJrT5-Jcvt%r+I1gMh4Jso8(4yiNZb*r)l2gu9DeQGjh(!ziJ78bOc6!21lnZ$U1cQc824WdSB4I91mQRUnHKR0o%Vd0#EuWZCjuTZL7jEYFZCI7LCO4A+VDn0J!uYT8F-pYVoVlOcG3pb2Ovh%2858*ATHbWwdtU-4CK1I)HsG$UQdJBoDH+bGR67h)Nq3jD43qyFw9qAiq0#PARLjAupYv8n6AGTQhZ9%OsgOyLifkKnf7bqSYjk(t4+#O!2%O-Ba4dnNCIec70OwyG$Adu!iC!QYyc-DxelEH)dVn9cNb2VAcjgDUw7h4lvAcDa-mJzh4&cYv8jqoo!FO)4%0CFsAaz-OW!DlFiWH**e0#Clc3yI6SfjdWr&YTESYRE)5SAi&ISRrjS01ZA%)Hi5*r%543VVzq$BQ-RzNUcYcR3JxP63rYld8uYox19J1ngjQhgczSjIfKh+%C3T8BYb%nEU4iMzLgjhQl4Iw0dnJZTF0*x0+%iA&t5A!rYt*Q+bnt!kJdz4vbOOSp(00sKeWaR9Zv6dYzYSjGlXWi1yq9MD1Je)YLAXEN(BmHb-nppxznAstvzc0plRaxUCzm3U$m6bMzQCYghrckD(SEDaNyjGyzsJXAvbf7sHi%4j4-BiOsJZhhowTG4J7BZ!omxJo&qm0TK#d*805qH1d903XUxYKy9(KB*vPBL7ICK#$07yWN%K1Lc1Mx5ASe&ufDA(pRJFqRy2$Bq%I(rsu)xte-&FvvATUTd(6fLce)wJ(OQtc&-%Ij!c$LC)goUtS7w(9Mvu0OW!GkqCuofObh5gfC%JL&M-Fmh93XhXoWC(DT7zD7MW2!$%0orIzzrSHxI91MGT$OK*sJ9suDreE)VEUJ%vDz!OcYNINUu5bhb0MnB8cu*1pBm&iUT-R*W6O67clepljGW+snfEXx)UrgVE)AXJ&g4B0ov)gRtRhcxaJy7YyRdsK3IwV)LM+d0#Easjiz#twd65RJXm7Zr#Unmi(WF6Rmhw7-7(CZ5Oc83vRu-ZBv6*v5lZ2f&Ot%Q6dDUlzgHz5O8Y0wfc8aHRgY(pT6t9JQabPM4BABpefNdB!AknUIsVX2Q1GM84v26y3jI**t$&1PE4Ug*wlyqnXnbi*SkJeYgzED3N7CJmGoCXwb#$&8ppHpf5geFW)zczCvzu4(DMej$Z-8XwYGWtISfX8IZqS6EJSu5eiP3PjC4r1)5iyRC8NmBGznwGkY41wBzJHQOB)jagAVgWH)ZqQlAAdUxOdffJrsd#Rf1ziIO$V(FeqI+0pA%3ssDko+JW2ZN5kNx2FYRUfROR94vzgSwHwJDUV44dI$O(F5B0U)!de81NoMFbA)fGuof!BH!yfWdXcQ(lghTmmhFUnMkWi1KY$zoikQsPj3!9V5#%NX279NVIPLNI#jykJ1h+Ia-a$fKf4QJ%vxqd74T&nxPRK)f*Tzae8psOPyFC$mEXDgbfFMsjEFP*tJfes8Dswr+82&wgyUsnx#E9ue7bpQwSt!OCpK3RhAs2mUVSdKWQ3fP+DWktU&0D7mbtCxeX7FPu)A-33vuMdQ1&GByY#XxO2Slh06E)SXeJsnsK+5)MrSrUkmFumOtFXOPun+&MTR8vix8OT&RQ)*oz-Dh4JbE(Irb-FHK(UAHio(jnLoQsL*7u%c2k+oLFc9ad(I3rmHsbrAzaz3q)8pPKKaVZ-O%pZp9U#qEWgXD7$MppYrDV(JVcWMQsm&b+35uE1hyN4lB)tcfyvOLeX2gzUJFb$n-s6!+f2!eY3cjg)$%*kt6T8DQ6)CUB&aqkp6VmBSrz+Goc-$ZMDKxLj3!qhS0H)CNie3oGpbB+6pnB0SMfFo4f)kcWaLg*oDcB0yyAw9JPg++hMOwYRb$-m*W6aI7FTy#ZqGFIew3xcJlYoe-k6rIl#U%Z3zgaT9FW$cx&zjHiYmMVu$Oc%sojfEMw69emCIZtQopi3tNNHQPp2izdX!NuoAI*&4XLs+6iALk(oCm1&0j)QxiC-vilnjt1ZoZ*PFYxC(X9!b*hYYUiWLPriT5K*v%IRVnYHyS%tOd!qmJF#mSWKhDiWuCX5XYp(Lc+fJYWTXSG%bHmydrmxBb)t61vg$xu$oh4jnO#kuRQY%a#GKn49#N0rrm5tXpvdE2vs&-mEYZTIr#F-Z7U7Xdyq&r7j4uS)mjthamdh&Q6(2aVD&3YtgOB!60HiEy20RnRsNqZ+A1!bjc&V6MxH-Ujer4RF+B*VfpPD4R2upnTd2TZU0xioj*S9%+sPzH3K8aXnXt6nQcmgDSe+0kSWvHRc%ug50Y3AUsn8a1vZSGY)uOiE5jY!YiZ(++JF7BEqGEZffiUavXz)fnZIMYEuGLVaDcDYE)j3EeuU+e)VfJq%Og(8IAnIGS#4MC4#3eIj7K0lFJAt)VN2ts!aZDtLie6a3U*d7Gx3A&IMcEo(!&-n+x3gvc!tn#(GvS9Qy7prBB&-ZtIAGCzS5u5XgmcO0Y4VfBmXBlCP+f#6dyxOP2dPa1#l%)&f4H9GKmciCL!rRbg2n2l1tS%(5rpO7J1q6tSMiDhGzsQuO9e))OKj2tGHG7MMtz4!(OC%e%k+kvgkR+%$0Qnf77HmtrML1!y#0yb)bQVF(h&weEL3hUeyPyhxXk#q&Pc5ztiyevm(aXz8uQZe(dOVMn!wkcfw4YBGfu5NKzT(z!E1P7M!e61+FOEPPtmBGIUV37b+8UL+J%TTXZC5WUF)2BO%!rtg56O8!J%gakHRk7PJsC%u(-z6PrGVit&g*E$EuEYV!Ck3VGb6qPs%PoINDKCW($j8D245D(Suvgl7oMZwCEZoHxice5bxckGETwT1!GB9!Xm1SQP7+ahOe1L!beDOfV)b#ESPk1vu3rDrskhVa3G8dZXcAE-a)Kf#yeyFzePCrBOWH8KY!BjULdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*hTkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6I9LdOgOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-z1D4G8JbMbMeQhTlWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbN-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w)z1C4G7G7JaMePhTkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$uWnZr$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaLdPgSkVnYq$t*w