arxiv: v1 [physics.flu-dyn] 4 Sep 2013

Transcription

1 THE THREE-DIMENSIONAL JUMP CONDITIONS FOR THE STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY, SINGULAR FORCES, AND AN INCOMPRESSIBLE INTERFACE PRERNA GERA AND DAVID SALAC arxiv: v1 physics.fl-dyn] 4 Sep 2013 Abstract. The three-dimensional jmp conditions for the pressre and velocity fields, p to the second normal derivative, across an incompressible/inextensible interface in the Stokes regime are derived herein. The flid viscosity is only piecewise continos in the domain while the embedded interface exerts singlar forces on the srrond flids. This gives rise to discontinos soltions in the pressre and velocity field. These jmp conditions are reqired to develop accrate nmerical methods, sch as the Immersed Interface Method, for the soltions of the Stokes eqations in sch sitations. Key words. Stokes eqations, discontinos viscosity, singlar force, immersed interface method, jmp conditions, incompressible interface, inextensible interface AMS sbject classifications. 65N06, 65N12, 76Z05, 76D07, 35R05 1. Introdction. Biological system containing cells can be considered as a mltiphase flid problem, with the cell membranes forming a bondary between flids of differing material properties. Unlike standard mltiphase flid problems, sch as droplets or bbbles, biological mltiphase systems, sch as red blood cells srronded by blood plasma, are characterized by inextensible membranes and singlar forces de to both a srface-tension like contribtion and a resistance to bending 16]. Varios techniqes have been developed to model sch systems. These techniqes can be split into two categories, based on how they track the location of the interface. Lagrangian methods, sch as the bondary integral method 24, 25] explicitly track the location of the interface. While these method can be highly accrate, they are difficlt to extend to three-dimensions and are not well sited when the membrane ndergoes large deformations. Elerian method, on the other hand, implicitly track the location of the interface. Examples of Elerian methods inclde the phase-field method 2, 13] and the level set method 18, 19]. Advantages of the Elerian methods inclde easy extensibility to three dimensions and the ability to handle large membrane deformations natrally. The major disadvantage is the added comptational cost of implicit membrane tracking. Elerian methods can be frther divided into two categories depending on how they treat the singlar forces arising at the embedded membrane. The first type, characterized by the immersed bondary method 15] or the continm srface force method 4], distribte (smooth) the singlar forces over a small region near the interface. This essentially trns the membrane forces into localized body-force terms, which are then inclded in the flid governing eqations. For most sitations this smoothing reslts in a first-order accrate method for the flid eqations 23]. The second type of Elerian method avoids the smearing of the interface by explicitly incorporating the singlar forces and embedded bondary conditions in the field eqations. These forces and immersed bondary conditions are inclded by explicitly considering the jmp in the soltion and derivatives of the soltion across the interface. Mayo 14] sed this idea to solve the Poisson and Biharmonic eqation This work was spported by the University at Bffalo SUNY and NSF Award # Department of Mechanical and Aerospace Engineering, University at Bffalo SUNY, Bffalo, NY,

2 on irreglar domains. Leveqe and Li called this techniqe the Immersed Interface Method (IIM) and sed it to solve Elliptic eqations 9]. They later extended this to Stokes flow with elastic bondaries or srface tension 8]. Later the immersed interface method was frther developed to model the Navier-Stokes eqations 7, 12, 27]. Interested readers are referred to the book by Li and Ito for more information abot the immersed interface method 11]. The application of the immersed interface method to a mltiphase flid system reqires that the jmps in the velocity and pressre field be explicitly calclated. For standard mltiphase systems these jmp conditions for a discontinos viscosity across the interface have been determined both in two dimensions 22, 23] and three dimensions 5, 28]. Only recently has work been done on extending the immersed interface method to mltiphase flows with inextensible membranes 21]. This recent work, thogh, is limited to the two-dimensional, constant-viscosity case and is not applicable to the more general discontinos viscosity sitation. In this work the three-dimensional jmp conditions for a mltiphase Stokes flow system with an inextensible interface, singlar forces, and a discontinos viscosity are for the first time derived. The jmp p to the second normal derivative of both the pressre and velocity fields will be developed. These jmp conditions will be sed to constrct an Immersed Interface solver. Several analytic test cases are also developed to verify the accracy of the jmp conditions. This paper is organized as follows. First, the basic governing eqations for the three-dimensional Stokes system with an inextensible membrane are presented. Next, in Sec. 3 the immersed interface method is briefly otlined. The jmp conditions for the velocity and pressre in the Stokes eqations are developed in Sec. 4. The derivation of analytic test cases is shown in Sec. 5 while reslts and convergence analysis follow in Sec. 6. A brief conclsion and ftre work follows in Sec Governing Eqations. Let a three-dimensional bonded domain, Ω, contain a closed and incompressible material interface Γ. For simplicity assme that the interface is wholly contained in the domain. The region enclosed by the interface is denoted as Ω while the region otside the interface is denoted by Ω +, Fig The flid in each domain is modeled sing the Stokes eqations: (2.1) (2.2) p ± =µ ± 2 ± + g ±, ± =0, with bondary conditions + = b on the bondary of Ω. The pressre, p, body force term, g, and viscosity, µ, are all piecewise constant in the domain, with a finite jmp occrring across the interface. The velocity field, = (, v, w), is assmed to be C 0 continos across the interface. It is also assmed that the velocity obey a no-slip condition on the interface. It is ths possible to state that the velocity of the interface Γ is given by q = + =. Application of the volme incompressibility constraint, Eq. (2.2), to the momentm eqations, Eq. (2.1), leads to a pressre-poisson eqation in each domain, (2.3) 2 p ± = g ±, with a bondary condition on the oter domain of (2.4) p + b = µ + n b n b g +, 2

3 n Γ Ω - n b Ω Fig Ω Sample comptational domain. where n b is the otward facing normal to the domain Ω. The incompressible interface reqires that the velocity on the interface be divergencefree: (2.5) s q = 0, where s = (I n n) is the srface gradient and s is the srface divergence operator. The presence of the interface reslts in a singlar force, f, exerted on the srronding flids. This singlar force is application specific. In vesicle/cell simlations this force wold be composed of a variable srface tension, γ, and other singlar forces, f 0 : (2.6) f = γhn + s γ + f 0, where n is the otward facing normal to the interface Γ (from Ω into Ω + ) and H is the total crvatre. An example of other singlar forces which might be applied are the bending forces in vesicle and red-blood cell simlations 19, 20] while the srface tension will be chosen to enforce the srface incompressibility constraint, Eq. (2.5). Ftre work will consider singlar forces of this form. The singlar force is balanced by a jmp in the normal component of the stress tensor: (2.7) ( σ + σ ) n = f, where the stress tensor is given by (2.8) σ ± = p ± I + µ ± ( ± + ( ±) T ). 3. Immersed Interface Method. The governing eqations as described in Sec. 2 reslt in two copled, bt discontinos, flid fields. In this section the scheme developed in Ref. 17] and sed to solve the two-dimensional Stokes eqations in Ref. 6] is reviewed. This is also the impets for deriving the normal jmp conditions across the interface. Begin by discritizing the domain Ω sing a niform Cartesian mesh with a grid spacing of h in all directions. While this section will focs on two-dimensional domains 3

4 the extension to three-dimensions is straight-forward and sed in the sbseqent sections. Now consider the soltion to a generic Poisson problem of the form (3.1) 2 = f, where both and f are (possibly) discontinos across the interface. Upon discritizing Eq. (3.1) it is possible to classify all grid points as either reglar or irreglar, Fig Reglar nodes are those nodes where the interface does not cross the discritized stencil. As the soltion is continos away from the interface the discretization at reglar nodes does not need to be modified. Irreglar nodes are defined as those where the interface does cross the discritized stencil. As the soltion may be discontinos across the interface it is not possible to se the standard discretization. Instead the system is modified at irreglar points to take into accont the discontinity. f x i,j G Fig Discritized Laplacian at a reglar node (circle) and irreglar node (sqare). Irreglar nodes are those where the interface intersects the stencil. The signed distance fnction from any grid point to the interface is given by φ. Consider the irreglar point x i,j shown in Fig Let this point exist in the Ω + domain. A standard second-order finite difference discretization reslts in (3.2) + i 1,j + + i+1,j + + i,j i,j+1 4+ i,j h 2 = f + i,j. The isse is the vales at x i 1,j and x i,j 1 exist in the Ω domain, not the Ω + domain. This fact mst be taken into accont in the discretization. Let a jmp across the interface be defined as (3.3) ] = lim ɛ 0 + (x Γ + ɛn) lim ɛ 0 + (x Γ ɛn), where x Γ is a point on the interface and n is the otward normal vector. Assme that the jmps in the soltion, ] Γ, the first normal derivative, /] Γ, and the second normal derivative 2 / 2 ] Γ are known on the interface for the problem given in Eq. (3.1). Using a Taylor Series expansion in the normal direction abot an interface point the jmps may be extended to a grid point by ] (3.4) ] i,j = ] Γ + φ i,j Γ 4 + φ2 i,j 2 2 ] 2 + O(h 3 ), Γ

5 where φ i,j is the signed-distance from the grid point x i,j to the interface, Fig Note that by sing the signed distance fnctions it is possible to accont for extension of the jmps into either domain. By extending the soltion jmps from the interface to the grid points it can be written that + i,j = i,j + ] i,j. Use this expression in the Eq. (3.2) to get (3.5) i 1,j + + i+1,j + i,j i,j+1 4+ i,j h 2 + ] i 1,j + ] i,j 1 h 2 = f + i,j. As the jmps ] i 1,j and ] i,j 1 are known they can be moved to the right-hand side as explicit corrections on the linear system, (3.6) i 1,j + + i+1,j + i,j i,j+1 4+ i,j h 2 = f + i,j + C i,j, where C i,j is the total correction needed at an irreglar node at location x i,j. In this case C i,j = (] i 1,j + ] i,j 1 )/h 2. This method is easily extended to irreglar nodes on either side of the interface and to three-dimensional systems. It shold be noted that the extension of the jmps mst be calclated p to third-order accracy to ensre that irreglar nodes have a local trncation error of O(h). Despite this higher local trncation error the overall method will retain the second-order accracy of the nderlying discretization 1, 8, 10]. If the second-normal derivative is not taken into accont the local trncation error for irreglar nodes will be redced to O(1) and the overall scheme wold only be first-order. 4. Derivation of Jmp Conditions. The governing eqations in Sec. 2 reslt in two copled Poisson problems. Using the Immersed Interface Method of Sec. 3 the two systems can be written over the entire domain by inclding the appropriate corrections: (4.1) (4.2) 2 (µ) = p g + C, 2 p = g + C p, where the corrections C and C p are non-zero only at irreglar grid points. To calclate these corrections at irreglar grid points the jmps p to the second-normal derivatives need to be derived for both the pressre and the agmented velocity, µ. This section will otline the derivation of these jmp conditions. It shold be noted that a C 0 -continos velocity indicates that ] = 0. It shold also be noted that the interface can be approached from either the Ω or Ω + domains. It is ths defined that (4.3) (4.4) f = lim f(x Γ + ɛn) ɛ 0 + f + = lim f(x Γ ɛn), ɛ 0 + where f is a qantity, sch as velocity, of interest. Let n be the otward (pointing into Ω + ) nit normal while t and b are two nit tangent vectors chosen sch that (n, t, b) form a Darbox Frame. Also chose t and b to lie along the principle directions of the interface, Fig Denote the directions 5

6 in the normal and two tangent vectors as n, t, and b, respectively. In the following sections all derivative in a particlar direction are denoted as /, / t, and / b. Note that a lack of bold-face indicates a direction and not a vector. Theorems 4.1 to 4.3 are non-trivial and non-obvios relations that will be sed to develop the jmp conditions for the pressre and velocity. n b t Fig The local vectors on an interface point. The nit vector n is the otward facing normal while t and b are nit tangent vectors which lie along the principle directions at the point on the interface. Theorem 4.1. The jmp in the viscosity times the normal derivative of the velocity in the normal direction is zero: (4.5) µ ] n = 0 Proof. For the flid in each domain the incompressibility conditions holds, µ ± ± = 0. Ths, the jmp in the incompressibility condition is zero across the interface. In terms of the local vectors this is (4.6) µ ( n + t t + b b )] = 0. The velocity as yo approach the interface mst also be srface-divergent free, µ ± s ± = 0. Therefore the the jmp in the srface divergence of the velocity mst also be zero, ( (4.7) µ t t + )] b b = 0. Combining Eqs. (4.6) and (4.7) reslts in ( ) ] (4.8) µ n As the normal vector is continos across the interface, n = n +, the jmp condition, Eq. (4.5), is obtained. Theorem 4.2. The following expression holds on the interface, ( 2 )] (4.9) µ 2 n + 2 t t + 2 b b = 0. 6 = 0.

7 Proof. The incompressibility constraint, = 0, mst hold identically in both domains. Ths, the normal derivative of this constraint mst also be zero, ( ) / = 0. It is ths possible to write µ ] ( ) = µ ( n + t t + )] b b ( 2 )] = µ 2 n + 2 t t + 2 b b (4.10) =0, where the fact that /] = 0, t/] = 0, and b/] = 0 are held to be tre. Theorem 4.3. The following expression holds on the interface, (4.11) µ 2 s ] ( 2 )] n = µ t b 2 n. Proof. The fll Laplacian in the local frame can be written as (4.12) 2 = t b 2, while the srface Laplacian can be written as 26] (4.13) 2 s = H. Ths the following holds, (4.14) µ 2 s ] n = ( 2 )] µ t b 2 n H µ ] n. Using the reslt of Theorem 4.1 the expression in Eq. (4.11) is shown to hold tre. The next three theorems otline the derivation of the pressre jmp conditions. Theorem 4.4. The jmp in the pressre is given by (4.15) p] = f n. Proof. The balance of forces on the interface, Eq. (2.7), is ( ) ] (4.16) pn + µ + ( ) T n = f. Taking the normal component of this balance reslts in (4.17) p] = 2 µ ] n f n. dn Using Eq. (4.5) reslts in the jmp of the pressre field, Eq. (4.15). 7

8 as (4.18) Theorem 4.5. The jmp in the normal derivative of the pressre may be written ] p = µ 2 ] n + g] n, or alternatively as (4.19) ] ( ) p =2µ] n 2 q sq + κ t t t + κ q b b b s f + (f n) H + g] n, where κ t and κ b are the principle crvatres along the t and b directions and recalling that q is the velocity on the interface. Proof. As Eq. (2.1) holds in each domain the jmp in the system mst be zero, p] = µ 2 ] + g]. Dotting this by n reslts in Eq. (4.18). As has been noted elsewhere 5] this form is not sefl for nmerical simlations de to the second-order partial derivatives of the velocity. To obtain the alternative form begin by taking the srface divergence of the force balance on the interface, s σ n] = s f, or in expanded form: ( ) ] (4.20) s pn] + s µ + ( ) T n = s f. Noting that the jmp operator commtes with differentiation along the interface 5, 28], individal portions of Eq. (4.20) can be written as (4.21) s pn] = s (pn)] = s p n + p s n] = p] s n = p]h = (f n) H. Next, write the srface divergence of a vector v on the interface as (4.22) s v = v t t + v b b. Then, (4.23) s µ n] = µ s ] ( = µ t ( ) t + b ( ) )] b When considering the derivatives of / = ( ) n along the tangent directions it is necessary to also take derivatives of the normal vector along the same directions. This reslts in the following expression along the t-direction (4.24) t ( ) t = 2 ( ) t t + ( ) T t dt, where the derivative along the b-direction is similar. In the Darbox Frame it can be written that /dt = κ t t + τ t b, where τ t is the geodesic torsion along the t direction. As the t-direction is a principle direction then the geodesic torsion is zero, τ t = 0 3]. Using this reslts in the following expression, (4.25) t ( ) t = 2 t t + κ t 8 t t.

9 Combining this with a similar reslt in the b-direction gives ( 2 )] (4.26) s µ n] = µ t t + κ t t t + 2 b b + κ b b b In a similar fashion the final portion of Eq. (4.20) can be written as ] ( (4.27) s µ ( ) T 2 )] n = µ t 2 n + κ t t t + 2 b 2 n + κ b b b Combining Eqs. (4.20), (4.21), (4.26), and (4.27) and sing Theorem 4.2 and Theorem 4.3 gives the expression (4.28) ( 2 )] ( )] µ 2 n = (f n) H s f + µ 2 s n + 2κ t t t + 2κ b b b. Next, se Eq. (4.12) to split the Laplacian term in the base jmp condition, Eq. (4.18), and combine with Eq to obtain ] ( )] p (4.29) = 2 µ n 2 s + κ t t t + κ b b b s f + (f n) H + g] n. Letting κ g t be the geodesic crvatre along the t direction allows for the following to hold, ] ( )] 2 ] = = t t t t t 2 κg t b κ t 2 ] ] ] = t 2 κ g t κ t b 2 ] ] = t 2 κ g ] t b κ t 2 ] ] (4.30) = t 2 κ t. De to the continity of tangential derivatives, / t] = ] / t = 0 holds on the entire interface. It is ths possible to state that / t / t] = 0. Therefore 2 ] ] (4.31) t 2 = κ t. When Eq. (4.31) is dotted by the nit normal vector and noting that /] n = 0, see Ref. 23], the expression 2 / t 2] n = 0 holds. A similar expression holds in the b direction. The exact form given in Eq. (4.19) can be obtained by writing jmps of the form a b] as a b] = a]b ± + a b]. This allows for the following to hold, (4.32) µ 2 t 2 n ] =µ] 2 q t 2 n + µ 2 t 2 ] n = µ] 2 q t 2 n. Replacing the appropriate terms in Eq. (4.29) reslts in the alternative jmp in the normal pressre derivative. Theorem 4.6. The jmp in the second-normal derivative of the pressre is given by (4.33) 2 ] p 2 = 2 s (f n) + g] 9 ] p H,

10 where the jmp in the first-normal derivative of the pressre is given in Eq. (4.19). Proof. As the pressre-poisson eqation, Eq. (2.3), holds in each domain the jmp is given by 2 p ] = g]. Splitting the Laplacian of the pressre into the local frame reslts in 2 p (4.34) sp + p ] H = g]. Solving for the jmp in the second derivative and sing the reslt of Theorem 4.4 reslts in Eq. (4.33). The final three theorems otline the derivation of the velocity jmp conditions. Theorem 4.7. The jmp in the agmented velocity is given by (4.35) µ] = µ] q. Proof. Expanding the jmp µ] reslts in (4.36) µ] = µ] q + µ ± ] = µ] q, de to ] = 0 and + = = q on the interface. Theorem 4.8. The jmp in the normal derivative of the agmented velocity is given by (4.37) µ ] (( ) ( ) ) q q = P f µ] t n t + b n b where P is the projection operator given by P = I n n. Proof. Dot the force balance on the interface, Eq. (2.7), by the tangent directions to obtain µ ] t =f t µ ] (4.38) n (4.39) µ ] b =f b t µ b ] n In conjnction with Eq. (4.5) a linear system can be written as 0 (4.40) n x n y n z t x t y t z µ ] = f t µ ] n t, b x b y b z f b µ ] n b where n = (n x, n y, n z ) T, t = (t x, t y, t z ) T, and b = (b x, b y, b z ) T are the components of the local vectors. As the inverse of the matrix in Eq. (4.40) is the transpose of the matrix the jmp can be obtained as (4.41) µ ] (( ) ( ) ) q q = (f t) t + (f b) b µ] t n t + b n b. 10

11 Realizing that the term (f t) t + (f b) b is simply the tangential projection of the force f onto the interface the jmp given in Eq. (4.37) is obtained. Theorem 4.9. The jmp in the second normal derivative of the agmented velocity is given by (4.42) ] µ 2 2 = µ] 2 sq µ ] H + ] p n s (f n) g]. Proof. Consider the jmp of the agmented momentm balance eqation, Eq. (4.1). Expand the Laplacian of the velocity and the gradient of the pressre into normal and tangential components: (4.43) µ µ 2 s + µ ] p H = n + p b b + p ] t t g]. De to continity of tangential derivatives, the nit normal, and the total crvatre, it can be stated that ] (4.44) µ 2 2 = µ 2 s ] µ ] ] p H + n + p] b b + p] t t g]. The jmp in the pressre is simply p] = f n while (4.45) µ 2 s] = µ] 2 sq + µ ± 2 s] = µ] 2 sq + µ ± 2 s] = µ] 2 sq. Finally, noticing that ( / t)t + ( / b)b = s is the srface gradient the jmp condition given in Eq. (4.42) can be obtained. 5. Sample Test Case. To verify the accracy of the derived jmp conditions a sample analytic test case has been created. To create this test case valid inner and oter velocity fields, and +, are determined for a pre-determined interface. These velocity fields are chosen to enforce the conditions = + and s = s + = 0 on the interface Γ and = + = 0 in the domain Ω. Next, arbitrary inner and oter pressre fields, p and p + are chosen. Using the velocity and pressre field the inner and oter body forces, g and g + can be calclated sing Eq. (2.1). The singlar force, f is obtained throgh the force balance on the interface, Eq. (2.7). The analytic body and singlar forces are then sed in conjnction with the known interface velocity q and appropriate bondary conditions to solve the discontinos Stokes eqations in the domain sing the Immersed Interface Method. The nmerical pressre and velocity fields are compared to their analytic versions to determine the overall spatial accracy of the method. In this example the interface is a nit sphere centered at the origin. The oter viscosity is taken to be 1 while the inner viscosity is given by µ. The exact inner velocity and pressre for the region x 2 + y 2 + z 2 1 are given by (5.1) (5.2) = 2yz xz xy, p = cos(x) sin(y + z). 11

12 This reslts in an inner body force of the form (5.3) g = sin(x) sin(y + z) cos(x) cos(y + z) cos(x) cos(y + z). The exact oter velocity and pressre for the region x 2 + y 2 + z 2 > 1 is 2yz + 1 (5.4) = xz x 2 + y 2 + z 2, xy (5.5) p + =0, which reslts in a body force of (5.6) g + = 1 12yz 6xz x 2 + y 2 + z 2 6xy A force balance on the interface reslts in a singlar force of 2(1 + µ)yz + x cos(x) sin(y + z) (5.7) f = (1 + µ)xz + y cos(x) sin(y + z) (1 + µ)xy + z cos(x) sin(y + z). Clearly on the interface q = = + holds on Γ and can be applied in the derived jmp conditions. 6. Reslts. In the following reslts the domain is a taken to be a 2, 2] 3 cbe. The pressre and velocity fields are obtained sing a split calclation where the pressre is calclated first by solving Eq. (2.3) and then obtaining the velocity field throgh the soltion of Eq. (2.1). It shold be noted that the jmp conditions allow for any Stokes field soltion techniqe, not only the one sed here. Sample pressre reslts for the z = 0 plane on a grid are shown in Fig. (6.1). The external viscosities shown are and It is clear that the sharp discontinity of the pressre is captred by the method and the reslt is not greatly affected by the viscosity ratio.. 12

13 (a) µ + = (b) µ + = 1000 Fig Example of the pressre distribtion on the z = 0 plane for a grid for oter viscositoes of and The inner viscosity is eqal to 1. 13

14 To demonstrate that the derived jmp conditions, when sed with a secondorder Immersed Interface Method, achieve the desired accracy the errors for both the pressre and velocity fields are now presented. In calclating the error external viscosities ranging from 10 3 to 10 3 have been considered. Grid sizes ranging from 48 3 to have been sed, corresponding to grid spacings of h to h Error e p Viscosity 1.0E E E E E E E h Grid Spacing (a) Error e p 2 Viscosity 1.0E E E E E E E h Grid Spacing (b) Fig Error of the pressre field verss grid spacing. The shown viscosities are the oter viscosities, µ +. The inner viscosity is taken to be µ = 1. 14

15 e Viscosity 1.0E E E E E E E+03 Error h Grid Spacing (a) Error e 2 Viscosity 1.0E E E E E E E h Grid Spacing (b) Fig Error of the velocity field verss grid spacing. The shown viscosities are the oter viscosities, µ +. The inner viscosity is taken to be µ = 1. 15

16 The error for the pressre field are presented in Fig. 6.2 while the corresponding velocity errors are shown in Fig For both the pressre and velocity fields the second-order accracy of the nderlying discretization is clearly preserved despite the presence of the discontinities in the soltion field. It shold be noted that the overall accracy of the method depends on the viscosity difference between the inner and oter flid, with matched viscosity prodcing the lowest error. This is to be expected as the jmp conditions are heavily dependent on this viscosity difference. Discretization errors in the srface derivative qantities, sch as the srface Laplacian of Eq. (4.19), are amplified as the viscosity difference increases. While more accrate srface derivative calclations wold redce this dependence, it will always remain and mst be taken into accont. 7. Conclding Remarks. In this paper the three-dimensional, normal jmp conditions for the pressre and velocity fields have been derived for an embedded incompressible/inextensible interface with a singlar force and in the presence of a discontinos viscosity field. The jmp conditions take into accont the additional constraint of interface incompressibility given by a srface-divergent free velocity field. The jmp conditions are applicable to any Stokes solver sing the Immersed Interface Method. As a demonstration of the accracy of the jmp conditions a sample analytic test case has been created. Error analysis sing this test case demonstrates that a second-order accrate discretization maintains the expected accracy despite the discontinity of the soltion fields. Ftre work will se these jmp conditions to constrct generalized Immersed Interface Methods which are capable of enforcing srface incompressibility while solve the Stokes eqations for mltiphase Stokes flow with discontinos viscosity. REFERENCES 1] Loyce Adams and Zhilin Li, The Immersed Interface/Mltigrid Methods for Interface Problems, SIAM J. Sci. Compt., 24 (2002), pp ] T Biben, K Kassner, and C Misbah, Phase-field approach to three-dimensional vesicle dynamics, Phys. Rev. E, 72 (2005), p ] MP Carmo, Differential geometry of crves and srfaces, Prentice-Hall, ] YC Chang, TY Ho, B Merriman, and S Osher, A level set formlation of elerian interface captring methods for incompressible flid flows, J. Compt. Phys., 124 (1996), pp ] K Ito, Z Li, and X Wan, Pressre Jmp Conditions for Stokes Eqations with Discontinos Viscosity in 2D and 3D, Methods Appl. Anal., 13 (2006), pp ] M.C. Lai and H.C. Tseng, A simple implementation of the immersed interface methods for Stokes flows with singlar forces, Compt. & Flids, 37 (2008), pp ] D.V. Le, B.C. Khoo, and J. Peraire, An immersed interface method for viscos incompressible flows involving rigid and flexible bondaries, J. Compt. Phys., 220 (2006), pp ] RJ Leveqe and ZL Li, Immersed interface methods for Stokes flow with elastic bondaries or srface tension, SIAM J. Sci. Compt., 18 (1997), pp ] Randall J Leveqe and Zhilin Li, The immersed interface method for elliptic eqations with discontinos coefficients and singlar sorces, SIAM J. Nmer. Anal., 31 (1994), pp ] Z. Li and K. Ito, Maximm Principle Preserving Schemes for Interface Problems with Discontinos Coefficients, SIAM J. Sci. Compt., 23 (2001), pp ], The Immersed Interface Method: Nmerical Soltions of PDEs Involving Interfaces and Irreglar Domains, Frontiers in Applied Mathematics, Society for Indstrial and Applied Mathematics, ] Zhilin Li and Ming-Chih Lai, The Immersed Interface Method for the Navier Stokes Eqations with Singlar Forces, J. Compt. Phys., 171 (2001), pp ] JS Lowengrb, A Ratz, and A Voigt, Phase-field modeling of the dynamics of mlticompo- 16

17 nent vesicles: Spinodal decomposition, coarsening, bdding, and fission, Phys. Rev. E, 79 (2009), p ] Anita Mayo, The Fast Soltion of Poisson s and the Biharmonic Eqations on Irreglar Regions, SIAM J. Nmer. Anal., 21 (1984), pp ] Charles S Peskin, Nmerical analysis of blood flow in the heart, J. Compt. Phys., 25 (1977), pp ] C Pozrikidis, Modeling and simlation of capsles and biological cells, CRC Press, 2003, ch. Shell theory for capsles and shells., pp ] D Rssell and ZJ Wang, A cartesian grid method for modeling mltiple moving objects in 2D incompressible viscos flow, J. Compt. Phys., 191 (2003), pp ] D. Salac and M. Miksis, A level set projection model of lipid vesicles in general flows, J. Compt. Phys., 230 (2011), pp ] David Salac and Michael J Miksis, Reynolds nmber effects on lipid vesicles, J. Flid Mech., 711 (2012), pp ] U Seifert, Configrations of flid membranes and vesicles, Adv. Phys., 46 (1997), pp ] Z Tan, D.V. Le, K.M Kim, and B.C. Khoo, An Immersed Interface Method for the Simlation of Inextensible Interfaces in Viscos Flids, Commn. Compt. Phys., 11 (2012), pp ] Z. Tan, D.V. Le, Z. Li, K.M. Lim, and B.C. Khoo, An immersed interface method for solving incompressible viscos flows with piecewise constant viscosity across a moving elastic membrane, J. Compt. Phys., 227 (2008), pp ] Z Tan, D Le, K Lim, and B Khoo, An Immersed Interface Method for the Incompressible Navier Stokes Eqations with Discontinos Viscosity Across the Interface, SIAM J. Sci. Compt., 31 (2009), pp ] SK Veerapaneni, D Geyffier, D Zorin, and G Biros, A bondary integral method for simlating the dynamics of inextensible vesicles sspended in a viscos flid in 2D, J. Compt. Phys., 228 (2009), pp ] SK Veerapaneni, A Rahimian, G Biros, and D Zorin, A fast algorithm for simlating vesicle flows in three dimensions, J. Compt. Phys., 230 (2011), pp ] Jian-Jn X and Hong-Kai Zhao, An Elerian Formlation for Solving Partial Differential Eqations Along a Moving Interface, J. Sci. Compt., 19 (2003), pp ] Sheng X and Z. Jane Wang, An immersed interface method for simlating the interaction of a flid with moving bondaries, J. Compt. Phys., 216 (2006), pp ] Sheng X and Z Jane Wang, Systematic derivation of jmp conditions for the immersed interface method in three-dimensional flow simlation, SIAM J. Sci. Compt., 27 (2006), pp