두 표면 사이의 비스코스 유체 흐름 모델 상호 상대적으로 이동
유체 역학 에서 쿠엣 흐름 은 두 표면 사이의 공간에 있는 점성 유체 의 흐름이며, 그 중 하나는 다른 표면과 비교하여 접선적 으로 움직인다. 표면의 상대적인 움직임은 유체에 전단 응력을 가하고 유동을 유도한다. 용어의 정의에 따라 흐름 방향에 적용된 압력 구배 가 있을 수도 있다.
쿠엣 구성은 지구의 맨틀 과 대기 와 같은 특정한 실제적인 문제를 모델링하며 가벼운 하중으로 저널 베어링 을 흐른다.[1] 또한 점도계 에 사용되며 가역성 의 근사치를 나타낸다.[2] [3]
19세기 말 프랑스 앵거스 대학 의 물리학과 교수인 모리스 쿠엣 의 이름을 따서 지어졌다.
플라나르 쿠에트 흐름 두 개의 무한 평면 플레이트를 사용한 간단한 쿠엣 구성. 쿠엣 흐름은 전단 구동 유체 운동을 설명하기 위해 학부 물리학과 공학 과정에서 자주 사용된다. 단순한 구성은 거리 h {\displaystyle h} 로 분리된 두 개의 무한 평행 판에 해당한다. 하나의 판은 자신의 평면에서 일정한 상대 속도 U {\displaystyle U} 로 번역된다. 압력 그라데이션 무시, Navier –스토크 방정식 은 다음과 같이 단순화된다.
d 2 u d y 2 = 0 , {\displaystyle {\frac {d^{2}u}{dy^{2}}=0,} 여기서 y {\displaystyle y} 은 플레이트에 정규적인 공간 좌표고 u ( ) {\displaystyle u(y)} 은 속도 필드다. 이 방정식은 흐름이 단방향 이라는 가정을 반영한다. 즉, 세 가지 속도 성분(u , v , w ) 중 하나만 비삼각형(u,v,w) 이다 . 하단 플레이트가 y = 0 {\displaystyle y=0 } 에 해당하는 경우 경계조건은 u ( 0 ) = 0 {\displaystyle u(0)=0 } 및 u ( h ) = U {\displaysty u(h)= 이다. U} . 정확한 해결책
u ( y ) = U y h u(y)= U{\frac{y}{h}} 두 번 통합하고 경계 조건을 사용한 상수에 대한 해결로 찾을 수 있다. 흐름의 주목할 만한 측면은 전단 응력 이 영역 전체에 걸쳐 일정하다는 것이다. 특히 속도의 첫 번째 파생상품인 U/ h {\displaystyle U/h } 은 일정하다. 뉴턴의 점도의 법칙 (뉴턴 액)에 따르면 전단 응력은 이 표현과 (정수) 유체 점도 의 산물이다.
시작 현실적으로 쿠엣 해법은 순간적으로 도달하지 못한다. 안정 상태에 대한 접근법을 설명하는 "시동 문제"는 다음에서 제시한다.
∂ u ∂ t = ν ∂ 2 u ∂ y 2 {\displaystyle {\frac {\fract t}}{\preason t}=\nu {\fract {\preason ^{2}u}{\preased y^{2}}: 초기의 조건에 따라.
u ( y , 0 ) = 0 , 0 < y < h , {\displaystyle u(y,0)=0,\pair 0(h)} 그리고 일정한 흐름과 동일한 경계 조건의 경우:
u ( 0 , t ) = 0 , u ( h , t ) = U , t > 0. {\displaystyle u(0,t)=0,\quad u(h,t)=U,\quad t>0. } 그 문제는 일정한 용액을 빼서 균일 하게 만들 수 있다. 그런 다음 변수 분리 를 적용하면 다음과 같은 해결책이 나온다.[4]
u ( y , t ) = U y h - 2 U π n = 1 ∞ 1 e - n n 2 ν t h 2 sin [ n π ( 1 - y h ) ] {\ displaysty u(y,t)= U{\frac {y}{h}-{\frac {2 } U}{\pi }{\pi }\ sum _{n=1 }^{n1}{n1}e^{ n1}{n^{2 }\pi ^{2}{\frac {\nu t}{h^{2}}:}\sin \left[n-{y}}}{h}}}}}}오른쪽 ]}}}. 안정 상태에 대한 이완을 설명하는 시간 척도는 그림에서와 같이 t ~ h 2 / ν {\displaystyle t\sim h^{2}/\nu } 이다. 안정 상태에 도달하는 데 필요한 시간은 플레이트 h {\displaystyle h} 과 유체의 동적 점성 사이의 간격에 따라 달라지지만 U {\displaystyle U} 에는 그렇지 않다.
압력 경사도가 있는 평면 흐름 보다 일반적인 쿠엣 흐름은 평판에 평행한 방향으로 일정 한 압력 그라데이션 G = - d p / d x = c o n s t {\ displaystyle G=-dp/dx=\mathrm {constant}} 을(를) 포함한다. 더 나비에 –스토크 방정식은
d 2 u d y 2 = − G μ , {\daystyle {\frac {d^{2}u}{dy^{2}}=-{\frac {G}{\mu }}},} 여기서 μ[\displaystyle \mu } 은 (는) 동적 점성 이다. 위의 방정식을 두 번 통합하고 경계 조건(압력 구배가 없는 쿠엣 흐름의 경우와 동일)을 적용하면 다음과 같은 결과를 얻을 수 있다.
u ( y ) = G 2 μ y ( h − y ) + U y h . {\displaystyle u(y)={\frac {G}{2\mu }}}y\,(h-y)++ U{\frac {y}{h}}. } 압력 구배는 양수(부압 구배) 또는 음수(우호 압력 구배)일 수 있다. 정지판의 제한 사례(U = 0 {\displaystyle U=0} )에서 유량을 평면 푸아세유 흐름 이라고 하며, 대칭(수평 중간 평면) 포물선 속도 프로파일을 참조한다.[5]
압축유동 M = 0 {\displaystyle \mathrm {M} =0} 에 대한 압축 가능한 쿠엣 흐름 M 2 P r = 7.5 {\displaystyle \mathrm {M} ^{2}\mathrm {Pr} = 7.5} 에 대한 압축식 쿠엣 흐름 압축 불가능한 흐름에서는 유체 온도가 일정하기 때문에 속도 프로파일이 선형이다. 상·하벽이 서로 다른 온도에서 유지될 때 속도 프로파일이 더 복잡해진다. 그러나 1950년 C. R. 일링워스에서 보듯이 정확한 암묵적 해결책을 가지고 있다.[6]
정지 상태에서 하부 벽이 있는 Couette 흐름과 등속 U {\displaystyle U} 을(를) 사용하여 상단 벽이 움직이는 평면 Couette 흐름을 고려하십시오. 하단 벽 의 유체 특성 은 첨자 w {\displaystyle w}, 상단 벽의 특성은 첨자 ∞ {\displaystyle \inft }. 속성과 압력. 벽은 기준 수량으로 규정되고 취해진다. l {\displaystyle l} 을(를) 두 벽 사이의 거리로 두십시오. 경계조건은
u = 0 , v = 0 , h = h w = c p w T w 에서 y = 0 , {\displaystyle u=0,\ v=0,\ h=h_{w}=c_{pw} T_{w}\{\text{at}\y=0,} u = U , v = 0 , h = h ∞ = c p ∞ T ∞ , p = p ∞ 에서 y = l {\displaystyle u=U,\ v=0,\h_{\\infit }=c_{p\infit }{\infit }\p=p_{\infit }\p=p_{\infit}\{\text{at}\}\y=l} 여기서 h {\displaystyle h} 은 (는) 특정 엔탈피 이고 c p {\ displaystyle c_{p} 는 특정 열 이다. 질량 및 y {\displaystyle y} -momentum을 보존하려면 흐름 영역 어디에나 v = 0, p = p p {\displaystyle v=0,\p=p_{\ft}}} 가 필요하다. 에너지 절약 및 x {\displaystyle x} -모멘텀 감소
d d y ( μ d u d y ) = 0 , ⇒ d τ d y = 0 , ⇒ τ = τ w {\displaystyle {\frac {d}{did}}\왼쪽(\mu {\frac {du}{dy}\오른쪽)=0,\quad \rightarrow \quad \frac {d\d\tau}}}=0,\quad \rigarrow \quau =\{w} 1 P r d d y ( μ d h d y ) + μ ( d u d y ) 2 = 0. {\displaystyle {\frac {1}{\mathrm {Pr}}}}}}{dy}}\좌(\mu {\frac {dh}{dy}}\우)+\mu \좌({\frac {du}{dy}}\우)^{2}=0. } 여기서 τ = τ w = 상수 {\ displaystyle \tau =\tau _{w}={\text{constant}}} 는 벽 전단 응력이다. The flow does not depend on the Reynolds number R e = U l / ν ∞ {\displaystyle \mathrm {Re} =Ul/\nu _{\infty }} , but rather on the Prandtl number P r = μ ∞ c p ∞ / κ ∞ {\displaystyle \mathrm {Pr} =\mu _{\infty }c_{p\infty }/\kappa _{\infty }} and the Mach number M = U / c ∞ = U / ( γ − 1 ) h ∞{\ displaystyle \ mathrm{M} =U/c_{\infit }=U/{\sqrt{(\gamma -1)h_{\ infit }}}}}, 여기서 where{\displaystystyle c} 는 음속 의 속도 , γ {\displaystytype \gamma. 비차원 변수 소개
y ~ = y l , T ~ = T T ∞ , T ~ w = T w T ∞ h ~ = h h ∞ h ~ w = h w h ∞ u ~ = u U μ ~ = μ μ ∞ τ ~ w = τ w μ ∞ U / l displaystyle {\y}={\frac}}}{\frac {l}},\quad {\tilde{\}}}}} T}={\frac {T}{\frac}}{\frac}}{{\frac}}}{{{p T_{\infty }}},\quad {\tilde {T}}_{w}={\frac {T_{w}}{T_{\infty }}},\quad {\tilde {h}}={\frac {h}{h_{\infty }}},\quad {\tilde {h}}_{w}={\frac {h_{w}}{h_{\infty }}},\quad {\tilde {u}}={\frac {u}{ U},\quad {\tilde{\mu }}={\frac {\mu }{\mu }{\mu }},\quad {\tilde {\\tau }}}={\frac {\w}{\mu }{\fl}}}}}}}}}}}}} 이러한 수량으로 볼 때 해결책은
h ~ = h ~ w + [ γ − 1 2 M 2 P r + ( 1 − h ~ w ) ] u ~ − γ − 1 2 M 2 P r u ~ 2 {\displaystyle {\tilde {h}}={\tilde {h}}_{w}+\left[{\frac {\gamma -1}{2}}\mathrm {M} ^{2}\mathrm {Pr} +(1-{\tilde {h}}_{w})\right]{\tild e {u}-{{\frac {\gamma -1}{2}}\mathrm {M}^{2}\mathrm {Pr} \,{\tilde {u}^{2},} y ~ = 1 τ ~ w ∫ 0 u ~ μ ~ d u ~ τ ~ w = ∫ 0 1 μ ~ d u ~ q w = − 1 P r τ w ( d h d u ) w {\displaystyle {\tilde {y}}={\frac {1}{{\tilde {\tau }}_{w}}}\int _{0}^{\tilde {u}}{\t ilde {\mu }}d{\tilde {u}},\quad {\tilde {\tau }}_{w}=\int _{0}^{1}{\tilde {\mu }}d{\tilde {u}},\quad q_{w}=-{\frac {1}{\mathrm {Pr} }}\tau _{w}\left({\frac {dh}{du}}\right)_{w},} 여기서 q w {\ displaystyle q_{w}} 는 하부 벽에서 단위 면적당 단위 시간 당 전달되는 열이다 . Thus h ~ , T ~ , u ~ , μ ~ {\displaystyle {\tilde {h}},{\tilde {T}},{\tilde {u}},{\tilde {\mu }}} are implicit functions of y {\displaystyle y} . One can also write the solution in terms of the recovery temperature T r {\displaystyle T_{r}} and recovery enthalpy h r {\displaystyle h_{r}} 절연벽의 온도에서 평가됨. 즉 , q = 0 {\displaystyle T_{w}} 및 h 값 {\ displaystyle h_{w }}.[clarification needed ] 그렇다면 해결책은
q w τ w U = T ~ w − T ~ r ( γ − 1 ) M 2 P r T ~ r = 1 + γ − 1 2 M 2 P r w}{w}{ww}}{w}}{\ {q_{w}}{\tau_{w}}} U}}={\frac {{\tilde {T}}_{w}-{\tilde {T}}_{r}}{(\gamma -1)\mathrm {M} ^{2}\mathrm {Pr} }},\quad {\tilde {T}}_{r}=1+{\frac {\gamma -1}{2}}\mathrm {M} ^{2}\mathrm {Pr} ,} h ~ = h ~ w + ( h ~ r − h ~ w ) u ~ − γ − 1 2 M 2 P r u ~ 2 . {\displaystyle {\tilde {h}}={\tilde {h}}_{w}+({\tilde {h}}_{r}-{\tilde {h}}_{w}){\tilde {u}}-{\frac {\gamma -1}{2}}\mathrm {M} ^{2}\mathrm {Pr} \,{\tilde {u}}^{2}. } 특정 열 이 일정하면 h ~ = T ~ {\ displaystyle {\tilde{h}={\tilde{{\t }}. M → 0 {\displaystyle \mathrm {M} \rightarrow 0} 및 T = t = q q = 0 {\displaystystytylease T_T} T_{\put },\Rightarrow q_{w}=0 }, T {\displaystyle T}, μ {\displaystyle \mu } 이 (가) 도처에 일정하게 나타나므로 , 압축 불가능한 쿠엣 플로우 용액을 회수한다. Otherwise, one must know the full temperature dependence of μ ~ ( T ~ ) {\displaystyle {\tilde {\mu }}({\tilde {T}})} . While there is no simple expression for μ ~ ( T ~ ) {\displaystyle {\tilde {\mu }}({\tilde {T}})} that is both accurate and general, there are several approximations for certain materials — 예 를 들어, 점도의 온도 의존성. When M → 0 {\displaystyle \mathrm {M} \rightarrow 0} and q w ≠ 0 {\displaystyle q_{w}\neq 0} , the recovery quantities become unity T ~ r = 1 {\displaystyle {\tilde {T}}_{r}=1} . For air, the values γ = 1.4 , μ ~ ( T ~ ) = T ~ 2 / 3 {\displaystyle \gamma =1.4,\ {\tilde {\mu }}({ \tilde{T}}={\tilde{T}^{2/3}}}} 이(가) 일반적으로 사용되며 , 이 경우에 대한 결과는 그림에 나와 있다.
분리 및 이온화 의 영향(즉 , c p {\ displaystyle c_{p }) 도 연구되었다. 이 경우 분자의 분리에 의해 회수 온도가 감소한다.[7]
직사각형 채널 1 couette flow with h/l=0.1 1차원 flow u( y ) {\displaystyle u(y)} 은 두 플레이트가 스트림 방향(x {\displaystyle x }) 및 스팬 방향(z {\displaysty z}) 에서 무한히 긴 경우에 유효하다 . 스팬웨이즈 길이가 유한하면 흐름이 2차원이 되고 u {\displaystyle u} 은 x {\displaystyle x} 와 z {\displaystyle z} 의 함수다. 다만 흐름의 단방향성을 확보하기 위해서는 스트림 방향의 무한 길이를 유지해야 한다.
예를 들어, 상단 벽이 일정한 속도 U {\displaystyle U} 을(를) 사용하여 이동하는 조건에 따라 가로 높이 h {\displaystyle h} 및 세로 폭 l {\displaystyle l} 을(를) 가진 무한히 긴 직사각형 채널을 고려해 보십시오. 가압 경사 없이 Navier –스토크 방정식은 다음으로 감소한다.
∂ 2 u ∂ y 2 + ∂ 2 u ∂ z 2 = 0 {\displaystyle {\frac {\fract ^{2}u}{\frac y^{2}}+{\frac {\frac ^{2}u}{\flash z^{2}}=0} .
u ( 0 z ) = 0 u ( h z ) = U {\displaystyle u(0,z)=0,\quad u(h,z)=U,} u ( y 0 ) = 0 u ( y l ) = 0. {\displaystyle u(y,0)=0,\display u(y,l)=0.} 변수의 분리 를 사용하여 용액은 다음과 같다.
u ( y z ) = 4 U π ∑ n = 1 ∞ 1 2 n − 1 징징거리다 ( β n y ) 징징거리다 ( β n h ) 죄를 짓다 ( β n z ) β n = ( 2 n − 1 ) π l . {\displaystyle u(y,z)={\frac {4} U}{\pi }}\sum _{n=1}^{\infty }{\frac {1}{2n-1}}{\frac {\sinh(\beta _{n}y)}{\sinh(\beta _{n}h)}}\sin(\beta _{n}z),\quad \beta _{n}={\frac {(2n-1)\pi }{l}}. } h / l ≪ 1 {\displaystyle h/l\ll 1}이( 가) 표시되면 그림과 같이 평면 쿠엣 흐름이 복구된다.
실린더 Taylor-Couette 흐름 은 무한히 긴 두 회전하는 동축 실린더 사이의 흐름이다.[8] 원래의 문제는 1845년 스톡스 에 의해 해결되었지만,[9] 제프리 잉그램 테일러 는 1923년 유명한 논문에서 그 안정성을 연구했기 때문에 그 흐름에 이름이 붙었다.[10]
The problem can be solved in cylindrical coordinates ( r , θ , z ) {\displaystyle (r,\theta ,z)} . Denote the radii of the inner and outer cylinders as R 1 {\displaystyle R_{1}} and R 2 {\displaystyle R_{2}} . Assuming the cylinders rotate at constant angular velocities Ω 1 {\displaystyle \Omega _{1}} and Ω 2 {\ displaystyle \Oomega _{2 }, 그러면 θ {\displaystyle \theta} -방향의 속도는[11]
v θ ( r ) = a r + b r a = Ω 2 R 2 2 − Ω 1 R 1 2 R 2 2 − R 1 2 b = ( Ω 1 − Ω 2 ) R 1 2 R 2 2 R 2 2 − R 1 2 . {\displaystyle v_{\theta }(r)=ar+{\frac {b}{r}},\qquad a={\frac {\Omega _{2}R_{2}^{2}-\Omega _{1}R_{1}^{2}}{R_{2}^{2}-R_{1}^{2}}},\quad b={\frac {(\Omega _{1}-\Omega _{2}) R_{1}^{2}R_{2}^{2}}: {R_{2}^{2}-R_{1}^{2}}. } 이 방정식은 곡률의 효과가 더 이상 흐름 영역에서 일정한 전단(shear)을 허용하지 않음을 보여준다.
길이가 유한한 동축 실린더 기존의 테일러-쿠엣 흐름 문제는 무한히 긴 실린더를 가정한다. 실린더에 불가결한 유한 길이 l {\displaystyle l} 이( 가) 있는 경우, 분석을 수정해야 한다(흐름이 여전히 단방향이지만). Ω 2 = 0 {\displaystyle \Oomega _{2}=0 } 의 경우 유한 길이 문제는 변수 또는 적분 변환 의 분리를 사용 하여 해결할 수 있으며,[12] 다음과 같은 장점이 있다.
v θ ( r , z ) = 4 R 1 Ω 1 π ∑ n = 1 ∞ 1 2 n − 1 I 1 ( β n R 2 ) K 1 ( β n r ) − K 1 ( β n R 2 ) I 1 ( β n r ) I 1 ( β n R 2 ) K 1 ( β n R 1 ) − K 1 ( β n R 2 ) I 1 ( β n R 1 ) 죄를 짓다 ( β n z ) , β n = ( 2 n − 1 ) π l , {\displaystyle v_{\theta }(r,z)={\frac {4R_{1}\Omega _{1}}{\pi }}\sum _{n=1}^{\infty }{\frac {1}{2n-1}}{\frac {I_{1}(\beta _{n}R_{2})K_{1}(\beta _{n}r)-K_{1}(\beta _{n}R_{2}) I_{1}(\beta _{n}r)}{ I_{1}(\beta _{n}R_{2})K_{1}(\beta _{n}R_{1}-K_{1}(\beta _{n}R_{2}) I_{1}(\beta _{n}R_{1}}})}\sin(\beta _{n}z),\quad \beta _{n}={\frac {(2n-1)\pi },},} 여기서 I( β n r ) , K ( β n r ) {\displaystyle I(\beta _{n}r),\ K(\beta _{n}r) 는 제1종과 제2종의 수정된 베셀 함수 다.
참고 항목 참조 ^ 질렌코 외(2018년) ^ 가이온 외 (2001), 페이지 136 ^ 헬러(1960년) ^ 포즈리키디스(2011), 페이지 338–339 ^ 쿤두 외 (2016), 페이지 415 ^ 라거스트롬(1996) ^ Lippmann 외 연구진(1956, 1957) ^ 란다우와 리프시츠 (1987년) ^ 스톡스 (1845) ^ 테일러(1923년) ^ 가이온 외 (2001년), 페이지 163–197 ^ 웬들 (1999년)
원천 Acheson, D.J. (1990). Elementary Fluid Dynamics . Oxford University Press. ISBN 0-19-859679-0 . Batchelor, G.K. (2000) [1967]. An Introduction to Fluid Dynamics . Cambridge University Press. ISBN 0-521-66396-2 . Guyon, Etienne; Hulin, Jean-Pierre; Petit, Luc; Mitescu, Catalin D. (2001). Physical Hydrodynamics . Oxford University Press. ISBN 0-19-851746-7 . Heller, John P. (1960). "An Unmixing Demonstration". American Journal of Physics . 28 (4): 348–353. Bibcode :1960AmJPh..28..348H . doi :10.1119/1.1935802 . ISSN 0002-9505 . Illingworth, C. R. (1950). "Some solutions of the equations of flow of a viscous compressible fluid". Mathematical Proceedings of the Cambridge Philosophical Society . 46 (3): 469–478. Bibcode :1950PCPS...46..469I . doi :10.1017/S0305004100025986 . ISSN 0305-0041 . Kundu, Pijush K.; Cohen, Ira M.; Dowling, David R. (2016). Fluid Mechanics (6th ed.). Elsevier. ISBN 978-0-12-405935-1 . Lagerstrom, Paco (1996). Laminar flow theory . Princeton University Press. ISBN 978-0691025988 . Landau, L. D.; Lifshitz, E.M. (1987). Fluid Mechanics (2nd ed.). Elsevier. ISBN 978-0-08-057073-0 . Lippmann, H. W., Z. O. Bleeviss. "분열과 이온화가 압축성 쿠엣 흐름에 미치는 영향." 더글러스 항공 파충류. SM-19831 130(1956년). 리프만, 한스 볼프강 , 아나톨 로슈코 . 가스 역학의 요소들. 1957년 택배회사 Pozrikidis, C. (2011). Introduction to Theoretical and Computational Fluid Dynamics . Oxford University Press. ISBN 978-0-19-975207-2 . 리처드 파인만 (1964) 파인만 물리학 강의: 주로 전자기 및 물질 , § 41–6 "쿠엣 흐름", 애디슨–웨슬리 ISBN 0-201-02117-X Stokes, George Gabriel (1880). "On the Theories of the Internal Friction of Fluids in Motion, and of the Equilibrium and Motion of Elastic Solids" . Mathematical and Physical Papers . Cambridge University Press: 75–129. doi :10.1017/CBO9780511702242.005 . ISBN 9780511702242 . Taylor, Geoffrey I. (1923). "Stability of a viscous liquid contained between two rotating cylinders" . Philosophical Transactions of the Royal Society of London . Series A, Containing Papers of a Mathematical or Physical Character. 223 (605–615): 289–343. Bibcode :1923RSPTA.223..289T . doi :10.1098/rsta.1923.0008 . JSTOR 91148 . Wendl, Michael C. (1999). "General solution for the Couette flow profile". Physical Review E . 60 (5): 6192–6194. Bibcode :1999PhRvE..60.6192W . doi :10.1103/PhysRevE.60.6192 . ISSN 1063-651X . PMID 11970531 . Zhilenko, Dmitry; Krivonosova, Olga; Gritsevich, Maria; Read, Peter (2018). "Wave number selection in the presence of noise: Experimental results". Chaos: An Interdisciplinary Journal of Nonlinear Science . 28 (5): 053110. Bibcode :2018Chaos..28e3110Z . doi :10.1063/1.5011349 . hdl :10138/240787 . ISSN 1054-1500 . PMID 29857673 . 외부 링크