ANALYSIS OF A PHASE FIELD NAVIER-STOKES VESICLE-FLUID INTERACTION MODEL. Qiang Du, Manlin Li and Chun Liu
|
|
- Solomon Malone
- 5 years ago
- Views:
Transcription
1 DISCRETE AND CONTINUOUS Website: DYNAMICAL SYSTEMS SERIES B Volume 8, Number 3, October 27 pp ANALYSIS OF A PHASE FIELD NAVIER-STOKES VESICLE-FLUID INTERACTION MODEL Qiang Du, Manlin Li an Chun Liu Department of Mathematics Pennsylvania Sate University University Park, PA 1682, USA Abstract. This paper is concerne with the ynamics of vesicle membranes in incompressible viscous fluis. Some rigorous theory are presente for the phase fiel Navier-Stokes moel propose in [7], which is base on an energetic variation approach an incorporates the effect of bening elasticity energy for the vesicle membranes. The existence an uniqueness results of the global weak solutions are establishe. 1. Introuction. The stuy of hyroynamical an rheological properties of fluis involving vesicle membranes an cells is of interest in many biological an physiological applications. Consierable research efforts have been evote to both experimental stuies [1, 11, 13] an the evelopment of mathematical moels an computational coes of various egrees of physical relevance an sophistication regaring the membrane properties/configurations an the flui constitution in recent years [2, 3, 4, 7, 14, 15, 16, 2, 21, 22, 24, 25, 26, 27]. Vesicle membranes are forme by lipi bilayers which play an essential role in biological functions. Their equilibrium shapes are often characterize by minimizing the bening elastic energy of the membrane [12, 19, 23, 28]: E = Γ k 2 (H c ) 2 S (1) where Γ is the surface of vesicle membrane, H is the mean curvature of Γ, c the spontaneous curvature an k the bening moulus. It is known that the behavior of these vesicles, in both static configurations an uner external flow fiels, ramatically iffers from that of those roplets whose shape is governe by the surface tension (with surface energy epening only on the surface area of the membranes). In this paper, we continue our earlier stuies [7] an consier the phase fiel Navier-Stokes moel for the vesicle shape ynamics, which is governe by the coupling of the hyroynamic flui flow an the bening elastic properties of the vesicle membrane. The resulting membrane configuration an the flow fiel reflect the competition an the coupling of the kinetic energy an membrane elastic energies. 2 Mathematics Subject Classification. Primary: 58F15. Key wors an phrases. Phase Fiel, Navier-Stokes, Well-poseness, Weak Solution, Existence, Uniqueness, Vesicle Membrane, Flui Interaction, Bening Elastic Energy. The research of Q. Du an M. Li is supporte in part by the NSF DMS-49297, NSF ITR an NSF CCF The research of C. Liu is supporte in part by the NSF DMS an DMS
2 54 Q. DU, M. LI AND C. LIU In the phase fiel Navier-Stokes moel, the escription of the membrane is given in terms of a phase fiel function φ (see [8] an the subsequent works [5, 6, 7, 9, 1, 32] for etails). The phase fiel function φ, roughly speaking, is a labeling function efine on computational omain. The function φ takes value nearly +1 insie the vesicle membrane an 1 outsie, with a thin transition layer of with characterize by a small positive parameter ǫ. The zero level surface of φ represents the surface of vesicle membrane. The avantage of introucing such a labeling function is to formulate the original Lagrangian escription of the membrane evolution in the Eulerian (observer s) coorinates. As in [8], we will approximate the elastic bening energy (1) by E ǫ (φ) = k ( ǫ φ + ( 1 2 2ǫ ǫ φ + c 2)(1 φ )) 2 x. For illustration purposes, the fluis both insie an outsie the vesicle are taken to be an incompressible viscous Newtonian flui an the elastic energy associate with the vesicle eformation mainly comes from the bening energy. We want to point out that it is easy to incorporate other physical consierations of the fluis an the membrane into our energetic variational approach. The vesicle eformation an the flui velocity fiel are then regare as the result of the competition between vesicle membrane bening energy an flui kinetic energy, subject to the constraints that the volume an surface area of the vesicle are preserve. The equations governing the ynamics of the phase fiel function φ an the flui velocity fiel u can be obtaine via the energetic variation approach [7]. To enforce the two constraints, one may either aopt the Lagrange multiplier approach or use a penalty formulation [7]. Here, we focus on the latter, that is, we a two penalty terms to the elastic bening energy E ǫ (φ) to enforce the volume an surface area constraints respectively. As in [9] an [7], the moifie energy is given by where E(φ) = E ǫ (φ) M 1(A(φ) α) M 2(B(φ) β) 2, (2) A(φ) = φx, B(φ) = ǫ 2 φ ǫ (φ2 1) 2 x. We efine the following corresponing action functional which illustrate the competition between ifferent part of the energies. T 1 A[x(t, X)] = 2 x t(t, X) 2 x E(φ(t, x(t, X)))t (3) where x(t, X) can be thought as the incompressible flui trajectory in the Lagrangian coorinate an u being the flui velocity fiel. The Least Action principle yiels the actual force balance (linear momentum) equation [7]. We are thus le to the following phase fiel Navier-Stokes equation for φ an u: u t + u u = p + µ u + φ in [, T], iv u = in [, T], φ t + u φ = γ in [, T], (4) u(, x) = ũ(x) in, φ(, x) = φ(x) in
3 where PHASE FIELD NAVIER-STOKES EQUATIONS 541, the so-calle chemical potential, enotes the variational erivative of E(φ) in the variable φ (its precise form is escribe later). The above equation is complemente by bounary conitions (BC). The particular BC consiere in this paper is of the Dirichlet type for the phase fiel function φ an the no-slip bounary conition for the velocity fiel u: u = on, (5) φ = 1, φ = on. (6) The main objective of this paper is to provie a rigorous mathematical founation to the above couple phase fiel Navier-Stokes (PFNS) equation. In particular, we present the proof of existence an uniqueness of weak solution to (4). Our results inicate that we can essentially control the coupling between the velocity fiel an the phase fiel so that the natural (energy) solution spaces for the PFNS equations remain the same as that for the ecouple conventional incompressible Navier- Stokes equation an the simple phase fiel graient flow for the bening elastic energy. We elect to only focus on the case c = in this paper, though the proof can be reaily extene to the non-zero spontaneous curvature case. We note that some formal analysis has been given in [7] on the sharp interface limit of the couple PFNS moel as the interfacial with parameter ǫ. In particular, it is seen that uner a general ansatz assumption, the extra term φ in the momentum equation leas to the well-known Willmore force acting between the backgroun flui an the membrane surface (see also [3]), though the corresponing well-poseness results for the limiting system are still open an uner investigation. 2. Main results an formal estimates. In this section, we state our main results concerning the well-poseness of the couple PFNS moel an the properties of their weak solutions. Throughout the iscussion, we use the space H () to enote the space of ivergence free vector fiels in H 1(), an L2 () for the closure of ivergence free subset of Cc () in L2 (). H 1 () enotes the ual space of H (). For notational convenience, for any given time T, we also use spaces like L p (, T; L q ()) for functions of both the time an space variables as efine in [29]. In aition, we use <, > to enote the inner prouct in (an uality pairing with respect to) L 2 (), an we also use the following trilinear form: B(u, v, w) = u v w x. (7) The main results of this paper are the following existence an uniqueness theorems. Theorem 1. Existence of Weak Solution. Let be an open, boune subset of R 3 either having a smooth bounary or being a convex polyhera. There exists a pair of functions φ an u with 1. u L 2 (, T; H ()) W 1, 4 3 (, T; H 1 ()) 2. φ L 2 (, T; H 2 ()) H 1 (, T; L 2 ()) which is a weak solution to equation (4) with bounary conition (5-6), that is,
4 542 Q. DU, M. LI AND C. LIU 1. for any δ(x) H (), ξ(x) L 2 (), an a.e. t [, T], we have < u t, δ > +B(u, u, δ) = µ < u, δ > + φ δ x, (8) < φ t, ξ > +B(u, φ, ξ) = γ <, ξ >. 2. u(, x) = ũ(x), φ(, x) = φ(x) where ũ L 2 () an φ + 1 H 2 (). Theorem 2. Uniqueness of Weak Solution. For the weak solutions to equation (4) iscusse in the previous existence theorem, if in aition we have the solution satisfying u L 8 (, T; L 4 ()), then the weak solution is unique. A few remarks are first in orer. First, we aopt suitable assumptions on the omain so that we can obtain the H 2 regularity for the Laplace operator with homogeneous bounary conition. Secon, ue to the stanar theory for the conventional Navier-Stokes equations without the membrane stress [29] an the simple L 2 graient flow of the elastic bening energy without the flui transport [31], it is easy to see that the main task at han is to analyze the coupling terms in the PFNS equation, which has similar spirits as that in the stuy of couple systems for flui an liqui crystal irector [17]. Therefore, we nee to consier (control) the contribution to the momentum equation of the aitional stress tensor ue to the membrane eformation an the contribution of the convection term to the phase fiel evolution. With the energy law establishe below, it turns out that the solution space L 2 (, T; H ()) W 1, 4 3 (, T; H 1 ()) for the velocity fiel u remains the same as that for the conventional three imensional incompressible Navier- Stokes equations. This reflects the fact that the membrane stress tensor oes not pose any extra limitation on the regularity of the weak solution of the velocity fiel. Meanwhile, the solution space L 2 (, T; H 2 ()) H 1 (, T; L 2 ()) for φ also coincies with the natural space for the solution of the simple L 2 graient flow of the elastic bening energy, again showing that the effect of flui transport on the phase fiel function can also be properly controlle. Thir, ue to the limite regularity in u, the issue of uniqueness of the weak solution remains open even for the conventional Navier-Stokes equations in three imension without the membrane effect, thus, we o not have the proof of uniqueness of the global weak solutions for the couple PFNS equation in general. However, with a better regularity assumption on the weak solutions, as in the case of the conventional Navier-Stokes equations [29], the uniqueness can be assure. Before we turn to the proofs, let us also mention that we have chosen to work with the penalty formulation to incorporate the volume an surface area conservation of the vesicle membrane in time. The results presente here are for given penalty constants, so the constraints are satisfie only approximately. Careful examination of our proofs inicates that much of the estimates erive in the paper can be mae to be inepenent of the penalty constants, thus allowing one to extract a suitable limit as the penalty constants approach infinity. Such a limiting process woul lea the existence of solutions in the Lagrange multiplier formulation with the constraints exactly satisfie. However, a etaile account of the epenence an inepenence of the estimates on the penalty constants woul require teious bookkeeping notation-wise. To avoi complication, we elect to ignore such epenence in the later presentation. The etaile proofs of the above theorems are ivie into several parts which are given in the later sections. The main steps inclue establishing solution estimates
5 PHASE FIELD NAVIER-STOKES EQUATIONS 543 from the energy law an passing to the limit via a moifie Galerkin proceure. For a etaile account of similar techniques, we refer to the iscussion on the conventional incompressible Navier-Stokes equation given in [29]. We note that though we have use the bounary conition (5-6), much of our analysis can be carrie out in other cases as well, such as the Neumann an perioic bounary conitions, an in particular, the case of an inhomogeneous velocity profile on the bounary which is often use in the stuy of cell eformation in shear an/or extensional flows Formal erivation of the energy law. The issipation of the kinetic energy is one of the most basic property of the conventional incompressible Navier-Stokes equations, an it can be easily erive. It is thus interesting to note that a similar energy law hols for the couple PFNS equation, with the membrane bening elastic energy being ae to the kinetic energy to prouce the total energy. For convenience, let us enote f(φ) = ǫ φ + 1 ǫ (φ2 1)φ, g(φ) = f(φ) + 1 ǫ 2 (3φ2 1)f(φ). Then, we may rewrite the energy as E(φ) = k f(φ) 2 x + 1 2ǫ 2 M 1(A(φ) α) M 2(B(φ) β) 2. The irect computation shows = kg(φ) + M 1(A(φ) α) + M 2 (B(φ) β)f(φ). (9) We now give a formal erivation of the energy law for smooth solutions u an φ of (4). Multiply u to the first equation in (4) an to the secon equation, then integrate over, we get the following issipative energy law 1 t 2 u 2 x + E(φ) = µ u 2 x γ 2 x. (1) Immeiately one can conclue that if u an φ are the solutions of the PFNS, we have the uniform bouns (with respect to any T > ) of the following type: u L (, T; L 2 ()) L 2 (, T; H ()); φ L (, T; H 2 ()); f(φ) L (, T; L 2 ()); L 2 ((, T) ()). For weak solutions of (4), we can erive a weaker version of the energy law rigorously via a Galerkin proceure outline below Formal estimates of u t H 1 () an φ t L2 (). Base on the bouns from the energy law an the PFNS equation, we may euce better estimates on the time erivatives. 1. Let v L 2 () an v L 2 () 1, from equation (4), φ t v x u φv x + γ v x.
6 544 Q. DU, M. LI AND C. LIU An, u φv x C v L 2 () φ L3 () u L6 () C φ H2 () u H 1 (), Therefore, φ t L 2 () C vx L 2 () v L 2 (). ( φ H 2 () u H 1 () + ) L 2 (). 2. We now take v H (). Also from equation (4), we get u t v x u u v x + φ v x +µ u L2 () v L2 (), u t H 1 () u u H 1 () + µ u L 2 () + φ H 1 (). By a well known interpolation result [29], we have, u u H 1 () C u 1 2 L 2 () u 3 2 L 2 (). As for φ H 1 (), we consier for any v H 1 (), with v H 1 () 1. By the Sobolev inequality v L 6 () C v H 1 (), we get φ v x L 2 () φ L 3 () v L 6 () K 1 L 2 () φ L 3 () K 2 L 2 () φ H 2 (). Hence, ( u t H 1 () C u 3 2 L2 () u 1 2 L2 () + u L 2 () + ) L 2 () φ H 2 (). 3. Proof the existence of weak solutions. This section is evote to the proof of existence of weak solutions of the couple PFNS equation (4) with bounary conition (5-6). We first outline a moifie Galerkin approximation, then we consier its weak limit an verify it as a weak solution of (4).
7 PHASE FIELD NAVIER-STOKES EQUATIONS A moifie Galerkin approximation. Let us choose {ω n } L 2 () to be the eigenfunctions of Stokes operator, such that {ω n } forms an orthonormal basis for L 2 (). Set W n = Span{ω 1, ω 2,..., ω n }. Apply the Galerkin approximation to the velocity fiel u, one can get the approximate equation for u W n an φ H 2 (): u t + P n (u u) = µ u + P n ( φ), φ t + u φ = γ, (11) u() = P n ũ(x), φ() = φ(x) where P n is the L 2 projection operator to W n. The following Lemma asserts the existence of a solution to the approximate equation (11). It also provies a uniform energy estimate on the solution (with respect to the imension n). Lemma 1. Existence of An Approximate Solution There exists a pair of functions u(t, x) W n, φ(t, x) H 2 () satisfying < u t, w > +B(u, u, w) = µ < u, w > + φ w x, < φ t, v > +B(u, φ, v) = γ <, v >, (12) u() = P n ũ(x) φ() = φ(x), w W n an v L 2 () for almost all t [, T]. Furthermore, for a.e. ˆT [, T], ˆT µ u 2 + γ 2 xt + with a constant M inepenent of W n. 1 2 u( ˆT, x) 2 x + E(φ( ˆT, x)) M, Proof. 1. Apply Galerkin approximation to φ Denote by {υ 1, υ 2,...} the eigenfunctions of operator uner the homogeneous Dirichlet bounary conition. They consist of an orthonormal basis of L 2 (). By the assumption on the omain, we have that the eigenfunctions of the Laplace operator have H 2 regularity. Set V m = Span{υ 1, υ 2,...,υ m }. Apply a moifie Galerkin metho to φ, we get an approximate equation to equation (12) as follows: fin u m (x, t) an φ m (x, t) of the form n m u m (x, t) = i (t)ω i (x) W n an φ m (x, t) + 1 = h j (t)υ j (x) V m i=1 such that for k = 1, 2,...,n an l = 1, 2,..., m, < u m, ω k > +B(u m, u m, ω k ) = µ < u m, w k > + π m ( δe(φm) ) φ m ω k x, δe(φ m) υ l x, < φ m, υ l > +B(u m, φ m, υ l ) = γ j=1 u m (, x) = P n ũ(x), φ m (, x) = π m ( φ(x) + 1) 1. (13) Here means ifferentiating in time, an π m enotes a L 2 (also H 2 if applicable) projection to V m. The solutions have natural epenence on the inex n, but for convenience, we have suppresse this epenence in the notation of u m an φ m.
8 546 Q. DU, M. LI AND C. LIU It is easy to see that the above finite imensional ODE system has a solution local in time. 2. Energy estimate In equation (13), replace w k with u m, υ l with π m ( δe(φm) ), we have, 1 u m 2 x = µ < u m, u m > + π m ( δe(φ m) ) φ m u m x t 2 < φ m, π m( δe(φ m) ) > + π m ( δe(φ m) ) φ m u m x Since = γ δe(φ m ) π m ( δe(φ m) )x. < φ m, π m( δe(φ m) ) >=< π m (φ m ), δe(φ m) >=< φ m m, δe(φ m) >= m t E(φ m), the summation of the two expressions above gives the energy law, t 1 2 It implies, u m 2 x + E(φ m ) = µ 1 2 u m 2 x γ ˆT u m ( ˆT, x) 2 x + E(φ m ( ˆT, x)) + π m( δe(φ m) 2 ) x. (14) µ u m 2 +γ π m( δe(φ 2 m) ) xt 1 u m (, x) 2 x + E(φ m (, x)). 2 We know that u m (, x) L 2 () ũ(x) L 2 (). By the construction of φ m (, x), we also have that φ m (, x) converges to φ(x) in H 2 () as m. Hence there exists some constant M inepenent of W n, such that, 1 2 ˆT + ˆT u m ( ˆT, x) 2 x + E(φ m ( ˆT, x)) + γ π m( δe(φ m) 2 ) xt M. µ u m 2 xt Note that such an energy law essentially gives the existence of local solutions for all time. 3. Compactness of {u m } an {φ m } The energy law ensures that u m (t) L 2 () an φ m (t) L 2 () are uniformly boune in m an t [, T]. Thus the solution of ODE (13) actually exists global in time. Furthermore, the energy law also inicates u m is uniformly boune in L (, T; L 2 ()) L 2 (, T; H ()). φ m is uniformly boune in L (, T; H 2 ()) (thus in L 2 (, T; H 2 ()))
9 PHASE FIELD NAVIER-STOKES EQUATIONS 547 π m ( δe(φm) ) is uniformly boune in L 2 ((, T) ). Similar to the previous formal erivation, we can rigorously obtaining estimates on u m H 1 () an φ m L 2 () : u m is uniformly boune L 4 3(, T; H 1 ()). φ m is uniformly boune L 2 ((, T) ). Therefore, using the Aubin-Lions type compact embeing results [29], there exist some φ, u, such that, u m has a subsequence u mk converging to u weakly in L 2 (, T; H ()) an strongly in L 2 (, T; L 2 ()). φ m has a subsequence φ ml converging to φ weakly in L 2 (, T; H 2 ()) an strongly in L 2 (, T; W 1,p ()) for 1 p < 6. For convenience, if there is no ambiguity, from now on we will ientify u m an φ m with their subsequences. 4. Passing weak limits of {u m } an {φ m } Choose w(t, x) = α(t)δ(x), v(t, x) = α(t)ξ(x) where α C([, T]), δ W n, an ξ V m, we have T < u m, w > +B(u m, u m, w)t = T In which, T T < φ m, v > +B(u m, φ m, v) = γ < π m ( δe(φ m) ) φ m, w > t T δe(φ m ) = k{ f(φ m ) + 1 ǫ 2 (3φ2 m 1)f(φ m )} µ < u m, w > t < δe(φ m), v > t. +M 1 (A(φ m ) α) + M 2 (B(φ m ) β)f(φ m ) = k{ [ ǫ φ m + 1 ǫ (φ2 m 1)φ m] + 1 ǫ 2 (3φ2 m 1)f(φ m)} +M 1 (A(φ m ) α) + M 2 (B(φ m ) β)f(φ m ) (15) = ǫk 2 φ m + L(φ m ) where L(φ m ) enotes the lower orer term. Note that the L(φ m ) L2 ((,T) ) is uniformly boune by the uniform boun on φ m in L (, T; H 2 ()). Then, π m ( δe(φ m) ) = π m (ǫk 2 φ m ) + π m (L(φ m )) = ǫk 2 φ m + π m (L(φ m )). Together with the energy estimate (14), we have ǫk 2 φ m L 2 ((,T) ) π m ( δe(φ m) ) L2 ((,T) ) + π m (L(φ m )) L2 ((,T) ) π m ( δe(φ m) ) L 2 ((,T) ) + L(φ m ) L 2 ((,T) ).
10 548 Q. DU, M. LI AND C. LIU Therefore, δe(φ m) L2 ((,T) ) π m ( δe(φ m) ) L2 ((,T) ) + 2 L(φ m ) L2 ((,T) ). (a) Weak limit of { δe(φm) }. We now prove that (a subsequence of) { δe(φm) } converges to weakly in L 2 ((, T) ). By the energy law (14), f(φ m ) is uniformly boune in L (, T; L 2 ()), hence uniformly boune in L 2 ((, T) ). It is enough to show for g C ([, T] ), T lim m T f(φ m )g xt = f(φ m )g xt = T T f(φ)g xt. [ ǫ φ m + 1 ǫ (φ2 m 1)φ m]g xt. It is sufficient to check on the nonlinear terms only. Since φ m is uniformly boune in L (, T; H 2 ()), we have φ m (t) C, 1 2 () C φ m(t) H2 () M for t [, T]. Hence, φ m (t, x) is uniformly boune in [, T]. Furthermore, φ m has as subsequence converging to φ strongly in L 2 (, T, W 1,p ) for 1 p < 6. Then a subsequence of φ m converges to φ almost everywhere in [, T]. By the Lebesgue-Dominate Theorem, lim n T φ 3 mg xt = Similarly, we nee to verify that T δe(φ m ) lim g xt = m T T T φ 3 g xt. g xt. (16) To give more etails, we follow (15) term by term. First, we have T T 2 φ m g xt = φ m g xt φ g xt = T 2 φg xt as m. Next, we have T (φ 3 m φ 3 )g xt = 3 T (φ 2 m φ m φ 2 φ) g xt T C (φ 2 m φ 2 ) φ m g xt + C T φ 2 ( φ φ m ) g xt C φ 2 m φ 2 L 2 ((,T) ) φ m L 2 ((,T) ) + C φ φ m L 1 ((,T) ).
11 PHASE FIELD NAVIER-STOKES EQUATIONS 549 Hence, as m, T (φ 3 m φ3 )g xt. T Now, consier (φ 2 mf(φ m ) φ 2 f(φ))g xt T (f(φ m ) f(φ))φ 2 g xt T + (φ 2 φ 2 m )gf(φ m)xt = I 1 + I 2. We have I 1 since f(φ m ) f(φ) weakly in L 2 ((, T) ). In aition, I 2 g L ((,T) ) f(φ m ) L 2 ((,T) ) φ 2 φ 2 m L 2 ((,T) ). It is also easy to show, T lim m B(φ m )f(φ m )g xt = T B(φ)f(φ)g xt. (b) Verifying the approximate equation. Choose g L 2 ((, T) ). Since δe(φm) weakly converges to π m (g) converges strongly to g in L 2 ((, T) ), we have T T = T, < π m ( δe(φ m) ), g > t, g >)t < δe(φ m), π m (g) > < < δe(φ m), g > t + T < δe(φ m), π m (g) g > t an as m, We can conclue that π m ( δe(φm) in L 2 ((, T) ). Now, one can let m to recover ) converges to weakly T T < u t, w > +B(u, u, w)t = T + < φ t, v > +B(u, φ, v) = γ T T µ < u, w > t φ w xt, <, v > t.
12 55 Q. DU, M. LI AND C. LIU An, ˆT ( 1 u( 2 x) 2 x + E(φ( ˆT, x)) + µ u 2 + γ xt M. Since α = α(t) is arbitrarily chosen in C([, T]), one can conclue for any δ W n an ξ V m, < u t, δ > +B(u, u, δ) = µ < u, δ > + φ δ x, < φ t, ξ > +B(u, φ, ξ) = γ ξ x. (17) By ensity argument, it is also true for any δ W n an ξ L 2 (). Setting α() = 1 an α(t) =, one can show that u(, x) = P n ũ(x) an φ(, x) = φ(x). Finally, since for i > n T u m (t, x)ζ(t)ω i (x)xt = for any ζ(t) C([, T]), we have by taking m that u(t, x)ω i (x)x = for almost all t [, T] when i > n. Therefore u W n. This completes the proof of this lemma Proof of the existence theorem. We now wrap up the proof of the existence theorem. Accoring to Lemma 1, for any positive time ˆT (, T), an for each W n, the equation (11) has a solution u n an φ n, such that ˆT ( µ u n 2 + γ δe(φ ) n) 2 xt+ n where M is inepenent of W n. Hence, Also, 1 2 u n( ˆT, x) 2 x+e(φ n ( ˆT, x)) M (18) u n is uniformly boune in L (, T; L 2 ()) L 2 (, T; H ()); φ n is uniformly boune in L (, T; H 2 ()) ( thus in L 2 (, T; H 2 ()) ); is uniformly boune in L 2 ((, T) ). δe(φn) n u n is uniformly boune L 3(, 4 T; H 1 ()); φ n is uniformly boune L2 (, T; L 2 ()). Therefore, there exist some φ an u, such that, u n has a subsequence u nk converging to u weakly in L 2 (, T; H ()) an strongly in L 2 ((, T) ); φ n has a subsequence φ nl converging to φ weakly in L 2 (, T; H 2 ()) an strongly in L 2 (, T; W 1,p ()) for 1 p < 6.
13 PHASE FIELD NAVIER-STOKES EQUATIONS 551 Similar to the claim in Lemma 1, (a subsequence of ) δe(φn) n converges weakly to in L 2 ((, T) ). Choose w(t, x) = α(t)δ(x), v(t, x) = α(t)ξ(x) where α C([, T]), δ W n, an ξ C(), then consier T < u n, w > +B(u n, u n, w)t = Let n, we get, T T < u t, w > +B(u, u, w)t = T T < δe(φ n) φ n, w > t n < φ n, v > +B(u n, φ n, v) = γ T T < T φ, w > t < φ t, v > +B(u n, φ, v)t = γ T T µ < u n, w > t <, v > t. µ < u, w > t, <, v > t. Because α(t) is an arbitrarily chosen function in C([, T]), one can conclue for any δ W n, ξ C(), < u t, δ > +B(u, u, δ) = µ < u, δ > + < φ t, ξ > +B(u, φ, ξ) = γ <, ξ > φ δ x, By a ensity argument, it is also true for any δ H 1 () an ξ L2 (). Set α() = 1 an α(t) =, one can show u(, x) = ũ(x), an φ(, x) = φ(x). This conclues the proof of the main existence theorem. 4. Uniqueness of weak solution. We now provie the proof of the uniqueness of the weak solution, uner the aitional regularity assumption on the velocity fiel. We first introuce a few notations to simplify the later iscussion. Define G(φ) = 1 (kǫ φ 2 + kǫ ) 2 φ 2 + φ 2 x, then it is easy to check that 1 C φ 2 H 2 () G(φ) C φ 2 H 2 (). Define also M(φ) = δg(φ) = kǫ 2 φ k φ + φ, N(φ) = ǫ M(φ). Assume u i, φ i an p i (i = 1, 2) are two weak solutions to equation (4) satisfying the assumptions given in the uniqueness theorem. Let û = u 1 u 2, ˆφ = φ 1 φ 2 an ˆp = p 1 p 2. We erive a Gronwall type inequality for û an ˆφ to prove the uniqueness Derivation of a Gronwall inequality. First, we have, û + û u 1 + u 2 û = ˆp + µ û + (M(φ 1 ) + N(φ 1 )) φ 1 (M(φ 2 ) + N(φ 2 )) φ 2, ˆφ + u 1 φ 1 u 2 φ 2 = γ(m(ˆφ) + N(φ 1 ) N(φ 2 )). (19)
14 552 Q. DU, M. LI AND C. LIU Multiply û to the first equation in (19) an M(ˆφ) to the secon one, integrate in space, we get 1 t 2 û 2 x + û u 1 ûx = µ û 2 x + (M(φ 1 ) φ 1 M(φ 2 ) φ 2 )û + (N(φ 1 ) φ 1 N(φ 2 ) φ 2 )ûx, t G(ˆφ) + (u 1 φ 1 u 2 φ 2 )M(ˆφ)x = γ M(ˆφ) 2 x Hence, γ M(ˆφ)(N(φ 1 ) N(φ 2 ))x 1 t 2 û 2 x + G(ˆφ) + γ M(ˆφ) 2 x + µ û 2 x + B(û, u 1, û) = (M(φ 1 ) φ 1 M(φ 2 ) φ 2 )û (u 1 φ 1 u 2 φ 2 )M(ˆφ)x + (N(φ 1 ) φ 1 N(φ 2 ) φ 2 )ûx γ Recall B(u, v, w) = u v wx. Since B(û, u 1, û) = B(û, û, u 1 ), we have B(û, u 1, û) = B(û, û, u 1 ) û L4 () û L2 () u 1 L4 () û 1/4 L 2 () û 7/4 L 2 () u 1 L 4 () M(ˆφ)(N(φ 1 ) N(φ 2 ))x û 2 L 2 () + C( ) û 2 L 2 () u 1 8 L 4 () where is an arbitrary small positive number. Direct calculation shows, (M(φ 1 ) φ 1 M(φ 2 ) φ 2 )û (u 1 φ 1 u 2 φ 2 )M(ˆφ)x an = (M(φ 1 ) ˆφ + M(ˆφ) φ 2 )û (u 1 ˆφ + û φ 2 )M(ˆφ)x = M(φ 1 ) ˆφ ûx u 1 ˆφ M(ˆφ)x = J 1 J 2, J 1 M(φ 1 ) L2 () ˆφ L4 () û L4 () û 2 H 1 () + C M(φ 1) 2 L 2 () ˆφ 2 H 2 (),
15 PHASE FIELD NAVIER-STOKES EQUATIONS 553 J 2 M(ˆφ) L 2 () u 1 L 4 () ˆφ L 4 () M(ˆφ) 2 L 2 () + C u 1 2 H 1 () ˆφ 2 H 2 (). with Moreover, (N(φ 1 ) φ 1 N(φ 2 ) φ 2 )ûx = (N(φ 1 ) ˆφ ûx + = K 1 + K 2, (N(φ 1 ) N(φ 2 )) φ 2 ûx K 1 N(φ 1 ) L 2 () ˆφ L 4 () û L 4 () C N(φ 1 ) L 2 () ˆφ H 2 () û H 1 () an û 2 C H 1 () + N(φ 1) 2 L 2 () ˆφ 2 H 2 (), K 2 N(φ 1 ) N(φ 2 ) L2 () φ 2 L4 () û L4 () Similarly, M(ˆφ)(N(φ 1 ) N(φ 2 ))x C N(φ 1 ) N(φ 2 ) L2 () φ 2 H2 () û H 1 () û 2 C H 1 () + N(φ 1) N(φ 2 ) 2 L 2 () φ 2 2 H 2 () û 2 C H 1 () + N(φ 1) N(φ 2 ) 2 L 2 (). We now use a claim (to be verifie later): M(ˆφ) L 2 () N(φ 1 ) N(φ 2 ) L 2 () M(ˆφ) 2 L 2 () N(φ 1) N(φ 2 ) 2 L 2 (). N(φ 1 ) N(φ 2 ) L 2 () C φ 1 φ 2 H 2 () = C ˆφ H 2 (). (2) Recall that ˆφ 2 H 2 () CG(ˆφ). Using (2) an putting everything together, we have 1 t 2 û 2 x + G(ˆφ) C 1 2 û 2 x + G(ˆφ) ) ( M(φ 1 ) 2 L 2 () + u 1 2H 1() + u 1 8 L 4 () + N(φ 1) 2 L 2 () + 1. Using the estimates alreay erive an the extra assumption on the velocity fiel, we have ( M(φ 1 ) 2 L 2 () + u 1 2 H 1() + u 1 8 L 4 () + N(φ 1) 2 L 2 () )
16 554 Q. DU, M. LI AND C. LIU is integrable in time. This implies if û(t) = an ˆφ(t) = at t =, then û = an ˆφ = for all time, which proves the uniqueness theorem Proof of claim (2). We now verify the claim (2) use in the above proof. By (9) an the efinitions of M an N, we get N(φ) = k ǫ φ3 + 2k 3k φ + ǫ ǫ 2 φ2 f(φ) k ǫ 2 f(φ) φ +M 1 (A(φ) α) + M 2 (B(φ) β)f(φ). This leas to N(φ 1 ) N(φ 2 ) L 2 () C{ φ 3 1 φ3 2 L 2 () + φ 1 φ 2 L 2 () + ˆφ L2 () + φ 2 1 f(φ 1) φ 2 2 f(φ 2) L2 () + f(φ 1 ) f(φ 2 ) L2 () + A(ˆφ) L 2 () + B(φ 1 )f(φ 1 ) B(φ 2 )f(φ 2 ) L 2 ()} Accoring to the energy estimate erive in the proof of existence theorem, we have φ i L (, T; H 2 ()), i = 1, 2. Therefore, for i, j = 1, 2, φ i (t, x) L ((,T) ) M, (φ i (t, x)φ j (t, x)) L ((,T);L 6 ()) M, (φ i (t, x)φ j (t, x)) L ((,T);L 2 ()) M for some constant M. Now we carefully estimate the iniviual terms respectively. A generic time-inepenent constant C is use. Let φ = φ φ 1φ 2 + φ 2 2, then φ 3 1 φ3 2 L 2 () ˆφ φ L 2 () + ˆφ φ L 2 () +2 ˆφ φ L2 () C ˆφ H2 (). f(φ 1 ) f(φ 2 ) L 2 () C( ˆφ L 2 () + ˆφ L 2 () + φ 3 1 φ3 2 L 2 ()) C ˆφ H2 (). φ 2 1 f(φ 1) φ 2 2 f(φ 2) L2 () (φ 2 1 φ2 2 )f(φ 1) L2 () + φ 2 2 (f(φ 1) f(φ 2 )) L2 () C ˆφ H2 (). B(φ 1 ) B(φ 2 ) C( (φ 1 + φ 2 ) L2 ()) ˆφ L2 () +C (φ φ2 2 2)(φ 1 + φ 2 )ˆφx C ˆφ H 2 (). B(φ 1 )f(φ 1 ) B(φ 2 )f(φ 2 ) L2 () B(φ 1 )(f(φ 1 ) f(φ 2 )) L2 () + (B(φ 1 ) B(φ 2 ))f(φ 2 ) L 2 () C ˆφ H 2 (). Summing together, the claim (2) is verifie.
17 PHASE FIELD NAVIER-STOKES EQUATIONS Conclusion. In this paper, a system of couple phase fiel Navier-Stokes equations moeling the eformation an evolution of three imensional vesicle membranes in a flui fiel is analyze. A major characteristic of the current moel stuie here, in comparison of many other similar moels stuie previously in the literature on membrane flui interactions, even in the phase fiel context, is the inclusion of the bening elastic energy which is crucial for vesicle bilayers. The resulting interaction mechanism is thus base on the competition between the elastic bening energy of the membrane, with prescribe bulk volume an surface area, an the kinetic energy in the surrouning flui velocity fiels. The variation of the elastic bening energy leas to an extra stress in the Navier-Stokes equation, which involves a nonlinear combination of higher orer spatial erivatives of the phase fiel function. Our main results illustrate that these aitional contributions can be properly controlle ue to the establishment of the energy law. Such an energy issipation mechanism is intrinsic in our erivation of the couple system. The results provie a rigorous mathematical founation to the couple phase-fiel Navier-Stokes equations an their numerical simulations [7]. The results can be extene to cases with the Neumann an perioic bounary conitions, as well as inhomogeneous Dirichlet conition for the velocity fiel. We note that our analysis largely relies on the amping term in the evolution of the phase fiel function an can not be reaily extene to the case of a pure transport. Such cases were treate in [18] for the systems for viscoelastic materials. However, the nonlinear coupling terms involve higher erivatives in this paper. In aition, the various estimates erive in the paper are not uniform with respect to the small interfacial with parameter ǫ an thus cannot be use in the stuy of the sharp interface limits. Such issues pose interesting challenges for future stuies. Extensions to interactions of vesicles with other types of fluis (possibly with ifferent types insie an outsie the vesicle) may also be consiere. Finally, it will be also interesting to stuy the existence of classical solutions of the system in this paper, in particular, the existence for the large viscosity situations. REFERENCES [1] M. Abkarian, C. Lartigue, an A. Viallat, Tank Treaing an Unbining of Deformable Vesicles in Shear Flow: Determination of the Lift Force, Phys. Rev. Lett., 88 (22), 6813 [2] J. Beaucourt, F. Rioual, T. Sion, T. Biben, an C. Misbah, Steay to unsteay ynamics of a vesicle in a flow, Phys. Rev. E, 69 (24), 1196 [3] T. Biben, K. Kassner an C. Misbah, Phase-fiel approach to three-imensional vesicle ynamics, Phys. Rev. E, 72 (25), [4] T. Biben an C. Misbah, Tumbling of vesicles uner shear flow within an avecte-fiel approach, Phys. Rev. E, 67 (23), 3198 [5] Q. Du, C. Liu, R. Ryham an X. Wang, A phase fiel formulation of the Willmore problem, Nonlinearity, 18 (25), , [6] Q. Du, C. Liu, R. Ryham an X. Wang, Moeling the Spontaneous Curvature Effects in Static Cell Membrane Deformations by a Phase Fiel Formulation, Communications in Pure an Applie Analysis, 4 (25), [7] Q. Du, C. Liu, R. Ryham an X. Wang, Moeling Vesicle Deformations in Flow Fiels via Energetic Variational Approaches, preprint, 26. [8] Q. Du, C. Liu, an X. Wang, A phase fiel approach in the numerical stuy of the elastic bening energy for vesicle membranes, Journal of Computational Physics, 198 (24), [9] Q. Du, C. Liu, an X. Wang, Retrieving Topological Information For Phase Fiel Moels, SIAM Journal on Applie Mathematics, 65 (25),
18 556 Q. DU, M. LI AND C. LIU [1] Q. Du, C. Liu, an X. Wang, Simulating the Deformation of Vesicle Membranes uner Elastic Bening Energy in Three Dimensions, Journal of Computational Physics, 212 (26), [11] K. e Haas, C. Bloom, D. van en Ene, M. Duits, an J. Mellema, Deformation of giant lipi bilayer vesicles in shear flow, Phys. Rev. E, 56 (1997), [12] W. Helfrich, Elastic properties of lipi bilayers: theory an possible experiments. Z. Naturforsch. C, 28 (1973), [13] V. Kantsler an V. Steinberg, Orientation an Dynamics of a Vesicle in Tank-Treaing Motion in Shear Flow, Phys. Rev. Lett., 95 (25), [14] S. Keller an R. Skalak, Motion of a tank-treaing ellipsoial particle in a shear flow, J. Flui Mech., 12 (1982), [15] M. Kraus, W. Wintz, U. Seifert, an R. Lipowsky. Flui Vesicles in Shear Flow, Phys. Rev. Lett., 77 (1996), [16] H. Li, H. Yi, X. Shan an H. Fang, Shape changes an motion of a vesicle in a flui using a lattice Boltzmann moel, Arxiv preprint physics/6774 [17] F. Lin an C. Liu, Nonparabolic issipative systems moeling the flow of liqui crystals, Comm. Pure. Appl. Math., XLVIII (1995), [18] F. Lin, C. Liu an P. Zhang, On Hyroynamics of Viscoelastic Fluis, Comm. Pure Appl. Math., Vol LVIII (25), [19] R. Lipowsky The morphology of lipi membranes, Current Opinion in Structural Biology, 5 (1995), [2] Y. Liu, L. Zhang, X. Wang an W.K. Liu, Coupling of Navier-Stokes equations with protein molecular ynamics an its application to hemoynamics, Inter. J. Numerical Methos in Fluis, 46 (24), [21] H. Noguchi an G. Gompper, Flui Vesicles with Viscous Membranes in Shear Flow, Phys. Rev. Lett., 93 (24), [22] H. Noguchi an G. Gompper, Meshless membrane moel base on the moving least-squares metho, Phys. Rev. E, 73 (26), 2193 [23] Z. Ou-Yang, J. Liu, an Y. Xie, Geometric Methos in the Elastic Theory of Membranes in Liqui Crystal Phases, Worl Scientific, Singapore, [24] C. Pozrikiis, Interfacial ynamics for Stokes flow, Journal of Computational Physics, 169 (21), [25] C. Pozrikiis, Numerical Simulation of the Flow-Inuce Deformation of Re Bloo Cells, Annals of Biomeical Engineering, 31 (23), [26] U. Seifert, Flui membranes in hyroynamic flow fiels: Formalism an an application to fluctuating quasispherical vesicles in shear flow, European Physical Journal B, 8 (1999), [27] U. Seifert, Hyroynamic Lift on Boun Vesicles, Phys. Rev. Lett., 83 (1999), [28] U. Seifert, K. Bernl an R. Lipowsky, Configurations of flui membranes an Vesicles, Physical Rev A, 44 (1991), , [29] R. Temam. Navier-Stokes Equations, Theory an Numerical Analysis, American Mathematical Society, 21 [3] X. Wang, Optimal Profile of the phase fiel moels for Willmore problems, preprint, 26. [31] X. Wang an Q. Du, Convergence of numerical approximations to a phase fiel bening elasticity moel of membrane eformations, to appear in Inter. J. Numer. Anal an Moeling, 26. [32] X. Wang an Q. Du, Moelling an Simulations of Multi-component Lipi Membranes an Open Membranes via a Diffuse Interface Approach, to appear in J. Mathematical Biology, 26. Receive August 26; revise February 27. aress: qu@math.psu.eu aress: li m@math.psu.eu aress: liu@math.psu.eu
The total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationGLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS
Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationAbstract A nonlinear partial differential equation of the following form is considered:
M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationOn Kelvin-Voigt model and its generalizations
Nečas Center for Mathematical Moeling On Kelvin-Voigt moel an its generalizations M. Bulíček, J. Málek an K. R. Rajagopal Preprint no. 1-11 Research Team 1 Mathematical Institute of the Charles University
More informationWUCHEN LI AND STANLEY OSHER
CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationProblem set 2: Solutions Math 207B, Winter 2016
Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationSINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationA transmission problem for the Timoshenko system
Volume 6, N., pp. 5 34, 7 Copyright 7 SBMAC ISSN -85 www.scielo.br/cam A transmission problem for the Timoshenko system C.A. RAPOSO, W.D. BASTOS an M.L. SANTOS 3 Department of Mathematics, UFSJ, Praça
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS. Gianluca Crippa
Manuscript submitte to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS Gianluca Crippa Departement Mathematik
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationOn a class of nonlinear viscoelastic Kirchhoff plates: well-posedness and general decay rates
On a class of nonlinear viscoelastic Kirchhoff plates: well-poseness an general ecay rates M. A. Jorge Silva Department of Mathematics, State University of Lonrina, 8657-97 Lonrina, PR, Brazil. J. E. Muñoz
More informationGlobal Solutions to the Coupled Chemotaxis-Fluid Equations
Global Solutions to the Couple Chemotaxis-Flui Equations Renjun Duan Johann Raon Institute for Computational an Applie Mathematics Austrian Acaemy of Sciences Altenbergerstrasse 69, A-44 Linz, Austria
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More information2.20 Marine Hydrodynamics Lecture 3
2.20 Marine Hyroynamics, Fall 2018 Lecture 3 Copyright c 2018 MIT - Department of Mechanical Engineering, All rights reserve. 1.7 Stress Tensor 2.20 Marine Hyroynamics Lecture 3 1.7.1 Stress Tensor τ ij
More informationWELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES
Electronic Journal of Differential Equations, Vol. 017 (017), No. 38, pp. 1 7. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationA new proof of the sharpness of the phase transition for Bernoulli percolation on Z d
A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition
More informationEuler Equations: derivation, basic invariants and formulae
Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the
More informationFLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction
FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More informationON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS
ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS ADRIAN CONSTANTIN AND BORIS KOLEV Abstract. Accoring to the principle of least action, the spatially perioic motions of one-imensional
More informationCentrum voor Wiskunde en Informatica
Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationIntroduction to variational calculus: Lecture notes 1
October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational
More informationA nonlinear inverse problem of the Korteweg-de Vries equation
Bull. Math. Sci. https://oi.org/0.007/s3373-08-025- A nonlinear inverse problem of the Korteweg-e Vries equation Shengqi Lu Miaochao Chen 2 Qilin Liu 3 Receive: 0 March 207 / Revise: 30 April 208 / Accepte:
More informationThe Generalized Incompressible Navier-Stokes Equations in Besov Spaces
Dynamics of PDE, Vol1, No4, 381-400, 2004 The Generalize Incompressible Navier-Stokes Equations in Besov Spaces Jiahong Wu Communicate by Charles Li, receive July 21, 2004 Abstract This paper is concerne
More informationarxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003
Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationTopological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem
Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationConvective heat transfer
CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More information3 The variational formulation of elliptic PDEs
Chapter 3 The variational formulation of elliptic PDEs We now begin the theoretical stuy of elliptic partial ifferential equations an bounary value problems. We will focus on one approach, which is calle
More informationarxiv: v1 [cond-mat.stat-mech] 9 Jan 2012
arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationCOUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND
More informationChapter 9 Method of Weighted Residuals
Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More informationSuppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane
Arch. Rational Mech. Anal. Digital Object Ientifier (DOI) https://oi.org/1.17/s25-18-1336-z Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane Siming He & Eitan Tamor Communicate
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationExponential asymptotic property of a parallel repairable system with warm standby under common-cause failure
J. Math. Anal. Appl. 341 (28) 457 466 www.elsevier.com/locate/jmaa Exponential asymptotic property of a parallel repairable system with warm stanby uner common-cause failure Zifei Shen, Xiaoxiao Hu, Weifeng
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More information7.1 Support Vector Machine
67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to
More informationTractability results for weighted Banach spaces of smooth functions
Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS
ON ISENTROPIC APPROXIMATIONS FOR COMPRESSILE EULER EQUATIONS JUNXIONG JIA AND RONGHUA PAN Abstract. In this paper, we first generalize the classical results on Cauchy problem for positive symmetric quasilinear
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationWELL-POSTEDNESS OF ORDINARY DIFFERENTIAL EQUATION ASSOCIATED WITH WEAKLY DIFFERENTIABLE VECTOR FIELDS
WELL-POSTEDNESS OF ORDINARY DIFFERENTIAL EQUATION ASSOCIATED WITH WEAKLY DIFFERENTIABLE VECTOR FIELDS Abstract. In these short notes, I extract the essential techniques in the famous paper [1] to show
More informationMany problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx
Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus
More informationThe continuity equation
Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationOn some parabolic systems arising from a nuclear reactor model
On some parabolic systems arising from a nuclear reactor moel Kosuke Kita Grauate School of Avance Science an Engineering, Wasea University Introuction NR We stuy the following initial-bounary value problem
More informationA Spectral Method for the Biharmonic Equation
A Spectral Metho for the Biharmonic Equation Kenall Atkinson, Davi Chien, an Olaf Hansen Abstract Let Ω be an open, simply connecte, an boune region in Ê,, with a smooth bounary Ω that is homeomorphic
More informationRegularity of Attractor for 3D Derivative Ginzburg-Landau Equation
Dynamics of PDE, Vol.11, No.1, 89-108, 2014 Regularity of Attractor for 3D Derivative Ginzburg-Lanau Equation Shujuan Lü an Zhaosheng Feng Communicate by Y. Charles Li, receive November 26, 2013. Abstract.
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationarxiv: v2 [math.dg] 16 Dec 2014
A ONOTONICITY FORULA AND TYPE-II SINGULARITIES FOR THE EAN CURVATURE FLOW arxiv:1312.4775v2 [math.dg] 16 Dec 2014 YONGBING ZHANG Abstract. In this paper, we introuce a monotonicity formula for the mean
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationMartin Luther Universität Halle Wittenberg Institut für Mathematik
Martin Luther Universität alle Wittenberg Institut für Mathematik Weak solutions of abstract evolutionary integro-ifferential equations in ilbert spaces Rico Zacher Report No. 1 28 Eitors: Professors of
More informationStable Poiseuille Flow Transfer for a Navier-Stokes System
Proceeings of the 26 American Control Conference Minneapolis, Minnesota, USA, June 14-16, 26 WeB2.2 Stable Poiseuille Flow Transfer for a Navier-Stokes System Rafael Vázquez, Emmanuel Trélat an Jean-Michel
More information6 Wave equation in spherical polar coordinates
6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.
More informationDelocalization of boundary states in disordered topological insulators
Journal of Physics A: Mathematical an Theoretical J. Phys. A: Math. Theor. 48 (05) FT0 (pp) oi:0.088/75-83/48//ft0 Fast Track Communication Delocalization of bounary states in isorere topological insulators
More informationBOUNDEDNESS IN A THREE-DIMENSIONAL ATTRACTION-REPULSION CHEMOTAXIS SYSTEM WITH NONLINEAR DIFFUSION AND LOGISTIC SOURCE
Electronic Journal of Differential Equations, Vol. 016 (016, No. 176, pp. 1 1. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu BOUNDEDNESS IN A THREE-DIMENSIONAL ATTRACTION-REPULSION
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationContinuum Mechanics Lecture 4 Fluid dynamics
Continuum Mechanics Lecture 4 Flui ynamics Prof. http://www.itp.uzh.ch/~teyssier Outline - Flui kinematics - Mass an momentum conservation laws - The energy equation - Real fluis - Ieal fluis - Incompressible
More information