Hướng dẫn chạy Module thuỷ lực (HD) trong MIKE_21 Flow Model thuộc bộ phần mềm MIKE

Hướng dẫn chạy Module thuỷ lực (HD) trong MIKE_21 Flow Model thuộc bộ phần mềm MIKE MỤC LỤC 1. GIỚI THIỆU VỀ MÔ HÌNH MIKE 21 FM ................................................................................................................2 2. CÁC DỮ LIỆU ĐẦU VÀO CỦA MIKE21 FM: ..........................................................................................................3 2.1 THIẾT LẬP LƯỚI TÍNH TOÁN. ......................................................................................................................................3 2.1.1 Tạo file .mdf từ file số liệu thô .xyz ............................................................................................................3 2.1.2 Chỉnh sửa lại đường bờ cho vùng cần tạo lưới tính...................................................................................8 2.1.3 Thiết lập lưới tính tam giác cho vùng tính. ................................................................................................9 2.2 SỬ DỤNG CÔNG CỤ TIME SERIES TRONG MIKE ZERO ĐỂ NHẬP MỰC NƯỚC THỰC ĐO TỪ DẠNG FILE.TXT:.................................13 2.3 THIẾT LẬP ĐIỀU KIỆN BIÊN: ......................................................................................................................................18 3. THIẾT LẬP MÔ PHỎNG TRONG MIKE21 FM:.....................................................................................................21 3.1 KHỞI ĐỘNG MIKE21 FM. ......................................................................................................................................21 3.2 VÙNG LÀM VIỆC (DOMAIN):....................................................................................................................................22 3.3 THỜI GIAN MÔ PHỎNG (TIME):................................................................................................................................23 3.4 ĐIỀU KIỆN BAN ĐẦU (INITIAL CONDITIONS):................................................................................................................23 3.5 ĐIỀU KIỆN BIÊN (BOUNDARY CONDITIONS):................................................................................................................24 3.6 KẾT QUẢ ĐẦU RA (OUTPUTS):..................................................................................................................................25 4. HIỆU CHỈNH CÁC THÔNG SỐ CỦA MÔ HÌNH:....................................................................................................28 Hướng dẫn sử dụng MIKE 21FM (HD) 2 1. Giới thiệu về mô hình MIKE 21 FM - Mike là bộ phần mềm của viện Thuỷ Lợi Đan Mạch được ứng dụng trong việc tính toán về thuỷ lực, tài nguyên và môi trường nước, bao gồm cả trong sông, vùng cửa sông, ven biển và biển Bộ phần mềm này đã được ứng dụng khá hiệu quả trong thực tế tại nhiều quốc gia trên thế giới. - Trong bộ phần mềm Mike có chia ra làm nhiều mô hình, như Mike Zero, Mike_11, Mike_21, Mike_3, Mike_She - Trong bài này ta quan tâm chủ yếu tới hai phần chính, đó là Mike Zero và Mike_21. 1) Mike Zero: Đây là 1 phần quan trọng trong mô hình Mike, tất cả các thông số đầu vào cũng như các phương án mô phỏng, các điều kiện biên đều được thiết lập thônng qua Mike Zero. 2) Mike_21: Là mô hình dòng chảy mặt 2D, mô hình Mike_21 được ứng dụng để mô phỏng các quá trình thuỷ lực và các hiện tượng về môi trường trong các hồ, các vùng cửa sông, vùng vịnh, vùng ven bờ và các vùng biển. - Mô hình Mike_21 bao gồm một số các Module sau: + Hydrodynamic (HD) + Transport (TR) + ECO Lab (EL) Hướng dẫn sử dụng MIKE 21FM (HD) 3 + Mud Transport (MT) + Sand Transport (ST) Trong các module của mô hình Mike_21 thì Module thuỷ động lực (Hydrodynamic HD) là module cơ bản nhất, nó cung cấp chế độ thuỷ động lực làm nền tảng cho quá trình tính toán của các module thuỷ lực khác. - Trong nghiên cứu này ta sẽ bước đầu khai thác module HD của Mike-21 Flow Model FM, và sau đây sẽ là một số bước cơ bản để có thể chạy module thuỷ động lực (HD) - Bài toán mẫu được đặt ra là sử dụng mô hình Mike_21 để quan trắc sự dao động mực nước và vận tốc dòng chảy của vùng biển nằm giữa Đan Mạch và Thụy Điển. 2. Các dữ liệu đầu vào của Mike21 FM: 2.1 Thiết lập lưới tính toán. Cửa sổ làm việc sau khi kích hoạt chương trình Mike 21 2.1.1 Tạo file .mdf từ file số liệu thô .xyz Ta sẽ bắt đầu việc tạo lưới tính toán bằng công cụ Mesh Generator: - Sau khi nhấn chuột vào New, chương trình Mike sẽ hiện ra cửa sổ như Hình I.1.1 ở dưới. - Chọn Mike Zero - Trong Mike Zero, chọn công cụ Mesh Generator Hướng dẫn sử dụng MIKE 21FM (HD) 4 (New Mike Zero  Mesh Generator) Hình I.1.1: Công cụ thiết lập lưới trong Mike Zero Sau khi khởi động Mesh Generator cửa sổ Workspace projection sẽ hiện ra như trong Hình I.1.2 và ta phải xác định hệ quy chiếu của khu vực nghiên cứu. Như trong ví dụ này ta sẽ chọn hệ quy chiếu của vùng tính là UTM-33. Hình I.1.2: Cửa sổ chọn hệ quy chiếu cho khu vực tính toán Sau khi chọn Projection ta sẽ có lưới ô toạ độ của vùng tính toán như trong hình dưới. Hướng dẫn sử dụng MIKE 21FM (HD) 5 Hình I.1.3: Lưới ô toạ độ của vùng làm việc Tiếp theo, ta cần có file số liệu thô (dạng file.xyz), trong đó mỗi điểm thuộc đường bờ hay các vị trí độ sâu sẽ được mô tả bởi 3 giá trị: x, y, và z. Trong đó x và y là toạ độ của điểm được xác định theo hệ kinh vĩ độ (Long/Lat), còn z là ID(đối với file đường bờ) hoặc là độ sâu (đối với file độ sâu). Các file có dạng như trong hình sau: File đường bờ với thứ tự các cột lần lượtlà: kinh độ, vĩ độ và ID của đường bờ (được quy định là 1): Hình I.1.4:File Shoreline.xyz Hướng dẫn sử dụng MIKE 21FM (HD) 6 File độ sâu với thứ tự các cột là: kinh độ, vĩ độ và độ sâu của điểm: Hình I.1.5:File Water.xyz Trên thanh công cụ phía trên, nhấp chuột vào Data và chọn Import Boundary, cửa sổ làm việc sẽ hiện ra, chỉ đường dẫn tới thư mục chứa file.xyz, lựa chọn file đường bờ và Open. Hình I.1.6: Cửa sổ nhập số liệu đường bờ Cần chú ý khi nhập số liệu đường bờ, trong phần lựa chọn Projection ta phải chọn hệ Long/Lat Hướng dẫn sử dụng MIKE 21FM (HD) 7 Hình I.1.7: Cửa sổ Boundary Properties Hình I.1.8: Dữ liệu số hoá đường bờ sau khi nhập, được thể hiện trên vùng làm việc Hướng dẫn sử dụng MIKE 21FM (HD) 8 2.1.2 Chỉnh sửa lại đường bờ cho vùng cần tạo lưới tính. Số liệu đường bờ sau khi được nhập có thể sẽ có 1 số các đường nối như đã thấy trong hình trên, việc đầu tiên cần làm là sửa lại đường bờ cho đúng: xoá các nốt và các Vertices không thuộc đường bờ trước khi tạo lưới tính. Để chỉnh sửa ta cần dùng đến 1 số các công cụ của thanh công cụ phía trên như: + Select points: Được dùng để xác định thuộc tính của các nodes hoặc các vertices như toạ độ của điểm, chuyển đổi giữa Nodes => Vertices và ngược lại + Select acrs: Được sử dụng để xác định thuộc tính của 1 cung như: Phân biệt giữa đường biên lỏng và đường biên cứng, chỉnh sửa lại sự phân bố của các Vertices, xuất số liệu của đường biên ra dạng file.xyz + Select polygons: Được dùng để xác định thuộc tính của 1 vùng khép kín. + Insert nodes: Được dùng chủ yếu để add thêm các nốt cần thiết khi lập lưới tính. + Draw arcs: Được dùng để nối liền đường bờ khi đường bờ bị đứt quãng, dùng để giới hạn đường biên lỏng của vùng nghiên cứu hoặc dùng để khoanh vùng 1 vùng tính với bước lưới khác biệt so với những vùng còn lại. + Insert polygons: được dùng để đánh dấu những vùng không phải là biển trong khu vực nghiên cứu như các đảo hoặc dùng để đánh dấu 1 vùng trong khu vực nghiên cứu có các thuộc tính khác với những vùng còn lại. + Move points: Được dùng để chỉnh sửa lại vị trí các nodes. + Delete points: Được dùng để xoá các đường nối, các điểm thừa không cần thiết. + Ngoài ra ta có thể dùng thêm 1 số các công cụ khác như Zoom In, Zoom Out, Pan để hỗ trợ trong việc chỉnh sửa đường biên. Việc tiếp theo là ta phải xác định các biên lỏng phía Bắc và phía Nam bằng cách kẻ các đường nối giữa Đan Mạch và Thuỵ Điển. Các đường biên này phải có toạ độ gần với các trạm thực đo đã cho. Lựa chọn các biên lỏng vừa thiết lập, nhấp chuột phải chọn Properties và thiết lập thuộc tính cho các biên này. Chọn thuộc tính 2 cho biên phía bắc và 3 cho biên phía nam. Đường biên giữa đất và nước sẽ được Mesh Generator mặc định và có thuộc tính là 1. Hướng dẫn sử dụng MIKE 21FM (HD) 9 Hình I.1.9: Cửa sổ thuộc tính của đường biên. Sau khi đã sửa lại đường bờ, giới hạn các biên lỏng và đánh dấu các đảo trong vùng, tiếp theo ta có thể thiết lập lưới tính toán cho vùng đã đóng kín. Nhưng trước tiên ta nên làm trơn lại đường bờ, Cần chú ý rằng quá trình tạo lưới sẽ được bắt đầu từ đường biên vì vậy số mắt lưới tính toán sẽ phụ thuộc rất nhiều vào số nốt và số Vertices của đường bờ. Ta có thể sử dụng công cụ Select arcs và chọn Redistribute Vertices để chỉnh sửa lại khoảng cách giữa các Vertices cho đều nhau, và với những vùng cần đặc biệt quan tâm như bến cảng, phá ta có thể chia lại khoảng cách giữa các Vertices với khoảng cách ngắn hơn. Sau khi đã chỉnh sửa lại đường biên, ta có thể ghi lại dưới dạng file.mdf 2.1.3 Thiết lập lưới tính tam giác cho vùng tính. Bước tiếp theo ta sẽ tạo lưới cho vùng tính. Đầu tiên, ta cần đánh dấu các vùng không cần chia lưới (như các đảo) bằng cách add thêm các nốt màu xanh lá cây. Hướng dẫn sử dụng MIKE 21FM (HD) 10 Hình I.1.10: Đường bờ sau khi đã được chỉnh sửa và khép kín vùng nghiên cứu Để tạo lưới tính, ta chọn Mesh trong thanh công cụ phía trên và chọn Triangulate. (Mesh Triangulate) Hình I.1.11: Cửa sổ thiết lập lưới đo cho vùng nghiên cứu. Cửa sổ Triangulate sẽ hiện ra như hình dưới và ta phải điền vào đó 1 số các giá trị như + Diện tích lớn nhất có thể của các ô lưới (m2) + Góc nhỏ nhất cho phép của các ô lưới. Hướng dẫn sử dụng MIKE 21FM (HD) 11 + Số mắt lưới tối đa của vùng tính. Cần chú ý rằng với những vùng đặc biệt cần quan tâm, ta có thể chia lưới mịn hơn(các ô lưới có diện tích nhỏ hơn) bằng cách khoanh vùng quan tâm, đánh dấu bằng công cụ Insert Polygons,sau đó chọn Select Polygons, phấp chuột phải chọn Properties và chọn Use local maximum area. Tiếp theo, ta có thể làm mịn lưới bằng cách chọn Mesh ở thanh công cụ phía trên và chọn Smooth Mesh (Mesh Smooth Mesh) Hình I.1.12: Vùng tính toán sau khi đã thiết lập lưới Bước tiếp theo ta sẽ nhập số liệu độ sâu: Nhấn chuột vào Data trên thanh công cụ và chọn Import Scatter Data, cửa sổ nhập số liệu độ sâu sẽ hiện ra như trong hình dưới: Hình I.1.13: Cửa sổ nhập file độ sâu địa hình vùng nghiên cứu Hướng dẫn sử dụng MIKE 21FM (HD) 12 Nhấn vào Add button và chỉ đường dẫn tới folder chứa file độ sâu (file.xyz). Khi add file độ sâu ta cũng cần nhớ chọn hệ quy chiếu Long/Lat khi được hỏi đến Map projection. (Data -> Import Scatter Data ->Add -> water.xyz) Hình I.1.14: Vùng tính toán khi đã nhập số liệu độ sâu mực nước Sau đó để nội suy độ sâu từ các điểm phân tán của file.xyz vào lưới tính ta chọn Mesh trên thanh công cụ phía trên và chọn Interpolate. (Mesh Interpolate) Hình I.1.15: Dạng file.mdf được tạo ra sau khi tiến hành nội suy Hướng dẫn sử dụng MIKE 21FM (HD) 13 Sau khi tiến hành nội suy ta đã có thể export ra dạng file.mesh để sử dụng trong quá trình mô phỏng: Nhấn vào Mesh trên thanh công cụ phía trên, lựa chọn Export Mesh, đặt tên file và Save. (Mesh Export Mesh) Hình I.1.16: Dạng file.mesh đã được tạo ra. 2.2 Sử dụng công cụ Time series trong Mike Zero để nhập mực nước thực đo từ dạng file.txt: Số liệu về mực nước thực đo của các trạm đo trong Mike_21 đựợc quản lý dưới dạng file.dfs0. Để tạo ra các file.dfs0 này ta cần có đầu vào là các file.txt có định dạng như trong hình dưới: Hình I.2.1: Dạng file.txt ghi mực nước thực đo tại các trạm đo Hướng dẫn sử dụng MIKE 21FM (HD) 14 Để tạo ra 1 file.dfs0 ta làm như sau: Trở về màn hình chính của Mike, nhấn vào New trên thanh công cụ phía trên, trên màn hình sẽ xuất hiện cửa sổ như hình dưới: Hình I.2.2:Công cụ Time Series trong Mike Zero Chọn Mike_Zero, chọn Time Series và Ok (FileNew TimeSeries) Trong cửa sổ New TimeSeries Chọn From ASCII file, trong cửa sổ làm việc mới hiện ra ta chọn Time description là Equidistant Calerdar Axis, chỉ đường dẫn tới chỗ đặt file mực nước dưới dạng file.txt và Ok. Nhấp chuột phải vào file dữ liệu mới được tạo ra, chọn Properties, chọn kiểu dữ liệu là “Water Level”, đặt tên file và save lại ta sẽ có được dạng file.dfs0. Làm tương tự cho các trạm đo mực nước còn lại. Hướng dẫn sử dụng MIKE 21FM (HD) 15 Hình I.2.3: Cửa sổ New TimeSeries Hình I.2.4: Cửa sổ nhập mực nước từ file.txt Hình I.2.5: Time seriesProperties Hướng dẫn sử dụng MIKE 21FM (HD) 16 Hình I.2.6: Dạng file mực nước (file.dfs0) sau khi được thiết lập Và để theo dõi quá trình thay đổi của mực nước theo thời gian của 1 trạm nào đó một cách trực quan hơn, ta có thể sử dụng công cụ Plot Composer trong Mike Zero để vẽ sự thay đổi của mực nước dưới dạng biểu đồ: Chọn New ở thanh công cụ phía trên, nhấn vào Mike Zero, chọn Plot Composer và nhấn Ok. New => Mike Zero => Plot Composer Hình I.2.7: Công cụ Plot Composer trong Mike Zero Hướng dẫn sử dụng MIKE 21FM (HD) 17 Cửa sổ làm việc sẽ hiện ra, nhấn vào Plot ở thanh công cụ phía trên và chọn Insert New Plot Object Sau đó, trong cửa sổ làm việc Insert Plot Object, ta chọn dạng biểu đồ là Time Series Plot Hình I.2.8: Cửa sổ Insert Plot Object ` Trong cửa sổ làm việc Time Series Plot Properties nhấn vào New Item và chỉ đường dẫn tới chỗ để file.dfs0 vừa mới tạo ra ở bước trên. Hình I.2.9: Cửa sổ Time Series Plot Properties Trong 1 file.dfs0 có thể chứa nhiều yếu tố khác nhau như mực nước, vận tốc ta sẽ tick vào yếu tố nào cần vẽ ra biểu đồ rồi nhấn Ok.. Hướng dẫn sử dụng MIKE 21FM (HD) 18 Hình I.2.10: Cửa sổ nhập file.dfs0 cần vẽ biểu đồ Hình I.2.11: Biểu đồ đường quá trình mực nước của 1 trạm đo 2.3 Thiết lập điều kiện biên: Để tạo ra file điều kiện biên dạng line, ta sử dụng công cụ Profile Series trong Mike Zero: Nhấn vào New, trong Mike_Zero ta chọn Profile Series. Hướng dẫn sử dụng MIKE 21FM (HD) 19 Hình I.3.1: Công cụ Profile Series trong Mike Zero Cửa sổ New File sẽ xuất hiện và ta chọn Blank T1 Document => Ok Hình I.3.2: Cửa sổ lựa chọn dạng file đầu vào cho file điều kiện. Hướng dẫn sử dụng MIKE 21FM (HD) 20 Hình I.3.3: Cửa sổ thiết lập các thuộc tính của biên Ta cần điền đầy đủ các thông tin sau cho từng biên: • Axis Type: ta sẽ chọn Equidistant Calendar Axis • Start Time: Ngày bắt đầu sẽ được xác định từ tài liệu thực tế. • Time step: khoảng thời gian giữa mỗi lần đo trong tài liệu thực đo mực nước của các trạm(ví dụ: mỗi lần đo cách nhau 1tiếng thì Time Step sẽ là 3600s). • No. of time steps: Số bước thời gian trong toàn bộ chuỗi tài liệu của các trạm đo. • No. of grid points: Số các điểm quan trắc trên đường biên, ví dụ trong bài này tại mỗi biên ta có hai trạm đo ở hai đầu thì chọn No. of grid points là 2 • Grid Step: Độ dài của biên lỏng (khoảng cách giữa 2 trạm đo ở 2 đầu mút của biên) tính bằng mét (m). Sau khi điền đầy đủ các thông tin trên và nhấn Ok, một cửa sổ mới sẽ hiện ra, ta sẽ coppy số liệu mực nước của hai trạm ở hai đầu mút của biên trên và paste vào hai cột số 0 và số 1 như trong hình dưới sau đó save lại. Mực nước tại các biên sẽ được nhập theo thứ tự ngược chiều kim đồng hồ. Như vậy ta đã tạo ra 1 file.dfs1 là file điều kiện biên dạng line. Làm tương tự như các bước trên cho biên còn lại. Hướng dẫn sử dụng MIKE 21FM (HD) 21 Hình I.3.4: Dạng fìle điều kiện biên dạng line (file.dfs1) sau khi được thiết lập 3. Thiết lập mô phỏng trong Mike21 FM: 3.1 Khởi động Mike21 FM. Chọn New ở thanh công cụ phía trên, khi cửa sổ New hiện ra ta chọn Mike_21 rồi chọn Flow Model FM => Ok Hình II.1.1: Module Flow Model FM trong Mike21. Hướng dẫn sử dụng MIKE 21FM (HD) 22 Hình II.1.2: Giao diện chính của Mike_21 Flow Model FM 3.2 Vùng làm việc (Domain): - Trong tab Mesh and Bathymetry ta sẽ phải nhập file.mesh đã tạo được ở bước trên: nhấp chuột vào nút [] và chỉ đường dẫn tới chỗ đặt file.mesh đã tạo ra ở trên. Sau khi nhập file.mesh vùng nghiên cứu sẽ hiện ra ngay ở cửa sổ phía dưới như trong hình II.2 - Tab Domain Specification trong Domain sẽ cho ta biết 1 số thông tin về vùng làm việc như: Map projection, Minimum depth cutoff - Tab Boundary names cho phép ta chỉnh sửa lại tên của các biên. Hình II.2: Cửa sổ Domain Hướng dẫn sử dụng MIKE 21FM (HD) 23 3.3 Thời gian mô phỏng (Time): Là khoảng thời gian ta muốn thiết lập sự mô phỏng cho vùng nghiên cứu. Thời gian mô phỏng sẽ được xác định từ chuỗi số liệu thực đo, nó phải nằm trong khoảng thời gian mà các trạm đo mực nước có số liệu đầy đủ. - Simulation start date: Ngày bắt đầu quá trình mô phỏng, được xác định dựa vào chuỗi tài liệu thực đo đã có và phụ thuộc vào yêu cầu nghiên cứu của người chạy mô hình. - Time step interval: Là khoảng thời gian 1 bước tính của mô hình, đơn vị được tính là giây (s). - No. of time steps: là số bước tính toán trong khoảng thời gian mô phỏng. Ta có thể tính ra bằng cách lấy khoảng thời gian mô phỏng chia cho thời gian của 1 bước tính => số bước tính toán. - Simulation end date: Từ Simulation start date, Time step interval và No. of time steps mô hình sẽ tự xác định ngày kết thúc quá trình mô phỏng. Hình II.3: Cửa sổ thiết lập thời gian cho quá trình mô phỏng 3.4 Điều kiện ban đầu (Initial Conditions): Điều kiện ban đầu bao gồm Surface elevation, U-Velocity và V-Velocity Trong trường hợp này ta chọn Surface elevation là -0,37m, các giá trị U và V Velocity là 0. Hướng dẫn sử dụng MIKE 21FM (HD) 24 Hình II.4: Cửa sổ điều kiện ban đầu 3.5 Điều kiện biên (Boundary Conditions): Trong ví dụ này ta sẽ có 2 biên: North và South - Type: Trong trường hợp này ta sẽ lựa chọn kiểu của biên là Specified level. “Specified level” có nghĩa là mực nước sẽ là nhân tố chính của biên ảnh hưởng tới quá trình mô phỏng. - Boundary Format phải được để ở dạng Varying in time and along boundary đối với kiểu file.dfs1 đã tạo ở trên. Sau đó nhấn vào nút Select và chỉ đường dẫn tới nơi đặt file điều kiện biên (file.dfs1) đã tạo được ở các bước trên. - Solf start là khoảng thời gian khởi động nhẹ mô hình khi điều kiện biên chưa ổn định. Ta có thể chọn Time interval là 2 tiếng (7200s) và Reference value bằng giá trị Initial Surface là -0,37(m). Hướng dẫn sử dụng MIKE 21FM (HD) 25 Hình II.5: Cửa sổ nhập điều kiện biên 3.6 Kết quả đầu ra (Outputs): Nhánh Outputs là nhánh kết quả đầu ra. Ta có thể có 1 hoặc nhiều kết quả đầu ra, bằng cách chọn nhánh Outputs và nhấn vào nút New Outputs. Hình II.6.1: Cửa sổ Outputs Tại mỗi cửa sổ đầu ra ta có thể lựa chọn 1 số các yếu tố sau: - Trong Tab Outputs specification có các lựa chọn sau: Hướng dẫn sử dụng MIKE 21FM (HD) 26 + Field type: ta sẽ chọn là 2D (horizontal) + Outputs format: Ta có thể lựa chọn kết quả đầu ra tại 1 điểm (Point series) , trên 1 đường (Line Series), hoặc mô phỏng cho toàn vùng nghiên cứu (Area series). + First time step: thường chọn là 0 + Last time step: là bước thời gian cuối cùng của khoảng thời gian mô phỏng + Frequency: là tần suất cho ta 1 giá trị kết quả đầu ra. Ví dụ bước thời gian đã chọn ở trên là 8s, tần suất là 15 thì sau mỗi 120s (2 phút) ta sẽ nhận được 1 giá trị mô phỏng của mô hình. + Nếu kết quả đầu ra tại 1 điểm thì ta cần xác định rõ vị trí, toạ độ của điểm quan trắc trên vùng nghiên cứu. Hình II.6.2: Cửa sổ lựa chọn các thuộc tính của file kết quả - Tab Output items cho phép ta lựa chọn các yếu tố đầu ra mà ta quan tâm bằng cách tick vào các ô tương ứng. Hướng dẫn sử dụng MIKE 21FM (HD) 27 Hình II.6.3: Cửa sổ lựa chọn các yếu tố kết quả đầu ra. Trong ví dụ này, các điều kiện khác ta có thể để ở chế độ mặc định: - Flood and dry : Drying depth: 0,005 Flooding depth: 0,05 Wetting depth: 0,1 - Density type: Barotropic - Eddy Viscosity: Được mô phỏng theo công thức Smagorinsky với hệ số là 0,28 - Bed resistance: Constan value: 32m1/3/s - Coriolis Forcing - Wind Forcing - Ice Coverage is not included - Tidal Potential is not included - Precipitation-Evaporation is not included - Wave Radiation is not included - Decoupling is not included Hướng dẫn sử dụng MIKE 21FM (HD) 28 Bây giờ mọi thứ đã sẵn sàng để chạy Mike21 Flow Model FM: Nhấn vào Run ở thanh công cụ phía trên và chọn Start simulation. (Run Start simulation) 4. Hiệu chỉnh các thông số của mô hình: Để hiệu chỉnh mô hình, ta cần có số liệu thực đo của 1 trạm đo nào đó nằm trong vùng nghiên cứu. Trong ví dụ mẫu này ta sẽ sử dụng giá trị mực nước thực đo của trạm Nrd_Roese để so sánh với kết quả đầu ra và hiệu chỉnh mô hình. Kết quả đầu ra của mô hình đối với Output dạng Point series sẽ cho ta một bảng giá trị về các yếu tố mà ta quan tâm trong khoảng thời gian mô phỏng dưới dạng file.dfs0. Hình III.1: File.dfs0 _Kết quả của Output dạng Point series Kết quả đầu ra của mô hình đối với Output dạng Area series được ghi lại dưới dạng file.dfsu. Nó cho phép ta quan sát một cách trực quan hơn về quá trình thay đổi của mực nước cũng như vận tốc dòng chảy trên toàn vùng nghiên cứu. Và ta cũng có thể lấy ra giá trị dao động mực nước hay vận tốc dòng chảy tại bất cứ thời điểm nào và tại bất kỳ vị trí nào trong vùng nghiên cứu. Kết quả lấy ra này được ghi lại dưới dạng file.dfs0 Hướng dẫn sử dụng MIKE 21FM (HD) 29 Hình III.2: File.dfs1 _Kết quả của Output dạng Area series Sau khi có kết quả đầu ra, ta có thể sử dụng công cụ Plot Composer để thể hiện các giá trị thực đo và kết quả tính toán của mô hình lên cùng 1 biểu đồ để so sánh một cách trực quan hơn. Khi nhận thấy kết quả đầu ra của mô hình và giá trị thực đo có sai số, ta cần hiệu chỉnh lại các thông số của mô hình. Hai thông số chính ảnh hưởng tới kết quả đầu ra của mô hình trong ví dụ này là hệ số Manning trong Bed Resistance và hệ số của công thức Smagrinsky trong Eddy viscosity. Ta sẽ hiệu chỉnh 2 giá trị này sao cho kết quả đầu ra của mô hình có sai số nhỏ nhất so với giá trị thực đo. Khi hiệu chỉnh, ta nên thay đổi từng thông số trong mỗi lần mô phỏng để đánh giá được mức độ ảnh hưởng của từng thông số tới kết quả đầu ra. Với số liệu của ví dụ của mẫu đã có, chương trình Mike_21 đã cho ta kết quả đầu ra khá sát với số liệu thực đo như đã thấy trong biểu đồ so sánh hình III.3 ở trên. Và trên cơ sở bước đầu chạy thành công Module thuỷ lực cho bài toán mẫu, ta sẽ tiến hành chạy và hiệu chỉnh mô hình cho vùng biển thuộc Vịnh Bắc Bộ để phục vụ cho các nghiên cứu khoa học tiếp theo. Hướng dẫn sử dụng MIKE 21FM (HD) 30 Hình III.3: Biểu đồ so sánh giữa giá trị thực đo tại trạm Nrd_Roese với kết quả đầu ra của mô hình dạng Point series