The Numerical Solution of Moving Boundary ...

European Symposium Symposium on on Computer Computer Arded Aided Process European Process Engineering Engineering –– 15 15 L. Puigjaner Puigjaner and and A. A. Espuña Espuña (Editors) (Editors) L. © 2005 2005 Elsevier Elsevier Science Science B.V. B.V. All All rights rights reserved. reserved. ©

7KH1XPHULFDO6ROXWLRQRI0RYLQJ%RXQGDU\3UREOHPV 8VLQJWKH0RYLQJ)LQLWH(OHPHQW0HWKRG 5XL5REDORD&DUORV6HUHQR0DULDGR&DUPR&RLPEUDEF DQG$OtULR5RGULJXHVEG D

'HSDUWDPHQWRGH0DWHPiWLFD8QLYHUVLGDGHGD%HLUD,QWHULRU 5XD0DUTXrV'¶ÈYLODH%RODPD&RYLOKm3RUWXJDO E

/DERUDWRU\RI6HSDUDWLRQDQG5HDFWLRQ(QJLQHHULQJ '(&)DFXOGDGHGH(QJHQKDULDGD8QLYHUVLGDGHGR3RUWR G '(4)DFXOGDGHGH(QJHQKDULDGD8QLYHUVLGDGHGR3RUWR 5XD'U5REHUWR)ULDVVQ3RUWR3RUWXJDO F

$EVWUDFW

$PRYLQJILQLWHHOHPHQW PHWKRG0)(0 LVGHYHORSHGIRUWKHQXPHULFDOVROXWLRQVRI WLPHGHSHQGHQWSDUWLDOGLIIHUHQWLDOHTXDWLRQVLQYROYLQJVROLGIOXLGUHDFWLRQV2XU0)(0 JHQHUDWHV DQ DGDSWLYH PHVK IRU HDFK GHSHQGHQW YDULDEOH VR WKURXJK VSDWLDO GRPDLQ GHFRPSRVLWLRQ LW FDQ EH PRGLILHG LQ RUGHU WR EH DQ HIILFLHQW VROYHU IRU D FODVV RI SUREOHPV VKRZLQJ D PRYLQJ LQWHUQDO ERXQGDU\ 7KH DOJRULWKP ZDV WHVWHG LQ WKH QXPHULFDO VLPXODWLRQ RI D QRQ FDWDO\WLF VROLGIOXLG UHDFWLRQ PRGHOOHG E\ WKH D  LVRWKHUPDOE QRQLVRWKHUPDOVKULQNLQJFRUHPRGHOZLWKQRQOLQHDUNLQHWLFV7KLVZRUN HVWDEOLVKHV WKH HIIHFWLYHQHVV DQG DSSOLFDELOLW\ RI WKH 0)(0 LQ VROYLQJ PRYLQJ ERXQGDU\SUREOHPV .H\ZRUGV 0RYLQJ )LQLWH (OHPHQWV 0RYLQJ %RXQGDU\ 3UREOHPV $GDSWLYH *ULGV 7LPH'HSHQGHQW3DUWLDO'LIIHUHQWLDO(TXDWLRQV6KULQNLQJ&RUH0RGHO

,QWURGXFWLRQ $ODUJHQXPEHURISUREOHPVLQVFLHQFHDQGHQJLQHHULQJDUHIRUPXODWHGLQWHUPVRIWLPH GHSHQGHQW SDUWLDO GLIIHUHQWLDO HTXDWLRQV LQYROYLQJ D PRYLQJ LQWHUIDFH 7KH DQDO\WLF VROXWLRQV RI PDWKHPDWLFDO PRGHOV RI G\QDPLF V\VWHPV DUH UDUHO\ NQRZQ 7KLV KDV PRWLYDWHG PDQ\ DXWKRUV DQG GLIIHUHQW WHFKQLTXHV KDYH EHHQ SURSRVHG LQ UHFHQW \HDUV &UDQN   SUHVHQWV D ZLGHUDQJLQJ GHVFULSWLRQ RQ PDWKHPDWLFDO DQG QXPHULFDO PHWKRGV IRU IUHH DQG PRYLQJ ERXQGDU\ SUREOHPV ,Q WKLV SDSHU ZH FRQVLGHU WKH GHYHORSPHQWRID0)(0WRFRPSXWHQXPHULFDOVROXWLRQRIRQHGLPHQVLRQDOSUREOHPV LQYROYLQJ WZR SKDVH ZLWK D PRYLQJ LQWHUIDFH 7KH 0)(0 6HUHQR HW DO  &RLPEUDHWDO JHQHUDWHVDQDGDSWLYHPHVKVRLWFDQEHPRGLILHGLQRUGHU WREHDQHIILFLHQWVROYHUIRUDFODVVRISUREOHPVVKRZLQJDPRYLQJLQWHUQDOERXQGDU\7R PRGLI\ WKH PRYLQJ ILQLWH HOHPHQW HTXDWLRQV ZH SHUIRUP WKH GHFRPSRVLWLRQ RI WKH VSDWLDO GRPDLQ E\ WKH LQWURGXFWLRQ RI D PRYLQJ QRGH GHVFULELQJ WKH SRVLWLRQ RI WKH LQWHUQDOLQWHUIDFH7KLVIRUPXODWLRQ5REDORHWDO LVH[SDQGHGLQWKLVZRUNWRD 

$XWKRUWRZKRPFRUUHVSRQGHQFHVKRXOGEHDGGUHVVHGPFRLPEUD#IHXSSW

ODUJHUFODVVRISUREOHPVE\DVVXPLQJWKHH[LVWHQFHRIQRQOLQHDUERXQGDU\FRQGLWLRQVDW WKH PRYLQJ ERXQGDU\ 7KH SDSHU LV RUJDQL]HG DV IROORZV ,Q VHFWLRQ  ZH SUHVHQW WKH SUREOHPWREHVROYHG7KHGHWDLOHGGHVFULSWLRQRIWKHPHWKRGLVSUHVHQWHGLQVHFWLRQ )LQDOO\ LQ VHFWLRQ  ZH DSSO\ WKH FRPSXWHU FRGH UHVXOWLQJ IURP WKH QXPHULFDO DOJRULWKP LPSOHPHQWDWLRQ WR WKH VLPXODWLRQ RI D QRQ FDWDO\WLF VROLGIOXLG UHDFWLRQ &DUDELQHWDO PRGHOOHGE\WKHD LVRWKHUPDOE QRQLVRWKHUPDOVKULQNLQJFRUH PRGHOZLWKQRQOLQHDUNLQHWLFV

7KH3UREOHP /HWXVGHILQHWKHPRYLQJERXQGDU\SUREOHPWKDWZHFRQVLGHULQWKLVZRUN7KHJHQHUDO PDWKHPDWLFDOPRGHORIWZRSKDVHV\VWHPLQRQHGLPHQVLRQDOIL[HGVSDFHGRPDLQZLWK DQ LQWHUQDO PRYLQJ LQWHUIDFH LV GHILQHG E\ D V\VWHP RI Q SDUDEROLF SDUWLDO GLIIHUHQWLDO HTXDWLRQV3'( ZKRVHPWKHTXDWLRQILOOWKHIRUP   ∂ \P  ∂ \ ⎞ ∂  \P  ∂\ ⎞ ⎛ ⎛    = I P ⎜ [ W  \  + J P ⎜ [ W  \  ∂W ∂[ ⎟⎠ ∂[  ∂[ ⎟⎠ ⎝ ⎝ RQ D ≤ [ ≤ E  W >   :H FRQVLGHU WKDW WKH VSDFH YDULDEOH VDWLVI\ D ≤ [ ≤ ; V ( W )  LI P ≤ Q, DQG ; V ( W ) ≤ [ ≤ E LI Q, < P ≤ Q  Q,  Q,, DUHWKHQXPEHURIGHSHQGHQWYDULDEOHV

LQ HDFK SKDVH ; V ( W )  LV WKH SRVLWLRQ RI WKH PRYLQJ ERXQGDU\ DW LQVWDQW W DQG 7 7      7 7 \ = ⎡⎣ \ ,  \ ,, ⎤⎦ = [ \  \Q ] ZKHUH \ , 7 = ⎡⎣ \,  \Q,, ⎤⎦ DQG \ ,, 7 = ⎡⎣ \,,  \Q,,,, ⎤⎦ /HW XVDVVXPHVRPHLQLWLDOFRQGLWLRQVDQGWKDWHTXDWLRQ LVVXEMHFWWR'LULFKOHWRU5RELQ ERXQGDU\FRQGLWLRQVDWWKHHQGSRLQWVRIVSDFHGRPDLQ$WWKHLQWHUIDFHZHVXSSRVHWKDW WKHERXQGDU\FRQGLWLRQVFDQEHH[SUHVVHGE\JHQHUDOLPSOLFLWFRQGLWLRQV

 ⎛  ∂\ KP ⎜ ; V  W  \ ;  V ⎜ ∂[ ⎝

 ;V

 ∂\, ∂[

 ; V−

 ∂ \ ,, ∂[

; V+

⎞ ⎟ =   ⎟ ⎠

 

7RFRPSOHWHWKHPRGHOOHWXVDVVXPHWKDWWKHLQWHUIDFHPRYHPHQWLVJRYHUQHGE\  ⎛ G; V  ∂\ = Z ⎜ ; V  W \ ;  V ⎜ GW ∂[ ⎝

 ;V

 ∂\, ∂[

 ; V−

 ∂ \ ,, ∂[

; V+

⎞ ⎟ ⎟ ⎠

 

7KLV FODVV RI SUREOHPV LV DQ H[WHQVLRQ RI WKRVH FRQVLGHUHG LQ 6HUHQR HW DO  &RLPEUD HW DO   ZKHUH DOO GLIIHUHQWLDO HTXDWLRQV KDYH WKH VDPH IL[HG RQH GLPHQVLRQDOVSDWLDOGRPDLQEXWHDFKGHSHQGHQWYDULDEOHKDVLWVRZQPRYLQJPHVK7KH PDLQ LGHD WR VROYH WKLV SUREOHP LV WR FRPELQH WKH GRPDLQ GHFRPSRVLWLRQ DQG WKHVH GLIIHUHQWPRYLQJILQLWHHOHPHQWPHVKHV

0HWKRG'HYHORSPHQW 7KH 0)(0 LV D GLVFUHWL]DWLRQ SURFHVV LQ WZR VWDJHV ILUVW VSDWLDO GLVFUHWL]DWLRQ XVLQJ ILQLWHHOHPHQWVDQGDOORZLQJWKHPRYHPHQWRIWKHVSDFHQRGHVLQZKLFKZHIRFXVRXU DWWHQWLRQ DQG VHFRQGO\ WKH WLPH LQWHJUDWLRQ RI WKH UHVXOWLQJ RUGLQDU\ GLIIHUHQWLDO

HTXDWLRQV2'( V\VWHP7RVROYHWKLVLPSOLFLWWLPHGHSHQGHQW2'(V\VWHPZHXVHWKH SDFNDJH/62',+LQGPDUVK  6SDWLDO'LVFUHWL]DWLRQ 7KHJULGFRQQHFWHGWRWKH PWKGHSHQGHQWYDULDEOH \P LVREWDLQHGE\SDUWLWLRQLQJWKH VSDWLDOGRPDLQLQ 1 P ILQLWHHOHPHQWVE\ 1 P −  LQWHULRUVHSDUDWLRQQRGHV 3P  ; P  ( W ) < ; P  ( W ) <  < ; P  1 P ( W ) < ; P  1P + ( W )  W ≥  

 

2QH RI WKH QRGHV ; P   RU ; P  1 P +  LV IL[HG DQG WKH RWKHU LV GHILQHG E\ WKH LQWHUIDFH SRVLWLRQDORQJWLPH,QHDFKRIWKHILQLWHHOHPHQWRIJULG 3P ZHGHILQHWKHDSSUR[LPDWH VROXWLRQE\DSRO\QRPLDOREWDLQHGE\/DJUDQJHLQWHUSRODWLRQ7KHUHODWLYHSRVLWLRQVRI WKHLQWHUSRODWLRQSRLQWVDUHRSWLPL]HGDVLQWKHRUWKRJRQDOFROORFDWLRQPHWKRG6HUHQR   7KH GHJUHH RI WKH SRO\QRPLDO DSSUR[LPDWLRQ PD\ EH FKRVHQ E\ WKH XVHU 7KLV ORFDODSSUR[LPDWLRQJXLGHVWRDJOREDODSSUR[LPDWLRQ

  =  + ∂W [ ∂[ ∂[

 

ZLWKLQLWLDOFRQGLWLRQVν  = ν  =  DQGERXQGDU\FRQGLWLRQV ∂ν  =  ∂[

∂ν  = − %LK &  [ =   ∂[

ν  = ν   −κ

 

∂ν  ∂ν  + − α & Q H − β +ν  =  [ = ; 6 W  ∂[ ∂[

 

%LK LVWKH%LRWQXPEHUIRUKHDWWUDQVIHU α LVWKHGLPHQVLRQOHVVKHDWRIUHDFWLRQDQG

κ =  )LJXUHVKRZVWKHWHPSHUDWXUHD DQGFRQFHQWUDWLRQE SURILOHVIRU Q =  

 )LJXUH7HPSHUDWXUHD DQGFRQFHQWUDWLRQE SURILOHVDWGLIIHUHQWWLPHV

7KLV UXQ WHVW ZDV SHUIRUPHG FRQVLGHULQJ β =   γ = × −  %LP = %LK =   α = − ×  DQG 'DW =  =  :HXVHILQLWHHOHPHQWVLQWKHJULGDVVRFLDWHGWR

νDQGLQWKHWZRRWKHUJULGVDVVRFLDWHGWR νDQG &8VLQJFXELFDSSUR[LPDWLRQVWKH &38 WLPH IRU WKLV UXQ LV  VHFRQGV )LJXUH D  DQG E  SUHVHQW WKH WUDMHFWRULHV RI QRGHV VHSDUDWLRQ DVVRFLDWHG WR WHPSHUDWXUH DQG FRQFHQWUDWLRQ UHVSHFWLYHO\ ,QLWLDOO\ QRGHVDUHFRQFHQWUDWHGQHDU[ 

 )LJXUH7HPSHUDWXUHD DQGFRQFHQWUDWLRQE DWQRGHVSRVLWLRQDQGLQWHUIDFH

&RQFOXVLRQV ,Q WKLV SDSHU ZH KDYH SUHVHQWHG D 0)(0 IRU PRYLQJ ERXQGDU\ SUREOHPV 7KH PDLQ DGYDQWDJHV RI WKH PHWKRG DUH LWV VLPSOLFLW\ DQG SHUIRUPDQFH ZKLFK DFFRXQW IRU LWV ZLGHVSUHDGXVH7KHFRPSOH[LW\RIWKHQXPHULFDOH[DPSOHVFRQVLGHUHGUHTXLUHVDUREXVW DQG HIILFLHQW DOJRULWKP ,Q IDFW WKH QXPHULFDO UHVXOWV REWDLQHG E\ WKLV IRUPXODWLRQ RI 0)(0H[KLELWYHU\JRRGDJUHHPHQWZLWKWKRVHREWDLQHGE\SUHYLRXVPHWKRGVVRWKLV ZRUN HVWDEOLVKHV WKH HIIHFWLYHQHVV DQG DSSOLFDELOLW\ RI WKH 0)(0 LQ VROYLQJ PRYLQJ ERXQGDU\SUREOHPV 5HIHUHQFHV

&DUDELQ 3 DQG ' %HUN  $QDO\VLV DQG PRGHOOLQJ RI WKH LVRWKHUPDO DQG QRQLVRWKHUPDO VKULQNLQJFRUHPRGHOZLWKQRQOLQHDUNLQHWLFV&KHP(QJQJ6FL &RLPEUD0&&$6HUHQRDQG$(5RGULJXHV$SSOLFDWLRQVRIDPRYLQJILQLWHHOHPHQW PHWKRG&KHPLFDO(QJLQHHULQJ-RXUQDO &RLPEUD0&&$6HUHQRDQG$(5RGULJXHV$(0RYLQJILQLWHHOHPHQWPHWKRG DSSOLFDWLRQWRVFLHQFHDQGHQJLQHHULQJSUREOHPV&RPS&KHP(QJQJ9ROSS &UDQN-)UHHDQG0RYLQJ%RXQGDU\3UREOHPV&ODUHQGRQ3UHVV2[IRUG +LQGPDUVK $ &  /62'( DQG /62', 7ZR 1HZ ,QLWLDO 9DOXH 2UGLQDU\ 'LIIHUHQWLDO (TXDWLRQ6ROYHUV$&06,*1801HZVOHWWHU9RO1RSS 0LOOHU.DQG510LOOHU0RYLQJ)LQLWH(OHPHQWV,6,$0-1XPHU$QDO9RO1R SS 5REDOR 5 0& &RLPEUD 6HUHQR & $ DQG 5RGULJXHV $ (  $SOLFDomR GR 0()0 D SUREOHPDV FRP IURQWHLUDV PyYHLV &RQJUHVVR GH 0pWRGRV &RPSXWDFLRQDLV HP (QJHQKDULD HGLWHGE\&0RWD6RDUHV$%DWLVWD*%XJHGD0&DVWDOHLUR-*RLFROHD-0DUWLQV&3LQD DQG+5RGULJXHV(GLWLRQV/1(&,6%1SS 6HUHQR&$$(5RGULJXHVDQG-9LOODGVHQ7KH0RYLQJ)LQLWH(OHPHQW0HWKRGZLWK 3RO\QRPLDO$SSUR[LPDWLRQRI$Q\'HJUHH&RPS&KHP(QJQJ9RO1RSJ 6HUHQR&$$(5RGULJXHVDQG-9LOODGVHQ6ROXWLRQRI3DUWLDO'LIIHUHQWLDO(TXDWLRQV 6\VWHPVE\0)(0&RPS&KHP(QJQJ9RO