BÁO CÁO MÔ TẢ HỌC THUẬT VÀ TỔNG QUAN NỘI DUNG GIÁO TRÌNH: DỰ BÁO THỜI TIẾT BẰNG PHƯƠNG PHÁP SỐ TRỊ
Thông tin thư mục xuất bản
- Tên giáo trình: Dự báo thời tiết bằng phương pháp số trị
- Tác giả: Trần Tấn Tiến
- Nhà xuất bản: Nhà xuất bản Đại học Quốc gia Hà Nội (NXB ĐHQGHN - 1997)
- Chịu trách nhiệm xuất bản: Giám đốc Nguyễn Văn Thỏa, Tổng biên tập Nghiêm Đình Vỳ
- Người nhận xét: PTS. Nguyễn Đăng Quế
- Biên tập và sửa bản in: Đinh Văn Vang
- Trình bày bìa: Ngọc Anh
- Mã số xuất bản: 01.97 s; Số xuất bản: 264/CXEB; Số trích ngang: 262 KH/XB; In 500 cuốn, nộp lưu chiểu tháng 10/1997.
Tổng quan về giáo trình
Dự báo thời tiết bằng phương pháp số trị là giáo trình chuyên ngành thuộc lĩnh vực Khí tượng Động lực học, được biên soạn nhằm phục vụ chương trình đào tạo đại học và sau đại học chuyên ngành Khí tượng học, Thủy văn và Vật lý Khí quyển. Môn học dự báo số trị khảo sát định lượng các quy luật vật lý của các quá trình khí quyển thông qua việc nghiên cứu các phương pháp tích phân số hệ phương trình thủy nhiệt động lực học.
Mục tiêu học tập của giáo trình hướng tới việc trang bị cho người học cơ sở lý thuyết hoàn chỉnh về hệ phương trình cơ bản mô tả chuyển động của khí quyển, phương pháp phân tích thứ nguyên và đánh giá bậc độ lớn của các quá trình quy mô synop và quy mô hành tinh, cũng như kỹ thuật chuyển đổi hệ phương trình sang các hệ tọa độ tựa Lagorăng phục vụ tích phân số.
Cấu trúc giáo trình tiếp cận theo hướng toán lý chặt chẽ: xuất phát từ các định luật bảo toàn vật lý (bảo toàn động lượng, bảo toàn khối lượng, bảo toàn năng lượng), thiết lập hệ phương trình khép kín cho chất lỏng lý tưởng và chất lỏng rối, sau đó áp dụng phép xấp xỉ tựa tĩnh học và gần đúng địa chuyển để xây dựng các sơ đồ dự báo hình thế khí áp và trường khí tượng.
Điểm đặc trưng của tài liệu là sự tích hợp chi tiết lịch sử phát triển của ngành dự báo số trị qua ba giai đoạn: giai đoạn hình thành với các công trình của V. Bjerknes, L. F. Richardson (1922), A. A. Fridman (1914); giai đoạn phát triển sơ đồ dự báo tựa tĩnh của I. A. Kibel (1940), C.-G. Rossby (1938–1940), E. N. Blinova (1943), N. I. Buleev và G. I. Marchuk (1951); và giai đoạn ứng dụng máy tính điện tử với các mô hình của J. G. Charney, N. A. Phillips (1950) và S. L. Belousov (1954).
Nội dung kiến thức cốt lõi
Các chương và chủ đề chính
Nội dung trọng tâm được xây dựng trên hệ thống các chủ đề toán lý khí tượng, tiêu biểu qua cấu trúc khảo sát của Chương I:
-
Hệ phương trình nhiệt động lực học cho chất lỏng lý tưởng trong hệ tọa độ Đề-các $(x, y, z)$:
Xây dựng hệ 5 phương trình cơ bản mô tả 7 đại lượng khí tượng ($u, v, w, P, \rho, T, q$):
- Phương trình chuyển động Navier-Stokes có tính đến lực khối (trọng lực $g = 9{,}8\text{ m/s}^2$ và lực Coriolis $f = 2\omega \sin\varphi$, $l = 2\omega \cos\varphi$).
- Phương trình liên tục biểu diễn định luật bảo toàn khối lượng: $\frac{\partial \rho}{\partial t} + \text{div}(\rho \mathbf{V}) = 0$.
- Phương trình trạng thái khí lý tưởng: $P = \rho R T$ ($R = 287\text{ m}^2/\text{s}^2\cdot\text{độ}$).
- Phương trình nhập nhiệt (nguyên lý thứ nhất nhiệt động lực học) với các nguồn nhiệt bức xạ $\varepsilon_{bx}$, rối $\varepsilon_R$ và chuyển pha hơi nước $\varepsilon_{ph} = Lm$.
- Phương trình vận chuyển ẩm: $\frac{dq}{dt} = -m$.
-
Hệ phương trình thủy nhiệt động lực học cho khí quyển rối:
Sử dụng quy tắc lấy trung bình Reynolds để tách biến thành giá trị trung bình $\bar{f}$ và xung lượng rối $f'$. Mô hình hóa các ứng suất rối và dòng nhiệt/ẩm rối thông qua hệ số trao đổi rối thẳng đứng $K$, nằm ngang $K'$ và gradien của đại lượng trung bình.
-
Chuyển đổi sang các hệ tọa độ tựa Lagorăng:
- Hệ tọa độ áp suất $(x, y, p, t)$: Sử dụng phương trình tĩnh học $\frac{\partial P}{\partial z} = -\rho g$ để thay biến $z$ bằng $p$. Vận tốc thẳng đứng được thay bằng đại lượng tương tự vận tốc thẳng đứng $\tau = \frac{dp}{dt}$. Khử ẩn số mật độ $\rho$ khỏi phương trình chuyển động và liên tục.
- Hệ tọa độ địa hình $\sigma$ $(x, y, \sigma, t)$: Với $\sigma = P / P_x$ ($P_x$ là áp suất mặt đất). Giúp chuẩn hóa điều kiện biên dưới tại mặt đất thành $\sigma = 1$ ($\dot{\sigma} = 0$).
- Hệ tọa độ đẳng Entropi $\theta$ $(x, y, \theta, t)$: Thiết lập trên các mặt nhiệt độ thế vị không đổi.
-
Phương trình xoáy và phương trình Đivecgiăng:
Khử các thành phần đạo hàm cục bộ có sai số lớn trong phương trình chuyển động bằng cách vi phân chéo, dẫn xuất phương trình biến đổi xoáy tuyệt đối $\Omega_a = \Omega + f$ ($\frac{d\Omega_a}{dt} = -(\Omega + f)D - \beta v$) và phương trình đivecgiăng $D = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y}$.
-
Phương pháp đánh giá bậc đại lượng và đơn giản hóa phương trình:
Trình bày 3 phương pháp xác định độ lớn: phương pháp thống kê nhiều năm của Fridman và Hesselberg; phương pháp giá trị đặc trưng của M. I. Yudin; phương pháp giải tích quy mô (scale analysis) sử dụng số Kibel $K = \frac{U}{f L} \approx 0{,}104$ ở vĩ độ trung bình. Trên cơ sở đó, thiết lập phương trình xoáy tựa địa chuyển và phương trình cân bằng phi tuyến dạng Monge-Ampère đối với hàm dòng $\psi$.
-
Hệ phương trình thủy nhiệt động lực học trong hệ tọa độ cầu $(\lambda, \theta, r)$:
Ứng dụng cho bài toán dự báo hạn dài và mô hình hoàn lưu khí quyển toàn cầu. Tích hợp định luật 4/3 Obukhov để xác định hệ số khuếch tán rối quy mô hành tinh: $(K_H)_0 = c \cdot D_0^{1/3} L_0^{4/3} \approx 8{,}75 \cdot 10^6\text{ m}^2/\text{s}$.
[Các định luật bảo toàn vật lý]
[Hệ phương trình cơ bản (x, y, z)]
[Mô hình hóa dòng rối & Lấy trung bình Reynolds]
[Chuyển đổi hệ tọa độ: (x, y, p, t) & (x, y, σ, t)]
[Phân tích quy mô (Scale Analysis) & Số Kibel K]
[Phương trình xoáy tựa địa chuyển & Phương trình cân bằng]
[Hệ tọa độ cầu cho mô hình toàn cầu & hạn dài]
Kiến thức nền tảng được xây dựng
- Lý thuyết thủy động lực học chất lưu địa cầu: Định luật bảo toàn mô men động lượng trên hệ quy chiếu quay, hiệu ứng tham số Coriolis $f = 2\omega \sin\varphi$ và biến thiên $\beta = \frac{\partial f}{\partial y}$.
- Nhiệt động lực học khí quyển: Khái niệm nhiệt độ thế vị $\theta = T (1000/p)^{R/c_p}$, nhiệt độ ảo $T_v$, và phương trình cân bằng tĩnh học.
- Lý thuyết nhiễu động và đồng dạng: Khung phân tích xấp xỉ địa chuyển ($u_g, v_g$), tính chất bất biến của xoáy tiềm năng, và điều kiện biên mô hình số.
Kỹ năng phát triển
- Kỹ năng toán lý chuyên sâu: Kỹ năng biến đổi tensor, vi phân hàm ẩn trong các hệ tọa độ cong không trực giao.
- Kỹ năng phân tích quy mô: Năng lực ước lượng bậc độ lớn của các số hạng đạo hàm để loại bỏ các sóng trọng trường hoặc nhiễu không mong muốn trong bài toán tích phân số.
- Kỹ năng thiết lập mô hình: Thiết lập bài toán biên ban đầu (Cauchy/Dirichlet) cho hệ phương trình khí quyển phân tầng.
Phương pháp giảng dạy và học tập
Phương pháp tiếp cận sư phạm
Giáo trình triển khai phương pháp diễn dịch toán lý: từ các định luật bảo toàn phổ quát đi đến các hệ phương trình vi phân đạo hàm riêng phi tuyến, sau đó áp dụng các giả thiết gần đúng vật lý (tựa tĩnh, tựa địa chuyển, tựa xô-lê-nô-ít) để đưa phương trình về dạng tích phân số khả thi. Mỗi bước biến đổi toán học đều đi kèm diễn giải ý nghĩa vật lý khí tượng cụ thể.
Bài tập và thực hành phân tích
- Phân tích thứ nguyên và bậc độ lớn: Sử dụng bảng số liệu đặc trưng của M. I. Yudin để tính toán định lượng tỷ số giữa lực quán tính và lực Coriolis, xác định giá trị số Kibel $K$ tại các vĩ độ khác nhau ($\varphi = 50^\circ, \varphi = 20^\circ$).
- Biến đổi giải tích: Thực hiện chứng minh toán học chuyển đổi hệ phương trình từ tọa độ Đề-các sang tọa độ $p$ và tọa độ địa hình $\sigma = P / P_x$, tính toán đại lượng vận tốc thẳng đứng $\omega = \frac{d\sigma}{dt}$.
- Khảo sát tính bảo toàn: Thiết lập nghiệm giải tích cho phương trình xoáy tuyệt đối $\frac{d(\Omega + f)}{dt} = 0$ và giải phương trình cân bằng Monge-Ampère $\nabla^2 \Phi = f \nabla^2 \psi + \dots$
Phương pháp kiểm tra đánh giá và tự học
- Đánh giá khả năng hiểu bản chất hệ phương trình khép kín và các tham số hóa vật lý (rối, bức xạ, ẩm).
- Kiểm tra năng lực phân tích điều kiện biên tại đáy khí quyển ($\sigma = 1 \Rightarrow \dot{\sigma} = 0$) và đỉnh tầng đối lưu.
- Hướng dẫn tự học: Sinh viên cần theo sát hệ thống ký hiệu hai chỉ số của công thức (ví dụ: công thức 1.15 là công thức thứ 15 trong mục 1) và thực hiện lại các bước đạo hàm trung gian để nắm vững logic toán học.
Điểm nổi bật và cập nhật học thuật
-
Khái quát hóa lịch sử phát triển khoa học dự báo số trị:
Hệ thống hóa toàn diện các mốc phát triển từ mô hình tính toán thủ công đầu tiên của Richardson (1922), qua bước nhảy vọt tìm ra nghiệm giải tích cho khí quyển tà áp của Buleev và Marchuk (1951), đến các mô hình máy tính điện tử của Charney, Phillips và Belousov.
-
Kỹ thuật xử lý địa hình bằng hệ tọa độ $\sigma$:
Tài liệu trình bày phương pháp giải quyết triệt để vấn đề biên dưới phức tạp của địa hình mặt đất thông qua biến đổi $\sigma = P / P_x$, cho phép đặt điều kiện biên vận tốc thẳng đứng $\dot{\sigma} = 0$ tại mặt đất chính xác hơn so với điều kiện biên trong hệ tọa độ $p$.
-
Mở rộng mô hình phi đoạn nhiệt và mô hình toàn cầu:
Tích hợp đầy đủ các dòng nhập nhiệt phi đoạn nhiệt (bức xạ $\varepsilon_{bx}$, rối $\varepsilon_R$, chuyển pha hơi nước $\varepsilon_{ph}$), mô tả phương trình cân bằng dòng solenoidal và thiết lập hệ phương trình thủy nhiệt động lực học trong hệ tọa độ cầu $(\lambda, \theta, r)$ ứng dụng cho dự báo hạn vừa và hạn dài.
Đối tượng sử dụng giáo trình
- Bậc đào tạo: Sinh viên đại học năm thứ 3, năm thứ 4 và học viên cao học, nghiên cứu sinh chuyên ngành Khí tượng học, Khí hậu học, Thủy văn, Hải dương học và Khoa học Môi trường.
- Kiến thức tiên quyết (Prerequisites):
- Giải tích toán học (phương trình vi phân đạo hàm riêng, giải tích vectơ và ten-xơ).
- Vật lý đại cương (cơ học chất lưu, nhiệt động lực học).
- Khí tượng động lực học đại cương và phương pháp tính/giải tích số.
- Mục đích sử dụng:
- Giảng viên: Sử dụng làm khung bài giảng chuẩn cho học phần Dự báo thời tiết số trị hoặc Động lực học khí quyển nâng cao.
- Cán bộ nghiên cứu và dự báo viên: Dùng làm tài liệu tra cứu các công thức lý thuyết nền tảng khi vận hành, hiệu chỉnh các mô hình số trị nghiệp vụ (như mô hình tựa địa chuyển, mô hình phương trình nguyên thủy).
Câu hỏi thường gặp
1. Giáo trình này phù hợp với đối tượng người học nào?
Tài liệu được thiết kế cho sinh viên chuyên ngành Khí tượng học từ năm thứ 3 trở lên, học viên cao học và các nghiên cứu sinh cần nắm vững bản chất toán lý của các mô hình dự báo thời tiết số trị.
2. Người học cần chuẩn bị kiến thức nền tảng nào trước khi tiếp cận?
Cần hoàn thành các học phần: Giải tích nhiều biến, Phương trình vi phân và Toán lý, Cơ học chất lỏng lý tưởng/chất lỏng thực, và học phần Khí tượng động lực học cơ sở.
3. Điểm khác biệt căn bản giữa hệ tọa độ áp suất $(p)$ và hệ tọa độ địa hình $(\sigma)$ trong giáo trình là gì?
Hệ tọa độ $p$ giúp loại bỏ ẩn số mật độ $\rho$ khỏi phương trình chuyển động và đơn giản hóa phương trình tĩnh học. Tuy nhiên, hệ tọa độ $\sigma = P / P_x$ chuẩn hóa mặt đất thành $\sigma = 1$, giúp xác định điều kiện biên dưới $\dot{\sigma} = 0$ chính xác tuyệt đối trên địa hình phức tạp.
4. Phương pháp đánh giá bậc đại lượng bằng số Kibel $K$ có vai trò gì trong việc xây dựng mô hình?
Số Kibel $K = \frac{U}{f L}$ (với $K \approx 0{,}104 \ll 1$ ở vĩ độ trung bình) là tiêu chuẩn không thứ nguyên cho phép lọc bỏ các thành phần có độ lớn nhỏ hơn $K$, dẫn xuất ra phương trình xoáy tựa địa chuyển nhằm loại trừ sóng âm và sóng trọng trường không mang thông tin dự báo synop.
5. Giáo trình có giới thiệu các mô hình dự báo thời tiết hạn dài không?
Có. Mục §6 của Chương I giới thiệu hệ phương trình thủy nhiệt động lực học trong hệ tọa độ cầu $(\lambda, \theta, r)$, tích hợp quy luật 4/3 Obukhov cho hệ số khuếch tán rối nằm ngang quy mô hành tinh phục vụ bài toán dự báo hạn dài toàn cầu.
Kết luận
Giáo trình Dự báo thời tiết bằng phương pháp số trị của tác giả Trần Tấn Tiến cung cấp hệ thống lý thuyết chuẩn mực về tích phân số hệ phương trình thủy nhiệt động lực học khí quyển. Tài liệu định hình lộ trình học tập logic từ việc nắm bắt các định luật bảo toàn cơ bản, kỹ thuật biến đổi hệ tọa độ tựa Lagorăng, cho đến phương pháp lọc sóng và phân tích quy mô. Đây là tài liệu học thuật cơ sở hỗ trợ việc nghiên cứu và phát triển các mô hình dự báo số trị nghiệp vụ trong ngành khí tượng thủy văn.
