BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI NCS LÊ THỊ THÁI NGHIÊN CỨU HIỆN TƯỢNG XÂM THỰC BAO QUANH CHÂN VỊT TÀU THỦY H Hà nội – 2013 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI NCS Lê Thị Thái NGHIÊN CỨU HIỆN TƯỢNG XÂM THỰC BAO QUANH CHÂN VỊT TÀU THỦY Chuyên ngành : Cơ học chất lỏng Mã số : 62442201 H PGS.TS Lê Quang TS Lê Thanh Tùng Hà nội - 2013 LỜI CẢM ƠN tài “Nghiên cứu tượng xâm thực bao quanh chân vịt tàu thủy” PGS.TS Lê Quang TS Lê Thanh Tùng TCc, Nghiên cứu sinh Lê Thị Thái LỜI CAM ĐOAN “Nghiên cứu tượng xâm thực bao quanh chân vịt tàu thủy” chép, gian Lê Thị Thái MỤC LỤC Phần mở đầu Chương 1: Tổng quan 1.1 5 1.1.1 1.1.2 1.1.3 1.1.4 13 13 16 18 23 1.3 25 1.4 27 Chương 2: Nghiên cứu, khảo sát dòng chảy bao chân vịt tàu thủy 29 2.1 Chân vc tính làm vic ca chân vt tàu thy 30 2.2 lý thuyt tính tốn 2.2.1 Mơ t v 32 32 2.2ng qt v dịng th 2.2u kin biên 34 35 2.2.4 Cu kin biên b mt xâm thc 37 2.3 Kho sát dòng chy bao chân vt bng p 2.3.1 Mơ hình tính tốn 2.3.2 Mơ hình hình hc min tính tốn 2.2.3 Kt qu tính tốn s 2.4 Kt lun Chương 3: Nghiên cứu, mơ dịng xâm thực qua chân vịt 3.1 Hing xâm thc thit b y 3.2 Nghiên cu s v mơ hình xâm thc 2D 3.2.1 S d 3.2.2 S d 45 45 48 48 55 56 56 57 57 65 3.3 Nghiên cu s v mô hình xâm thc 3D 3.3.1 Tính tốn dịng qua chân vt tàu thy 3.3.2 Tính tốn dịng qua h thng chân vt bánh lái 3.4 Kt lun Chương 4: Khảo sát thực nghiệm đặc tính xâm thực chân vịt 75 75 80 84 85 85 85 88 92 94 94 95` 4.3 P 96 4.3 96 4.3 4.3 99 103 Kết luận kiến nghị 104 106 PH L : Code chng trình tính tốn s v mơ hình xâm th 2D b phng pháp LES 109 DANH MỤC CÁC KÝ HIỆU A Ac Cw C x, Cy, Cz C L, CD C l, C d Cp D F(s1,s2,s3 ,t) F G(x,x) J KT KQ M : Vector mômen O Q = Mx R C = SB SC SBC : Sw theo Swc T (u,v) V = (vx ,vy ,vz VC VO VS : Vg Vm Vw(x,y,z) VRef VS1, VS2 , VS3 VU1, VU2 , VU3 X = (X,Y,Z) g g : Hàm ngnh cho b mt siêu xâm thc (i,j,k) : Vector vn t xâm thc tham chiu cho vt th c nh lC : Chiu dài xâm thc n T vòng quay P : Áp sut PV : Áp su PRef : Áp sut tham chiu r : Vector khong cách rf : Khong cách gim trng cong tm tuyn tính tham chiu SO : Chiu dài cung tm tách xâm thc (S1 ,S 2,S 3) : H t không trc giao t : Thi gian * t : Th (t1 ,t2,t3 ) (u1,u 2,u3 ) : H t tham chiu trc giao gn vi b mt dng cong tuyn tính V : Vn tc nhiu lon X = (x,y,z) : H t tham chiu cho vt th c nh : H t cc : Bin, loi kích c tm (2D) - : Hồn l ng dịng bên ngồi : V 0 t2 : Trng ng riêng ca c Vl Vv : Th tích pha lng : Th tích pha hi GIẢI THÍCH CÁC TỪ VIẾT TẮT CFD : Computational Fluid Dynamics ITTC : Internationl Towing Tank Committee FDM FEM SEM FVM BEM : Finite Difference Method : Finite Element Method : SpEctral Method : Finite Volume Method : Boundary Element Method RANS LES DES DNS VM : Reynolds Average NavieStokes : Large Eddy Simulation : Detached Eddy Simulation : Direct Numerical Simulation : Vortex Method AFD EFD VOF KBC DBC : Analysis Fluid Dynamics : Experimental Fluid Dynamics : Volume Of Fluid : Kinetics Boundary Conditions : Dynamics Boundary Conditions “UBlendingFactor", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE fvc::interpolate(gamma) # include "createPhi.H" Inforho2(); const dimensionedScalar& pSat = twoPhaseProperties->pSat(); // Need to store rho for ddt(rho, U) volScalarField rho IOobject "rho", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE gamma*rho1 + (scalar(1) - gamma)*rho2, gamma.boundaryField().types() rho.oldTime(); label pdRefCell = 0; scalar pdRefValue = 0.0; setRefCell(pd, mesh.solutionDict().subDict("PISO"), pdRefCell, pdRefValue); Info