Xây dựng mô hình mô phỏng dòng chảy và truyền chất hai chiều, ứng dụng cho vùng vịnh hạ long

DAI HOC QUOC GIA HA NQI Trudng Dai hoc Khoa hoc Tir nhien XkY DUNG MO HINH MO PHONG DONG CHAY VA TRUYEN CHAT HAI CHKU, VTNG DUNG CHO YUNG VINH HA LONG MQ SO: Q T Chu Iri 3S kcii: TS Phung Dang Hie'u CQC thcinh vien thorn gio: ThS Phqm Hodng Lam CN Trdn Dvcc Trvt OAI H O C Q U O C GIA HA NO: TRUNG T.Ar/ THONG TIN THU ' " ^ ^ DT / 64^ HA NOI - 2005 BAO CAO TOM TAT D £ TAI Ten de tai: Xay dung m6 hinh mo phdng ddng chay va truydn chat hai chieu, Ihig dung cho viing vinh Ha Long Ma sd: QT.06.34 Chu tri de tai: TS Phiing Dang Hieu Cac can bo tham gia: 1) ThS Pham Hoang Lam 2) CN Trdn Dixc Trur Muc tieu va not dung nghien curu: 1) Muc tieu: - Xay dung dugc mo hinh so m6 phdng ddng chay va truyen chat dang hai chidu sir dung h6 phuong trinh nude nong phi tuyen va phuong trinh truydn tai khuech tan - LSig dung tinh loan thir nghiem cho viing vinh Ha Long 2) Npi dung: a) xay dung mo hinh so mo phong ddng chay va truydn chat hai chieu - Nghien cihi thuat toan thich hop sai phan giai he phuong trinh nude nong hai chi6u Lap chuong trinh may tinh - Nghien cuu thuat toan thich hop sai phan giai he phuong trinh truyen chat Lap chuong trinh may tinh - Tinh toan va kiem nghiem mo hinh so theo cac di6u kien ly thuyet va thuc te b) Ung dung mo hinh tinh toan thir nghiem cho viing vinh Ha Long - Thu thap va so hoa trudng dp sau viing vinh Ha Long; thu thap cac sd lieu v6 thuy trieu va gid viing bien ven bd Quang Ninh - ung dung mo hinh tinh toan thir nghiem cho viing vinh Ha Long theo cac kich ban - Phan tich ket qua, danh gia mo hinh va viet bao cao long ket, nop san pham, nghiem thu Cac ket qua dat dugc: 1) Da utig dung phuong phap CIP (Constraint Interpolation Profile) hien dai cd dp chinh xac bac ba va phat trie'n mot phuong phap sai phan an vao viec sai phan hoa giai 11 h6 phuong trinh nude ndng phi tuy^n cho bai toan ddng chay Viet cong chuong trinh may tinh m6 phdng ddng chay sinh gid va thuy trilu Tinh toan thu nghidm va ki^m chung md hinh cho cac dieu kien ly thuyet cd nghiem giai tich va kilm chiing trudng hop dieu kien thuc te ciia bai toan truydn trieu Bien Dong K^t qua cho tha^y phuong phap CIP sir dung d day cho bai toan mo phdng ddng chay la r^t hop ly va cd d6 chinh xac kha cao Mo hinh chiing to kha nang ung dung tot cho cac viing bi^n nong ven bd 2) Da xay dung cong chuong trinh truyin tai va khuech tan vat cha't su dung phuong phap th^ tich huu han vdi so TVD (Total Variation Diminishing) cd dp chinh xac bac hai Viet cong chuong trinh may tinh va kiem chung mo hinh cho cac bai toan ly thuyet cd nghiem giai tich Ket qua cho thay phuong phap TVD siJ dung d day cd dp chinh xac cao hon nhieu so vdi mot sd phuong phap sai phan truyen thdng hay dugc sir dung nhu UPWIND bac mot hay Lax-Wendrof Mo hinh cho phep md phdng qua trinh truydn tai khuech tan vat cha't dudi tac dpng cua ddng chay trung binh 3) Da thu thap va sd hoa dia hinh viing vinh Ha Long tren ludi vuong gdc toa dp D6 Cac vdi budc ludi 250 m x 250 m Thu thap va phan tich nhung sd lieu v6 tdc dp gid va hudng gid viing ven bd Quang Ninh, sd lieu thiiy tridu LSig dung mo hinh m6 phdng he thdng ddng chay gay gid trung binh mija Dong va miia He va tac dpng ciaa thuy trieu cac chu ky trilu len va xudng Idn Mo phdng Ian truyen vat cha't dugc thuc hien tren n6n ddng chay long hop ca gid va thuy tri^u vdi gia thiet chat chi thi bin vung khong phan huy va cd sir cd mo trudng thai tai mot sd di^m gia dinh Cac ket qua cho cac bu'c tranh phan bd trudng ddng chay va phan bd nong dp cha't nhilm cd tac dpng ciia cac yeu td dpng lire la gid va thuy trieu 4) Ndi chung md hinh xay dung di tai da dap ihig dugc yeu cau v6 mo phdng ddng chay va truyen chat vdi dp chinh xac dang tin cay va cd kha nang urng dung de nghien cihi cho cac viing bien ven bd khac Ket qua ung dung tinh toan cho viing vinh Ha Long, Quang Ninh cho thay ro kha nang ung dung ciia md hinh da phat trien 5) Dang dugc mot bai bao tren tap chi khoa hpc Dai hoc Qudc Gia Tinh hinh kinh phi cua de tai Tdng kinh phi dugc cap: 20.000.000 ddng Da nhan: 20.000.000 ddng Da toan: 20.000.000 dong 111 XAC NHAN CHU TRI D £ CUA BAN CHU NHlfiM KHOA TAX (^>^>t^-^^ TS Phung Ddng Hieu PGS TS Pham Van Hudn XAC NHAN CUA TRU^ONG OHO HIEU TRUONU GS.TS.^^rfin/ ^ / ^ / ^ IV ABSTRACT Project title: Development of a numerical 2-D model based on the shallow water equation and convection-diffusion equation for the application to the Halong Bay Code number: QT.06.34 Project Leader: Dr Phung Dang Hieu 4- Members: 1) M Sc Pham Hoang Lam 2) B.Sc Tran Due Tru Aims and contents of project: 1) Aims: - Development of a numerical 2-D model based on the shallow water equation and convection-diffusion equation - Application of the model to the study of tide-induced currents and windinduced currents, the transportation of tracer in the Halong Bay 2) Contents: a) Development of a numerical 2-D model - Development and application of a finite difference scheme for solving the shallow water equation Development of a computer program - Development and application of a finite volume TVD scheme for solving the conveetion-diffution equation Development of a computer program - Verification of the numerical model with some analytical and practical conditions b) Application of the model to the study of tide-induced currents and windinduced currents, the transportation of tracer in the Halong Bay - Collection and digitization of the bathymetry of Halong Bay; Collection of wind and tidal data in the near-shore area of Quang Ninh - Application of the model to the simulation of tide-induced, wind-induced currents in the Halong Bay Simulation of the propagation and distribution of pollution (tracer) based on some scenarios - Analysing and assessing the simulated results in order to give conclusions and writing the final report Results: 1) The CIP (Constraint Interpolation Profile) method and a implicit finite difference scheme have been applied to the developement of a numerical model based on the shallow water equation A numerical computer program has been successfully developped for the simulation of wind and tide induced currents Computation has been carried out using the numerical model for the caseses of analiytical and practical problems Computed results show that the application of CIP method for the simulation of shallow water currents is fairly accurate The numerical model can be applied to the near-shore shalow seas 2) A numerical model based on the convection-diffusion equation using the TVD (Total Variation Diminishing) scheme has been successfully developed A computer program has developed Some numerical tests with analytical conditions show that the TVD method used here is more accurate than the conventional methods such as Upwind first order or Lax-Wendrof The numerical model is capable of simulation of the advection and diffusion processes of suspended materials under the driven action of the depth averaged currents 3) Collection and digitization of the bathymetry in the Halong have been done for an orthogonal mesh with the same mesh size in x and y direction 250 m The wind and tidal data in the near shore area of Quang Ninh also were collected Then, the numerical model was applied to the simulation of wind induced currents resulted by the averaged wind in the sommer and winter; and by the tide Simulation of propagation and distribution of suspended materials initiated at some locations under the depth averaged currents induced separately by wind and tide, and by both wind and tide was also carried out Simulated resuhs showed resonable pictures of the water current and pollution concentration distribution 4) In general, the numerical model developed in this study meets the requirement on the simulation of depth averaged currents and pollution propogation with an acceptable accuracy and the model can be applied to other areas of shallow sea 5) Results of the study were published as an article on the Journal of Science, VNU VI MUC LUC Mdd^u Chuong Md hinh md phdng ddng chay hai chilu 1.1 H6 phuong trinh ca ban va cac di6u kien bien 1.2 Phuong phap sai phan theo so dd CIP cho he phuong trinh nude nong 1.3 Kiem chung md hinh sd'theo cac bai toan ly thuyet va thuc te 1.4 Nhan xet 4 10 16 Chuong Mo hinh mo phdng truydn chat hai chi6u 2.1 He phuong trinh coban vacac didu kien bien 2.2 He phuong trinh dudi dang rdi rac sir dung phuong phap TVD 2.3 Kiem chung mo hinh sd theo cac bai toan ly thuyet 17 17 17 18 Chuong 3: Ung dung md hinh cho vung vinh Ha Long 3.1 Gidi thieu viing nghien cuxi 3.2 Sdlieu sijrdung nghien cuu 3.3 Ket qua mo phdng 26 26 35 35 Ket luan va kien nghi 51 Tai lieu tham khao 53 Phuluc 54 Tom tat cac cong trinh NCKH cua ca nhan (Miu 1) 55 Scientific Project (Mlu 2) 56 Phieu dang ky ket qua nghien cuu KH-CN 57 MODAU O nhidm moi trudng bien hien dang la va'n dl dugc quan tarn ciia hSu het cac qudc gia cd bien tren the gidi Vdi nhu c^u phat trien kinh te hudng vao tiem nang ciia bi^n da cang gia tang ap lire v6 nhilm moi trudng bien ndi chung va d nhilm mdi trudng vung ven bien ndi rieng Chinh vi le dd, viec phat trien kinh te ciia mdt viing ludn dugc ddi hoi dat bdi canh phat tri^n ben vung Trong dd viec danh gia tac ddng ciia boat dpng kinh te dan sinh d6'n thay doi moi trudng ddi hoi phai dugc thuc hidn mdt each nghiem tiic de cd nhirng chien luge phat trien cung nhu qui hoach phat trien dudi sir quan ly tdng th^ cho vira dap ung dugc nhu cau phat trien tat yeu cua kinh te xa hoi, vira dam bao giam thieu dugc tac hai cho moi trudng tir nhien vd'n cd Hiu het cac qudc gia phat trien va dang phat trien cd bien (nhu: My, Nhat ban, Anh, Ha Ian, Mexico, Trung Qudc, v.v.) d6u da va dang thuc hien cac dir an khao sat, danh gia va du bao bien dpng moi trudng nham dua cac canh bao, cac thong tin trp giiip cho he thdng chinh sach, phue vu cho qui hoach phat trien kinh te cua cac viing bi^n va ven bien quan (Peterlin va nnk, 2005 [10]; Lozano va nnk, 2005 [7]; Lau, 2005 [8]) Trong dd viec md phdng va danh gia muc dp Ian truyen chat d nhi6m mot viing nude cu the la mot nhung mat xich quan trpng cac chuong trinh quan ly tich hgp phuc vu cho viec xir dung hop ly tai nguyen va phat trien ben vimg d h^u het cac qudc gia tren the gidi dd c6 Viet Nam Cac cong cu md phdng thuy dpng lire va Ian truydn vat cha't noi chung va cac chat d nhilm ndi rieng (COD, BOD, TSS, ) da va dang dupe rat nhieu trung tam nghien cull cd uy tin tren the gidi phat trien va cho cac san phfa thuong mai cd gia tri cao dap ung dugc cac nhu cau ung dung quan ly va sir dung hgp ly ngudn nude cung nhu trg giiip cho qui hoach phat trien kinh te Cac san pham cong cu phan mem dai dien nhu SMS (ciia My), MINLAKE (cua My), MIKE 21 (cua Dan mach) Dd la nhung cong cu hiru ich de ting dung tren thuc te nude ta cac chuyen gia ve mo hinh hoa va md phong loan hpc tren may tinh cac qua trinh thuy dpng lire hpc va Ian truyen cac cha't d nhidm da va dang dan tiep can dugc nhimg tuu khoa hpc hien dai tien tien tren the gidi He phuong trinh nude nong phi tuyen hai chi6u da dugc nhi^u tac giai sir dung de md phdng bai toan ddng chay hai chi^u trung binh theo dp sau nhu bai toan thuy trieu, bai toan ddng chay gio Tuy nhien cac tac gia chu yeu don gian hoa he phuong trinh dang tuyen tinh hoa nhu Nguyln Ngpc Thuy (1969), Dang Cong Minh (1975) [1, 2] hoac giai phuong trinh nude nong d^y dii nhung vdi cac so dd toan don gian nhu UPWIND bac hoac sai phan trung tam cd dp chinh xac khong cao Ddi vdi bai toan truyen tai vat chat thi va'n dt khd hon la vdi cac so dd toan don gian thi viec ton tai cac nghiem sd khong that (thi du nghiem am, hoac nhiiu sd) la va'n de thudng gap Do vay viec phat trien mdt md hinh md phdng ddng chay va chuydn cha't hai chi6u dang tin cay cho viing ven bd la v^n di cin thiet phuc vu cho nghien cuu ciing nhu qui hoach moi trudng cac viing nude ven bd Muc tieu ciia dd tai la phat trien va ung dung cong cac so dd sd tien tien cd dd chinh xac cao vao xay dung mo hinh sd md phdng dugc ddng chay va truyen chat hai chi^u viing ven bd dudi tac dpng ciia gid va thuy tridu Sau dd ung dung tinh toan thii nghiem cho viing vinh Ha Long Ngoai phSn md ddu, ket luan va kien nghi, tai li6u tham khao, npi dung bao cao cua d^ tai dugc trinh bay ba chuong: Chuong : Md hinh ddng chay hai chidu Trong chuong trinh bay co sd ly thuyet, he phuong trinh xua^t phat, cac dieu kien bien ciia md hinh md phdng ddng chay hai chi6u cd tinh den anh hudng ciia gid Trinh bay tdm luge phuong phap sai phan cd dd chinh xac cao CIP (Constrained Interpolation Profile); phuong phap sai phan an dang SMAC va each ung dung rdi rac hoa giai he phuong trinh nude nong phi tuyen Cac bai toan kiem nghiem md hinh dudi dang ly thuyet va thuc tiln cung dupe trinh bay nham khang dinh tinh on dinh, dung dan cua md hinh so Chuong : Md hinh truyin cha't hai chidu Trong chuong trinh bay he phuong trinh md ta qua trinh truyen tai va khuech tan vat chat hai chieu dudi tac dpng cua ddng chay nen Md ta cac dilu kien bien Trinh bay phuong phap sai phan dang TVD cd dp chinh xac bac hai va ling dung rdi rac hoa giai he phuong trinh xuat phat Md hinh sd dugc kiem chung thdng qua cac bai toan co ban cd nghiem giai tich nham kiem nghiem kha nang mo phdng qua trinh truydn tai vat chat ciia md hinh Chuong 3: Ung dung md hinh cho viing vinh Ha Long Chuong trinh bay tom lupe viing nghien cun, di6u kien gid, thuy treu, dia hinh Lira chpn cac sd lieu dau vao de md phdng Cac ket qua mo phdng Ian truyen vat chat gia dinh khong bi phan huy theo thdi gian theo cac kich ban dudi tac dpng ciia gid va thuy trieu Cac nhan xet v6 ket qua nghien cun De tai dupe hoan vdi sir hd trp kinh phi tir phia Dai hpc Qudc gia (DHQG) Ha Ndi, sir giiip ciia Ban Khoa hpc & Cong nghe, DHQG, Phong Khoa hpc Cdng nghe, trudng Dai hpc Khoa hpc Tir nhien, sir ung ho nhiet tinh ciia Hoi ddng Khoa hpc Trai da't, DHQG HN, ciia Ban chu nhiem khoa Khi tugng Thuy van va Hai duong hpc, cung nhu su ddng gop y kien eiia cac nha khoa hpc, cac ddng nghiep tham gia di tai va tap the can bd cua Bp mdn Hai duong hpc Nhan day chiing tdi xin chan cam on nhimg giiip qui bau dd CHUONG 1: MO HINH MO PHONG DONG CHAY HAI CHltv LI H £ PHUONG TRINH CO BAN VA CAC Difiu KI$N BifiN He phuong trinh xuat phat diroc sir dung phat trien md hinh la he phuong trinh nude ndng phi tuyen, bao gom : Phuong trinh bao toan dpng lugng du du du dC — + u— + v fv = -S^Adt dx dy dx p{h + dv dv dv ^ dC — + u — + v — + fu = -2 -^- + dt dx dy dy (LI) k'rlh^^^ (1.2) p{h + Phuong trinh bao toan khdi lupng dC , d(h + Ou , d(h + Ov dx dt dd dy ( du^ I dx.1 dy dx d ( +— dy dx (1.3) du] (1.4) dv] (1.5) ( •' ^ 'J, — w,v: cac phan van tdc trung binh dp sau theo cac hudng ox va oy\ f: tham sd Coriolis; (^ : dao dpng muc nude; p : mat dp nude; h : dp sau nude; u^: he sd nhdt rdi T^Ty la cac phan ung suat gio tren mat theo true x va true y, chung dupe xac dinh theo cdng thuc sau: T: =P^C,^\W\W/, T].=P^C\W\W^,^ (1.6) dd p^ la mat dp khong khi; w la vector van tdc gid He sd C, dugc xac dinh theo Smith & Banke (1975) nhu sau : (0.63 + 0.066|W^|)xlO^ C (2.28+ 0.033(1^^^1-20)) X 10" |fF|2t)mls Cac phan ung suat ma sat day dupe xac dinh theo cdng thuc sau: CLU^JU + V h r^ I 2 Tj = p Cf,vyju + v^ He sd ma sat Q dugc xac dinh theo cdng thuc C, (h + cy Ai^ 0.0264 (1.7) Phitog Dang Hitfo 46 Bing Cac bien long va cac hing so dieu hoa Vido M2 Kinh S2 K1 01 H(cm) G" H(cm) G' H(cni) G- H(cm) G- Ddi Loan 23*33' 119^3' 87 326 23 15 25 278 21 236 Basi 2r58' 120°45' 26 188 11 214 23 262 18 224 108*'59* 26 125 199 14 66 16 10 Malaea ^''\v Hinh trinh bay tnl6ng v4n toe dong trung binh tai luc 10 gid ng4y thang nam 2004 Hinh Ket qua mo ph&ng tn/dng van toe trung binh tren biln Dong luc 10 gid 1/1/2004 Ke-t qua mo phong dao dong m^c nUdc tai Hon Dau dUdc so sanh vdi so heu thUc tai tram Hon Dau thang nim 1988 Hinh trinh bay dao dong mUc n u ^ trilu tai khu vuc Hon Dau thang nim 1988 Tren hinh ve ta th^y dao dpng mUc nuac tinh toan kha p h h ^ v6i so^heu thuc Dao ^^^ ^^.^^'^^^J^^ ^ ' ^ ^ nhat trilu diu va xu the^ thay do'i cua dao dong mUc nU6c trnh toan r&t phU hop ven thu J ' T r e n M n ^ vg so sanh ta cOng tha'y c6 sU sai khac nho vl bi.n d tneu g^Ua mo rh6ng so-va so^Ueu thUc do, nguyen nhan co thi giai thich v, t n lay ket qua m h ng sai lech kha I6n so v6i vi tri c6 so^ heu vi lua la 1/5 dp va hdn nUa lUP mo Phuong phiip sai phan Ai cho bai toto thuy tri^u bidn ven 47 phong kha tho nen chUa mo ta dUdc het nhiJng anh hudng phi tuy^n cua viing nUdc nong va sU thay doi dia hinh phtic tap len dao dong mUc nUdc trieu •T£nh to6n Thuc • I -1 -2 50 100 150 200 250 300 350 400 Thdi ^an (gijy) Hinh Dao dong mUc nUdc tai khu vi/c Hdn Dau (b^t diu tCrO gid ng^y 01/1/1988) Mo hinh cung dUdc tinh toan thuy trieu thii cho vtog vinh Ha Long Gidi han vtog vinh dUdc chon la toan bo vung vinh Ha Long va Bai T i Long Vung ap dung tfnh toan thtt nghiem dUdc chia eac lUdi vuong vdi bUdc lUdi 250m Cac bien long dUdc sii dung la eac dao dong mUc nUdc trich tii mo hinh Idn tinh cho viing vinh Bac Bo vdi ludi la 1/30 Tai cac diem bien tren mien tinh nho khong triing vdi dilra lUdi Idn c6 so* lieu trich thi dUdc npi suy tuyen tinh theo cac so Heu tren cac diem ciia lUdi 16n Tren hinh 7, trinh bay ddng chay tai hai thdi diem nUdc xuong va nUdc len vung vinh Ha Long (Delta X-250 m) Hinh Trudng van toe d6ng chiy vung vjnh Ha Long luc pha trieu xuong Phiing Dang Hi^u imJk X (DeKaz-250m) Hinh TrUdng van toe ddng chSy vung vjnh Ha Long luc pha trieu len Tren cac hinh va ta tha'y toe dong chay trieu tai eac khe hep giiJa cac dao va tai nhiing vung nUdc nong la kha Idn K^t qua nhin chung kha phu hdp vdi qui luat cua dong chay tac dong ciia thuy trieu Do nguon so' Ueu thUc ve dong chay tai cac eo khe dao d tai viing Ha Long khong c6 nen viec so sanh giiia ket qua mo phong va so lieu thUc te chUa thUc hien dUdc Tuy nhien ket qua bude d^u cho tha'y phUdng phap sai phin ^n trinh bay bai bao la ra't hieu qua va de dang ap dung Dieu chiing to mtic on dinh ciing nhU kha nang ap dung vao thUe te ciia phUdng phap sai phan trinh bay cl tren Ket luan • PhUdng phap sai phan in trinh bay tren ra't ddn gian, de sii dung va co chinh xac ciing nhu on dinh kha tin cay Ket qua sii dung phUdng phap cho mot so bai toan 1^ thuyet ciing nhu thUc te trinh bay d tren da chiing to kha nang ting dung thuc te ciia mo hinh toan phat triln nghien ciiu la ra't hien thUc Tuy nhien de di den hoan thiSn ph^n mim iing dung mo hinh toan cin dUdc kiem chiing them niia vdi cac so heu thUc te Ldi cam an Noi dung bai bao trinh bay cl day la mot phfin cua de tai nghien ciiu khoa hoc ca'p Dai hoc Quoc gia tac gia chii tri Nhan day tac gia giii ldi cam dn tdi phong khoa hoc trudng Dai hoc KHTN va Dai hoc Quoc gia Ha Noi da trd giiip kinh phi va tao dieu kien cho tac gia thUc hien de tai ma so QT-06-34 phuong phip sai phan fc cho Mitodnthuy tri^u bi^n ven q TAI L I E U T H A M KHAO [1] Pham v a n Ninh va nnk, Bien Dong Tap 11: Khi tuang thuy van dong lUc bien, ChUdng trinh dieu tra nghien ctiu bien c^p nha nu6c KHCN-06 (1996-2000), NXB Dai hoc Quoc gia Ha N6i, Ha N o i , 2003, 565 tr [2] Nguyin Thi Viet Lien, M6 hinh so tri bai toan thuy trieu bien Dong, LuAn an Ph6 tien si khoa hoc, Vien Cd hoc, 1997 [3] Yabe T., T Aoki., A universal solver for hyperbolic equations by cubic-polynomial interpolation, Comp, Phys Communications, Vol 66, Elsevier, 1991, p 219-232 [4] Hieu P D., K Tanimoto, C T Vu., Numerical simulation of breaking waves using a two-phase flow model Int J Applied Math Modeling, Elsevier, Vol 28, No.ll(2004) p.983-1005 [5] Admiraty Tidal Tables, Pacific Ocean And Adjacent Seas, Vol 3, 1971 VNU JOURNAL OF SCIENCE Nat, Sci., & Tech., T.XXII No2B AP 2006 AN IMPLICIT FINITE DIFFERENCE SCHEME APPLIED TO THE TIDAL SIMULATION MODEL FOR SHALLOW SEAS P h u n g D a n g Hieu Department of Hydro-Meteorology & Oceanography College of Science, VNU The p a p e r p r e s e n t s a n impHcit finite difference scheme for t h e shallow water equations a n d a n application for the tidal simulation of a shallow sea The implicit finite difference method used here is somewhat similar to the SMAC method, which was employed for t h e discretisation of t h e Navier-Stokes equations for the simulation of breaking waves in t h e previous studies of the author The advantage of the method is t h a t it is stable a n d accurate enough for a simulation of discontinuity when it combined with C I P scheme Firstly, a numerical model is developed using t h e implicit scheme and verified a g a i n s t the analitical solutions in t h e case of dam-break problem Then, t h e model is applied to the simulation of tidal waves in the Bien Dong Sea and in the Halong bay The n u m e r i c a l results show t h a t t h e proposed model is capable of simulation of tide in shallow seas TRircWVG DAI HQC KHOA HOC TlT NHIEN TRAN Dire TRir Sir DUNG PHirONG PHAP CIP vA HE PHl/ONG TRINH N i r d c NONG PHI TUYEN MO PHONG NUCiC D A N G DO B A O Chuyen nganh: Hai dirang hoc Ma s6: 1.07.07 LUAN VAN THAC SY KHOA HOC NGLTCII Hl/dNG DAN KHOA HOC: TS PHUNG D A N G HIEU Ha Noi - 2006 TOM TAT CAC CONG TRINH NCKH CUA CA NHAN (LIEN QUAN DEN DE T A I ) Nganh: Hai duorng hoc; Chuyen nganh: Vat ly va moi truong bien Ho vd ten: Phung Dang Hieu Nam: 2006 Ten bdi bao: Phuong phap sai phan an cho bai toan thuy trieu bien ven Ten tap chi: Tap chi Khoa hoc, T.XXII, No 2B AP, 2006, DHQG Ha Noi tr 39- 49 Tom tat: Bai bao trinh bay mot phirong phap sai phan an cho he phuong trinh nude nong hai chieu va irng dung cho bai toan thuy trieu bie'n ven Phuong phap sai phan in su dung Of day dugc thuc hien tuong tu nhu phuong phap SMAC da dugfc ung dung cho he phuong trinh Navier-Stokes cong trinh nghien cim song ngan cua tac gia Dac die'm cua phuong phap la lam tang tfnh 6n dinh cua mo hinh sd va 6k dang giai bang cac phuong phap giai lap hien dai Thanh phan binh luu phi tuyen he phuong trinh nude nong dugc xa'p XI bang so CIP (constraint interpolation profile) co chinh xac bac ba va phuong phap giai lap BiCGSTAB dugc ap dung Trudc tien mo hinh sd dugc kiem nghiem vdi mot sd trudng hgp co nghiem giai tich vdi bai toan va dap, sau mo hinh sd dugc ap dung tinh toan thu cho truyen trieu tren bien Dong va tren vung bie'n vinh Ha Long Quang Ninh Ket qua ban dau cho thay phuong phap sai phan an trinh bay d day co tinh 6n dinh cao, dt ap dung va co kha nang ap dung tot cho bai toan thuy trieu bie'n ven Title: An Implicit Finite Difference Scheme Applied to the Tidal Simulation Model for Shallow Seas Journal: Journal of Science T.XXII, No.2B AP, 2006, Vietnam National University Hanoi, pp 39-49 Summary: The paper presents an implicit finite difference scheme for the shallow water equations and an application for the tidal simulation of a shallow sea The implicit finite difference method used here is somewhat similar to the SMAC method, which was employed for the discretisation of the Navier-Stokes equations for the simulation of breaking waves in the previous studies of the author The advantage of the method is that it is stable and accurate enough for a simulation of discontinuity when it combined with CIP scheme Firstly, a numerical model is developed using the implicit scheme and verified against the analitical solutions in the case of dam-break problem Then, the model is applied to the simulation of tidal waves in the Bien Dong Sea and in the Halong bay The numerical results show that the proposed model is capable of simulation of tide in shallow seas 55 SCIENTIFIC PROJECT Branch: Oceanography; Project category: Ministry level Title: Development of a numerical 2-D model based on the shallow water equation and convection-diffusion equation for the application to the Halong Bay Code: QT.06.34 Managing Institution: Vietnam National University, Hanoi Implementing Institution: University of Science, VNU, Hanoi Collaborating Institutions: None Coordinator: None Key implementors: Dr Phung Dang Hieu; M.Sc Pham Hoang Lam; B.Sc Tran Due Tru Duration: Jan 2006-Dec 2006 Budget: 20.000.000 VND 10 Main results: • Results in science and technology: Contribution to the improvement of knowledge on marine hydrodynamie and environmental modeling Provide a good tool for the researches on hydrodynamie processes resulted by tide, wind and on distribution of pollution in shallow water • Resuhs in practical application: Will be used as a good reference material for the further researches in the field of coastal hydrodynamics and environmental modeling • Results in training: master thesis • Publications: article published in the Journal of Science, VNU Hanoi 11 Evaluation grade (if the project has been evaluated by the evaluation committee: excellent, good, fair): 56 PHIEU DANG KY KET QUA NGHIEN CUU KH-CN Ten d^ tai: My dung mo hinh mo phong ddng chay va truyen chat hai chieu I'mg dung cho viing vinh Hq Long Ma so: QT.06.34 Ca quan chu tri de tai: Trucmg Dai hoc Khoa hoc Tu nhien Dja chi: 334 Nguyln Trai, Thanh Xuan, Ha Noi Dien thoai: (04) 8581419 Co quan quan ly de tai: Dai hoc Qudc gia Ha Noi Dia chi: 144 Duofng Xuan Thiiy, Cau Giay, Ha Npi Dien thoai: Tong kinh phi thuc chi: 20.000.000 dong (Hai muoi trieu dong) Trong do: - Tir ngan sach Nha nude: 20.000.000 dong - Kinh phi ciia trucmg: 0,0 - Vay tin dung: 0,0 - Vontuco: 0,0 - Thu h6i: 0,0 Then gian nghien curu: 12 thang Thdi gian bat dau: 1-2006 Then gian ket thuc: 12-2006 Ten cac can bo phdi hup nghien curu: 1) ThS Pham Hoang Lam 2) CN Tran Du-c Trii: (hoc vien cao hoc) S6' dang ky de tai So chumg nhan dang ky ket qua nghien cihi: Ngay: 57 Bao mat: a Pho bien rong rai: V b Pho bien han che: c Bao mat: Tom tdt k^t qua nghidn curu: 1) Da ling dung phucmg phap CIP (Constraint Interpolation Profile) hien dai cd chinh xac bac ba va phat tri^n mot phuong phap sai phan an vao viec sai phan hoa giai he phucfng trinh nude ndng phi tuyen cho bai toan dong chay Viet cong chuong trinh may tinh md phong dong chay sinh gio va thuy trieu Tinh toan thu nghiem va kiem chiing m6 hinh cho cac dieu kien ly thuyet co nghiem giai tich va kie'm chung tnrofng hgp dieu ki^n thirc te ciia bai toan truyen trieu bie'n Dong Ket qua cho thay phuong ph^p CIP sir dung d day cho bai toan mo phong dong chay la rat hap ly va co chinh xac kha cao Mo hinh chiing to kha nang ling dung tot cho cac vung bien nong ven bd 2) Da xay dung cong chuong trinh truyan tai va khuech tan vat chat su dung phuong phap the tich him han vdi so TVD (Total Variation Diminishing) co chinh xac bac hai Viet cong chuong trinh may tinh va kiem chiing mo hinh cho cac bai toan ly thuyet cd nghiem giai tich Ket qua cho thay phuong phap TVD sir dung d day co chinh xac cao hOn nhieu so vdi mot sd phuong phap sai phan truyen thdng hay duoc su dung nhu UPWIND bac mot hay Lax-Wendrof Mo hinh cho phep mo phong qua trinh truyen lai khue'ch tan vat cha't dudi tac dong ctia dong chay trung binh 3) Da thu thap va sd hoa dia hinh vung vinh Ha Long tren ludi vuong gdc toa d6 De Cac vdi budc ludi 250 m x 250 m Thu thap va phan tich nhung sd lieu ve toe gio va hudng gio vung ven bd Quang Ninh, sd lieu hang sd dieu hoa vung vinh Ha Long, ling diing mo hinh mo phong he thdng dong chay gay gio trung binh miia Dong va mua He va tac ddng cua thuy trieu cac chu ky trieu len va xudng Idn Mo phong lan truyen vat cha't dugc thuc hien tren nen dong chay t6ng hgp ca gio va thuy trieu vdi gia thiet cha't chi thi ben vung khong phan huy va co sir cd mo trudng thai tai mot sd diem gia dinh Cac ke't qua cho cac biic tranh phan bd trudng dong chay va phan bo ndng chat nhiem co tac dong ciia cac yeu td dpng lire la gio va thuy trieu 4) Noi chung mo hinh xay dung de tai da dap ling dugc yeu cau ve mo phong dong chay va truyen chat vdi dp chinh xac dang tin cay va co kha nang ling dung de nghien ciiu cho cac vung bien ven bd khac Ket qua ling dung li'nh toan cho vung vinh Ha Long, Quang Ninh cho thay ro kha nang iing dung ciia mo hinh da phat trien 5) Dang dugc mot bai bao tren tap chi khoa hpc Dai hpc Quoc Gia Kien nghi ve qui mo va ddi tugng ap dung nghien cuu: Nhimg ket qua nghien ciiu ciia de tai budc dau cho tha'y mo hinh toan phat trien de tai mo phong bai toan dong chay va truyin chat hai chieu co kha nang iing dung tren thuc tiln Ke't qua cua de tai co the la tai lieu tham khao cho cae nha nghien cuu vl dpng lire hpc va moi truong bien cOng nhu phuc vu tot cho cong tac dao tao dai hpc va thac sy San pham ciia dl tai cung tao mot phucmg phap va cong cu him ich cho viec mo phong va nghien ciiu dong chay, lan truyen Gjia't.6 fihu^' viing ven bien Tcri • • ',; -'^r C-';Gf;GHt Chii nhiem de tai Ho ten Thii trutfngcaquan chu t r i de tai Chii tich Hoi dong danh gia chinh thiJc Thii trirong co ;*