Латтис Больцманы Арга

Чөлөөт нэвтэрхий толь — Википедиагаас
Jump to navigation Jump to search

Латтис Больцманы арга (LBM эсвэл ЛБА)' (эсвэл Дулааны Латтис Больцманы арга (TLBM)) нь тооцон бодох шингэний динамик (CFD) дахь шингэний симуляцийг гүйцэтгэх нэг арга юм. Навиер-Стокесийн тэгшитгэлийг шийдэхийн оронд тухайлсан (дискрет) Больцманы тэгшитгэл нь Ватнагар-Гросс-Крүүк (BGK-ВГК) гэх мэт мөргөлдөөний загваруудын тусламжтайгаар Ньютоны шингэний урсгалыг симуляцлах боломжтой. Тодорхой тооны эгэл хэсгүүдийн хувьд урсгал болон мөргөлдөөний процессыг симуляци хийснээр бодит шингэний эгэл хэсгүүдийн хоорондын харилцан үйлчлэл нь массаас илүүтэйгээр урсгалын шинж чанарт бичил-орчны зунгааралт нөлөөлдөг гэдгийг батладаг.

Алгоритм[засварлах | edit source]

ЛБА бол комплекс шингэний системд харьцангуй шинэ тооцооллын арга бөгөөд тооцон бодох шинжлэх ухаанд ялангуяа физик, шингэний динамикийн салбарт судлаачдийн анхаарлыг маш их татаж байна. Хадгалагдах хуулийн Макро-орчны шинж чанаруудыг (ж.н масс, момент, энерги) тоон аргаар шийдвэрлэдэг уламжлалт ТБШД-ын аргуудтай (хэлхээст арга) адилгүйгээр ЛБА нь шингэнийг төсөөллийн эгэл хэсгүүдээс бүрдсэн мөнөөх эгэл хэсгүүд нь дискрет сүлжээ (Латтис гэдэг нь сүлжээ юм) бүхий хэлхээсийн дагуу дэс дараалсан түгэлт болон мөргөлдөөнийг бий болгоно гэж загварчилдаг. Энэхүү хэсгийн чухал байдал, орчны динамикийн ачаар ЛБА нь бусад нийтлэг хэрэглэгдэх ШДТБАргуудаас илүү гарах хэд хэдэн давуу талийг бий болгох ба ялангуяа комплекс хязгаартай бодлогууд, макро орчны харилцан үйлчлэлийн хамаарал, алгоритмыг параллелчлах нөхцөлд маш тохиромжтой давуу талуудыг бий болгоно. Сүлжээний буюу Латтис Больцманы тэгшитгэлийн өөр нэг илрэл нь тухайлсан-хурдны Больцманы тэгшитгэл юм. Тухайлсан тэгшитгэл гэдэг нь дискрет орчинд бичигдсэн эсвэл хөрвүүлэгдсэн тэгшитгэл гэж ойлгоно. Тухайн дифференциал тэгшитгэлийн системийг шийдэх тоон арга нь төсөөллийн эгэл хэсгүүдийн түгэлт ба мөргөлдөөнийг харуулсан дискрет дүрслэлийг бий болгоно.

Комьпютер алгоритмд, мөргөлдөөн болон урсгалын алхамууд нь дараах байдлаар илэрхийлэгдэнэ:

Мөргөлдөх алхам:

Урсах алхам:

Энд i моментийн чиглэлийг заах ба түгэлт гэдэг нь төсөөллийн эгэл хэсгүүдийн урсах алхамыг илэрхийлж байгаа юм.

Хийн сүлжээний автоматын (LGA-ХСА) аргаас (ЛБА) хөгжсөн нь[засварлах | edit source]

ЛБА нь орон зай хугацаанд төсөөллийн молекулын динамикийг хялбарчилан авч үздэг хийн сүлжээний автомат (LGA) гэдэг аргаас үүсэлтэй ба хийн автоматад эгэл хэсгийн хурд нь үргэлж дискрет байдаг. Тухайлбал 2 хэмжээст FHP загварт сүлжээний зангилаа бүр гурвалжин байдлаар хөрш сүлжээний зангилаатайгаа 6-н шилбээр холбогдох ба ямарч эгэл хэсэг нь уг шилбээр дамжин өгөгдлөн сүлжээний хурдаар хөдлөнө. Хугацааны интервалын дараа, эгэл хэсэг бүр өгөгдсөн чиглэлээрээ хөрш зангилааруу нүүх ба энэ процессийг түгэлтийн буюу урсгалын шат гэж нэрлэдэг. Нэгээс илүү олон эгэл хэсгүүд өөр өөр чиглэлээс ижил зангилаанд ирэх тохиолдолд хоорондоо мөргөлдөж чиглэл болон хурдаа тусгай мөргөлдөөний дүрмийн дагуу өөрчилнө. Урсгалын болон мөргөлдөх шатууд тусгай хувилбартай байдаг. Хамгийн тохиромжтой мөргөлдөөний дүрэм нь мөргөлдөөний өмнөх ба дараах үеийн масс эгэл хэсгийн тоо, момент (хөдөлгөөний тоо хэмжээ), энерги зэргийг хадгалсан байх ёстой. ХСА нь гидродинамикийн симуляци хийхэд зарим нэг засагдашгүй сул талуудтай. Жишээ нь: хурд ихтэй урсгалд Галилейн инвариант алдагдах, сүлжээний хэмжээнээс хамаарч статистик шуугиан гарах эсвэл Рейнольдсын тоо хязгаарлагдах гэх мэт. Хэдий тийм боловч ХСА нь реактив тархалт болон молекул динамикийн загварчлалд хүрэхэд маш тохиромжтой шилжүүлэгддэг.

ХСА-аас ЛБА руу шилжих гол үндэслэл нь сүлжээний чиглэл дахь Боолиан эгэл хэсгийн тоо ба түүний олонлогийн дунджийг нягтын түгэлтийн функцаар сольж статистик шуугианыг үгүй хийх хүсэл байсан юм. Энэ солилтыг хийснээр дискрет мөргөлдөөний дүрэм нь мөргөлдөөний оператор гэж нэрлэгдэх тасралтгүйн функцаар бас солигдсон юм. ЛБА-ын хөгжлийн үе дахь чухал нэг хялбарчилгаа бол мөргөлдөөний операторыг Ватнагар-Гросс-Крүүк (BGK-ВГК) гэх сулралтын (заримдаа тайвшрах) илэрхийлэлээр ойролцоолох явдал байв. Энэхүү сүлжээний ВГК (LBGK) загвар нь симуляцийг илүү үр ашигтай болгох ба шилжүүлгийн коэффициентийн уян чанарыг бий болгоно. Өөрөөр хэлбэл ЛБА нь тасралтгүй Больцманы тэгшитгэлийн онцгой тухайлсан хэлбэрийг авч үздэг гэсэн үг юм. Чапман-Энскогийн онолоор ЛБА-ын алгоритмаас Навиер-Стокесийн тэгшитгэл болон тасралгүйн тэгшитгэлийг сэргээж болох ба даралтын орон нь уламжлалт ТБШД-ын аргууд шиг нэмэлт тэгшитгэл эсвэл Пауссоны тэгшитгэл тэгшитгэлийг шийдэлгүйгээр нягтын хуваарилалтаас шууд олдох боломжтой байдаг.

Сүлжээ ба DnQm ангилал[засварлах | edit source]

Латтис Больцманы загвар нь куб болон гурвалжин хэлбэртэй эсвэл хослосон байхаар хэд хэдэн төрлийн сүлжээ болон дискрет түгэлтийн функцийн хөдөлгөөнтэй ба хөдөлгөөнгүй эгэл хэсгүүдийг ашиглан ажилладаг. Сүлжээний өөр байдлаас бий болж байгаа аргуудыг ангилах хамгийн зөв зам нь DnQm схем юм. Энд "Dn" нь "n хэмжээстэй" гэдгийг илтгэж байхад "Qm" нь "m хурдтай" гэдгийг илтгэнэ. Жишээлбэл, D3Q15 бол хөдөлгөөнгүй эгэл хэсгүүд бүхий куб торон дээрх гурван хэмжээст Латтис Больцманы загварыг илэрхийлж байна. Зангилаа бүр нь кристал хэлбэрийг бий болгох ба гадаргаа хуваалцаж байгаа 6, хүрээгээр нийлж байгаа 8-н хөршдөө тус бүр эгэл хэсгээ түгээх боломжтой юм.[1] (D3Q15 загвар нь хүрээгээр нийлж байгаа 12 хөршрүүгээ эгэл хэсгүүдийг түгээх боломж байхгүй ба хэрэв боломжтой болчихвол энэ нь "D3Q27" загвар болно.)

Симуляцийг хийхийн тулд орон зай ба хугацаа гэх мэт бодит тоо хэмжээнүүд сүлжээний нэгжрүү шилжүүлэгдсэн байх хэрэгтэй. Рейнольдсын тоо гэх мэт хэмжээсийн бус хэмжигдэхүүнүүд хэвээр байж болно.

Сүлжээний нэгжийн шилжүүлэг[засварлах | edit source]

Ихэнх Латтис Больцманы симуляцид нь сүлжээний хэмжээсийн үндсэн нэгж байх ба хэрэв бодлогийн хүрээний урт -ийн турш тооны шилбэ оршино гэж үзвэл хэмжээсийн нэгж нь хялбар байдлаар гэж тодорхойлогдоно. Латтис Больцманы симуляци дахь хурд нь ерөнхийдөө дууны хурдны илэрхийллээр өгөгддөг. Иймд дискрет хугацааны нэгж нь гэж илэрхийлэгдэх ба үүний хуваагч нь дууны физик хурд юм.[2]

Бичил хэмжээстэй урсгалд (механик ба гидравликт сүвэрхэг орчин дундуурх урсгал) бодит дууны хурдыг ашиглах нь үл зөвшөөрөх жижиг хугацааны алхамыг өгдөг. Иймд сүлжээний Махын тоо гэж бий болох ба энэ нь бодит Махын тооноос их утгатай байна. Мөн Рейнольдсын тоог хэвээр хадгалахын тулд зунгааралтыг нэмэгдүүлж тэнцвэржүүлдэг.[3]

Холимог шингэний симуляци[засварлах | edit source]

Хөдлөх болон хэв гажих харилцах хязгаартай учир олон фазат/холимог урсгалын симуляци нь ТБШДинамик дахь уламжлалт аргуудын томоохон шалгуур байсаар ирсэн. Илүү нарийвчилбал ялгаатай фазууд (шингэн ба уур) эсвэл бүрэлдэхүүнүүд (тухайлбал тос ба ус) хоорондын харилцах хязгаар нь шингэний молекулуудын хооронд тодорхой харилцан үйлчлэлийг бий болгодог. Иймд макро-орчны шинжийг агуулах Навиер-Стокесийн тэгшитгэлрүү энэ мэт харилцан үйлчлэлийн микро-орчны шинж чанаруудыг боловсруулж хийх нь хүндрэлтэй байдаг. Гэхдээ ЛБА-д хэсгийн кинетик нь мөргөлдөөний операторыг засварласанаар суурь болсон микро-орчны шинж чанаруудтай холбогдох аргыг тохиромжтой харьцангуй хялбар болгож өгдөг. Хэд хэдэн ЛБА-ын олон фазат/холимог урсгалын загварууд хөгжсөөр байна. Эдгээрт фазийн ялгарал нь эгэл хэсгийн динамикаас автоматаар бий болох ба уламжлалт хэлхээст аргууд шиг ямар нэгэн тусгай арга шаардлагагүй. Олон фазат/холимог ЛБА нь өөр өөр комплекс шингэний системд тухайлбал харилцах хязгаарын тогтворгүйжилттэй урсгал, хөөс/дусалын динамик, хатуу гадаргуугын норолт,шингэний дотоод үе хоорондын гулсалт, дуслын электрогидродинамикийн деформаци гэх мэтэд амжилттай хэрэглэгдэж байна.

Бага Махын тооны горимтой үед нягтын өөрчлөлтийн ач холбогдлыг авч үзсэн хийн холимогийн шаталт дахь Латтис Больцманы загвар нь саяхан дэвшүүлэгдээд байна.[4]

Эндээс үзвэл дан ЛБА нь томоохон физик оронг (тохиромжит аргатай харьцуулбал) авч үзэж болох ба реактив хийн холимогийн симуляци нь том хэмжээний шаталтын механизмыг нарийвчлан авч үзэж байгаа мэт санах ойн хэрэгцээний хувьд зарим нэг нэмэлт сайжруулалт хэрэгтэй гэдгийг харуулдаг юм. Эдгээр асуудлууд нь "загварын системчилсэн өөрчлөлтийн техник"-т дахин ангилагдаж судлагдаж, боловсруулагдаж байна.[5][6][7]

Дулааны Латтис Больцманы арга[засварлах | edit source]

Саяхан (2009 онд), Дулааны Латтис-Больцманы арга (TLBM - ДЛБА) нь Олон-хурдат аргачлал /multi-speed approach/,[8] Идэвхгүй скаляр аргачлал /passive scalar approach/,[9] Дулааны энергийн түгэлтийн аргачлал /thermal energy distribution/[10] гэх гурван ангилалын сүүлчийнх нь болсон юм.

Хязгаарлалтууд[засварлах | edit source]

Комплекс шингэний динамикийн симуляцид ЛБА-ын хэрэглээний өсөлтыг үл анхаарвал энэхүү шинэ аргачлалд зарим хязгаарлалтууд байдаг юм. Одоогийн байдлаар аэродинамик дахь өндөр Махын тоотой урсгалыг симуляцилахад хүндрэлтэй байгаа ба тохиромжит термо-гидродинамикийн схем олдоогүй байна. Гэвч Навиер-Стокесийн тэгшитгэл дээр суурилсан ТБШД-ын аргуудтай ЛБА нь амжилттай хосолж дулааны шилжилттэй (хатуу биеын чанарт суурилсан дулааны дамжилт, нэвчилт ба ялгаралт) симуляциуудыг чадварлаг шийдвэрлэж байна. Олон фазат/холимог урсгалын загварын хувьд харилцах үеийн зузаан нь ихэнхдээ том байх ба энэ орчмийн нягтын харьцаа нь бодит шингэнтэй харьцуулахад бага байдаг. Саяхан энэ асуудал Юуан болон Счаефер нараар шийдэврлэгдлээ. Тэд Шан ба Чен нарын загвар дээр сайжруулалт хийж энэ үр дүнд хүрсэн ба төлөвийн тэгшитгэлийг хялбарчилсанаар 1000:1 -ын нягтын харьцаанд хүрэх боломжтой гэж харуулжээ.

Гэсэн хэдий ч, энэ аргын сүүлийн хорин жилийн хурдан дэвшил, өргөн хэрэглээ нь энэ аргыг микро-урсгалыг багтаагаад тооцон бодох шинжлэх ухаанд ялангуяа физикийн салбарт чадамжтай гэдгийг баталж байна. ЛБА нь өндөр Кнүдсэний тоотой урсгалд хамгийн нарийвчлалтай үр дүнг өгдөг нь туршилтаар батлагдсан.

Тухайлсан ЛБТэгшитгэлээс Навиер-Стокесийн тэгшитгэлийг гаргах[засварлах | edit source]

Тухайлсан Латтис Больцманы тэгшитгэлээс (мөргөлдөөний операторыг хэрэглэсэн ЛВГК (LBGK) тэгшитгэл гэж нэрлэнэ) эхэлье. Бид эхлээд ЛБТэгшитгэлийн зүүн талыг -р эрэмбийн Тэйлорын цуваанд задлая. Хялбарчилсан эрэмбийн Тэйлорын цувааг сонгосон бол ЛБТэгшитгэл нь утгаа алдана. - р эрэмбийн Тэйлорын цуваанд задлахад баруун талаас 0 болон нэгдүгээр эрэмбийн уламжлалтай хэсэг нь орхигдож, тэйлорын цувааны 1 ба 2-р эрэмбийн уламжлал болон мөргөлдөөний операторууд тус тус үлдээнэ.

Хялбарчлахын тулд - ыг гэж бичье. Тэйлорын цувааг хөнгхөн хялбарчилж диад хоорондын үржвэрийг тодорхойлох цэг ":" -ээр тэмдэглэе.

Тэнцвэрт ба тэнцвэрт бус бүрдүүлэгчрүү эгэл хэсгийн түгэлтийн функцыг өргөтгөж мөн Күдсэний тоо бүхий Чапман-Энсогийн задаргааг хэрэглэн Тэйлорын задаргаа бүхий ЛБТ нь тодорхой тасралтгүй тэгшитгэлийг олохын тулд Кудсэний тооны хувьд өөр өөр утгын хэмжээрүү задрах болно.

Тэнцвэрт ба тэнцвэрт бус түгэлт нь макро-түвшний хэмжигдэхүүнүүдтэй дараах байдлын холбоог хангах ёстой. Эгэл хэсгийн хуваарилалт нь эгэл хэсгээр макро-түвшинг хэмжихийн тулд 'зөв хэлбэрт' байх ёстой.

Чэпман-Энскогийн задаргаа нь:

.

Тэйлорын цуваанд тэнцвэрт ба тэнцвэрт бус түгэлтийн задаргааг орлуулж, өөр өөр эрэмбэд ялгасанаар тасралтгүй орчны тэгшитгэл илэрхийлэгдэж эхэлнэ.

Эрэмбэ -ийн хувьд:

Эрэмби -ийн хувьд:

Хоёр дахь тэгшитгэлд зарим хувиргалт хийж эхний тэгшитгэлд орлуулбал дараах байдалтай болно.

Дээрхээс макро-орчны шинж чанар болон эгэл хэсгийн хоорондын хамаарлыг хэрэглэн масс болон моментийн тэгшитгэлийг гарган авч болно.

Моментийн цүлхэлтийн тенсор нь дараах хэлбэртэй байна.

Үүнд: бол -ын (ж.н ) бүх бүрдүүлэгчийн нийлбэрийн квадратын хувьд товч тэмдэглэгээ ба Навиер-Стокесийн тэгшитгэлтэй харьцуулж болохуйц хоёрдугаар эрэмбийн тэнцвэржүүлэгч эгэл хэсгийн түгэлтийн тэгшитгэл нь:

.

Тэнцвэржүүлэгч буюу тэнцвэрт түгэлт нь зөвхөн бага хэмжээний Махийн тоотой буюу бага хурдтай үед л хүчинтэй. Тэнцвэрт түгэлтийг цүлхэлтийн тенсорлуу оруулж ирсэнээр:

Эцэст нь Навиер-Стокесийн тэгшитгэлийг нягтын маш бага өөрчлөлттэй гэсэн нөхцөл дор гарган авч байна.

Энэ гаргалгаа бол Чен болон Доолений ажил юм.[11]

Симульци дахь математик тэгшитгэлүүд[засварлах | edit source]

Тасралтгүй Больцманы тэгшитгэл нь дан эгэл хэсгийн магадлалт хуваарилалтын /түгэлт/ функц -ын хувьд үнэлэмжийн тэгшитгэл нь бөгөөд дотоод энергийн нягтын хуваарилалтын функцууд (Хе нар.) нь тус бүр:

Үүнд: нь дараах байдлаар -тэй холбогдоно.

бол гадаад хүч, бол мөргөлдөөний интеграл ба (зарим эх сурвалжид гэж тэмдэглэгдэнэ) бол микро түвшний хурд юм. Гадаад хүч температураас хамааралтай гадаад хүч болох -тэй дараах байдлаар холбогдоно. Хүчний харилцааны загварын нэг шалгалт нь -ийн хувьд Раелейг-Венардын нэвчилт юм.

Нягт , хурд болон температур гэх мэт макро-орчний хэмжигдэхүүнүүд нь нягтын түгэлтийн функцын момент шиг тодорхойлогдох боломжтой ба илэрхийлбэл:

Латтис Больцманы арга нь сүлжээрүү хязгаарлах хэмжээсийг оруулж мөн хурдны оронг микро орчны хурдны тухайлсан бүрдэлрүү (ж.н ) шилжүүлсэнээр дээрх тэгшитгэлийг тухайлж дискрет тэгшитгэл болгоно. Тухайлбал D2Q9, D3Q15, болон D3Q19 дахь микро-орчны хурдууд нь дараах байдлаар илэрхийлэгдэнэ.

Массын нягт, дотоод энергийн нягтын хувьд дан фазад тухайлагдсан Больцманы тэгшитгэл нь:

Мөргөлдөөний оператор нь зарим нэг нөхцөлд тухайлбал хадгалагдах хуулийг хангаж байхын тулд үргэлж ВГК мөргөлдөөний оператороор ойролцоологдоно.

Мөргөлдөөний операторт байх нь дискрет ба үүнийг тэнцвэрт эгэл хэсгийн магадлалт түгэлтийн функ гэж үзнэ. D2Q9 ба D3Q19-ын хувьд операторыг тасралтгүй ба дискрет хэлбэрт үл шахагдах шингэний хувьд харуулсан ба үүнд D, R ба T нь хэмжээснүүд, хийн тогтмол болон абсольют температурууд юм. Тасралтгүйгээс дискрет хэлбэрийн хувьд тухайн уламжлал нь 2-р эрэмбийн тохиролтойгоор /нарийвчлал/ хийгдэв.

гэж авч үзвэл эцсийн дүнд:

Дан бүрэлдэхүүнт урсгалын хувьд энэ төрлийн ажлууд маш их хийгдсэн байдаг ба одоо Дулааны ЛБА-ын тухай авч үзье. Олон бүрэлдэхүүнт/олон фазат ДЛБА нь нэг бүрдэлт урсгалаас илүү хялбар ашиглагдахуйц, сонирхол татдаг юм. Орчин үеийн судалгаатай хөл нийлүүлэхийн тулд элементтэй системийн бүх бүрдүүлэгчийг (ж.н сүвэхрхэг орчны хана, олон төрлийн шингэний хий холимог) -аар тодорхойлох явдал юм.

Тайвшрах буюу сулралтын параметр нь кинематик зунгааралт - тай дараах холбоогоор холбогдоно.

-ын моментууд нь хадгалагдсан орчны хэмжигдэхүүнүүдийг өгнө. Энэ үед нягт нь:

Жинлэсэн дундаж хурд болон орчны момент /хөдөлгөөний тоо хэмжээ/ нь дараах байдлаар илэрхийлэгдэнэ.

Тэнцвэрт хурд -ийн хувьд өгөгдсөн дээрх тэгшитгэлд гишүүн нь бүрэлдэхүүн бусад бүрэлдэхүүнтэйгээ харилцан үйлчлэлцэх хүч юм. Энэ нь шингэн шингэнтэй, хий шингэнтэй, шингэн хатуу биетэй хэрхэн харилцахыг тодорхойлдог ердийн тохируулах параметр боловч одоог хүртэл хэлэлцээний гол сэдэв болж байна. Франк ба бусад эрдэмтэд энэ хүчний хувьд нэгэн загвар жагсаалт гаргаж ирсэн. Түгээмэл ашиглагддаг гаргалгаа нь Гүнстенсений хромодинамикийн загвар ба энэ нь шингэн уурын систем, бинарын шингэнд маш тохиромжтой Свифтийн чөлөөт энерги дээр суурилсан аргачлал юм. Үүнд дотоод молекулын харилцаанд суурилсан загвар, Инамурогийн арга, Лее ба Линий аргачлалч[12] мөн хамаарна.

-ийн ерөнхий хэлбэр нь доорх байдлаар хэд хэдэн эрдэмтэдээр өгөдсөн байдаг.[13][14]

бол үр ашигтай масс ба нь хөрш хэсэгтэйгээ - ээр харилцах дотоод эгэл хэсгийг харуулах Греений функц юм. гэсэн нөхцөл хангаж байвал болох ба эгэл хэсгүүд хоорондоо түлхэлцэх хүчийг илэрхийлнэ. D2Q9 болон D3Q19 хувьд энэ нь:

Шан ба Чен нараар дэвшүүлэгдсэн идэвхит масс нь дан ба олон бүрэлдэхүүнт системд үр ашигт масс гэдгээрээ доорхи байдлаар хэрэглэгдэнэ. Төлөвийн тэгшитгэл нь дан болон олон фазат урсгалын нөхцөлд өгөгдөх боломжтой.

Цаашид ба нь төлөвийн тэгшитгэл (EOS-ТТ) - ийн системрүү орсоноор тохируулах чөлөөт тогтмол шиг харагддаг боловч болон гэх мэт критик нөхцөл шиг гидродинамикийн харилцан холбоог хангаж байх ёстой. Төлөвийн тэгшитгэлд бол D2Q9 ба D3Q19 хувьд 3.0 байхад D3Q15[15] хувьд 10.0-тэй тэнцүү байна.

Олон фазат урсгалын симуляцийг илүү тохиролтой буюу нарийвчлалтай байлгахын тулд үр ашигт массын нягт нь өөрчлөгдсөн байх хэрэгтэй гэдгийг Юан ба Счаефер[16] нар харуулсан юм. Тэд өөрсдийн үр дүнг Шан ба Чен (SC), Карнахан-Старлинг (C–S), Ван дер Ваал (vdW), Редлич-Квонг (R–K), Редлич-Квонг Соаве (RKS), болон Пенг-Робинсон (P–R) нарын төлөвийн тэгшитгэлтэй харьцуулсан юм. Тэдний үр дүнд SC ТТ нь маш сул байсан ба C–S, P–R, R–K, болон RKS ТТэгшитгэлүүд нь олон ба дан фазат урсгалд тохиромжтой илүү тохиролыг өгч байсан байна.

Алдарт изотерм ЛБА-ын хувьд эдгээр нь зөвхөн тоо хэмжээнүүдийг л хадгална. Дулааны загвар нь мөн энергийг хадгалж байх ёстой учир нэмэлт хадгалагдах тоо хэмжээ шаардлагатай болно.

Програм хангамжууд[засварлах | edit source]

Нээлттэй үүсвэр / чөлөөт програмууд[засварлах | edit source]

  • Advanced Simulation Library: График Процессийн нэгж (OpenCL) бүхий Нээлттэй ЛБА-ны код (AGPL). Өөр өөр сүлжээний хувьд хүчинтэй.
  • Taxila LBM: Лос Аламос Үндэсний лаборатороос гаргасан 2 ба 3н хэмжээст параллелчилсан нээлттэй үүсвэр (LGPL) бүхий Фортран код. Олон фазат урсгалд хчүинтэй.
  • LIMBES: Гаусс-Хермитийн квадратур дээр суурилсан 2 хэмжээст нээлттэй үүсвэр (GPL) код, Фортран 90, параллелчилсан (openmp)
  • OpenLB: C++, параллелчилсан ЛБА дээр суурилсан нээрттэй үүсвэр (GPLv2) бүхий сан
  • Palabos: Нээлттэй үүсвэр (AGPL) параллел C++ код, дулааны, олон фазат, чөлөөт гадаргуутай урсгалд хүчинтэй
  • Sailfish: График Процессийн нэгж (CUDA/OpenCL) бүхий Нээлттэй ЛБА-ны код (LGPL). Өөр өөр сүлжээний хувьд хүчинтэй.
  • Blender: ЛБА-ыг хэрэглэсэн 3н хэмжээст загварчлагч (http://wiki.blender.org/index.php/Doc:2.6/Manual/Physics/Fluid/Technical_details)

олон фазат болон турбулент урсгалт хүчинтэй.

  • El'Beem: ЛБА-ыг хэрэглэсэн чөлөөт ТБШД-ын код (GPL)
  • J-Lattice-Boltzmann: Жава апплеттэй ЛБА-ын хамаарал.
  • C examples: C (GPL) хэл дээрх зарим ЛБА-ыг ашигласан бодлогын жишээнүүд.
  • LBM-C: Нээлттэй (GPLv2) CUDA код
  • waLBerla: Хэрэглээ өргөнтэй ЛБА
  • elbe: Бүтээмж өндөртэй ЛБА

Чөлөөт програм[засварлах | edit source]

  • LBSim: C++ хэл дээр бичигдсэн үнэгүй програм

Арилжааны програм[засварлах | edit source]

  • LaBS: LaBS эсвэл Даттис Больцманы Солвер гэдэг нь шахагдах шингэнийг загварчилах арилжааны ТБШД-ын код юм.
  • LBHydra: LBHydra бол Миннесотагийн Их сургуулиас гаргасан 2 ба 3н хэмжээст турбулент, олон ба дан фазат урсгалыг загварчлах өргөтгөж болон Латтис Больцманы симулятор юм.
  • PowerFLOW: Экса Корпорацийн ТБШД-ын код.
  • FlowKit: Нээлттэй үүсвэр болох Palabos програмын мэргэжлийн туслах програм.
  • XFlow: Некст Лимит Технологи компанийн бүтээл ЛБА-д суурилсан ТБШД-ын хамгийн сүүлчийн үеийн код.
  • ultraFluidX: Мульт-GPU -ээр дэмжигдсэн ЛБА дээр суурилсан ТБШД-ын код, FluiDyna GmbH FluiDyna.

Цаашид унших[засварлах | edit source]

Тэмдэглэл[засварлах | edit source]

  1. Succi, p 68
  2. Succi, Appendix D (p. 261-262)
  3. Succi, chapter 8.3, p. 117-119
  4. (2012) "Lattice Boltzmann model for reactive flow simulations". EPL 98. DOI:10.1209/0295-5075/98/34001.
  5. (2010) "Coupling of the model reduction technique with the lattice Boltzmann method for combustion simulations". Combust. Flame 157: 1833-1849.
  6. (2012) "Efficient simulations of detailed combustion fields via the lattice Boltzmann method". International Journal of Numerical Methods for Heat & Fluid Flow 21.
  7. (2009) "Combustion simulation via lattice Boltzmann and reduced chemical kinetics". Journal of Statistical Mechanics: Theory and Experiment. DOI:10.1088/1742-5468/2009/06/P06013.
  8. McNamara, G., Garcia, A., and Alder, B., "A hydrodynamically correct thermal lattice boltzmann model", Journal of Statistical Physics, vol. 87, no. 5, pp. 1111-1121, 1997.
  9. Shan, X., "Simulation of rayleigh-b'enard convection using a lattice boltzmann method", Physical Review E, vol. 55, pp. 2780-2788, The American Physical Society, 1997.
  10. He, X., Chen, S., and Doolen, G.D., "A novel thermal model for the lattice boltzmann method in incompressible limit", Journal of Computational Physics, vol. 146, pp. 282-300, 1998.
  11. Chen, S., and Doolen, G.D., "Lattice Boltzmann Method for Fluid Flows", Annual Review of Fluid Mechanics, vol. 30, pp. 329-364, 1998.
  12. Frank, X., Almeida, G., Perre, P., "Multiphase flow in the vascular system of wood: From microscopic exploration to 3-D Lattice Boltzmann experiments", International Journal of Multiphase Flow, vol. 36, pp. 599-607, 2010.
  13. Yuan, P., Schaefer, L., "Equations of State in a lattice Boltzmann model", Physics of Fluids, vol. 18, 2006.
  14. Harting, J., Chin, J., Maddalena, V., Coveney, P., "Large-scale lattice Boltzmann simulations of complex fluids: advances through the advent of computational Grids", Philosophical Transactions of the Royal Society A, vol. 363, pp. 1895–1915 2005.
  15. Yuan, P., Schaefer, L., "A Thermal Lattice Boltzmann Two-Phase Flow Model and its Application to Heat Transfer Problems-Part 1. Theoretical Foundation", Journal of Fluid Engineering 142-150, vol. 128, 2006.
  16. Yuan, P., Schaefer, L., "Equations of State in a lattice Boltzmann model", Physics of Fluids, vol. 18, 2006.

Гадаад линк[засварлах | edit source]