• 検索結果がありません。

2成分層流の時間的大域解について(流体とプラズマの諸現象の数学解析)

N/A
N/A
Protected

Academic year: 2021

シェア "2成分層流の時間的大域解について(流体とプラズマの諸現象の数学解析)"

Copied!
4
0
0

読み込み中.... (全文を見る)

全文

(1)

2 成分層流の時間的大域解について

YASUSHI HATAYA 京大理幡谷 泰史 YOSHIAKI TERAMOTO 京大・理 寺本 恵昭 ABSTRACT. 2成分流体の平面 Couette 流について考察する。粘性が異なる場合について、 適当な仮定の下で大域解が存在し、指数的に減衰することを示す。

1.

INTRODUCTION

上壁は水平方向に–定速度で動き、

下壁は固定されている

2

枚の平行平板間の

2

成分流 体を考える。

2

流体は粘性のみ異なるものとし、 それ自身未知関数である自由界面で分けら

れているものとする。それぞれの流体の運動は

Navier-Stokes

方程式、界面では流速の連

続性・応力と表面張力のつりあい・運動力学的境界条件、壁では粘着境界条件によって記述

されているものとする。

[1]

における無次元化を採用すれば方程式、境界条件は次のように

なる。 方程式

$\mathrm{u}_{t}^{I}+(\mathrm{u}^{I}\cdot\nabla)\mathrm{u}^{I}+\nabla p^{I}=\mu_{1}.\Delta \mathrm{u}^{I}$

,

$t>0,0<x_{2}<l+h(t, x_{1})$

$\nabla\cdot \mathrm{u}^{I}=0$

,

$\mathrm{u}_{t}^{II}+,(\mathrm{u}^{II}\cdot\nabla)\mathrm{u}^{IIII}+\nabla p=\mu_{2}\Delta \mathrm{u}^{II}$

,

$t>0,$ $l+h(t,x_{1})<x_{2}<1$

.

$\nabla\cdot \mathrm{u}^{II}=0$

,

初期境界条件 $\mathrm{u}^{j}|_{t=0}=\mathrm{u}_{0}^{j}$

,

$h|_{t=0}=h_{0}$ $\mathrm{u}^{I}|_{x_{2}=0}=(0,0)$

,

$\mathrm{u}^{II}|_{x_{2}1}==(1,0)$ . $h_{t}+u_{11}^{I}h_{x}=u_{2}^{I}$ $u_{1}^{I}=u_{1}^{II}$

,

$u_{2}^{I}=u_{2}^{II}$ デ.$T^{I}\vec{n}=\tilde{\tau}\cdot T^{II}\tilde{n}$

$\vec{n}\cdot T^{I}\vec{n}-\vec{n}\cdot T^{II}\vec{n}=\frac{Sh_{x_{1^{X}1}}}{(1+(hx_{1})2)\frac{3}{2}}$

.

ここで、水平方向に $L-$ 周期的、

$hd$

r=0

、正数$S$ は表面張力係数、 $\vec{n},\vec{\tau}$ はそれぞ

れ界面における法ベクトル、接ベクトルとする。

数理解析研究所講究録

(2)

以下では次のような定常解

(COuette 流) の周りで考える。 定常解

$u_{1}(x_{1}, x_{2},.t)=U_{0}(x_{2})=\{$ $\frac{\frac{x_{2}}{l+m(1-l)m(x2-1)}}{l+m(1-^{\iota})},+1$

,

$\mathrm{i}\mathrm{n}l<X<1\mathrm{i}\mathrm{n}0<x_{2}2<\iota$

,

$u_{2}^{J}\equiv 0$

,

$h\equiv 0$

,

$p^{I}\equiv p^{II}\equiv p0$

(const).

ここで、 $m=\mu_{\underline{1}}$ とした。この

2

成分流体の線形安定性は、 いくつかの特殊な極限につい て詳しく調べられている

(Joseph

and Renardy [1])

。また、

2

成分流体の解の存在を扱っ た論文として

Tanaka

[3]

などがある。

2.

LOCAL

SOLUTION

自由界面を持つ領域の問題から、

定常解が占める固定領域へ問題を書き換える。

参照領

域への変換写像を次の $\ominus$ の逆写像として定義する。 $\tilde{h}(x_{1}, x_{2}, t):=\sum_{\xi\in \mathbb{Z}^{\frac{\hat{h}(\xi,x_{2},t)e^{x_{1}\epsilon}}{1+\xi^{2}(x_{2}-\iota)^{2}}}}$ とお

くと $(x_{1}, x_{2}, t)\in\Omega(t)\cross(\mathrm{O}, T)$

,

$(x_{1}’, x_{2}, t\prime\prime)\in\Omega_{0}\cross(0, T)$ に対して

$\ominus:\Omega_{i}^{0}\cross(0,\tau)arrow\Omega_{i}(t)\cross(0,T)$

,

$x_{1}=x_{1}’$

,

$x_{2}=X_{2}’+ \frac{x_{2}’(_{X_{2}-}/1)}{l(l-1)}\tilde{h}(x_{1’ 2};’/x,t)$

,

$t=t’$

.

この変換は $h(x_{1}, t)$

が十分小さければ微分同相写像を定めるので、 以下では固定領域で解を

構成しよう。方程式境界条件は次のように書き換わる。

$\mathrm{u}_{t}^{I}-\mu_{1}\Delta \mathrm{u}^{I}+\nabla p^{I}=f(\mathrm{u}^{I}, h)+b(h)\nabla p^{I}$

,

. $\nabla\cdot \mathrm{u}^{I}=0$

,

$t>0,0<X_{2}<l$

.$\mathrm{u}_{t}^{II}-\mu 2\Delta.\mathrm{u}^{II.I}+\nabla p=If(\mathrm{u}, hII)+b(h)\nabla p^{I}I$

,

$\nabla\cdot \mathrm{u}^{II}=0$

,

(2.1)

$t>0,$ $l<x_{2}<1$ $\mathrm{u}^{I}|_{x_{2}=0}=(0,0)$

,

$\mathrm{u}^{II}|_{x_{2}}=1=(\mathrm{o}, 0)$ $h_{t}-u_{2}^{I}=\partial x_{1}- a\mathrm{o}(h)$

(2.2)

$(u_{1}^{I}-u_{1})II+’. \frac{1-mj}{l+m(1-l)}h=a_{1}(h)$ $(u_{2}^{I}-u^{II})2x1a=\partial 2(h)$ $\mu_{1}(u_{2,x1}^{I}+uI1,x_{2})-\mu 2(u_{2,1,x}+u^{II})II=a_{3}(x_{1}2\mathrm{u}, h)$ $(2\mu 1u_{2,x2}^{I}-p^{I})-(2\mu_{2}u^{II}p^{II}2,x2^{-})-shx_{1}x1=a_{4}(\mathrm{u}, h)$

.

VASOlonnikov

[6]

の手法を用いて、次の局所解を得ることができる。

120

(3)

$\mathrm{p}\mathrm{r}\mathrm{o}\mathrm{p}_{\mathrm{o}\mathrm{S}}\mathrm{i}\mathrm{t}\mathrm{i}\mathrm{o}\mathrm{n}2.1$

.

$\frac{3}{2}<r<2$ とする。 $\epsilon>0$ が存在して $|m-1|<\epsilon$ ならば、ある適合

条件を満たす初期条件

$(u_{0}, h\mathrm{o})$ と任意に与えられた時間 $T>0$ に対して、 正数$\delta$ が存在し

て、

$|h_{0}|_{H^{r-}}1/2(\Omega^{0})+|u_{0}|_{H\mathrm{t}^{\mathrm{R}})}n-1\leq\delta$

ならば、初期値境界値問題 $(2.1)_{f}(2.2)$ はつぎのような解 $(h, u,p)$ を持つ。

$u\in H^{r,\frac{f}{2}}(\Omega_{0}\cross(0, T))$

,

$h\in H^{r+\frac{1}{2},\frac{\tau+1/2}{2}}(\Gamma_{0}\cross(0, T))$

,

$\nabla p\in H^{r-2}’\frac{\tau-2}{2}(\Omega_{0}\cross(0, T))$

,

$p|_{\Gamma}\in H^{r-\frac{3}{2},\frac{r-\theta/2}{2}}(\mathrm{r}_{0}\cross(0, T))$

.

3. GLOBAL SOLUTION

TheOrem 3.1.

正数 $\epsilon_{1}$ に対して、 $|m-1|<\epsilon_{1}$ であるとする。 ある正数 $\gamma$

が存在して、

任意の$t_{0}>0$ に対して正数 $M$ が存在して、任意の$T>0$ に対して局所解 $(u, h,p)$ が

$(3.1)$ $\sup$ $\{||u(t)||_{H^{2}(\Omega_{0}^{j})}+|h(t)|_{H^{5/}(\Gamma)}20+|p(t)|_{H^{1/2}}(\mathrm{r}0)\}\leq\epsilon_{1}$

$t_{0}\leq t\leq t\mathit{0}+T$ を満たすのならば、

(3.2)

$||u(t)||_{2}+|h(t)|_{3}\leq Me^{-\gamma t}\{||u(t\mathrm{o})||_{2}+|h(t_{0})|_{3}\}$

,

$1\leq t\leq T$

,

が成り立つ。

証明の方針

$E= \sum_{j=0,12},||\partial_{x}^{j}\mathrm{u}(t)||^{2}+||\mathrm{u}_{t}(t)||^{2}+\frac{S}{2}(\sum_{j=0,12},|h_{x}(t)|^{2}+|h_{xt}(t)|^{2})$

,

$F= \sum$ $||\nabla\partial_{x}^{j}\mathrm{u}(t)||^{2}+||\nabla \mathrm{u}_{t}(t)||2$

,

$j=0,1,2$

$G= \int_{\Gamma}hxuxu_{1},z+\int_{\Gamma}(hh_{x})_{t}p+\frac{m-1}{\mathcal{R}_{2}(l+m(1-l))}\int_{\Gamma}hu_{1,z}+\cdots$

.

とおく。 $G$ について、

(3.1)

と $|m-1|<\epsilon_{1}$ を用いて

(3.3)

$|G(t)|\leq c\epsilon_{1}E(t)$

,

を満たすように$\epsilon_{1}$ を十分小さくとれる。 また、ある $\gamma’>0$がとれて、 $\frac{d}{dt}(E+G)+\gamma^{l}F\leq 0$

,

を導くことができる。

POincare’ の不等式を用いれば

$E+G\leq \mathrm{c}F$ であるので、

(3.4)

$\frac{d}{dt}(E+G)+\gamma(E+G)\leq 0$

,

for

$\gamma>0$

.

また、垂直方向の微分についての評価は、

(4)

Proposition 3.2.

定常

Stokes

方程式 $-\mu_{j}\Delta \mathrm{u}+\nabla p=f$

in

$\Omega_{j}$ $\nabla u=0$

(3.5)

$u=0$

on

$S_{B}$ $u_{2}=u_{2}$ $u_{1,x_{2}}+u_{2,x_{1}}=u_{1,x_{2}}+u_{2,x_{1}}$

on

$x_{2}=l$

,

の解 $(u,p)$ について、次の評価が成り立つ。 $||u||s+2^{+}||\nabla p||_{s}\leq C||f|_{s}$ $s=0,1$

.

この

Proposition

(2.1), (2.2)

を用いれば、 次の評価をうる。

(3.6)

$\sum_{s=0,1}(||\partial_{x_{2}}S\mathrm{u}||2\nabla+||p||_{s})\leq cF$

.

(3.3), (3.4), (3.6)

をあわせれば、

(3.1)

を得る。

REFERENCES

1. D. D.JOSEPH AND Y.RENARDY, Fundamentals

of two-fluid

$dynamiCs:Mathematical$theory and its

applications, Springer, New York, IAM series 3, 1992

2. C. S.YIH, Instability due to viscosity strtification, J.Fluid Mech.26,337 (1967).

3. N.TANAKA

,

$Globa\iota$ existence

of

two phase homogeneous viscous incompressible

fluid

flow,Comm.P.D.E.,vo118

,

(1&2), pp.41-81, (1993)

4. T.NISHIDA AND Y.TERAMOTO AND HtayAung Win, Navier-Stokes

flow

down an inclined plane:

Downward periodic motion, J.Math.Kyoto.Vniv. vol.33, $\mathrm{N}\mathrm{o}.3,(1993)$

.

5. J. T.BEALE, Large time behavior

of

viscous

surface

waves, Arch.Rat.Mech.Anal. 84, pp.304-352

(1984).

6. V. A.SOLONNIKOV, On aninitial-boundary value problem

for

theStoke systems arising in the study

of

problem with a

free

boundary,Proc.Steklov.Inst.Math. 3, pp.191-239 (1991).

7. J. L.LIONS AND E.MAGENES, Nonhomogeneous boundary value problem and its applications,

Springer, NewYork, 1969

8. O. $\mathrm{A}.\mathrm{L}\mathrm{A}\mathrm{D}\mathrm{Y}\check{\mathrm{Z}}\mathrm{E}\mathrm{N}\mathrm{S}\mathrm{K}\mathrm{A}\mathrm{j}\mathrm{A}$,V. A.SOLONNIKOV

AND N.

N.URALCS’VA,

Linearand quasilinear equations

of

parabolic type, Translations ofMathematical Monographs$\mathrm{V}\mathrm{o}\mathrm{l}23,\mathrm{A}\mathrm{M}\mathrm{S}$(1968).

9. S. UKAI, A Solution Formula

for

the Stokes Equation in$R_{+}^{n}$

,

Comm.Pure.Apll.Math.$\mathrm{X}\mathrm{L}$

,

pp.611-621 (1987)

10. R. S. ADAMS, Sobolev Spaces, Academic Press, NewYork, 1975

11. M. S. AGRANOVICHAND M. S. VISHIK, Elliptic problems with aparameter and parabolec problems

of

general type,Russian Math. Surveys,19, pp.53-161 (1964)

DEPARTMENTOF MATHEMATICS, KYOTO UNIVERSITY, KYOTO

参照

関連したドキュメント

2 解析手法 2.1 解析手法の概要 本研究で用いる個別要素法は計算負担が大きく,山

そのため本研究では,数理的解析手法の一つである サポートベクタマシン 2) (Support Vector

金沢大学は学部,大学院ともに,人間社会学分野,理工学分野,医薬保健学分野の三領域体制を

振動流中および一様 流中に没水 した小口径の直立 円柱周辺の3次 元流体場 に関する数値解析 を行った.円 柱高 さの違いに よる流況および底面せん断力

非自明な和として分解できない結び目を 素な結び目 と いう... 定理 (

Supersingular abelian varieties and curves, and their moduli spaces 11:10 – 12:10 Tomoyoshi Ibukiyama (Osaka University).. Supersingular loci of low dimensions and parahoric subgroups

3 Numerical simulation for the mteraction analysis between fluid and

Mochizuki, Topics Surrounding the Combinatorial Anabelian Geometry of Hyperbolic Curves III: Tripods and Tempered Fundamental Groups, RIMS Preprint 1763 (November 2012).