Modeli ya Benjamin-Aaron-Mayumbo (BAM) ni modeli ya makroskopiki ya mtiririko wa trafiki ya daraja la pili inayofanya kazi kwa vigeu vya jumla kama vile wiani wa trafiki na kasi ya wastani, na inayopanua modeli ya kawaida ya Payne-Whitham (PW) kwa mwitikio wa dereva unaotegemea nafasi inayopatikana hadi gari lililo mbele. Ubunifu mkuu wa modeli ni kutumia nafasi ya mbele kati ya magari si tu kubainisha kasi inayotakiwa, bali pia kubadilisha muda faafu wa ulegezaji unaobainisha jinsi dereva anavyokaribia kwa kasi kutoka kasi ya sasa hadi kasi inayotakiwa. Utafiti unalinganisha BAM na modeli ya kawaida ya PW kwa uchanganuzi wa uthabiti wa kianaliti na uigaji wa kinambari unaodhibitiwa.

Modeli hutumia uhusiano \(\Delta x \approx 1/\rho\) kubadilisha wiani wa makroskopiki kuwa makadirio ya nafasi kati ya magari. Njia hii inategemea dhana kwamba magari katika mtiririko wa trafiki yamesambazwa kwa nafasi takribani zinazofanana. Nafasi inayopatikana kisha huingia katika mifumo miwili tofauti: kasi inayotakiwa inayodhibitiwa na nafasi \(v_d\), ambayo dereva analenga kuifikia, na muda faafu wa ulegezaji \(\tau_{eff}\), unaobainisha kipimo cha muda cha urekebishaji wa kasi. Hivyo, wiani wa trafiki unapoongezeka, umbali unaopatikana mbele hupungua; modeli huunganisha kasi lengwa na tabia ya kurekebisha kasi na hali ileile ya trafiki.

Katika ulinganisho mkuu wa utafiti, modeli ya kawaida ya PW ilizalisha mitetemo mikubwa sana ya kasi ya kinambari katika hali zilizojaribiwa, ilhali BAM iliweka kasi ndani ya masafa yaliyokusudiwa kwa modeli na kuzalisha mawimbi laini zaidi ya trafiki. Chanzo kinaripoti kuwa katika uigaji wa PW kasi inaweza kupanda hadi takribani \(320\,m/s\) na kushuka hadi takribani \(-180\,m/s\), huku BAM ikibaki ndani ya masafa ya uundaji modeli yaliyowekwa kuwa 0–20 \(m/s\). Hii inamaanisha kwamba BAM ilionyesha tabia iliyodhibitiwa zaidi katika usanidi wa kinambari uliochunguzwa; haimaanishi kwamba inafanya utabiri sahihi zaidi katika trafiki halisi.

Katika uchanganuzi wa uthabiti wa mstari, sharti la kimuundo lililotolewa na utafiti lina umbo \(C_0 \geq -\rho_0 v_d'(\rho_0)\). Derivative ya funksi ya kasi inayotakiwa inayodhibitiwa na nafasi inapowekwa, mpaka huwa \(C_0 \geq \frac{s\rho_0v_{max}}{(1+s\rho_0)^2}\). Kwa \(\rho_0=0.30\,veh/m\), \(v_{max}=20\,m/s\) na \(s=8\,m\), utafiti ulihesabu thamani muhimu kuwa takribani \(C_{0,critical}=4.15\,m/s\). \(C_0=2.0\,m/s\) ilijaribiwa kwa kinambari kama hali isiyo thabiti, \(4.15\,m/s\) kama mpaka muhimu, na \(20.0\,m/s\) kama hali thabiti.

Kikomo muhimu zaidi cha tafsiri ni hiki: BAM bado haijakalibishwa kwa uchunguzi wa trafiki halisi wala kuthibitishwa kwa kujitegemea. Matokeo yanatokana na majaribio ya kinambari yaliyodhibitiwa yaliyofanywa kwenye barabara za mzunguko uliofungwa, chini ya usumbufu wa awali uliobainishwa na vigezo vilivyochaguliwa. Kwa hiyo, utafiti hauonyeshi kwamba BAM hufanya utabiri sahihi zaidi kuliko PW katika mitandao halisi ya barabara; unaonyesha kwamba inaweza kuzalisha tabia ya kinambari iliyodhibitiwa na thabiti zaidi.

Modeli ya Mtiririko wa Trafiki ya BAM ni Nini?

Muundo wa Benjamin-Aaron-Mayumbo (BAM) ni modeli ya makroskopiki ya mtiririko wa trafiki ya daraja la pili inayohifadhi muundo wa uhifadhi wa magari na shinikizo la trafiki/utarajiaji wa modeli ya kawaida ya Payne-Whitham, lakini badala ya tabia ya ulegezaji isiyobadilika hutumia kasi inayotakiwa na muda faafu wa ulegezaji unaotegemea nafasi ya mbele kati ya magari inayotolewa kutokana na wiani wa trafiki.

Modeli ya makroskopiki ya trafiki inawakilisha nini?

Wakati modeli za mikroskopiki za trafiki hufuatilia kila gari na dereva mmoja mmoja, modeli za makroskopiki huichukulia trafiki kama mazingira endelevu yanayotiririka. Kwa hiyo, badala ya kuunda modeli tofauti kwa amri za usukani, kuongeza kasi au breki za kila gari, hutumia vigeu vya uga kama wiani, kasi ya wastani na mtiririko kando ya barabara.

Vigeu vikuu katika utafiti huu ni hivi:

Modeli ya LWR: uhifadhi wa magari

Mojawapo ya sehemu za kuanzia za kihisabati za utafiti ni modeli ya Lighthill-Whitham-Richards (LWR). Mlinganyo (1) katika chanzo:

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

Mlinganyo huu ni sheria ya uhifadhi. Zao \(\rho v\) linawakilisha mtiririko wa trafiki. Ikiwa wiani katika sehemu fulani ya barabara unabadilika kwa muda, sababu ni magari kuingia na kutoka katika sehemu hiyo. Mlinganyo wenyewe hausemi jinsi dereva anavyoongeza au kupunguza kasi; unaeleza tu uhifadhi wa idadi ya magari.

Kwa nini modeli ya Payne-Whitham inahitaji mlinganyo wa pili?

Modeli ya LWR haielezi mabadiliko ya mienendo ya kasi kwa mlinganyo tofauti. Kwa hiyo modeli ya Payne-Whitham huongeza mlinganyo wa kasi/momentum kwenye mlinganyo wa uhifadhi:

\[ \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} \]

Katika Mlinganyo (3), \(\frac{\partial v}{\partial t}\) inawakilisha mabadiliko ya kasi kwa muda, \(v\frac{\partial v}{\partial x}\) mabadiliko ya kasi yanayosafirishwa na mtiririko, na \(\frac{C_0^2}{\rho}\frac{\partial\rho}{\partial x}\) athari ya shinikizo/utarajiaji inayosababishwa na gradienti ya wiani wa trafiki. \(\frac{v_e(\rho)-v}{\tau}\) upande wa kulia ni sehemu ya ulegezaji inayoeleza urekebishaji wa kasi ya sasa ya dereva kuelekea kasi ya msawazo.

Tatizo katika modeli ya kawaida ya PW ni kwamba \(\tau\) huchukuliwa kuwa isiyobadilika. Mabadiliko ya msingi ya BAM yanalenga hasa sehemu hii ya chanzo.

Sehemu ya mwitikio unaodhibitiwa na nafasi

Mlinganyo (4) wa utafiti:

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

Sehemu ya kwanza ya namba juu inawakilisha kasi ambayo dereva anataka kufikia chini ya nafasi iliyopo, huku namba chini ikiwakilisha kipimo cha muda ambacho kasi inakaribia thamani hiyo. Hivyo BAM haiunganishi tu swali “ninataka kwenda kwa kasi gani?” na nafasi kati ya magari, bali pia swali “nitakaribia kasi hiyo kwa haraka kiasi gani?”.

Uhusiano wa kasi-wiani wa Greenshields

Chanzo hutumia modeli ya Greenshields kwa uhusiano wa kawaida wa msawazo:

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

Wiani unapokaribia sifuri, kasi ya msawazo hukaribia \(v_{max}\). Wiani unapokaribia \(\rho_m\), kasi ya msawazo katika mlinganyo hukaribia sifuri.

Mlinganyo (6) ni derivative ya uhusiano huu kwa wiani:

\[ 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} \]

Alama hasi inaonyesha kihisabati kwamba kasi ya msawazo hupungua kadiri wiani unavyoongezeka.

Je, Muda Faafu wa Ulegezaji ni Ucheleweshaji Halisi wa Mwitikio wa Dereva?

Hapana. \(\tau_{eff}\) katika BAM si ucheleweshaji wa muda wa moja kwa moja unaowakilisha dereva kusubiri muda fulani baada ya kuona onyo; ni kipimo cha muda kinachotegemea hali ya trafiki na kinachobainisha kasi ambayo kasi ya sasa inakaribia kasi inayotakiwa inayodhibitiwa na nafasi. Chanzo kinasisitiza tofauti hii na kueleza kwamba katika modeli halisi ya ucheleweshaji, mfumo ungepaswa kutegemea hali ya trafiki ya wakati uliopita kama \(t-\tau\).

Kasi inayotakiwa inayodhibitiwa na nafasi

Mlinganyo (7):

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

\(\Delta x\) inapoongezeka, uwiano \(\frac{\Delta x}{\Delta x+s}\) hukaribia 1 na \(v_d\) hukaribia kasi ya mtiririko huru \(v_{max}\). Magari yanapokaribiana, uwiano hupungua na kasi inayotakiwa hushuka.

Muda faafu wa ulegezaji unaotegemea nafasi

Mlinganyo (8):

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

Fomula hii ni sehemu ya pili muhimu ambayo BAM hutofautiana na mkabala wa kawaida wa \(\tau\) isiyobadilika. \(\tau_{eff}\) hubadilika kulingana na hali ya trafiki.

Katika chanzo, \(s_0\) inafafanuliwa kama “spacing scaling parameter” na thamani 0.30 inatolewa katika Jedwali 5. Hata hivyo, kitengo cha kigezo hiki hakijaripotiwa wazi katika jedwali. Kwa kuwa katika muundo wa mlinganyo kinatokea katika jumla \(\Delta x+s_0\), upungufu huu ni muhimu kwa tafsiri ya vipimo; Verianla haiongezi kitengo ambacho hakipo katika chanzo.

Mlinganyo wa uongezaji kasi

Mlinganyo (9):

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

Ikiwa kasi inayotakiwa ni kubwa kuliko kasi ya sasa, namba juu ni chanya na modeli hutoa mwitikio wa kuongeza kasi; ikiwa kasi ya sasa ni kubwa kuliko kasi inayotakiwa, namba juu ni hasi na modeli hutoa mwitikio wa kupunguza kasi. \(\tau_{eff}\) hupima jinsi urekebishaji huu unavyotokea kwa haraka.

Maelezo ya Mlinganyo (7) na Mlinganyo (8) yanapounganishwa katika Mlinganyo (9), muundo wa uongezaji kasi unaodhibitiwa na nafasi unaotolewa katika chanzo kama Mlinganyo (10) hupatikana:

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

Dokezo la uorodheshaji wa milinganyo ndani ya chanzo: Maandishi ya maelezo katika PDF yanatoa marejeo mtambuka yasiyo sahihi kwa namba za milinganyo katika baadhi ya sehemu za utoaji huu. Hapa fomula zimepangwa kulingana na usemi wa kihisabati uliochapishwa katika chanzo; makosa ya marejeo mtambuka hayajarudiwa kimya kimya kana kwamba ni sehemu ya maandishi ya chanzo.

Mpito kutoka nafasi ya mikroskopiki hadi wiani wa makroskopiki

Mkabala muhimu wa utafiti:

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

Uhusiano huu hutumia kinyume cha wiani kama makadirio ya makroskopiki ya nafasi ya wastani ya mbele katika mtiririko wa trafiki ambamo magari yamesambazwa takribani kwa usawa. \(\rho\) inapoelezwa kwa veh/m, \(1/\rho\) hutoa kipimo cha urefu kwa kila gari. Chanzo pia kinaeleza wazi kwamba ubadilishaji huu unategemea dhana ya nafasi kati ya magari kuwa takribani sawa.

Mlinganyo (11) unapotumika katika mlinganyo wa kasi inayotakiwa:

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

Vivyo hivyo, muda faafu wa ulegezaji:

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

Hivyo usemi wa uongezaji kasi wa BAM unaweza kuandikwa moja kwa moja kwa wiani wa trafiki:

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

Umbo lililopangwa upya:

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

Mlinganyo (15) unaonyesha kiini cha BAM kwa uwazi sana. Wiani \(\rho\) huathiri kasi lengwa \(\frac{v_{max}}{1+s\rho}\) na kipimo cha urekebishaji \(\frac{1+s_0\rho}{\tau_0}\).

Mlinganyo wa mwisho wa kasi wa BAM

Sehemu ya ulegezaji ya kawaida ya PW inapobadilishwa na sehemu mpya ya chanzo inayotegemea wiani:

\[ \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} \]

Mlinganyo huu, pamoja na mlinganyo wa uhifadhi wa magari, huunda mfumo wa BAM:

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

Tofauti kuu ya kihisabati kati ya PW na BAM

Utenganishaji wa Roe na umbo la uhifadhi

Milinganyo ya BAM na PW imediscretishwa kwa suluhisho la kinambari katika umbo la ujazo wenye ukomo kwa kutumia Roe approximate Riemann solver. Chanzo hutoa sheria ya mizani ya hyperboliki katika umbo la jumla lifuatalo:

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

Uwakilishi wa kwanza wa vekta katika chanzo:

\[ 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} \]

Hapa \(G\) inawakilisha hali ya trafiki, \(f(G)\) usafirishaji wa hali hiyo katika nafasi, na \(S(G)\) athari za chanzo kama urekebishaji wa dereva.

Katika chanzo, sehemu ya kwanza ya derivative katika Mlinganyo (22) inaonekana katika uchapishaji kwa umbo hili:

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

Usemi huu unafafanuliwa katika maelezo yaliyo karibu nao kama “quasilinear form”. Ukilinganishwa na noteshini ya kawaida ya mfumo wa hyperboliki, mwelekeo wa derivative ya kwanza unaonekana kuandikwa kinyume katika uchapishaji; hili ni kutokulingana kwa upangaji/noteshini ya derivative ndani ya chanzo. Verianla haisahihishi kimya kimya kana kwamba ni fomula ya chanzo, bali inaashiria kutokulingana huku wazi.

Mpito kwenda katika umbo la momentum lililohifadhiwa: Milinganyo (23)–(36)

Chanzo kinaunda mnyororo wa mabadiliko ya kialjebra ili kufikia muundo wa uhifadhi unaoweza kutumiwa katika kisuluhishi cha Roe. Mlinganyo (23):

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

Mlinganyo (24):

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

Kutokana na haya, Mlinganyo (25) huandikwa:

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

Mlinganyo wa kasi wa PW unapozidishwa kwa wiani, chanzo kinatoa Mlinganyo (26):

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

Kwa sehemu ya usafirishaji wa nafasi, Mlinganyo (27):

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

Mahusiano haya yanapounganishwa:

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

Chanzo kinatoa kama Mlinganyo (29) badiliko la kati linaloeleza upya mtiririko wa momentum kwa kutumia kigeu kilichohifadhiwa \(\rho v\), na kisha hufikia umbo la momentum lililohifadhiwa la PW ya kawaida:

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

Katika Mlinganyo (31) wa chanzo, kuna kutokulingana kwa noteshini kati ya uwasilishaji wa vekta ya hali na maelezo yaliyo mara moja chini yake kwamba “vigeu vilivyohifadhiwa ni wiani na momentum”. Mlinganyo wa baadaye (36) unaonyesha wazi vekta ya hali iliyohifadhiwa kama \([\rho,\rho v]^T\). Kwa hiyo, mawasilisho haya mawili ndani ya chanzo hayachanganywi hapa.

Mlinganyo wa kasi unaofanana na momentum kwa BAM, tena:

\[ 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} \]

Umbo lililozidishwa kwa wiani:

\[ \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} \]

Baada ya kutumia kanuni ya zao na mlinganyo wa uhifadhi wa magari:

\[ (\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} \]

Kwa kutumia kasi \(v=(\rho v)/\rho\), kwa vigeu vilivyohifadhiwa kikamilifu:

\[ (\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} \]

Uwasilishaji wa mwisho wa vekta ya uhifadhi katika chanzo:

\[ 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 na kasi za mawimbi ya tabia

Jacobian ya muundo wa mtiririko wa PW:

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

Tatizo la eigenvalue:

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

Eigenvalue za tabia:

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

Kwa kuwa BAM hubadilisha sehemu ya chanzo bila kubadilisha muundo msingi wa mtiririko, utafiti unahitimisha kuwa PW na BAM zina kasi zilezile za mawimbi ya tabia. Eigenvalue zinapokuwa halisi na tofauti, mfumo huhifadhi muundo wake wa hyperboliki madhubuti. Chanzo pia kinaeleza kuwa katika nukta za soniki ambako eigenvalue hukaribia sifuri, kisuluhishi cha Roe hutumia Harten entropy fix ili kuzuia mishtuko ya upanuzi isiyo ya kifizikia.

Sharti la Uthabiti la Modeli ya BAM ni Lipi?

Katika uchanganuzi wa usumbufu wa mstari wa utafiti, sharti la kimuundo lililopatikana ili BAM iweze kudhoofisha usumbufu mdogo kwa njia ya asymptotic ni \(C_0\geq-\rho_0v_d'(\rho_0)\); derivative ya kasi inayotakiwa inayodhibitiwa na nafasi ya modeli inapowekwa badala yake, kikomo hiki hubadilika kuwa \(C_0\geq\frac{s\rho_0v_{max}}{(1+s\rho_0)^2}\).

Ufafanuzi wa usumbufu mdogo

Karibu na hali sare ya usawa:

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

Hapa inachukuliwa kuwa \(\delta\rho\ll\rho_0\) na \(\delta v\ll v_0\). Yaani, uchanganuzi wa mstari huchunguza usumbufu mdogo karibu na usawa, si mabadiliko makubwa ya trafiki.

Mlinganyo wa uhifadhi uliowekwa mstari:

\[ \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} \]

Mlinganyo wa kasi uliowekwa mstari:

\[ \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} \]

Muda faafu wa ulegezaji katika hali ya usawa:

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

Chanzo, baada ya kuondoa vigeu, kinatoa mlinganyo wa diferensheli sehemu wa daraja la pili kwa usumbufu wa wiani katika umbo lifuatalo:

\[ \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} \]

Ufumbuzi wa wimbi na uhusiano wa dispersion

Chanzo kinashughulikia usumbufu katika umbo la wimbi la harmoniki:

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

Hapa \(k\) ni namba ya wimbi ya nafasi, na \(\omega=\omega_r+i\omega_i\) ni masafa changamano ya angular. Dhana hii inapotumika katika Mlinganyo (43), chanzo kinatoa Mlinganyo (44):

\[ (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} \]

Kwa usumbufu wenye urefu mrefu wa wimbi, yaani \(k\) ndogo, upanuzi unapofanywa kasi ya uenezaji wa wimbi la kinematiki ni:

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

Vipengele vya daraja la pili hutoa kigezo cha ukuaji/kudhoofika:

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

Sharti la uthabiti wa mstari lililopatikana na chanzo:

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

Derivative ya wiani ya kasi inayotakiwa inayodhibitiwa na nafasi:

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

Hivyo kikomo cha mwisho cha uthabiti ni:

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

Uthibitishaji wa namba wa kikomo cha uthabiti

Kwa \(\rho_0=0.30\,veh/m\), \(v_{max}=20\,m/s\) na \(s=8\,m\), chanzo kinatoa:

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

Katika Kielelezo 1, chanzo kinaonyesha profaili za wiani za hali hizo tatu katika nafasi mbalimbali za barabara kwa nyakati tofauti. Katika hali thabiti ya \(C_0=20\,m/s\), amplitude ya wimbi la sinusoidal la awali hupungua kadiri muda unavyopita, ilhali katika hali ya \(C_0=2\,m/s\) usumbufu hukua kwa dhahiri. Grafu haijatumiwa tena hapa kwa nakala halisi; matokeo ya kisayansi yameundwa kwa kujitegemea kwa kutumia thamani za chanzo zilizo hapo juu.

Mbinu na Matokeo ya Utafiti

Mbinu ya ufumbuzi wa namba

Mifumo ya PW na BAM imegawanywa kwa njia ya ujazo wenye ukomo kwa kutumia Roe approximate Riemann solver. Jukumu la mbinu ya Roe katika utafiti ni kufuatilia kwa njia ya namba mabadiliko makali na mawimbi yanayofanana na mshtuko yanayoweza kutokea katika mtiririko wa trafiki, ndani ya mfumo wa sheria za uhifadhi za hyperboliki. Kwa kuwa Jacobian ya mtiririko ya BAM hubaki sawa na ya PW, mgawanyo uleule wa msingi wa Roe unaweza kutumiwa kwa modeli zote mbili.

Chanzo kinaeleza kwamba Harten entropy fix hutumika katika nukta za soniki pale ambapo mojawapo ya eigenvalue hukaribia sifuri. Madhumuni ya marekebisho haya ni kuzuia kisuluhishi cha namba kuzalisha mishtuko ya upanuzi isiyo ya kifizikia.

Vigezo msingi vya uigaji

Kutokulingana kwa hatua ya muda ndani ya chanzo: Jedwali 5 linaandika \(\Delta t=0.005\,s\), na kwa kuwa \(1200\times0.005=6\,s\), jedwali linawiana ndani yake na muda wa jumla. Hata hivyo, katika maandishi ya maelezo kwenye ukurasa wa 18 imeandikwa \(\Delta t=0.006\,s\). Verianla haichagui kimya kimya moja ya thamani hizi mbili kwa niaba ya mtafiti; tofauti kati ya jedwali na simulizi inaachwa wazi. Katika majaribio ya baadaye ya uimara ya sekunde 50, chanzo kinaeleza tena wazi kwamba \(\Delta t=0.005\,s\) ilitumika.

Mlinganyo wa wiani wa awali

Chanzo pia kinafafanua wiani wa awali kwa vipande katika Mlinganyo (49) kama ifuatavyo:

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

Dokezo la sharti la awali ndani ya chanzo: Profaili hii ya vipande si ufafanuzi uleule moja kwa moja na usemi “wiani wa wastani wa awali = 0.30 veh/m” uliotolewa katika Jedwali 5. Chanzo kinatoa taarifa hizi mbili pamoja lakini hakielezi kwa kina uhusiano kati yake. Kwa hiyo, Verianla haizipatanishi kwa njia ya kubuni.

Thamani ya sekunde 0.052 ina maana gani?

Ingawa chanzo kinatumia \(\tau_0=0.52\,s\) kama thamani ya marejeleo ya mwitikio wa binadamu, hakiwasilishi \(0.052\,s\) kama muda wa kawaida wa mwitikio wa dereva wa binadamu. Thamani hiyo ni kigezo cha unyeti kilichodhibitiwa, kilichoundwa kwa kupunguza thamani msingi kwa oda moja ya ukubwa. Vivyo hivyo, nafasi ya standstill ya 2 m haitumiwi kama thamani ya ulimwengu mzima ya tabia ya kawaida ya dereva, bali kama sharti la kulinganisha unyeti.

Ulinganisho wa tabia ya kasi ya PW na BAM

Katika ulinganisho msingi kwenye barabara ya mzunguko uliofungwa yenye urefu wa 500 m, modeli ya PW imeonyesha misukosuko mikali ya namba ambayo chanzo kinaieleza kuwa “haiwezekani kifizikia” chini ya mabadiliko ya ghafla ya trafiki. Chanzo kinaripoti excursion za kasi ya juu zaidi ya takribani \(320\,m/s\) na ya chini zaidi ya \(-180\,m/s\).

Katika masharti yale yale ya majaribio, BAM imeweka kasi ndani ya masafa ya modeli yaliyowekwa ya 0–20 m/s na kutoa mabadiliko ya kasi yaliyo laini zaidi. Hata hivyo, waandishi hawatafsiri ulinganisho huu kama ubora wa usahihi wa kiempiria. Matokeo yanaonyesha tu kwamba, chini ya masharti ya namba yaliyojaribiwa, BAM ilitoa tabia ya kasi ya namba yenye mipaka zaidi na inayokubalika zaidi kifizikia.

Sehemu ya wiani imebadilikaje?

Katika uigaji wa wiani wa sekunde 50 wa BAM kwenye barabara ya mzunguko wa 500 m, chanzo kinaripoti karibu na \(x=150\,m\) upeo wa ndani wa takribani \(0.45\,veh/m\), na karibu na \(x=375\,m\) kiwango cha chini cha takribani \(0.15\,veh/m\). Sehemu hiyo hubaki laini kwa ujumla na haizalishi kutokuwepo mwendelezo kunakoonekana wazi.

Katika sehemu ya wiani ya PW ya kawaida, chanzo kinaripoti karibu na \(150\,m\) upeo wa \(0.42\,veh/m\), na karibu na \(350\,m\) kiwango cha chini cha \(0.17\,veh/m\). Kadiri muda unavyopita, profaili laini ya awali hubadilika kuwa miundo mikali zaidi na mawimbi ya stop-and-go yanayoenea kurudi nyuma. Tafsiri ya utafiti ni kwamba, chini ya masharti yaliyojaribiwa, PW huongeza usumbufu wa awali, ilhali BAM hutoa uenezaji wa mawimbi ulio laini zaidi.

Matokeo ya mtiririko wa trafiki

Kwa PW, katika uchunguzi wa mtiririko wa barabara ya 500 m kwa sekunde 6, chanzo kinaripoti karibu na \(150\,m\) upeo wa \(4.55\,veh/s\), na karibu na \(350\,m\) kiwango cha chini cha \(3.4\,veh/s\). Umbo la wimbi hupoteza ulinganifu wake kadiri muda unavyopita na huonyesha deformation.

Kwa BAM, upeo wa mtiririko umeripotiwa karibu na \(125\,m\) kuwa \(4.55\,veh/s\), na kiwango cha chini karibu na \(375\,m\) kuwa \(3.4\,veh/s\). Ni muhimu kwamba ukubwa wa mtiririko uko karibu: kama waandishi wanavyosema pia katika sehemu ya hitimisho, tofauti kuu ya BAM katika majaribio si kuzalisha “mtiririko mkubwa zaidi wa trafiki”, bali ni maendeleo yaliyodhibitiwa zaidi ya hali za trafiki katika muda na nafasi.

Majaribio ya uimara yalibuniwaje?

Utafiti kisha ulichunguza BAM kwa sekunde 50 chini ya masharti manne kwenye barabara za mzunguko uliofungwa zenye urefu wa 250, 500 na 1000 m. Katika majaribio haya, hatua ya nafasi ilihifadhiwa kuwa 1 m, hatua ya muda 0.005 s, \(C_0=20\,m/s\), \(s_0=0.30\) na \(v_{max}=20\,m/s\).

Katika Senario D, kikomo cha uwezo wa ndani kinachoweza kurudiwa kiliundwa kwa kupunguza kasi ya juu inayopatikana kwa %25 katika sehemu ya %10 iliyo katikati ya barabara.

Uwiano wa usumbufu wa RMS unaonyesha nini?

Utafiti hutumia uwiano wa usumbufu wa mwisho wa wiani wa RMS kwa thamani ya awali ya RMS ili kutathmini ni kiasi gani cha usumbufu wa wiani wa awali kinabaki mwishoni mwa uigaji. Uwiano unapokuwa chini ya 1, usumbufu kwa ujumla umedhoofika; unapozidi 1, usumbufu wa mwisho umekuwa mkubwa kuliko wa awali.

Katika Senario A, kadiri urefu wa barabara ulivyoongezeka, uwiano wa mwisho wa RMS ulipanda kutoka 0.002 hadi 0.213 na 0.681. Yaani, ingawa usumbufu wa awali karibu ulidhoofishwa kabisa kwenye barabara fupi, barabara ndefu zaidi ilibeba mabadiliko zaidi ya nafasi yaliyobaki.

Katika Senario B, licha ya wiani wa juu zaidi wa awali na usumbufu wenye nguvu zaidi, ufumbuzi ulibaki wenye mipaka katika urefu wote mitatu wa barabara uliojaribiwa. Uwiano wa RMS ni 0.002, 0.219 na 0.686, kwa mtiririko huo.

Katika Senario C, kipimo msingi cha muda wa ulegezaji kilicho kirefu zaidi na nafasi kubwa zaidi ya standstill vilitumika. Ingawa maandishi ya chanzo yanaeleza uwiano wa RMS wa takribani 0.000001 kwa 250 m, Jedwali 7 linaonyesha 0.000 kwa nafasi tatu za desimali. Tofauti hii iko katika kiwango cha kuzungusha; hapa thamani ya jedwali imehifadhiwa jinsi ilivyo ndani ya jedwali.

Usumbufu uliobaki wenye nguvu zaidi ulitokea katika Senario D yenye sehemu finyu amilifu. Kwa barabara za 500 m na 1000 m, uwiano wa RMS ni 1.332 na 1.923, kwa mtiririko huo; yaani usumbufu wa mwisho wa wiani ni mkubwa kuliko thamani yake ya awali. Chanzo kinatafsiri hili kama ugawaji upya wa kudumu wa nafasi na ukuaji wa usumbufu unaosababishwa na kikomo cha uwezo wa ndani kilichotumika.

Katika maandishi ya chanzo, kasi ya juu zaidi ya Senario D kwa 1000 m imeandikwa kuwa 6.728 m/s, ilhali thamani katika Jedwali 7 ni 6.729 m/s. Verianla inahifadhi thamani ya 6.729 m/s ya jedwali la chanzo katika jedwali la namba na haifichi tofauti ya 0.001 m/s kati ya simulizi na jedwali.

Ulinganisho wa barabara za mita 250, 500 na 1000

Katika Vielelezo 8–15, chanzo kinatoa profaili za mwisho za wiani na kasi kwa senario nne katika urefu wa barabara wa 250, 500 na 1000 m. Mwelekeo wa jumla ni kwamba barabara fupi hudhoofisha usumbufu wa awali kwa nguvu zaidi, ilhali barabara ndefu huruhusu mabadiliko zaidi ya nafasi kubaki. Matokeo haya yanaonekana wazi hasa katika uwiano wa RMS.

Hata hivyo, hali hii haiwezi kujumlishwa kuwa “barabara fupi huwa thabiti zaidi kila wakati katika trafiki halisi”. Athari ya urefu wa barabara katika utafiti ni matokeo ya namba yaliyopatikana chini ya mzunguko uliofungwa, sharti la mpaka la kipindi, vigezo vya BAM vilivyochaguliwa na usumbufu wa awali uliofafanuliwa.

Unyeti wa muda wa mwitikio na nafasi ya standstill

Vielelezo 16 na 17 vya chanzo vinalinganisha seti mbili za vigezo kwenye barabara ya mzunguko wa 500 m. Katika hali ya kwanza, \(\tau=0.52\,s\) na \(l_s=8\,m\); katika hali ya pili, kipimo kifupi zaidi cha muda cha \(0.052\,s\) na nafasi ya 2 m zimetumiwa. Chanzo kinaripoti kwamba mawimbi ya wiani ni laini zaidi kwa kipimo kikubwa zaidi cha muda na nafasi ya standstill, ilhali mabadiliko makubwa zaidi ya kasi na wiani hutokea kwa thamani ndogo.

Thamani za 0.052 s na 2 m hazipaswi kuchukuliwa kuwa zinawakilisha uendeshaji wa kawaida wa binadamu. Makala inatumia thamani hizi waziwazi kama masharti ya unyeti yaliyodhibitiwa.

Matokeo yanayoungwa mkono na utafiti

• Inaungwa mkono kwamba, chini ya masharti ya namba yaliyojaribiwa, BAM inaweza kutoa thamani za kasi zenye mipaka zaidi na maendeleo laini zaidi ya mawimbi ya trafiki kuliko PW ya kawaida.
• Matumizi ya pamoja ya kasi inayotakiwa inayodhibitiwa na nafasi na kipimo cha muda wa ulegezaji kinachotegemea wiani hurekebisha sehemu ya chanzo ya modeli kulingana na hali ya trafiki.
• Uchanganuzi wa mstari unatoa kikomo cha wazi cha uthabiti kwa BAM.
• Majaribio ya \(C_0=2.0\), \(4.15\) na \(20.0\,m/s\) yameonyesha tabia ya namba inayolingana na kikomo cha kinadharia cha uthabiti.
• Katika majaribio ya uimara ya 250–1000 m, ufumbuzi wa BAM ulibaki wenye mipaka ndani ya nafasi ya vigezo iliyochunguzwa.
• Sehemu finyu amilifu ilisababisha usumbufu mkubwa zaidi wa wiani uliobaki kati ya masharti yaliyochunguzwa.
• Ukubwa wa mtiririko wa trafiki wa BAM na PW unaweza kufanana; tofauti kuu katika utafiti ilijitokeza zaidi katika uthabiti wa hali na maendeleo yake katika muda na nafasi.

Matokeo yasiyoungwa mkono na utafiti

• Haijaonyeshwa kwamba BAM ni sahihi zaidi kwa uhakika kuliko PW katika mitandao halisi ya barabara.
• Modeli haijakalibishwa kwa sensa halisi za trafiki, loop detector au data za trajectory za magari.
• Utendaji wa utabiri haujapimwa katika seti huru ya data ya ulimwengu halisi.
• Urefu tofauti wa magari, aina za magari na sifa za madereva binafsi hazitatuliwi waziwazi.
• Mbinu ya \(\Delta x\approx1/\rho\) katika utafiti haiwakilishi kikamilifu nafasi zote zisizo sawia kati ya magari.
• Kigezo cha sekunde 0.052 si matokeo ya muda halisi wa mwitikio wa dereva wa binadamu.
• Nafasi ya standstill ya 2 m haiwasilishwi kama kanuni ya jumla ya trafiki halisi.
• Uthabiti wa namba katika uigaji uliodhibitiwa si ushahidi wa usahihi wa utabiri wa trafiki halisi.

Mapungufu msingi ya utafiti

Kizuizi muhimu zaidi ambacho waandishi wameeleza wazi ni kutokuwepo kwa ukalibishaji na uthibitishaji wa kiempiria. Matokeo ya sasa yanategemea majaribio ya namba yaliyodhibitiwa kwenye barabara ya mzunguko uliofungwa na usumbufu wa awali uliobainishwa mapema.

BAM pia hutoa nafasi ya mbele ya gari kutoka kwa wiani wa makroskopiki kwa kutumia mbinu ya \(\Delta x\approx1/\rho\). Mbinu hii hubeba dhana ya nafasi zilizo takribani sawia na haitatui waziwazi utofauti unaotokana na aina tofauti za magari, urefu wa magari au sifa za madereva binafsi.

Kwa kazi ya baadaye, waandishi wanapendekeza kukalibisha vigezo vya BAM kwa data halisi za trafiki; kuthibitisha kwa vipimo huru kama mtiririko wa trafiki, kasi, wiani, headway na muda wa safari; na kuchunguza hali za trafiki zisizo sawia, unyeti wa vigezo na miundo ya mtandao iliyo halisi zaidi.

Dokezo la Chanzo na Mbinu

Kichwa asilia: A macroscopic traffic flow model with gap-controlled driver response and adaptive reaction dynamics

Waandishi na mpangilio: Benjamin Simwinga; Aaron Zimba; Mayumbo Nyirenda.

Mwandishi wa mawasiliano: Benjamin Simwinga.

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

Jarida: Discover Networks.

Mchapishaji/jukwaa: Springer Nature / Springer Nature Link.

Juzuu na namba ya makala: Volume 2, Article 39 (2026).

DOI: 10.1007/s44354-026-00053-y

Tarehe ya kuchapishwa: 8 Oktoba 2026.

Tarehe za mchakato: Ilipokelewa: 10 Agosti 2026; ilifanyiwa marekebisho: 19 Septemba 2026; ilikubaliwa: 28 Septemba 2026.

Aina ya chanzo: Research article; uundaji wa modeli ya kihisabati, uchanganuzi wa uthabiti wa mstari na utafiti wa uigaji wa namba uliodhibitiwa.

Leseni: Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0). Leseni inaruhusu ufikiaji wazi lakini hairuhusu usambazaji wa nyenzo zilizobadilishwa; kwa sababu hiyo, simulizi ya Verianla imeandaliwa kwa usanifu wa taarifa za kisayansi ulio huru na maandishi na vielelezo vya chanzo.

Michango ya waandishi: Chanzo kinaeleza kwamba BS aliandika makala; AZ na MN walitoa mapitio na mwongozo.

Ufadhili: Utafiti unaripoti kwamba haukupokea ufadhili.

Upatikanaji wa data: Imeelezwa kwamba seti za data zilizozalishwa na/au kuchanganuliwa wakati wa utafiti zinaweza kupatikana kutoka kwa mwandishi wa mawasiliano kwa ombi linalofaa.

Mgongano wa maslahi: Waandishi wanaripoti kwamba hakuna mgongano wa maslahi.

Idhini ya kimaadili na ridhaa ya washiriki: Imeripotiwa kuwa haitumiki.

Hali ya uthibitishaji wa kiempiria: Hakuna. BAM bado haijakalibishwa kwa data za trafiki zilizozingatiwa wala kuthibitishwa kwa vipimo huru vya uwanjani.

Uhitimishaji wa kitakwimu: Utafiti hauripoti thamani ya p, muda wa kujiamini, takwimu za sampuli za kliniki/majaribio au jaribio la nadharia ya kitakwimu. Matokeo ni masharti ya uchanganuzi ya uthabiti na matokeo ya uigaji wa namba.

Kutokulingana kwa kiufundi ndani ya chanzo: Ingawa Jedwali 5 linatoa \(\Delta t=0.005\,s\), sehemu moja ya maandishi ya maelezo inaandika 0.006 s. Thamani ya 0.005 s ya Jedwali 5 inawiana ndani yake na hatua 1200 na muda wa jumla wa sekunde 6; hata hivyo, Verianla haijafuta kimya kimya usemi wa 0.006 s katika maandishi ya chanzo kwa sababu hiyo. Uhusiano kati ya “wiani wa wastani wa awali” wa 0.30 veh/m uliotolewa katika Jedwali 5 na profaili ya awali ya vipande 0.01/0.30/0.10 katika Mlinganyo (49) pia hauelezwi kwa kina katika chanzo. Pia kuna kutokulingana kwa upangaji/noteshini katika marejeo ya mlinganyo na baadhi ya uwasilishaji wa vekta.