Luận án Tiến sĩ: Bài toán ngược trong lý thuyết nhiệt - Phạm Hoàng Quân

Theo Alifanov, có bốn loại bài toán nhiệt ngược cơ bản. Bài toán ngược thời gian xác định nhiệt độ ban đầu từ dữ liệu thời điểm cuối. Bài toán biên ngược tìm phân bố nhiệt độ hoặc thông lượng nhiệt trên biên. Bài toán xác định hệ số khôi phục các tham số như hệ số dẫn nhiệt hay nguồn nhiệt. Bài toán hình học xác định đặc trưng hình học như lỗ hổng trong vật dẫn nhiệt. Luận án tập trung vào ba loại đầu tiên.

Bài toán gọi là chỉnh theo Hadamard khi nghiệm tồn tại, duy nhất và ổn định. Bài toán ngược nhiệt thường không thỏa ít nhất một điều kiện này. Tính không chỉnh xuất hiện do bản chất của phương trình khuếch tán. Nghiệm không ổn định nghĩa là sai số nhỏ trong dữ liệu gây biến động lớn trong nghiệm. Điều này làm bài toán khó giải bằng phương pháp số trực tiếp.

Bài toán ngược nhiệt có nhiều ứng dụng trong công nghiệp và khoa học. Xác định nhiệt độ bề mặt động cơ từ cảm biến bên trong là một ví dụ. Phát hiện khuyết tật trong vật liệu thông qua phân tích nhiệt là ứng dụng khác. Trong y học, phương pháp này giúp chẩn đoán bệnh qua hình ảnh nhiệt. Ngành hàng không vũ trụ sử dụng để thiết kế lá chắn nhiệt.

Mollification biến đổi dữ liệu thô thành hàm trơn hơn. Phương pháp sử dụng tích chập với hạt nhân đặc biệt. Hạt nhân Dirichlet và de la Vallée Poussin là hai lựa chọn phổ biến. Dữ liệu sau mollification thuộc không gian hàm nguyên với cấp hữu hạn. Trong không gian này, bài toán ban đầu trở thành chỉnh. Kỹ thuật này đặc biệt hiệu quả với phương trình parabolic.

Tham số chính quy hóa cân bằng giữa độ chính xác và ổn định. Tham số quá lớn cho nghiệm ổn định nhưng thiếu chính xác. Tham số quá nhỏ làm nghiệm nhạy cảm với nhiễu. Nguyên lý chọn hậu nghiệm của Morozov là một phương pháp phổ biến. Phương pháp L-curve cũng được sử dụng rộng rãi. Lựa chọn đúng tham số quyết định chất lượng nghiệm xấp xỉ.

Đánh giá sai số là bước quan trọng trong chính quy hóa. Sai số tổng gồm sai số xấp xỉ và sai số nhiễu. Định lý hội tụ chứng minh nghiệm chính quy hóa tiến tới nghiệm chính xác khi nhiễu giảm. Tốc độ hội tụ phụ thuộc vào độ trơn của nghiệm chính xác. Các đánh giá tiên nghiệm và hậu nghiệm cung cấp ước lượng sai số.

Bài toán Cauchy cho phương trình khuếch tán là không chỉnh nghiêm trọng. Nghiệm không ổn định liên tục theo dữ liệu trong mọi chuẩn thông thường. Sai số nhỏ trong dữ liệu có thể gây sai số lớn tùy ý trong nghiệm. Tính không chỉnh này do bản chất làm trơn của phương trình parabolic. Đặc trưng phổ của toán tử elliptic liên quan làm rõ vấn đề này.

Quasi-reversibility thay phương trình parabolic bằng phương trình nhiễu loạn. Phương trình mới có cấp cao hơn và tính chỉnh tốt hơn. Tham số nhiễu loạn đóng vai trò tham số chính quy hóa. Nghiệm của phương trình nhiễu loạn xấp xỉ nghiệm bài toán gốc. Phương pháp này được Lattès và Lions phát triển từ những năm 1960.

Phương pháp quasi-boundary value thay đổi điều kiện biên của bài toán. Điều kiện biên mới được chọn để bài toán trở nên chỉnh. Tham số trong điều kiện biên mới là tham số chính quy hóa. Phương pháp này đặc biệt hiệu quả cho bài toán ngược thời gian. Đánh giá sai số cho thấy tốc độ hội tụ logarit hoặc Hölder.

Nguồn nhiệt dạng q(x,t)f(x) với q chưa biết là trường hợp phổ biến. Hàm f(x) giả sử đã biết từ thông tin vật lý. Bài toán tìm q(x) từ phép đo nhiệt độ tại các điểm hoặc thời điểm. Điều kiện f(x) khác không trong miền là cần thiết. Isakov chỉ ra điều kiện này đảm bảo tính duy nhất. Không có điều kiện này, bài toán có thể có vô số nghiệm.

Hệ số dẫn nhiệt đặc trưng cho khả năng truyền nhiệt của vật liệu. Bài toán xác định hệ số này từ phép đo nhiệt độ là bài toán ngược phi tuyến. Dữ liệu cần phải đủ nhiều để bù đắp cho tính phi tuyến. Phương pháp gradient và Newton thường được áp dụng. Ổn định của nghiệm phụ thuộc vào miền biến thiên của hệ số.

Tính duy nhất đòi hỏi dữ liệu đủ phong phú và điều kiện tương thích. Đánh giá ổn định kiểu Lipschitz hoặc Hölder là mục tiêu mong muốn. Điều kiện trên hệ số như bị chặn trên dưới là cần thiết. Không gian hàm chứa nghiệm cần được chọn cẩn thận. Các kết quả ổn định có điều kiện (conditional stability) thường đạt được.

Phương trình truyền nhiệt được cho trong miền với điều kiện ban đầu. Điều kiện biên cho trên một phần biên có thể tiếp cận. Dữ liệu bổ sung là nhiệt độ đo tại các điểm bên trong hoặc trên biên khác. Cần tìm điều kiện biên trên phần biên còn lại. Bài toán này là bài toán không chỉnh do tính chất của phương trình parabolic.

Phương pháp phần tử hữu hạn rời rạc hóa miền không gian. Sai phân hữu hạn hoặc phương pháp đường thẳng xử lý biến thời gian. Hệ phương trình đại số thu được thường có ma trận bệnh. Chính quy hóa Tikhonov áp dụng cho hệ rời rạc. Phương pháp gradient liên hợp tìm nghiệm tối ưu.

Ngành luyện kim sử dụng để giám sát nhiệt độ lò. Công nghiệp điện tử áp dụng cho thiết kế tản nhiệt. Hàng không vũ trụ cần xác định nhiệt độ bề mặt tàu vũ trụ. Y học sử dụng trong điều trị ung thư bằng nhiệt. Mỗi ứng dụng có yêu cầu riêng về độ chính xác và tốc độ tính toán.

Luận án thiết lập các định lý duy nhất mới cho bài toán nguồn nhiệt. Đánh giá ổn định có điều kiện được cải thiện so với kết quả trước. Phương pháp chính quy hóa mới kết hợp mollification và quasi-boundary value. Thuật toán số hiệu quả được đề xuất và lập trình. Thí nghiệm số minh họa tính khả thi của các phương pháp.

Kết quả mở rộng các công trình của Isakov về bài toán nguồn. Cải tiến phương pháp mollification so với công trình Seidman. Bổ sung điều kiện yếu hơn cho tính duy nhất. Đánh giá sai số chính xác hơn trong một số trường hợp đặc biệt. Phương pháp số ổn định hơn các thuật toán trước đó.

Mở rộng cho phương trình truyền nhiệt phi tuyến là hướng tự nhiên. Bài toán trong miền có biên không trơn cần quan tâm. Hệ số gián đoạn mô tả vật liệu nhiều lớp là thách thức. Kết hợp phương pháp học sâu để tăng tốc tính toán. Ứng dụng vào bài toán thực tế phức tạp hơn cần được thử nghiệm.