Benjamin-Aaron-Mayumbo (BAM) modeli, trafik yoğunluğu ve ortalama hız gibi toplu değişkenlerle çalışan, klasik Payne-Whitham (PW) modelini öndeki araçla kullanılabilir mesafeye bağlı sürücü tepkisiyle genişleten ikinci dereceden bir makroskopik trafik akış modelidir. Modelin temel yeniliği, ileri yönlü araç aralığını yalnız istenen hızı belirlemek için değil, sürücünün mevcut hızdan istenen hıza ne kadar hızlı yaklaşacağını tanımlayan etkin gevşeme süresini de değiştirmek için kullanmasıdır. Çalışma BAM'ı analitik kararlılık incelemesi ve kontrollü sayısal simülasyonlarla klasik PW modeliyle karşılaştırmaktadır.

Model makroskopik yoğunluktan yaklaşık araç aralığına geçmek için \(\Delta x \approx 1/\rho\) ilişkisini kullanır. Bu yaklaşım, trafik akımındaki araçların yaklaşık homojen aralıklarla dağıldığı varsayımına dayanır. Elde edilen aralık daha sonra iki ayrı mekanizmaya girer: sürücünün ulaşmayı hedeflediği boşluk kontrollü hız \(v_d\) ve hız ayarlamasının zaman ölçeğini belirleyen etkin gevşeme süresi \(\tau_{eff}\). Böylece trafik yoğunluğu arttığında kullanılabilir ileri mesafe küçülür; model hem hedef hızı hem de hız ayarlama davranışını aynı trafik durumuna bağlar.

Çalışmanın temel karşılaştırmalarında klasik PW modeli test edilen koşullarda çok büyük sayısal hız salınımları üretirken, BAM hızları model için öngörülen aralık içinde tutmuş ve daha düzgün trafik dalgaları oluşturmuştur. Kaynak, PW simülasyonunda hızın yaklaşık \(320\,m/s\)'ye kadar yükselip yaklaşık \(-180\,m/s\)'ye kadar düşebildiğini bildirirken BAM'ın 0–20 \(m/s\) olarak belirlenen modelleme aralığının içinde kaldığını belirtmektedir. Bu, BAM'ın incelenen sayısal düzenekte daha kontrollü davranması anlamına gelir; gerçek trafikte daha doğru tahmin yaptığı anlamına gelmez.

Doğrusal kararlılık analizinde çalışmanın türettiği yapısal koşul \(C_0 \geq -\rho_0 v_d'(\rho_0)\) biçimindedir. Boşluk kontrollü istenen hız fonksiyonunun türevi yerine konulduğunda sınır \(C_0 \geq \frac{s\rho_0v_{max}}{(1+s\rho_0)^2}\) olur. \(\rho_0=0.30\,veh/m\), \(v_{max}=20\,m/s\) ve \(s=8\,m\) için çalışma kritik değeri yaklaşık \(C_{0,critical}=4.15\,m/s\) olarak hesaplamıştır. \(C_0=2.0\,m/s\) kararsız, \(4.15\,m/s\) kritik sınır ve \(20.0\,m/s\) kararlı durum olarak sayısal biçimde sınanmıştır.

En önemli yorum sınırı şudur: BAM henüz gerçek trafik gözlemleriyle kalibre edilmemiş veya bağımsız biçimde doğrulanmamıştır. Sonuçlar kapalı çevrim yollar, tanımlanmış başlangıç pertürbasyonları ve seçilmiş parametreler altında gerçekleştirilen kontrollü nümerik deneylerin sonuçlarıdır. Dolayısıyla çalışma BAM'ın gerçek yol ağlarında PW'den kesin olarak daha doğru tahmin yaptığını göstermemekte; daha kontrollü ve kararlı sayısal davranış üretebildiğini göstermektedir.

BAM Trafik Akış Modeli Nedir?

Benjamin-Aaron-Mayumbo (BAM) modeli, klasik Payne-Whitham modelinin araç korunumu ve trafik-basıncı/öngörü yapısını koruyan, ancak sabit gevşeme davranışı yerine trafik yoğunluğundan türetilen ileri araç aralığına bağlı istenen hız ve etkin gevşeme süresi kullanan ikinci dereceden makroskopik bir trafik akış modelidir.

Makroskopik trafik modeli neyi temsil eder?

Mikroskopik trafik modelleri tek tek araçları ve sürücüleri takip ederken, makroskopik modeller trafiği akan sürekli bir ortam gibi ele alır. Bu nedenle her aracın direksiyon, gaz veya fren komutunu ayrı ayrı modellemek yerine yol boyunca yoğunluk, ortalama hız ve akış gibi alan değişkenleri kullanılır.

Bu çalışmada temel değişkenler şunlardır:

LWR modeli: araçların korunumu

Çalışmanın matematiksel başlangıç noktalarından biri Lighthill-Whitham-Richards (LWR) modelidir. Kaynaktaki Denklem (1):

\[ \frac{\partial \rho}{\partial t}+\frac{\partial(\rho v)}{\partial x}=0 \tag{1} \]

Bu denklem bir korunum yasasıdır. \(\rho v\) çarpımı trafik akışını temsil eder. Belirli bir yol kesimindeki yoğunluk zamanla değişiyorsa bunun nedeni araçların kesime girip çıkmasıdır. Denklem kendi başına sürücünün nasıl hızlandığını veya yavaşladığını söylemez; yalnız araç sayısının korunmasını ifade eder.

Payne-Whitham modeli neden ikinci bir denkleme ihtiyaç duyar?

LWR modeli hızın dinamik değişimini ayrı bir denklemle açıklamaz. Payne-Whitham modeli bu nedenle korunum denklemine bir hız/momentum denklemi ekler:

\[ \frac{\partial \rho}{\partial t}+\frac{\partial(\rho v)}{\partial x}=0 \tag{2} \]

\[ \frac{\partial v}{\partial t} +v\frac{\partial v}{\partial x} +\frac{C_0^2}{\rho}\frac{\partial\rho}{\partial x} = \frac{v_e(\rho)-v}{\tau} \tag{3} \]

Denklem (3)'te \(\frac{\partial v}{\partial t}\) hızın zamana göre değişimini, \(v\frac{\partial v}{\partial x}\) akımla taşınan hız değişimini ve \(\frac{C_0^2}{\rho}\frac{\partial\rho}{\partial x}\) trafik yoğunluğu gradyanının oluşturduğu basınç/öngörü etkisini temsil eder. Sağ taraftaki \(\frac{v_e(\rho)-v}{\tau}\) ise sürücünün mevcut hızını denge hızına doğru ayarlamasını ifade eden gevşeme terimidir.

Klasik PW modelinde sorun, \(\tau\)'nun sabit alınmasıdır. BAM'ın temel müdahalesi tam olarak bu kaynak terimine yöneliktir.

Boşluk kontrollü reaksiyon terimi

Çalışmanın Denklem (4)'ü:

\[ A_{gap} = \frac{ v_{max}\frac{\Delta x}{\Delta x+s}-v }{ \tau_0\frac{\Delta x}{\Delta x+s_0} } \tag{4} \]

Payın ilk kısmı sürücünün mevcut boşluk altında ulaşmak istediği hızı, payda ise hızın bu değere hangi zaman ölçeğinde yaklaşacağını temsil eder. Böylece BAM yalnız “ne kadar hızlı gitmek istiyorum?” sorusunu değil, “o hıza ne kadar hızlı yaklaşacağım?” sorusunu da araçlar arasındaki boşluğa bağlar.

Greenshields hız-yoğunluk ilişkisi

Kaynak, klasik denge ilişkisi için Greenshields modelini kullanır:

\[ v_e(\rho)=v_{max}\left(1-\frac{\rho}{\rho_m}\right) \tag{5} \]

Yoğunluk sıfıra yaklaştığında denge hızı \(v_{max}\)'a yaklaşır. Yoğunluk \(\rho_m\)'ye yaklaştığında denklemde denge hızı sıfıra yaklaşır.

Denklem (6), bu ilişkinin yoğunluğa göre türevidir:

\[ v_e'(\rho) = \frac{d}{d\rho} \left[ v_{max}\left(1-\frac{\rho}{\rho_m}\right) \right] = -\frac{v_{max}}{\rho_m} \tag{6} \]

Negatif işaret, yoğunluk yükseldikçe denge hızının azaldığını matematiksel olarak gösterir.

Etkin Gevşeme Süresi Gerçek Sürücü Reaksiyon Gecikmesi midir?

Hayır. BAM'daki \(\tau_{eff}\), sürücünün bir uyarıyı gördükten sonra belirli bir süre beklemesini temsil eden açık bir zaman gecikmesi değildir; mevcut hızın boşluk kontrollü istenen hıza yaklaşma hızını belirleyen, trafik durumuna bağlı bir zaman ölçeğidir. Kaynak bu ayrımı özellikle vurgular ve gerçek bir gecikme modelinde sistemin \(t-\tau\) gibi geçmiş zamandaki trafik durumuna bağlı olması gerektiğini belirtir.

Boşluk kontrollü istenen hız

Denklem (7):

\[ v_d = v_{max}\frac{\Delta x}{\Delta x+s} \tag{7} \]

\(\Delta x\) büyüdükçe \(\frac{\Delta x}{\Delta x+s}\) oranı 1'e yaklaşır ve \(v_d\) serbest akış hızı \(v_{max}\)'a yaklaşır. Araçlar birbirine yaklaştıkça oran küçülür ve istenen hız düşer.

Boşluğa bağlı etkin gevşeme süresi

Denklem (8):

\[ \tau_{eff} = \tau_0\frac{\Delta x}{\Delta x+s_0} \tag{8} \]

Bu formül BAM'ın klasik sabit \(\tau\) yaklaşımından ayrıldığı ikinci önemli noktadır. \(\tau_{eff}\) trafik durumuyla değişir.

Kaynakta \(s_0\) “spacing scaling parameter” olarak tanımlanmakta ve Tablo 5'te 0.30 değeri verilmektedir. Ancak bu parametrenin birimi tabloda açık biçimde raporlanmamıştır. Denklem yapısında \(\Delta x+s_0\) toplamında yer aldığı için boyutsal yorum açısından bu eksiklik önemlidir; Verianla bu boşluğa kaynakta olmayan bir birim eklememektedir.

İvme denklemi

Denklem (9):

\[ a=\frac{v_d-v}{\tau_{eff}} \tag{9} \]

İstenen hız mevcut hızdan büyükse pay pozitiftir ve model hızlanma yönünde; mevcut hız istenen hızdan yüksekse pay negatiftir ve yavaşlama yönünde tepki üretir. \(\tau_{eff}\) ise bu uyarlamanın ne kadar hızlı gerçekleşeceğini ölçekler.

Denklem (7) ve Denklem (8)'deki ifadeler Denklem (9)'da birleştirildiğinde kaynakta Denklem (10) olarak verilen boşluk kontrollü ivme yapısı elde edilir:

\[ a= \frac{ v_{max}\frac{\Delta x}{\Delta x+s}-v }{ \tau_0\frac{\Delta x}{\Delta x+s_0} } \tag{10} \]

Kaynak içi numaralandırma notu: PDF açıklama metni bu türetimin bazı yerlerinde denklem numaralarına hatalı çapraz referanslar vermektedir. Burada formüller kaynakta basılan matematiksel ifadeye göre sıralanmıştır; çapraz referans hataları sessizce kaynak metninin parçasıymış gibi yeniden üretilmemiştir.

Mikroskopik aralıktan makroskopik yoğunluğa geçiş

Çalışmanın kilit yaklaşımı:

\[ \Delta x\approx\frac{1}{\rho} \tag{11} \]

Bu ilişki, araçların yaklaşık homojen biçimde dağıldığı bir trafik akımında yoğunluğun tersini ortalama ileri mesafe için makroskopik bir yaklaşık olarak kullanır. \(\rho\) veh/m olarak ifade edildiğinde \(1/\rho\), araç başına düşen uzunluk ölçeğini verir. Kaynak da bu dönüşümün nispeten homojen araç aralığı varsayımına dayandığını açıkça belirtmektedir.

Denklem (11), istenen hız denklemine uygulanınca:

\[ v_d = v_{max} \frac{1/\rho}{1/\rho+s} = \frac{v_{max}}{1+s\rho} \tag{12} \]

Benzer şekilde etkin gevşeme süresi:

\[ \tau_{eff} = \tau_0 \frac{1/\rho}{1/\rho+s_0} = \frac{\tau_0}{1+s_0\rho} \tag{13} \]

Böylece BAM'ın ivme ifadesi doğrudan trafik yoğunluğuyla yazılabilir:

\[ a= \frac{ \frac{v_{max}}{1+s\rho}-v }{ \frac{\tau_0}{1+s_0\rho} } \tag{14} \]

Yeniden düzenlenmiş biçimi:

\[ a= \frac{1+s_0\rho}{\tau_0} \left( \frac{v_{max}}{1+s\rho}-v \right) \tag{15} \]

Denklem (15) BAM'ın özünü çok açık biçimde gösterir. Yoğunluk \(\rho\), hem hedef hız \(\frac{v_{max}}{1+s\rho}\) hem de uyarlama ölçeği \(\frac{1+s_0\rho}{\tau_0}\) üzerinde etkilidir.

Nihai BAM hız denklemi

Klasik PW gevşeme terimi yeni yoğunluk-bağımlı kaynak terimiyle değiştirildiğinde:

\[ \frac{\partial v}{\partial t} + v\frac{\partial v}{\partial x} + \frac{C_0^2}{\rho} \frac{\partial\rho}{\partial x} = \frac{1+s_0\rho}{\tau_0} \left( \frac{v_{max}}{1+s\rho}-v \right) \tag{16} \]

Bu denklem, araç korunumu denklemiyle birlikte BAM sistemini oluşturur:

\[ \frac{\partial\rho}{\partial t} + \frac{\partial(\rho v)}{\partial x} = 0 \tag{17} \]

PW ile BAM arasındaki temel matematiksel fark

Roe ayrıştırması ve korunum biçimi

BAM ve PW denklemleri sayısal çözüm için Roe yaklaşık Riemann çözücüsüyle sonlu hacim biçiminde ayrıklaştırılmıştır. Kaynak hiperbolik denge yasasını şu genel biçimde verir:

\[ G_t+f(G)_x=S(G) \tag{18} \]

Kaynakta ilk vektör gösterimi:

\[ G= \begin{bmatrix} \rho\\ v \end{bmatrix} \tag{19} \]

\[ f(G)= \begin{bmatrix} \rho v\\ \frac{v^2}{2}+C_0^2\ln(\rho) \end{bmatrix} \tag{20} \]

\[ S(G)= \begin{bmatrix} 0\\ \frac{1+s_0\rho}{\tau_0} \left( \frac{v_{max}}{1+s\rho}-v \right) \end{bmatrix} \tag{21} \]

Burada \(G\) trafik durumunu, \(f(G)\) bu durumun uzay boyunca taşınmasını ve \(S(G)\) sürücü uyarlaması gibi kaynak etkilerini ifade eder.

Kaynakta Denklem (22)'nin ilk türev terimi baskıda şu biçimde görünmektedir:

\[ \frac{\partial t}{\partial G} + A(G)\frac{\partial G}{\partial x} = 0 \tag{22, kaynakta basılı biçim} \]

Bu ifade, hemen çevresindeki açıklamada “quasilinear form” olarak tanımlanmaktadır. Standart hiperbolik sistem notasyonuyla karşılaştırıldığında ilk türev yönünün baskıda ters yazılmış görünmesi bir kaynak içi dizgi/türev gösterimi tutarsızlığıdır. Verianla bunu sessizce kaynak formülüymüş gibi düzeltmemekte, tutarsızlığı açıkça işaretlemektedir.

Korunumlu momentum biçimine geçiş: Denklem (23)–(36)

Kaynak, Roe çözücüsünde kullanılabilecek korunumlu yapıya ulaşmak için cebirsel dönüşüm zinciri kurmaktadır. Denklem (23):

\[ v\rho_t+v(\rho v)_x=0 \tag{23} \]

Denklem (24):

\[ v\rho_t=(\rho v)_t-\rho v_t \tag{24} \]

Bunlardan Denklem (25) yazılır:

\[ \rho v_t=(\rho v)_t+v(\rho v)_x \tag{25} \]

PW hız denklemi yoğunlukla çarpıldığında kaynak Denklem (26)'yı verir:

\[ \rho v_t+\rho vv_x+C_0^2\rho_x = \rho\frac{v_e(\rho)-v}{\tau} \tag{26} \]

Uzaysal taşıma terimi için Denklem (27):

\[ \rho vv_x=(\rho v^2)_x-v(\rho v)_x \tag{27} \]

Bu ilişkiler birleştirildiğinde:

\[ (\rho v)_t+(\rho v^2)_x+C_0^2\rho_x = \rho\frac{v_e(\rho)-v}{\tau} \tag{28} \]

Kaynak, momentum akısını korunan değişken \(\rho v\) ile yeniden ifade eden ara dönüşümü Denklem (29) olarak verir ve ardından klasik PW'nin korunumlu momentum biçimine ulaşır:

\[ (\rho v)_t+ \left( \frac{(\rho v)^2}{\rho}+C_0^2\rho \right)_x = \rho\frac{v_e(\rho)-v}{\tau} \tag{30} \]

Kaynağın Denklem (31)'inde durum vektörü gösterimi ile hemen altındaki “korunan değişkenler yoğunluk ve momentumdur” açıklaması arasında notasyon uyumsuzluğu görülmektedir. Daha sonraki Denklem (36), korunan durum vektörünü açık biçimde \([\rho,\rho v]^T\) olarak vermektedir. Bu nedenle burada kaynak içindeki iki gösterim birbirine karıştırılmamaktadır.

BAM için momentum-benzeri hız denklemi tekrar:

\[ v_t+vv_x+\frac{C_0^2}{\rho}\rho_x = \frac{1+s_0\rho}{\tau_0} \left( \frac{v_{max}}{1+s\rho}-v \right) \tag{32} \]

Yoğunlukla çarpılmış biçimi:

\[ \rho v_t+\rho vv_x+C_0^2\rho_x = \rho\frac{1+s_0\rho}{\tau_0} \left( \frac{v_{max}}{1+s\rho}-v \right) \tag{33} \]

Çarpım kuralı ve araç korunum denklemi kullanıldıktan sonra:

\[ (\rho v)_t+ (\rho v^2+C_0^2\rho)_x = \rho\frac{1+s_0\rho}{\tau_0} \left( \frac{v_{max}}{1+s\rho}-v \right) \tag{34} \]

Hız \(v=(\rho v)/\rho\) kullanılarak tamamen korunan değişkenlerle:

\[ (\rho v)_t+ \left( \frac{(\rho v)^2}{\rho}+C_0^2\rho \right)_x = \rho\frac{1+s_0\rho}{\tau_0} \left( \frac{v_{max}}{1+s\rho}-v \right) \tag{35} \]

Kaynağın nihai korunumlu vektör gösterimi:

\[ G= \begin{bmatrix} \rho\\ \rho v \end{bmatrix}, \qquad f(G)= \begin{bmatrix} \rho v\\ \frac{(\rho v)^2}{\rho}+C_0^2\rho \end{bmatrix} \]

\[ S(G)= \begin{bmatrix} 0\\ \rho\frac{1+s_0\rho}{\tau_0} \left( \frac{v_{max}}{1+s\rho}-v \right) \end{bmatrix} \tag{36} \]

Jacobian ve karakteristik dalga hızları

PW akı yapısına ait Jacobian:

\[ A(G)= \begin{bmatrix} 0 & 1\\ C_0^2-v^2 & 2v \end{bmatrix} \tag{37} \]

Özdeğer problemi:

\[ A-\lambda I= \begin{bmatrix} -\lambda & 1\\ C_0^2-v^2 & 2v-\lambda \end{bmatrix} \tag{38} \]

Karakteristik özdeğerler:

\[ \lambda_1=v+C_0, \qquad \lambda_2=v-C_0 \tag{39} \]

BAM kaynak terimini değiştirip temel akı yapısını değiştirmediği için çalışmada PW ve BAM'ın aynı karakteristik dalga hızlarına sahip olduğu sonucuna varılmaktadır. Özdeğerler gerçek ve ayrı olduğunda sistem sıkı hiperbolik yapısını korur. Kaynak ayrıca özdeğerlerin sıfıra yaklaştığı sonik noktalarda Roe çözücüsünün fiziksel olmayan genişleme şoklarını önlemek amacıyla Harten entropy fix kullandığını belirtmektedir.

BAM Modelinin Kararlılık Koşulu Nedir?

Çalışmanın doğrusal pertürbasyon analizinde BAM'ın küçük bozulmaları asimptotik olarak sönümleyebilmesi için elde edilen yapısal koşul \(C_0\geq-\rho_0v_d'(\rho_0)\)'dır; modelin boşluk kontrollü istenen hız türevi yerine konulduğunda bu sınır \(C_0\geq\frac{s\rho_0v_{max}}{(1+s\rho_0)^2}\) biçimine dönüşür.

Küçük pertürbasyonların tanımlanması

Düzgün denge durumu çevresinde:

\[ \rho(x,t)=\rho_0+\delta\rho(x,t), \qquad v(x,t)=v_0+\delta v(x,t) \tag{40} \]

Burada \(\delta\rho\ll\rho_0\) ve \(\delta v\ll v_0\) kabul edilir. Yani doğrusal analiz büyük trafik değişimlerini değil, denge etrafındaki küçük bozulmaları inceler.

Doğrusallaştırılmış korunum denklemi:

\[ \frac{\partial(\delta\rho)}{\partial t} + v_0\frac{\partial(\delta\rho)}{\partial x} + \rho_0\frac{\partial(\delta v)}{\partial x} = 0 \tag{41} \]

Doğrusallaştırılmış hız denklemi:

\[ \frac{\partial(\delta v)}{\partial t} + v_0\frac{\partial(\delta v)}{\partial x} + \frac{C_0^2}{\rho_0} \frac{\partial(\delta\rho)}{\partial x} = \frac{1+s_0\rho_0}{\tau_0} \left[ v_d'(\rho_0)\delta\rho-\delta v \right] \tag{42} \]

Denge durumundaki etkin gevşeme süresi:

\[ \tau_{eff}=\frac{\tau_0}{1+s_0\rho_0} \]

Kaynak, değişkenler elimine edildiğinde yoğunluk pertürbasyonu için ikinci dereceden kısmi diferansiyel denklemi şu biçimde verir:

\[ \left( \frac{\partial}{\partial t} + v_0\frac{\partial}{\partial x} \right)^2 \delta\rho - C_0^2 \frac{\partial^2(\delta\rho)}{\partial x^2} = -\frac{1}{\tau_{eff}} \left( \frac{\partial}{\partial t} + v_0\frac{\partial}{\partial x} \right) \delta\rho - \frac{\rho_0v_d'(\rho_0)}{\tau_{eff}} \frac{\partial(\delta\rho)}{\partial x} \tag{43} \]

Dalga çözümü ve dispersiyon ilişkisi

Kaynak, pertürbasyonu harmonik dalga biçiminde ele alır:

\[ \delta\rho(x,t)=\hat{\rho}e^{i(kx-\omega t)} \]

Burada \(k\) uzaysal dalga sayısı, \(\omega=\omega_r+i\omega_i\) karmaşık açısal frekanstır. Bu varsayım Denklem (43)'e uygulandığında kaynak Denklem (44)'ü verir:

\[ (v_0k-\omega)^2 - C_0^2k^2 - \frac{i}{\tau_{eff}}(v_0k-\omega) - \frac{i\rho_0v_d'(\rho_0)k}{\tau_{eff}} = 0 \tag{44} \]

Uzun dalga boylu, yani küçük \(k\) bozulmaları için açılım yapıldığında kinematik dalga yayılma hızı:

\[ c_k=v_0+\rho_0v_d'(\rho_0) \tag{45} \]

İkinci derece terimler büyüme/sönüm parametresini verir:

\[ \alpha= \tau_{eff} \left[ C_0^2- \left( \rho_0v_d'(\rho_0) \right)^2 \right] \tag{46} \]

Kaynağın elde ettiği doğrusal kararlılık koşulu:

\[ C_0\geq-\rho_0v_d'(\rho_0) \tag{47} \]

Boşluk kontrollü istenen hızın yoğunluk türevi:

\[ v_d'(\rho_0) = -\frac{s\,v_{max}}{(1+s\rho_0)^2} \]

Böylece nihai kararlılık sınırı:

\[ C_0 \geq \frac{s\rho_0v_{max}} {(1+s\rho_0)^2} \tag{48} \]

Kararlılık sınırının sayısal doğrulanması

Kaynak \(\rho_0=0.30\,veh/m\), \(v_{max}=20\,m/s\) ve \(s=8\,m\) için:

\[ C_{0,critical} = \frac{8\times0.30\times20} {(1+8\times0.30)^2} \approx 4.15\,m/s \]

Kaynak Şekil 1'de üç durumun yol konumu boyunca yoğunluk profilleri farklı zamanlarda gösterilmektedir. Kararlı \(C_0=20\,m/s\) durumunda başlangıç sinüzoidal dalgasının genliği zamanla küçülürken, \(C_0=2\,m/s\) durumunda bozulma belirgin biçimde büyür. Grafik burada birebir yeniden kullanılmamıştır; bilimsel sonucu yukarıdaki kaynak değerleriyle bağımsız biçimde yapılandırılmıştır.

Çalışmanın Yöntemi ve Bulguları

Sayısal çözüm yöntemi

PW ve BAM sistemleri Roe approximate Riemann solver kullanılarak sonlu hacim yöntemiyle ayrıklaştırılmıştır. Roe yaklaşımının çalışmadaki rolü, trafik akışında oluşabilen keskin değişimleri ve şok benzeri dalgaları hiperbolik korunum yasaları çerçevesinde sayısal olarak izlemektir. BAM'ın akı Jacobian'ı PW ile aynı kaldığından iki model için aynı temel Roe ayrıştırması kullanılabilmektedir.

Kaynak, sonik noktalarda özdeğerlerden birinin sıfıra yaklaşması hâlinde Harten entropy fix uygulandığını belirtir. Bu düzeltmenin amacı sayısal çözücünün fiziksel olmayan genişleme şokları üretmesini engellemektir.

Temel simülasyon parametreleri

Kaynak içi zaman adımı tutarsızlığı: Tablo 5'te \(\Delta t=0.005\,s\) yazmaktadır ve \(1200\times0.005=6\,s\) olduğundan tablo kendi içinde toplam süreyle uyumludur. Ancak 18. sayfadaki açıklama metninde \(\Delta t=0.006\,s\) yazmaktadır. Verianla bu iki değerden birini araştırmacının yerine sessizce seçmemekte; tablo ve anlatım arasındaki farkı açık bırakmaktadır. Daha sonraki 50 saniyelik sağlamlık deneylerinde kaynak yeniden \(\Delta t=0.005\,s\) kullanıldığını açıkça belirtmektedir.

Başlangıç yoğunluğu denklemi

Kaynak ayrıca Denklem (49)'da başlangıç yoğunluğunu parça parça şöyle tanımlar:

\[ \rho(x,0)= \begin{cases} 0.01, & x\leq30\\ 0.3, & 30<x\leq100\\ 0.1, & x>100 \end{cases} \tag{49} \]

Kaynak içi başlangıç koşulu notu: Bu parça-parça profil, Tablo 5'te “mean initial density = 0.30 veh/m” olarak verilen ifadeyle doğrudan aynı tanım değildir. Kaynak iki bilgiyi birlikte vermektedir ancak aralarındaki ilişkiyi ayrıntılı biçimde açıklamamaktadır. Bu nedenle Verianla bunları yapay biçimde uzlaştırmamaktadır.

0.052 saniyelik değer ne anlama geliyor?

Kaynak \(\tau_0=0.52\,s\)'yi referans insan tepki değeri olarak kullanırken \(0.052\,s\)'yi tipik insan sürücü reaksiyon süresi olarak sunmamaktadır. Değer, temel değerin bir mertebe azaltılmasıyla oluşturulan kontrollü bir duyarlılık parametresidir. Benzer biçimde 2 m standstill aralığı da normal sürücü davranışının evrensel değeri olarak değil, karşılaştırmalı duyarlılık koşulu olarak kullanılır.

PW ve BAM hız davranışının karşılaştırılması

500 m kapalı çevrim yoldaki temel karşılaştırmada PW modeli, ani trafik değişimleri altında kaynakta “fiziksel olarak imkânsız” olarak nitelendirilen şiddetli sayısal salınımlar göstermiştir. Kaynak yaklaşık \(320\,m/s\) maksimum ve \(-180\,m/s\) minimum hız excursion'ları raporlamaktadır.

BAM aynı test koşullarında hızı belirlenmiş 0–20 m/s model aralığında tutmuş ve çok daha düzgün hız evrimi üretmiştir. Ancak yazarlar bu karşılaştırmayı ampirik doğruluk üstünlüğü olarak yorumlamaz. Bulgular yalnız test edilen sayısal koşullarda BAM'ın daha sınırlı ve fiziksel bakımdan daha makul sayısal hız davranışı ürettiğini gösterir.

Yoğunluk alanı nasıl değişmiştir?

BAM için 500 m çevrim yol üzerindeki 50 saniyelik yoğunluk simülasyonunda kaynak, yaklaşık \(x=150\,m\) civarında yaklaşık \(0.45\,veh/m\) yerel maksimum ve yaklaşık \(x=375\,m\) civarında yaklaşık \(0.15\,veh/m\) minimum raporlamaktadır. Alan genel olarak düzgün kalmakta ve belirgin süreksizlikler oluşturmamaktadır.

Klasik PW yoğunluk alanında ise kaynak yaklaşık \(150\,m\) civarında \(0.42\,veh/m\) maksimum ve yaklaşık \(350\,m\) civarında \(0.17\,veh/m\) minimum bildirir. Başlangıçtaki düzgün profil zaman ilerledikçe daha keskin yapılara ve geriye yayılan stop-and-go dalgalarına dönüşmektedir. Çalışmanın yorumu, test edilen koşullarda PW'nin başlangıç bozulmasını büyütmesi, BAM'ın ise daha düzgün dalga yayılımı üretmesidir.

Trafik akışı sonuçları

PW için kaynak, 500 m yol ve 6 saniyelik akış incelemesinde yaklaşık \(150\,m\)'de \(4.55\,veh/s\) maksimum ve yaklaşık \(350\,m\)'de \(3.4\,veh/s\) minimum raporlamaktadır. Dalga biçimi zamanla simetrisini kaybetmekte ve deformasyon göstermektedir.

BAM için akış maksimumu yaklaşık \(125\,m\)'de yine \(4.55\,veh/s\), minimum ise yaklaşık \(375\,m\)'de \(3.4\,veh/s\) olarak raporlanmıştır. Akış büyüklüklerinin birbirine yakın olması önemlidir: yazarların sonuç bölümünde de belirttiği gibi BAM'ın başlıca farkı testlerde “daha fazla trafik akışı” üretmek değil, trafik durumlarının zaman ve uzaydaki gelişiminin daha kontrollü olmasıdır.

Sağlamlık deneyleri nasıl tasarlandı?

Çalışma daha sonra BAM'ı 250, 500 ve 1000 m uzunlukta kapalı çevrim yollar üzerinde dört koşul altında 50 saniye boyunca incelemiştir. Bu deneylerde uzaysal adım 1 m, zaman adımı 0.005 s, \(C_0=20\,m/s\), \(s_0=0.30\) ve \(v_{max}=20\,m/s\) olarak korunmuştur.

D senaryosunda yolun merkezindeki %10'luk kesimde erişilebilir maksimum hız %25 azaltılarak tekrarlanabilir bir yerel kapasite kısıtı oluşturulmuştur.

RMS bozulma oranı neyi gösterir?

Çalışma, başlangıçtaki yoğunluk pertürbasyonunun simülasyon sonunda ne kadar kaldığını değerlendirmek için son RMS yoğunluk pertürbasyonunun başlangıç RMS değerine oranını kullanır. Oran 1'in altında olduğunda bozulma genel olarak sönümlenmiş, 1'in üzerine çıktığında ise son bozulma başlangıçtakinden daha büyük hâle gelmiştir.

Senaryo A'da yol uzunluğu arttıkça son RMS oranı 0.002'den 0.213 ve 0.681'e yükselmiştir. Yani kısa yolda başlangıç bozulması neredeyse tamamen sönümlenirken daha uzun yol daha fazla artık uzaysal değişim taşımıştır.

Senaryo B'de daha yüksek başlangıç yoğunluğu ve daha güçlü pertürbasyona rağmen çözümler test edilen üç yol uzunluğunda da sınırlı kalmıştır. RMS oranları sırasıyla 0.002, 0.219 ve 0.686'dır.

Senaryo C'de daha uzun temel gevşeme zaman ölçeği ve daha büyük standstill aralığı kullanılmıştır. Kaynak metni 250 m için RMS oranını yaklaşık 0.000001 olarak açıklarken Tablo 7 bunu üç ondalıkla 0.000 olarak göstermektedir. Bu fark yuvarlama düzeyindedir; burada tablo değeri tablo içinde aynen korunmuştur.

En güçlü artık bozulma aktif darboğaz Senaryo D'de ortaya çıkmıştır. 500 m ve 1000 m yollar için RMS oranları sırasıyla 1.332 ve 1.923'tür; yani son yoğunluk pertürbasyonu başlangıç değerinden daha büyüktür. Kaynak bunu uygulanan yerel kapasite kısıtının oluşturduğu kalıcı uzaysal yeniden dağılım ve bozulma büyümesi olarak yorumlamaktadır.

Kaynak metninde Senaryo D, 1000 m için maksimum hız 6.728 m/s olarak yazılırken Tablo 7 değeri 6.729 m/s'dir. Verianla sayısal tabloda kaynak tablonun 6.729 m/s değerini korumakta ve anlatı ile tablo arasındaki 0.001 m/s farkı gizlememektedir.

250, 500 ve 1000 metrelik yolların karşılaştırılması

Kaynak Şekil 8–15'te dört senaryonun 250, 500 ve 1000 m yol uzunlukları için son yoğunluk ve hız profilleri verilmektedir. Genel eğilim, kısa yolların başlangıç bozulmasını daha güçlü sönümlemesi; uzun yolların ise daha fazla uzaysal değişimin kalmasına izin vermesidir. Bu sonuç özellikle RMS oranlarında açık biçimde görülmektedir.

Ancak bu durum “kısa yollar gerçek trafikte her zaman daha kararlıdır” biçiminde genellenemez. Çalışmanın yol uzunluğu etkisi kapalı çevrim, periyodik sınır koşulu, seçilen BAM parametreleri ve tanımlanan başlangıç bozulmaları altında elde edilmiş nümerik bir bulgudur.

Reaksiyon zamanı ve standstill aralığı duyarlılığı

Kaynak Şekil 16 ve 17, 500 m çevrim yol üzerinde iki parametre kümesini karşılaştırmaktadır. Birinci durumda \(\tau=0.52\,s\) ve \(l_s=8\,m\); ikinci durumda daha kısa \(0.052\,s\) zaman ölçeği ve 2 m aralık kullanılmıştır. Kaynak, daha büyük zaman ölçeği ve standstill aralığında yoğunluk dalgalarının daha düzgün olduğunu; daha küçük değerlerde ise daha güçlü hız ve yoğunluk değişimleri oluştuğunu bildirmektedir.

0.052 s ve 2 m değerlerinin normal insan sürüşünü temsil ettiği düşünülmemelidir. Makale bu değerleri açıkça kontrollü duyarlılık koşulları olarak kullanmaktadır.

Çalışmanın desteklediği sonuçlar

• BAM'ın test edilen sayısal koşullarda klasik PW'ye göre daha sınırlı hız değerleri ve daha düzgün trafik dalgası evrimi üretebildiği desteklenmektedir.
• Boşluk kontrollü istenen hız ile yoğunluk-bağımlı gevşeme zaman ölçeğinin birlikte kullanılması, modelin kaynak terimini trafik durumuna uyarlamaktadır.
• Doğrusal analiz BAM için açık bir kararlılık sınırı vermektedir.
• \(C_0=2.0\), \(4.15\) ve \(20.0\,m/s\) testleri teorik kararlılık sınırıyla uyumlu sayısal davranış göstermiştir.
• 250–1000 m sağlamlık deneylerinde BAM çözümleri incelenen parametre uzayında sınırlı kalmıştır.
• Aktif darboğaz, incelenen koşullar arasında en büyük artık yoğunluk bozulmasını oluşturmuştur.
• BAM ve PW'nin trafik akış büyüklükleri benzer olabilmektedir; temel fark çalışmada daha çok durumların kararlılığı ve zaman-uzay evrimi olarak ortaya çıkmıştır.

Çalışmanın desteklemediği sonuçlar

• BAM'ın gerçek yol ağlarında PW'den kesin olarak daha doğru olduğu gösterilmemiştir.
• Model gerçek trafik sensörü, loop detector veya araç yörünge verisiyle kalibre edilmemiştir.
• Bağımsız bir gerçek dünya veri setinde tahmin başarısı ölçülmemiştir.
• Farklı araç uzunlukları, araç tipleri ve bireysel sürücü özellikleri açıkça çözülmemektedir.
• Çalışmadaki \(\Delta x\approx1/\rho\) yaklaşımı heterojen araç aralıklarının tamamını temsil etmez.
• 0.052 saniyelik parametre gerçek insan sürücü reaksiyon süresi sonucu değildir.
• 2 m standstill aralığı genel gerçek trafik normu olarak sunulmamaktadır.
• Kontrollü simülasyondaki sayısal kararlılık, gerçek trafik tahmin doğruluğunun kanıtı değildir.

Çalışmanın temel sınırlılıkları

Yazarların açıkça belirttiği en önemli sınırlılık ampirik kalibrasyon ve doğrulamanın bulunmamasıdır. Mevcut sonuçlar kapalı çevrim yol ve önceden tanımlanmış başlangıç pertürbasyonları üzerindeki kontrollü nümerik deneylere dayanmaktadır.

BAM ayrıca ileri araç aralığını \(\Delta x\approx1/\rho\) yaklaşımıyla makroskopik yoğunluktan türetmektedir. Bu yaklaşım nispeten homojen aralık varsayımı taşır ve farklı araç tipleri, araç uzunlukları veya bireysel sürücü özelliklerinden doğan heterojenliği açıkça çözmez.

Gelecek çalışma olarak yazarlar BAM parametrelerini gerçek trafik verisiyle kalibre etmeyi; trafik akışı, hız, yoğunluk, headway ve seyahat süresi gibi bağımsız ölçümlerle doğrulamayı; heterojen trafik koşullarını, parametre duyarlılığını ve daha gerçekçi ağ yapılarını incelemeyi önermektedir.

Kaynak ve Yöntem Notu

Özgün başlık: A macroscopic traffic flow model with gap-controlled driver response and adaptive reaction dynamics

Yazarlar ve sıra: Benjamin Simwinga; Aaron Zimba; Mayumbo Nyirenda.

Sorumlu yazar: Benjamin Simwinga.

Kurumlar: Benjamin Simwinga ve Aaron Zimba — School of Computing, Technology and Applied Sciences, ZCAS University, Lusaka, Zambia. Mayumbo Nyirenda — Department of Computer Science, University of Zambia, Lusaka, Zambia.

Dergi: Discover Networks.

Yayınevi/platform: Springer Nature / Springer Nature Link.

Cilt ve makale numarası: Volume 2, Article 39 (2026).

DOI: 10.1007/s44354-026-00053-y

Yayın tarihi: 8 Ekim 2026.

Süreç tarihleri: Alındı: 10 Ağustos 2026; revize edildi: 19 Eylül 2026; kabul edildi: 28 Eylül 2026.

Kaynak türü: Research article; matematiksel model geliştirme, doğrusal kararlılık analizi ve kontrollü nümerik simülasyon çalışması.

Lisans: Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0). Lisans açık erişim sağlar ancak uyarlanmış materyalin paylaşımına izin vermez; bu nedenle Verianla anlatımı kaynak metin ve figürlerden bağımsız bilimsel bilgi mimarisiyle hazırlanmıştır.

Yazar katkıları: Kaynakta BS'nin makaleyi yazdığı; AZ ve MN'nin gözden geçirme ve yönlendirme sağladığı belirtilmektedir.

Finansman: Çalışma finansman almadığını bildirmektedir.

Veri erişilebilirliği: Çalışma sırasında üretilen ve/veya analiz edilen veri setlerinin makul talep üzerine sorumlu yazardan alınabileceği belirtilmektedir.

Çıkar çatışması: Yazarlar çıkar çatışması olmadığını bildirmektedir.

Etik onay ve katılımcı onamı: Uygulanabilir değil olarak bildirilmiştir.

Ampirik doğrulama durumu: Yoktur. BAM henüz gözlenen trafik verileriyle kalibre edilmemiş ve bağımsız saha ölçümleriyle doğrulanmamıştır.

İstatistiksel çıkarım: Çalışmada p değeri, güven aralığı, klinik/deneysel örneklem istatistiği veya istatistiksel hipotez testi raporlanmamaktadır. Bulgular analitik kararlılık koşulları ve sayısal simülasyon çıktılarıdır.

Kaynak içi teknik tutarsızlıklar: Tablo 5'te \(\Delta t=0.005\,s\) verilmesine karşılık açıklama metninin bir yerinde 0.006 s yazmaktadır. Tablo 5'in 0.005 s değeri 1200 adım ve 6 saniyelik toplam süreyle kendi içinde tutarlıdır; ancak Verianla bu nedenle kaynak metnindeki 0.006 s ifadesini sessizce silmemiştir. Tablo 5'te verilen 0.30 veh/m “ortalama başlangıç yoğunluğu” ile Denklem (49)'daki 0.01/0.30/0.10 parça-parça başlangıç profili arasındaki ilişki de kaynakta ayrıntılı açıklanmamaktadır. Denklem çapraz referanslarında ve bazı vektör gösterimlerinde de dizgi/notasyon uyuşmazlıkları bulunmaktadır.