求解交通流LWR模型方程的高效數(shù)值方法研究

馮娟娟　楊苗苗

摘 要 對(duì)于交通流LWR模型的數(shù)值求解，單元交界面處采用五階WENO-Z+重構(gòu)，時(shí)間方向的推進(jìn)采用強(qiáng)穩(wěn)定的三步三階Runge-Kutta方法。從而得到了一種高精度、高分辨率以及數(shù)值穩(wěn)定的WENO-Z+格式。結(jié)果顯示，對(duì)激波有良好的捕捉效果，在間斷區(qū)域沒(méi)有非物理振蕩，是模擬交通流LWR模型的理想方法。

關(guān)鍵詞 交通流 高精度 高分辨率

中圖分類號(hào)：TP301.6 文獻(xiàn)標(biāo)識(shí)碼：A

1概述

對(duì)于交通流理論的研究始于上世紀(jì)三十年代，最初是應(yīng)用概率論方法分析交通流量和車速的關(guān)系，1933年Kinzer首次將Poisson分布應(yīng)用于交通流的分析。在Greenshields的研究基礎(chǔ)上，Lighthill和Whitham與Richards各自獨(dú)立提出了交通流連續(xù)介質(zhì)理論（trafficflowcontinuumtheory），即現(xiàn)在通常所指的LWR模型。LWR模型有易于實(shí)現(xiàn)計(jì)算的特點(diǎn)，同時(shí)從實(shí)際的應(yīng)用來(lái)看，對(duì)于交通管理和信號(hào)燈控制有良好的模擬效果。因此LWR模型具有重要的理論和實(shí)際意義，本文以該模型為研究對(duì)象。

2 LWR模型

交通流LWR模型的質(zhì)量守恒微分形式以及速度公式如下：

其中

為時(shí)間，為距離，為車流密度，為最大限速，為阻塞密度。交通流LWR模型的偏微分方程為雙曲形式，具有明確的質(zhì)量守恒意義。

3五階WENO-Z+重構(gòu)

WENO方法，是近年來(lái)廣泛使用的一種高分辨率數(shù)值計(jì)算方法，主要用于雙曲守恒律方程的數(shù)值求解。WENO方法是在ENO格式的基礎(chǔ)上加權(quán)構(gòu)造的，其性能更穩(wěn)定，對(duì)定常問(wèn)題收斂性更好，它能保證在解的光滑區(qū)精度更高，在解的間斷區(qū)保持陡峭的間斷過(guò)渡和本質(zhì)無(wú)振蕩的優(yōu)良性質(zhì)。然而，WENO格式在極值點(diǎn)處精度會(huì)降低（例如五階WENO格式在極值點(diǎn)處會(huì)降為三階）。所以，為了提高WENO格式在極值點(diǎn)處的精度，Borges等通過(guò)構(gòu)造全局高階光滑因子的方式，提出五階WENO-Z格式。通過(guò)增大非光滑模板所對(duì)應(yīng)的非線性權(quán)重的方式構(gòu)造了WENO-Z+格式，該格式具有比五階WENO-Z格式更低的耗散性。此格式在解的光滑區(qū)域有更好的分辨率，在激波和間斷區(qū)域可以保持之前的穩(wěn)定性。

它的新奇之處是在經(jīng)典WENO格式下添加了全局光滑因子，這樣可以提高精度，全局光滑因子是之前光滑因子的一個(gè)線性組合。例如，一般情況下，對(duì)于五階WENO-Z格式，全局光滑因子為。

為了減少WENO-Z的耗散，在非光滑模板上增加一個(gè)可控的量，理想條件下，它需要滿足的要求是：在非光滑模板上，其數(shù)值隨著的增大而增大以提高分辨率，但是如果在模版內(nèi)解出現(xiàn)間斷時(shí)，為防止振蕩的出現(xiàn)以及數(shù)值不穩(wěn)定，其數(shù)值隨著的增大而減小。定義。其中是非光滑模板權(quán)重的一個(gè)細(xì)微增量，實(shí)驗(yàn)結(jié)果表明取時(shí)，效果最好。因此得到了WENO-Z+的非歸一化非線性權(quán)重為：

在時(shí)間方向的推進(jìn)采用三階的Runge-Kutta方法：

參考文獻(xiàn)

[1] Acker，F(xiàn).&R.; B. D. R. Borges&B.; Costa.An improved WENO-Z scheme[J]. J. Comput. Phys.，2016（313）：726-753.