Пісок та сипкі матеріали

1. Чим особливі сипкі середовища?

2. Кут відкосу

Насипте пісок на поверхню — і утвориться конічна купа. Кут схилу θᵣ, за якого купа стабілізується, — це кут відкосу. Він визначається суто тертям між частинками:

Насправді тут є гістерезис: купа може стояти під максимальним кутом стійкості θₘ > θᵣ, доки її не збурять, а потім лавиною осідає назад до θᵣ. Ця бістабільність зумовлює переривчасті лавини з розподілом за розміром.

3. Силові ланцюги

Напруження в сипкому середовищі розподіляється не рівномірно, як у рідині, — воно поширюється вздовж силових ланцюгів: квазілінійних шляхів зерен, що несуть навантаження значно вище за середнє. Експерименти з фотопружними дисками роблять ці ланцюги видимими як яскраві лінії в поляризованому світлі.

Силові ланцюги — це мікромасштабний механізм, що стоїть за ефектом Янссена: у високому циліндричному силосі тиск біля дна насичується з висотою (на відміну від тиску рідини, що зростає лінійно), бо зерна утворюють арки й перенаправляють навантаження на стінки.

4. Клітинний піщаний алгоритм

Клітинний автомат «падіння піску» представляє світ як 2D-сітку пікселів. Кожна клітина містить тип матеріалу (повітря, пісок, вода, вогонь, …). На кожному такті до кожної клітини застосовуються правила оновлення:

Правило недетерміноване на кроці 2/3 (випадковий вибір ліворуч чи праворуч), щоб запобігти артефактам ґратки. Кут відкосу 45° у КА — це артефакт квадратної сітки; реальний пісок із більшою кількістю діагональних варіантів (8 напрямків) може наближатися до 30–35°.

5. Лавини та самоорганізована критичність

Бак, Танг і Візенфельд (1987) показали, що повільно керована піщана купа самоорганізується до критичного стану, де розміри лавин підкоряються степеневому закону — без характерного масштабу довжини чи часу. Це самоорганізована критичність (SOC):

SOC пропонують як пояснення для шуму 1/f, розподілів розмірів землетрусів (Гутенберга-Ріхтера), розмірів лісових пожеж, подій вимирання видів та крахів фінансових ринків — усі вони демонструють степеневу статистику розмірів подій.

6. Масштабування Бегнольда та щільний потік

Ральф Бегнольд (1941) охарактеризував текучі сипкі середовища за допомогою безрозмірного числа Бегнольда, що порівнює напруження від зіткнень зерен із в'язким напруженням:

У глибокій воді (мале Ba) зерна піску поводяться як в'язка суспензія. У повітрі (велике Ba) домінують міжзеренні зіткнення, і зсувне напруження масштабується з квадратом швидкості зсуву — ключовий фактор у формуванні еолових дюн.

7. Симуляція падіння піску на JavaScript

// Симуляція падіння піску на клітинних автоматах
const W = 200, H = 150;
const EMPTY = 0, SAND = 1, WALL = 2;
const grid = new Uint8Array(W * H);
const idx  = (x, y) => y * W + x;
const get  = (x, y) => (x < 0 || x >= W || y < 0 || y >= H) ? WALL : grid[idx(x, y)];
const set  = (x, y, v) => { if (x >= 0 && x < W && y >= 0 && y < H) grid[idx(x, y)] = v; };

function stepSand() {
  // Ітеруємо знизу вгору, випадковий вибір ліво-право, щоб уникнути упередження
  const lr = Math.random() < 0.5;
  for (let y = H - 2; y >= 0; y--) {
    for (let x = 0; x < W; x++) {
      if (grid[idx(x, y)] !== SAND) continue;
      if (get(x, y + 1) === EMPTY) {
        set(x, y + 1, SAND); set(x, y, EMPTY);
      } else {
        const dl = get(x - 1, y + 1) === EMPTY;
        const dr = get(x + 1, y + 1) === EMPTY;
        if (dl && dr) {
          const dx = lr ? -1 : 1;
          set(x + dx, y + 1, SAND); set(x, y, EMPTY);
        } else if (dl) {
          set(x - 1, y + 1, SAND); set(x, y, EMPTY);
        } else if (dr) {
          set(x + 1, y + 1, SAND); set(x, y, EMPTY);
        }
      }
    }
  }
}

// Палітра кольорів: колір піщаного пікселя за висотою для ілюзії глибини
function sandColour(y) {
  const t = y / H;
  const r = Math.round(194 + t * 20);
  const g = Math.round(160 - t * 20);
  return `rgb(${r},${g},80)`;
}

// SOC-модель піщаної купи BTW для навчального порівняння
function btwTopple(z, W, H) {
  let avalanche = 0;
  let active = true;
  while (active) {
    active = false;
    for (let y = 1; y < H - 1; y++)
      for (let x = 1; x < W - 1; x++) {
        if (z[y * W + x] >= 4) {
          z[y * W + x] -= 4;
          z[(y-1) * W + x]++; z[(y+1) * W + x]++;
          z[y * W + (x-1)]++;  z[y * W + (x+1)]++;
          avalanche++; active = true;
        }
      }
  }
  return avalanche;
}

8. Застосування

• Цивільна та гірнича інженерія: проєктування силосів, аналіз стійкості схилів, будівництво насипів та підземне тунелювання — усе це вимагає розуміння потоку сипких матеріалів та динаміки силових ланцюгів.
• Фармацевтичне виробництво: пресування таблеток, змішування порошків та наповнення капсул залежать від надійного потоку порошку — сегрегація за розміром зерна є постійною проблемою під час змішування.
• Еолова геоморфологія: міграція пустельних дюн, сальтація піску та хвилясті патерни постають із взаємодії зерен і повітря в режимі Бегнольда — досліджуються для планетарної науки про поверхні (Марс, Титан).
• Прогнозування лавин: моделі снігових лавин поєднують фізику кута відкосу з даними про погоду та геометрією схилу, щоб видавати попередження.
• Ігрові рушії: ігри про падіння піску (Powder Toy, Noita) реалізують розширені клітинні автомати з вогнем, водою, кислотою та електричними взаємодіями — популярний жанр, що виник із цієї простої фізики.