Pressure estimate for Navier-Stokes equation in bounded domains ...

6 downloads 3980 Views 601KB Size Report
Helmholtz projection onto divergence-free vector fields. L. 2. (Ω,R ..... Stability analysis: dot with −∆ u n+1 u .... Theorem Assume Ω is a bounded domain in R. N.
1

Pressure estimate for Navier-Stokes equation in bounded domains Jian-Guo Liu & Jie Liu (U Maryland College Park) Bob Pego (Carnegie Mellon) • Stability and convergence of efficient Navier-Stokes solvers via a commutator estimate, to appear in Comm. Pure Appl. Math. • Error estimates for finite element Navier-Stokes solvers with explicit time stepping for pressure (submitted)

2

~ut + ~u · ∇~u + ∇p = ν∆~u + f

in Ω

∇ · ~u = 0

in Ω

~u = 0

on Γ = ∂Ω

Phenomena modeled by Navier-Stokes dynamics: Lift, drag, boundary-layer separation, vortex shedding, . . .

3

~ut + ~u · ∇~u + ∇p = ν∆~u + f

in Ω

∇ · ~u = 0

in Ω

XXX   X  XXX

~u = 0

on Γ = ∂Ω

Phenomena modeled by Navier-Stokes dynamics: Lift, drag, boundary-layer separation, vortex shedding, . . .

4

 Z Z = ν∆~ ~ut + ~u · ∇~u + Z∇p u+f Z

in Ω

∇ · ~u = 0

in Ω

XXX   X  XXX

XXX  XX X

~u = 0

on Γ = ∂Ω

Phenomena modeled by Navier-Stokes dynamics: Lift, drag, boundary-layer separation, vortex shedding, . . .

5

Helmholtz projection onto divergence-free vector fields L2(Ω, RN ) = PL2(Ω, RN ) ⊕ ∇H 1(Ω) Given ~v ∈ L2(Ω, RN ), there exists q ∈ H 1(Ω) so that ~v = P~v − ∇q satisfies

hP~v , ∇φi = h~v + ∇q, ∇φi = 0 for all φ ∈ H 1(Ω).

5

Helmholtz projection onto divergence-free vector fields L2(Ω, RN ) = PL2(Ω, RN ) ⊕ ∇H 1(Ω) Given ~v ∈ L2(Ω, RN ), there exists q ∈ H 1(Ω) so that ~v = P~v − ∇q satisfies

hP~v , ∇φi = h~v + ∇q, ∇φi = 0 for all φ ∈ H 1(Ω).

Then ∇ · (P~v ) = 0 Note so

in Ω,

~n · P~v = 0 on Γ.

P~v ∈ H(div; Ω) = {f~ ∈ L2(Ω, RN ) : ∇ · f~ ∈ L2}, ~n · P~v ∈ H −1/2(Γ)

by a standard trace theorem.

6

Traditional unconstrained formulation of NSE

~ut + P(~u · ∇~u − f~) = νP∆~u, • Formally

~u|Γ = 0

∂t(∇ · ~u) = 0

• Perform analysis and computation in spaces of divergence-free fields (unconstrained Stokes operator P∆u is incompletely dissipative). • Inf-Sup/LBB condition(Ladyzhenskaya-Babuˇska-Brezzi)

7

Alternative unconstrained formulation of NSE ~ut + P(~u · ∇~u − f~) = νP∆~u + ν∇(∇ · ~u), • Formally, ∂t(∇ · ~u) = ν∆(∇ · ~u),

~u|Γ = 0

∂n(∇ · ~u)|Γ = 0

7

Alternative unconstrained formulation of NSE ~ut + P(~u · ∇~u − f~) = νP∆~u + ν∇(∇ · ~u), • Formally, ∂t(∇ · ~u) = ν∆(∇ · ~u),

~u|Γ = 0

∂n(∇ · ~u)|Γ = 0

• Equivalent to a ’reduced’ formulation of Grubb & Solonnikov (1991)

7

Alternative unconstrained formulation of NSE ~ut + P(~u · ∇~u − f~) = νP∆~u + ν∇(∇ · ~u), • Formally, ∂t(∇ · ~u) = ν∆(∇ · ~u),

~u|Γ = 0

∂n(∇ · ~u)|Γ = 0

• Equivalent to a ’reduced’ formulation of Grubb & Solonnikov (1991) • Lemma For all ~u ∈ L2(Ω, RN ), ∇(∇ · ~u) = ∆(I−P)~u

7

Alternative unconstrained formulation of NSE ~ut + P(~u · ∇~u − f~) = νP∆~u + ν∇(∇ · ~u), • Formally, ∂t(∇ · ~u) = ν∆(∇ · ~u),

~u|Γ = 0

∂n(∇ · ~u)|Γ = 0

• Equivalent to a ’reduced’ formulation of Grubb & Solonnikov (1991) • Lemma For all ~u ∈ L2(Ω, RN ), ∇(∇ · ~u) = ∆(I−P)~u • NSE: ~ut +P(~u · ∇~u − f~) = ν∆~u +ν(P∆ − ∆P)~u,

~u|Γ = 0

7

Alternative unconstrained formulation of NSE ~ut + P(~u · ∇~u − f~) = νP∆~u + ν∇(∇ · ~u), • Formally, ∂t(∇ · ~u) = ν∆(∇ · ~u),

~u|Γ = 0

∂n(∇ · ~u)|Γ = 0

• Equivalent to a ’reduced’ formulation of Grubb & Solonnikov (1991) • Lemma For all ~u ∈ L2(Ω, RN ), ∇(∇ · ~u) = ∆(I−P)~u • NSE: ~ut +P(~u · ∇~u − f~) = ν∆~u +ν(P∆ − ∆P)~u,

~u|Γ = 0

• commutator/gradient: [∆, P]~u = (I − P)(∆~u − ∇∇ · ~u) := ∇pS

7

Alternative unconstrained formulation of NSE ~ut + P(~u · ∇~u − f~) = νP∆~u + ν∇(∇ · ~u), • Formally, ∂t(∇ · ~u) = ν∆(∇ · ~u),

~u|Γ = 0

∂n(∇ · ~u)|Γ = 0

• Equivalent to a ’reduced’ formulation of Grubb & Solonnikov (1991) • Lemma For all ~u ∈ L2(Ω, RN ), ∇(∇ · ~u) = ∆(I−P)~u • NSE: ~ut +P(~u · ∇~u − f~) = ν∆~u +ν(P∆ − ∆P)~u,

~u|Γ = 0

• commutator/gradient: [∆, P]~u = (I − P)(∆~u − ∇∇ · ~u) := ∇pS • For ~u ∈ H 2 ∩ H01(Ω, RN ), Stokes pressure satisfies ∆pS = 0 in Ω with BC: ∂npS = ~n ·(∆ − ∇∇·)~u = −~n ·∇×∇×~u in H −1/2(Γ) first used by Orszag (1986) for consistency in a projection step.

8

Pressure Poisson equation • Unconstrained reformulation of Navier-Stokes equations: ut + ~u · ∇~u + ∇p = ν∆~u + f~,

~u|Γ = 0

Total pressure p is determined by the weak form h∇p, ∇φi = hf~ − ~u·∇~u, ∇φi + νh∆ − ∇∇·)~u, ∇φi

∀φ ∈ H 1(Ω).

• In computation, we use ∆pn = ∇ · (f~n − ~un · ∇~un) in Ω with BC:

~n ·∇pn = ~n · f~n − ν~n · (∇ × ∇ × ~un) on Γ

• ∇pS arises from tangential vorticity at the boundary: Z Z 3D weak form: ∇pS·∇φ = (∇ × ~u)·(~n×∇φ) ∀φ ∈ H 1(Ω) Ω

Γ

9

Space of Stokes pressures Sp = {p ∈ H 1(Ω)/R : ∆p = 0 in Ω, ~n ·∇p ∈ SΓ}, Z SΓ = {f ∈ H −1/2(Γ) : f = 0 ∀components G of Γ}. G

• ∃ a bounded right inverse ∇pS 7→ ~u from Sp → H 2 ∩ H01(Ω, RN ) • In R3, ∇Sp is the space of simultaneous gradients and curls: ∇Sp = ∇H 1(Ω) ∩ ∇ × H 1(Ω, R3)

10

The commutator term is strictly controlled by viscosity If ∇ · ~u = 0, then k[∆, P]~uk = k((I − P)∆~u − ∇∇ · ~uk ≤ k∆~uk. For period box [P, ∆]~u = 0

10

The commutator term is strictly controlled by viscosity If ∇ · ~u = 0, then k[∆, P]~uk = k((I − P)∆~u − ∇∇ · ~uk ≤ k∆~uk. For period box [P, ∆]~u = 0 Main Theorem Let Ω ⊂ RN (N ≥ 2), bounded, ∂Ω ∈ C 3. Then, ∀ε > 0, ∃C ≥ 0, s.t. for all ~u ∈ H 2 ∩ H01(Ω, RN ), Z Z Z 1 2 |(∆P − P∆)~u| ≤ ( + ε) |∆~u|2 + C |∇~u|2 2 Ω Ω Ω

10

The commutator term is strictly controlled by viscosity If ∇ · ~u = 0, then k[∆, P]~uk = k((I − P)∆~u − ∇∇ · ~uk ≤ k∆~uk. For period box [P, ∆]~u = 0 Main Theorem Let Ω ⊂ RN (N ≥ 2), bounded, ∂Ω ∈ C 3. Then, ∀ε > 0, ∃C ≥ 0, s.t. for all ~u ∈ H 2 ∩ H01(Ω, RN ), Z Z Z 1 2 |(∆P − P∆)~u| ≤ ( + ε) |∆~u|2 + C |∇~u|2 2 Ω Ω Ω

Hence our unconstrained NSE is fully dissipative: ~ut + P(~u · ∇~u − f~) + ν[∆, P]~u = ν∆~u,

~u|Γ = 0

11

The commutator term is strictly controlled by viscosity If ∇ · ~u = 0, then k[∆, P]~uk = k((I − P)∆ − ∇∇·)~uk ≤ k∆~uk. For period box [P, ∆]~u = 0. Main Theorem Let Ω ⊂ RN (N ≥ 2), bounded, ∂Ω ∈ C 3. Then, ∀ε > 0, ∃C ≥ 0, s.t. for all ~u ∈ H 2 ∩ H01(Ω, RN ), Z Z Z 1 2 |(∆P − P∆)~u| ≤ ( + ε) |∆~u|2 + C |∇~u|2 2 Ω Ω Ω Hence our unconstrained NSE is fully dissipative: XXX   XX XXX XXX XXX   X  X

XXX   X  ~ X X XXu ~ut + P(~u · ∇~u − f ) + ν[∆, u,  P]~ XX = ν∆~

NSE as perturbed heat equation!

~u|Γ = 0

12

Proof of the theorem – estimate on commutator [∆, P] Decompose ~u ∈ H 2 ∩ H01 into parts parallel and normal to Γ: Let Φ(x) = dist(x, Γ), ~n(x) = −∇Φ(x), ξ a cutoff = 1 near Γ. ~u = ~uk + ~u⊥,

~uk = ξ(I − ~n~nt)~u.

Boundary identities on Γ:

∇ · ~uk = 0,

Stokes pressure satisfies

∆pS = 0,

~n ·∇×∇×~u⊥ = 0. ∂npS = −~n ·∇×∇×~u

Hence ∇pS = (∆P − P∆)~u = (I − P)(∆ − ∇∇·)(~uk + 0). h∇pS, ∇pSi = h∆~uk − ∇∇ · ~uk, ∇pSi = h∆~uk, ∇pSi h∇pS − ∆~uk, ∇pSi = 0 k∆~ukk2 = k∇pSk2 + k∇pS − ∆~ukk2

13

D2N/N2D bounds on tubes Ωs = {x ∈ Ω | Φ(x) < s} Lemma For s0 > 0 small ∃C0 such that whenever ∆p = 0 in Ωs0 and 0 < s < s0 then Z Z 2 t 2 2 |~ n ·∇p| − |(I − ~ n ~ n )∇p| ≤ C s |∇p| 0 Φ 0, x ∈ Ω),

∇ · ~u = h

(t ≥ 0, x ∈ Ω),

~u = ~g

(t ≥ 0, x ∈ Γ),

~u = ~uin

(t = 0, x ∈ Ω).

32

Non-homogeneous side conditions ∂t~u + ~u ·∇~u + ∇p = ν∆~u + f~

(t > 0, x ∈ Ω),

∇ · ~u = h

(t ≥ 0, x ∈ Ω),

~u = ~g

(t ≥ 0, x ∈ Γ),

~u = ~uin

(t = 0, x ∈ Ω).

Unconstrained formulation involves an inhomogeneous pressure pgh: ∂t~u + P(~u · ∇~u − f~ − ν∆~u) + ∇pgh = ν∇∇ · ~u h∇pgh, ∇φi = −h~n · ∂t~g , φiΓ + h∂th, φi + hν∇h, ∇φi

∀φ ∈ H 1(Ω).

∇ · ~u − h satisfies a heat equation with Neumann BCs as before.

33

Summary: NSE is a perturbed heat equation! ~ut + P(~u · ∇~u − f~) = νP∆~u

33

Summary: NSE is a perturbed heat equation! ~ut + P(~u · ∇~u − f~) = νP∆~u + ν∇(∇ · ~u)

33

Summary: NSE is a perturbed heat equation! ~ut + P(~u · ∇~u − f~) = νP∆~u + ν∇(∇ · ~u) = νP∆~u + ν∆(I−P )~u

33

Summary: NSE is a perturbed heat equation! ~ut + P(~u · ∇~u − f~) = νP∆~u + ν∇(∇ · ~u) = νP∆~u + ν∆(I−P )~u = ν∆~u − ν[∆, P]~u

33

Summary: NSE is a perturbed heat equation! ~ut + P(~u · ∇~u − f~) = νP∆~u + ν∇(∇ · ~u) = νP∆~u + ν∆(I−P )~u = ν∆~u − ν[∆, P]~u = ν∆~u − ν∇pS • ∇ · ~u satisfies a heat equation with Neumann BCs • The Stokes pressure term is strictly controlled by viscosity • Stability, existence, uniqueness theory is greatly simplified • Promise of enhanced flexibility in design of numerical schemes

34

Summary: NSE is a perturbed heat equation! XX  XXX ~ut + P(~uX· X∇~ u − f~) = νP∆~u + ν∇(∇ · ~u)

= νP∆~u + ν∆(I−P )~u) = ν∆~u − ν[∆, P]~u   HH H = ν∆~u − νH∇p S H

• ∇ · ~u satisfies a heat equation with Neumann BCs • The Stokes pressure term is strictly controlled by viscosity • Stability, existence, uniqueness theory is greatly simplified • Promise of enhanced flexibility in design of numerical schemes

35

References • G. Grubb and V. A. Solonnikov, Boundary value problems for the nonstationary Navier-Stokes equations treated by pseudodifferential methods, Math. Scand. 69 (1991) 217–290. • J. L. Guermond and J. Shen, A new class of truly consistent splitting schemes for incompressible flows, J. Comp. Phys. 192 (2003) 262–276. • H. Johnston and J.-G. Liu, Accurate, stable and efficient Navier-Stokes solvers based on explicit treatment of the pressure term, J. Comput. Phys. 199 (1) (2004) 221–259 • N. A. Petersson, Stability of pressure boundary conditions for Stokes and Navier-Stokes equations, J. Comp. Phys. 172 (2001) 40-70. • L. J. P. Timmermans, P. D. Minev, F. N. Van De Vosse, An approximate projection scheme for incompressible flow using spectral elements, Int. J. Numer. Methods Fluids 22 (1996) 673–688.

36

Existence and uniqueness theorem Assume ~uin ∈ Huin := H 1(Ω, RN ), f~ ∈ Hf := L2(0, T ; L2(Ω, RN )), ~g ∈ Hg := H 3/4(0, T ; L2(Γ, RN )) ∩ L2(0, T ; H 3/2(Γ, RN )) ∩{~g | ∂t(~n · ~g ) ∈ L2(0, T ; H −1/2(Γ))}, h ∈ Hh := L2(0, T ; H 1(Ω)) ∩ H 1(0, T ; (H 1)0(Ω)). and the compatibility conditions ~g = ~uin

(t = 0, x ∈ Γ),

h~n · ∂t~g , 1iΓ = h∂th, 1iΩ.

Then ∃T ∗ > 0 so that a unique strong solution exists, with ~u ∈ L2([0, T ∗], H 2) ∩ H 1([0, T ∗], L2) ,→ C([0, T ∗], H 1).