6.6.2: Безперервність осаду
Для отримання концентрації осаду необхідно вирішити рівняння балансу маси для осаду. Загальна заява про збереження осаду говорить:
∂c∂t⏟change insedimentconcentration+∂uc∂x+∂vc∂y⏟net import of sedimentby the horizontalfluid velocity+∂wc∂z⏟net upward transportof sediment by thevertical fluid velocity−∂wsc∂z⏟net downward transportwith the fall velocity=0
або використовуючи рівняння неперервності для рідини:
∂c∂t+u∂c∂x+v∂c∂y+w∂c∂z−∂wsc∂z=0
Зверніть увагу, що в Eqs. ???а??? швидкість і концентрація - це сумарні сигнали, що складаються із середньої, коливальної та турбулентної частини.
Горизонтальні адвективні терміни часто (але не завжди) менші за вертикальні адвективні терміни (див. Інтермеццо 6.6). Давайте для простоти знехтуємо горизонтальними адвективними термінами Eq. ???так що нам залишилося:
∂c∂t+∂wc∂z−∂wsc∂z=0
Аналогічно дрейфу води (Дрейф Стокса, Сект. 5.5.1), також спостерігається дрейф зваженого осаду. Під гребенем хвилі концентрація зваженого осаду розтягується, тоді як під жолобом вона стискається. Дрейф Стокса пояснювався більш високим рівнем води під час прямого орбітального руху, ніж під час руху назад. Аналогічно відбувається дрейф підвішеного осаду. Цей дрейф враховується лише тоді, коли горизонтальні адвективні члени включені в рівняння адвекції-дифузії.
Як було сказано раніше, швидкість і концентрація складаються з середньої, коливальної і турбулентної частини. Значить,w=W+˜w+w′ іc=C+˜c+c′. Ми, звичайно, не збираємося вирішувати турбулентний рух і хотіли б усереднити цей рух. Це називається усереднення Рейнольдса. Якщо ми усереднюємо Eq. ???за шкалою турбулентності більшість термінів з турбулентними коливаннями усереднюють, крім одного терміну, а саме∂⟨c′,w′⟩/∂z. Тут дужки позначають усереднення за часовою шкалою турбулентності і⟨c′,w′⟩ є (висхідним) потоком опадів по турбулентності. Крім того, усереднену Рейнольдсом вертикальну швидкість води можна вважати незначною порівняно зі швидкістю падіння осаду. При спрощенні постійної швидкості падіння (в реальності швидкість падіння залежить від концентрації, див. 6.2.3) ми отримуємо таку часто використовувану форму рівняння адвекції-дифузії:
∂c∂t−ws∂c∂z⏟sediment net going downwardwith its fall velocity+∂⟨c′,w′⟩∂z⏟sediment net going upwardwith fluid turbulence=0
Зверніть увагу, що в цьому рівнянні концентраціяc тепер позначає турбулентність усередненої концентрації:c=C+˜c.
Для моделювання потоку осаду через турбулентність ми робимо аналогічне припущення про висхідний транспорт за рахунок турбулентної дифузії, як і для рідини (див. Наприклад, Eqs. 5.5.5.14 і 5.5.5.17):
−⟨c′,w′⟩=vt,sδcδz
в якійvt,s турбулентна дифузійність осадової маси вm2/s і теперc визначається як об'ємна концентрація (див. 6.2.2). Цей висхідний транспорт за допомогою турбулентної дифузії можна зрозуміти наступним чином. Турбулентний обмін гарантує, що посилка рідини, завантаженої осадом, піднімається вгору до рівня з меншою концентрацією осаду. Посилка осаду, що йде вниз, містить менше осаду, ніж середня посилка на рівні, куди вона надходить. Оскільки більше частинок осаду переноситься вгору, ніж вниз, чистий ефект - це транспорт вгору.
Іноді турбулентнуvt,s дифузійність осадової маси приймають рівноюvt турбулентної в'язкості води. Однак також можна стверджувати, що змішування води і осаду - це дві різні речі. Який би підхід не був прийнятий, зазвичай враховується загасання турбулентності через високі концентрації осаду. Мається на увазі вплив, який надають частинки осаду на турбулентну структуру рідини. Цей ефект стає все більш важливим для високих концентрацій осаду, які призводять до розшарування і, отже, демпфування турбулентності. Це впливає як на рух води, так і на розподіл осаду. Іноді використовуються емпіричні склади, які знижують вихрову в'язкість в залежності від концентрації осаду.
Нестаціонарне рівняння адвекції-дифузії тепер говорить:
∂c∂t−ws∂c∂z−∂∂zvt,sδcδz=0
В якості нижньої граничної умови часто прописується концентрація осаду на певному рівні біля шару, так звана еталонна концентрація. Ця концентрація призначається на певному еталонному рівні, наприкладza=2D50, також залежно від рецептури навантаження на ліжко. Довідкова концентрація часто вважається функцією напруги зсуву шару (подібно до складів транспортування навантаження на ліжко). Отже, на ліжку реакція квазістійка, тоді як вища по вертикалі концентрація осаду відстає від напруги зсуву на шарі. Альтернативною граничною умовою є функція підбору, яка прописує вертикальний градієнт концентрації замість концентрації.
Детальні моделі, які також вирішують хвильовий прикордонний шар (Sect. 5.4.3) вирішити залежну від часу концентрацію у вертикалі з метою моделювання транспортування листового потоку та транспортування підвісних вантажів. Залежну від часу концентрацію зваженого осаду можна визначити за допомогою нестаціонарного рівняння адвекції-дифузії, такого як Eq. ???. Цей підхід передбачає моделювання турбулентності, яка збільшується і зменшується протягом хвильового періоду, і є складним і трудомістким (див. Також Intermezzo 6.1).