Số biot là gì? Các bài báo nghiên cứu khoa học liên quan

Giới thiệu chung về số Biot

Số Biot, ký hiệu là Bi, là một đại lượng vô thứ nguyên giữ vai trò nền tảng trong các bài toán nhiệt động học và truyền nhiệt không ổn định. Đại lượng này biểu thị tỷ số giữa trở kháng nhiệt dẫn truyền bên trong vật thể rắn và trở kháng nhiệt đối lưu tỏa ra môi trường xung quanh tại bề mặt tiếp xúc. Nói cách khác, số Biot phản ánh mức độ chênh lệch nhiệt độ tương đối giữa tâm lõi của vật thể và lớp bề mặt khi có sự trao đổi năng lượng nhiệt với môi trường bao quanh.

Trong tính toán kỹ thuật nhiệt, việc đánh giá chính xác giá trị của số Biot cho phép các kỹ sư và nhà nghiên cứu xác định phương pháp mô hình hóa nhiệt phù hợp nhất, từ phương pháp đơn giản hóa với giả định phân bố nhiệt độ đồng đều trong toàn khối đến việc giải hệ phương trình vi phân đạo hàm riêng mô tả trường nhiệt độ phụ thuộc không gian và thời gian.

Các lĩnh vực kỹ thuật và công nghệ ứng dụng then chốt của số Biot bao gồm:

• Kỹ thuật nhiệt, năng lượng và công nghệ điện lạnh
• Thiết kế và tối ưu hóa hiệu suất của các thiết bị trao đổi nhiệt
• Kỹ thuật hàng không, vũ trụ và bảo vệ nhiệt khí động học
• Cơ sinh học, y sinh học và kỹ thuật nhiệt trị trong y khoa
• Khoa học vật liệu, luyện kim và các quy trình xử lý nhiệt tôi phôi

Biểu thức toán học của số Biot

Số Biot được xác định theo công thức định lượng:

$Bi = \frac{hL_c}{k}$

Trong đó:

• h: hệ số truyền nhiệt đối lưu tại bề mặt tiếp xúc giữa vật thể và môi chất xung quanh (đơn vị: W/m²·K)
• $L_c$: chiều dài đặc trưng về mặt hình học của vật thể rắn (đơn vị: m)
• k: hệ số dẫn nhiệt của vật liệu cấu thành vật thể rắn (đơn vị: W/m·K)

Biểu thức trên chứng minh rằng số Biot là thước đo tương đối giữa tốc độ truyền nhiệt qua ranh giới bề mặt và tốc độ vận chuyển nhiệt bên trong cấu trúc khối vật liệu. Khi vật thể có khả năng dẫn nhiệt nội tại rất cao nhưng trao đổi nhiệt bề mặt bị hạn chế, số Biot sẽ nhận giá trị rất nhỏ. Ngược lại, khi vật liệu có tính cách nhiệt cao (hệ số dẫn nhiệt nhỏ) trong khi trao đổi nhiệt cưỡng bức tại bề mặt diễn ra mãnh liệt, số Biot sẽ nhận giá trị lớn.

Ý nghĩa vật lý và tiêu chuẩn phân loại theo số Biot

Giá trị của số Biot đóng vai trò là tiêu chuẩn quyết định để xác định xem gradient nhiệt độ bên trong một cấu trúc vật thể rắn có thể được bỏ qua một cách an toàn hay không:

• Khi $Bi \le 0.1$: Trở kháng nhiệt nội tại nhỏ hơn rất nhiều so với trở kháng đối lưu bề mặt. Gradient nhiệt độ bên trong khối vật liệu là không đáng kể, do đó toàn bộ vật thể có thể xem như đồng nhất về nhiệt độ tại mỗi thời điểm. Sai số tính toán khi áp dụng mô hình khối lượng đồng nhất (lumped capacitance method) được duy trì ở mức dưới 5%.
• Khi số Biot xấp xỉ bậc đơn vị hoặc lớn hơn 0.1: Sự chênh lệch nhiệt độ giữa vùng trung tâm và vùng biên bắt đầu tăng cao. Việc xem vật thể là đẳng nhiệt không còn chính xác, đòi hỏi phải giải phương trình vi phân dẫn nhiệt nhiều chiều phụ thuộc tọa độ không gian.

Khoảng giá trị của Bi	Mức độ kháng nhiệt tương đối	Đặc điểm trường nhiệt độ	Phương pháp tính toán phù hợp
Số Biot nhỏ hơn hoặc bằng 0.1	Trở kháng dẫn nhiệt nội tại không đáng kể so với trở kháng đối lưu bề mặt	Nhiệt độ xem như đồng đều trong toàn bộ thể tích vật thể	Mô hình nhiệt dung gộp đơn giản (sai số dưới 5%)
Số Biot lớn hơn 0.1	Trở kháng dẫn nhiệt nội tại đáng kể hoặc vượt trội trở kháng đối lưu	Tồn tại gradient nhiệt độ rõ rệt giữa tâm và bề mặt	Nghiệm giải tích chuỗi một số hạng hoặc phương pháp phần tử hữu hạn số

Chiều dài đặc trưng và vai trò hình học

Chiều dài đặc trưng $L_c$ là thông số hình học cốt lõi trong công thức tính số Biot. Đối với một vật thể bất kỳ, chiều dài đặc trưng thường được định nghĩa bằng tỷ số giữa thể tích và diện tích bề mặt trao đổi nhiệt:

$L_c = \frac{V}{A}$

Trong đó:

• V: tổng thể tích hình học của vật thể rắn
• A: tổng diện tích bề mặt tiếp xúc trực tiếp tham gia trao đổi nhiệt đối lưu

Đối với các dạng hình học cơ bản trong kỹ thuật nhiệt, chiều dài đặc trưng $L_c$ được suy biến thành các biểu thức giải tích tường minh:

• Khối cầu bán kính R: $L_c = \frac{R}{3}$
• Khối lập phương cạnh a: $L_c = \frac{a^3}{6a^2} = \frac{a}{6}$
• Khối trụ dài bán kính R và chiều dài L: $L_c = \frac{\pi R^2 L}{2\pi R L + 2\pi R^2}$

Việc thiết lập chuẩn xác chiều dài đặc trưng theo đúng hình dạng tiếp xúc nhiệt bảo đảm cho kết quả tính toán số Biot phản ánh trung thực bản chất vật lý của quá trình trao đổi nhiệt.

Mô hình nhiệt dung gộp (Lumped Capacitance Method)

Khi thỏa mãn điều kiện số Biot nhỏ hơn hoặc bằng 0.1, các kỹ sư có thể áp dụng mô hình nhiệt dung gộp (còn gọi là mô hình khối lượng đồng nhất) để mô tả quá trình biến thiên nhiệt độ theo thời gian. Mô hình này giả thiết rằng nhiệt độ tại mọi điểm trong vật thể là như nhau tại cùng một thời điểm, loại bỏ sự phụ thuộc không gian của trường nhiệt độ.

Sự suy giảm hoặc gia tăng nhiệt độ của vật thể theo thời gian được mô tả bằng hàm mũ:

$T(t) = T_{\infty} + (T_i - T_{\infty}) e^{-\frac{hA}{\rho c_p V}t}$

Trong đó:

• T(t): nhiệt độ của vật thể tại thời điểm t
• T_i: nhiệt độ ban đầu đồng nhất của vật thể rắn
• T_∞: nhiệt độ của môi chất xung quanh ở xa bề mặt
• ρ: khối lượng riêng của vật liệu cấu tạo
• c_p: nhiệt dung riêng đẳng áp của vật liệu
• V, A: thể tích và diện tích bề mặt trao đổi nhiệt của cấu trúc

Mô hình này giúp đơn giản hóa việc giải phương trình vi phân đạo hàm riêng phi tuyến thành phương trình vi phân thường bậc nhất tuyến tính, hỗ trợ đắc lực cho việc đánh giá thời gian nung nóng, thời gian làm nguội phôi kim loại và phản ứng của các đầu đo nhiệt điện ngâm trong dòng lưu chất.

So sánh số Biot với các số vô thứ nguyên nhiệt học

Trong phân tích truyền nhiệt quá độ, số Biot luôn được khảo sát đồng thời cùng các đại lượng vô thứ nguyên quan trọng khác:

Tên đại lượng	Ký hiệu	Ý nghĩa vật lý	Biểu thức định nghĩa
Số Biot	Bi	Tỷ số giữa trở kháng nhiệt dẫn truyền nội tại và trở kháng nhiệt đối lưu bề mặt	$Bi = \frac{hL_c}{k}$
Số Fourier	Fo	Thước đo không thứ nguyên về tốc độ khuếch tán nhiệt độ theo thời gian	$Fo = \frac{\alpha t}{L_c^2}$
Số Nusselt	Nu	Tỷ số giữa cường độ truyền nhiệt đối lưu và dẫn nhiệt thuần túy trong lớp lưu chất	$Nu = \frac{hL}{k_f}$

Cần lưu ý sự phân biệt rõ ràng giữa số Biot và số Nusselt dù có biểu thức hình thức tương tự nhau: trong số Biot, hệ số dẫn nhiệt k là của vật thể rắn; trong khi ở số Nusselt, hệ số dẫn nhiệt k_f là của môi chất lưu chất bao quanh. Khi số Fourier lớn hơn 0.2, các nghiệm chuỗi vô hạn trong truyền nhiệt quá độ qua các kết cấu phẳng, trụ và cầu hội tụ nhanh chóng và có thể được xấp xỉ chính xác bằng số hạng giải tích đầu tiên kết hợp với số Biot.

Ứng dụng số Biot trong mô phỏng số và kỹ thuật tính toán

Trong các gói phần mềm mô phỏng trường nhiệt vật lý như COMSOL Multiphysics, ANSYS Fluent hoặc OpenFOAM, số Biot được sử dụng làm tiêu chí để phân loại bài toán và lập chiến lược chia lưới. Khi số Biot được xác định là nhỏ hơn 0.1, hệ thống có thể kích hoạt mô hình zero chiều (0D lumped model) để tiết kiệm tài nguyên tính toán mà vẫn đảm bảo sai số dưới 5%. Ngược lại, khi số Biot vượt quá ngưỡng này, phần mềm buộc phải thiết lập lưới không gian chi tiết để giải phương trình vi phân truyền nhiệt tổng quát:

$\rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + Q$

với điều kiện biên hỗn hợp loại ba (điều kiện biên Robin) thể hiện qua số Biot tại bề mặt vật thể:

$-k \frac{\partial T}{\partial n} = h(T - T_{\infty})$

Quy trình đánh giá sơ bộ số Biot trước khi khởi chạy mô phỏng 3D giúp tối ưu hóa thời gian tính toán và hạn chế các lỗi xấp xỉ trong thiết kế tản nhiệt điện tử, van nhiệt lò phản ứng và cấu kiện hàng không.