Các phương trình được sử dụng
Phƣơng trình cho phần trăm thể tích hơi nƣớc (water vapor)
(3.5)
Phƣơng trình thể hiện tỷ lệ thay đổi pha (1 ) 3 2( ) 3 v l v B l p p R R p (3.6)
với α tỷ lệ phần trăm pha hơi
33 3 4 3 4 1 3 b B b B n R n R
Trong đó: khối lƣợng riêng hỗn hợp (1 ) l v
RB bán kính bọt khí, nb là mật độ bọt khí, các đại lƣợng cụ thể khác có trong [11] Các phƣơng trình 3.1, 3.2, 3.3, 3.4 đƣợc sử dụng để giải bài toán mô hình không xâm thực, với mô hình xâm thực còn có thêm phƣơng trình 3.5 và 3.6.
Thông số vật lý chất lỏng mô phỏng
Thông số Water-liquid Water-vapor
Khối lƣơng riêng (kg/m3
) 998.2 0.5542
Độ nhớt động lực học (kg/m.s) 0.001003 1.34e-05
Bảng 3.1: Thông số vật lý của pha nước và pha hơi
( ) .( ) .( ) v v l v D u t Dt
Bộ môn Máy và Tự động Thủy khí-Viện Cơ khí Động lực-Đại học BKHN 45
Hình 3.5: Mô hình lưới 2D (a)Mô hình lưới trên TD3
(b)Lưới tại sát cánh công tác
Thiết đặt cơ bản trong Ansys Fluent
Thiết đặt Chế độ không xâm thực Chế độ xâm thực
Thuật giải Pressure-Base, phụ thuộc thời gian: Transient
Pressure-Base, Transient (unsteady)
Mô hình rối K-Eps Realizable, không trao đổi nhiệt
K-Eps Realizable, không trao đổi nhiệt
Đa pha Không Multiphase-Mixture
Áp suất bão hòa Pbh=0.03493at, nb=10e13
Biên Pressure outlet (Pout=0 at) (theo áp suất dƣ)
Non-slip wall (tƣờng biên không trƣợt)
Biên Periodic tại mặt cắt
Pressure outlet (Pout=1 at)(theo áp suất tuyệt đối) Non-slip wall
Biên Periodic tại mặt cắt
Nội suy Áp suất: standard
Mô men, K,E… chọn second order up wind.
Áp suất: standard
Mô men, K,E… chọn QUICK
(a)
(b) )
Bộ môn Máy và Tự động Thủy khí-Viện Cơ khí Động lực-Đại học BKHN 46
Điều kiện mô phỏng
Tính theo áp suất dƣ Operating Pressure=1 at
Theo áp suất tuyệt đối Operating Pressure=0 at
Độ chính xác 1e-05 Vapor 1e-04, còn lại 1e-05
Time step 1e-03 s 1e-03 s