In honor of Prof. Hokee Minn
|
|
- Jerome Carter
- 5 years ago
- Views:
Transcription
1 In honor of Prof. Hokee Minn It is my greatest honor to be here today to deliver the nd Minn Hokee memorial lecture. I would like to thank the Prof. Ha for the nomination and the selecting committee for giving me this honor. I
2 On incompressible Navier-Stokes equation (NSE): a new approach for analysis and computation Jian-Guo Liu The title of my talk is: on... Collaborators: Weinan E, Hans Johnston, Jie Liu, Bob Pego, Cheng Wang The nd MINNHOKEE Memorial Lecture Seoul National University May 30, 007
3 Flow phenomena modeled by incompressible NSE Lift, viscous drag, boundary layers & separation, vortex shedding,... The NSE is a simplest mathematical eqn which describes a wide range of phenomena in fluid dynamics, such as lift, viscous drag, vprtex shedding... These are all fundamental phenomena in nature which have important applications in real life.
4 Here I show a movie which simulates fluid passing a cylinder with a Reynold s number of 50,000. Solving such a problem accurately is a very hard problem. We solved this using numerical method we developed. Most recently, we solved the case with the Reynold number of 1000,000. method: vorticity-stream function, 4th order compact FD, RK4, fast FFT Poisson solver, FGP (Johnston, Liu & E)
5 3 Navier-Stokes equations for incompressible flow u t + u u + p = ν u + f in Ω momentum u = 0 in Ω incompressibility u = 0 on Ω no slip Claude-Louis Navier George Gabriel Stokes Pressure is like a Lagrange multiplier to enforce incompressibility. Q: How should one compute the pressure? (stably, accurately, fast)
6 4 Navier-Stokes equations for incompressible flow u t + u u + p = ν u + f in Ω momentum u = 0 in Ω incompressibility u = 0 on Ω no slip Claude-Louis Navier George Gabriel Stokes Pressure is like a Lagrange multiplier to enforce incompressibility. Q: How should one compute the pressure? (stably, accurately, fast) A: Solve p = ( f u u) n p = n ν( ) u + n f in Ω on Ω
7 5 Navier-Stokes equations for incompressible flow u t + (l. l. o. t.) = ν u + f in Ω momentum u = 0 in Ω incompressibility u = 0 on Ω no slip Claude-Louis Navier George Gabriel Stokes Pressure is like a Lagrange multiplier to enforce incompressibility. Q: How should one compute the pressure? (stably, accurately, fast) A: Solve p = ( f u u) n p = n ν( ) u + n f in Ω on Ω Upshot: Well-posed IBVP like a perturbed diffusion equation!
8 6 Traditional implicit Stokes/explicit advection discretization (Temam s book, many mixed FEM): u n+1 u n t + p n+1 ν u n+1 = f n u n u n u n+1 = 0 u n+1 = 0 on Ω New: proof of unconditional stability and convergence for u n+1 u n t ν u n+1 = f n u n u n p n u n+1 = 0 on Ω p n = ( f n u n u n ) n p n = n ν( ) u n + n f n on Ω
9 7 Outline Derive an unconstrained formulation of NSE in terms of the Laplace-Leray commutator P P = p Stokes Why viscosity dominates pressure a new estimate ( Simple theory of existence & uniqueness for strong solutions) New interpretation of existing projection methods and gauge methods (Kim & Moin, Orszag et al, Timmermans, Henshaw- Petersson, Johnston-Liu, E-Liu, Guermond-Minev-Shen) Simple C 0 finite element solvers apparently needing no velocity/pressure compatibility restrictions for stability C 1 finite element solvers with full theory of unconditional stability, convergence, error estimates, in C 3 domains
10 Projection Method, Chorin 68, Temam 69 u u n t + (u n )u n = ν u, u Γ = 0 u n+1 u + p n+1 = 0, u n+1 = 0, u n+1 n Γ = 0 t Efficient, Stable, Numerical boundary layer of order ν t in u, p (Orszag etc 86, E-Liu, 95) nd order projections, Kim, Moin 85, van Kan 86, Orszag etc 86, Bell, Colella, Glaz 89, reduce numerical boundary layer Recent consistent projections: Timmermans etc 96, Gauge method 97, Brown, Cortez, Minion 01, Henshaw, Petersson 01, Guermond, Shen 03, Johnston, Liu 04, Liu, Liu, Pego
11 Impulse Density Formulation Oseledets (1989) rewritten NSE as { t M + (u )M + M ( u) = Re 1 M, u = PM, where u is the velocity field, P is the Helmholtz-Leray projection operator to the space of divergence-free vector fields. Buttke/Chorin 1993: Velicity methods: Lagrangian numerical methods which preserve the Hamiltonian structure of incompressible fluid flow; Magnetization Roberts 197 A Hamiltonain theory for weakly interacting vortices Summers & Chorin 1995, Cortez 1995 etc
12 Impulse density formulation E-Liu, JCP Fourier Collocation Method Finite Difference Method
13 4 Impulse Density Formulation, E-Liu, JCP 1997 Linearized system: { t ˆM + i(ūk1 + vk ) ˆM + ik 1 ( Mû + N ˆv) = 0 t ˆN + i(ūk1 + vk ) ˆN + ik ( Mû + N ˆv) = 0 ( û ˆv ) = ( I k k ) ( ˆM k ˆN Take N = 0, v = 0 and k 1 = 0, we can show that û = ˆM. Hence ( ) ( ) ( ) ˆM 0 0 ˆM t + ik M = 0 ˆN 1 0 ˆN ˆM(k, t) = ˆM(k, 0), ˆN(k, t) = ˆN(k, 0) + ik t M ˆM(k, 0) marginally ill-posed. or, weakly well posed )
14 5 Impulse Density Formulation, E-Liu, JCP 1997 A Simplified Velocity-Impulse Density Formulation/gauge method: { t M + (u )u + q = 1 Re M u = PM If we write M = u + φ, then pressure is given by p = q + t φ 1 Re φ The variable q is arbitrary and can be chosen to suit our numerical purposes. In computation we set q = 0. And we renamed this case as gauge method.
15 Impulse density formulation E-Liu, JCP Re=150, t=/pi, finite difference method Re=10000, t=1000, 18X18 Staggered grid
16 3 Gauge Method, E & Liu 97, 03 Set a = u φ, The momentum equation becomes ( ) a t + (u )u + t φ ν φ + p = ν a. Incompressibility constraint becomes φ = a, impose BC φ n leads to gauge choice = 0, a n = 0, a τ = φ τ, on Ω, t φ ν φ = p,
17 5 Gauge Method, E & Liu 97, 03 Backward Euler, u n+1 = a n+1 + φ n+1 a n+1 a n t + (u n )u n = a n+1 a n+1 n Γ = 0, a n+1 τ Γ = φn τ, φ n+1 = a n+1, φ n+1 n Γ = 0 Theorem (Wang & Liu 00). Unconditional Stability for Stokes Flow u n+1 u n t + û n + ( φ n+1 φ n ) 0.
18 6 Geometric gauge 6 6 EL gauge J. Fluid Mech. (1999), vol. 391, pp Printed in the United Kingdom c 1999 Cambridge University Press JFM, 1999 Impulse formulation of the Euler equations: general properties and numerical methods t= By G I O V A N N I R U S S O 1 AND P E T E R S M E R E K A 0 G. Russo and P. Smereka (a) t= Geometric gauge (b) t= EL gauge vorticity pressure pressure 4 6 Figure 4. PhiD results. The left-hand column shows contours of ω at times t = 15, 1, and 7 from top to bottom, the contours are separated by 0.1 The middle column shows contours of p at times t = 15, 1, and 7 from top to bottom, in the geometric gauge; the contours are separated by 0.5. The right-hand column shows the contours of p in the EL gauge; the contours are separated by Figure 5. The L, L, and L 1 norms are shown for the results in figure 4, corresponding to the geometric gauge (a) and EL gauge (b). In each plot, the upper curve is the L norm, the middle one is the L norm, and the lower curve is the L 1 norm. The norms have been normalized to be equal to 1 at t = 0.
19 4 Gauge Method, E & Liu 97, 03 Gauge formulation of NSE with u = a + φ a t + u u = ν a, a n Γ = 0, a τ Γ = τ φ, φ φ = a, n Γ gauge choice becomes: t φ ν φ = p, = 0 a t + P(u )u = ν a, a n Γ = 0, a τ Γ = τ φ, φ φ = a, = 0 n Γ gauge choice becomes: t φ ν φ = p s. Equivalent with the Pressure Poisson formulation.
20 Pressure Poisson formulation, Johnston & Liu 04 Navier-Stokes equation (NSE): u t + (u )u + p = ν u, u Γ = 0, Incompressibility u = 0 is replaced by p = (u u), n p = νn ( u) Γ, BC first used by Orszag (1986), motivated by identity u = u + ( u). t ( u) = ν ( u), ( u) t=0 = 0, n ( u) Γ = 0.
21 Pressure Poisson formulation, Johnston & Liu 04 Navier-Stokes equation (NSE): u t + (u )u + p = ν u, u Γ = 0, Incompressibility u = 0 is replaced by p = (u u), n p = νn ( u) Γ, BC first used by Orszag (1986), motivated by identity u = u + ( u). t ( u) = ν ( u), ( u) t=0 = 0, n ( u) Γ = 0. Recast as, u t + P( u u) = νp u + ν ( u), u Γ = 0
22 5 Pressure Poisson formulation, Johnston & Liu 04 Galerkin formulation (good for C 0 finite element method): u t, v + u u, v + p, v = ν u, v, v H 1 0 p, ψ = u u, ψ + ν u, n ψ Γ, ψ H 1 Finite element spaces: X h H0(Ω, 1 R N ), v h X h and q h Y h, require Y h H 1 (Ω)/R. For all 1 t ( u n+1 h, v h u n h, v h ) + ν u n+1 h, v h = p n h, v h + f n u n h u n h, v h. p n h, q h = f n u n h u n h, q h + ν u n h, n q h Γ,
23 !0.5 Driven cavity flow!4! !3 1 3!5! !1!0.5!0.5! 0!1!1!!4!!3!5!3!3!5!4!3! !!!0.5!0.5!1! ! !!1!! !!1! !1! Re= P4 elements. CN/AB time stepping
24 Backward facing step flow 0.5 0! Re= P4 elements. CN/AB time stepping 0.5 X 1 X 0 S! X 3 Re= P4 elements. CN/AB time stepping
25 Hermann von Helmholtz W.V.D. Hodge Hermann von Helmholtz Hermann von Helmholtz Helmholtz projection onto divergence-free vector fields 5 L (Ω, R N ) = PL (Ω, R N ) H 1 (Ω) Given v L (Ω, R N ), there exists q H 1 (Ω) so that v = P v q satisfies P v, φ = v + q, φ = 0 for all φ H 1 (Ω).
26 5 Helmholtz projection onto divergence-free vector fields L (Ω, R N ) = PL (Ω, R N ) H 1 (Ω) Given v L (Ω, R N ), there exists q H 1 (Ω) so that v = P v q satisfies P v, φ = v + q, φ = 0 for all φ H 1 (Ω). Then (P v) = 0 in Ω, n P v = 0 on Γ. Note P v H(div; Ω) = { f L (Ω, R N ) : f L }, so n P v H 1/ (Γ) by a standard trace theorem.
27 Jean Leray Olga A. Ladyzhenskaya Traditional unconstrained formulation of NSE mgrestosio Kato ( ) Formally t ( u) = 0 u t + P( u u f) = νp u, 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ška-brezzi)
28 10 Unconstrained formulation of NSE u t + P( u u f) = νp u+ν u No difference! as long as u = 0. Formally, w = u satisfies a heat equation with Neumann BCs: t w = ν w in Ω, n w = 0 on Ω. w = n u = 0 = exponential decay Ω Ω
29 10 Unconstrained formulation of NSE u t + P( u u f) = νp u+ν u No difference! as long as u = 0. Formally, w = u satisfies a heat equation with Neumann BCs: t w = ν w in Ω, n w = 0 on Ω. w = n u = 0 = exponential decay Ω Ω Equivalent to a reduced formulation of Grubb & Solonnikov (1991) studied as a nondegenerate parabolic pseudodifferential IBVP
30 11 Equivalent formulation with commutator Lemma For all u L (Ω, R N ), u = u P u
31 11 Equivalent formulation with commutator Lemma For all u L (Ω, R N ), u = u P u U-NSE: u t + P( u u f) = ν u ν( P P ) u
32 11 Equivalent formulation with commutator Lemma For all u L (Ω, R N ), u = u P u U-NSE: u t + P( u u f) = ν u ν( P P ) u P u = ( ) u ( P P ) u = (I P)( ) u =: p S Stokes pressure BVP!! For u H H0(Ω, 1 R N ), we get p S H(div; Ω), with p S = 0 in Ω, n p S = n ( ) u in H 1/ ( Ω). Stokes pressure arises from tangential vorticity at the boundary: 3D weak form: p S φ = ( u) ( n φ) φ H 1 (Ω) Ω Ω
33 1 Unconstrained form pressure decomposition U-NSE: u t + P( u u f) = ν u ν( P P ) u u = 0 on Ω Pressure-Poisson form: u t + u u + p = ν u + f u = 0 on Ω p = ( f u u) in Ω n p = n ν( ) u + n f on Ω The pressure p consists of two parts: p E is nonlinear but lower order p S is nd order but a commutator p = p E + ν p S
34 13 Viscosity dominates Stokes pressure gradient Theorem Let Ω R N (N ) be a bounded domain with C 3 boundary. Then, for every ε > 0 there exists C 0 such that for all u H H 1 0(Ω, R N ), Ω ( P P ) u ( ) 1 + ε Ω u + C Ω u,
35 13 Viscosity dominates Stokes pressure gradient Theorem Let Ω R N (N ) be a bounded domain with C 3 boundary. Then, for every ε > 0 there exists C 0 such that for all u H H 1 0(Ω, R N ), Ω ( P P ) u ( ) 1 + ε Ω Ingredients of the (elementary!) proof: u + C Ω u, Decompose u H H 1 0 into parts parallel and normal to Ω: Let Φ(x) = dist(x, Ω), n(x) = Φ(x), ξ a cutoff = 1 near Ω. u = u + u, u = ξ(i n n t ) u.
36 14 1. Boundary identities (i) u = 0 and (ii) u = 0 on Ω Hence p S = ( P P ) u = (I P)( )( u + 0) u = u + u, u + u (1 ε) u C u goal: ( ε) p S C u
37 15. Orthogonal decomposition of u u = p S + p S u p S = (I P)( ) u = u u + P v p S u H0(Ω) 1 + PL, since u = 0 on Ω p S PL, and p S = 0 p S, p S u = p S, u = p S, u = 0 It remains only to show p S u (1 ε) p S C u Break in parallel, perpendicular parts, and use somewhat novel...
38 16 3. Neumann-to-Dirichlet estimates on tubes Ω s = {x Ω Φ(x) < s} Lemma For s 0 > 0 small C 0 such that whenever p = 0 in Ω s0 and 0 < s < s 0 then p p C 0 s p Φ<s Φ<s 0 Morally: p Φ<s p Φ<s + junk In a half space: p = p exactly ( nˆp = k ˆp) Now: p S u s n ( p S u ) s (1 ε) p s C u ( 1 ε) p S s C u
39 17 Stability of time-difference scheme with pressure explicit (related: Ti96,Pe01,GS03,JL04) u n+1 u n t ν u n+1 = f n u n u n p n, u n+1 = 0 on Ω, p n ν u n + ν u n + u n u n f n, φ = 0 φ H 1 (Ω). p n ν p n S + u n u n f n, φ = 0 φ H 1 (Ω). Taking φ = p n gives the pressure estimate p n ν p n S + f n u n u n
40 Unconditional stability theorem for N =, 3 Liu, Liu & Pego, CPAM Theorem Take f L (0, T ; L (Ω, R N )), u 0 H H 1 0(Ω, R N ). Then positive constants T and C depending only upon Ω, ν and M 0 := u 0 + ν t u 0 + so that whenever n t T we have T 0 f, sup u k + 0 k n ( n 1 k=0 n u k t C, k=0 u k+1 u k t + u k u k ) t C.
41 Unconditional stability theorem for N =, 3 Liu, Liu & Pego, CPAM Theorem Take f L (0, T ; L (Ω, R N )), u 0 H H 1 0(Ω, R N ). Then positive constants T and C depending only upon Ω, ν and M 0 := u 0 + ν t u 0 + so that whenever n t T we have T 0 f, sup u k + 0 k n ( n 1 k=0 n u k t C, k=0 u k+1 u k t + u k u k ) t C. u in H 1 0 strong solution u L (0, T ; H ) H 1 (0, T ; L )
42 3 C 0 finite element schemes Space-continuous projection method after Orszag et al (1986): u = P u n = u n q n p n ν u n + ν u n + u u f n, φ = 0 u n+1 u t φ H 1 (Ω). ν u n+1 = f n u u p n, u n+1 Ω = 0. C 0 finite element method: v h X h H 1 0(Ω) N, φ h Y h H 1 (Ω), u n+1 h u h, v h t u n h q h, φ h = 0, u h := u n h q h p n h, φ h = f n u h u h, φ h + ν u n h, n φ h Ω, + ν u n+1 h, v h = p n h, v h + f n u h u h, v h.
43 4 Analogous C 1 finite element schemes Finite element spaces: X h H H 1 0(Ω, R N ), Y h H 1 (Ω)/R. For all v h X h and φ h Y h, require u h, v h = P u n, v h = u n h, v h u n h, v h p n h, φ h = f n u h u h, φ h + ν u n h, n φ h Ω u n+1 h u h, v h t + ν u n+1 h, v h = p n h + u h u h f n, v h
44 4 Analogous C 1 finite element schemes Finite element spaces: X h H H 1 0(Ω, R N ), Y h H 1 (Ω)/R. For all v h X h and φ h Y h, require u h, v h = P u n, v h = u n h, v h u n h, v h p n h, φ h = f n u h u h, φ h + ν u n h, n φ h Ω, u n+1 h u h, v h t Pressure estimate: Take φ h = p n h, then + ν u n+1 h, v h = p n h + u h u h f n, v h p n h f n u h u h + ν p S ( u n h) Unconditional stability and convergence follow as before!
45 5 No compatibility required between velocity and pressure spaces In the proof, we do not require the spaces X 0,h and Y h satisfy the inf-sup stability condition for weak solutions: ( vh )q h inf sup q h Y h v h X 0,h v h q h β > 0 The inf-sup condition is necessary and sufficient for stability for weak solutions to the mixed formulation of steady Stokes problems It seriously constrains kinds of discretizations that are weakly stable Our proofs are for strong solutions in C 3 domains and we compute pressure stably from a well-posed BVP. In standard test problems, the simple C 0 scheme works well. Nevertheless, the situation is not fully understood.
46 8 Stability check: C 0 FE scheme with smooth solution 10 0 u!u L p!p L u!u L p!p L t v.s. t 10 1 (u!u) L (p!p) L u x +v y L ω!ω L t v.s. t !1 10!1 10! t 10! t 1: BE/FE, T=1000, fixed mesh, t 8
47 9 Spatial accuracy check for C 0 FE scheme !1 10!1 10! 10!3 u!u L p!p L u!u L 10! (u!u) L (p!p) L u +v x y L p!p L ω!ω L 10!4 h v.s. h 10! 10!1 h 10!3 h v.s. h 10! 10!1 h : P1/P1 finite elements
48 30 Spatial accuracy check for C 0 FE scheme ! 10! 10!4 10!6 u!u L 10!4 (u!u) L (p!p) L p!p L u x +v y L 10!8 u!u L p!p L h v.s. h 5 10!6 ω!ω L h v.s. h 4 h v.s. h 5 10!10 10!1 h 10!8 10!1 h 3: P4/P4 finite elements
49 31 Stability check: C 1 FE scheme with smooth solution !1 u!u L p!p L u!u L p!p L t v.s. t (u!u) L (p!p) L (u!u) L (p!p) L u x +v y L t v.s. t 10! 10!1 10!3 10! 10! t 10! 10! 10! t 4: BE/FE, T=1000, fixed mesh, t 8
50 3 Spatial accuracy check: C 1 FE scheme with smooth solution 10! 10!1 10!3 10! 10!4 10!3 10!5 10!4 10!6 10!7 10!8 u!u L p!p L u!u L p!p L h v.s. h 4 10!5 10!6 10!7 (u!u) L (p!p) L (u!u) L (p!p) L u x +v y L h v.s. h 3 h v.s. h !9 10! 10!1 h 10!8 10! 10!1 h
51 33 Numerics: driven cavity and backward facing step (u,v)=(1,0) (u,v) given (u,v) given
52 C 0 FEM. Driven Cavity. Re=1000. P4/P elements. 34
53 35 ʯ ۿ ޱ ԏݏ ȶԏ ੴԏ ӹۿſ ϻ ϻ ÿ ۿ ۿ ۿ
54 36 C 0 FEM. Backward facing step. Re=600. P4/P elements 0.5 X 1 X 0 S! X 3
55 C 1 FEM. Driven Cavity. Re= C 1 elements 37
56 38 C 1 FEM. Backward facing step. Re= ! (with a few recycled C 1 elements near the step)
57 39 C 1 FEM. Backward facing step. Re= C 1 elements X 1 X S! X 3 (with a few recycled C 1 elements near the step)
58 33 Summary: NSE is a perturbed heat equation! u t + P( u u f) = νp u
59 33 Summary: NSE is a perturbed heat equation! u t + P( u u f) = νp u + ν ( u)
60 33 Summary: NSE is a perturbed heat equation! u t + P( u u f) = νp u + ν ( u) = νp u + ν (I P ) u
61 33 Summary: NSE is a perturbed heat equation! u t + P( u u f) = νp u + ν ( u) = νp u + ν (I P ) u = ν u ν[, P] u
62 33 Summary: NSE is a perturbed heat equation! u t + P( u u f) = νp u + ν ( u) = νp u + ν (I P ) u = ν u ν[, P] u = ν u ν p S 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
63 34 Summary: NSE is a perturbed heat equation! u t + P( u u f) = νp u + ν ( u) = νp u + ν (I P ) u) = ν u ν[, P] u = ν u ν p S 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
64 Thank You!
65 3 Non-homogeneous side conditions t u + u u + p = ν u + f u = h u = g u = u in (t > 0, x Ω), (t 0, x Ω), (t 0, x Γ), (t = 0, x Ω).
66 3 Non-homogeneous side conditions t u + u u + p = ν u + f u = h u = g u = u in (t > 0, x Ω), (t 0, x Ω), (t 0, x Γ), (t = 0, x Ω). Unconstrained formulation involves an inhomogeneous pressure p gh : t u + P( u u f ν u) + p gh = ν u p gh, φ = n t g, φ Γ + t h, φ + ν h, φ φ H 1 (Ω). u h satisfies a heat equation with Neumann BCs as before.
67 36 Existence and uniqueness theorem Assume u in H uin := H 1 (Ω, R N ), f H f := L (0, T ; L (Ω, R N )), g H g := H 3/4 (0, T ; L (Γ, R N )) L (0, T ; H 3/ (Γ, R N )) { g t ( n g) L (0, T ; H 1/ (Γ))}, h H h := L (0, T ; H 1 (Ω)) H 1 (0, T ; (H 1 ) (Ω)). and the compatibility conditions g = u in (t = 0, x Γ), n t g, 1 Γ = t h, 1 Ω. Then T > 0 so that a unique strong solution exists, with u L ([0, T ], H ) H 1 ([0, T ], L ) C([0, T ], H 1 ).
68 6 Kim Moin s Method vs Gauge Method Kim-Moin: u u n t u n = 0, + (u u) n+1 = ν u +u n, u τ = t pn τ, p n+1 p n+1 = u, n Γ = 0 u n+1 = u + t p n+1 Gauge Method: u n+1 = a n+1 + φ n+1 a n+1 a n + (u u) n+1/ = ν an+1 + a n t a n+1 n Γ = 0, a n+1 τ Γ = (φn p n 1 ) τ φ n+1 = a n+1, φ n+1 n Γ = 0,
69 3 Non-homogeneous side conditions t u + u u + p = ν u + f u = h u = g u = u in (t > 0, x Ω), (t 0, x Ω), (t 0, x Γ), (t = 0, x Ω).
70 3 Non-homogeneous side conditions t u + u u + p = ν u + f u = h u = g u = u in (t > 0, x Ω), (t 0, x Ω), (t 0, x Γ), (t = 0, x Ω). Unconstrained formulation involves an inhomogeneous pressure p gh : t u + P( u u f ν u) + p gh = ν u p gh, φ = n t g, φ Γ + t h, φ + ν h, φ φ H 1 (Ω). u h satisfies a heat equation with Neumann BCs as before.
Analysis and Computation of Navier-Stokes equation in bounded domains
1 Analysis and Computation of Navier-Stokes equation in bounded domains Jian-Guo Liu Department of Physics and Department of Mathematics Duke University with Jie Liu (National University of Singapore)
More informationGauge finite element method for incompressible flows
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 34: 701 710 Gauge finite element method for incompressible flows Weinan E a, *,1 and Jian-Guo Liu b,2 a Courant Institute
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationDivorcing pressure from viscosity in incompressible Navier-Stokes dynamics
Divorcing pressure from viscosity in incompressible Navier-Stokes dynamics Jian-Guo Liu 1, Jie Liu 2, and Robert L. Pego 3 February 18, 2005 Abstract We show that in bounded domains with no-slip boundary
More informationImplementation of 3D Incompressible N-S Equations. Mikhail Sekachev
Implementation of 3D Incompressible N-S Equations Mikhail Sekachev Navier-Stokes Equations The motion of a viscous incompressible fluid is governed by the Navier-Stokes equations u + t u = ( u ) 0 Quick
More informationCONVERGENCE OF GAUGE METHOD FOR INCOMPRESSIBLE FLOW CHENG WANG AND JIAN-GUO LIU
MATHEMATICS OF COMPUTATION Volume 69, Number 232, Pages 135{1407 S 0025-571(00)0124-5 Article electronically published on March 24, 2000 CONVERGENCE OF GAUGE METHOD FOR INCOMPRESSIBLE FLOW CHENG WANG AND
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationSome recent results for two extended Navier-Stokes systems
Some recent results for two extended Navier-Stokes systems Jim Kelliher UC Riverside Joint work with Mihaela Ignatova (UCR), Gautam Iyer (Carnegie Mellon), Bob Pego (Carnegie Mellon), and Arghir Dani Zarnescu
More informationFinite Difference Schemes for Incompressible Flows in the Velocity Impulse Density Formulation
JOURNAL OF COMPUTATIONAL PHYSICS 130, 67 76 (1997) ARTICLE NO. CP965537 Finite Difference Schemes for Incompressible Flows in the Velocity Impulse Density Formulation Weinan E 1 Courant Institute of Mathematical
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationFourth Order Convergence of Compact Finite Difference Solver for 2D Incompressible Flow
Fourth Order Convergence of Compact Finite Difference Solver for D Incompressible Flow Cheng Wang 1 Institute for Scientific Computing and Applied Mathematics and Department of Mathematics, Indiana University,
More informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
More informationPROJECTION METHOD III: SPATIAL DISCRETIZATION ON THE STAGGERED GRID
MATHEMATICS OF COMPUTATION Volume 71, Number 237, Pages 27 47 S 0025-5718(01)01313-8 Article electronically published on May 14, 2001 PROJECTION METHOD III: SPATIAL DISCRETIZATION ON THE STAGGERED GRID
More informationProjection Method with Spatial Discretizations
Projection Method with Spatial Discretizations Weinan E School of Mathematics Institute for Advanced Study Princeton, NJ 08540 and Jian-Guo Liu Courant Institute of Mathematical Sciences New York, NY 10012
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationAccurate Projection Methods for the Incompressible Navier-Stokes Equations
Accurate Projection Methods for the Incompressible Navier-Stokes Equations David L. Brown 1 Center for Applied Scientific Computing Lawrence Livermore National Laboratory, Livermore, CA 94551, USA E-mail:
More informationWorkshop on PDEs in Fluid Dynamics. Department of Mathematics, University of Pittsburgh. November 3-5, Program
Workshop on PDEs in Fluid Dynamics Department of Mathematics, University of Pittsburgh November 3-5, 2017 Program All talks are in Thackerary Hall 704 in the Department of Mathematics, Pittsburgh, PA 15260.
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationMinimal stabilization techniques for incompressible flows
Minimal stabilization tecniques for incompressible flows G. Lube 1, L. Röe 1 and T. Knopp 2 1 Numerical and Applied Matematics Georg-August-University of Göttingen D-37083 Göttingen, Germany 2 German Aerospace
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationNOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7 Number Pages 3 39 c Institute for Scientific Computing and Information NOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL
More informationLeast-Squares Finite Element Methods
Pavel В. Bochev Max D. Gunzburger Least-Squares Finite Element Methods Spri ringer Contents Part I Survey of Variational Principles and Associated Finite Element Methods 1 Classical Variational Methods
More informationVector Penalty-Projection Methods for the Solution of Unsteady Incompressible Flows
Author manuscript, published in "5th International Symposium on Finite Volumes for Complex Applications, Aussois : France 2008)" Vector Penalty-Projection Methods for the Solution of Unsteady Incompressible
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationERROR ESTIMATES FOR SEMI-DISCRETE GAUGE METHODS FOR THE NAVIER-STOKES EQUATIONS : FIRST-ORDER SCHEMES
MATHEMATICS OF COMPUTATION Volume, Number, Pages S 5-578XX- ERROR ESTIMATES FOR SEMI-DISCRETE GAUGE METHODS FOR THE NAVIER-STOKES EQUATIONS : FIRST-ORDER SCHEMES RICARDO H. NOCHETTO AND JAE-HONG PYO Abstract.
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationHypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th
Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th Department of Mathematics, University of Wisconsin Madison Venue: van Vleck Hall 911 Monday,
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationThe Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations
The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech PDE Conference in honor of Blake Temple, University of Michigan
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationBorel Summability in PDE initial value problems
Borel Summability in PDE initial value problems Saleh Tanveer (Ohio State University) Collaborator Ovidiu Costin & Guo Luo Research supported in part by Institute for Math Sciences (IC), EPSRC & NSF. Main
More informationFinite Difference Methods for 3D Viscous Incompressible Flows in the Vorticity Vector Potential Formulation on Nonstaggered Grids
JOURNAL OF COMPUTATIONAL PHYSICS 138, 57 82 (1997) ARTICLE NO. CP975815 Finite Difference Methods for 3D Viscous Incompressible Flows in the Vorticity Vector Potential Formulation on Nonstaggered Grids
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationHamiltonian aspects of fluid dynamics
Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW AT Re τ = 5200 DATA SET Data provenance: M. Lee 1 & R. D. Moser 1 Database ingest and Web Services: Z. Wu 2, G. Lemson 2, R. Burns 2, A. Szalay
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationBlow-up or No Blow-up? Fluid Dynamic Perspective of the Clay Millennium Problem
Blow-up or No Blow-up? Fluid Dynamic Perspective of the Clay Millennium Problem Thomas Y. Hou Applied and Comput. Mathematics, Caltech Research was supported by NSF Theodore Y. Wu Lecture in Aerospace,
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationREMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID
REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID DRAGOŞ IFTIMIE AND JAMES P. KELLIHER Abstract. In [Math. Ann. 336 (2006), 449-489] the authors consider the two dimensional
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationHIGH-ORDER NUMERICAL METHODS FOR PRESSURE POISSON NAVIER-STOKES EQUATIONS
HIGH-ORDER NUMERICAL METHODS FOR PRESSURE POISSON EQUATION REFORMULATIONS OF THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS DOCTOR OF PHILOSOPHY Dong Zhou August, 214 Benjamin Seibold, Advisory Chair, Mathematics
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
More informationSecond-order finite-element projection method for 3D flows
Second-order finite-element projection method for 3D flows J.-L. Guermond I and L. Quartapelle 2 1 LIMSI-CNRS, BP 133, 91403 Orsay, France, (guermond@limsi.fr) 2 Dipartimento di Fisica del Politecnico
More informationFast Solvers for Unsteady Incompressible Flow
ICFD25 Workshop p. 1/39 Fast Solvers for Unsteady Incompressible Flow David Silvester University of Manchester http://www.maths.manchester.ac.uk/~djs/ ICFD25 Workshop p. 2/39 Outline Semi-Discrete Navier-Stokes
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationRunge-Kutta-Chebyshev Projection Method
Runge-Kutta-Chebyshev Projection Method Zheming Zheng Linda Petzold Department of Mechanical Engineering, University of California Santa Barbara, Santa Barbara, CA 93106, USA June 8, 2006 This work was
More informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationOf The Navier Stokes Equations Nar Associates
We have made it easy for you to find a PDF Ebooks without any digging. And by having access to our ebooks online or by storing it on your computer, you have convenient answers with of the navier stokes
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationAn evaluation of a conservative fourth order DNS code in turbulent channel flow
Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationA Review of Solution Techniques for Unsteady Incompressible Flow
Zeist 2009 p. 1/57 A Review of Solution Techniques for Unsteady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/57 PDEs Review : 1966 1999 Update : 2000
More informationEnergy dissipation caused by boundary layer instability at vanishing viscosity
Energy dissipation caused by boundary layer instability at vanishing viscosity Marie Farge, Ecole Normale Supérieure, Paris Kai Schneider, Université d Aix-Marseille in collaboration with Romain Nguyen-Nouch
More informationVector penalty-projection methods for outflow boundary conditions with optimal second-order accuracy
Vector penalty-projection methods for outflow boundary conditions with optimal second-order accuracy Philippe Angot, Rima Cheaytou To cite this version: Philippe Angot, Rima Cheaytou. Vector penalty-projection
More informationBlow-up or No Blow-up? the Role of Convection in 3D Incompressible Navier-Stokes Equations
Blow-up or No Blow-up? the Role of Convection in 3D Incompressible Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech Joint work with Zhen Lei; Congming Li, Ruo Li, and Guo
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationThe enigma of the equations of fluid motion: A survey of existence and regularity results
The enigma of the equations of fluid motion: A survey of existence and regularity results RTG summer school: Analysis, PDEs and Mathematical Physics The University of Texas at Austin Lecture 1 1 The review
More informationA Projection FEM for Variable Density Incompressible Flows
Journal of Computational Physics 65, 67 88 (000) doi:0.006/jcph.000.6609, available online at http://www.idealibrary.com on A Projection FEM for Variable Density Incompressible Flows J.-L. Guermond and
More informationOn a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions
Proceedings of the 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August -, 5 pp36-41 On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability
More informationDivergence-free wavelet projection method for incompressible viscous flow on the square
Divergence-free wavelet proection method for incompressible viscous flow on the square Souleymane Kadri Harouna, Valérie Perrier To cite this version: Souleymane Kadri Harouna, Valérie Perrier. Divergence-free
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationIMA Preprint Series # 2143
A BASIC INEQUALITY FOR THE STOKES OPERATOR RELATED TO THE NAVIER BOUNDARY CONDITION By Luan Thach Hoang IMA Preprint Series # 2143 ( November 2006 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationThe Navier-Stokes Equations with Time Delay. Werner Varnhorn. Faculty of Mathematics University of Kassel, Germany
The Navier-Stokes Equations with Time Delay Werner Varnhorn Faculty of Mathematics University of Kassel, Germany AMS: 35 (A 35, D 5, K 55, Q 1), 65 M 1, 76 D 5 Abstract In the present paper we use a time
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW DATA SET Data provenance: J. Graham 1, M. Lee 2, N. Malaya 2, R.D. Moser 2, G. Eyink 1 & C. Meneveau 1 Database ingest and Web Services: K. Kanov
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationWeierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, Berlin, Germany,
Volker John On the numerical simulation of population balance systems Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, 10117 Berlin, Germany, Free University of Berlin, Department
More informationIssues for a mathematical definition of LES
Issues for a mathematical definition of LES Jean-Luc Guermond 1 and Serge Prudhomme 2 1 Texas A&M University, College Station TX 77843, USA, and LIMSI, CNRS UPR 3251, BP 133 Orsay Cedex, France, guermond@math.tamu.edu
More informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationCONNECTIONS BETWEEN A CONJECTURE OF SCHIFFER S AND INCOMPRESSIBLE FLUID MECHANICS
CONNECTIONS BETWEEN A CONJECTURE OF SCHIFFER S AND INCOMPRESSIBLE FLUID MECHANICS JAMES P. KELLIHER Abstract. We demonstrate connections that exists between a conjecture of Schiffer s (which is equivalent
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationThe Kolmogorov Law of turbulence
What can rigorously be proved? IRMAR, UMR CNRS 6625. Labex CHL. University of RENNES 1, FRANCE Introduction Aim: Mathematical framework for the Kolomogorov laws. Table of contents 1 Incompressible Navier-Stokes
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More information