一維守恒型方程（組）的熵耗散格式.pdf

1、本文考慮的一維雙曲守恒律方程(組)。我們給出了一種設(shè)計(jì)高分辨率格式的方法。此格式是Godunov型的，用的是分片線性重構(gòu)。與傳統(tǒng)的Godunov型格式不同的是此格式在計(jì)算過程中不僅計(jì)算數(shù)值解還計(jì)算了數(shù)值熵。在每個(gè)網(wǎng)格中線性重構(gòu)函數(shù)的斜率是根據(jù)熵耗散得到的，它要求此網(wǎng)格上重構(gòu)函數(shù)熵的網(wǎng)格平均應(yīng)與此網(wǎng)格上的數(shù)值熵相等。數(shù)值解是根據(jù)有限體積法(finite-volume)求得的，同時(shí)數(shù)值熵計(jì)算的時(shí)候格式中有一熵耗散項(xiàng)，它在計(jì)算過程中耗散熵.通

2、過這種方式為數(shù)值計(jì)算引入了適當(dāng)?shù)恼承?，穩(wěn)定了計(jì)算。 在[25]中我們對(duì)于特殊的熵函數(shù)分析了格式的精度。本文中對(duì)于一般的熵函數(shù)，我們分析了格式的精度，格式在遠(yuǎn)離極值點(diǎn)附近為二階精度，在極值點(diǎn)附近為一階精度。因?yàn)楦袷脚c傳統(tǒng)的守恒型格式不同，所以我們給出了它的相容性定義和Lax-Wendroff定理，定理說如果用我們的數(shù)值格式求得的數(shù)值解收斂，則它一定收斂到方程的滿足熵條件的弱解。 設(shè)計(jì)這樣一種格式的一個(gè)重要?jiǎng)訖C(jī)是期望用此類格

3、式來克服守恒型方程組的數(shù)值模擬中，線性特征場(chǎng)上的數(shù)值耗散問題。為此對(duì)于線性傳輸方程，我們用不帶熵耗散函數(shù)的格式進(jìn)行了一些數(shù)值實(shí)驗(yàn)，以研究數(shù)值熵對(duì)消除線性耗散問題的作用。在此研究的基礎(chǔ)上，我們?cè)O(shè)計(jì)了兩種格式，一種格式中我們?nèi)越o線性方程或方程組的線性特征場(chǎng)以一定的熵耗散以穩(wěn)定計(jì)算。另一種格式中我們保持線性方程或方程組的線性特征場(chǎng)上的熵守恒，同時(shí)為數(shù)值解的重構(gòu)設(shè)計(jì)了一種所謂的“極大值減少極小值增”(Minimums-Increase-and-