From d9abd3df0b3c0f0536d29c1b98b7e8923c77539c Mon Sep 17 00:00:00 2001 From: danielque Date: Mon, 23 Jul 2018 15:20:44 +0200 Subject: [PATCH] [DOC] update model options --- doc/Eqs/voidfractionModel_divided_pic1.png | Bin 0 -> 1757 bytes doc/Eqs/voidfractionModel_divided_pic2.png | Bin 0 -> 4812 bytes doc/Eqs/voidfractionModel_divided_pic3.png | Bin 0 -> 5124 bytes doc/Eqs/voidfractionModel_divided_pic4.png | Bin 0 -> 15150 bytes doc/Eqs/voidfractionModel_divided_pic5.png | Bin 0 -> 4217 bytes doc/Eqs/voidfractionModel_divided_pic6.png | Bin 0 -> 1945 bytes doc/dataExchangeModel_twoWayFiles.txt | 6 ++- doc/liggghtsCommandModel_execute.txt | 12 +++--- doc/liggghtsCommandModel_readLiggghtsData.txt | 37 ++++++++++++++++-- doc/liggghtsCommandModel_writeLiggghts.txt | 10 +++-- doc/momCoupleModel_implicitCouple.txt | 6 ++- doc/smoothingModel_constDiffSmoothing.txt | 2 + doc/voidFractionModel_GaussVoidFraction.txt | 4 +- ...dFractionModel_bigParticleVoidFraction.txt | 6 +-- doc/voidFractionModel_dividedVoidFraction.txt | 36 +++++++++++++++-- 15 files changed, 94 insertions(+), 25 deletions(-) create mode 100755 doc/Eqs/voidfractionModel_divided_pic1.png create mode 100755 doc/Eqs/voidfractionModel_divided_pic2.png create mode 100755 doc/Eqs/voidfractionModel_divided_pic3.png create mode 100755 doc/Eqs/voidfractionModel_divided_pic4.png create mode 100755 doc/Eqs/voidfractionModel_divided_pic5.png create mode 100755 doc/Eqs/voidfractionModel_divided_pic6.png diff --git a/doc/Eqs/voidfractionModel_divided_pic1.png b/doc/Eqs/voidfractionModel_divided_pic1.png new file mode 100755 index 0000000000000000000000000000000000000000..b455f0679be08052c569f2bccd3a8e6771fb73fb GIT binary patch literal 1757 zcmds&{X5f(8^=Gxe6uOfQm36n9`ZDzVR2+8^U!k5!h~oJtLA$MLu8XuYsaE_s(FYe z)>-zo@|0sG+ewInR!G(+N-fD$G>)@#o&Vzdx~}*2zFxoF*LB~&-1Javz-GvH2mk<^ z12Jfvy8G3^Z_)$+P{E%hAGLteZ~-Wwd1%MHTC|dU5BUN>%N6KaJXq~dLPPNW>i+-W zwwDc7y@_8s4jTsaMsFJENYEc6n=QOdpSfcDL0RjeYCOQlP}(g*A4FZvpGG95KY+6Mw8bRlSpbJgM92m&7nT_@Z*+Qc>x? z(8l!dmuB^H$@K}hpA=e-lia*LXM(Nz&W|0lGQnzXOa#>~H5eMSG0xO9wIqy$UNuYI z)30DdjMuAI&-LFx)tmqC8q#RYC)X`y%l*;yl8yko>WsnvwCy+A&^(h-4AY4?8$t}I zzgGcT81mrK#OqNfNyH6V%cdajEGv`^ip_dk0$$?k*3>S7bGeVcX*o)K-1+3CX<{G|N(j?zI>yprHtZgP?65x3Z)C%HGm zqN{mzSMjWR+5bvO5I;2{;&<^TE@#qgBR=1QSvPDea3ykFaCUX#E=)|*p+xtkkVp2a zsyf(m>nl^IXKc8I?U$BgToOaiMDr;^u=NM&qDimxF!~+s)Hm(b28Q1d>$t;M%|nE5 zkWg}UTISez?v}H}?L~Hc?Vw@#AF6L^s2;27{_PUZxiej{BPsb-_BUIO4HzYcL_s7XBM=LMC z8ORObhMFO!qk3fdR!dC#`Ye>zDf$R3kzA=KRdg_$q{KXnK8Qko%z&pwUs_o(9m>X? z&F%}Y8$uwcJR`iT;b8XuJTJrij`CUhc4%Wms~|P%*US9mPjBmP^N`)SMH9gW%?q$v z<|v(8c6J*%>3m7+g-^vE+lc1OYPs!~(L4Fe`E#pS`HvNvK+d8_C_2&tsBaMK zNY{itMfH8{MD|OK)9boU$wzxHKD^L?Nyc<~o}OiM1Fya5D&QA&?MmbiQ;kyNP=|?6 z(~pC%6;=7gLiBptGSKTlqFq+X+>h`63n?tyMWCOStm(Eoa=^7WQtu7}S+Hb(2SQ+m zuwcIV9mo)eG(lz!-pbes$~A6*gTzt#f(JS^BdLNdyq(bGqxJ%Gf=;ki3WQ286>%u5 zdRukf7Z2{XuKjDc1kTf{yfFaHWPS$EMTUSnV#v{!(HeQ#d%u8kO>}zPW5&jV$5)69 z-6uzx(YG`>^3TC9HS+o7i(JA9_@L_daQ&FOmy%K(aMbXsIJ~I)JZx?JA$a)Cb%BM{op=Rm%Pz5Gg*Y9`d)On7=m0pHV> zqdI)#o328Lo{X3xhdr1+{a8fR3om4>_kU%t+gMv$Zviz3L4D8f^jN@GQUCm*Kz}T{ I85MW#ZxjIJ=l}o! literal 0 HcmV?d00001 diff --git a/doc/Eqs/voidfractionModel_divided_pic2.png b/doc/Eqs/voidfractionModel_divided_pic2.png new file mode 100755 index 0000000000000000000000000000000000000000..9b331fb5ec10f764ea139676bd382da2c42b9e14 GIT binary patch literal 4812 zcmdUzI%46lvn@&09Q#-RvQ37>i?I)7%2ZzK;~@3KSA=;R*(kN zj8Px_6I5F%H7Nj~9)kU7iT2O?G*xxw{>}gINV-+M`FDeqm$teNpwARz<`s; z-RsN2o`1_Q-4u;H0RS?>{}$3O<_ie`fKX0JR!Y~`>^LvR_)wT2B)kFxYoD+#HkQ-n zOYE1Tt@~#UZ(HURft3Z*zo4h5x`1rp2&a|T^R z*t?1`WQzN}V&3d>v+O^xl79(Vs*A_5*1LeNusA0tmJFG3?)G)(Uq@|dj-9oHNoazh ztg3#*EoyK2-0WOpV;g{Rx;nGFo0{n_?Dmn?csv*9W&&dH$v%?LqjFxuS(tdpjwPr& zr0k6BmrToN$>)i-!+(4D_ch|{55BR zAk3-D0`FFRRY(qWycF)NpQ>3vhy_3sRaPch{ABmdp}{5-=hPQ0K)R!zPC=n!`PTF- zBLbB=UCr8vE6D^JUJ>C~^YI`ID<{4lM7jY8t`OGq+G#Nf(OdTS*lG}M^Gyjft+_}K zRL~yiLMH0DmBEyZ4GC0mnxFWS^)+nPS?q(vMHyMccsd-$p{qTG#@o~`MNfjXkRq$; zF{UpV8`8QuxEHnqiCoqc+ctKv4UWeTM&lPWX#AUSFZdTpVRrHHTLOfZ4)T-P*Yu*$-c;BJS@|#gJFWd zU%@aGcXL2)jSv>QMad@iXTT)=oepM?4C`|hv@>2fiktTYOmuVkB?7@tYmV8BSvXB~ zH9?k%#>%PWOTZt;x4*`+_29zvXvWM!(;`ald!YD=WNmk#B?_?O( z4yEFpT540ETuaZlkyJPv@wpwLtfnq(SV5%Rb^eZQc6Eghk&5*!Ef_}%#Z8K2;4{2^ z`*v%|4qQ2Wa?=1tR))|3v;q0iL9VJG3)G0p!S4e-n@2$6mQw3%`j?lOYifp;N~z$t?TRWmNV7E7B&W? zk7$RndJFT*eILsdH;G`*c2~T&u(_@LU`$Pj|pk1~ISrDMKV=+i+Z^gVm4zy~x z3ibI>o!c*tqwRVb4xTbnLr69hbC31#)PBZ$ZjO3C^CL0|*3TE^d*X>FmM;IZiD(zN z=h{w}tOq41>2ag}<^g>Y@Kw|+Xq2l zwCSGv6zX;(xq4_Jv*O_!XCEu}VJy0YQ-pD%dB!*^OaArorL)GeP;=#P`2TP#*a{VR z9NX!J)E_=IMu~BVZ+Q=8^-9g>^twD-`b##Osz6?@K`M#v&^-@tJ$)SIpMZXZWw}+UIV|RH z)|Zlrs+ybTER4NMH9WRyo+NudxgT-*a_v}o62D1Erk-i@gFIh`u%}khigvm=;k)`+ ziYu?CDUqde@Ix2vgxcK!12|C1q13+? zTN%#rUTK>?qM-b}!n^dOe$NAeM?R||H)MBWEPH648NCv>X9nZ>YxhF5i7sJJePkWc zNZ_kOJ$zuk1nbQ&D$3U_WyG8!Thibiqj!~pon;D+gPmgpN853Fd1BSKw$QROK zB~0$1Ts>HE5a0X;D64o5wug(vr1A2)q~ug3#TPvs=?@R~_n7Y~`#DMEbDAv&J>*v` z4$8)nI^5?mzz>F4S44;=RJ`lLJ8u_@maUu1OZ=f&XmF6WM+(e8>MB*Z(oJInFXx6k zGlxgeJy1$q928^bGyV=2!DS*`Z$z@g_ikGlZ$PeL)saEY$B zAe!m1`ao688`!62m1ZMHU~Rh5+h$1QYWH1)en|^X5l=FjSFzufoKFe#x4F(WUEz`O zTPH{4eXLVPUS`Z@y5RxMI|VcWL%ceN-tH~H3Go+P-zu0Qi^wVJ2!6zItJ}dgC_~C- zbhsEF-G7fGS2R9S_}5nX=xZ`xoBb$HGG9uRkUrwO0(=I{)<0b@Ix*Rzc7CsUpQuVS zib28$rKcw=FXoWJO-|JYS*NZRc=v{L66l9U`%kBviV+1m5KpOjQfN|uivKDpGt(;p zX<&F*DwU68M>O_2O#yeuBMV?Rm*SbzK~hxry2TtXpyqa#C#+B|P!BSSsQJ|H=s$kg zI5m9FKprzm(x;M(N|uYs{!D@6)dRsO6ehP&0;Wi5GVddy8GLu1qFtRv5wdIdrz{U6 zCKe7v5o&LZD3HF1QiB%bE|OuIdik`L&N2vxDebtW;W^`4m&1i@4&sLRR+2Dj^l z<)25(6S(AFGJUqKhH@q?#dbMo9po+K0=g+N_v%2&j~I25%xH;t@iJ8P9k0Zczl1&X zweW!O4=-;2M!)XMni&MoXgiKJ=@+lNsS3I!To z)xTQ5kDq)J&<5SAOVZfSg(kpDi=9pC(|D$K_e!!f;7IFmWYRy|W2 zF3-w`C+5VJtFZV4jCr4bE2frcsy$&xmlzw?HV-NPle6@qY-TkbQ9p1r8Sz`FGGWhx zU!)}TL;a;;Idj^hODtE%N#B{HR6gko8Ir#?)`e5*Lw6ioTxZ<>h3OPPwzPx>Hxjd= zSc98cW7afL`apMld))^awp$^2LouS~P?7jfLmIM8n)s7a|K11Ivfb;q&!d2;Q|dL> zON|Mn#61}L>X-dIVeb|ecO0e*SK*V9US@CvOt|^x|8__rnnGusZd~N`Z30& zaCs%VA>_Sjoa4<`)*tm-x&eo1=|Z=UFRElEEL)>glt$hn+x|ea1VAZ_x3ru5`n9nELL$J%aE_lwp{4V zZ-PmEcW~}Uw{9sA{@xE(wFd62w~P8)b}t_E!`UB;G++5|sB~$cFAl2t?7o^r=jro7 zA-Vfy#Hti(L~6uL3D_}rcf#KrUodB+G^FxdtLWxPw4|Bxz^P1Ouz5Qo+Xs+uNpSTV z4re+N9}U=_zHpGF=--`|2;z;;<#1TC?VHvxQ0s9W?GL1V>>^jPw(M>EhVty*RDpnJ zYgg@B)=q`UnKSDj>-%8+iAZXl*BJ7JFn!o04y4-{qDW&`NC{(FcGTIwRw+Xw_?*>{ zA1i0!{KIN*q~ErpJTn<7yIG*!jrlM2cVPWbGp7!&d;eCs1(|?W^z`$f%PDJttP7>a zReRgbD!$z6!3QH{wMV9DGj@=p@s1;u#l5&b1kBHW41wFLnlaA`iJ=bRWj*SgZ zEDq<7v-s6ko7v3J@@IBqKD>67?1;RZi+Awub5NpfpPMWAI*lvQNp@Qo#mBIdrxr3w zlMQ@XHlEjpg{_ zm~|sh5SA&Qkw~xN2YU<1rJo@M_$x?SC|XY`!d_<0dlgeAepyeG<^DLfn|yIiNX<~iKNHdX&O8nz~^*go`~ecT#8HN)Mw<+{1W=<p-ewRte=J+>{tBF@ZQssnI6(S+W|Pph zeSDC$SJv2iW^&M@>V!`7JP=wXT?6pa$7>tia1)@5>bE0gEFsz6_TBFI*A2IUvK&dI ztbY4_owQ5qzM~t*jS06zf;b%OvMoKGF~K%^Jj9+srD)ZFhexwxLbVg`X}XU?Hy1fM zI)`BFI34mF_QYL>CL)QPb&J{s!FG~l+5CviEx}3@ z%!jzV14!6@W9hg+P*XnN6KN=PsI&ew<-o>Pa4%AeV?Z?v3TgZP(aqI(mXW{O%CUuG zReZKH1T&mEU z?aalz*NOYH#U+`_tMIXll)rLHqj$_6dNp**BJZe~AW-0d3C*&yQp9DDL1rdU${ARe+;?~;1^=9%> z%}Jzenb>7lU16!?obq*+E)k!q)iqw|(bKYWd*oU<9LP!1P+I1jY>vo(2nEwh{uP%V z$i!?UNEThYE1s2>(bE2j-qQpXq;7-ijks{HM4NjDpr$P%uLW3^7cxHa-xX@@ewi>0 z`)rWm2wY;6&-lv$4~i8J47e$VKNWTC7T2RqZ_#tR ztV=-JGyCb4oj?X#&kh$zBS9j;LNb6$$mKg{1q_{&tt>bbH(%%BL5+O6UMZOKV}r<> zxtP~GBe%V$MZqihJ~JifqV7x{4;Vu8>|A67u`j_xiKN}cDYRz${RHYYhnr-td%bW1h$%~wAEs!b5B8mHg zxN*DZWbDGOBH_V=jc+)?skf*n)nX64l$axJPfvrpaa$)}<$uV5`-6XE!SNBsB<4@` zmxYTw6VeZ4FcjCzmOnse+dEwTJ*wG<(*twZ?$6k6be1$1ns1V|)ePlKy0!uu5wvDb z1R~84dch2ynZMj6iJ}Tte1Dh{&p7TPOicw|`{d0QmAMbzoUVG%7M-8{Y9i$QUxEOo a&nP-eiga-PqTv6&os{I%Wox7@!u|sjnC3tL literal 0 HcmV?d00001 diff --git a/doc/Eqs/voidfractionModel_divided_pic3.png b/doc/Eqs/voidfractionModel_divided_pic3.png new file mode 100755 index 0000000000000000000000000000000000000000..5430f275cfe65d7255d5e487699463cd3075ffb1 GIT binary patch literal 5124 zcmdUz)mIb%x5WqP4goh;~i~qC*4*EZp#{w<+XV4(J zDhdGD5bf?igW({rDGvbD#Npq+!~EyFwKepW|Mmax$bGD2{I`SrCtWQ)K({##(w1c0 zg=#+H6XSXulw&_l;4D|-u2imCwNbmn7^BykU>Hn39rBTOIog|bJ1L0gC{yIBQ2MD- zsZOI+|MDF1``<9!kE$jR0N@4be+#Xj{Y(}BAoEvKlGpdPJj}*5Thk)$LPhuI_yPC@Vv3(O_yBbLXBMXtEM#P7o~lVkYP zOXXQU&CNDzlZy1!zV*WbyK8Izv8`~~jhC8B2n*;3pRu2$Kt$7YEkrCanR_Pz`1Axs|y9EB(nw@4zVPka_NwgdA!+j z4nYPrJQL5$Xm5bm!wjt~#E&0bebSzuHNXN(wtNOYVVuHUu3N*JZx%E~~VuV1A?t;|MQ9b@ua zf(#$@;#;z2&}l-lZsN6ZUzXz+eFIt>jrOWJ=xVA1=ymz>yK86pDjQ&lC?uFja}u`F zayMn&jupz$rQ|8a?iAe4N z^UsEovl6xuIPy%^lxmdjr-}20{4&tf_$627o}gj5o>=_JIqr2V3Ott_hUBsbzt6Ej zrPjmVo!<`_>wI(obFJsnP!FX@#NIXD6_`gwT6kf{ z!mGhXCprf^{%0>c5q(NXc~*#c?Luz82ff*lj-#t{@~(Dc*P0QilYi z+(~n)r=SU=_eZ^m{DQs=X+-<R7z)XR>c}i$waluB`0EK1=ETlb>G;F=?81**j)o>HO7Jo z+f-6}iX9f+5{Ky)_H(L!ExdwlnGVIo@GJ44lJr}HskxcdFKsUjc$-7Fd|z8JC^Cso zb9G=S(YjYcRxQLSL$VUwy<2`;{7fTD);i!wv(vS~7~LIri87O`?DRMZmo<^BLKlla z$N163_%mbDHTtiA5)f~GtAnCZ?ujAUN{YkIRU~z0ZKD6WE2^D}ThsQD9L*%J1K6vN zV%-{ksoqgK%#eZ)i!E8j^;Lf3{7Kd3*<1#gM5?08u%s|ySTGAo%~=X=FHVmA(0Q`= zS|v&9`U-r(v;Rl+wp?F9q(j{Nrccxo*OsrcMgWVRemwq6VgackqKC}g~pwDCi87V-EAbWFPl5EZ-XoRg9pd1<0t z#&$<>Lz~RRca=(|b)3&C+OD>v7&AP)eT(ovR@ozGWD-qW49kZNhGeRSPbWI5;e=Y` zq5;WCRJ=#SY(VU?)e(uOnYXk9?iKuq@Dhn%I+mJjzbY<+6$zH&_G(xM!Xh%SFvKJN0_L$9|+C_oo5Vl4@Uk zD-vG(_F&X3M6-mVd?mv2Kde6m*{GhtpRO1iBM%;6@6s0Ki?E3ve!~ z<-=2edgg%UCYToj_fj8NU>`@X|7h0$a~hDzmqpkk;A{v$`*Fcbhd}{9NZ=e@ICEH zkogV6tMw`}cEIRh@#l_^1N8GB+}{-}L(t!VcG-Q>tw>}~PgiN4vbj4RM}x%GQjGNM zUg@pW9PBf<)Bea-vaP_Aarb3$ZDIZ*XEUSM%*)0vX+Pyw7!9c;d$oBqe7^frXL3LF zJZFzHdu^onHAn;JLq;2GX9UK2}BecR8u38vH1Pw=Xe_XEHQ+Nrl~w-=hZ{QLY5 ztQQoe6m#s}aj#JQP&Fr|j$T zOEH86l{&4e^`SI8iQ%eYdZJG(!hYPkN#4O)QG9F;JS0skJ0)T+f$j9vxS2i#{zjWv zRV2<3_Ez7tmSZ}99K^jr-}o9nvocj5wL4Ox;1`&Bx;i|(Y{!PaW?yS!L0E4pdL7u> zcIBmA?+dGMz!jyvznn7eqz}IHLV~8kFGKB z#1$p}5kli9ohMN}Tt8lq^-VZYUFlSavfd*+ewY*L+?Mq^P55rpbZYB0gH0hjx7{pL zv0PeZf<=q&a7X)X1h3h*BD3EUjS3vuKCg6)*g$c;OwOI^KRjU{J*&a5DLJ@Y%C_vS zb{3f{(C*cn-r;#ghvI1T{zUjvQP70%*+c(mB($TP{0^6#oDzmcU)u4+a3cAd@vZC9 z%T!yF7h3qZ1`Bkz-_V||StiaEpeTrQgxNSlc*S?8yomi#x<(=0tSL8&#mj`VS!h{h zwUt^5CT~d7aKfInJ4X)PLkMeA$tCEvl8DFNkIHhk-OM8gQ;2mdE4Z-$l+h6hx+#Lk0y1OKS@C?aSqrgbi*7R(+C{Y{ z7)hbV&-W~!i`eWGWs3lJCoH^3fxaw-TFHuJ@1wg}gtoQ#iX?SDsO`PFDrtCO<5MWcaaUf{n4p{w|5dlka#aZD^hxJgxuzxmT6V(fkk)-&7Ur@=fkzVs z>wPkLU@@0YQISWJtwfas#y8Yifg?nIVUZQZ<`jrO;-OQX6ZB0V z4)5-h8iAH64vq3>Y?$lz7?{2JvM6F|)rdDi7^st6Gr+?7w+lk{_Hln{O{WBaIyk3H zSG^CprpzSzIwkE^WVi*t>FDzV3)U&MJw-Uuz{K2P_O}0rDywrP z9bII@f}fgcj#1EBv}tOr0V@R)@L^1XvdW2 z>%mpdCK-CJjw=6r5taE})c@wz15=6)yiN zX_Xx+z?p%2_55M6=61=2M8lyxAGJd0TwFUn-#7#!z1#Hg=+u3=I`?o8^ZV}mq&saI z&r4$YHm#W5?+QZW#?{Y2^~jRK0_a zqT*s7WFz;eUyW%=G&EgE_iuGP>=cYp${ihAq5DFqPZm%0@6ut@3D8~o~nou3K<@Zm5-FIB!my1 zTmh>fMUITS3GZ#V-lUM9d_bt3L%YEDx^$6L(^7$acIjqR{Fl21eCxo*-s6_qtYr?Z zMI605YTwe|4Iiay)@gOFta&slG{c?S9!~p~z6RGDEEwdyZvs09zk5pp*D){cjEG+%KB|uvJ?GgWLR-<&7`%$ZQ}GYF`R_Pb1S*X z`PTurmI}LYsgcv2%X8mn-NOd97Sb=m)Hu$sjLH6tUD$vGcg`87cR#^2O^%O7NeRy= z*wvqQVgj@1uoYd$84#zvBI#81na{T!jiHm)1D8)Y6%*E3Fsk?t7f~#fT$SqJHKF;J zssXt{VMCuir^WMqTy9?d)*9cVc31AUB4fks(h`__L4?UJmPKo;20gTKkcJqTJu z*=VU6jqJ$H%X6mNglQ8Y_#Uz|Gb<_zlbAlee<&|L{`1*bDECsj+W+}jncwJgCor^y z?2iqoTv%ueK`BqxPUKMv?ECw-&j%8p^gLl)`%eFeUm!h*&h0uk{#5->v-z`6t%^Z~ z7)H%Bjv}uQ@z4eQ{atwj!XDn`xyz!*KrrtvI1$lF@&`uDvR}ixi?64Su$I~QU!c?= zoIg+2CRn^c<8=eweVOkZoYQCl?WQ-hv2q?B?lwMNYSTyTvrd1ZvS&69gomEawfKhw zU|i?Utc8_b>Cw+PzdsP)&(T?h+%gWDV~IHQG|Mv1=z_Y;x>Ubg9vWb#R$n|hi6IQ30Ku0TeB4(Yh4v_&t#_*r{bJI zR}g|WiAr0YlFKdA8BuSoQA_*V6`fYa0PNmJauk31PD|UID8_W9#$xteHF;eQ^BD*~ zA?Yz#v9heX(JEWFDoC|}6S0}8DoOj)ge_4#>+Y)K^f_RCa#iw>1)!%ZVPSupXB^k~ zSFAg2sxEpxijt%U_u(=tC@MdrujhAn_4wErpNF%Rl!dkh#YKlPjs7cc%~#%B0!0;p z^?>?3ijWVU##qcW9_w49TBao~PJ_wQ^FK5-`#0UmcIH9j$PVn&(=tAImteMU1m*v;dL!P&0IWO5ts6SdC9k{Ix+45WM5PEAV95llF$yY_hUHO1zf z=$DgQAB4)G>sF7@K$$H~Wg+Cqs654wqueuGoezY{4@gVc$2hJV7dvx zTTcrXu@+UdRbX5fv8|y#-{oUl;B!c;Oy?IukHAx))I0BiIQDLeR8KpX7xh!Nm=L!A hy)jtUWn=6fh&#glTB@?|r|*o$Sf1mF&!9XJz&e zvnE_#Rtz2n8wLmn2wp;5SP=*a^z2`J1_cTP1k9a)u>CIqc2pD-1ge?AJ^z=0nF+`U z00GskZCfn@r6J`DRqBBZ=F|1N_DGbf&2+(=w~mfB_8@p_EzCd z{uZ_on(0YLd{kQ5qzA54-uG@n(fG`E5+3Y{cYRzLT-7#{&?4>hzyQ;R0>J5T+uO~# z*6FSj$_DzNJGXBaLw-KO*+yi`iZqEHNZ=3NCh_6V->mvAN72@RCxYl#VjSFwx`gG< z6sBGBke;6_#yq+t`Y{$w_Zq!jtumJ{$p5}Av)fm;>8Xnjm3Q#=nrYo;&jC~u6CfTR z)|QmC?w7r@j{@rg@pc?%#R8a;sy9XUfQUPFpYZo7lBT>_#_aUQB4)o3$Tt>C*sP?g}A=W~tsp6YARWP>ZNr>$!8+bXD+A`!6!o{H@ zkz2t#BoJ&$Qv;lRVn;fUEYiJ=_8M5?cE|jDO7FmaKS$a?9hx+4sB*rBuVEor4Oo2< zZG-M?TPd~uT*Nw28sojp-kqhLk(L#xTmkw)AL73y2&e8YVOZNJ@vLcOr%86es@?IR zG4W;TSLXYx7&GhOJA7v0TH!Q)2Jf(ZDkiX7;me8?+orbYF~qRte0il+Se!l*$lsfa zpHPCi?}QZ_XZS4{#)ErldXP)?VN~ZM>hBX`)eqNydzDuqyWax`%*uiBp`C@p3zqTS zaz%tawDmTRFEbZwb~~n3$=H{2t;iZy(-`!&yzzx7)XZ9lHMEHo9~wa8evZ0!wVuk~ zyUpglz;@Rp0;A10>u45_%QmMkDL=5=3-&_q_5d9ZYqKXWg0|d{bo=UOCqSoEOavOthY#Vo2t%Xiy0$aJt zPt50AZ8-Yi@NI4v_2kBItD;TNZCjaL77M#KmV4hXt&6g5$cIOUC)$`3o{s10qU1l_ zBU8Y?oyf%BooY=E^6BTwf6zWJ`@Ltu(2K?bM9hOiv3o>^fol<=Opv|@=d8Af&OpMv ze#e5xFREFf)32h+CpD5cIpo^HTNHl^gRd3GYo#kqaMf$6T|-D;--w{+hl+SInojVZl+J3klZ- z)B=;6nQ_(*Mu*im!~d*Uq7j-;}s60hoEB+-i8lamq+m| z`>LM`5bF6}YLv|v{xP2n)RFXB3)^({JF3R!_eWx%9+2_y^+)x0t}V}ti@A?AOX$xv zA4*0;LU)iDx`6FVEYEf|e5K^IaG&+LLxxT8-Jq`r3+Yq{sI18;h$j=Y4UI1glB?cc z>H8P^L1ibt?~Zyi+=cUsP>kp^aYts77SQWkkU4TO)pno{(N(b=f9}U*J|e_?GRYbW zS4qF9P1f0w-mcu9T&dnI8wBkz3*b%*%LW~hbl}zJUESv0NdqLtKBuzbl!jh)R`zk{ zm%(@r>xHB~bCN1QK2<>|$V6fnVN;XtFT@386xc%Jl4vQI_YE7O1`WMExESm}$j^1! z#RifgB3nHcOz*-qyk4K&7gGlazYrh|Ui!5-Wr;0(cn!A71 z7SIb5l;R1!L%c+P;210R2eU*a$OzqLUQoln1#{E{eLHjtGFnV+3qr*RhT*}5^U^Dy z5OcqqmdD9d_C3d;ICHhgfIhm^K0|wFi$eHY;!V4QzJ`!?fNC@Twtq!&|OEM8$ry_)ewm}&&N8VJ?p9h=4 zx95~A;=E%WBX;5QYmkl{C)IMq1om&-Vo;YlcuPQtLHKT zr&;a0_gD0+O~|lSH+Rr~J%AH4Y4qk6)IK72~LqU|uZyj@ALyV;Fe-qELq zS7JjNdWHDq{FT_+BXl5+?9x442f84&z~^OOk1sQ(oFT_6?k8HsEVmS^6VT7_*$>+C z0J68(=S;#Uu1R>tp6ISB)NQTI{O%SNtM84E_U3komi@d>7P>@v!FR_P;rc`T9o>Y{ zX&cv%cxkSqc9v8fz=J%usrG(}F4$hT{wjmy%~JkX1a>RoB=Wj60PT5Z`OGZ9?QyIt z6CV15EJ#LgDF#%6o1*R~WOFlGYQKUFBB{mx=@U;o(m@&&k+yG*rbRo}Y)7f`>lzkD zLEq$iIXC8W>>PZM0dA0CsutK}x=%Xc4V@I5n0i7&;er0J-=2L)@`Zpv?f3q3=9DTw zXMV(^P~T5>YE|}+G*#r1@4PgY>TBB37B7V@PkVLf1;ayd-<>BNM>+z@iwV{p!mNI| zP@|>EgN){r-?K#GLcf%GP+2D@!Ec5*S4QZRx+HzYZ{p$eo*y*bEA zOmB1QwJ|8x5h5$|7VahD3Jff`vo|}ysAehJ5RkXhL!Ikldf)A~0papn!UJ|Ee4Q;_ zFTN|m7<0_BIU^gYR9y1JQJ;8DMm5Yz=w<987R8*Fk-LOAf9>=E`5pqQ%_xefe_sDFNE>ZLXh={`Jqnz*;6e+ zQ8w7)hs2m`eRJz5?m*k(ojo~p&(9fGoRfZM2%_o~aVDFiO1opm}iTP*sj0$*mWE8NBy@>0- zZ8S6@WSY*}l!6^KZpoMyaW`x2+I~K(ih&PZbE~Z-xL2-OUfO%#pRL5p_U)ZUSlx96Vn#9+>Lnnq`W`_Xtsg; z11R=N^gl)jB!N{RF3 zfJaEPUDax)I9<$%eX6vKkFiS}?i=C$94}06H@1zQI`kK-+xT^`;FB9_#!XTN>)a8` zF>{S9X3w>b&_SDHA$a1HZ?B5{9-U*46G?}v@`|}cau-=u*R=V-9Uw3e6wI+pzSh2q z3lxcVMIaiBQKSrngL57a7C(46RwB>Zx5Q zW@VE~`9&E|a!FmaoIymLK`QTXSbV%(Vg1@6e$=YZ%ReSr9-`aF_jB3W&vx{A^O%~A zj6mD2Nh{%7t|_y;1>}_1DcZ8>H4! zhZ2J0u_|6A_6|eOG2cK15mD(f6{bzSu0hB|B4vLs3UEbwAUiTet4*c*xtXe59Wi2e z#Rv|M*|IRu&!myHAj?_Tj8Vb)GrVTj>k(GyC{9rOuFa#O-;yX~gS=UPowfZ3VMt?6 zwuDHf3vSR0gT7gK5~D5jZJ2bj=&k_SjlzYy+s)<&ww z(!%9wfwqw`Yh=rswHVLpvWOId!fbP7D)@JLxn4Heu>a_XJ#_W(y==jA+o!W>VUYHx z&}=;x!gFU5SP_Pb+!2zHFG6IIpUHgb>1%s0?b-UT{c4)vQEGpfG3T0NP%j0_I~gR% zJhf$^F`%$2ZGl8xv#A%PgU!km1{*>HD0X)=*B)m^L>47_<3K*$>=vttkz$hiD)T!Z1`1@nDOUgZRvTpPZ&zH^FB1 zl!$vvc{a6B>aR@&mb-8d0H_^UH_k+Zq%X6Km`oBxmLTaR&{1kd2(-J}_~BSDCC|>` ze*3C_Ot^^Edfb=uQ9U;7f*e+{AUV@G-A{1~>uFj}RapK)nevd+0wLZoXuB;TW z6M+KR0v$Xu1=>O2uQw|ul5j_DHO{DkD91GPsQtFr^ivN zg-^D+))70Uz7N8{&HiTA`GrNPOLmYW#!AEavY4E3F;Mv%PL9zGU-^)$TR~Pnwrliq zd*y1zT|qfkCKd47oHM>!*=c3{cGUh&3il_xFm1F|RvV#^oXvw3) z&b>GrUX((ZN6G%idrmUZnk~gZ`cti_im=Rb{K!;!#%CrsvKy-I^GQAP({XX6C@3?K z5UZarC>;|ettc?ttDx2oN0u!pJPl&MRBA^#P#Za1RL8V#`_ZT!4jhVo$*F3!2)#%i zKJqzzZyHDj2!g1;W&!034NU*={(r`%Y~z{CUH_*KKDOb60j&=ubG&Oi3A-QCDp5;$ zc>9>0K3?jZU;<%7Lx2^p)~S&r`yXavbo(H`ioYhSCM{vC+F|#YR#~{{+w&%s_Ovkz zmh~&$;=85hxla0XkHX(Q)^f@NawnyWU3|Cl(Jw;}94jUpE(+ZoG%CYZGFyo%*H`_z zRkV6|OTv~E!X+o|7(-ejMH~@F7)q5rB=@ZqtqR`nF5=5DhQ{C&3G`hu$1O~#F3uB& zkW{Mid{<7J%y2XhN*v(##t(m~WQ>YdL2hk`pq642&#mes)6+5ZA0hU7Z00f^>SazidNDLHTUsTv^UsP)MmG-1ZYk#nwqF z`;L@IFItCTTUT38TnO^5UgN8rY_zsnZ0r|isqeKemfxFp>|mr{Q6D?HD<=CS`?01B z0SJD6QIJt>6m<5=73ey&R>|shqOm!y^Pxu%X-|gQz)0jp?&&O9O1-H52AjOx1Z^i7 zi+VPomz*ikO&9HAcoc;c%K(w|ECg=n)1aFW97B_PY^4{ykfzQ1AwdH*KJe;*#b7H? zi|?6Cegb>&uV`*W zG|QB-8{4FmB)G(UIzRhHqQ<*L)M)yPW$cI&5vty{JhQ*jd{j#`{PRl2m}D#Y73*n#0l8|mhVKCs zAC#;7z}jGRl$^isamdBE)GONnUlyA29YkGZnRuu)SdFAwtlL)zDTzR7voH#_D12Ic zGCthz<%gLiVNfk{?FxFsulT9@#nkI;uqHDk9fyJ;9aQp2b`BgYH8)|kMjhimm~WG{ z_us)rAu%!-G}--uO`Zw z4F6Jt1Z(tZ7Y3ehR5&7B!-~zh{6pnqS~KPXX?R{tGHdWK88=;BFnDb^EnbKXsORGaEbh%9I(?pq}E^E@XH)TdzC0m>&gHFH z=x?Fci`HDi)~r;ZzCh#ke)jSpshNG^eQ3r#@Aq;q&TU)m`0*c8sq|!LUt=AXVn+&` z%e#K{UyLK2k^`?F71&{e*4^IhFpm!~b}#lM^QPe{y0@F$C1unKW~^>|8&AugVK{S) z>j)}gUyw_a(Vrra&zkIx=%%^y@PhDmL;IU#(k$6t(}1kO%u>3X;P%^<_1@7Ql8_?* z)9tU^w)xG+NtZ`4Les(2@AO#vI*;r1@<%+9+`bfXm-+GINrO;VXN^N<#$h)830*TT z8C-G55z=u~({WWKp3lwUi4h!9wZB`iulde1WnWKg&T)MTS<<89efGg330`uY0LHlR z20};U;NNld+O9)1ltU)#cUp;~l#40NcJmUE3f8?r=iGT(DMs;_a@F)mib2qZ%l;L7 zpOyZY^s~*CHc?!0gjQ(_fo=8DCQGnSoYj@M-0OI~6Zil^HGL2(y3-s>MYAye-`5X> zvsp5*ix&$4B58yxDzf^4iM)Akc&WPXmmxC-_l!838Ceb}v3I$V^<{nC0jFI9x< zO<|{*2kOer3zs#37QQIwjQgTXOnwMj}T3V8QA(g!|Qku!bIR}=cI4&1|DSiFT`s$BVW*D52y&ldx zY3{$_RzyOr?p~>CS*j#eDyNEs!fyN6*&*w9M+ETsq}<1ZBEfj>#k&G_D9IFln~uJM zkKcuW$r36Y0*Fh%Ir`9COmYix_0BRYiHofctG%YF`c2L%sTb;gUKxm##B31prr&$4 zi8%S)F6v(;`o~4B|blkxi*qSZS`nx`{5KI#)KOBK4(5fY7)&IseRS6JY_;f9m{8 z%ve8971-7_XpJOP(|TUBlWQ!DZkx08{Zd`T^6I4t*qE9oQkY-AK)ilU8ZkLN1r)zC z{W8tY3dL(Sezr0D4aBm6n4eUU7yav>O&e(Hcw#Y9td$W3D-K0Hw1+(g#jg9Ab3Ju3 z4m^GMVzTTBX_uG6H@v^8rEp5<3E^)47GML_Dj;Qx4W4!jbO|7aVRzx1-ppy9@5Er7 zIT%_`YfoKAw!HqRVp%s`n9UxQMbJ?^DCJ%(jTb%EpfjD6I#5+QidiapIs05|vg>s&BvMN| zsA}jj;T(+sv@)d%$=H6X1iJTcK~;%St%L{TUwcuxK3m!CSd+y(OQ{l`K<}B`l3gZR z&0NEM1*1xcA1)dlB1q>obk7}E#(YRb2@&RfFh+)lgJ0S7#){{Lv0Ahm()_wp9n4F& zZ35_Rw~cH%niej%7J3uoKB)#0I~SgH(x#C# z&poIm^o3H((UZq8f|K=R5~hw}Li6)U(bKgN9ISg^ZGx0NMoF6EGu2a0F3qGsN96BG z9Yo!WDT>FMkHs$X;bX~;N6}h{C%uq~(vzNfOfT%|BsH+}NR(Qk6g1B`X`;!)$>w30 z4-DRas~roFMvX_ooOM-Cp#l zf6(I?W3Tck!}K62$M)fe?82@rBXul1eF&ZTfh}yg|4OkJ6BwmGRj6gs>jThIGE>6y zmY^4(bk^c-Rb51tsk&Q@RfSKv))#8FRRtmzadG)~m3Ki@=LHTzwT5`3^ThoVa^ngM zhXZD+!j8dTBx>)ij>e;YucR^f0Ed<+liX(57{A<&(uC;r$23chMj&VHu%d*(!l*RV@yqT{x3l!WxWZm?TIa-^Wvux3TAIP?s!t;#!^oPsO;DRceVv@-0dBBK=Eu`q&3x%BYnyX1NfCFClmwyak&ZzjE3@8I_Z znTM&-*al~Uz@3(6clnuic%)9Z@2TAwuit>9XMDmlifUH;I-w<_^3!DwCprGr)5@`m zGEYkaC$6b5_%RhzvKnL;Y3-fqfCl@_c-um6ELUj>ysv$4LN?&otK39}s&=VVP#XIr z;Lm=yM48s&Yne`#JZH#yO-s2^PoTGF%@=!iMTmDJ^~jr9Ld>xeK=DM5-4D3^qlKB^yAdhy@z9)`Z{~A06UOIP|@!bOUqKWZYq(u=Dd8GHE1HMuH zLUnm4UiXq0?Q|NyyDHSio}RM8eXM12p{_|!7?SgOu$W~b%OqC$1tq;1>2%2z9!;@^ za~jsr$d=3Rln^~ zVFGfOyGl#_>KL)Uo=7}(N@gJ@k!0X@17TH~&FXq3m-k{t;+f_#xCj_l?!28-d{sQJ_sZ9E?XEzrTnBKqwi~4DImB9 z8z;SJ41o175i(ot$Mgx!R~W#3Wa<4dGA{gu9V&v&ID?*?UM10Z^;T3!1e@PHhqzj2nz;~K#fl6jzk z-bx%tQH+ZcK1;rh51#YiswzJZ z^5qz$HMc`fKnXd?U`8szv_C@3-j5+GM~?C=vO$!dplRqT1}k}4XHcy;Pp;YIuuA)O zM8V2yh{j{6pYaqmaf5=u|4=zzo85et_S@%1e}XX-g9~PFXsf-)|)@jmeZYavO;L!DCRZ{Hzv?Ex9PV z__k$R(ZDptQo%g$K6Y}YT=X>!DeR>@Om~sq${q?cxGAaWIvyLK=O5yP9i#!a zyB+hNn&d^Z2&9#0UB+NowK`d<42o`@#PrWf8!bVsR1!3O|MJcv-;#p^$~chQmk+CC zZgu-qrXx_rSE-$+tw@KrL3EAMMf_kZKU!?6nHG`SLLjHNbGM?QgWoTCevHz|pK)pH zkcZ)i49W?{Xk(qJqdr@i=6FTNf@=)Gt7X8OF7N5AJ1!5&$w+>Z!OcBcSZg-e8vgRZ zR-$@H4OHwyh5I;dTUDTazJH$Ib`lYF^(7;dQKy3QV?wu@-V@*@*Mv1M(B?_N9#)cv zmeezWBpwP%L5gP&h5=%Vys-z#`;EN=LYu^PZ~q*1>V62-JC%X!7Hl!kIdEr0*a5p$ zMhey{bGAD95Uv_wFiWVSt@g0<{xI@{!#M$DwY|*rHf<@}r@kSzKWfU=@^}A_Q?m`{ zFW&L+E@|{-_?*hhw+d%Oim1uoIF)e{4&d?2*PD9$F31WTxR|DtYkAx7-JYq)c` zT2)=Y7@<^N^wt+^ec00|qgPj>9AJRyiuOWb-;sSE5k#5TO2Y1H9uS82Ms#^>e8L=Q zwrnlmVB)uB49tI+CI;f6)|FhgLtdl|2_x?r8RE3~rdccz!cti_;3Y zRLt9N{U*?n)R}2g_pmK%h8)4UF;@}=H%yDR zx6m-M$-qsEEcw+HYp!i89@!PX;1bN_Ut)5ms|@nly@DC!05hU9yR>qRQS&x>y4m^0 zP9O;|;;4sV{09#|Okc80X5$}5<$*hFfUs3i)`+d(z(v@fh}1(8U6qO6s&dslrcklK zLq8RBHTP<%cMn=3N3_v^5Cb2^xazIcU4y54gPA=L@@??>f^R4LO?-R5NfYKh6Pqoz$Y2x;k$xJjDW>`(e!qN1j#~8+IsEb*Y|T{grIM>d z@v(>L6Bfizh!7sMjp^PS(`EbK(5L?8Uc6K!@x_-pVW1Z zN%vI6zYuRwfbtzBDVkJLBT>+s(c0z)nsNIvZfvJYiKxZ1__fU4)tkS1K@fxu9of(RSi_r*gzRXbgw%U?yH7jw-<}P*>~b;v;*gA>GS0FDv{5Z zJs|sQ$TEdr%ei^!1S6Gv8#qm~zMMKTI(H^0Q7KpbzKeEA1J{U+e!Xl^Wp+YB&4|hr z+IKo|~+dUOCa)6)%PPTDkDBjek_lzJ_5^m=5GMTFqP6F$ebRzHDrx z-j)Y^#M{qZjob7Nw$3o#9Yd(Mf(o`f$}N5CWls>T!69uJ=SsZi*(bx`@LI8Eb*fu6 zrrwg5Ha~BUPh7t+K7OHda9^of6^DI6)T*T9#bUXhBVF?-&GWb&DR!hwFftf5$=RhT ziU%Zp7E7*>-YQCv;u_Z_tRR_xMQe{yW;rQ`8|H3`+2zDry}~)Pd;1$N1}C|3yF7;= z2hK)h!&9R-zYwlYRVem+(f=znkt&~g`RsHINjieN7J%JhX;70SeIRUVb3SvN=UEYO zp&NLC?U*wDk0heDO-f{NB1-Skpke4CzQaA>+u9hVO((n(!Sise8c>p8sVUgzn#21o z&MmJj&jBn2B2n2aZT5sA+PWdTTkMOHqX94mj;bi$LZy}PS;2H6TqnGPRRYwZFMG|6 z7tIRxE#)r<0q{H%#`I|R#UhcWVBae347TJQYY8Dfgj^{#l#u5!mW835Y$UeHK18zl z89jfDlwPd;8|7DAe1)w;!X=>hnr%klrd}07?i$iM>@DVF5$!m6v%Qe4n)Nvv(Lvvp zp(jfjQVGm0cB`r@^Iac`Y6#-F_%Gp}i_iNbrGZks5El9I zg(D1BjTU#^^?)tbZtIR$+A3A6f!cixRG>wuDO3XRpOYx?rlI7S)Q_gY(JP<+t~*Zb zoQ#dMRA&|dj(t5*172G_U{Z3gN{kr2?t{pl6IF-(g! zm=(IKDGo)WOipE14!KTsTu>`qt290px745b?L}FO-{3FYOY%^!Vbk?tuBXjjuwU2squF0*4sosld8{rpA z+KJ-|)2d+S5Ch(HCT=Ro<;}Mm>#dg)0m;YIDma=GujHjLYHv8dj1Y-F8W7z~H zCdR3kIZ2#j-loZsX6$4!hu3TRac?@=snEuD2QCY&qnr#EjTZ2Eixjf0IjoaT@phFt z6Ko$?PTMI-g!1s%X_M^$NZT{G2OZotmAP7Z3hmzh@tWNo-Xw@eF|*?PqcW^PF)rsW zvl*@zJu}L_|MyKfSFFTWDRi)P`<$ED6?TxJH{g$SQ!&iHAo8Ot6 z-N}R!VOYoQD{)B~4-fIL55qyy2=hb>GYG=nBa;^K(VoMnVqlKo@?24bZ;JvE03t1% z$h;D5G0f2G-x!w{QT|z^S^eCNybHPo(Gdo!7PmAVm&x)EWq;twx8~A&$L0`Nuw_bY zrt;||lLY*T8&$@!idE_ye>cZW?LQu1>5+vdjyj|Vh~K8Yax3Oyh+1h^3U_9s{@j?*Ekm%6&ya#!lh(z=PfO=n8^HLYajB6GGUs&*g898Lbw-P(G}cje+0 zpQ=WTF7WdERu3V;J>c|(;1_P^(crgN_*f`5jc$G}g5iAndB9s|f4A8J(h`)jvpuxO z^!jLgjYYd*e@~vn2|@cCl3H=#(kV>q^?*Fu%im%?sI^D=uzr&3xgpZ?#~a$rsGd+@ zy#rpv*{jW%bwcw0_&H$lG{1yQUs7fJ%*ub3!)@6#8j|~z%}%Lf`&{UA9Yo6+ z9$Oox>5T^n*8Og5hT-N$q>JD=yOaL=1j0sq*->i(hB{kCnDrAVogwVGEO?#9fP7l`0TgLt=lE*X*FFl z0%-(OoK_?FQs>VFi?wf-1np@Xja3;;CSNx%OBOXw%zrPIqb*0-LD-uNyf=6?xQWUa z(5>Cs_?Ui!PU`s0F<{&$!ErxDVCy(OYY~wwFc|A(@x1-Yz{gm~H?i^I=<)+&l5S^p zj-K#isX``T#p%?24eyV!(ljPRuV;tMp(q{KQPEfSzzp>v6}<&qe>VIZ4541-G6YXr zHDxQqROf3r{8xE-c>OO-(ZULKu38I))}^Y#-qg>)p2nWsz?KchB!K?3!->Rn)v|z` zJCW$g17I|9p&M*t!W<^W8FA1$YNVI7Z`|q7_Lwh1b+{78-mPg0?zxNl9(q%GcCd%d z(bA!y8%S{Su72i|a;CVde^rC#
j3@z4rhX4^j~<>k*~N|zl|^!@|ZO>W5!(a znyHx_zbYiAdLjKF*QmA57ogWK@ivt1q-UC6=8En;r=#KKE+-@hx>UJ543H4%)M7a|n8<+IS{v1XH-rHF)8{~kN|8fdCbW)wwy0NW#m_d7Jrkz>d@X{$qoyB(L^54S@*qD6&b zO|jLqIOGtv8s+nHbZz@{6NcS(wanxcZEA0a2S@rcXQ%D1BKJ~rDlakxeSe58eJx8f zyk5M@%~Z*_@gk%=4R_bH&skYO;fZ<_>Mg4m)_?)iI`CqBy^ahU@O8~uj=$yyHdQ`f zAPN1%xf)L#>gSShI7N_)ORc2B)_|I>!y6pK(l&zZ!E&tRO#W#$5X01I{CIi!{$8hD z#sLrHNA%*eqT0>m$0mGP^p3r_wo9h#?c#b#&P`9~ zmeg{asWZj{28Wa2)KwdKLF_yuTB*GYCEd#x_(-8f!r*3hj<5|hL{bz>tHbSzy1@)-|nMqpj5JLz`XVR?+z?_H0G*30Ovopt_KGSb? zTsxc*{3dnVpW_!5Dbhu%U&aCq&>!#K*#uzFKOtrl+RI+j#O&?)%VY z`U3X=+wL^XuOUdr`QVP?gmeO0*?)-CONH$Tmd@IjZ3UNeR@^JL9w==@mJ9#@mo&k$ zTN50%CWT_2$R#*->WMxeZI4u;A+d|B4w_lr#Kb=WyJ$=Kx^3jQn{r20P3~mja6Y8= zhu~Xa9PY(-Uj1eXJuPxePcmKGgFVV*^m~JF+5~ko6J!Htb2_09JOo4;+!wSDjTJd( z)9>z2y*u^#(n9W8J_yM$B|UaXXnxuYnSC=^6jS!{WRwVBo}ZkVQI*p_IQaV@Y2U}f z&dj@mXTK$qO7w1*K$&c>-Z5rn@q~+*j}6|Q%7)litL|F-Yg^s@2L8|vSy$y`=$ypf zFm*m^05U5@40#hN(w^5ec4Pp{sZ#I5F-cx(1fPFL{R+WZc{^Z?SJ$CyELBH@cqjG@iQ}i1n9QC zG33HD?ggBC?Z>x_4i(*gD}7w3{HIc*LGOdq@MgHI4Q=ui)y+&%E*y(Qdf|LfO#M55R#^Nwm#!lQ`&G zj;%dqzz;aj``Ruy!}mpL{hNGJSji5{noA+T0;&3G;yz9UW#35f*Kyd5{UJ+z^)e|? zZlrxA_nA>qqhpublJD=7ys#i~mzCaphB?8LYcAGI{RHCb)P+}wY$F(SHy|{Cr>0MZ zhr;)&vm)SRk~>l5JgGbuegrJyN0%JcB~w|;MNC4&kps>^i_jw=+Uor&q-@XfeXIN_ z-u3bMLZ9W^jsx-dL0A142HkG!f)$!~yaFZaHeO(7Z22so}VI@v+ zJOqm93DZ6e=n#TJk(E7}{FgPG@46LHZp-s#yCUry0~n9j$6YfrVX=PRvbKQ-7)y)jp;CU!0H;Dvyc3UQr{CiQ79C1$v_T1_P>0 zQN#XTMUhWf^0^=GJ}t4$Br$$VXK>7tmMbM)LqCj*DQ`f)@%Lf#qGn^^^FOy)9mL!Y zhLElR4$dL@Rc4hM0a>oDFgllq57V2Xi+dGcd?WZ>z83I9tLI+vdPl7k*p;@xi^oD| zTkBrnIPhe@6meQq4(aH0ZIK)!7c~i-$$36)jdJ28<_S806htjGBaqgfy5f{w^d!cF zPf!GKzf3fWbfMe7X*`GJ(5|ua@R(e`HY*}HexUp|SxwE$mL()5+vq(*nFoKU8GTZ% z)Tkl`LvFMOWFI~3+S1$>x3XY}Il6;p1>-2YcV@U$#I})#SDJ=~hUuttC?2W%{w^9C z{G$a@ed^ix~J3@Hnpm zB>(1V1TW*GbR7`mRpU6E5-EZCV7Ni;sYztk2p7 zw)o%N12b{TC%)C!j3OTQ*!JG)82`@7ErumII~gzOfjF5b9S2o+>q#LPr* zTwpTAnJcLx6?**u6N&WA2dUt^Aq751f`_rTGGV!gDv#ol@&Z1}d*3X?Ec>n60Cfwxf)f2h8S6B;Q^>L=?8a~ z>1DuuAwt3LS3ZWYI2xJ_x9Ts0v0np|Zz_0`1>;Q|^UL(}J|;lUh3%!R`u4^pUJmgx z?}Ws+#i|X~>5Oe@+QaaebU{lNW6&E#-cySHF}{ve_k@FBy~rV*HQL@RG{8-f+ptI= zT#2UPw%E8hBRIEfknhW=(998`=RaMlaaE3S`2I9Tq&qapR+h)Hdf?(q z2+<22Jej|AbU@H;-kAN!!v8S?u_K}=GAy5!Lc58b4&mP5b`Rai0;psXng5_9iF&Xz z#ZX(Rcuesc1!|?;fTjg(xN(c)mqP>>)DTzytodV11YJvbIkpSpvYX|UgRdDVShZdO z4n&YH``HJr^Yd;Adat;b<;N=g1d~aa!|pXx|DFgPJ_u%4{}K%^(Tfow^W@I7ju$Rj z_oEi+M7>xSCPBFAXsS{Qmqor|VieAbcW2MN{9M1Ay!DJ0ZxZ&&`NYzZyC7b)IUiZ? zmAB5N^+9$m>--_9&A@0c6>3v!oAQBS!q)2|7^TEZGoxM86K&nyq zLTuz_H5rO~s^@;bEc$r6$U3E~YyB4V5x6_HTpC1OAUpgexSYhL3Xxy?e7&z&MCA5i zZ(1@nm84Z6o`S*pRjc+D;Hz2NYI04BpOOr~4~2Nr$x#+HRkpewgaga;QA^@oHf79c?m?oxF;M{+73<*Fdokr)uxR?boSC^-9@T(FOF%JDdV$s4}}r2s6XVC#=I-6R)3};HrBr(Dzin zTC1Jdb6d!riyt=kkZ7u?los5jXcS{9{i%In<{dcxmdb|!n(rw~5@d35bW$jznG_u5B6zT^Aw_(Rgzhsv^Jn3H>Wt>^!DA<7ie@}AQ&X{K^w!~XkTX#;~8;f<=Hqf0i#CKu(@?U ztn%kBKg%esfF?6*mxPAQ;Y*6R86Jx&BvD$v zq`d|B%|IMo0@bfam*`c0YtGc<8G8Ihvei5HPk`U_s`@KSj1u=#4X<#UY%AtzRk6Ry zB^6$0f@40N3R&?P!h2w558IczLH~Q<|JN=p43`^MO(sX)Psykq)w(i>C3?F8*?TMR zo#;%8hyAwF@HOwBeKJy*`Da{RPgXL&lZRbLSyE%Fi*!IsxRFPZTBDwAK%9vam?+LJ z|Kam9MA;6hvpMiMie+9(H!wQNJG(*D-jpn7MyIh^`r1}miO!jH_d%?U|BZ9k?u|%i zU|3jT2pp~-jH46PUjJM zvcqyD`>~3a=}Ws;*8IB00aO{%RJX!7$qHA8uNfbxDeDpL`tChTh_nrVerH%?!+Lw_;$bHeMM-W>0VD&#VOmjJ!WQ z);rD&cv2JT;eXro+p}qjyWeqOgSH50ZF4LB(uX3INBQ-B8JjO%YFtOc(&kR}38J#r zyl%bV?%d4st-wz6;?z~S50;O5dy<-eNs3%?6nk%fQ5XA@3TIOqJ_#$b$n6fp)-N`% zKFHkNe2gR+8h>RIaG#R6mh#sUD0C5oF{-n}is!buklswM@}=*=Sey6UERMAc*S(B@ zQ;Sq*Ew(Aii7P?=tdfyuz6;DZ&O$G%{?3`vVp1zIc=DO-D9nQuS3D+dfbsT$USFBF z{pbvznRx0yQ)I!)X@dak)hi$@H9E^;e@aCz+^m^K^^};{86`EA(}Kyjd(;o*49xwb z!A!^1M>IiXsccS?I-`vx35>e9v&-G;E)?Za-?bHxX2-g!oLboWF7nO(%AFh^nl+|G zxUkJm%&uK4z^Igu*=M13XRz?ASh%GIAtiKzw9A*q=PzLG%*oXw`q^BJV^cD&J<`u? z_lEieFHWVd1Uo7JbrjzvI;u7bv5{06HCd^rYaWi1LU+#`q-@Qj==ygVdp-KMCbpb} z4`G>z$iaZ&HA26&Dq{xi-`zp)a8DT<3rRbM)7(INclZZ9XWZqe-NIDwu_|$aCN=mI zxhY5ta*tx9Ra}(SHZ`DoY=WUMku{_~*Uk_iT)4hJ$5?9?997$a@;cXI*_R6mC{HK) z^6jnP$$cnGI%S1{wIb7K(LTP_FPRb&WyhoA&Q;8{BJx4vPUX#viYz5fRdJk29Xg!# z>E}*TZ^`9qf(ImRbBDJmlhArKRTx-U2(}-f%)WVL-Fuo#xdc`biJI+n%jZ_4wxvfD z#IJzEW2RD8r5H_|0%5Pnvl#Kg-_|_O^}Zfna&==_#|mW*DYR_ly=YbO%20}J_3cyW znxyHA*p+!Nr)f(URcn$xp+YtZ3An3)(V&8PS+fhY)Ds5%X^3ZwC!J$SO;7m_jHEBC zs-X-gUM@NG6tZ8aM9&A>ZA4V0&(DyYh!(b4g`QB96(m4>A5-}d33sPRQLutlJbKeMV9TVD*#&d!=K0Ne$zTnfMv$1eX1`FI-iI#DYkCNcvZ;V%` z%F%2B@^A=_+C4KMX};@L0@890&j0Ir%ImE6(zU%f0J~sTDN3TVcTp?!I#Bs+bYkBw z=Gt=)Fr*$`!1XR+HiK0mdfo=l(P_g+z_ zV@!PecokwVUdbdmB(b@9GLNE)YDqthQ#PS92^Md8%9)#ZwIkrnGB|}EqWyJ71rlgr zd>>YEQ|+GatFEz7=&%DJX)fVivX9r1jwtjQ?P(a(YHRaKBkPx#1ZP1_K>z@;j|==^1poj5Gf+%aMF0Q*000000000000000000000000000000 z000000000000000000000000000000qXkCE0000HbW%=J0L!@Ub|P4oFG}1z|Al9$ z6UX0Sr2qf`32;bRa{vGr5dZ)e5dq33^FIIp2GU7HK~!i%?VE|3vnmXRSzIA1_kXqL zOF-*rwbiy$|2e$ROreShNr13aeGlKJ$Z5!>2%=OsF{BbfUk$`x0f&atw_ZL;OHfcf zCQ_=F&*`Z>!ZQUQQg|Iq3it3Sd<;EE3Q1N}#WC)Q_^FO7|3N8jPEPqM;3K`kk7OD zWpTetNm~O{Ci88cny364()nD4243j3x{`|(r{5?}X_k)fw|OkLdOr9NJ)D?Jj%U^3 zHpas((tCCsKHMqiqp4wMfl1#d9B0Dw6C0@JctOkY2^Ymh#eanv3W7MIT1=NXqN^d^ z-xd<7erM22MUzb7#Kda#ODELC`-)0LOMY`rWMs<5D+c^_c?LJe;%-4YAWhSftnBQ; zlQ&t9*v!p!g&GQ$wS|$x=R#8MC>j(GUuur~1}>q@F1HzyH1ANkSKANIpN~y5sI0a< zPZ&vxoCCrcc5+|S?bYAgi#b&7Ytnfqgc-o_rs-TJlgVWM9z0H52YadF^!=rZV=tvF zxH=aqPB(Z@^&!{t(_E`%U9LDATV5U6JB#jL146rwPtR13Cxkx=yq@?k2(szpcqQOC z>Da>_KWCpZUk2ON5Ic3xpBnhdL1nqLauUT*0VW51h?%hd5z)RUyR%c_-kB{+uubBeHLL>rIKcSt?8N;(rJM z2@TsizqUb9BrWNZtG^@E&M5!n3Pn+=T9U3d8Pg&dFD+>NLvXtWoj5;T186YiA9jO(%#Z()(A+k{n|} z;^!K?e8A0UB?4q0g(iz~Tp0=Z)_r7qx3)?@++LHVL5{38{3bZBSPxDDYoRIxDNe^1 z@Q)JL9Ih)Hi5&wpE8!lzzD8U@*QF~f#D+x^7V<^lq6L%g@qVrknkRpM@G?%5q6@-6 zr+cMH$#x!%?n!nw9&kCF>$68B@vTlQDj@=m`1o7HO#{jKzl%ymkW$z(E_OeT}LEa6&vJC~v*%EH;e*fI;} zKBeXLp)8zx6tWtag>w%bcNbVV7do+a4S0EeVr%L8|4wdJJ7nQ(QgfqvxcfANkEDnC z&?EmhWMk~YbkdW?qx8&T6Z(Vr+Pdno7FwmAr=AH_vat?s;ZPo()}&MR>Aag?`U5nk-t! z5LJ!sr?SQ)gZTUFK7~wUZ58-Rl0ipZfIw)}xjwNun1nNI?nzI-iyGOflF*)R1s^RVH-VmNq*lo1VkEcG zc+o6R9M|AN?6$_k%?fDw{ot&_g+h;(ji{~I{tWfyIO;DP@!Jbn#UQxx=bh2;ZiZ8w z@;D(Jc|JJbQc?A|a3-4{q>>=x&_%h|$A@i0R{9sT`^j0X$)fH}otg+_RceW%7PTq| zEvLm87v+ZQZNrD%q?PrDU8|d)+DFr$Z1pS{ltf*Xr9l=nxObd4nWrdHLu}`!3g?Jy zJsM$asZM=m30Zbw@@3HA;+rAbc%`tUqfwp1id5P|N`^*E>nn_zhim^f(3q0As@JW4 zgGyA?9BTk^SP~(f=}*! fOeT~0`zZAv@-;Y>!#kM~00000NkvXXu0mjfm6(xc literal 0 HcmV?d00001 diff --git a/doc/dataExchangeModel_twoWayFiles.txt b/doc/dataExchangeModel_twoWayFiles.txt index 9561b12d..d2cb2a38 100644 --- a/doc/dataExchangeModel_twoWayFiles.txt +++ b/doc/dataExchangeModel_twoWayFiles.txt @@ -16,11 +16,13 @@ dataExchangeModel twoWayFiles; twoWayFilesProps \{ couplingFilename "filename"; - maxNumberOfParticles number; + maxNumberOfParticles scalar1; + DEMts scalar2; \} :pre {filename} = filename of the VTK file series :ulb,l -{number} = maximum number of particles in DEM simulation :l +{scalar1} = maximum number of particles in DEM simulation :l +{scalar2} = DEM time step width :l :ule [Examples:] diff --git a/doc/liggghtsCommandModel_execute.txt b/doc/liggghtsCommandModel_execute.txt index e008f7a9..4d404c6d 100644 --- a/doc/liggghtsCommandModel_execute.txt +++ b/doc/liggghtsCommandModel_execute.txt @@ -23,18 +23,18 @@ executeProps0 run $couplingInterval ); - runFirst switch1; + runFirst switch1; runLast switch2; runEveryCouplingStep switch3; runEveryWriteStep switch4; verbose; \} :pre -{command} = LIGGGHTS command to be executed. Each word in a new line, numbers and symbols need special treatment (e.g. $couplingInterval will be replaced by correct coupling interval in the simulation) :ulb,l -{switch1} = switch (choose on/off) if the command is executed only at first time step :l -{switch2} = switch (choose on/off) if the command is executed only at last time step :l -{switch3} = switch (choose on/off) if the command is executed at every coupling step :l -{switch4} = switch (choose on/off) if the command is executed at every writing step :l +{command} = LIGGGHTS command to be executed. Each word in a new line, numbers and symbols need special treatment (e.g. $couplingInterval will be replaced by correct coupling interval in the simulation) :ulb,l +{switch1} = (optional, default off) if the command is executed only at first time step :l +{switch2} = (optional, default off) if the command is executed only at last time step (requires {switch1} to be off) :l +{switch3} = (optional, default off) if the command is executed at every coupling step (requires {switch1} and {switch2} to be off) :l +{switch4} = (optional, default off) if the command is executed at every writing step (requires {switch1} to {switch3} to be off) :l {verbose} = (normally off) for verbose run :l :ule diff --git a/doc/liggghtsCommandModel_readLiggghtsData.txt b/doc/liggghtsCommandModel_readLiggghtsData.txt index e2fc588d..29e9d71c 100644 --- a/doc/liggghtsCommandModel_readLiggghtsData.txt +++ b/doc/liggghtsCommandModel_readLiggghtsData.txt @@ -18,9 +18,26 @@ liggghtsCommandModels ); readLiggghtsDataProps0 \{ - ??? + startIndex scalar1; + verbose; + exactTiming; + filePath + ( + "word" + ); + startTime scalar2; + endTime scalar3; + timeInterval scalar4; \} :pre +{scalar1} = start index of data file to be read; the index is appended to {filePath} :ulb,l +{verbose} = (default off) flag for verbose run :l +{exactTiming} = flag indicating that start time should be kept even during a coupling interval :l +{filePath} = path to LIGGGHTS data file. Each word starts in a new line; special characters, i.e. dots and slashes, need special treatment (e.g. dotdot will be replaced by "..") :l +{scalar2} = start reading at this time :l +{scalar3} = end reading at this time :l +{scalar4} = repeat reading at this time interval while increasing the data file index :l +:ule [Examples:] @@ -30,13 +47,27 @@ liggghtsCommandModels ); readLiggghtsDataProps0 \{ - ??? + startIndex 0; + exactTiming; + filePath + ( + dotdot + slash + DEM + slash + packing.data + ); + startTime 0.002; + endTime 0.012; + timeInterval 0.001; \} :pre [Description:] The {readLiggghtsData} liggghtsCommandModel can be used to read LIGGGHTS data -files into LIGGGHTS during runtime of a coupled simulation. +files into LIGGGHTS during runtime of a coupled simulation. This corresponds to +the {read_data} command in LIGGGHTS with the {add} option, i.e. read in +particles are added to existing particles. IMPORTANT NOTE: Model is outdated. diff --git a/doc/liggghtsCommandModel_writeLiggghts.txt b/doc/liggghtsCommandModel_writeLiggghts.txt index f3f2f05b..53bad2c9 100644 --- a/doc/liggghtsCommandModel_writeLiggghts.txt +++ b/doc/liggghtsCommandModel_writeLiggghts.txt @@ -20,15 +20,17 @@ liggghtsCommandModels writeLiggghtsProps \{ writeLast switch1; + path "path"; writeName "name"; overwrite switch2; verbose; \} :pre -{switch1} = switch (choose on/off) to select if only last step is stored or every write step (default on). :ulb,l -{name} = name of the restart file to be written in /$caseDir/DEM/ default (default "liggghts.restartCFDEM") :l -{switch2} = switch (choose on/off) to select if only one restart file $name or many files $name_$timeStamp are written (default off):l -{verbose} = (default off) for verbose run :l +{switch1} = (optional, default on) select if only last step is stored or every write step :ulb,l +{path} = (optional, default "../DEM") alternative path (relative to execution directory) for saving the restart file :l +{name} = (optional, default "liggghts.restartCFDEM") name of the restart file to be written in ../DEM/ :l +{switch2} = (optional, default off) select if only one restart file $name or multiple files $name_$timeStamp are written :l +{verbose} = (optional, default off) for verbose run :l :ule [Examples:] diff --git a/doc/momCoupleModel_implicitCouple.txt b/doc/momCoupleModel_implicitCouple.txt index aa5da194..f746477b 100644 --- a/doc/momCoupleModel_implicitCouple.txt +++ b/doc/momCoupleModel_implicitCouple.txt @@ -21,13 +21,15 @@ implicitCoupleProps velFieldName "U"; granVelFieldName "Us"; voidfractionFieldName "voidfraction"; - minAlphaP number; + KslLimit scalar1; + minAlphaP scalar2; \} :pre {U} = name of the finite volume fluid velocity field :ulb,l {Us} = name of the finite volume granular velocity field :l {voidfraction} = name of the finite volume voidfraction field :l -{number} = minimum value for local particle volume fraction to calculate the exchange filed (default SMALL) :l +{scalar1} = (optional, default 1e10) limit implicit momentum exchange field :l +{scalar2} = (optional, default SMALL) minimum value for local particle volume fraction to calculate the exchange field :l :ule [Examples:] diff --git a/doc/smoothingModel_constDiffSmoothing.txt b/doc/smoothingModel_constDiffSmoothing.txt index c6d83f2e..1bc69bde 100644 --- a/doc/smoothingModel_constDiffSmoothing.txt +++ b/doc/smoothingModel_constDiffSmoothing.txt @@ -19,12 +19,14 @@ constDiffSmoothingProps upperLimit number2; smoothingLength lengthScale; smoothingLengthReferenceField lengthScaleRefField; + verbose; \} :pre {number1} = scalar fields will be bound to this lower value :ulb,l {number2} = scalar fields will be bound to this upper value :l {lengthScale} = length scale over which the exchange fields will be smoothed out :l {lengthScaleRefField} = length scale over which reference fields (e.g., the average particle velocity) will be smoothed out. Should be always larger than lengthScale. If not specified, will be equal to lengthScale. :l +{verbose} = (optional, default false) flag for debugging output :l :ule [Examples:] diff --git a/doc/voidFractionModel_GaussVoidFraction.txt b/doc/voidFractionModel_GaussVoidFraction.txt index f382bd06..0c792001 100644 --- a/doc/voidFractionModel_GaussVoidFraction.txt +++ b/doc/voidFractionModel_GaussVoidFraction.txt @@ -23,8 +23,8 @@ GaussProps {number1} = maximum number of cells covered by a particle (search will fail when more than {number1} cells are covered by the particle) :ulb,l {number2} = minimum limit for voidfraction :l -{number3} = (optional) scaling of the particle volume to account for porosity or agglomerations. :l -{number4} = (optional) diameter of the particle's representation is artificially increased according to {number2} * Vparticle, volume remains unaltered! :l +{number3} = (optional, default 1.0) scaling of the particle volume to account for porosity or agglomerations. :l +{number4} = (optional, default 1.0) diameter of the particle's representation is artificially increased according to {number2} * Vparticle, volume remains unaltered! :l :ule [Examples:] diff --git a/doc/voidFractionModel_bigParticleVoidFraction.txt b/doc/voidFractionModel_bigParticleVoidFraction.txt index 495d1714..724b0489 100644 --- a/doc/voidFractionModel_bigParticleVoidFraction.txt +++ b/doc/voidFractionModel_bigParticleVoidFraction.txt @@ -23,8 +23,8 @@ bigParticleProps {number1} = maximum number of cells covered by a particle (search will fail when more than {number1} cells are covered by the particle) :ulb,l {number2} = minimum limit for voidfraction :l -{number3} = (optional) scaling of the particle volume to account for porosity or agglomerations. :l -{number4} = (optional) diameter of the particle's representation is artificially increased according to {number2} * Vparticle, volume remains unaltered! :l +{number3} = (optional, default 1.0) scaling of the particle volume to account for porosity or agglomerations. :l +{number4} = (optional, default 1.0) diameter of the particle's representation is artificially increased according to {number2} * Vparticle, volume remains unaltered! :l :ule [Examples:] @@ -34,7 +34,7 @@ bigParticleProps \{ maxCellsPerParticle 1000; alphaMin 0.10; - weight 1.; + weight 1.0; porosity 5.0; \} :pre diff --git a/doc/voidFractionModel_dividedVoidFraction.txt b/doc/voidFractionModel_dividedVoidFraction.txt index d87cda43..b4e5906a 100644 --- a/doc/voidFractionModel_dividedVoidFraction.txt +++ b/doc/voidFractionModel_dividedVoidFraction.txt @@ -19,12 +19,16 @@ dividedProps interpolation; weight number2; porosity number3; + procBoundaryCorrection switch1; + verbose; \} :pre {number1} = minimum limit for voidfraction :ulb,l {interpolation} = flag to interpolate voidfraction to particle positions (normally off) :l {number2} = (optional) scaling of the particle volume to account for porosity or agglomerations. :l {number3} = (optional) diameter of the particle's representation is artificially increased according to {number2} * Vparticle, volume remains unaltered! :l +{switch1} = (optional, default false) allow for correction at processor boundaries. This requires the use of engineIB and vice versa. :l +{verbose} = (optional, default false) flag for debugging output :l :ule [Examples:] @@ -41,12 +45,38 @@ The {divided} voidFraction model is supposed to be used when a particle (or its representation) is in the size range of a CFD cell. Satellite points are used to divide the particle's volume to the touched cells. +The particle has radius R and its volume is divided into 29 non-overlapping +regions of equal volume. The centroids of these volumes are then used to +reproduce each volume. The first volume is a sphere with the center coinciding +with the particle center. The radius of this subsphere can be found as follows: + +:c,image(Eqs/voidfractionModel_divided_pic2.png) + +The rest of the volume is a spherical layer that is divided into 2 layers of +equal volume. Position of the border between these two spherical layers in +radial direction can be easily obtained: + +:c,image(Eqs/voidfractionModel_divided_pic3.png) + +Each of these spherical layers is later divided into 14 elements of equal volume. +Position of the centroid point in radial direction of each volume in the first +spherical layer is as follows + +:c,image(Eqs/voidfractionModel_divided_pic4.png) + +Similarly, for the second spherical layer remembering that the external radius +is the particle radius: + +:c,image(Eqs/voidfractionModel_divided_pic5.png) + The region of influence of a particle can be increased artificially by -"porosity", which blows up the particles, but keeps their volume (for -voidfraction calculation) constant. +{porosity}, which blows up the particles, but keeps their volume (for +void fraction calculation) constant. The particle volume occupied in the CFD domain can be adjusted by the parameter -"weight", using Vparticle=dsphere^3*pi/6*weight. +{weight}, using + +:c,image(Eqs/voidfractionModel_divided_pic6.png). In the basic implementation of solvers, the void fraction is calculated based on all particles. Depending on the solver used, the void fraction calculation is