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

We consider a viscous, incompressible fluid whose motion is subject to a random force. The container of the fluid is supposed to be the torus T

N/A
N/A
Protected

Academic year: 2021

シェア "We consider a viscous, incompressible fluid whose motion is subject to a random force. The container of the fluid is supposed to be the torus T"

Copied!
3
0
0

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

全文

(1)

Stochastic Power Law Fluid Equations

寺澤祐高

(Yutaka Terasawa),

東大数理、日本学術振興会特別研究員

PD

吉田伸生氏(京大・理)との共同研究

We consider a viscous, incompressible fluid whose motion is subject to a random force. The container of the fluid is supposed to be the torus T

d

= ( R / Z )

d

= [0, 1)

d

as a part of idealization. For a differentiable vector field v : T

d

R

d

, which is interpreted as the velocity field of the fluid, we denote the rate of strain tensor by:

e(v) =

(

i

v

j

+

j

v

i

2

)

: T

d

R

d

R

d

. (0.1) We assume that the extra stress tensor

τ (v) : T

d

R

d

R

d

depends on e(v) polynomially. More precisely, for ν > 0 (the kinematic viscosity) and p > 1,

τ (v) = 2ν(1 + | e(v) |

2

)

p−22

e(v). (0.2) Given an initial velocity u

0

: T

d

R

d

, the dynamics of the fluid is described by the following SPDE:

div u = 0, (0.3)

t

u + (u · ∇ )u = −∇ Π + div τ (u) +

t

W, (0.4) where

u · ∇ =

d

j=1

u

j

j

and div τ(u) = (

d

j=1

j

τ

ij

(u) )

d

i=1

. (0.5)

The unknown process in the SPDE are the velocity field u = u(t, x) = (u

i

(t, x))

di=1

and the pressure Π = Π(t, x). The Brownian motion W = W (t, x) = (W

i

(t, x))

di=1

with values in L

2

( T

d

R

d

) (the set of vector fields on T

d

with L

2

components) is added as the random force. Physical interpretation of (0.3) and (0.4) are the mass conservation, and the motion equation, respectively. We note that the SPDE (0.3)–(0.4) for the case p = 2 is the stochastic Navier-Stokes equation [1, 2].

Our motivation comes from works by J. M´ alek, J. Ne˘ cas, M. Rokyta, and M.

R u˘

zi˘ cka [3], where the deterministic equation (where the colored noise

t

W in (0.3)–(0.4) is dropped) is investigated.

We need the following definition.

Definition 1 Let H be a Hilbert space, and Γ : H H be a self-adjoint, non- negative definite operator of trace class. A random variable (W

t

)

t≥0

with values in C([0, ) H) is called a H-valued Brownian motion with the covariance operator Γ (abbreviated by BM(H, Γ) below) if, for each ϕ H and 0 s < t,

E [exp (i h ϕ, W

t

W

s

i ) | (W

u

)

us

] = exp (

t s

2 h ϕ, Γϕ i )

, a.s.

(2)

We also need some function spaces. Let V be the set of R

d

-valued divergence free, mean-zero trigonometric polynomials. We equip the torus T

d

with the Lebesgue measure. For p [1, ) and α R , we introduce:

V

p,α

= the completion of V with respect to the norm k · k

p,α

, (0.6) where

k v k

pp,α

=

Td

| (1 ∆)

α/2

v |

p

. (0.7) Definition 2 Suppose that

I Γ : V

2,0

V

2,0

is a bounded self-adjoint, non-negative definite operator of trace class;

I µ

0

is a Borel probability measure on V

2,0

.

I (X, Y ) = ((X

t

, Y

t

))

t0

is a process defined on a probability space (Ω, F , P ) such that:

X L

p,loc

([0, ) V

p,1

) L

,loc

([0, ) V

2,0

) C([0, ) V

2p0,−β

), (0.8) for some β > 0, and (Y

t

)

t≥0

is a BM(V

2,0

, Γ) (cf. Definition 1).

Then, the process (X, Y ) is said to be a weak solution to the SDE (stochastic differential equation)

X

t

= X

0

+

t 0

b(X

s

)ds + Y

t

(0.9)

with the initial law µ

0

if the following conditions are satisfied;

P (X

0

∈ · ) = µ

0

;

Y

t+·

Y

t

and {h ϕ, X

s

i ; s t, ϕ ∈ V} are independent for any t 0;

h ϕ, X

t

i = h ϕ, X

0

i +

t 0

h ϕ, b(X

s

) i ds + h ϕ, Y

t

i , for all ϕ ∈ V and t 0.

We can now state our existence result.

Theorem 3 Let Γ and µ

0

be as in Definition 2 and suppose additionally that p (

32

, ) (d = 2), p (

95

, ) (d = 3) · · · ,

I ∆Γ = Γ∆ and both Γ, ∆Γ are of trace class;

I µ

0

is a probability measure on V

2,1

and m

α

=

k ξ k

22,α

µ

0

(dξ) < for α = 0, 1. (0.10)

hold. Then, there exists a weak solution to the SDE (0.9) with the initial law µ

0

(cf. Definition 2) such that (0.8) holds with β = β(p) > 0. Moreover, for any T > 0,

E [

sup

t≤T

k X

t

k

22

+

T 0

k X

t

k

pp,1

dt ]

(1 + T )C < , (0.11)

where C = C(d, p, Γ, m

0

) < .

(3)

The proof is done by using Galerkin approximation. We obtain Stochastic ODEs by this. By combining the deterministic method from [3] and Ito’s formula, we get their uniform apriori estimates needed for passing to the limit.

References

[1] Flandoli, Franco : An introduction to 3D stochastic fluid dynamics. SPDE in hy- drodynamic: recent progress and prospects, 51–150, Lecture Notes in Math., 1942, Springer, Berlin, 2008.

[2] Flandoli, Franco ; Gatarek, Dariusz: Martingale and stationary solutions for stochastic Navier-Stokes equations. Probab. Theory Related Fields 102 (1995), no.

3, 367–391.

[3] M´ alek, J.; Ne˘ cas, J.; Rokyta, M.; R u˘

zi˘ cka, M.: Weak and measure-valued solutions to evolutionary PDEs. Applied Mathematics and Mathematical Computation, 13.

Chapman & Hall, London, 1996. xii+317 pp. ISBN: 0-412-57750-X

[4] Yutaka Terasawa and Nobuo Yoshida, Stochastic Power Law Fluids: Existence and

Uniqueness of Weak Solutions, Annals of Applied Probability 21 (2011), 1827–1859.

参照

関連したドキュメント

The analysis of the displacement fields in elastic composite media can be applied to solve the problem of the slow deformation of an incompressible homogen- eous viscous

We consider the Cauchy problem for nonstationary 1D flow of a compressible viscous and heat-conducting micropolar fluid, assuming that it is in the thermodynamical sense perfect

This article is devoted to establishing the global existence and uniqueness of a mild solution of the modified Navier-Stokes equations with a small initial data in the critical

Wheeler, “A splitting method using discontinuous Galerkin for the transient incompressible Navier-Stokes equations,” Mathematical Modelling and Numerical Analysis, vol. Schotzau,

The numerical tests that we have done showed significant gain in computing time of this method in comparison with the usual Galerkin method and kept a comparable precision to this

In this section, we shall prove the following local existence and uniqueness of strong solutions to the Cauchy problem (1.1)..

Then (v, p), where p is the corresponding pressure, is the axisymmetric strong solution to problem (1.1) which is unique in the class of all weak solutions satisfying the

A mathematical formulation of well-posed initial boundary value problems for viscous incompressible fluid flow-through-bounded domain is described for the case where the values