Trójwymiarowe modelowanie do dokładnego przewidywania krytycznej średnicy oparty na trajektoriach w deterministycznej lateralnej dyspersji mikroprzepływów

Sortowanie maleńkich cząstek w maleńkich labiryntach

Wyobraź sobie test laboratoryjny potrafiący w kilka minut wyłowić rzadkie komórki nowotworowe albo wirusy z próbki krwi, używając tylko kropli płynu i przezroczystej plastikowej płytki. Ten artykuł bada jedną z kluczowych technologii stojących za takimi testami — mikroskopijny „tor przeszkód” dla cząstek zwany deterministyczną lateralną dyspersją (DLD) — i przedstawia nową metodę przewidywania, które cząstki zostaną rozdzielone, a które prześlizgną się dalej.

Jak mikro-tory przeszkód sortują według rozmiaru

Urządzenia DLD to płaskie kanały mikrofluidyczne wypełnione regularnie rozmieszczonymi filarami, niczym uporządkowany las słupków. Płyn przepływa równomiernie przez szczeliny. Małe cząstki podążają za liniami prądu płynu i przechodzą prosto w zygzakowatym wzorze. Większe cząstki nie mieszczą się jednak w najwęższych ścieżkach przepływu; wielokrotnie zderzają się z filarami i są odsuwane na boki, ostatecznie trafiając do innego wylotu. Graniczny rozmiar, który decyduje, czy cząstka zygzakuje, czy zostaje odrzucona, nazywany jest krytyczną średnicą. Znajomość tej wartości z wyprzedzeniem jest niezbędna do projektowania chipów, które niezawodnie rozdzielają komórki, krople czy nanocząstki do zastosowań diagnostycznych i badawczych.

Dlaczego dotychczasowe zasady projektowania zawodzą

Do tej pory większość zasad projektowania urządzeń DLD traktowała cząstki jako idealne punkty, a kanał jako idealnie dwuwymiarowy. Proste wzory lub modele komputerowe szacowały krytyczną średnicę używając tylko rozstawu filarów w płaszczyźnie. Tymczasem rzeczywiste urządzenia mają skończoną wysokość, a płyn zwalnia w pobliżu ścianek górnej i dolnej. Filary mogą nie być idealnie okrągłe, być rozmieszczone nierównomiernie lub powstać z drobnymi wadami produkcyjnymi. Wcześniejsze symulacje trójwymiarowe opierały się albo na empirycznych współczynnikach dopasowania, które zmieniały się z projektem, albo były tak obciążające obliczeniowo, że nie nadawały się do rutynowego stosowania. W efekcie przewidywania rozdziału rozmiarów często były niedokładne, zwłaszcza dla bardziej zaawansowanych kształtów filarów lub ściśle dostrojonych urządzeń.

Trójwymiarowa mapa sił działających na każdą cząstkę

Autorzy prezentują nowe trójwymiarowe ramy modelowania, które podchodzą do problemu z punktu widzenia cząstki. Najpierw obliczają szczegółowe pole przepływu 3D w małym, reprezentatywnym bloku czterech sąsiednich filarów przy użyciu oprogramowania elementów skończonych. Następnie, zamiast traktować cząstkę jako punkt, dzielą powierzchnię sferycznej cząstki na wiele drobnych fragmentów. Dla każdego fragmentu obliczają, jak lokalne prędkości i ciśnienia płynu działają na cząstkę — uwzględniając lepki opór, siły ciśnieniowe oraz subtelne siły unoszenia tworzone przez gradienty prędkości i pobliskie ścianki. Te lokalne siły są łączone, by krok po kroku aktualizować ruch cząstki. Śledząc wiele rozmiarów cząstek przez powtarzające się kopie tego samego bloku czterech filarów, metoda ujawnia, czy dany rozmiar podąża ścieżką zygzakową, ścieżką odbić, czy czymś pośrednim.

Ukryte trzecie zachowanie w wymiarze pionowym

Stosując podejście 3D, badacze odkryli, że krytyczna średnica nie jest jedną stałą liczbą, lecz zmienia się w zależności od wysokości w kanale. Tworzy ona w rzeczywistości krzywą w kształcie litery U: cząstki na środku wysokości są separowane przy najmniejszym rozmiarze, natomiast te blisko górnej i dolnej ścianki wymagają większej różnicy rozmiaru, by zostać odchylone. Między tymi skrajnościami znajduje się strefa przejściowa, w której cząstka o danym rozmiarze może przełączać się między zygzakiem a odbiciem, gdy subtelnie oscyluje w górę i w dół. To mieszane zachowanie tworzy „zmieniony zygzak” o netto przesunięciu bocznym słabszym i bardziej zmiennym niż w przypadku czystego odbijania. Symulacje zespołu zgadzają się z opublikowanymi eksperymentami i nowymi testami na specjalnie wykonanych chipach — zmierzone trajektorie cząstek pokrywają się w granicach około mikrometra.

Projektowanie ostrzejszych i mądrzejszych układów sortujących

Dla nie-specjalistów kluczowa konkluzja jest taka, że pionowa struktura przepływu — a nie tylko układ filarów widziany z góry — w dużym stopniu wpływa na to, jak dobrze urządzenie DLD potrafi rozróżniać cząstki o podobnych rozmiarach. Poprzez jawne modelowanie sił 3D nowe ramy potrafią przewidzieć i wyjaśnić niejednoznyste trajektorie, które wcześniej rozmywały wydajność separacji. Pokazują też, jak pewne kształty filarów, na przykład odwrócone trójkątne, mogą zmniejszyć strefę przejściową i wyostrzyć rozdzielczość urządzenia. Ponieważ metoda wykorzystuje tylko jedno ustalone rozwiązanie przepływu 3D, a następnie efektywnie je ponownie używa, oferuje praktyczne narzędzie do szybkiego badania nowych geometrii chipów. W dłuższej perspektywie autorzy przewidują połączenie tego modelu opartego na fizyce z automatyzacją, tak aby separatory mikrofluidyczne mogły być projektowane na żądanie do zadań od izolacji rzadkich komórek po diagnostykę przy łóżku pacjenta.

Cytowanie: Chen, J., Huang, X., Xuan, W. et al. A 3D modeling framework for accurate trajectory-based prediction of critical diameter in deterministic lateral displacement microfluidics. Microsyst Nanoeng 12, 78 (2026). https://doi.org/10.1038/s41378-025-01139-3

Słowa kluczowe: separacja mikroprzepływowa, deterministyczna lateralna dyspersja, sortowanie cząstek, lab-on-a-chip, analiza komórek i nanocząstek