Đồ án kỹ thuật hạt nhân

Đồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhânĐồ án kỹ thuật hạt nhân

LỜI CẢM ƠN Đầu tiên, em xin chân thành cảm ơn đến TS Nguyễn Văn Thái giúp đỡ đưa lời khuyên sâu sắc để em thực đồ án tốt nghiệp Nhờ hướng dẫn tận tình thầy, em thực mở mang kiến thức lĩnh vực thủy nhiệt hạt nhân, chuẩn bị cho trình nghiên cứu tới Tiếp theo, em xin cảm ơn tới thầy cô Viện Kỹ thuật hạt nhân Vật lý môi trường góp ý động viên em q trình em thực đồ án Cuối cùng, em xin gửi lời cảm ơn tới gia đình bạn bè – người ln ủng hộ nhiệt tình để em vượt qua thời điểm khó khăn LỜI CAM ĐOAN Tôi xin cam đoan luận văn tốt nghiệp tơi tự nghiên cứu, thiết kế thu hoạch kết hướng dẫn TS Nguyễn Văn Thái Để hồn thành đồ án này, tơi sử dụng tài liệu ghi phần tài liệu tham khảo, ngồi tơi khơng sử dụng tài liệu khác mà không ghi Nếu sai, tơi xin chịu hình thức kỷ luật theo quy định Sinh viên thực Nguyễn Ngọc Quý MỤC LỤC DANH MỤC HÌNH VẼ DANH MỤC KÝ HIỆU DANH MỤC BẢNG BIỂU MỞ ĐẦU Cơ sở lý thuyết động học dòng hai pha kênh dẫn đứng 1.1 Dòng chảy hai pha kênh dẫn đứng 1.2 Phương pháp mơ hình hóa q trình động học dòng hai pha 1.3 Mơ hình hai dòng chảy cho phương trình bảo tồn 11 1.4 Các mơ hình tương tác động học 13 1.4.1 Drag force 13 1.4.2 Lift force 17 1.4.3 Wall lubrication force 21 1.4.4 Turbulence dispersion force 23 1.4.5 Virtual mass force 25 Phương pháp mơ động học dòng hai pha kênh dẫn đứng 27 2.1 Phương pháp thể tích hữu hạn (Finite Control Volume) 27 2.1.1 Rời rạc miền không gian thời gian 29 2.1.2 Sắp xếp lưu trữ biến 30 2.1.3 Rời rạc phương trình 31 2.1.4 Rời rạc theo khoảng thời gian 33 2.1.5 Phương pháp giải hệ phương trình đại số tuyến tính 34 2.2 Chương trình tính tốn động học dòng chảy mã nguồn mở OpenFOAM 36 Áp dụng phương pháp giải vào trường hợp thực nghiệm trường đại học PURDUE sử dụng chương trình tính tốn mã nguồn mở OpenFOAM 39 3.1 Hệ thực nghiệm đo đạc q trình động học bọt điều kiện khơng trao đổi nhiệt đại học PURDUE 39 3.2 Tiền xử lý tốn mơ hệ thực nghiệm 41 3.3 Sơ đồ thuật giải tốn động học dòng hai pha kênh dẫn đứng chương trình OpenFOAM 43 3.3.1 Vòng lặp hiệu chỉnh tăng tiến theo thời gian 43 3.3.2 Vòng lặp hiệu chỉnh áp suất-vận tốc 44 Kết mô nhận xét 47 Kết luận 52 TÀI LIỆU THAM KHẢO 53 PHỤ LỤC 54 DANH MỤC HÌNH VẼ Hình 1: Các chế độ dòng chảy truyền nhiệt kênh dẫn đứng hình ống [3] Hình 2: Đường cong sơi chế độ trao đổi nhiệt Hình 3: Dòng phân tách (dispersed flow) dịch chuyển lên phía kênh dẫn đứng (pha phân tách pha khí) Hình 4: Sơ đồ phương pháp chung mơ tả tượng vật lý sử dụng mơ hình tốn học Hình Mối liên hệ tượng vật lý, mơ hình vật lý liệu thực nghiệm 10 Hình 6: Các tương tác lực tác dụng lên bọt khí dòng chảy có shear thấp [2] 15 Hình 7: Đường cong hệ số drag xậy dựng qua thực nghiệm Newton, Stoke qua phương pháp số[2] 16 Hình 8: Lift force gây shear dòng chảy lên khối cầu [2] 18 Hình 9: Sự xuất xốy nước khơng bọt khí di chuyển dòng chảy có vận tốc Ub [2] 19 Hình 10: Lift force tác dụng lên bọt khí cầu bọt khí bị biến dạng theo hướng khác [2] 19 Hình 11: Đường cong lift force với (chấm đen) (tam giác đen) [2] 21 Hình 12: (a) Steamline dòng chất lưu qua hai trụ với lưu thơng đặt trước (b) Đóng góp wall lubrication force vào phân bố void fraction [6] 22 Hình 13: Ảnh hưởng lực turbulence dispersion lên phân bố mật độ pha khí [6] 24 Hình 14: Thí nghiệm Serizawa đo độ lệch áp suất qua khe thắt ống venturi [8] 26 Hình 15: Định nghĩa lực virtual mass 26 Hình 16: Thể tích điều khiển hữu hạn khơng gian [10] 28 Hình 17: Sơ đồ phương pháp rời rạc thể tích hữu hạn 29 Hình 18: Rời rạc miền giải [2] 30 Hình 19: Hai cách xếp lưu trữ biến: a) trường thể tích; b) trường mặt [2] 31 Hình 20: Sơ đồ lặp Jacobi 35 Hình 21: Tổng quan cấu trúc OpenFOAM 36 Hình 22: Cấu trúc thư mục “case” OpenFOAM 38 Hình 23: Hệ đo đạc thực nghiệm trường đại học PURDUE (2001) [11] 40 Hình 24: Lưới hóa hình học không gian dạng hai chiều 41 Hình 25: Sơ đồ vòng lặp giải cho bước thời gian 44 Hình 26: Phân bố bọt dọc theo đường kính ống vị trí z/D = 6.00 53.5 49 Hình 27: Sự thay đổi vận tốc riêng pha lỏng theo bán kính vị trí vị trí z/D = 6.00 53.5 50 Hình 28: Sự thay đổi vận tốc riêng pha khí theo bán kính vị trí vị trí z/D = 6.00 53.5 51 DANH MỤC KÝ HIỆU Ký hiệu Clift Cdrag Cwall CTD Đơn vị - g m N m/s - M N / kg  kg / m2 / s Cwm DS F j Ý nghĩa Hệ số lực lift Hệ số lực drag Hệ số lực tường Hệ số lực turbulent dispersion Hệ số lực virtual mass Đường kính Sauter Lực Vận tốc riêng pha Gia tốc trọng trường Lực tác dụng lên mặt p S - q" W / m2  kg / s3 t T U z/D s K °C m/s -   -  kg / m3  N / m2  kg / m / s Pa kg / m3 / s biên đơn vị thể tích Áp suất Nguồn sinh thể tích hữu hạn Thơng lượng nhiệt Thời gian Nhiệt độ Vận tốc Tỷ lệ vị trí đo đạc độ dài ống Hệ số thể tích pha Tốc độ trao đổi khối lượng đơn vị thể tích Mật độ khối Tensor độ nhớt Ký hiệu số dưới: Ý nghĩa Pha khí Mặt biên phân cách hai pha Pha lỏng Chỉ số pha Chảy rối Ký hiệu g i l k t DANH MỤC BẢNG BIỂU Bảng 1: Các số hạng rời rạc 33 Bảng 2: Thông số ban đầu hệ thực nghiệm PURDUE (2001) 41 Bảng 3: Điều kiên nguồn biên cho đại lượng 42 Bảng 4: Thông số vận chuyển (transport) nước khí đặt cho tốn mơ 42 Bảng 5: Kết so sánh giá trị trung bình void fraction theo diện tích thể tích điều khiển dọc theo bán kính thực nghiệm mơ 51 MỞ ĐẦU Ngày nay, ngành lượng hạt nhân phải đối mặt với nhiều thách thức sau cố tai nạn nhà máy điện hạt nhân Fukushima năm 2011 Đặc biệt lĩnh vực thủy nhiệt hạt nhân, yêu cầu đặt cần phải có hiểu biết sâu tượng xảy lò phản ứng, từ tạo bước tiến vấn đề nghiên cứu phân tích an toàn nhà máy điện hạt nhân Một số tượng quan trọng phân tích an tồn thủy nhiệt q trình sơi đối lưu (convective boiling) động học bọt dịch chuyển kênh dẫn đứng bó nhiên liệu nhà máy điện hạt nhân [3, 4] Quá trình quan trọng hoạt động bình thường hay việc dự đoán chuyển tiếp xảy nhà máy điện hạt nhân, đồng thời phức tạp chia thành nhiều chế độ khác (hình 2) Hiện tượng sôi kênh dẫn bắt đầu với chế độ sơi bọt bão hòa, nước làm mát tiếp xúc với bó nhiên liệu nhận nhiệt bay tạo bọt nhỏ bám bề mặt ống (hình 2) Tiếp tục gia nhiệt, chế độ sơi bọt bão hòa thiết lập lúc bọt sinh nhiều Nếu bọt sinh nhiều đến mức tạo lớp màng khí bọc lên nhiên liệu, thông lượng nhiệt đạt đến tới hạn (Critical Heat Flux) đánh dấu chuyển tiếp chế độ sôi từ sôi bọt sang sôi màng Lúc màng khí có hệ số trao đổi nhiệt nhỏ nhiều so với nước, nhiên liệu tải lượng nhiệt lớn mà làm nóng nó, gây nứt vỡ nóng chảy Ngồi ra, lò nước áp lực PWR (Pressure Water Reactor), chế độ sơi bão hòa xảy phía vùng hoạt ảnh hưởng đến phân bố cơng suất Còn lò nước sơi BWR ( Boiling Water Reactor ), phân bố không đồng di chuyển bọt ảnh hượng đến độ phản ứng Vì lý quan trọng trên, em lựa chọn đề tài đồ án tốt nghiệp là: “Nghiên cứu mơ dòng hai pha kênh dẫn đứng sử dụng chương trình tính tốn mã nguồn mở OpenFOAM” Hình 1: Các chế độ dòng chảy truyền nhiệt kênh dẫn đứng hình ống [3] Lưới chia thành thể tích điều khiển với 10 đoạn theo trục x (theo bán kính), 150 đoạn theo trục y không chia nhỏ trục z Điều kiện biên điều kiện nguồn đặt bảng bảng Đoạn code cho phần lưới hóa trình bày phần phụ lục cuối luận văn Bảng 3: Điều kiên nguồn biên cho đại lượng Đại lượng Giá trị nguồn Giá trị vào Giá trị Giá trị tường Ub 0.517 0.517 0.517 Ua 1.112 1.112 1.112  0.05 zeroGradient zeroGradient p 20000 zeroGradient 20000 zeroGradient (zeroGradient dạng điều kiện biên mà giá trị gradient biên đại lượng [2]) Bảng 4: Thơng số vận chuyển (transport) nước khí đặt cho tốn mơ Đại lượng Giá trị b 1027kg / m3 C pb 4195 j / kg / K b Prb 1004e  06Pa.s C pa 1007 j / kg / K a Pra 1.84e  05Pa.s 7.02 0.7 42 3.3 Sơ đồ thuật giải toán động học dòng hai pha kênh dẫn đứng chương trình OpenFOAM Sau chia miền khơng gian, lưới hóa đặt điều kiền biên điều kiện đầu toán, ta lựa chọn solver OpenFOAM để giải phương trình Solver lựa chọn CompressibleTwoPhaseEulerFoam, tích hợp sẵn OpenFOAM Tuy nhiên để phù hợp với tốn mơ phỏng, đặc biệt khảo sát ảnh hưởng lực lên trình di chuyển bọt, solver bổ sung thêm mô hình loại lực 3.3.1 Vòng lặp hiệu chỉnh tăng tiến theo thời gian Các bước giải chương trình trình bày giống luận án thạc sỹ Ghione [3] Thuật giải trình lặp tăng tiến bước thời gian qua thể tích điều khiển chia Trước bước vào vòng lặp, số Courant Co  U t tính với t bước thời gian rời rạc, L kích thước L lưới, U vận tốc chênh lệnh hai pha Số Courant phải nhỏ để tránh bất ổn định tính tốn số có chuyển tiếp, khơng nhỏ q để giải xác cho thể tích điều khiển Thuật giải bao gồm hai vòng lặp hiệu chỉnh cho bước thời gian chuyển tiếp q trình tính tốn Vòng lặp ngồi outer corrector loop nhằm mục đích tính tốn số hạng quan trọng cần cho việc giải hai phương trình bảo tồn khối lượng bảo tồn momen, đồng thời giải phương trình momen chưa thỏa mãn phương trình liên tục Các phương trình solver gồm có:  “alphaEqn.H” chứa phương trình liên tục cho pha; 43  “kEpsilon.H” chứa mơ hình k-ε tính tốn cho tensor căng Reynolds phương trình momen;  “liftDragCoeffs.H” tính tốn loại lực tác dụng lên bọt;  “Ueqn.H” chứa phương trình momen cho pha rời rạc;  “pEqn.H” chứa phương trình áp suất để hiệu chỉnh trường vận tốc; Hình 25: Sơ đồ vòng lặp giải cho bước thời gian 3.3.2 Vòng lặp hiệu chỉnh áp suất-vận tốc Lý solver cần thêm vòng lặp PISO (Pressure-Implicit with Splitting of Operators) mối quan hệ hai đại lượng áp suất vận tốc (pressurevelocity coupling) Hai đại lượng phụ thuộc lẫn nhau, dẫn đến phương trình momen (1.2) khơng tuyến tính Vòng lặp PISO nhắm mục đích hiệu chỉnh qua lại hai đại lượng áp suất vận tốc để đạt trường vận tốc thỏa mãn 44 phương trình liên tục (1.1) phương trình momen (1.2) Thuật giải vòng lặp PISO trình bày giống [3] tổng quát sau: Đầu tiên trường vận tốc đoán cho thơng lượng đốn đưa vào vòng lặp PISO tương tự sơ đồ Jacobi: (2.9) Với A hệ phương trình tuyến tính tổng qt từ rời rạc phương trình momen pha mà khơng có thành phần gradient áp suất gradient trọng trường Sau solver, thông lượng hiệu chỉnh thêm số hạng ảnh hưởng trọng trường: k**  k*  ( g.S )( AD1,k ) f (2.10) Do trường vận tốc cần thỏa mãn phương trình liên tục, đồng thời đảm bảo tính bị chặn số pha     nên sử dụng định nghĩa thông lượng hỗn hợp: (2.11) Tương tự cho tốc hỗn hợp: (2.12) (1.12) ) Bây giờ, đại lượng áp suất giải qua phương trình riêngx ây dựng từ phương trình liên tục phương trình động lượng sử dụng thông lượng hỗn hợp 45 Theo Weller , phương trình liên tục (1.2) sử dụng dạng thông lượng hỗn hợp: .    a D (  a )  b D ( b )  a Dt b Dt (2.13) Khi đó, hiệu chỉnh giá trị thông lượng xét thêm gradient áp suất:   A  k D ,k k  k**      S f  f p f (2.14) Thay (2.11) (2.14) vào (2.13):  .      af   **  AD1,a      bf  a  f  AD1,b     b  f    D (  a )  b D ( b )  p    a  (2.15)   a Dt b Dt   Với Thay giá trị mật độ biểu diễn qua áp suất hệ số nén k   k p kí hiệu cho hệ số nén (trong mơ phỏng, chất lỏng khơng nén (2.16) ) Từ đó, sau biến đổi ta có phương trình áp suất thỏa mãn (1.1) (1.2):          gf      AD1,k   g      lf  f  AD1,k   l     f     p    .k**  g   g      l    l p    l prevIter lf   l .l prevIter  l  t      a p       prevIter    . prevIter    g gf  g g  t    (2.17)  46 (prevIter giá trị xác định trước vòng lặp) Giải phương trình áp suất, ta hiệu chỉnh lại phương trình vận tốc bao gồm ảnh hưởng trọng trường gradient áp suất Phương trình liên tục (1.1) với giá trị mật độ cần tính lại với trường vận tốc mới:   U k  U k*  reconstruct   g.S   AD1,k     A   k D ,k      S f  f p   f    k k     k kU k    k t (2.18) (2.19) Tổng kết, vòng lặp PISO trình bày tóm tắt sau: 1) Dự đốn trường thơng lượng mặt sử dụng (2.9) hiệu chỉnh lại (2.10); 2) Xây dựng giải phương trình áp suất (2.17); 3) Hiệu chỉnh thông lượng qua (2.14) vận tốc qua (2.18) nghiệm áp suất vừa giải được; 4) Cập nhật lại giá trị mật độ pha giải (2.19); Kết mô nhận xét Kết mơ thể qua việc dự đốn hệ số bọt (void fraction) vận tốc riêng pha, cụ thể định nghĩa sau:  Trong thể tích điều khiển (control volume), void fraction định nghĩa tỷ lệ thể tích pha khí thể tích tồn bộ: 47  Vg V  1 Vl V (3.1)  Vận tốc riêng pha định nghĩa sau: j g   vg jl  (1   )vl (3.2) Do mục đích nghiên cứu khảo sát thay đổi khoảng giới hạn mơ hình lực wall lubrication lift tác dụng lên bọt, code mô sử dụng hai mơ hình hệ số lực wall lubrication với ba trường hợp hệ số lực khác nhau:  Mơ hình hệ số lực tường Antal thay đổi với hai trường hợp: Antal 01 ( CW  0.011 CW  0.012 ứng với Cl  0.1 ) Antal 02 ( CW  0.01 CW  0.05 ứng với Cl  0.27 );  Mơ hình hệ số lực tường Frank khảo sát với CWC  10 , CWD  6.8 , p  1.7 ứng với Cl  0.27 ); Các hệ số wall hệ số lift cần phải lựa chọn hợp lý mơ hình mơ xác mặt vật lý Đối với mơ hình lực drag, code mơ sử dụng mơ hình Schiller-Naumann với lực turbulent dispersion, em sử dụng mơ hình Drew Lahey (RPI) (được nêu trang 17 25 luận văn) Kết mô từ OpenFOAM (hình 26) cho ba trường hợp thấy xuất peak void fraction Sự xuất peak cân lực wall lubrication lift lực turbulent dispersion đóng vai trò làm mịn phân bố bọt theo bán kính Tuy nhiên vị trí độ lớn peak thể mơ hình có khác biệt: 48  Với z/D = 6.00, hai trường hợp sử dụng mơ hình Antal ứng với hai hệ số Cl riêng thể profile void fraction giống với thực nghiệm, vị trí peak gần tường so với peak thực nghiệm Điều lý giải mơ hình lực tường Antal xây dựng cho dòng chảy tầng, khoảng cách ym từ bọt tới tường nhỏ nên lực tường mô hình Antal yếu so với hai hệ số Cl  0.27 Cl  0.1 Còn với trường hợp sử dụng mơ hình Frank, thấy rõ bọt xuất mật độ lớn gần trục ống Nguyên nhân dự đoán khuếch tán dòng chảy rối Sau kết mô phân bố bọt (void fraction) theo bán kính ống hai vị trí z/D = 6.00 z/D =53.5: Hình 26: Phân bố bọt dọc theo đường kính ống vị trí z/D = 6.00 53.5  Với z/D = 53.5, không cân lực mơ hình Antal thể rõ peak bị đẩy tường tạo dao động vật lý phân bố bọt gần tường với vai trò tốn tử gradient mơ hình turbulent dispersion Còn với trường hợp sử dụng mơ hình Frank, lúc dòng 49 chảy ổn định so với vị trị z/D = 6.00 nên khơng thấy bọt tập trung nhiều gần trục ống Lực tường thể mạnh lực lift nên peak bị đẩy xa peak thực nghiệm Ngồi ra, việc mơ hình chưa đầy đủ ảnh hưởng đến việc dự đốn vận riêng pha (hình 27 28) Mơ hình Antal dự đốn xu hướng vận tốc riêng thay đổi theo bán kính so với mơ hình Frank Hình 27: Sự thay đổi vận tốc riêng pha lỏng theo bán kính vị trí vị trí z/D = 6.00 53.5 50 Hình 28: Sự thay đổi vận tốc riêng pha khí theo bán kính vị trí vị trí z/D = 6.00 53.5 Bảng 5: Kết so sánh giá trị trung bình void fraction theo diện tích thể tích điều khiển dọc theo bán kính thực nghiệm mơ Alpha dự đốn 0.05 Antal01 với Antal02 với Frank với Cl = 0.1 Cl = 0.27 Cl = 0.25 z/D Thực nghiệm 0.053 0.075 0.077 0.083 53.5 0.086 0.081 0.082 0.07 Từ bảng ta thấy việc dự đoán giá trị void fraction vị trí z/D = 53.5 cho kết tốt so với vị trí z/D = Tuy nhiên, để đánh giá sai số tồn q trình mơ tốn, ngồi việc đánh giá sai số mơ hình lực cần phải đánh giá sai số xuất phát từ trình giải phương trình đưa nghiệm cho biến, cụ thể đánh giá độ xác phương 51 pháp tính tốn số sử dụng mơ Sai số đánh giá phương pháp tính tốn số sử dụng mơ gồm có: sai số trình rời rạc phương trình sai số trình giải ma trận đại số Kết thay đổi phương pháp nội suy rời rạc cho thấy sai số phương pháp tốn tối ưu nhất, sai số chủ yếu q trình mơ đến từ việc mơ hình mơ tả tượng động học dòng hai pha chưa đầy đủ Đó nhược điểm mơ hình hai dòng chảy lĩnh vực nghiên cứu dòng hai pha Kết luận Luận văn trình bày sở lý thuyết động học phương pháp chung mơ dòng hai pha kênh dẫn đứng (trường hợp chưa có trao đổi nhiệt chuyển trạng thái pha), phương trình bảo tồn theo mơ hình hai dòng chảy mơ hình tương tác lực qua mặt biên bọt Hướng nghiên cứu sử dụng chương trình tính tốn động học dòng chảy OpenFOAM mơ tốn dòng hai pha hướng hoàn toàn phù hợp Việt Nam, giúp người dùng nâng cao hiểu biết tượng thủy nhiệt dòng hai pha, đồng thời có khả xây dựng hệ thí nghiệm ảo sử dụng sở liệu thực nghiệm có sẵn từ nhóm nghiên cứu giới Kết nghiên cứu luận án mô cho trường hợp đơn giản, chưa tính tốn đến hợp tách bóng trao đổi nhiệt chuyển pha Đây mục nghiên cứu em thời gian tới 52 TÀI LIỆU THAM KHẢO [1] Mamoru Ishii and Takashi Hibiki Thermo-Fluid Dynamics of Two-Phase Flow, Second Edition [2] Henrik Rusche Computational Fluid Dynamics of Dispersed Two-Phase Flows at High Phase Fractions [3] MSThesis Ghione Development and validation of two-phase CFD model using OpenFOAM [4] MSThesis Michita OpenFOAM_NucleateBoiling_KTH_MS [5] Antal Analysis of phase distribution in fully developed laminar bubbly two-phase flow [6] Th Frank, J.-M Shi, E Krepper Non-drag Forces in Gas Liquid Bubbly Flows and Validation ofExisting Formulations [7] Jun-Mei Shi, Ulrich Rohde, Horst-Michael Prasser Turbulent Dispersion of Bubbles in Poly-dispersed Gas-Liquid Flows in a Vertical Pipe [8] R Maliska, E Paladino The Role of Virtual Mass, Lift and Wall Lubrication Forces in Accelerated Bubbly Flows [9] Nguyen VT, Song C-H, Bae BU, Euh DJ The Dependence of Wall Lubrication Force on Liquid Velociy in Turbulent Bubbly Two-phase Flow [10] John D Anderson, Jr Computational Fluid Dynamics – The Basics with Applications Takashi Hibiki,Mamoru Ishii, Zheng Xiao Axial interfacial area transport of vertical bubbly flows http://openfoam.org The Open Source CFD Toolbox – User Guide 4.0 and Programmers Guide 53 [13] T.R.Auton The lift force on a spherical rotational flow PHỤ LỤC Code OpenFOAM cho việc chia lưới miền khơng gian tính tốn convertToMeters 1; vertices ( (0 0) (0.025348 -0.001107) (0 3.061 0) (0.025348 3.061 -0.001107) (0.025348 0.001107) (0.025348 3.061 0.001107) ); blocks ( hex (0 2) (10 150 1) simpleGrading (1 1) //edgeGrading (1 1 1 1 1 1 1) ); edges ( ); boundary ( axis { type empty; faces 54 ( (0 2 0) ); } inlet { type patch; faces ( (0 0) ); } outlet { type patch; faces ( (2 2) ); } walls { type wall; faces ( (1 4) ); } front { type wedge; faces ( (0 4) ); } back { 55 type wedge; faces ( (0 1) ); } ); mergePatchPairs ( ); // ****************************************************** ******************* // 56 ... đoan luận văn tốt nghiệp tơi tự nghiên cứu, thiết kế thu hoạch kết hướng dẫn TS Nguyễn Văn Thái Để hoàn thành đồ án này, sử dụng tài liệu ghi phần tài liệu tham khảo, ngồi tơi không sử dụng tài... liệu khác mà không ghi Nếu sai, tơi xin chịu hình thức kỷ luật theo quy định Sinh viên thực Nguyễn Ngọc Quý MỤC LỤC DANH MỤC HÌNH VẼ DANH MỤC KÝ HIỆU DANH MỤC BẢNG BIỂU ... nước sơi BWR ( Boiling Water Reactor ), phân bố không đồng di chuyển bọt ảnh hượng đến độ phản ứng Vì lý quan trọng trên, em lựa chọn đề tài đồ án tốt nghiệp là: “Nghiên cứu mơ dòng hai pha kênh