Lỗ nhỏ và dòng nghẹt là hai điều kiện khác nhau
Một khe rất nhỏ có thể tạo lưu lượng thấp, nhưng kích thước nhỏ tự nó chưa chứng minh dòng nghẹt. Tỷ số áp suất, hình dạng đường rò, ma sát và chế độ liên tục của khí cùng quyết định mô hình. Trong bài này, “khe” được thay bằng một lỗ ngắn tương đương có diện tích và hệ số lưu lượng đã giả định. Kết quả nhằm nghiên cứu độ nhạy, không suy ra kích thước vết nứt thực từ một phép đo duy nhất.
Phân tích của Kandula và Nufer xem xét cả hình học đường rò và sự thay đổi chế độ ở kích thước rất nhỏ; hiệu ứng khí loãng trở nên đáng kể trong miền phù hợp. [2] Vì vậy, mô hình lỗ lý tưởng dưới đây không được kéo dài tự động sang khe dài, mạng mao quản hoặc rò phân tử trong chân không.
Công thức dùng áp suất tuyệt đối
Gọi p_u và T_u là áp suất tuyệt đối, Pa, và nhiệt độ, K, trong bình lớn phía trước lỗ; vận tốc trong bình coi nhỏ. Áp suất ngoài p_d là tuyệt đối. Đặt β=p_d/p_u, γ=c_p/c_v, R_s là hằng số khí riêng, J/(kg·K), A là diện tích lỗ, m², và C_d không thứ nguyên. Tỷ số tới hạn của mô hình khí lý tưởng, nhiệt dung không đổi:
β_c=[2/(γ+1)]^[γ/(γ−1)].
Khi β≤β_c, lưu lượng khối lượng ra khỏi bình là:
m_dot=C_d A p_u√[γ/(R_sT_u)] [2/(γ+1)]^[(γ+1)/(2(γ−1))].
Khi β_c<β<1, dùng nhánh chưa nghẹt:
m_dot=C_d A p_u√{2γ/[(γ−1)R_sT_u] × [β^(2/γ)−β^((γ+1)/γ)]}.
Hai nhánh có cùng giá trị tại β_c. Tài liệu mô hình lỗ của MathWorks trình bày các quan hệ và nêu C_d phụ thuộc hình học. [1] Đây là tài liệu mô hình, không phải lựa chọn một giá trị C_d được chứng nhận cho khe đang khảo sát. Hạ p_d thêm khi đã nghẹt không làm tăng lưu lượng nếu p_u, T_u và lỗ giữ nguyên; tăng p_u vẫn có thể làm tăng lưu lượng.
Trạng thái ban đầu của một bình giả định
Cho V=0,100 m³, p_u0=800.000 Pa tuyệt đối, T_0=300 K, p_d=100.000 Pa tuyệt đối, R_s=287 J/(kg·K), γ=1,40. Lỗ tròn tương đương d=0,00100 m có A=7,8539816×10⁻⁷ m²; chọn C_d=0,650 chỉ cho ví dụ.
Tính được β_c=0,528281788. Ban đầu β=0,125 nên nhánh nghẹt áp dụng trong mô hình. Hệ số m_dot/p_u bằng 1,191302195×10⁻⁹ kg/(s·Pa), cho m_dot0=0,0009530418 kg/s. Khối lượng khí ban đầu m_0=p_u0V/(R_sT_0)=0,92915215 kg.
Không lấy 700.000 Pa áp suất tương đối thay cho p_u0. Cũng không lấy tỷ số chênh áp/p_u để kiểm β_c: β ở đây là áp suất hạ lưu chia áp suất thượng lưu. Hai lỗi có thể cùng tạo một con số trông hợp lý nhưng làm sai cả quyết định chọn nhánh lẫn lưu lượng.
Tích trữ làm áp suất và lưu lượng cùng giảm
Giả định trước tiên thành bình cấp nhiệt đủ nhanh để khí trong bình giữ T_0. Trong đoạn nghẹt, viết m_dot=Cp với C như trên. Cân bằng dm/dt=−Cp và m=pV/(R_sT_0) dẫn đến:
dp/dt=−p/τ; τ=V/(R_sT_0C); p(t)=p_0 exp(−t/τ).
τ=974,933305 s. Sau 600 s, áp suất còn 432.328,196 Pa và vẫn ở miền nghẹt. Lưu lượng tại thời điểm đó giảm cùng tỷ lệ áp suất. Lấy lưu lượng ban đầu nhân 600 s sẽ tính lượng mất lớn hơn nghiệm tích phân, vì nó giả định nguồn áp không suy giảm.
Dòng hết nghẹt khi p đạt p_c=p_d/β_c=189.292,916 Pa. Thời điểm này t_c=τ ln(p_0/p_c)=1.405,187062 s. Không tiếp tục hàm m_dot=Cp xuống áp suất khí quyển: sau t_c, tỷ số β đã vượt β_c và lưu lượng phải tính bằng nhánh chưa nghẹt.
Để tìm thời gian từ p_c xuống 120.000 Pa, tích phân dt=−V dp/[R_sT_0 m_dot(p)] với công thức nhánh chưa nghẹt. Tích phân số cho 478,296860 s, nên thời gian từ đầu là 1.883,483921 s. Lưu lượng tại 120.000 Pa còn 0,0001092516 kg/s. Các số này là nghiệm mô hình đẳng nhiệt và C_d không đổi; không phải lịch giảm áp vận hành được khuyến nghị.
Đẳng nhiệt của bình không đồng nghĩa không trao đổi nhiệt
Một lỗ lý tưởng có thể được mô hình như dòng giãn nở rất nhanh, trong khi bình phía trước trao đổi nhiệt với thành và môi trường. Hai thể tích kiểm soát khác nhau nên không thể gọi toàn hệ “đoạn nhiệt” chỉ vì công thức lỗ bắt đầu từ quan hệ giãn nở.
Với bình cứng và khí lý tưởng có T không đổi, cân bằng năng lượng cho Q_dot=m_dot R_sT nếu bỏ công trục và động năng trong bình. Ban đầu cần Q_dot=82,056895 W truyền vào khí để giữ đúng giả thiết 300 K. Đây là nhu cầu nhiệt của mô hình; chưa có hệ số truyền nhiệt hay nhiệt dung thành bình để chứng minh thiết bị thực đáp ứng nó. Mô hình buồng khí của MathWorks giữ riêng cân bằng khối lượng và năng lượng vì cả p lẫn T có thể thay đổi. [3]
Ở giới hạn khác, một bình khí đồng nhất, cứng, đoạn nhiệt và không có dòng vào cho T/T_0=(p/p_0)^[(γ−1)/γ]. Tại cùng áp suất 189.292,916 Pa, nhiệt độ mô hình này là 198,736112 K. Không được lấy thời điểm 1.405,187062 s của nghiệm đẳng nhiệt gán cho trường hợp đoạn nhiệt: lưu lượng phụ thuộc nhiệt độ và lịch biến đổi cũng khác. Khí thật còn có thể đổi tính chất hoặc trạng thái pha trước khi đạt điều kiện này.
Độ nhạy theo lỗ và nhiệt độ đo
Trong cùng nhánh nghẹt và cùng điều kiện bình, m_dot tỷ lệ C_d d² p/√T. Tăng d 10% làm diện tích và lưu lượng tăng 21%, còn τ giảm từ 974,933305 xuống 805,730004 s nếu mọi giả thiết khác giữ nguyên. Hệ số C_d và diện tích thường không thể nhận dạng riêng chỉ từ một đường suy giảm áp: dữ liệu chủ yếu xác định tích C_dA trong mô hình này.
Một kiểm tra khác xét bình hoàn toàn kín nhưng nguội từ 300 xuống 294 K, giữ V và m. Áp suất giảm từ 800.000 xuống 784.000 Pa. Nếu người xử lý giữ nhiệt độ 300 K, họ sẽ suy ra mất 0,018583043 kg, dù khối lượng thực không đổi. Nếu quá trình kéo dài 60 s, “lưu lượng rò” giả sẽ là 0,0003097174 kg/s. Cần đồng bộ p và T, và phải biết nhiệt độ đo đại diện cho khí hay chỉ cho thành.
Phạm vi của kết luận
Diện tích tương đương không xác định được hình dáng vết rò; hai khe có cùng diện tích nhưng chiều dài và bề mặt khác có thể cho lưu lượng khác. Dòng nhiều pha, biến dạng khe theo áp, trao đổi nhiệt mạnh trong mao quản hoặc khí loãng đòi hỏi mô hình bổ sung. Bài NASA dùng các chế độ khác nhau chính vì một quan hệ không bao phủ mọi kích thước và điều kiện. [2]
Giá trị của phép tính là kiểm nhất quán áp suất tuyệt đối, thời điểm đổi nhánh và khối lượng tích trữ. Một đánh giá rò đáng tin còn cần trạng thái nhiệt, thể tích thực nối với phép đo, tính chất khí và dữ liệu hình học hoặc hiệu chuẩn lưu lượng tương ứng. Nó không thay thế quy trình thử kín và tiêu chí chấp nhận riêng của hệ thống.
Nguồn tham khảo
[1] Constant Area Pneumatic Orifice — Sharp-edged orifice in pneumatic systems — MathWorks, tài liệu mô hình Pneumatic legacy. Truy cập ngày 09/11/2026.
[2] On the Equivalence of Gas and Liquid Leak Flow Rates for Application to Propulsion and Propellant Transfer Systems — Max Kandula và Brian M. Nufer; NASA Kennedy Space Center. Truy cập ngày 09/11/2026.
[3] Constant Volume Chamber (G) — Chamber with fixed volume of gas and variable number of ports — MathWorks. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một khe rất nhỏ có thể tạo lưu lượng thấp, nhưng kích thước nhỏ tự nó chưa chứng minh dòng nghẹt. Tỷ số áp suất, hình dạng đường rò, ma sát và chế độ liên tục của khí cùng quyết định mô hình. Trong bài này, “khe” được thay bằng một lỗ ngắn tương đương có diện tích và hệ số lưu lượng đã giả định. Kết quả nhằm nghiên cứu độ nhạy, không suy ra kích thước vết nứt thực từ một phép đo duy nhất.
Phân tích của Kandula và Nufer xem xét cả hình học đường rò và sự thay đổi chế độ ở kích thước rất nhỏ; hiệu ứng khí loãng trở nên đáng kể trong miền phù hợp. [2] Vì vậy, mô hình lỗ lý tưởng dưới đây không được kéo dài tự động sang khe dài, mạng mao quản hoặc rò phân tử trong chân không.
Công thức dùng áp suất tuyệt đối
Gọi p_u và T_u là áp suất tuyệt đối, Pa, và nhiệt độ, K, trong bình lớn phía trước lỗ; vận tốc trong bình coi nhỏ. Áp suất ngoài p_d là tuyệt đối. Đặt β=p_d/p_u, γ=c_p/c_v, R_s là hằng số khí riêng, J/(kg·K), A là diện tích lỗ, m², và C_d không thứ nguyên. Tỷ số tới hạn của mô hình khí lý tưởng, nhiệt dung không đổi:
β_c=[2/(γ+1)]^[γ/(γ−1)].
Khi β≤β_c, lưu lượng khối lượng ra khỏi bình là:
m_dot=C_d A p_u√[γ/(R_sT_u)] [2/(γ+1)]^[(γ+1)/(2(γ−1))].
Khi β_c<β<1, dùng nhánh chưa nghẹt:
m_dot=C_d A p_u√{2γ/[(γ−1)R_sT_u] × [β^(2/γ)−β^((γ+1)/γ)]}.
Hai nhánh có cùng giá trị tại β_c. Tài liệu mô hình lỗ của MathWorks trình bày các quan hệ và nêu C_d phụ thuộc hình học. [1] Đây là tài liệu mô hình, không phải lựa chọn một giá trị C_d được chứng nhận cho khe đang khảo sát. Hạ p_d thêm khi đã nghẹt không làm tăng lưu lượng nếu p_u, T_u và lỗ giữ nguyên; tăng p_u vẫn có thể làm tăng lưu lượng.
Trạng thái ban đầu của một bình giả định
Cho V=0,100 m³, p_u0=800.000 Pa tuyệt đối, T_0=300 K, p_d=100.000 Pa tuyệt đối, R_s=287 J/(kg·K), γ=1,40. Lỗ tròn tương đương d=0,00100 m có A=7,8539816×10⁻⁷ m²; chọn C_d=0,650 chỉ cho ví dụ.
Tính được β_c=0,528281788. Ban đầu β=0,125 nên nhánh nghẹt áp dụng trong mô hình. Hệ số m_dot/p_u bằng 1,191302195×10⁻⁹ kg/(s·Pa), cho m_dot0=0,0009530418 kg/s. Khối lượng khí ban đầu m_0=p_u0V/(R_sT_0)=0,92915215 kg.
Không lấy 700.000 Pa áp suất tương đối thay cho p_u0. Cũng không lấy tỷ số chênh áp/p_u để kiểm β_c: β ở đây là áp suất hạ lưu chia áp suất thượng lưu. Hai lỗi có thể cùng tạo một con số trông hợp lý nhưng làm sai cả quyết định chọn nhánh lẫn lưu lượng.
Tích trữ làm áp suất và lưu lượng cùng giảm
Giả định trước tiên thành bình cấp nhiệt đủ nhanh để khí trong bình giữ T_0. Trong đoạn nghẹt, viết m_dot=Cp với C như trên. Cân bằng dm/dt=−Cp và m=pV/(R_sT_0) dẫn đến:
dp/dt=−p/τ; τ=V/(R_sT_0C); p(t)=p_0 exp(−t/τ).
τ=974,933305 s. Sau 600 s, áp suất còn 432.328,196 Pa và vẫn ở miền nghẹt. Lưu lượng tại thời điểm đó giảm cùng tỷ lệ áp suất. Lấy lưu lượng ban đầu nhân 600 s sẽ tính lượng mất lớn hơn nghiệm tích phân, vì nó giả định nguồn áp không suy giảm.
Dòng hết nghẹt khi p đạt p_c=p_d/β_c=189.292,916 Pa. Thời điểm này t_c=τ ln(p_0/p_c)=1.405,187062 s. Không tiếp tục hàm m_dot=Cp xuống áp suất khí quyển: sau t_c, tỷ số β đã vượt β_c và lưu lượng phải tính bằng nhánh chưa nghẹt.
Để tìm thời gian từ p_c xuống 120.000 Pa, tích phân dt=−V dp/[R_sT_0 m_dot(p)] với công thức nhánh chưa nghẹt. Tích phân số cho 478,296860 s, nên thời gian từ đầu là 1.883,483921 s. Lưu lượng tại 120.000 Pa còn 0,0001092516 kg/s. Các số này là nghiệm mô hình đẳng nhiệt và C_d không đổi; không phải lịch giảm áp vận hành được khuyến nghị.
Đẳng nhiệt của bình không đồng nghĩa không trao đổi nhiệt
Một lỗ lý tưởng có thể được mô hình như dòng giãn nở rất nhanh, trong khi bình phía trước trao đổi nhiệt với thành và môi trường. Hai thể tích kiểm soát khác nhau nên không thể gọi toàn hệ “đoạn nhiệt” chỉ vì công thức lỗ bắt đầu từ quan hệ giãn nở.
Với bình cứng và khí lý tưởng có T không đổi, cân bằng năng lượng cho Q_dot=m_dot R_sT nếu bỏ công trục và động năng trong bình. Ban đầu cần Q_dot=82,056895 W truyền vào khí để giữ đúng giả thiết 300 K. Đây là nhu cầu nhiệt của mô hình; chưa có hệ số truyền nhiệt hay nhiệt dung thành bình để chứng minh thiết bị thực đáp ứng nó. Mô hình buồng khí của MathWorks giữ riêng cân bằng khối lượng và năng lượng vì cả p lẫn T có thể thay đổi. [3]
Ở giới hạn khác, một bình khí đồng nhất, cứng, đoạn nhiệt và không có dòng vào cho T/T_0=(p/p_0)^[(γ−1)/γ]. Tại cùng áp suất 189.292,916 Pa, nhiệt độ mô hình này là 198,736112 K. Không được lấy thời điểm 1.405,187062 s của nghiệm đẳng nhiệt gán cho trường hợp đoạn nhiệt: lưu lượng phụ thuộc nhiệt độ và lịch biến đổi cũng khác. Khí thật còn có thể đổi tính chất hoặc trạng thái pha trước khi đạt điều kiện này.
Độ nhạy theo lỗ và nhiệt độ đo
Trong cùng nhánh nghẹt và cùng điều kiện bình, m_dot tỷ lệ C_d d² p/√T. Tăng d 10% làm diện tích và lưu lượng tăng 21%, còn τ giảm từ 974,933305 xuống 805,730004 s nếu mọi giả thiết khác giữ nguyên. Hệ số C_d và diện tích thường không thể nhận dạng riêng chỉ từ một đường suy giảm áp: dữ liệu chủ yếu xác định tích C_dA trong mô hình này.
Một kiểm tra khác xét bình hoàn toàn kín nhưng nguội từ 300 xuống 294 K, giữ V và m. Áp suất giảm từ 800.000 xuống 784.000 Pa. Nếu người xử lý giữ nhiệt độ 300 K, họ sẽ suy ra mất 0,018583043 kg, dù khối lượng thực không đổi. Nếu quá trình kéo dài 60 s, “lưu lượng rò” giả sẽ là 0,0003097174 kg/s. Cần đồng bộ p và T, và phải biết nhiệt độ đo đại diện cho khí hay chỉ cho thành.
Phạm vi của kết luận
Diện tích tương đương không xác định được hình dáng vết rò; hai khe có cùng diện tích nhưng chiều dài và bề mặt khác có thể cho lưu lượng khác. Dòng nhiều pha, biến dạng khe theo áp, trao đổi nhiệt mạnh trong mao quản hoặc khí loãng đòi hỏi mô hình bổ sung. Bài NASA dùng các chế độ khác nhau chính vì một quan hệ không bao phủ mọi kích thước và điều kiện. [2]
Giá trị của phép tính là kiểm nhất quán áp suất tuyệt đối, thời điểm đổi nhánh và khối lượng tích trữ. Một đánh giá rò đáng tin còn cần trạng thái nhiệt, thể tích thực nối với phép đo, tính chất khí và dữ liệu hình học hoặc hiệu chuẩn lưu lượng tương ứng. Nó không thay thế quy trình thử kín và tiêu chí chấp nhận riêng của hệ thống.
Nguồn tham khảo
[1] Constant Area Pneumatic Orifice — Sharp-edged orifice in pneumatic systems — MathWorks, tài liệu mô hình Pneumatic legacy. Truy cập ngày 09/11/2026.
[2] On the Equivalence of Gas and Liquid Leak Flow Rates for Application to Propulsion and Propellant Transfer Systems — Max Kandula và Brian M. Nufer; NASA Kennedy Space Center. Truy cập ngày 09/11/2026.
[3] Constant Volume Chamber (G) — Chamber with fixed volume of gas and variable number of ports — MathWorks. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
