Lax – Wendroff yöntemi - Lax–Wendroff method - Wikipedia

Lax – Wendroff yöntemi, adını Peter Lax ve Burton Wendroff, bir sayısal çözümü için yöntem hiperbolik kısmi diferansiyel denklemler, dayalı sonlu farklar. Hem uzay hem de zamanda ikinci dereceden doğrudur. Bu yöntem bir örnektir açık zaman entegrasyonu yönetim denklemini tanımlayan fonksiyon mevcut zamanda değerlendirilir.

Tanım

Denklemin aşağıdaki formda olduğunu varsayalım:

${ displaystyle { frac { kısmi u (x, t)} { kısmi t}} + { frac { kısmi f (u (x, t))} { kısmi x}} = 0}$

nerede x ve t bağımsız değişkenlerdir ve başlangıç ​​durumu, u (x, 0) verilir.

Doğrusal durum

Doğrusal durumda, nerede f (u) = Au , ve Bir sabittir^[1]

${ displaystyle u_ {i} ^ {n + 1} = u_ {i} ^ {n} - { frac { Delta t} {2 Delta x}} A sol [u_ {i + 1} ^ { n} -u_ {i-1} ^ {n} right] + { frac { Delta t ^ {2}} {2 Delta x ^ {2}}} A ^ {2} left [u_ { i + 1} ^ {n} -2u_ {i} ^ {n} + u_ {i-1} ^ {n} sağ].}$

Bu doğrusal şema, farklı şekillerde genel doğrusal olmayan duruma genişletilebilir. Bunlardan biri izin veriyor

${ displaystyle A (u) = f '(u) = { frac { kısmi f} { kısmi u}}}$

Doğrusal olmayan durum

Doğrusal olmayan genel bir denklem için Lax-Wendroff'un muhafazakar formu şu şekildedir:

${ displaystyle u_ {i} ^ {n + 1} = u_ {i} ^ {n} - { frac { Delta t} {2 Delta x}} sol [f (u_ {i + 1} ^ {n}) - f (u_ {i-1} ^ {n}) sağ] + { frac { Delta t ^ {2}} {2 Delta x ^ {2}}} sol [A_ { i + 1/2} left (f (u_ {i + 1} ^ {n}) - f (u_ {i} ^ {n}) sağ) -A_ {i-1/2} left (f (u_ {i} ^ {n}) - f (u_ {i-1} ^ {n}) sağ) sağ].}$

nerede ${ displaystyle A_ {i pm 1/2}}$ Jacobian matrisi değerlendirilir mi ${ displaystyle { frac {1} {2}} (u_ {i} ^ {n} + u_ {i pm 1} ^ {n})}$.

Jacobian içermeyen yöntemler

Jacobian değerlendirmesinden kaçınmak için iki aşamalı bir prosedür kullanın.

Richtmyer yöntemi

Aşağıda Richtmyer iki adımlı Lax – Wendroff yöntemi yer almaktadır. Richtmyer iki adımlı Lax – Wendroff yöntemindeki ilk adım, f (u (x, t)) yarı zamanlı adımlarda, t_n + 1/2 ve yarım ızgara noktaları, x_ben + 1/2. İkinci adımda değerler t_n + 1 verileri kullanılarak hesaplanır t_n ve t_n + 1/2.

İlk (Lax) adımlar:

${ displaystyle u_ {i + 1/2} ^ {n + 1/2} = { frac {1} {2}} (u_ {i + 1} ^ {n} + u_ {i} ^ {n} ) - { frac { Delta t} {2 , Delta x}} (f (u_ {i + 1} ^ {n}) - f (u_ {i} ^ {n})),}$

${ displaystyle u_ {i-1/2} ^ {n + 1/2} = { frac {1} {2}} (u_ {i} ^ {n} + u_ {i-1} ^ {n} ) - { frac { Delta t} {2 , Delta x}} (f (u_ {i} ^ {n}) - f (u_ {i-1} ^ {n})).}$

İkinci adım:

${ displaystyle u_ {i} ^ {n + 1} = u_ {i} ^ {n} - { frac { Delta t} { Delta x}} sol [f (u_ {i + 1/2} ^ {n + 1/2}) - f (u_ {i-1/2} ^ {n + 1/2}) sağ].}$

MacCormack yöntemi

Aynı türden başka bir yöntem de MacCormack tarafından önerildi. MacCormack'in yöntemi önce ileri farklılaştırmayı ve ardından geriye doğru farklılaşmayı kullanır:

İlk adım:

${ displaystyle u_ {i} ^ {*} = u_ {i} ^ {n} - { frac { Delta t} { Delta x}} (f (u_ {i + 1} ^ {n}) - f (u_ {i} ^ {n})).}$

İkinci adım:

${ displaystyle u_ {i} ^ {n + 1} = { frac {1} {2}} (u_ {i} ^ {n} + u_ {i} ^ {*}) - { frac { Delta t} {2 Delta x}} left [f (u_ {i} ^ {*}) - f (u_ {i-1} ^ {*}) sağ].}$

Alternatif olarak, İlk adım:

${ displaystyle u_ {i} ^ {*} = u_ {i} ^ {n} - { frac { Delta t} { Delta x}} (f (u_ {i} ^ {n}) - f ( u_ {i-1} ^ {n})).}$

İkinci adım:

${ displaystyle u_ {i} ^ {n + 1} = { frac {1} {2}} (u_ {i} ^ {n} + u_ {i} ^ {*}) - { frac { Delta t} {2 Delta x}} left [f (u_ {i + 1} ^ {*}) - f (u_ {i} ^ {*}) sağ].}$

Referanslar

• P.D Lax; B. Wendroff (1960). "Koruma yasaları sistemleri". Commun. Pure Appl. Matematik. 13 (2): 217–237. doi:10.1002 / cpa.3160130205.
• Michael J. Thompson, Astrofiziksel Akışkanlar Dinamiğine Giriş, Imperial College Press, Londra, 2006.
• Basın, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). "Bölüm 20.1. Akışı Muhafazakar Başlangıç ​​Değer Sorunları". Sayısal Tarifler: Bilimsel Hesaplama Sanatı (3. baskı). New York: Cambridge University Press. s. 1040. ISBN  978-0-521-88068-8.