Benjamin-Aaron-Mayumbo (BAM) модели — транспорттун тыгыздыгы жана орточо ылдамдык сыяктуу жыйынды өзгөрмөлөр менен иштеген, классикалык Payne-Whitham (PW) моделин алдындагы унаага чейин жеткиликтүү болгон аралыкка көз каранды айдоочу реакциясы менен кеңейтүүчү экинчи тартиптеги макроскопиялык транспорт агымынын модели. Моделдин негизги жаңылыгы — алдыдагы унаалар аралык аралыкты каалаган ылдамдыкты аныктоо үчүн гана эмес, айдоочунун учурдагы ылдамдыктан каалаган ылдамдыкка канчалык тез жакындай турганын аныктаган эффективдүү релаксация убактысын да өзгөртүү үчүн колдонушу. Изилдөө BAMды аналитикалык туруктуулук талдоосу жана көзөмөлдөнгөн сандык симуляциялар аркылуу классикалык PW модели менен салыштырат.

Модель макроскопиялык тыгыздыктан болжолдуу унаалар аралык аралыкка өтүү үчүн \(\Delta x \approx 1/\rho\) катышын колдонот. Бул ыкма транспорт агымындагы унаалар болжол менен бир тектүү аралыктар менен жайгашкан деген божомолго негизделет. Алынган аралык андан кийин эки өзүнчө механизмге кирет: айдоочу жетүүгө умтулган аралыкка жараша башкарылган ылдамдык \(v_d\) жана ылдамдыкты жөндөөнүн убакыт масштабын аныктаган эффективдүү релаксация убактысы \(\tau_{eff}\). Ошентип, транспорттун тыгыздыгы жогорулаганда алдыда жеткиликтүү болгон аралык кыскарат; модель максаттуу ылдамдыкты да, ылдамдыкты жөндөө жүрүм-турумун да ошол эле транспорт абалына байланыштырат.

Изилдөөнүн негизги салыштырууларында классикалык PW модели сыналган шарттарда өтө чоң сандык ылдамдык термелүүлөрүн жаратса, BAM ылдамдыктарды модель үчүн каралган диапазондун ичинде кармап, кыйла жылмакай транспорт толкундарын түзгөн. Булак PW симуляциясында ылдамдык болжол менен \(320\,m/s\) чейин көтөрүлүп, болжол менен \(-180\,m/s\) чейин түшүшү мүмкүн экенин билдирет, ал эми BAM 0–20 \(m/s\) деп белгиленген моделдөө диапазонунда калганын белгилейт. Бул BAM каралган сандык түзүлүштө көбүрөөк көзөмөлдөнгөн жүрүм-турум көрсөткөнүн билдирет; бул анын реалдуу транспортто такыраак божомол берерин билдирбейт.

Сызыктуу туруктуулук талдоосунда изилдөө чыгарган структуралык шарт \(C_0 \geq -\rho_0 v_d'(\rho_0)\) түрүндө. Аралыкка жараша башкарылган каалаган ылдамдык функциясынын туундусу ордуна коюлганда чек \(C_0 \geq \frac{s\rho_0v_{max}}{(1+s\rho_0)^2}\) болот. \(\rho_0=0.30\,veh/m\), \(v_{max}=20\,m/s\) жана \(s=8\,m\) үчүн изилдөө критикалык маанини болжол менен \(C_{0,critical}=4.15\,m/s\) деп эсептеген. \(C_0=2.0\,m/s\) туруксуз, \(4.15\,m/s\) критикалык чек жана \(20.0\,m/s\) туруктуу абал катары сандык түрдө сыналган.

Эң маанилүү чечмелөө чеги төмөнкүдөй: BAM али реалдуу транспорт байкоолору менен калибрленген эмес же көз карандысыз түрдө текшерилген эмес. Натыйжалар жабык шакек жолдордо, аныкталган баштапкы пертурбациялар жана тандалган параметрлер астында жүргүзүлгөн көзөмөлдөнгөн сандык эксперименттердин жыйынтыктары. Ошондуктан изилдөө BAM реалдуу жол тармактарында PWге караганда сөзсүз такыраак божомол берерин көрсөтпөйт; ал көбүрөөк көзөмөлдөнгөн жана туруктуу сандык жүрүм-турум жарата аларын көрсөтөт.

BAM Транспорт Агымынын Модели Деген Эмне?

Benjamin-Aaron-Mayumbo (BAM) модели классикалык Payne-Whitham моделинин унааларды сактоо жана транспорт-басым/алдын ала эсепке алуу түзүмүн сактаган, бирок туруктуу релаксация жүрүм-турумунун ордуна транспорт тыгыздыгынан чыгарылган алдыдагы унаалар аралык аралыкка көз каранды каалаган ылдамдыкты жана эффективдүү релаксация убактысын колдонгон экинчи тартиптеги макроскопиялык транспорт агымынын модели.

Макроскопиялык транспорт модели эмнени чагылдырат?

Микроскопиялык транспорт моделдери ар бир унааны жана айдоочуну өзүнчө көзөмөлдөсө, макроскопиялык моделдер транспортту агып жаткан үзгүлтүксүз чөйрө катары карайт. Ошондуктан ар бир унаанын рул, газ же тормоз буйругун өзүнчө моделдөөнүн ордуна жол боюнча тыгыздык, орточо ылдамдык жана агым сыяктуу талаа өзгөрмөлөрү колдонулат.

Бул изилдөөдөгү негизги өзгөрмөлөр төмөнкүлөр:

LWR модели: унаалардын сакталуусу

Изилдөөнүн математикалык баштапкы чекиттеринин бири Lighthill-Whitham-Richards (LWR) модели. Булактагы Теңдеме (1):

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

Бул теңдеме сакталуу мыйзамы. \(\rho v\) көбөйтүндүсү транспорт агымын билдирет. Белгилүү бир жол кесиндисиндеги тыгыздык убакыт өткөн сайын өзгөрсө, анын себеби унаалардын ошол кесиндиге кирип-чыгышы. Теңдеме өзүнчө айдоочу кандай тездетерин же жайлатарын айтпайт; ал унаалардын санынын сакталуусун гана билдирет.

Payne-Whitham модели эмне үчүн экинчи теңдемеге муктаж?

LWR модели ылдамдыктын динамикалык өзгөрүшүн өзүнчө теңдеме менен түшүндүрбөйт. Ошондуктан Payne-Whitham модели сакталуу теңдемесине ылдамдык/импульс теңдемесин кошот:

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

Теңдеме (3)тө \(\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}\) болсо айдоочунун учурдагы ылдамдыгын тең салмак ылдамдыкка карай жөндөөсүн билдирген релаксация мүчөсү.

Классикалык PW моделиндеги көйгөй — \(\tau\) туруктуу деп алынат. BAMдын негизги өзгөртүүсү дал ушул булак мүчөсүнө багытталган.

Аралыкка жараша башкарылган реакция мүчөсү

Изилдөөнүн Теңдеме (4)ү:

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

Алымдын биринчи бөлүгү айдоочунун учурдагы аралык шартында жетүүнү каалаган ылдамдыгын, бөлүмү болсо ылдамдык бул мааниге кайсы убакыт масштабында жакындай турганын билдирет. Ошентип, BAM “канчалык тез жүргүм келет?” деген суроону гана эмес, “ошол ылдамдыкка канчалык тез жакындайм?” деген суроону да унаалар аралык аралыкка байланыштырат.

Greenshields ылдамдык-тыгыздык катышы

Булак классикалык тең салмак катышы үчүн Greenshields моделин колдонот:

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

Тыгыздык нөлгө жакындаганда тең салмак ылдамдык \(v_{max}\)ка жакындайт. Тыгыздык \(\rho_m\)ге жакындаганда теңдемедеги тең салмак ылдамдык нөлгө жакындайт.

Теңдеме (6) — бул катыштын тыгыздык боюнча туундусу:

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

Терс белги тыгыздык жогорулаган сайын тең салмак ылдамдык азаярын математикалык түрдө көрсөтөт.

Эффективдүү Релаксация Убактысы Чыныгы Айдоочу Реакциясынын Кечигүүсүбү?

Жок. BAMдагы \(\tau_{eff}\) айдоочу эскертүүнү көргөндөн кийин белгилүү убакыт күтүүсүн билдирген ачык убакыт кечигүүсү эмес; ал учурдагы ылдамдыктын аралыкка жараша башкарылган каалаган ылдамдыкка жакындаш ылдамдыгын аныктаган, транспорт абалына көз каранды убакыт масштабы. Булак бул айырманы өзгөчө баса белгилейт жана чыныгы кечигүү моделинде система \(t-\tau\) сыяктуу өткөн убактагы транспорт абалына көз каранды болушу керек экенин айтат.

Аралыкка жараша башкарылган каалаган ылдамдык

Теңдеме (7):

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

\(\Delta x\) чоңойгон сайын \(\frac{\Delta x}{\Delta x+s}\) катышы 1ге жакындайт жана \(v_d\) эркин агым ылдамдыгы \(v_{max}\)ка жакындайт. Унаалар бири-бирине жакындаган сайын катыш кичирейип, каалаган ылдамдык төмөндөйт.

Аралыкка көз каранды эффективдүү релаксация убактысы

Теңдеме (8):

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

Бул формула BAMдын классикалык туруктуу \(\tau\) ыкмасынан айырмаланган экинчи маанилүү жери. \(\tau_{eff}\) транспорт абалына жараша өзгөрөт.

Булакта \(s_0\) “spacing scaling parameter” деп аныкталып, 5-таблицада 0.30 мааниси берилген. Бирок бул параметрдин бирдиги таблицада ачык көрсөтүлгөн эмес. Теңдеменин түзүлүшүндө \(\Delta x+s_0\) суммасына киргендиктен, өлчөмдүк чечмелөө жагынан бул жетишпестик маанилүү; Verianla бул боштукка булакта жок бирдикти кошпойт.

Ылдамдануу теңдемеси

Теңдеме (9):

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

Каалаган ылдамдык учурдагы ылдамдыктан чоң болсо, алымы оң болуп, модель ылдамдануу багытында; учурдагы ылдамдык каалаган ылдамдыктан жогору болсо, алымы терс болуп, жайлоо багытында реакция жаратат. \(\tau_{eff}\) болсо бул ыңгайлашуу канчалык тез болорун масштабдайт.

Теңдеме (7) жана Теңдеме (8)деги туюнтмалар Теңдеме (9)да бириктирилгенде, булакта Теңдеме (10) катары берилген аралыкка жараша башкарылган ылдамдануу түзүмү алынат:

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

Булак ичиндеги номерлөө эскертүүсү: PDF түшүндүрмө тексти бул чыгаруунун айрым жерлеринде теңдеме номерлерине туура эмес кайчылаш шилтемелерди берет. Бул жерде формулалар булакта басылган математикалык туюнтмага ылайык иреттелген; кайчылаш шилтеме каталары булак текстинин бөлүгү катары унчукпай кайра чыгарылган эмес.

Микроскопиялык аралыктан макроскопиялык тыгыздыкка өтүү

Изилдөөнүн негизги ыкмасы:

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

Бул катыш унаалар болжол менен бир тектүү жайгашкан транспорт агымында тыгыздыктын тескери маанисин орточо алдыдагы аралык үчүн макроскопиялык жакындоо катары колдонот. \(\rho\) veh/m менен берилгенде \(1/\rho\) бир унаага туура келген узундук масштабын берет. Булак да бул өзгөртүү унаалар аралык аралык салыштырмалуу бир тектүү деген божомолго негизделерин ачык айтат.

Теңдеме (11) каалаган ылдамдык теңдемесине колдонулганда:

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

Ушундай эле эффективдүү релаксация убактысы:

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

Ошентип, BAMдын ылдамдануу туюнтмасын түздөн-түз транспорт тыгыздыгы менен жазууга болот:

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

Кайра иреттелген түрү:

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

Теңдеме (15) BAMдын маңызын абдан ачык көрсөтөт. Тыгыздык \(\rho\) максаттуу ылдамдык \(\frac{v_{max}}{1+s\rho}\)ка да, ыңгайлашуу масштабы \(\frac{1+s_0\rho}{\tau_0}\)на да таасир этет.

BAMдын акыркы ылдамдык теңдемеси

Классикалык PW релаксация мүчөсү жаңы тыгыздыкка көз каранды булак мүчөсү менен алмаштырылганда:

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

Бул теңдеме унааларды сактоо теңдемеси менен бирге BAM системасын түзөт:

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

PW менен BAMдын негизги математикалык айырмасы

Roe ажыратуусу жана сакталуу формасы

BAM жана PW теңдемелери сандык чечим үчүн Roe approximate Riemann solver колдонулуп, чектүү көлөмдөр түрүндө дискреттештирилген. Булак гиперболалык баланс мыйзамын төмөнкү жалпы формада берет:

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

Булактагы биринчи вектордук көрсөтүлүш:

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

Бул жерде \(G\) транспорт абалын, \(f(G)\) бул абалдын мейкиндик боюнча ташылышын жана \(S(G)\) айдоочунун ыңгайлашуусу сыяктуу булак таасирлерин билдирет.

Булакта Теңдеме (22)нин биринчи туунду мүчөсү басмада төмөнкүдөй көрүнөт:

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

Бул туюнтма дароо жанындагы түшүндүрмөдө “quasilinear form” деп аныкталат. Стандарттык гиперболалык система нотасиясы менен салыштырганда, биринчи туундунун багыты басмада тескери жазылгандай көрүнөт; бул булак ичиндеги терүү/туунду көрсөтүү дал келбестиги. Verianla муну булак формуласы катары унчукпай оңдобойт, дал келбестикти ачык белгилейт.

Сакталуучу импульс формасына өтүү: Теңдеме (23)–(36)

Булак Roe чечүүчүсүндө колдонууга боло турган сакталуучу түзүмгө жетүү үчүн алгебралык өзгөртүүлөрдүн тизмегин түзөт. Теңдеме (23):

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

Теңдеме (24):

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

Булардын негизинде Теңдеме (25) жазылат:

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

PW ылдамдык теңдемеси тыгыздыкка көбөйтүлгөндө булак Теңдеме (26)ны берет:

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

Мейкиндиктик ташуу мүчөсү үчүн Теңдеме (27):

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

Бул катыштар бириктирилгенде:

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

Булак импульс агымын сакталуучу өзгөрмө \(\rho v\) менен кайра туюнткан аралык өзгөртүүнү Теңдеме (29) катары берет, андан кийин классикалык PWнин сакталуучу импульс формасына келет:

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

Булактын Теңдеме (31)индеги абал векторунун көрсөтүлүшү менен анын дароо астындагы “сакталуучу өзгөрмөлөр — тыгыздык жана импульс” деген түшүндүрмөнүн ортосунда нотасиялык дал келбестик бар. Кийинки Теңдеме (36) сакталуучу абал векторун ачык түрдө \([\rho,\rho v]^T\) деп берет. Ошондуктан бул жерде булактын ичиндеги эки көрсөтүлүш бири-бирине аралаштырылбайт.

BAM үчүн импульска окшош ылдамдык теңдемеси кайрадан:

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

Тыгыздыкка көбөйтүлгөн түрү:

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

Көбөйтүндү эрежеси жана унааларды сактоо теңдемеси колдонулгандан кийин:

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

Ылдамдык \(v=(\rho v)/\rho\) колдонулуп, толугу менен сакталуучу өзгөрмөлөр менен:

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

Булактын акыркы сакталуучу вектордук көрсөтүлүшү:

\[ 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 жана мүнөздөмөлүү толкун ылдамдыктары

PW агым түзүмүнө тиешелүү Jacobian:

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

Өздүк маани маселеси:

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

Мүнөздөмөлүү өздүк маанилер:

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

BAM булак мүчөсүн өзгөртүп, негизги агым түзүмүн өзгөртпөгөндүктөн, изилдөөдө PW менен BAM бирдей мүнөздөмөлүү толкун ылдамдыктарына ээ деген жыйынтык чыгарылат. Өздүк маанилер чыныгы жана бири-биринен айырмалуу болгондо система катуу гиперболалык түзүмүн сактайт. Булак ошондой эле өздүк маанилер нөлгө жакындаган соникалык чекиттерде Roe чечүүчүсүнүн физикалык эмес кеңейүү шокторун алдын алуу үчүн Harten entropy fix колдонулганын билдирет.

BAM Моделинин Туруктуулук Шарты Кандай?

Изилдөөнүн сызыктуу пертурбация талдоосунда BAMдын кичине бузулууларды асимптотикалык түрдө басаңдата алышы үчүн алынган структуралык шарт \(C_0\geq-\rho_0v_d'(\rho_0)\); моделдин аралыкка жараша башкарылган каалаган ылдамдык туундусу ордуна коюлганда бул чек \(C_0\geq\frac{s\rho_0v_{max}}{(1+s\rho_0)^2}\) түрүнө өтөт.

Кичине пертурбацияларды аныктоо

Бир тектүү тең салмак абалынын айланасында:

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

Бул жерде \(\delta\rho\ll\rho_0\) жана \(\delta v\ll v_0\) деп кабыл алынат. Башкача айтканда, сызыктуу талдоо чоң транспорт өзгөрүүлөрүн эмес, тең салмактын айланасындагы кичине бузулууларды изилдейт.

Сызыктууланган сакталуу теңдемеси:

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

Сызыктууланган ылдамдык теңдемеси:

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

Тең салмак абалындагы эффективдүү релаксация убактысы:

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

Булак өзгөрмөлөр жок кылынганда, тыгыздык пертурбациясы үчүн экинчи тартиптеги жеке туундулуу дифференциалдык теңдемени төмөнкүдөй берет:

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

Толкун чечими жана дисперсиялык катыш

Булак пертурбацияны гармониялык толкун түрүндө карайт:

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

Бул жерде \(k\) — мейкиндиктик толкун саны, \(\omega=\omega_r+i\omega_i\) — комплекстүү бурчтук жыштык. Бул божомол Теңдеме (43)кө колдонулганда булак Теңдеме (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} \]

Узун толкундуу, башкача айтканда кичине \(k\) бузулуулар үчүн кеңейтүү жүргүзүлгөндө кинематикалык толкундун таралуу ылдамдыгы:

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

Экинчи тартиптеги мүчөлөр өсүү/басаңдоо параметрин берет:

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

Булак алган сызыктуу туруктуулук шарты:

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

Аралыкка жараша башкарылган каалаган ылдамдыктын тыгыздык боюнча туундусу:

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

Ошентип, акыркы туруктуулук чеги:

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

Туруктуулук чегин сандык текшерүү

Булак \(\rho_0=0.30\,veh/m\), \(v_{max}=20\,m/s\) жана \(s=8\,m\) үчүн:

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

Булак 1-сүрөттө үч абалдын жол боюнча тыгыздык профилдерин ар башка убакыттарда көрсөтөт. Туруктуу \(C_0=20\,m/s\) абалында баштапкы синусоидалык толкундун амплитудасы убакыт өткөн сайын азайса, \(C_0=2\,m/s\) абалында бузулуу ачык өсөт. График бул жерде түз көчүрүлүп колдонулган эмес; илимий жыйынтык жогорудагы булак маанилери менен өз алдынча түзүмдөштүрүлгөн.

Изилдөөнүн Ыкмасы жана Жыйынтыктары

Сандык чечүү ыкмасы

PW жана BAM системалары Roe approximate Riemann solver колдонулуп, чектүү көлөмдөр ыкмасы менен дискреттештирилген. Roe ыкмасынын изилдөөдөгү ролу — транспорт агымында пайда болушу мүмкүн болгон кескин өзгөрүүлөрдү жана шокко окшош толкундарды гиперболалык сакталуу мыйзамдарынын алкагында сандык түрдө көзөмөлдөө. BAMдын агым Jacobian-ы PW менен бирдей калгандыктан, эки модель үчүн бирдей негизги Roe ажыратуусун колдонууга болот.

Булак соникалык чекиттерде өздүк маанилердин бири нөлгө жакындаганда Harten entropy fix колдонуларын билдирет. Бул оңдоонун максаты — сандык чечүүчүнүн физикалык эмес кеңейүү шокторун жаратышын алдын алуу.

Негизги симуляция параметрлери

Булак ичиндеги убакыт кадамынын дал келбестиги: 5-таблицада \(\Delta t=0.005\,s\) деп жазылган жана \(1200\times0.005=6\,s\) болгондуктан таблица өз ичинде жалпы убакыт менен дал келет. Бирок 18-беттеги түшүндүрмө текстте \(\Delta t=0.006\,s\) деп жазылган. Verianla бул эки маанинин бирин изилдөөчүнүн ордуна унчукпай тандабайт; таблица менен баяндоонун ортосундагы айырманы ачык калтырат. Кийинки 50 секунддук робасттуулук эксперименттеринде булак кайрадан \(\Delta t=0.005\,s\) колдонулганын ачык билдирет.

Баштапкы тыгыздык теңдемеси

Булак ошондой эле Теңдеме (49)да баштапкы тыгыздыкты бөлүктөр боюнча төмөнкүдөй аныктайт:

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

Булак ичиндеги баштапкы шарт эскертүүсү: Бул бөлүктүк профиль 5-таблицада “mean initial density = 0.30 veh/m” деп берилген туюнтма менен түздөн-түз бирдей аныктама эмес. Булак эки маалыматты тең берет, бирок алардын өз ара байланышын кеңири түшүндүрбөйт. Ошондуктан Verianla аларды жасалма түрдө шайкеш келтирбейт.

0.052 секунддук маани эмнени билдирет?

Булак \(\tau_0=0.52\,s\)ти эталондук адам реакциясынын мааниси катары колдонгону менен, \(0.052\,s\)ти типтүү адам айдоочусунун реакция убактысы катары сунуштабайт. Бул маани негизги маанини бир иреттик чоңдукка азайтуу менен түзүлгөн көзөмөлдөнгөн сезгичтик параметри. Ошол сыяктуу 2 m standstill аралыгы да кадимки айдоочу жүрүм-турумунун универсалдуу мааниси катары эмес, салыштырма сезгичтик шарты катары колдонулат.

PW жана BAM ылдамдык жүрүм-турумун салыштыруу

500 m жабык шакек жолдогу негизги салыштырууда PW модели кескин транспорт өзгөрүүлөрүнүн шартында булакта “физикалык жактан мүмкүн эмес” деп мүнөздөлгөн күчтүү сандык термелүүлөрдү көрсөткөн. Булак болжол менен \(320\,m/s\) максималдуу жана \(-180\,m/s\) минималдуу ылдамдык excursion-дарын билдирет.

BAM ошол эле сыноо шарттарында ылдамдыкты белгиленген 0–20 m/s модель диапазонунда кармап, кыйла жылмакай ылдамдык эволюциясын жараткан. Бирок авторлор бул салыштырууну эмпирикалык тактыктын артыкчылыгы катары чечмелешпейт. Жыйынтыктар сыналган сандык шарттарда гана BAM чектелүүрөөк жана физикалык жактан көбүрөөк негиздүү сандык ылдамдык жүрүм-турумун түзгөнүн көрсөтөт.

Тыгыздык талаасы кандай өзгөргөн?

BAM үчүн 500 m шакек жолдогу 50 секунддук тыгыздык симуляциясында булак болжол менен \(x=150\,m\) аймагында болжол менен \(0.45\,veh/m\) локалдык максимум жана болжол менен \(x=375\,m\) аймагында болжол менен \(0.15\,veh/m\) минимум билдирет. Талаа жалпысынан жылмакай бойдон калып, айкын үзүлүүлөрдү түзбөйт.

Классикалык PW тыгыздык талаасында болсо булак болжол менен \(150\,m\) аймагында \(0.42\,veh/m\) максимум жана болжол менен \(350\,m\) аймагында \(0.17\,veh/m\) минимум билдирет. Баштапкы жылмакай профиль убакыт өткөн сайын курчураак түзүмдөргө жана артка тараган stop-and-go толкундарына айланат. Изилдөөнүн чечмелөөсү боюнча, сыналган шарттарда PW баштапкы бузулууну күчөтөт, BAM болсо жылмакайыраак толкун таралышын жаратат.

Транспорт агымынын жыйынтыктары

PW үчүн булак 500 m жол жана 6 секунддук агым талдоосунда болжол менен \(150\,m\)де \(4.55\,veh/s\) максимум жана болжол менен \(350\,m\)де \(3.4\,veh/s\) минимум билдирет. Толкундун формасы убакыт өткөн сайын симметриясын жоготуп, деформацияланат.

BAM үчүн агымдын максимуму болжол менен \(125\,m\)де кайра эле \(4.55\,veh/s\), минимуму болсо болжол менен \(375\,m\)де \(3.4\,veh/s\) деп берилген. Агым чоңдуктарынын бири-бирине жакын болушу маанилүү: авторлор жыйынтык бөлүмүндө да белгилегендей, BAMдын сыноолордогу негизги айырмасы “көбүрөөк транспорт агымын” жаратуу эмес, транспорт абалдарынын убакыт жана мейкиндик боюнча өнүгүүсүнүн көбүрөөк көзөмөлдөнүшү.

Робасттуулук эксперименттери кандай түзүлгөн?

Изилдөө андан кийин BAMды узундугу 250, 500 жана 1000 m болгон жабык шакек жолдордо төрт шарт астында 50 секунд бою изилдеген. Бул эксперименттерде мейкиндик кадамы 1 m, убакыт кадамы 0.005 s, \(C_0=20\,m/s\), \(s_0=0.30\) жана \(v_{max}=20\,m/s\) катары сакталган.

D сценарийинде жолдун борборундагы %10 бөлүктө жеткиликтүү максималдуу ылдамдык %25 азайтылып, кайталануучу локалдык өткөрүү жөндөмүнүн чектөөсү түзүлгөн.

RMS бузулуу катышы эмнени көрсөтөт?

Изилдөө баштапкы тыгыздык пертурбациясынын симуляциянын аягында канчасы калганын баалоо үчүн акыркы RMS тыгыздык пертурбациясынын баштапкы RMS маанисине катышын колдонот. Катыш 1ден төмөн болгондо бузулуу жалпысынан басаңдаган, 1ден жогору болгондо болсо акыркы бузулуу баштапкыдан чоң болгон.

A сценарийинде жол узарган сайын акыркы RMS катышы 0.002ден 0.213 жана 0.681ге жогорулаган. Башкача айтканда, кыска жолдо баштапкы бузулуу дээрлик толугу менен басаңдаса, узунураак жол көбүрөөк калдык мейкиндиктик өзгөрмөлүүлүктү сактаган.

B сценарийинде жогорку баштапкы тыгыздык жана күчтүүрөөк пертурбацияга карабастан, чечимдер сыналган үч жол узундугунда тең чектелген бойдон калган. RMS катыштары тиешелүүлүгүнө жараша 0.002, 0.219 жана 0.686.

C сценарийинде узагыраак негизги релаксация убакыт масштабы жана чоңураак standstill аралыгы колдонулган. Булак тексти 250 m үчүн RMS катышын болжол менен 0.000001 деп түшүндүрөт, ал эми 7-таблица муну үч ондук белгиге чейин 0.000 деп көрсөтөт. Бул айырма тегеректөө деңгээлинде; бул жерде таблицадагы маани таблицанын ичинде так сакталган.

Эң күчтүү калдык бузулуу активдүү тар жер болгон D сценарийинде чыккан. 500 m жана 1000 m жолдор үчүн RMS катыштары тиешелүүлүгүнө жараша 1.332 жана 1.923; башкача айтканда, акыркы тыгыздык пертурбациясы баштапкы мааниден чоң. Булак муну колдонулган локалдык өткөрүү жөндөмүнүн чектөөсү жараткан туруктуу мейкиндиктик кайра бөлүштүрүү жана бузулуунун өсүшү катары чечмелейт.

Булак текстинде D сценарийи, 1000 m үчүн максималдуу ылдамдык 6.728 m/s деп жазылса, 7-таблицадагы маани 6.729 m/s. Verianla сандык таблицада булак таблицасынын 6.729 m/s маанисин сактайт жана баяндоо менен таблицанын ортосундагы 0.001 m/s айырманы жашырбайт.

250, 500 жана 1000 метрлик жолдорду салыштыруу

Булак 8–15-сүрөттөрдө төрт сценарийдин 250, 500 жана 1000 m жол узундуктары үчүн акыркы тыгыздык жана ылдамдык профилдерин берет. Жалпы тенденция — кыска жолдор баштапкы бузулууну күчтүүрөөк басаңдатат, узун жолдор болсо көбүрөөк мейкиндиктик өзгөрмөлүүлүктүн калып калышына мүмкүнчүлүк берет. Бул жыйынтык айрыкча RMS катыштарында ачык көрүнөт.

Бирок бул абалды “кыска жолдор реалдуу транспортто дайыма туруктуураак” деп жалпылоого болбойт. Изилдөөдөгү жол узундугунун таасири жабык шакек, периоддук чектик шарт, тандалган BAM параметрлери жана аныкталган баштапкы бузулуулар астында алынган сандык жыйынтык.

Реакция убактысы жана standstill аралыгынын сезгичтиги

Булак 16 жана 17-сүрөттөрдө 500 m шакек жолдо эки параметр топтомун салыштырат. Биринчи учурда \(\tau=0.52\,s\) жана \(l_s=8\,m\); экинчи учурда кыскараак \(0.052\,s\) убакыт масштабы жана 2 m аралык колдонулган. Булак чоңураак убакыт масштабы жана standstill аралыгы шартында тыгыздык толкундары жылмакайыраак болгонун; кичирээк маанилерде болсо күчтүүрөөк ылдамдык жана тыгыздык өзгөрүүлөрү пайда болгонун билдирет.

0.052 s жана 2 m маанилери кадимки адам айдоосун билдирет деп эсептелбеши керек. Макала бул маанилерди ачык түрдө көзөмөлдөнгөн сезгичтик шарттары катары колдонот.

Изилдөө колдогон жыйынтыктар

• BAM сыналган сандык шарттарда классикалык PWге салыштырмалуу чектелүүрөөк ылдамдык маанилерин жана жылмакайыраак транспорт толкун эволюциясын жарата алары колдоого алынат.
• Аралыкка жараша башкарылган каалаган ылдамдык менен тыгыздыкка көз каранды релаксация убакыт масштабын бирге колдонуу моделдин булак мүчөсүн транспорт абалына ыңгайлаштырат.
• Сызыктуу талдоо BAM үчүн ачык туруктуулук чегин берет.
• \(C_0=2.0\), \(4.15\) жана \(20.0\,m/s\) сыноолору теориялык туруктуулук чеги менен шайкеш сандык жүрүм-турум көрсөткөн.
• 250–1000 m робасттуулук эксперименттеринде BAM чечимдери изилденген параметр мейкиндигинде чектелген бойдон калган.
• Активдүү тар жер изилденген шарттардын ичинен эң чоң калдык тыгыздык бузулуусун жараткан.
• BAM менен PWнин транспорт агымынын чоңдуктары окшош болушу мүмкүн; негизги айырма изилдөөдө көбүрөөк абалдардын туруктуулугу жана убакыт-мейкиндик эволюциясы катары көрүнгөн.

Изилдөө колдобогон жыйынтыктар

• BAM реалдуу жол тармактарында PWге караганда сөзсүз такыраак экени көрсөтүлгөн эмес.
• Модель реалдуу транспорт сенсору, loop detector же унаа траекториясынын маалыматтары менен калибрленген эмес.
• Көз карандысыз реалдуу дүйнө маалымат топтомунда божомол тактыгы өлчөнгөн эмес.
• Ар түрдүү унаа узундуктары, унаа түрлөрү жана жеке айдоочу өзгөчөлүктөрү ачык чечилбейт.
• Изилдөөдөгү \(\Delta x\approx1/\rho\) ыкмасы гетерогендүү унаалар аралык аралыктардын баарын чагылдырбайт.
• 0.052 секунддук параметр реалдуу адам айдоочусунун реакция убактысынын натыйжасы эмес.
• 2 m standstill аралыгы жалпы реалдуу транспорт нормасы катары берилбейт.
• Көзөмөлдөнгөн симуляциядагы сандык туруктуулук реалдуу транспорт божомолунун тактыгынын далили эмес.

Изилдөөнүн негизги чектөөлөрү

Авторлор ачык белгилеген эң маанилүү чектөө — эмпирикалык калибрлөөнүн жана текшерүүнүн жоктугу. Азыркы жыйынтыктар жабык шакек жолдо жана алдын ала аныкталган баштапкы пертурбациялар боюнча жүргүзүлгөн көзөмөлдөнгөн сандык эксперименттерге негизделген.

BAM ошондой эле алдыдагы унаалар аралык аралыкты \(\Delta x\approx1/\rho\) ыкмасы менен макроскопиялык тыгыздыктан чыгарат. Бул ыкма салыштырмалуу бир тектүү аралык божомолун камтыйт жана ар түрдүү унаа түрлөрү, унаа узундуктары же жеке айдоочу өзгөчөлүктөрүнөн келип чыккан гетерогендүүлүктү ачык чечпейт.

Келечектеги иштер катары авторлор BAM параметрлерин реалдуу транспорт маалыматтары менен калибрлөөнү; транспорт агымы, ылдамдык, тыгыздык, headway жана сапар убактысы сыяктуу көз карандысыз өлчөөлөр менен текшерүүнү; гетерогендүү транспорт шарттарын, параметр сезгичтигин жана реалдуураак тармак түзүмдөрүн изилдөөнү сунушташат.

Булак жана Ыкма Жөнүндө Эскертүү

Түпнуска аталыш: A macroscopic traffic flow model with gap-controlled driver response and adaptive reaction dynamics

Авторлор жана ирети: Benjamin Simwinga; Aaron Zimba; Mayumbo Nyirenda.

Жооптуу автор: Benjamin Simwinga.

Мекемелер: Benjamin Simwinga жана Aaron Zimba — School of Computing, Technology and Applied Sciences, ZCAS University, Lusaka, Zambia. Mayumbo Nyirenda — Department of Computer Science, University of Zambia, Lusaka, Zambia.

Журнал: Discover Networks.

Басмачы/платформа: Springer Nature / Springer Nature Link.

Том жана макаланын номери: Volume 2, Article 39 (2026).

DOI: 10.1007/s44354-026-00053-y

Жарыяланган күнү: 8 Октябрь 2026.

Процесс даталары: Алынды: 10 Август 2026; кайра каралды: 19 Сентябрь 2026; кабыл алынды: 28 Сентябрь 2026.

Булактын түрү: Research article; математикалык моделди иштеп чыгуу, сызыктуу туруктуулук талдоосу жана көзөмөлдөнгөн сандык симуляция изилдөөсү.

Лицензия: Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0). Лицензия ачык жеткиликтүүлүктү камсыздайт, бирок ыңгайлаштырылган материалды бөлүшүүгө уруксат бербейт; ошондуктан Verianla баяндоосу булак текстинен жана фигураларынан көз карандысыз илимий маалымат архитектурасы менен даярдалган.

Авторлордун салымы: Булакта BS макаланы жазганы; AZ жана MN карап чыгуу жана багыт берүү камсыз кылганы айтылат.

Каржылоо: Изилдөө каржылоо албаганын билдирет.

Маалыматтардын жеткиликтүүлүгү: Изилдөө учурунда түзүлгөн жана/же талданган маалымат топтомдорун негиздүү суроо-талап боюнча жооптуу автордон алууга болору билдирилет.

Кызыкчылыктардын кагылышы: Авторлор кызыкчылыктардын кагылышы жок экенин билдиришет.

Этикалык бекитүү жана катышуучунун макулдугу: Колдонулбайт деп билдирилген.

Эмпирикалык текшерүү абалы: Жок. BAM али байкалган транспорт маалыматтары менен калибрленген эмес жана көз карандысыз талаа өлчөөлөрү менен текшерилген эмес.

Статистикалык тыянак: Изилдөөдө p мааниси, ишеним аралыгы, клиникалык/эксперименттик тандалма статистикасы же статистикалык гипотеза тести берилген эмес. Жыйынтыктар аналитикалык туруктуулук шарттары жана сандык симуляция натыйжалары.

Булак ичиндеги техникалык дал келбестиктер: 5-таблицада \(\Delta t=0.005\,s\) берилгенине карабай, түшүндүрмө тексттин бир жеринде 0.006 s деп жазылган. 5-таблицадагы 0.005 s мааниси 1200 кадам жана 6 секунддук жалпы убакыт менен өз ичинде шайкеш; бирок Verianla ушул себептен булак текстиндеги 0.006 s туюнтмасын унчукпай өчүргөн эмес. 5-таблицадагы 0.30 veh/m “орточо баштапкы тыгыздык” менен Теңдеме (49)дагы 0.01/0.30/0.10 бөлүктүк баштапкы профилдин ортосундагы байланыш да булакта кеңири түшүндүрүлгөн эмес. Теңдемелердин кайчылаш шилтемелеринде жана айрым вектордук көрсөтүлүштөрдө да терүү/нотация дал келбестиктери бар.