浅谈反应扩散方程.docxVIP

PAGE PAGE #数學傳播39卷4期民104年12月 PAGE PAGE #数學傳播39卷4期民104年12月 数學傳播 数學傳播39卷4期.pp. PAGE #?26 淺談反應擴散方程 倪維明 從物理、化學到生物，反應擴敞方程在近代科學中已被廣泛地用來描述各種現象。披散，如 我們下面第一節所談的，似乎是一種「平凡化」的過程：任何初始狀態，經過艮時間擴散.總是 達到處處常數的境界。因此，在1952年，Alan Turing用反應擄散方程來描述自然界中的種種 pattern formation,是個劇時代的創舉！ 原始形態學中最早的一個例子，應該是A. TYembley在1744年發現的水蛔｛hydra)的 再生現象：當水蝗的頭部被切掉48小時後，又能長出一個新的頭。爲了瞭解這一奇特的現象, Turing在1952年提出了一個關於activator-inhibitor的數學模型。在這個理論中，Turing 假設activator擴散的很慢但inhibitor擴散的很快、而他充分利用這兩者之間擴散速度的差 異.給出一個「直觀J的論證，來說明水媳頭部的再生現象。他的觀點在1972年經過數值模擬, 終於被formulated成一個非線性方程組，並引發了大量數學上的研究。 1998年在Notices of AMS上筆者有一篇文章[13],介紹了在「自律(autonomous)J系 統中，反應與搪散的互動及其所產生的各種不同的機制,包括了有趣的「凝聚(condenSat.ion)J 現象，例如：尖峰解(peak solutions),也簡要說明了上述Alan Turning原創性觀點的數學 證明。 在這篇短文中，我們將由擴散的機制談起，但我們的後績討論著重於「非自律」系統中，反 應、擴散輿空間異質性(spatial heterogeneity)之間的交互作用。我們將以生態學中古老的 logistic equation爲例，來闡明這三者之間的互動關係。 要特別一提的是，雖然本文以引發讀者的興趣爲主，但文中所選的題材，都是基本、簡單而 又重要的，同時包含了當前研究前沿中較受關注的open problems.感興趣的讀者，歡迎與筆 者聯繫，或直接嘗試解決這些問題。 1.首先我們由熱方程(heat equation)談起 假設在空間Rn中我們有一有界區域Q,它光滑的邊界是由良好的絕熱材料製成。若將Q 17 內部在點与時間tQ的部在,記爲u(x.t),則u(x.t)滿足下列方程 ut ut = △?/ 0) = Uo(x) dyu = 0 in Q x (0)oo), in Q. on dQ X (0, oo), 此處△ = 讀是Laplwe運算元，。代表知上的單位外法線方向，邊值條件 M = 0 即代表了這是個封爾系統：無熱流入也無熱流出。爲了簡單起見，在本文中，我們將始終假設 s N 0 且 “0 第 0。 顕而易見的(我們也已提過)，當我們讓時間t T 8時，不論初始值心(了)如何分佈, “(/,￡)必定趨於一個常散布一即u0(x)在Q上的平均值。數學上來說，我們有 u(x, t)=丽 + ( Z UQ物)物(])+ ... 其中0 = 1 V 2 W ”3 4 ...是△的一串固有值(eigenvalues),它們對應的正規化固有函 教(normalized eigenfunctions)爲 物三常晚 响,洌… ( △我+四四=0 in Q, %Wk =。 on dQ. 事實上，(2)式不但告訴了我們解u(x.t)趨於一個常數uo (初始值心(时的平均值)，它 還指出趨於可的速度(rate)是由△的第一個非零的固有值“2所決定。 幸或不辛，有趣的是，到目前爲止，數學家們還沒有找到決定M2大小的有效方法，但是, 許多人甚至以爲決定Neumann固有值网大小的方式與決定Dirichlet |S1有值大小的方式是 一樣或類似的° ([15, pp.308-309], |10, p.239]) Dirichlet固有值問題： (4)in Q, on 00., (4) 滿足一個簡單的性質: 「区域單訶性(Domain-Monotonicity Property)J:若 Qi C Q2｝则入(%) Afc(fi2),封所有A; = 1,2,...都成立。 在2007年筆者與王學鋒的文章[14]裡，對這類問題有一些初步的討論。當然，上述敎科 書([15], [10|)在[14] 一文發表後，都做了適度的更正。 淺談反應损散方程 淺談反應损散方程19 淺談反應损散方程 淺談反應损散方程19 在［14］中，作者提出了幾個相關的幾何量：例如:「內在直徑(intrinsic diametcrjj, ... 等。我們先來看個例子。設