Нелинейный фриктиофорез

редактировать

Нелинейный фриктиофорез - это однонаправленный дрейф частицы в среде, вызванный периодической движущей силой с нулевым средним. Эффект возможен из-за нелинейной зависимости силы трения-сопротивления от скорости частицы. Он был открыт теоретически и известен в основном как нелинейный электрофриктиофорез. На первый взгляд, периодическая движущая сила с нулевым средним значением способна вовлечь частицу в колебательное движение без однонаправленного дрейфа, потому что интегральный импульс, передаваемый частице силой, равен нулю. Возможность однонаправленного дрейфа можно распознать, если учесть, что частица сама теряет импульс, передавая его дальше среде, в которой она движется. Если трение нелинейно, то может случиться так, что потеря количества движения при движении в одном направлении не будет равна потере импульса в противоположном направлении, и это вызовет однонаправленный дрейф. Чтобы это произошло, зависимость движущей силы от времени должна быть более сложной, чем в случае одной синусоидальной гармоники.

Содержание
  • 1 Простой пример - пластик Бингема
    • 1.1 Нелинейное трение
    • 1.2 Периодическая движущая сила
    • 1.3 Однонаправленный дрейф
  • 2 Математический анализ
  • 3 Приложения
    • 3.1 Разделение ДНК по длине
    • 3.2 Разделение по удельной массе
  • 4 Удлинения
    • 4.1 Вращение с помощью круговой волны
    • 4.2 Модификация потенциальной функции
    • 4.3 Другая нелинейность
  • 5 Ссылки
Простой пример - пластик Бингема

Нелинейное трение

Простейшим случаем зависимости трения от скорости является закон Стокса :

(1) F dr ( v) знак равно λ v, {\ displaystyle (1) \ qquad F_ {dr} (v) = \ lambda v,}{\ displaystyle (1) \ qquad F_ {dr} (v) = \ lambda v,}

где F dr (v) {\ displaystyle F_ {dr} (v)}{\ displaystyle F_ {dr} (v)} - сила трения / сопротивления, приложенная к частице, движущейся со скоростью v {\ displaystyle v}v в среде. Закон трения-скорости (1) наблюдается для медленно движущейся сферической частицы в ньютоновской жидкости.

Рис. 1 Линейное трение

Оно линейное, см. Рис. 1, и не подходит для проведения нелинейного фриктиофореза. Характерным свойством закона (1) является то, что любая, даже очень небольшая движущая сила способна заставить частицу двигаться. Это не относится к таким материалам, как пластик Bingham. Для этих сред необходимо применить некоторую пороговую силу, d {\ displaystyle d}d , чтобы заставить частицу двигаться. Этот вид закона трения-скорости (сухое трение) имеет скачкообразный разрыв в v = 0 {\ displaystyle v = 0}v = 0 :

(2) F dr (v) = λ v + d ⋅ sign (v). {\ displaystyle (2) \ qquad F_ {dr} (v) = \ lambda v + d \ cdot \ mathrm {sign} (v).}{ \ Displaystyle (2) \ qquad F_ {dr} (v) = \ лямбда v + d \ cdot \ mathrm {sig п} (v).}
Рис. 2 Пример нелинейного трения

Он нелинейный, см. Рис. 2, и используется в этом примере.

Периодическая движущая сила

Пусть T>0 {\ displaystyle T>0}T>0 обозначают период движущей силы. Выберите значение времени t 1 {\ displaystyle t_ {1} }t_ {1} такой, что 0 < t 1 < T {\displaystyle 0{\ displaystyle 0 <t_ {1} <T} и два значения силы, F +>0 {\ displaystyle F ^ {+}>0}{\displaystyle F^{+}>0} , F - < 0 {\displaystyle F^{-}<0}{\ displaystyle F ^ {-} <0} таким образом, чтобы выполнялись следующие отношения :

F +>d, | F - | < d, {\displaystyle \qquad \qquad F^{+}>d, \ quad | F ^ {-} | {\displaystyle \qquad \qquad F^{+}>d, \ quad | F ^ {-} | <d,}
(3) F + t 1 + F - (T - t 1) = 0. {\ displaystyle (3) \ qquad F ^ {+} t_ {1} + F ^ {-} (T-t_ {1}) = 0.}{\ displaystyle (3) \ qquad F ^ {+} t_ {1} + F ^ {-} (T-t_ {1}) = 0.}

Периодическая движущая сила f ( t) {\ displaystyle f (t)}е (т) , используемый в этом примере, выглядит следующим образом:

(4) f (t) = {F +, if 0 < t ≤ t 1, F −, if t 1 < t ≤ T, f ( t + T) = f ( t). {\displaystyle (4)\qquad f(t)={\begin{cases}F^{+},{\text{ if }}0{\ displaystyle (4) \ qquad f (t) = {\ begin {cases} F ^ {+}, {\ text {if} } 0 <t \ leq t_ {1}, \\ F ^ {-}, {\ text {if}} t_ {1} <t \ leq T, \ quad f (t + T) = f (t). \ end { case}}}

Очевидно, что из-за to (3), f (t) {\ displaystyle f (t)}е (т) имеет нулевое среднее значение:

∫ 0 T f (t) dt = F + t 1 + F - ( T - t 1) = 0. {\ displaystyle \ int \ limits _ {0} ^ {T} f (t) dt = F ^ {+} t_ {1} + F ^ {-} (T-t_ {1 }) = 0.}{\ displaystyle \ int \ limits _ { 0} ^ {T} f (t) dt = F ^ {+} t_ {1} + F ^ {-} (T-t_ {1}) = 0.}
Рис. 3 Пример с нулевой средней движущей силой

См. Также рис. 3.

Однонаправленный дрейф

Для простоты мы рассматриваем здесь физическую ситуацию, когда инерцией можно пренебречь. Последнее может быть достигнуто, если масса частицы мала, скорость мала и трение велико. Эти условия должны гарантировать, что τ ≪ t 1 {\ displaystyle \ tau \ ll t_ {1}}{\ displaystyle \ tau \ ll t_ {1}} , где τ {\ displaystyle \ tau}\ tau равно время релаксации. В этой ситуации частица, движущаяся силой (4), немедленно начинает двигаться с постоянной скоростью v + = 1 λ (F + - d) {\ displaystyle v ^ {+} = {\ frac {1} {\ lambda }} (F ^ {+} - d)}{\ displaystyle v ^ {+} = {\ frac {1} {\ lambda}} (F ^ {+} - d)} в течение интервала 0 < t ≤ t 1 {\displaystyle 0{\ displaystyle 0 <t \ leq t_ {1}} и немедленно прекратит движение в течение интервала t 1 < t ≤ T {\displaystyle t_{1}{\ displaystyle t_ {1} <t \ leq T} , см. Рис. 4.

Рис. 4 Скорость с ненулевым средним

Это приводит к положительной средней скорости однонаправленного дрейфа:

v (t) ¯ = 1 T ∫ 0 T v (t) dt = t 1 λ T (F + - d)>0. {\ displaystyle \ qquad \ qquad {\ overline {v (t)}} = {\ frac {1} {T}} \ int \ limits _ {0} ^ {T} v (t) dt = {\ frac { t_ {1}} {\ lambda T}} (F ^ {+} - d)>0.}{\displaystyle \qquad \qquad {\overline {v(t)}}={\frac {1}{T}}\int \limits _{0}^{T}v(t)dt={\frac {t_{1}}{\lambda T}}(F^{+}-d)>0.}
Математический анализ

Анализ возможности получения ненулевого дрейфа с помощью периодической силы было сделано в. Безразмерное уравнение движения частицы, движущейся под действием периодической силы f (t) {\ displaystyle f (t)}е (т) , f (t + 1) = f (t) {\ displaystyle f ( t + 1) знак равно е (t)}{\ displaystyle f (t + 1) = f (t)} , ∫ 0 1 f (t) dt = 0 {\ displaystyle \ int _ {0} ^ {1} f (t) dt = 0}{\ displaystyle \ int _ {0} ^ {1} е (t) dt = 0} выглядит следующим образом:

(5) v ˙ + λ v - ϵ g (v) = f (t), {\ displaystyle (5) \ qquad {\ dot {v}} + \ lambda v- \ epsilon g (v) = f (t),}{\ displaystyle (5) \ qquad {\ dot {v}} + \ lambda v- \ epsilon g (v) = f (t),}

где сила трения / сопротивления F dr (v) = λ v - ϵ g (v) {\ displaystyle F_ {dr} (v) = \ lambda v - \ epsilon g (v)}{\ displaystyle F_ {dr} (v) = \ lambda v- \ epsilon g (v)} удовлетворяет следующему:

F dr (- v) = - F dr (v), ddv F dr ( v) ≥ 0. {\ displaystyle \ qquad \ qquad F_ {dr} (- v) = - F_ {dr} (v), \ quad {\ frac {d} {dv}} F_ {dr} (v) \ geq 0.}{\ displaystyle \ qquad \ qquad F_ {dr} (- v) = - F_ {dr} (v), \ quad {\ frac {d} {dv}} F_ {dr} (v) \ geq 0.}

Доказано, что любое решение (5) переходит в периодический режим v ∗ (t) {\ displaystyle v ^ {*} (t)}{\ displaystyle v ^ {*} (t)} , v ∗ (t + 1) = v ∗ (t) {\ displaystyle v ^ {*} (t + 1) = v ^ {*} (t)}{\ displaystyle v ^ {*} (t + 1) = v ^ {*} (t)} , что имеет ненулевое среднее значение:

v ∗ ( t) ¯ знак равно ∫ 0 1 v ∗ (t) dt ≠ 0, {\ displaystyle \ qquad \ qquad {\ overline {v ^ {*} (t)}} = \ int \ limits _ {0} ^ {1} v ^ {*} (t) dt \ neq 0,}{\ displaystyle \ qquad \ qquad {\ overline {v ^ {*} (t)}} = \ int \ limits _ {0} ^ {1} v ^ {*} (t) dt \ neq 0,}

почти наверняка при условии, что f (t) {\ displaystyle f (t)}е (т) не является антипериодическим.

Для g (v) = v 3 {\ displaystyle g (v) = v ^ {3}}{\ displaystyle g (v) = v ^ {3}} , два случая f (t) {\ displaystyle f (t)}е (т) были учтены явно:

1. Пилообразная движущая сила, см. Рис. 5:

Рис. 5 Пилообразная движущая сила
f (t) = a t, t ∈ [- 1/2; 1/2], f (t + 1) = f (t), t ∈] - ∞; ∞ [. {\ Displaystyle \ qquad \ qquad е (т) = ат, \ четырехъядерный т \ в [-1/2; 1/2], \ четырехъядерный е (т + 1) = е (т), \ четырехъядерный т \ в] - \ infty; \ infty [.}{\ displaystyle \ qquad \ qquad f (t) = at, \ quad t \ in [-1/2; 1/2], \ quad f (t + 1) = f (t), \ quad t \ in] - \ infty; \ infty [.}

В данном случае найдено в первом порядке в приближении ϵ {\ displaystyle \ epsilon}\ epsilon к v ∗ (t) {\ displaystyle v ^ {*} (t)}{\ displaystyle v ^ {*} (t)} , v 1 ∗ (t) {\ displaystyle v_ {1} ^ {*} (t)}{\ displaystyle v_ {1} ^ {*} (t)} , имеет следующее среднее значение:

| v 1 ∗ (t) ¯ | ≡ | ∫ 0 1 v 1 ∗ (t) d t | ≥ 2 3 ϵ а 3 / λ 5. {\ Displaystyle \ qquad \ qquad {\ Big |} {\ overline {v_ {1} ^ {*} (t)}} {\ Big |} \ Equiv {\ Big |} \ int \ limits _ {0} ^ {1} v_ {1} ^ {*} (t) dt {\ Big |} \ geq {\ frac {2} {3}} \ epsilon a ^ {3} / \ lambda ^ {5}.}{\ displaystyle \ qquad \ qquad {\ Big |} {\ overline {v_ { 1} ^ {*} (t)}} {\ Big |} \ Equiv {\ Big |} \ int \ limits _ {0} ^ {1} v_ {1} ^ {*} (t) dt {\ Big |} \ geq {\ frac {2} {3}} \ epsilon a ^ {3} / \ lambda ^ {5}.}

Эта оценка сделана с ожиданием λ ≫ 1 {\ displaystyle \ lambda \ gg 1}{\ displaystyle \ lambda \ gg 1} .

2. Движущая сила двух гармоник,

f (t) = a cos ⁡ (2 π t) + b cos ⁡ (4 π t + ψ). {\ displaystyle \ qquad \ qquad f (t) = a \ cos (2 \ pi t) + b \ cos (4 \ pi t + \ psi).}{\ displaystyle \ qquad \ qquad f (T) знак равно a \ соз (2 \ пи T) + б \ соз (4 \ пи t + \ фунт / кв. дюйм).}

В этом случае первый порядок в ϵ {\ displaystyle \ epsilon}\ epsilon аппроксимация имеет следующее среднее значение:

| v 1 ∗ (t) ¯ | = | ∫ 0 1 v 1 ∗ (t) d t | ≥ 2 27 ϵ λ (M λ) 3 1 (1 + 4 π 2 λ 2) (1 + 16 π 2 λ 2) 1/2. {\ displaystyle \ qquad \ qquad {\ Big |} {\ overline {v_ {1} ^ {*} (t)}} {\ Big |} = {\ Big |} \ int \ limits _ {0} ^ { 1} v_ {1} ^ {*} (t) dt {\ Big |} \ geq {\ frac {2} {27}} {\ frac {\ epsilon} {\ lambda}} \ left ({\ frac { M} {\ lambda}} \ right) ^ {3} {\ frac {1} {(1 + {\ frac {4 \ pi ^ {2}} {\ lambda ^ {2}}}) (1+ { \ frac {16 \ pi ^ {2}} {\ lambda ^ {2}}}) ^ {1/2}}}.}{\ displaystyle \ qquad \ qquad {\ Big |} {\ overline {v_ {1} ^ {*} (t)}} {\ Big |} = {\ Big |} \ int \ limits _ {0} ^ {1} v_ {1} ^ {*} (t) dt {\ Big |} \ geq {\ frac {2} {27}} {\ frac {\ epsilon} {\ lambda}} \ left ({\ frac {M} {\ lambda}} \ right) ^ {3} {\ frac { 1} {(1 + {\ frac {4 \ pi ^ {2}} {\ lambda ^ {2}}}) (1 + {\ frac {16 \ pi ^ {2}} {\ lambda ^ {2}) }}) ^ {1/2}}}.}

Это значение максимизируется в ψ {\ displaystyle \ psi}\ psi , a, b {\ displaystyle a, b}a, b , сохраняя константу a + b = M {\ displaystyle a + b = M}{\ displaystyle a + b = M} . Интересно, что значение дрейфа зависит от ψ {\ displaystyle \ psi}\ psi и дважды меняет свое направление по мере того, как ψ {\ displaystyle \ psi}\ psi проходит в интервале [0; 2 π] {\ displaystyle [0; 2 \ pi]}{\ displaystyle [0; 2 \ pi]} . Другой тип анализа, основанный на нарушении симметрии, также предполагает, что нулевое среднее значение движущей силы способно вызвать направленный дрейф.

Применения
Рис.6 (а): сплошная линия - сила сопротивления на заряд на одиночном б.п. от скорости, пунктирная линия - линейная аппроксимация для сравнения. (b): то же, что и (a), но в мелком масштабе

В приложениях характер силы f (t) {\ displaystyle f (t)}е (т) в (5), обычно электрическая, подобно силам, действующим во время стандартного электрофореза. Единственное отличие состоит в том, что сила периодическая и без постоянной составляющей.

Чтобы эффект проявился, зависимость силы трения / сопротивления от скорости должна быть нелинейной. Так обстоит дело с многочисленными веществами, известными как неньютоновские жидкости. Среди них гели и дилатантные жидкости, псевдопластические жидкости, жидкие кристаллы. Специализированные эксперименты определили F d r (v) {\ displaystyle F_ {dr} (v)}{\ displaystyle F_ {dr} (v)} для стандартной лестницы ДНК длиной до 1500 п.н. в 1,5% агарозном геле. Найденная зависимость, см. Рис. 6, подтверждает возможность нелинейного фриктиофореза в такой системе. На основе данных на рис. 6 был найден оптимальный временной ход для возбуждения электрического поля с нулевым средним, E (t) {\ displaystyle E (t)}E (t) , в котором обеспечивается максимальное дрейф на 1500 б.п. длинный фрагмент, см. рис. 7.

Рис. 7 Временной ход оптимального электрического поля, E (t) {\ displaystyle E (t)}E (t) , E (t ± 1 s) = E (t) {\ displaystyle E (t \ pm 1s) = E (t)}{\ displaystyle E (t \ pm 1s) = E (t)} , установившаяся скорость, v ∗ (t) {\ displaystyle v ^ {*} (t)}{\ displaystyle v ^ {*} (t)} , и смещение, X (т) {\ Displaystyle X (т)}X (t) . Значение сноса v ∗ (t) ¯ = 0,31 μ {\ displaystyle {\ overline {v ^ {*} (t)}} = 0,31 \ mu}{\ displaystyle {\ overline {v ^ {*} (t)}} = 0,31 \ mu} м / с.

Эффект однонаправленного дрейфа, вызванный периодической силой с нулевым интегральным значением, имеет своеобразную зависимость от времени действия приложенной силы. Примеры см. В предыдущем разделе. Это открывает новое измерение для ряда проблем разделения.

Разделение ДНК по длине

При разделении фрагментов ДНК в электрофорезе с нулевым интегральным полем (ZIFE) используется периодическое электрическое поле с нулевым средним периодом, где зависимость поля от времени аналогична показанный на рис. 3. Это позволяет разделять длинные фрагменты в агарозном геле, не разделяемые стандартным электрофорезом в постоянном поле. Длинная геометрия ДНК и способ ее движения в геле, известный как рептация, не позволяют напрямую применять соображения, основанные на уравнении. (5), выше.

Разделение по удельной массе

Было замечено, что при определенных физических условиях механизм, описанный в разделе «Математический анализ» выше, можно использовать для разделения по удельной массе, например, частицы изготовлены из изотопов того же материала.

Расширения

Идея организации направленного дрейфа с периодическим возбуждением с нулевым средним получила дальнейшее развитие для других конфигураций и других физических механизмов нелинейности.

Вращение с помощью круговой волны

электрический диполь, свободно вращающийся вокруг z {\ displaystyle z}z -оси в среде с нелинейным трением можно управлять, применяя электромагнитную волну, поляризованную по кругу вдоль ± z {\ displaystyle \ pm z}{\ displaystyle \ pm z} и состоящую из двух гармоник. Уравнение движения для этой системы выглядит следующим образом:

ω ˙ + λ ω - ϵ g (ω) = f (t, θ (t)), θ ˙ = ω, {\ displaystyle \ qquad \ qquad {\ точка {\ omega}} + \ lambda \ omega - \ epsilon g (\ omega) = f (t, \ theta (t)), \ quad {\ dot {\ theta}} = \ omega,}{\ displaystyle \ qquad \ qquad {\ dot {\ omega }} + \ lambda \ omega - \ epsilon g (\ omega) = f (t, \ theta (t)), \ quad {\ dot {\ theta}} = \ omega,}

где f (t, θ (t)) {\ displaystyle f (t, \ theta (t))}{\ displaystyle f (t, \ theta (t))} - крутящий момент, действующий на диполь из-за круговой волны:

(6) f (t, θ (t)) ∼ | p | (A 1 соз ⁡ (ω T ± θ) + A 2 соз ⁡ (2 ω T ± θ + ψ)), {\ displaystyle (6) \ qquad f (t, \ theta (t)) \ sim | p | (A_ {1} \ cos (\ omega t \ pm \ theta) + A_ {2} \ cos (2 \ omega t \ pm \ theta + \ psi)),}{\ displaystyle (6) \ qquad f (t, \ theta (t)) \ sim | p | (A_ {1} \ cos (\ omega t \ pm \ theta) + A_ {2} \ cos (2 \ omega t \ pm \ theta + \ psi)),}

где p {\ displaystyle p}p - составляющая дипольного момента, ортогональная оси z {\ displaystyle z}z , а θ {\ displaystyle \ theta}\ theta определяет направление диполя в плоскости XY {\ displaystyle XY}XY . Выбрав правильный фазовый сдвиг ψ {\ displaystyle \ psi}\ psi в (6), можно ориентировать диполь в любом желаемом направлении, θ 0 {\ displaystyle \ theta _ {0 }}\ theta _ {0} . Направление θ 0 {\ displaystyle \ theta _ {0}}\ theta _ {0} достигается за счет направленного углового дрейфа, который становится равным нулю, когда θ = θ 0 {\ displaystyle \ theta = \ theta _ {0}}\ theta = \ theta _ {0} . Небольшая расстройка между первой и второй гармониками в (6) приводит к непрерывному вращательному дрейфу.

Модификация потенциальной функции

Рис. 8 Пример модификации потенциальной функции U (x) {\ displaystyle U (x)}U (x) из-за нелинейного фриктиофореза. (a) начальный U (x) {\ displaystyle U (x)}U (x) , (b) модифицированный U (x) {\ displaystyle U (x)}U (x) .

Если a частица совершает направленный дрейф при свободном перемещении в соответствии с формулой. (5), то он дрейфует аналогично, если наложено достаточно мелкое потенциальное поле U (x) {\ displaystyle U (x)}U (x) . Уравнение движения в этом случае:

v ˙ + λ v - ϵ g (v) = f (t) - ϕ (x), x ˙ = v, {\ displaystyle \ qquad \ qquad {\ dot {v) }} + \ lambda v- \ epsilon g (v) = f (t) - \ phi (x), \ quad {\ dot {x}} = v,}{\ displaystyle \ qquad \ qquad {\ dot { v}} + \ lambda v- \ epsilon g (v) = f (t) - \ phi (x), \ quad {\ dot {x}} = v,}

где ϕ (x) { \ displaystyle \ phi (x)}\ phi (x) - сила, создаваемая потенциальным полем. Дрейф продолжается до тех пор, пока не будет достигнута достаточно крутая область в течение U (x) {\ displaystyle U (x)}U (x) , которая может остановить дрейф. Такое поведение, как показывает строгий математический анализ, приводит к модификации U (x) {\ displaystyle U (x)}U (x) путем добавления линейного в x {\ displaystyle x}x срок. Это может изменить U (x) {\ displaystyle U (x)}U (x) качественно, например, изменение количества точек равновесия, см. рис. 8. Эффект может быть существенным при воздействии высокочастотного электрического поля на биополимеры.

Другая нелинейность

Для электрофореза коллоидных частиц в электрическом поле малой напряженности сила f (t) {\ displaystyle f (t)}е (т) в правой части уравнения. (5) линейно пропорциональна прочности E (t) {\ displaystyle E (t)}E (t) приложенного электрического поля. При высокой прочности линейность нарушается из-за нелинейной поляризации. В результате сила может нелинейно зависеть от приложенного поля:

f (t) ∼ E (t) + α (E (t)) 3. {\ displaystyle \ qquad \ qquad f (t) \ sim E (t) + \ alpha (E (t)) ^ {3}.}{\ displaystyle \ qquad \ qquad f (t) \ sim E (t) + \ альфа (E (t)) ^ {3}.}

В последнем выражении, даже если примененное поле, E (t) {\ displaystyle E (t)}E (t) имеет нулевое среднее значение, приложенная сила f (t) {\ displaystyle f (t)}е (т) может иметь постоянная составляющая, которая может вызвать направленный дрейф. Как и выше, для этого E (t) {\ displaystyle E (t)}E (t) должен иметь более одной синусоидальной гармоники. Такой же эффект для жидкости в трубке может использоваться в электроосмотическом насосе, приводимом в действие с нулевым средним электрическим полем.

Ссылки
Последняя правка сделана 2021-05-31 12:12:13
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте