Analysis and Computation of Navier-Stokes equation in bounded domains
|
|
- Georgia Andrews
- 5 years ago
- Views:
Transcription
1 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) Bob Pego (Carnegie Mellon)
2 2 Outline: The Laplace-Leray commutator and Stokes pressure A formula/decomposition for Navier-Stokes pressure Extended/unconstrained NSE Estimate for the Laplace-Leray commutator Analysis and numerical analysis for (extended) NSE Finite element method with equal order polynomials Stable third order time-split/projection methods Numerics: stability accuracy check, driven cavity, backward-facing step, flow past cylinder
3 3 Navier-Stokes equations for incompressible flow u t + u u + p = ν u + f in Ω momentum u = in Ω incompressibility u = on Ω = Γ no slip Ways to regard and treat pressure: (i) Pressure is like a Lagrange multiplier to enforce incompressibility (ii) Pressure can be eliminated by projection on divergence-free fields or taking curl to get a vorticity equation (iii) Pressure can be found from u and f by solving Poisson equations
4 4 Leray-Helmholtz projection P onto divergence-free fields L 2 (Ω, R N ) = H 1 (Ω) PL 2 (Ω, R N ) v = φ + w, w, ψ = v φ, ψ = for all ψ H 1 (Ω). φ = v, n φ = n v Γ w = in Ω, n w = on Ω. Notation: w = Pv, φ = Qv.
5 5 The Laplace-Leray commutator and Stokes pressure For all v H 2 (Ω, R N ) we have (I P)v = φ = φ = v, Pv = ( )v = v, P v = P( )v ( P P )v = (I P)( )v. We define the Stokes pressure for u H 2 (Ω, R N ) via p S (u) = Q( )u, then p S (u) = ( P P )u. Note ( )u H(div; Ω), hence p S satsifies the BVP p S = in Ω, n p S = n ( )u in H 1/2 (Γ)
6 6 p S arises from tangential vorticity at the boundary Johnston-Liu p S, q = n ( u), q Γ, q H 1 (Ω)
7 Johnston-Liu 7
8 8 Space of Stokes pressures p S = in Ω, n p S = n ( )u = f in H 1/2 (Γ) S p = {p H 1 (Ω)/R : p = in Ω, n p S Γ }, S Γ = {f H 1/2 (Γ) : f = components G of Γ}. G a bounded right inverse p S u from S p H 2 H(Ω, 1 R N ) In R 3, S p is the space of simultaneous gradients and curls: S p = H 1 (Ω) H 1 (Ω, R 3 )
9 9 A formula/decomposition for the Navier-Stokes pressure Suppose that u = Pu is a strong solution of NSE: u t + p = ν u + f u u u = u = in Ω in Ω on Ω Apply P and note P p =, P u = u p S (u): u t + ν p S (u) = ν u + P(f u u) Subtracting, we find that (up to spatial constant) necessarily p = νp S (u) + Q(f u u). p = (f u u), n p = n f n u Γ
10 1 Pressure formula with inflow/outflow Suppose we require: u = g on Ω, where to conserve volume, n g = (t ). Ω Let R(g) solve R(g) = in Ω, n R(g) = n g on Ω. If u = then u R(g) = P(u R(g)) = Pu. Applying P to NSE we now find t (u R(g)) + ν p S (u) = ν u + P(f u u), p = νp S (u) R( t g) + Q(f u u). For q H 1 (Ω) such that: p, q = ν u, n q Γ n g t, q Γ + f u u, q
11 11 Extended/unconstrained NSE u t + u u + p = ν u + f, u Ω = g p = νp S (u) R( t g) + Q(f u u) p = (f u u), n p = n (f u u g t ) n u Γ p, q = ν u, n q Γ n g t, q Γ + f u u, q p, q = ν ( )u, q n g t, q Γ + f u u, q n g t, q Γ + u u + p, q = ν u + f, q ν u, q ( u) t, q + ν ( u), q = ( u) t = ν ( u) in Ω, n ( u) Γ = If ( u) =, then u. Equivalent to the original NSE t=
12 12 Comparing with traditional formulation of NSE u t + P(u u f) = νp u Formally t ( u) = Perform analysis in spaces of divergence-free fields (unconstrained Stokes operator P u is incompletely dissipative). Inf-Sup/LBB condition(ladyzhenskaya-babuška-brezzi) u t + P(u u f) = ν u ν p S = νp u + ν ( u) Improved enforcement of the divergence constraint Formally t ( u) = ν ( u), n ( u) = Equivalent to a reduced formulation of Grubb & Solonnikov (1991) studied as a pseudodifferential IBVP
13 13 Local-time well posedness for extended NSE u t + u u + p = ν u + f, u Ω = p = νp S (u) R( t g) + Q(f u u) Assume u in H uin := H 1 (Ω, R N ), f H f := L 2 (, T ; L 2 (Ω, R N )), g H g := H 3/4 (, T ; L 2 (Γ, R N )) L 2 (, T ; H 3/2 (Γ, R N )) {g t (n g) L 2 (, T ; H 1/2 (Γ))}, g = u in (t =, x Γ), n t g, 1 Γ = Then T > so that a unique strong solution exists, with u L 2 ([, T ], H 2 ) H 1 ([, T ], L 2 ) C([, T ], H 1 ).
14 14 Estimate for the Laplace-Leray commutator Recall [, P]u = ( P P )u = p S For period box [P, ]u = If u =, then [, P]u = ((I P) )u u. Theorem Let Ω R N (N 2) be a bounded C 3 domain Then ε > C such that for all u H 2 H(Ω, 1 R N ), ( ) 1 ( P P )u ε u 2 + C u 2 Ω Hence our unconstrained NSE is fully dissipative: u t + P(u u f) + ν p S = ν u NSE as perturbed heat equation! Ω Ω
15 15 An example in half period strip < x < 2π, y u(x, y) = sin(kx)ye ky, any v(x, y) with v(x, ) = Stokes pressure: p(x, y) = cos(kx)e ky p =, y p(x, ) = k cos(kx) = x ( y u x v) y= u = 2k sin(kx)e ky = 2 x p, (u, ) =, (, v) = at y = (boundary properties ) p x 2 = p y 2 = πk 2 (equal partition) u x p y p, x p = p x 2 p y 2 = (orthogonality) y p This implies: p 2 = 1 2 u ( u 2 + v 2 )
16 16 Basic ingredients in the (elementary) proof: Let Φ(x) = dist(x, Γ), n(x) = Φ(x) Ω s = {x Ω Φ(x) < s}, f, g s := f, g Ωs Choose a cutoff function ξ so ξ(x) = 1 (Φ < 1 2s), ξ(x) = (Φ > s). Decompose velocity field into parts parallel and normal to Γ: u = u + u, u = ξnn t u + (1 ξ)u, u = ξ(i nn t )u. (i) A parallel-normal decomposition u = u + u satisfying on Ω: u =, u = (ii) An orthogonality identity: u 2 = p S 2 + p S u 2 (iii) Sharp N2D estimates on tubes: p S 2 s p S 2 s
17 17 p S = ( P P )u = (I P)( )(u + u ) (due to u = on Ω) p S = u u P( )u p S, p S = p S, u p S, u p S, Pv = p S, u (due to u = on Ω, p S =, Pv) u 2 = p S 2 + p S u 2 (orthogonality)
18 18 p u 2 = p 2 Ω c + ( p u s ) 2 s + ( p u ) 2 s p 2 Ω c + (1 ε) p 2 s s + ( p u ) 2 s + = ( p u ), p s + ( p u ), p s + Hence Φ>s p 2 2 ( p u ), p s 2 ( p u ), p s ( p u ) 2 s + p 2 s ( p u ) 2 s p 2 s + 2( p 2 s p 2 s) junk p u 2 (1 ε) p 2 + 2( p 2 s p 2 s) junk
19 19 D2N/N2D bounds on tubes Ω s = {x Ω Φ(x) < s} Lemma For s > small C such that whenever p = in Ω s and < s < s then n p 2 (I nn t ) p 2 C s p 2 Φ<s Φ<s In the limit s, it reduce to n p 2 (I nn t ) p 2 C Γ Γ Ω p 2 In a half space: n p 2 = (I nn t ) p 2. Known as Rellich identity in 2D circular disk.
20 2 u 2 = u u, u + u 2 (1 ε) u 2 C u 2 = (1 ε)( p S 2 + p S u 2 ) C u 2 (2 3ε) p S 2 C u 2
21 21 Planar domains with corners a negative result Theorem (Cozzi & Pego, 21) Let Ω R 2 be a bounded domain with a locally straight corner. Given any β < 1 and C >, there exists u H 2 H 1 (Ω, R N ) such that Ω p S (u) 2 > β Ω u 2 + C Ω u 2 Open Q: Is extended Navier-Stokes dynamics well-posed in Lipschitz (or polygonal) domains?
22 22 Difference scheme with pressure explicit (Johnston-L, 24) u n+1 u n t p n = (f n u n u n ) in Ω n p n = n f n νn ( u n ) on Ω ν u n+1 = f n u n u n p n u n+1 = on Ω No Stokes solver! u n is not forced to be zero Demonstrated numerical stability, full time accuracy Proved linear stability of tangential velocity in slab geometry by energy estimates Equivalent to gauge method (E & L). Unconditional stability for schemes with time lag BC for Stokes eq (Wang & L, 2)
23 23 Time-discrete scheme with pressure explicit (related to projection methods: Ti96,Pe1,GS3,JL4) u n+1 u n t ν u n+1 = f n u n u n p n, u n+1 = on Ω, p n = ν p n + (I P)(f n u n u n ). S We have the pressure estimate, for fixed constant β < 1 p n ν p n + f n u n u n S νβ u n + f n u n u n + C u n
24 24 Stability analysis in H 1 : dot with u n+1 u n+1 u n t ν u n+1 = f n u n u n p n u n+1 u n 2 + u n+1 2 u n 2 + ν u n t u n+1 (2 f n u n u n + ν p n ) S ( ε ν 2 ) un ε 1 f n u n u n 2 + ν 2 p S 2 This gives u n+1 2 u n 2 t + (ν ε 1 ) u n ε 1 ( f n 2 + u n u n 2 ) + ν p n S 2
25 25 Handling the pressure u n+1 2 u n 2 + (ν ε 1 ) u n+1 2 t 8 ( f n 2 + u n u n 2 ) + ν p n 2 S ε 1 Use the theorem (β = ε): ν pn S 2 νβ u n 2 + νc u n 2 u n+1 2 u n 2 t + (ν ε 1 )( u n+1 2 u n 2 ) + (ν ε 1 νβ) u n 2 8 ε 1 ( f n 2 + u n u n 2 ) + νc u n 2
26 26 Handling the nonlinear term Use Ladyzhenskaya s inequalities to get for N = 2, 3 (standard) u n u n 2 ε 2 u n 2 + 4C ε 2 u n 6 Take ε 1, ε 2 small so ε := (ν ε 1 νβ ε 2 ) > to get: u n+1 2 u n 2 + ε u n 2 C( f n 2 + u n 2 + u n 6 ) An easy discrete Gronwall inequality leads to:
27 27 Stability theorem for N = 2, 3 Theorem Take f L 2 (, T ; L 2 (Ω, R N )), u H 2 H 1 (Ω, R N ). Then positive constants T and C depending only upon Ω, ν and M := u 2 + ν t u 2 + T f 2, so that whenever n t T we have sup u k 2 + k n n 1 k= ( u k 2 + u k+1 u k 2) t C. t
28 28 Classic inf-sup condition for steady Stokes flow FEM Find u X H(Ω, 1 R N ), p Y L 2 (Ω)/R so that u, v + p, v = f, v v X, u, q = q Y. Well-known: inf sup p Y v X p, v p v c >. (1) Condition (1) is necessary and sufficient for stability of (p, u) in L 2 H 1 (the weak NS solution space). Many simple finite elements fail this condition (e.g., equal-order Lagrange)
29 29 Steady Stokes flow in terms of Stokes pressure Find u X H 2 H(Ω, 1 R N ), p Y H 1 (Ω)/R so that u, v + p, v = f, v v X, p, q u, n q Γ = f, q q Y. Stability of (p, u) H 1 (H 2 H 1 (Ω, R N )) (strong solution space) is not subject to standard inf-sup velocity-pressure compatibility. Still a mixed formulation computationally for steady flow Possible gain in simplicity of discretization with C elements (e.g. for complex fluid-coupled systems) nonconforming for H 2
30 3 C Finite element schemes (Johnston-Liu, 4) C finite element scheme: X h H(Ω, 1 R N ), all v h X h and q h Y h, we require require Y h H 1 (Ω)/R. For p n h, q h = f n u n h u n h, q h + ν u n h, n q h Γ, un+1 h u n h t, v h + ν u n+1 h, v h = f n u n h u n h p n h, v h. This scheme has fully accuracy for both velocity and pressure when there are enough resolution and regularity for the solution. We have no theorem on why it works. For C 1 FEM without inf-sup compatibility, however, we have fully theory: stability and error estimates (both for velocity and pressure).
31 31 Error estimates of C 1 FE scheme u n h X h H 2 H 1 (Ω, R N ), p n h Y h H 1 (Ω)/R. For all v h X h and q h Y h, require 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 Γ, Theorem Assume Ω is a bounded domain in R N (N=2,3) with C 3 boundary. Let M, >, and let T > be given by the stability theorem. Let m 2, m 1 be integers, and assume For any v H m+1 H 1 (Ω, R N ) and q H m (Ω), inf v h X,h (v v h ) C h k 1 v H k+1 for 2 k m,
32 32 inf (q q h ) C h m 1 q q h Y H m, h where C > is independent of v, q and h. Then there exists C 1 > with the following property. Whenever u h X h, < h < 1, < n t min(t, T ), and n u h 2 + t u h 2 + f(t k ) 2 t M, then e n = u(t n ) u n h, rn = p(t n ) p n h of C1 finite element scheme satisfy sup e k 2 + k n k= n ( e k 2 + r k 2) t k= C 1 ( t 2 + h 2m 2 + h 2m 2 + e 2 + e 2 t).
33 33 3. Split-step time discretizations of NSE Chorin/Temam projection method Given velocity U n = PU n, let F n+1 = f n+1 U n U n and update velocity by solving u n+1 U n = ν u n+1 + F n+1 t in Ω, u n+1 = on Ω, U n+1 = u n+1 φ n+1 = Pu n+1 Why does it work? Apply P. Stokes pressure reappears implicitly: U n+1 U n + ν p S (U n+1 ) = ν U n+1 + PF n+1, t U n+1 = φ n+1 on Ω (slip error)
34 34 Subtract Formally we get a consistent 1st-order scheme, since φ n+1 t ν φ n+1 = νp S (U n+1 ) + Q(f n+1 U n U n ), n φ n+1 = on Ω. thus φ n+1 Ω = O( t). The correct corresponding pressure p n+1 = νp S (U n+1 ) + Q(f n+1 U n+1 U n+1 ) = φn+1 t ν φ n+1 + O( t) The (old) pressure approximation φ n+1 / t has a boundary layer with O(1) gradient error on Ω, since n φ n+1 = but n p n+1 = n (ν U n+1 + f n+1 U n+1 U n+1 ).
35 35 Higher-order accurate time discretization Chorin-Temam method: 1st-order in time, boundary layers for p Brown, Cortez, Minion 1: Explain 2nd-order time accuracy Guermond, Minev, Shen 6: Review many projection methods, FEM use higher-order backward time difference formulas, and improve boundary-layer accuracy using various strategies, e.g.: Pressure approximation (Karniadakis et al 91): extrapolate curl-curl BCs to approximate p n+1 in a RHS F n+1 Pressure update (van Kan 86, BCG 89): predict then correct pressure, avoiding curl-curl BCs Slip correction (Kim-Moin 86): extrapolate to improve BCs for U n+1
36 36 A stable (3,2) slip-corrected projection method Use BD3 formula D 3 u n+1 = 1 t ( α u n+1 + α k u n+1 k) k 1 Write H j = U j U j & extrapolate: E 3 H n+1 = 3H n 3H n 1 + H n 2, E 2 φ n+1 = 2φ n φ n Given U j u(j t) (j n), solve (without any pressure!) 1 ( α u n+1 t + k 1 α k U n+1 k) = ν u n+1 + f n+1 E 3 H n+1, u n+1 = g n+1 + E 2 φ n+1 on Γ. 2. U n+1 = u n+1 φ n+1, φ n+1 = u n+1, n φ n+1 =.
37 37 Formal accuracy: Apply P and recall Pu j = PU j = U j R(g j ), P u n+1 = Pu n+1 p S (u n+1 ) = U n+1 p S (U n+1 ) : D 3 U n+1 + p n+1 = ν U n+1 + f n+1 E 3 H n+1, p n+1 = νp S (U n+1 ) R(D 3 g n+1 ) + Qf n+1 QE 3 H n+1, U n+1 = g n+1 (φ n+1 E 2 φ n+1 ) on Γ. ( α ) t ν (φ n+1 E 2 φ n+1 ) = p n+1 2 p n + p n 1 = O( t 2 ), hence U n+1 g n+1 = O( t 3 ) formally and the scheme is formally 3rd-order accurate overall.
38 38 A stable (3,2) pressure approximation method Update intermediate velocity u n via (H n = u n u n ) F = f n+1 E 3 H n+1 1 t j 1 α j u n+1 j, P = νp S (E 2 u n+1 ) α t R(gn+1 ) + QF α t un+1 ν u n+1 = P + F u n+1 = g n+1 on Γ The div-free velocity and consistent pressure are U n+1 = u n+1 φ n+1, φ n+1 = u n+1, n φ n+1 =, p n+1 = P ν u n+1. (Improves over previous methods)
39 39 (3,3) pressure update method (cf. BCG 89, Ren et al 5) 1. Find u n+1 to solve 1 ( α u n+1 t + k 1 α k U n+1 k) + E 3 p n+1 = ν u n+1 + f n+1 E 3 H n+1, u n+1 = g n+1 on Γ. 2. U n+1 = u n+1 φ n+1, φ n+1 = u n+1, n φ n+1 =. 3. Update pressure via p n+1 = E 3 p n+1 + ( α t ν ) φ n+1
40 4 Single mode Stokes test problem in strip 1 < x < 1 t u + p S (u) = u, u x=±1 = u(t, x, y) = e iξy σt (u(x, ξ), iv(x, ξ)) For the case of odd symmetry: p = e iξy σt sinh ξx, ( ) cosh ξx cos µx u(x) = A cosh ξ cos µ ( ) sinh ξx sin µx v(x) = B sinh ξ sin µ, A = ξ cosh ξ ξ 2 + µ 2,, B = ξ sinh ξ ξ 2 + µ 2. ξ tanh ξ + µ tan µ =, σ = ξ 2 + µ 2. Usually take ξ = 1, µ , σ
41 41 Single-mode time-discrete stability tests A u n+1 + A 1 u n A k u n+1 k =. Look for u n = κ n u with I I A 3 A 2 A 1 u κu κ 2 u = κ I I A u κu κ 2 u. Solve for eigenvalues κ using Matlab eigs. Plot max κ vs t for various combinations of finite-element velocity/pressure pairs and orders of accuracy.
42 42 Slip-correction schemes, single-mode stability test (3,2) P1/P1 (3,2) P2/P1 (2,2) P1/P1 (2,2) P2/P1 (4,3) P1/P1 (4,3) P2/P t 1: max κ vs t. 3 elements for each var. Solid lines: space-continuous theory. (2,2), (3,2): unconditional stability. (cf. Leriche et al 6 w/pa) (3,3), (4,3): stability for t < t c independent of h, ξ. (smooth?)
43 43 Pressure update schemes, single-mode stability test (3,2) P1/P1 (3,2) P2/P1 (2,2) P1/P1 (2,2) P2/P1 (3,3) P1/P1 (3,3) P2/P1 (4,3) P1/P1 (4,3) P2/P t 2: max κ vs t. 3 elements for each var. Solid lines: interpolated. All P1/P1 PU schemes: unstable. P2/P1 velocity-pressure pairs: (2,2), (3,2): unconditional stability. (3,3),(4,3): stability window
44 44 2D test problem: stability & accuracy Ω = square [ 1, 1] 2 hole, or ring. ν = 2, t = 2 u(t) = g(t) ( ) cos 2 (πx/2) sin(πy) sin(πx) cos 2, (πy/2) p(t) = g(t) cos(πx/2) sin(πy/2).
45 45 2D test: stability in domains with & without corners t c \ N 1 2 (3,3) P1/P (3,3) P2/P (4,3) P1/P (4,3) P2/P Largest time step for linear stability of PA schemes in a square with hole. N = # of refinements from grid in figure. Diffusive constraint on t. t c \ N 1 2 (3,3) P1/P (3,3) P2/P (4,3) P1/P (4,3) P2/P Largest time step for linear stability of PA schemes in annulus. N = # of refinements from grid in figure. No diffusive constraint on t.
46 46 Accuracy: slip-corrected P4 FEM, square with hole log 1 E (& local order α) vs. t E \ t p ˆp h L (2.71) 5.8 (2.98) 5.99 (3.2) (p ˆp h ) L (3.73) 4.39 (3.4) 5.17 (2.62) u u h L 5.75 (2.73) 6.65 (3.2) 7.57 (3.5) (u u h ) L 3.94 (3.22) 4.95 (3.35) 5.83 (2.92) Top: (3,2) SC scheme. Bottom: (2,1) SC scheme. E \ t p ˆp h L (1.97) 2.64 (2.2) 3.26 (2.4) (p ˆp h ) L (1.96) 2.15 (2) 2.76 (2.2) u u h L 3.61 (1.96) 4.22 (2.2) 4.84 (2.5) (u u h ) L 2.19 (1.97) 2.8 (2.1) 3.41 (2.4)
47 47 Pressure error for square with hole, (3,2) schemes x 1 5 x Pressure error for the (3,2) PA (left) and SC (right) schemes. t = P4 elements for each variable. Only values at vertices of triangles are used in the plots.
48 48 Benchmark tests with equal-order C finite elements (u,v)=(1,) (u,v) given (u,v) given 3: Mesh used for backward facing step flow when ν = 1/6 and for flow past a cylinder when ν = 1/1.
49 .5 49 Driven Cavity, Re=1, t= stretched rectangular grid 8192 piecewise linear C elements for each variable. h min =.594, t =.6, (3,2) slip-corrected scheme. Left: vorticity contours. Right: pressure contours.
50 5 Backward-facing step.5 S X (3,2) SC scheme, 664 P1 elements for each variable (dof=3487). ν = 1/1, t = 2, h min =.783, t =.6..5 X 1 X 2 S X (3,2) SC scheme, 17 P2 elements for each variable (dof=3925). ν = 1/6, t = 12, h min =.186, t =.3.
51 51 Flow past cylinder, ν = 1/
52 52 Drag, lift, pressure drop t(c d,max )=3.9348, c d,max = t t(c l,max )=5.6932, c l,max = t p max =2.3219, p(8)= t (3,2) SC scheme 763 isoparametric P4 elements (dof=6322) for each variable.
53 53 Summary The Navier-Stokes pressure for strong solutions in C 3 domains is determined by the current velocity and data by p = νp S (u) R( t g) + Q(f u u). The Laplace-Leray commutator is p S (u) and is strictly controlled by u at leading order, when u = on the boundary. Improved accuracy (esp. for pressure) in numerical computations Improved flexibility of finite-element methods inf-sup conditions need not be respected. Short-time wellposedness extends to u with simple proof
54 Thank you! 54
In honor of Prof. Hokee Minn
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
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 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 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 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 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 informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
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 informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationProjection Methods for Rotating Flow
Projection Methods for Rotating Flow Daniel Arndt Gert Lube Georg-August-Universität Göttingen Institute for Numerical and Applied Mathematics IACM - ECCOMAS 2014 Computational Modeling of Turbulent and
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
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 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 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 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 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 informationNavier-Stokes equations in thin domains with Navier friction boundary conditions
Navier-Stokes equations in thin domains with Navier friction boundary conditions Luan Thach Hoang Department of Mathematics and Statistics, Texas Tech University www.math.umn.edu/ lhoang/ luan.hoang@ttu.edu
More informationON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS
Submitted to Math. Comp. ON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS J.L. GUERMOND 1 AND JIE SHEN 2 Abstract. In this paper we study the rotational form of the pressure-correction
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 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 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 informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
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 informationINCOMPRESSIBLE FLUIDS IN THIN DOMAINS WITH NAVIER FRICTION BOUNDARY CONDITIONS (II) Luan Thach Hoang. IMA Preprint Series #2406.
INCOMPRESSIBLE FLUIDS IN THIN DOMAINS WITH NAVIER FRICTION BOUNDARY CONDITIONS II By Luan Thach Hoang IMA Preprint Series #2406 August 2012 INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY OF
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
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 informationON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS
MATHEMATICS OF COMPUTATION Volume 73, Number 248, Pages 1719 1737 S 0025-5718(03)01621-1 Article electronically published on December 19, 2003 ON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
More informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
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 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 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 informationMathematical analysis of the stationary Navier-Stokes equations
Mathematical analysis of the Department of Mathematics, Sogang University, Republic of Korea The 3rd GCOE International Symposium Weaving Science Web beyond Particle Matter Hierarchy February 17-19, 2011,
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 informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
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 informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More informationKey words. Incompressible magnetohydrodynamics, mixed finite element methods, discontinuous Galerkin methods
A MIXED DG METHOD FOR LINEARIZED INCOMPRESSIBLE MAGNETOHYDRODYNAMICS PAUL HOUSTON, DOMINIK SCHÖTZAU, AND XIAOXI WEI Journal of Scientific Computing, vol. 40, pp. 8 34, 009 Abstract. We introduce and analyze
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
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 informationR. M. Brown. 29 March 2008 / Regional AMS meeting in Baton Rouge. Department of Mathematics University of Kentucky. The mixed problem.
mixed R. M. Department of Mathematics University of Kentucky 29 March 2008 / Regional AMS meeting in Baton Rouge Outline mixed 1 mixed 2 3 4 mixed We consider the mixed boundary value Lu = 0 u = f D u
More informationSTABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS
STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS RAYTCHO LAZAROV AND XIU YE Abstract. In this paper, we derive two stabilized discontinuous finite element formulations, symmetric
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 informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
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 informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More information1. Introduction. Consider nonstationary Navier Stokes equations with a nonslip boundary condition
SIAM J. SCI. COMPUT. Vol. 3, No., pp. 398 49 c 8 Society for Industrial and Applied Mathematics LONG TIME NUMERICAL SOLUTION OF THE NAVIER STOKES EQUATIONS BASED ON A SEQUENTIAL REGULARIZATION FORMULATION
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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
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 informationGlobal regularity of a modified Navier-Stokes equation
Global regularity of a modified Navier-Stokes equation Tobias Grafke, Rainer Grauer and Thomas C. Sideris Institut für Theoretische Physik I, Ruhr-Universität Bochum, Germany Department of Mathematics,
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationOn Some Variational Optimization Problems in Classical Fluids and Superfluids
On Some Variational Optimization Problems in Classical Fluids and Superfluids Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
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 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 informationProjection Method II : Godunov-Ryabenki Analysis
Projection Method II : Godunov-Ryabenki Analysis Weinan E 1 and Jian-Guo Liu 2 Abstract: This is the second of a series of papers on the subject of projection methods for viscous incompressible flow calculations.
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 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 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 informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
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 informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationThe Virtual Element Method: an introduction with focus on fluid flows
The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P.
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 informationVALÉRIE PERRIER AND. Key words. Divergence-free wavelets, Navier-Stokes simulation, physical boundary conditions
DIVERGENCE-FREE WAVELET PROJECTION METHOD FOR INCOMPRESSIBLE VISCOUS FLOW SOULEYMANE KADRI HAROUNA AND VALÉRIE PERRIER Abstract. We present a new wavelet numerical scheme for the discretization of Navier-Stokes
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 informationDIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER
Conference on Turbulence and Interactions TI26, May 29 - June 2, 26, Porquerolles, France DIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER E. Leriche, Laboratoire d Ingénierie
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationRecent progress on the vanishing viscosity limit in the presence of a boundary
Recent progress on the vanishing viscosity limit in the presence of a boundary JIM KELLIHER 1 1 University of California Riverside The 7th Pacific RIM Conference on Mathematics Seoul, South Korea 1 July
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 informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationAdvances in the mathematical theory of the finite element immersed boundary method
Advances in the mathematical theory of the finite element immersed boundary method Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 12, 2014
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 informationTrefftz-Discontinuous Galerkin Methods for the Time-Harmonic Maxwell Equations
Trefftz-Discontinuous Galerkin Methods for the Time-Harmonic Maxwell Equations Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationVALÉRIE PERRIER AND. Key words. Divergence-free wavelets, Navier-Stokes simulation, physical boundary conditions
DIVERGENCE-FREE WAVELET PROJECTION METHOD FOR INCOMPRESSIBLE VISCOUS FLOW SOULEYMANE KADRI HAROUNA AND VALÉRIE PERRIER Abstract. We present a new wavelet numerical scheme for the discretization of Navier-Stokes
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
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 informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
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 informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
MATHEMATICS OF COMPUTATION Volume 75, Number 256, October 2006, Pages 1659 1674 S 0025-57180601872-2 Article electronically published on June 26, 2006 ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED
More informationDivergence Boundary Conditions for Vector Helmholtz Equations with Divergence Constraints
Divergence Boundary Conditions for Vector Helmholtz Equations with Divergence Constraints Urve Kangro Department of Mathematical Sciences Carnegie Mellon University Pittsburgh, PA 15213 Roy Nicolaides
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
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 informationASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 11, Number 1, July 004 pp. 189 04 ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS Tian Ma Department of
More informationGeneralized Newtonian Fluid Flow through a Porous Medium
Generalized Newtonian Fluid Flow through a Porous Medium V.J. Ervin Hyesuk Lee A.J. Salgado April 24, 204 Abstract We present a model for generalized Newtonian fluid flow through a porous medium. In the
More informationJ.-L. Guermond 1 FOR TURBULENT FLOWS LARGE EDDY SIMULATION MODEL A HYPERVISCOSITY SPECTRAL
A HYPERVISCOSITY SPECTRAL LARGE EDDY SIMULATION MODEL FOR TURBULENT FLOWS J.-L. Guermond 1 Collaborator: S. Prudhomme 2 Mathematical Aspects of Computational Fluid Dynamics Oberwolfach November 9th to
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 informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
More informationTime domain boundary elements for dynamic contact problems
Time domain boundary elements for dynamic contact problems Heiko Gimperlein (joint with F. Meyer 3, C. Özdemir 4, D. Stark, E. P. Stephan 4 ) : Heriot Watt University, Edinburgh, UK 2: Universität Paderborn,
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 informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More information