IFE for Stokes interface problem

Transcription

1 IFE for Stokes interface problem Nabil Chaabane Slimane Adjerid, Tao Lin Virginia Tech SIAM chapter February 4, 24 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 / 2

2 Problem statement Consider the Stokes problem where the Stress tensor is S(u, p) = f, in Ω, (a) u =, in Ω, (b) u = g, on Ω, (c) S lj (u, p) = ν(d l u j + D j u l ) pδ lj, l, j =, 2. Ω is occupied by two fluids ν + and ν separated by an interface Γ. The jump conditions are [S(u, p) n] Γ =, (2) [u] Γ =. (3) In this work, we assume that the pressure is continuous. Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 2 / 2

3 Ω Γ Ω + Ω Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 3 / 2

4 Ω Γ Ω + Ω Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 4 / 2

5 IFE basis functions We approximate the interface with a straight line. We map the physical element to the reference element using appropriate mapping. Then, we construct a piecewise vector-valued function ˆΦ s.t: { ˆΦ ˆΦ + if (ˆx, ŷ) ˆT + = ˆΦ if (ˆx, ŷ) ˆT, Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 5 / 2

6 IFE basis functions ˆΦ (ˆx, ŷ) = ( a + b ˆx + c ŷ + d ˆxŷ a 2 + b 2 ˆx + c 2 ŷ + d 2 ˆxŷ ˆΦ + (ˆx, ŷ) = ( a + + b + ˆx + c+ ŷ + d + ˆxŷ a b+ 2 ˆx + c+ 2 ŷ + d + 2 ˆxŷ The IFE basis functions satisfy the following jump constraints ), (4a) ), (4b) ˆΦ (Ê) = ˆΦ + (Ê), ˆΦ j (Â i ) = δ ij, i, j =, 2,..., 8, (5) ˆΦ ( ˆD) = ˆΦ + ( ˆD), 2 ˆΦ ˆx ŷ = 2 ˆΦ + ˆx ŷ, (6) [S( ˆΦ, p) n ]ds =. (7) ˆDÊ ˆDÊ The jump conditions yield 6 constraints that are used to solve for the unknowns. Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 6 / 2

7 Properties of the bilinear IFE space On every interface element T the IFE basis: are well defined and are unique. are continuous. form a partition of unity i.e 4 Φ i (x, y) = i= ( ), 8 Φ i (x, y) = i=4 ( ), (x, y) t T. Furthermore, If the viscosity parameters have no discontinuity i.e ν + = ν, then Φ i becomes the standard vector-valued bilinear nodal basis Ψ i. If min{ T +, T } shrinks to, then Φ i becomes the standard vector bilinear nodal basis Ψ i. The interpolation error converges with optimal rate. Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 7 / 2

8 IFE basis functions Figure: First and second components of the IFE basis functions with ν =, ν + = 5.3 on the left. The first and second components of the standard nodal finite element basis on the right. Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 8 / 2

9 Interpolation error h u I h u,ω Order u 2 I h u 2,Ω Order.6844e-.6844e e e e e e e e e e e-3.37 Table: L 2 error of the interpolation error with ν = and ν + = 5.3 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 9 / 2

10 Interpolation error h u (I h u),ω Order u 2 (I h u) 2,Ω Order e e e e e e e e e e Table: L 2 error of the interpolation error with ν = and ν + = 5.3 h u (I h u),ω Order u 2 (I h u) 2,Ω Order e e e e e e e e e e Table: H error of the interpolation error with ν = and ν + = 5.3 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 / 2

11 IFE spaces using Q 2 iso Q /Q elements Consider an arbitrary uniform cartesian mesh T h and its refinement T h/2 and define the following on T h/2 { span{φ i (x, y), i =, 2,..., 8}, T is an interface element S h/2 (T ) = span{ψ i (x, y), i =, 2,..., 8}, otherwise, and define S h/2 (Ω) = T Th/2 S h/2 (T ), and S h/2, (Ω) = {v S h/2 (Ω) : v e =, e Eh/2 B \Ei h/2 }, where Eh/2 i is the set of interface edges. On the mesh T h, define M h (Ω) = {q C L 2 : q T Q }. Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 / 2

12 IFE formulation The finite element formulation is: Find (u h, p h ) S h (Ω) M h (Ω) s.t 2 S(u h, p h ) : v h dx = f v h dx, v h S h,(ω) (8a) 2 Ω Ω u h q h dx =, q h M h (Ω), (8b) Ω Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 2 / 2

13 Numerical results The domain Ω is the square [, ] 2, the interface is a half-circle with radius r =.3, and center x =.4327, y =.329 that separates Ω into two regions Ω + = {(x, y) t : (x x ) 2 + (y y ) 2 > r 2 } and Ω = {(x, y) t : (x x ) 2 + (y y ) 2 < r 2 }. The function u is defined as { u + u + (x, y) = = ((x x ) 2 + (y y ) 2 r 2 )(y y ), u 2 + = ((x x ) 2 + (y y ) 2 r 2, (9a) )(x x ), u = ν+ ν u+. (9b) Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 3 / 2

14 Numerical results h u u,h,ω u 2 u 2,h,Ω p p h,ω e e e e e e e e e e e e e e e e e e-2 Table: L 2 errors with ν = and ν + = 5.3 for the non-penalized method. u u,h,ω Ch.587, u 2 u 2,h,Ω Ch.6484, p p h,ω Ch Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 4 / 2

15 Figure: Standard Lagrange finite element basis Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 5 / 2

16 Figure: Immersed finite element basis Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 6 / 2

17 IFE formulation The finite element formulation is: Find (u h, p h ) S h (Ω) M h (Ω) s.t 2 S(u h, p h ) : v h dx Ω e Eh i + e E i h \EB h e e E i h \EB h [u h ] [v h ]ds = Ω Ω e {(S(u h, p h ) n)} [v h ]ds e {(S(v h, ) n)} [u h ]ds f v h dx, (a) v h S h,(ω) (b) 2 u h q h dx =, q h M h (Ω), (c) Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 7 / 2

18 Numerical results The domain Ω is the square [, ] 2, the interface is a half-circle with radius r =.3, and center x =.4327, y =.329 that separates Ω into two regions Ω + = {(x, y) t : (x x ) 2 + (y y ) 2 > r 2 } and Ω = {(x, y) t : (x x ) 2 + (y y ) 2 < r 2 }. The function u is defined as { u + u + (x, y) = = ((x x ) 2 + (y y ) 2 r 2 )(y y ), u 2 + = ((x x ) 2 + (y y ) 2 r 2, (a) )(x x ), u = ν+ ν u+. (b) Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 8 / 2

19 Numerical results h u u,h,ω Order u 2 u 2,h,Ω Order.592e e e e e e e e e e e e-6.99 Table: L 2 error and order of convergence with ν = and ν + = 5.3 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 9 / 2

20 Numerical results h p p h,ω Order 5.363e e e e e e-4.45 Table: L 2 error and order of convergence with ν = and ν + = 5.3 Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 2 / 2

21 THANK YOU! Nabil Chaabane (VT) IFE for Stokes interface problem February 4, 24 2 / 2