OPTIMAL CONTROL OF DISTRIBUTED PARAMETER SYSTEM GIVEN BY CAHN-HILLIARD EQUATION
|
|
- Barbara Gordon
- 5 years ago
- Views:
Transcription
1 Nonlinear Functional Analysis and Applications ol. 19, No. 1 (214), pp Copyright c 214 Kyungnam University Press KUPress OPTIMAL CONTROL OF DISTRIBUTED PARAMETER SYSTEM GIEN BY CAHN-HILLIARD EQUATION Quan-Fang Wang 1 and Shin-ichi Nakagiri 2 1 Chinese University of Hong Kong N. T. Hong Kong wang quanfang@yahoo.co.jp 2 Kobe University Kobe, , Japan nakagiri@kobe-u.ac.jp Abstract. In this work, optimal control problem is addressed for distributed systems described by Cahn-Hilliard equation by the means of distributed control, initial control and Neumann boundary control. The existence and uniqueness is provided for weak solution using variational method. Further, existence of optimal control is proved completely, and optimality conditions is established for integral cost and quadratic cost, respectively. Lastly, Bang-Bang principle is deduced. 1. Introduction In investigating optimal control problem for distributed parameter system described by Cahn-Hilliard (C-H) equation, let be an open bounded set R n of x with a piecewise smooth boundary Γ =. Let Q = (, T ) and Σ = (, T ) Γ with T >. Consider the system described by Cahn-Hilliard equation as y t + γ 2 y λ f(y) = B v in Q, y =, ( y) = B 1 v 1 on Σ, (1.1) y() = Bw on, Received July 4, 213. Revised January 7, Mathematics Subject Classification: 35B37, 35A4, 35A15, 35K35, 35K55, 35Q93, 49J2, 49J4, 93C2. Keywords: Cahn-Hilliard equation, optimal control, necessary optimality condition.
2 2 Q. F. Wang and S. Nakagiri where γ, λ > are constants. denote Laplacian in R n. y is the concentration function of binary components. B, B 1, B are control operators, v, v 1, w are control variables, is a unit outer normal derivative. f(y) is a polynomial of order 2p 1 as follows f(y) = 2p 1 j=1 a j y j, p N, p 1. (1.2) Here N is set of natural numbers, and the leading coefficient of f is positive, i.e., a 2p 1 >. Cahn-Hilliard equation is based on a continuum model for phase transition in binary systems such as alloy, glasses, and polymer-mixtures (cf. [3]). Since Cahn and Hilliard proposed the equation in 1958, contributing research had been reported in many related works (cf. [5, 7, 15]). Paper [4] using free energy method for C-H equation. Article [8] show finite difference scheme for C-H equation with numerical experiment. Further, [11] is a paper on discontinuous finite element method for C-H equation with demonstrations using color graphics visualization. The discontinuous Galerkin method is executed in [22] for Cahn-Hilliard equation to get good experimental agreement. All of these works are focused either on system (analytic, numerical) solution or on system properties. There are few results (cf. [23]) on optimal control problem theoretically and computationally, even nonlinear term in [23] is limited to 3 order polynomial. To develop optimal control theory as [12] for system (1.1) with integral cost in the form of J(v) = Φ(y(T, v)) + L(t, v, y(t, v))dt, v U, (1.3) where Φ, L are continuous functions on t [, T ], v = (v, v 1, w) is a control vector and U is a Hilbert space of v, U = 1 W and, 1, W are Hilbert spaces of control variables v, v 1, w. Our goal is to find such optimal control v and characterize optimality conditions on v such that inf J(v) = J(v ) for an admissible set U ad U. v U ad The paper is organized as follows. Section 2 gives the mathematical setting and some preparation. Section 3 states existence and uniqueness of weak solution for free system. Section 4 considers the optimal control problem of system (1.1) with integral performance index (1.3). Section 5 contains conclusions.
3 Optimal control of distributed parameter C-H equation Mathematical setting If using [6, 17, 18, 19, 2, 21] to give mathematical setting. Consider the free system with B v = g, B 1 v 1 = k 1, Bw = y in (1.1). Introduce Hilbert space H = L 2 () as usual. Define Hilbert space = H 2 ( ; ) = { φ φ L 2 (), φ L 2 (), with inner product (φ, ϕ) = φ φ(x)ϕ(x)dx + φ(x) ϕ(x)dx. = on Γ} is equipped with the norm φ = φ L 2 () + φ L 2 (). Neglecting φ C R, it is clear that the norm φ is equivalent to φ L 2 () and φ H 2 () in /R (cf. [19, 2, 21]). Let be the dual space of, denote its norm by. The symbol, denotes the dual pair of and. Denote φ = dφ dt to define the solution space by W (, T ;, ) = { φ φ L 2 (, T ; ), φ L 2 (, T ; ) }. According to Neumann boundary condition in (1.1), Hilbert space H 3 2 (Γ) is introduced with norm H 3 2 (Γ) defined by ( φ 3 H 2 = (Γ) Γ φ(s) 2 ds + Γ α =1 D α φ(s) 2 φ(s) φ( s) 2 ds + Γ Γ s s n+1 dsd s Its dual space H 3 2 (Γ) with norm denoted by H 2 3 [1], the symbol (Γ), H 3 2 (Γ),H 3 denotes the dual pair of 2 (Γ) H 3 2 (Γ) and H 3 2 (Γ). Definition 2.1. A function y is a weak solution (cf. [21]) of (1.1) if y W (, T ;, ) satisfies weak form d (y, v) + γ( y, v) λ(f(y), v) dt = g, v + γk 1, v Γ H 3 2 (Γ),H 2 3, v, (2.1) (Γ) y() = y H, in the sense of D (, T ), where D (, T ) denotes the space of distributions on (, T ). In order to obtain weaker restrict on nonlinear term as in [21], the assumption on p and n is needed to assure f(φ) L 2 () for all φ. Define ) 1 2.
4 22 Q. F. Wang and S. Nakagiri the arbitrary p, if 1 n 4. [ n 2 ] Assumption (A) : p <, if 5 n 6. n 4 p = 1, if n 7. Here [ ] denote the Gauss symbol. Citing Gagliardo-Nirenberg inequality in [9, 14, 15]. Lemma 2.2. Let p, r, q be integers. If 1 r < q and p q. Then for φ(t) W m,p () L r () to have φ(t) C φ(t) W k,q () holds with some C > and θ = ( k n + 1 r 1 q θ W m,p () φ(t) 1 θ L r () ) ( m n ) 1 r p (2.2) provided that < θ 1 (assume that < θ < 1 if q = + ). Here in (2.2), k, m are the derivatives orders of φ(t) in space W k,q (), W m,p (), respectively. k < m. Lemma 2.3. Let assumption (A) be satisfied, then f(φ(t)) H C ( φ(t) 4p 2 n(p 1) 2 H ) φ(t) n(p 1) for φ(t) /R, where C > is constant independent of φ(t). Taking r = p = m = 2, q = 4p 2, k = in (2.2) of Lemma 2.2, then n(p 1) θ = < 1 to verify Lemma p 2 Lemma 2.4. Let assumption (A) be satisfied, then f(y(t)), as a nonlinear mapping: H, satisfies the local Lipschitz continuity. That is for some C >, and φ, ϕ to have f(φ(t)) f(ψ(t)) 2 H C( φ(t) 4p 4 One can refer (6.17) in p. 175 of [16] for above Lemma. + ψ(t) 4p 4 ) φ(t) ψ(t) 2.
5 Optimal control of distributed parameter C-H equation Existences and uniqueness of weak solution Consider y L 2 (), y(t) for t [, T ], then theoretical study is closely related to n and p of the nonlinearity (cf. [21]). If using Gauss symbol [ ], define [ p n = ], n N. n Consider y L 2 () and y(t) for t [η, T ], where η > arbitrary, then the following theorem depends on the conditions of n and p (cf. [21]). One can propose arbitrary [ p, if 1 n 2. n ] Assumption (B) : p, if 3 n 4. n 2 p = 1, if n 5. Theorem 3.1. Let assumption (A) be satisfied. Assume that y L 2 () and g L 2 (, T ; ), then there exists a unique weak solution y W (, T ;, ) for free system of (1.1), which is belonging to For C > to get estimate L (, T ; H) L 2 (, T ; ), T >. (3.1) y 2 L (,T ;H) + y 2 L 2 (,T ; ) ( C 1+ y 2 H+ g 2 L 2 (,T ; ) + ) k 1 2 L 2 (,T ;H 3. (3.2) 2 (Γ)) Further, if p p n, then we obtain estimation y 4p 2 C L 4p 2 (,T ;L 4p 2 ()) ( 1+ y 2 H + g 2 L 2 (,T ; ) + ) k 1 2 k, L 2 (,T ;H 2 3 (3.3) (Γ)) where C and k are constants independent of y, g and k 1. In addition, for any < η < T, if n and p satisfy assumption (B), then the estimation can be given by y 4p 2 L 4p 2 (η,t ;L 4p 2 ()) ( C(η) 1+ y 2 H + g 2 L 2 (,T ; ) + ) k 1 2 k, L 2 (,T ;H 2 3 (3.4) (Γ)) where C(η) and k are independent of y, g and k 1. Proof. Since the boundary Γ of is piecewise smooth, by trace theorem (cf. [13]), it is easy to verify (3.1) and (3.2) (cf. g = in [16]). To prove (3.2) and
6 24 Q. F. Wang and S. Nakagiri (3.3), from Lemma 2.3 to deduce that y(t) 4p 2 dx C ( 1 + y 2 H + g 2 L 2 (,T ; ) + ) k 1 2 k L 2 (,T ;H 2 3 y(t) n(p 1) (Γ)), (3.5) where k = 4p 2 n(p 1). If p p n and p 1, using Hölder inequality to get that ( y n(p 1) T ) n(p 1) dt y 2 2 dt T 2 n(p 1) 2. (3.6) By (3.2), (3.5) and (3.6) to deduce the estimate (3.3). To derive (3.4), let F denote the primitive of f vanishing at y =, i.e. F (y) = jb j = a j 1 and 2 j 2p. Define Lyapunov function by Y (y(t)) = γ 2 y(t) 2 H + λ F (y(t))dx. Let us quote a result at footnote of p.155 in [16], for arbitrary η >, for t [η, T ] to have Y (y(t)) Y (y(η)) y L (η, T ; H 1 ()) L (η, T ; L 2p ()). Further, if y H 1 () L 2p (), then Y (y(t)) Y (y ) for all t. Setting r = 2p, p = m = 2, q = 4p 2, k = in Lemma 2.2, then take n(p 1) θ 1 = 1 under assumption (B). Applying Gagliardo- (2p 1)(4p n(p 1)) Nirenberg inequality (cf. [9, 14]) to y(t) H 2 () L 2p (), by the equivalence of norm y(t) H 2 (), y(t) and y(t) H in /R to deduce that 2p j=2 b j y j y(t) L 4p 2 () C y(t) θ 1 y(t) 1 θ 1 L 2p () < (3.7) for y(t) L 2p () and some C >. Therefore, by (3.7) to find that η y(t) 4p 2 L 4p 2 () dt C sup y(t) 1 θ 1 L 2p () t [,T ] under assumption (B). It provides (3.4). existence by [6] in a routing way. η y(t) θ 1(4p 2) dt η y(t) 2 dt for Above a priori estimate to prove The regularity results can be proved by deriving an energy-type equation.
7 Optimal control of distributed parameter C-H equation 25 In order to prove uniqueness for y L 2 (), we use the following inequalities. Let f(s) be a polynomial defined by (1.2). Then for b, c R, there exists C such that f(b) f(c) C b c (1 + b 2p 2 + c 2p 2 ). (3.8) By Young inequality (cf. [12]), for a, b >, it is clear that ab εp p ap + 1 q b q ε q, where 1 p + 1 q = 1. (3.9) Let y 1 (t) and y 2 (t) be two weak solutions of (1.1) on [, T ] with respect to inputs g(t), k 1 (t) and y. Set ȳ(t) = y 1 (t) y 2 (t), then ȳ(t) satisfies dȳ(t) dt ȳ() =, + γ 2 ȳ(t) λ (f(y 1 (t)) f(y 2 (t))) =, (3.1) in weak sense. Consider y 1 (t), y 2 (t), ȳ(t), and eak form of (2.1), then we have 1 d 2 dt ȳ(t) 2 H + γ ȳ(t) 2 H = λ(f(y 1 (t)) f(y 2 (t)), ȳ(t)). (3.11) Let us utilize (3.8) to estimate the right hand side of (3.11) as follows (f(y 1 (t)) f(y 2 (t)), ȳ(t)) L 2 () C ε ȳ(t) 2 L 2 () + Cε ȳ(t) 2 H + C +C ȳ(t) ȳ(t) y 2 (t) 2p 2 dx. ȳ(t) ȳ(t) y 1 (t) 2p 2 dx For the exponent of y i (i = 1, 2), using Hölder inequality, taking p = 2(4p 2) 4p 2 + n(p 1), q = such that 1 p + 1 q = 1 in (3.9), then deduce that ȳ(t) ȳ(t) y i (t) 2p 2 dx C εp p ȳ(t) 2 + 2(4p 2) 4p 2 n(p 1) C yi q (t) (2p 2)(8p 4) 4p 2 n(p 1) ȳ(t) 2 ε q L 4p 2 () H (3.12)
8 26 Q. F. Wang and S. Nakagiri for C >. Setting k i (t) = C q y i (t) (2p 2)(8p 4) 4p 2 n(p 1) and C ε q L 4p 2 () i (ε) = C εp p, then obtain ȳ(t) ȳ(t) y i (t) 2p 2 dx C i (ε) ȳ(t) 2 + k i (t) ȳ(t) 2 H. If y L 2 (), one can prove that k i (t) is integrable on [, T ]. In fact y i (t) (2p 2)(8p 4) 4p 2 n(p 1) dt C y i 2p 2 L (,T ;H) L 4p 2 () y i (t) 2n(p 1) 2 4p 2 n(p 1) 4p 2 n(p 1) for C >. Since n(p 1) 2 > 1 as p p n, by Hölder inequality and for some k >, one has that yi (t) 2n(p 1) 2 4p 2 n(p 1) ( k dt C y i (t) 2 dt) as p pn. This implies that k i (t) is integrable on [, T ] as p p n for y L 2 (). By (3.12), one can rewrite (3.11) as follows 1 d 2 dt ȳ(t) 2 H + (γ Cε CλC(ε)) ȳ(t) 2 H ( 1 ) C ε + λk(t) ȳ(t) 2 H, (3.13) where C(ε) = C 1 (ε) + C 2 (ε) and K(t) = k 1 (t) + k 2 (t). Finding ε in (3.13) such that γ Cε CλC(ε) >, by Bellman-Gronwall theorem (cf. [16]) to have ȳ(t) H = in H for all t [, T ] as p p n. It completes the proof of uniqueness in Theorem 3.1. dt 4. Optimal control problems As pointed out early, U = 1 W denotes Hilbert space of control vector v = (v, v 1, w). Let B L(, L 2 (, T ; )), B 1 L( 1, L 2 (, T ; H 3 2 (Γ))), B L(W, H). For control vector v = (v, v 1, w), by virtue of Theorem 3.1, the solution mapping v y(v) from U into W (, T ;, ) is well defined. Here y(v) is called the state of control system (1.1). The integral cost associated with the control system (1.1) is given by (1.3), and we regard Φ : H R +, L(t, v, y(t, v)) : [, T ] U H R +. Let U ad = ad ad 1 W ad be a closed convex (bounded) subset of U = 1 W, which is called the admissible set. Integral cost optimal control problem subject to (1.1) and (1.3) is:
9 Optimal control of distributed parameter C-H equation 27 (i) find an element v = (v, v 1, w ) U ad such that (ii) characterize such v. inf J(v) = J(v ); v U ad Such a v is called optimal control for C-H problem (cf. [12]) Existence of optimal control. Assume U ad is bounded. Suppose assumptions on Φ and L: H(Φ) Function Φ : H R + is continuous and convex. H(L) Function L : [, T ] U H R + is an integrand such that (a) For arbitrary (v, y) U H, L(t, v, y) is measurable in t [, T ]; (b) For a.e. t [, T ], L(t, v, y) is continuous and convex for each (v, y) U H; (c) For arbitrary bounded set E U, there exists an m E L 1 (, T ) s.t. sup L(t, v, y(t, v)) m E (t), a.e. in [, T ]; (4.1) v E (d) L(t, v, y) is locally uniformly Lipschitz continuous with respect to y, i.e., for bounded set K = E F U H, there exists an m K L 2 (, T ) s.t. L(t, φ, y) L(t, φ, z) m K (t) y z H. (4.2) Theorem 4.1. Let H(Φ) and H(L) be satisfied and assume y L 2 (). Consider y(t), t [, T ] as p p n and y(t), t [η, T ] under assumption (B). If U ad is closed convex (bounded) subset, then there exists at least one optimal control v U ad such that v minimizes the cost (1.3). Proof. By virtue of Theorem 3.1, there exists a weak solution y(v) of equation (1.1) with B v L 2 (, T ; ), B 1 v 1 L 2 (, T ; H 3 2 (Γ)) and Bw H. y(v) is uniformly bounded for v E, i.e., sup{ y(t, v) H : v E, t [, T ]} < +. By (c) of H(L) to know that J(v) is bounded on E, it means that J(v) makes sense for any v U. Since U ad is non-empty, there exists a sequence {v n }, v n = (v n, vn 1, wn ), in U such that inf J(v)= lim J(v n) = J. Here y n is the v U ad n trajectory corresponding to v n, that is y n = y(t, v n ). Because U ad is bounded, convexity and closed, one can choose a subsequence {v m }, v m = (v m, vm 1, wm ), of {v n } and find a v = (v, v 1, w ) U ad such that v m v weakly in U, as m. (4.3)
10 28 Q. F. Wang and S. Nakagiri Analogous to the proof of Theorem 2.1, there exists a subsequence (rewritten as {y(v m )}) of {y(v m )} and z W (, T ;, ) such that y(v m ) z weakly in L (, T ; H), (4.4) y(v m ) z weakly in W (, T ;, ), (4.5) y(v m ) z strongly in L 2 (, T ; H). (4.6) Because of convergences (4.4), (4.5) and Gagliardo-Nirenberg inequality (cf. [9, 14]), one sees that y(v m ) belongs to L 4p 2 (, T ; L 4p 2 ()) for any m as p p n. This implies that f(y(v m )) belongs to L 2 (Q) as p p n. By (4.6) and the continuity of f(y) with respect to y, we know that f(y(v m )(t, x)) f(z(t, x)) a.e. in Q as p p n. Noticing the boundedness of f(y(v m )) for any m, by the same argument in subsection 4.2 to deduce that f(y(v m )(t, x)) f(z(t, x)) weakly in L 2 (Q) (4.7) as m. The weak form is expressed as d dt (y(t, v m), v) + γ( y(t, v m ), v) λ f(y(t, v m )) vdx = B v m (t), v + γb 1 v m 1 (t), v Γ H 3 2 (Γ),H 3 2 (Γ), (4.8) for v. By (4.3), (4.4), (4.5), (4.6) and (4.7), taking the limit in (4.8) yields (2.1). One can obtain z = y(v ) via uniqueness, that is, function y(v ) is the state corresponding to control variable v. Therefore from (4.5) to have y(t, v m ) y(t, v ) weakly in H, a.e. in [, T ]. (4.9) It is well known that continuity plus convexity imply weak lower semi-continuity (cf. [12]), then from H(Φ) and (4.9) with t = T to get lim Φ(y(T, v m )) Φ(y(T, v )). (4.1) m Since v m v weakly in U, there exists a constant M > such that v m U M, this means there exists a bounded set E = {v : v m U M} U. Then by the estimate in Theorem 3.1, we have { } N =sup y(t, v m ) H : v m U M, t [, T ] <+. This deduces that there exists a bounded set F = {y : y(t, v m ) H N} H. For E and F, there exists a bounded set K = E F = { v : v m U M } { y : y(t, v m ) H N } U H. (4.11)
11 Optimal control of distributed parameter C-H equation 29 By local uniformly Lipschitz continuity assumption (d) with (4.2) of H(L), for the bounded set K there exists an m K L 2 (, T ) such that lim m L(t, v m, y(t, v m )) L(t, v m, y(t, v )) dt m K L 2 (,T ) lim m y(t, v m) y(t, v ) L 2 (,T ;H)=. (4.12) By assumption (c) with (4.1) of H(L), for the bounded set E, there exists m E L 1 (, T ) s.t. sup L(t, v m, y(t, v )) dt m E (t)dt < +. m Also because continuity plus convexity imply weak lower semi-continuity, the assumption (b) of H(L) and (4.3), it follows from the Lebesgue-Fatou lemma (cf. [13]) that lim m L(t, v m, y(t, v ))dt L(t, v, y(t, v ))dt. (4.13) Therefore, from (4.1), (4.12) and (4.13) to get J = lim inf J(v m) J(v ), m then J(v ) = inf J(v), i.e., v is an optimal control for cost (1.3). It finish v U ad proof of Theorem Necessary conditions for optimality. It is well known that the optimality condition in [12] for v is given by the variational inequality J (v )(v v ), v U ad, where J (v ) denotes the Gâteaux derivative of J(v) in (1.3) at v. In order to derive the optimality conditions, the following assumptions are posed. A(Φ) The function Φ : H R + is Gâteaux differentiable, and Φ (y) is continuous on H. A(L) The function L : [, T ] U H R + satisfies (a) For fixed (t, v) [, T ] U, there exists Gâteaux derivative L y(t, v, y), which is continuous in (v, y) U H. For bounded set K = E F U H, there exists m 1 k (t) L1 (, T ) s.t. sup L y(t, v, y) L(H) m 1 k (t) a.e. in [, T ]; (v,y) K (b) For fixed (t, y) [, T ] H, there exists Gâteaux derivative L v(t, v, y), which is continuous in v U. For bounded set K = E F U H, there exists an m 2 k (t) L2 (, T ) s.t. sup L v(t, v, y) L(U) m 2 k (t) a.e. in [, T ]. (v,y) K
12 3 Q. F. Wang and S. Nakagiri The Gâteaux derivative of the cost J(v) can be calculated by assumption A(Φ) and A(L). It means that the cost J(v) is weak Gâteaux differentiable at v in the direction v v. Therefore optimality condition J (v )(v v ) can be rewritten as Φ (y(t, v ))z(t ) + L y(t, v, y(t, v ))z(t)dt + L v (t, v, y(t, v ))(v v )dt (4.14) for all v U ad. Here z(t) = Dy(v )(v v ) is Gâteaux differential of y(v) at v in the direction v v. In general, solve problem (ii) by introducing adjoint system, forming (4.14) to derive theorem on optimality condition. Theorem 4.2. Let H(Φ), H(L), A(Φ) and A(L) be satisfied, and assume that y L 2 (). If supposing y(t), t [, T ] as p p n, or y(t), t [η, T ] under assumption (B). Then the optimal control v = (v, v 1, w ) U ad for (1.3) is characterized by optimality system y t + γ 2 y λ f(y) = B v in Q, y =, ( y) = B 1 v1 on Σ, y() = Bw on. p(v ) + γ 2 p(v ) λf t y(y(v )) p(v ) = L y(t, v, y(v )) in Q, p(v ) ( p(v )) =, = on Σ, p(t, v ) = Φ (y(t, v )) on. B (v v ), p(v ) dt + +(B(w w ), p(, v )) H + for all v = (v, v 1, w) U ad. γb 1 (v 1 v 1), p(v ) Γ L v(t, v, y(v ))(v v )dt, H 3 2 (Γ),H 3 2 (Γ) dt Consider the distributed observation z 1 (v) = C 1 y(v) and terminal value observation z 2 (v) = C 2 y(t, v), where C 1 L(W (, T ;, ), M 1 ) and C 2 L(H, M 2 ) are operators, also called the observers, and M 1, M 2 are observation
13 Optimal control of distributed parameter C-H equation 31 spaces. The cost function associated with (1.1) is given by J(v) = C 1 y(v) z d 2 M 1 + C 2 y(t, v) z T d 2 M 2 + (Nv, v) U, v U. (4.15) Here z d M 1, z T d M 2 are desired values of z 1 (v) and z 2 (v) at t and T, respectively. N = (N, N 1, N) L(U, U) is symmetric and positive operator, (Nv, v) U = (N v, v ) + (N 1 v 1, v 1 ) 1 + (Nw, w) W. Theorem 4.3. Let the assumption in Theorem 4.1 be satisfied. Then the optimal control v = (v, v 1, w ) U ad for cost function (4.15) is characterized by the following equalities and inequality, optimality system y t + γ 2 y λ f(y) = B v in Q, y =, ( y) = B 1 v1 on Σ, y() = Bw on, p(v ) + γ 2 p(v ) λf (y(v )) p(v ) t = C1 Λ M 1 (C 1 y(v ) z d ) in Q, p(v ) =, ( p(v )) = on Σ, p(t, v ) = C2 Λ M 2 (C 2 y(t, v ) zd T ) on. + for all v = (v, v 1, w) U ad. B (v v ), p(v ) dt γb 1 (v 1 v 1), p(v ) Γ H 3 2 (Γ),H 3 2 (Γ) dt +(B(w w ), p(, v )) H + (Nv, v v ) U, (4.16) If specifying that U = 1 W = L 2 (Q) L 2 (Σ) L 2 (), and U ad = ad 1 ad W ad = {v v a v v b, a.e. on Q} {v 1 v a 1 v 1 v b 1, a.e. on Σ} {w w a w w b, a.e. on } with v a, vb L (Q), v a 1, vb 1 L (Σ) and w a, w b L (). Assume that N =, B = B = B 1 = I and U ad is nonempty. Since U ad is closed and convex
14 32 Q. F. Wang and S. Nakagiri in U, then from the optimality condition (4.16), we have v (t) v (t), p(t, v ) dt + (w w, p(, v )) H + γ(v1 (t) v 1(t)), p(t, v ) Γ H 3 2 (Γ),H 3 2 (Γ) dt (4.17) for all v U ad, where v = (v, v 1, w ) U ad. By setting (v, v 1, w) U ad in (4.17), we get (p(, v ), w w ) L 2 (), w W ad. By Lebesgue convergence theorem in [12] and (4.17), we have for a.e. t [, T ], (p(t, v ), v (t) v (t)) L 2 (), v ad, (p(t, v ) Γ, γ(v 1 (t) v 1(t))) L 2 (Γ), v 1 1 ad. Then one can convert the following property of v : i) w (x) = w a (x) if p(, v, x) > ; w (x) = w b (x) if p(, v, x) < for x. ii) v (t, x) = va (t, x) if p(t, v, x) > ; v (t, x) = vb (t, x) if p(t, v, x) < for (t, x) Q. iii) v1 (t, ξ) = va 1 (t, ξ) if p(t, v, ξ)> ; v1 (t, ξ) = vb 1 (t, ξ) if p(t, v, ξ) < for (t, ξ) Σ. As is well known, it is Bang-Bang principle (cf. [12]) of optimal control v = (v, v 1, w ). 5. Conclusions In this paper, surveyed optimal control problem of distributed parameter system given by Cahn-Hilliard equation by the means of distributed control, initial control and Neumann boundary control. Further, for integral and quadratic cost function, existences of optimal control is straitforwardly proved under rational assumptions, and necessary optimality conditions are considered, respectively. Acknowledgments The authors would like to express full of gratitude to Japan Scientific Research Grant (C) (2) (No ) from year 1999 to 22. The first author thank to Dr. C. Cao and E. Bernard (University of Connecticut). References [1] R. Adams, Sobolev Spaces, Academic Press, New York (1975).
15 Optimal control of distributed parameter C-H equation 33 [2] H.T. Banks, Control and Estimation in Distributed Parameter Systems, Society for Industrial and Applied Mathematics (1992). [3] J.W. Cahn and J.E. Hilliard, Free energy of a nonuniform system-i: Interfacial free energy, The Journal of Chemical Physics, 28 (1958), [4] A. Cohen, Energy method for the Cahn-Hilliard equations, Quart. Appl. Math., (1988), [5] M. Copetti and C. Elliott, Numberical analysis of the Cahn-Hilliard equation with a logarithmic free energy, Numerische Mathematik, 63 (1992), [6] R. Dautary and J.L. Lions, Mathematical Analysis and Numerical Methods for Science and Technology, Springer-erlag, Berlin-Heidelberg-New York (1992). [7] C. Elliot and S. Zheng, On the Cahn-Hilliard equations, Arch. Rat. Mech. Anal., 96 (1986), [8] D. Furihata, A stable and conservative finite difference scheme for the Cahn-Hilliard equation, Numerische Mathematik, 87 (21), [9] E. Gagliardo, Proprietà di alcune classi di funzioni più variabili, Ricerche Mat., 7 (1958), [1] I. Gelfand and G. Silov, Functions généralisees 1,2,3: Moscow (1958). [11] D. Kay,. Styles and E. Süli, Discontinuous Galerkin finite element approximation of the Cahn-Hilliard equation with convection, SIAM Journal on Numerical Analysis, 47 (29) [12] J.L. Lions, Optimal Control of Systems Governed by Partial Differential Equations, Springer-erlag, Berlin Heidelberg New York (1971). [13] J.L. Lions and E. Magenes, Non-Homogeneous Boundary alue Problems and Applications I, II, Springer-erlag, Berlin-Heidelberg-New York (1972). [14] L. Nirenberg, On ellliptic partial differential equations, Ann. Sc. Norm. Pisa, 13 (1959), [15] J.Y. Park and J.J. Bae, On the existence of solution of strongly damped nonlinear wave equations, Internat. J. Math. and Math. Sci., 23 (2), [16] R. Temam, Infinite-Dimensional Dynamical Systems in Mechanics and Physis, Applied Mathematica Sciences 68, Springer-erlag, Berlin-Heidelberg-New York (1988). [17] Q.F. Wang, Practical Application of Optimal Control Theory: computational approach: Lambert Academic Publishing, Germany (211). [18] Q.F. Wang, Optimal Control for Nonlinear Parabolic Distributed Parameter Systems: with numerical analysis: Lambert Academic Publishing, Germany (211). [19] Q.F. Wang, Optimal distributed control of nonlinear Cahn-Hilliard systems with computational realization, Journal of Mathematical Sciences, 177 (211), [2] Q.F. Wang and S. Nakagiri, Weak solution of Cahn-Hilliard equations having forcing term and optimal control problems, Research Institute for Mathematical Sciences, Kyoto University, Japan, Kokyuroku, 1128 (2), [21] Q.F. Wang, Parameters identification problem of distributed system given by Cahn- Hilliard equations, International Journal of Pure and Applied Mathematics, 22 (25), [22] G. Well, E. Kuhl and K. Garikipati, A discontinuous Galerkin method for the Cahn- Hilliard equation, Journal of Computational Physics, 218 (26), [23] J. Yong and S. Zheng, Feedback stabilization and optimal control for Cahn-Hilliard equation, Nonlinear Analysis TMA, 17 (1991),
Junhong Ha and Shin-ichi Nakagiri
J. Korean Math. Soc. 39 (22), No. 4, pp. 59 524 IDENTIFICATION PROBLEMS OF DAMPED SINE-GORDON EQUATIONS WITH CONSTANT PARAMETERS Junhong Ha and Shin-ichi Nakagiri Abstract. We study the problems of identification
More informationA CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS 1. INTRODUCTION
Trends in Mathematics Information Center for Mathematical Sciences Volume 4, Number 2, December 21, Pages 39 44 A CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS
More informationThe Navier-Stokes Equations with Time Delay. Werner Varnhorn. Faculty of Mathematics University of Kassel, Germany
The Navier-Stokes Equations with Time Delay Werner Varnhorn Faculty of Mathematics University of Kassel, Germany AMS: 35 (A 35, D 5, K 55, Q 1), 65 M 1, 76 D 5 Abstract In the present paper we use a time
More informationWeak Solutions to Nonlinear Parabolic Problems with Variable Exponent
International Journal of Mathematical Analysis Vol. 1, 216, no. 12, 553-564 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ijma.216.6223 Weak Solutions to Nonlinear Parabolic Problems with Variable
More informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationWELL-POSEDNESS OF WEAK SOLUTIONS TO ELECTRORHEOLOGICAL FLUID EQUATIONS WITH DEGENERACY ON THE BOUNDARY
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 13, pp. 1 15. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu WELL-POSEDNESS OF WEAK SOLUTIONS TO ELECTRORHEOLOGICAL
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationSobolev Spaces. Chapter Hölder spaces
Chapter 2 Sobolev Spaces Sobolev spaces turn out often to be the proper setting in which to apply ideas of functional analysis to get information concerning partial differential equations. Here, we collect
More informationELLIPTIC EQUATIONS WITH MEASURE DATA IN ORLICZ SPACES
Electronic Journal of Differential Equations, Vol. 2008(2008), No. 76, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ELLIPTIC
More informationANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE INTO CARTOON PLUS TEXTURE. C.M. Elliott and S.A.
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 6, Number 4, December 27 pp. 917 936 ANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE
More informationNonlinear Analysis 72 (2010) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 72 (2010) 4298 4303 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na Local C 1 ()-minimizers versus local W 1,p ()-minimizers
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 155 A posteriori error estimates for stationary slow flows of power-law fluids Michael Bildhauer,
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
More informationEXISTENCE AND REGULARITY RESULTS FOR SOME NONLINEAR PARABOLIC EQUATIONS
EXISTECE AD REGULARITY RESULTS FOR SOME OLIEAR PARABOLIC EUATIOS Lucio BOCCARDO 1 Andrea DALL AGLIO 2 Thierry GALLOUËT3 Luigi ORSIA 1 Abstract We prove summability results for the solutions of nonlinear
More informationON NONHOMOGENEOUS BIHARMONIC EQUATIONS INVOLVING CRITICAL SOBOLEV EXPONENT
PORTUGALIAE MATHEMATICA Vol. 56 Fasc. 3 1999 ON NONHOMOGENEOUS BIHARMONIC EQUATIONS INVOLVING CRITICAL SOBOLEV EXPONENT M. Guedda Abstract: In this paper we consider the problem u = λ u u + f in, u = u
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More informationb i (µ, x, s) ei ϕ(x) µ s (dx) ds (2) i=1
NONLINEAR EVOLTION EQATIONS FOR MEASRES ON INFINITE DIMENSIONAL SPACES V.I. Bogachev 1, G. Da Prato 2, M. Röckner 3, S.V. Shaposhnikov 1 The goal of this work is to prove the existence of a solution to
More informationTHEOREMS, ETC., FOR MATH 515
THEOREMS, ETC., FOR MATH 515 Proposition 1 (=comment on page 17). If A is an algebra, then any finite union or finite intersection of sets in A is also in A. Proposition 2 (=Proposition 1.1). For every
More informationGLOBAL EXISTENCE AND ENERGY DECAY OF SOLUTIONS TO A PETROVSKY EQUATION WITH GENERAL NONLINEAR DISSIPATION AND SOURCE TERM
Georgian Mathematical Journal Volume 3 (26), Number 3, 397 4 GLOBAL EXITENCE AND ENERGY DECAY OF OLUTION TO A PETROVKY EQUATION WITH GENERAL NONLINEAR DIIPATION AND OURCE TERM NOUR-EDDINE AMROUN AND ABBE
More informationAnalysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient
Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient Yifeng Xu 1 Jun Zou 2 Abstract Based on a new a posteriori error estimator, an adaptive finite element method is proposed
More informationMixed exterior Laplace s problem
Mixed exterior Laplace s problem Chérif Amrouche, Florian Bonzom Laboratoire de mathématiques appliquées, CNRS UMR 5142, Université de Pau et des Pays de l Adour, IPRA, Avenue de l Université, 64000 Pau
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationWeak solutions for the Cahn-Hilliard equation with degenerate mobility
Archive for Rational Mechanics and Analysis manuscript No. (will be inserted by the editor) Shibin Dai Qiang Du Weak solutions for the Cahn-Hilliard equation with degenerate mobility Abstract In this paper,
More informationNONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS
Fixed Point Theory, Volume 9, No. 1, 28, 3-16 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS GIOVANNI ANELLO Department of Mathematics University
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More information2 Statement of the problem and assumptions
Mathematical Notes, 25, vol. 78, no. 4, pp. 466 48. Existence Theorem for Optimal Control Problems on an Infinite Time Interval A.V. Dmitruk and N.V. Kuz kina We consider an optimal control problem on
More informationA Product Property of Sobolev Spaces with Application to Elliptic Estimates
Rend. Sem. Mat. Univ. Padova Manoscritto in corso di stampa pervenuto il 23 luglio 2012 accettato l 1 ottobre 2012 A Product Property of Sobolev Spaces with Application to Elliptic Estimates by Henry C.
More informationPERTURBATION THEORY FOR NONLINEAR DIRICHLET PROBLEMS
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 28, 2003, 207 222 PERTURBATION THEORY FOR NONLINEAR DIRICHLET PROBLEMS Fumi-Yuki Maeda and Takayori Ono Hiroshima Institute of Technology, Miyake,
More informationSYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS. M. Grossi S. Kesavan F. Pacella M. Ramaswamy. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 12, 1998, 47 59 SYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS M. Grossi S. Kesavan F. Pacella M. Ramaswamy
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationOptimal stopping time formulation of adaptive image filtering
Optimal stopping time formulation of adaptive image filtering I. Capuzzo Dolcetta, R. Ferretti 19.04.2000 Abstract This paper presents an approach to image filtering based on an optimal stopping time problem
More informationEXISTENCE RESULTS FOR QUASILINEAR HEMIVARIATIONAL INEQUALITIES AT RESONANCE. Leszek Gasiński
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 2007 pp. 409 418 EXISTENCE RESULTS FOR QUASILINEAR HEMIVARIATIONAL INEQUALITIES AT RESONANCE Leszek Gasiński Jagiellonian
More informationOn the bang-bang property of time optimal controls for infinite dimensional linear systems
On the bang-bang property of time optimal controls for infinite dimensional linear systems Marius Tucsnak Université de Lorraine Paris, 6 janvier 2012 Notation and problem statement (I) Notation: X (the
More informationFriedrich symmetric systems
viii CHAPTER 8 Friedrich symmetric systems In this chapter, we describe a theory due to Friedrich [13] for positive symmetric systems, which gives the existence and uniqueness of weak solutions of boundary
More informationSobolev Spaces 27 PART II. Review of Sobolev Spaces
Sobolev Spaces 27 PART II Review of Sobolev Spaces Sobolev Spaces 28 SOBOLEV SPACES WEAK DERIVATIVES I Given R d, define a multi index α as an ordered collection of integers α = (α 1,...,α d ), such that
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 207 The elastic-plastic torsion problem: a posteriori error estimates for approximate solutions
More informationViscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces
Viscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces YUAN-HENG WANG Zhejiang Normal University Department of Mathematics Yingbing Road 688, 321004 Jinhua
More informationMULTIPLE SOLUTIONS FOR CRITICAL ELLIPTIC PROBLEMS WITH FRACTIONAL LAPLACIAN
Electronic Journal of Differential Equations, Vol. 016 (016), No. 97, pp. 1 11. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu MULTIPLE SOLUTIONS
More informationON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN
ON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN XIAN LIAO Abstract. In this work we will show the global existence of the strong solutions of the inhomogeneous
More informationl(y j ) = 0 for all y j (1)
Problem 1. The closed linear span of a subset {y j } of a normed vector space is defined as the intersection of all closed subspaces containing all y j and thus the smallest such subspace. 1 Show that
More informationHamilton-Jacobi Equation and Optimality conditions for control systems governed by semilinear parabolic equations with boundary control
Hamilton-Jacobi Equation and Optimality conditions for control systems governed by semilinear parabolic equations with boundary control Marc Quincampoix and Catalin Trenchea Laboratoire de Mathématiques,
More informationON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM
Internat. J. Math. & Math. Sci. Vol. 22, No. 3 (999 587 595 S 6-72 9922587-2 Electronic Publishing House ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR
More informationON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS
Chin. Ann. Math.??B(?), 200?, 1 20 DOI: 10.1007/s11401-007-0001-x ON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS Michel CHIPOT Abstract We present a method allowing to obtain existence of a solution for
More informationt y n (s) ds. t y(s) ds, x(t) = x(0) +
1 Appendix Definition (Closed Linear Operator) (1) The graph G(T ) of a linear operator T on the domain D(T ) X into Y is the set (x, T x) : x D(T )} in the product space X Y. Then T is closed if its graph
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationA TWO PARAMETERS AMBROSETTI PRODI PROBLEM*
PORTUGALIAE MATHEMATICA Vol. 53 Fasc. 3 1996 A TWO PARAMETERS AMBROSETTI PRODI PROBLEM* C. De Coster** and P. Habets 1 Introduction The study of the Ambrosetti Prodi problem has started with the paper
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 3. Fundamental properties IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Example Consider the system ẋ = f
More informationEXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem
Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS
More informationINFINITE TIME HORIZON OPTIMAL CONTROL OF THE SEMILINEAR HEAT EQUATION
Nonlinear Funct. Anal. & Appl., Vol. 7, No. (22), pp. 69 83 INFINITE TIME HORIZON OPTIMAL CONTROL OF THE SEMILINEAR HEAT EQUATION Mihai Sîrbu Abstract. We consider here the infinite horizon control problem
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More informationA VARIATIONAL INEQUALITY RELATED TO AN ELLIPTIC OPERATOR
Proyecciones Vol. 19, N o 2, pp. 105-112, August 2000 Universidad Católica del Norte Antofagasta - Chile A VARIATIONAL INEQUALITY RELATED TO AN ELLIPTIC OPERATOR A. WANDERLEY Universidade do Estado do
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationNonlinear elliptic systems with exponential nonlinearities
22-Fez conference on Partial Differential Equations, Electronic Journal of Differential Equations, Conference 9, 22, pp 139 147. http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu
More informationON THE EXISTENCE OF THREE SOLUTIONS FOR QUASILINEAR ELLIPTIC PROBLEM. Paweł Goncerz
Opuscula Mathematica Vol. 32 No. 3 2012 http://dx.doi.org/10.7494/opmath.2012.32.3.473 ON THE EXISTENCE OF THREE SOLUTIONS FOR QUASILINEAR ELLIPTIC PROBLEM Paweł Goncerz Abstract. We consider a quasilinear
More informationPDE Constrained Optimization selected Proofs
PDE Constrained Optimization selected Proofs Jeff Snider jeff@snider.com George Mason University Department of Mathematical Sciences April, 214 Outline 1 Prelim 2 Thms 3.9 3.11 3 Thm 3.12 4 Thm 3.13 5
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationLocal null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain
Local null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain Nicolás Carreño Université Pierre et Marie Curie-Paris 6 UMR 7598 Laboratoire
More informationOn Asymptotic Variational Wave Equations
On Asymptotic Variational Wave Equations Alberto Bressan 1, Ping Zhang 2, and Yuxi Zheng 1 1 Department of Mathematics, Penn State University, PA 1682. E-mail: bressan@math.psu.edu; yzheng@math.psu.edu
More informationStability for solution of Differential Variational Inequalitiy
Stability for solution of Differential Variational Inequalitiy Jian-xun Zhao Department of Mathematics, School of Science, Tianjin University Tianjin 300072, P.R. China Abstract In this paper we study
More informationASYMPTOTIC BEHAVIOUR OF NONLINEAR ELLIPTIC SYSTEMS ON VARYING DOMAINS
ASYMPTOTIC BEHAVIOUR OF NONLINEAR ELLIPTIC SYSTEMS ON VARYING DOMAINS Juan CASADO DIAZ ( 1 ) Adriana GARRONI ( 2 ) Abstract We consider a monotone operator of the form Au = div(a(x, Du)), with R N and
More informationBLOW-UP AND EXTINCTION OF SOLUTIONS TO A FAST DIFFUSION EQUATION WITH HOMOGENEOUS NEUMANN BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 016 (016), No. 36, pp. 1 10. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu BLOW-UP AND EXTINCTION OF SOLUTIONS TO A FAST
More informationOBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 27 (27, No. 6, pp. 2. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu OBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS
More informationOn the Solvability Conditions for a Linearized Cahn-Hilliard Equation
Rend. Istit. Mat. Univ. Trieste Volume 43 (2011), 1 9 On the Solvability Conditions for a Linearized Cahn-Hilliard Equation Vitaly Volpert and Vitali Vougalter Abstract. We derive solvability conditions
More informationA generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem
A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem Dave McCormick joint work with James Robinson and José Rodrigo Mathematics and Statistics Centre for Doctoral Training University
More informationEXISTENCE OF WEAK SOLUTIONS FOR A NONUNIFORMLY ELLIPTIC NONLINEAR SYSTEM IN R N. 1. Introduction We study the nonuniformly elliptic, nonlinear system
Electronic Journal of Differential Equations, Vol. 20082008), No. 119, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu login: ftp) EXISTENCE
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationTHE CAUCHY PROBLEM AND STEADY STATE SOLUTIONS FOR A NONLOCAL CAHN-HILLIARD EQUATION
Electronic Journal of Differential Equations, Vol. 24(24), No. 113, pp. 1 9. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) THE CAUCHY
More informationUNIFORM DECAY OF SOLUTIONS FOR COUPLED VISCOELASTIC WAVE EQUATIONS
Electronic Journal of Differential Equations, Vol. 16 16, No. 7, pp. 1 11. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu UNIFORM DECAY OF SOLUTIONS
More informationA POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS
A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS Z. LU Communicated by Gabriela Marinoschi In this paper, we discuss a
More informationConvergence Rates in Regularization for Nonlinear Ill-Posed Equations Involving m-accretive Mappings in Banach Spaces
Applied Mathematical Sciences, Vol. 6, 212, no. 63, 319-3117 Convergence Rates in Regularization for Nonlinear Ill-Posed Equations Involving m-accretive Mappings in Banach Spaces Nguyen Buong Vietnamese
More informationExistence Results for Multivalued Semilinear Functional Differential Equations
E extracta mathematicae Vol. 18, Núm. 1, 1 12 (23) Existence Results for Multivalued Semilinear Functional Differential Equations M. Benchohra, S.K. Ntouyas Department of Mathematics, University of Sidi
More informationarxiv: v1 [math.oc] 22 Sep 2016
EUIVALENCE BETWEEN MINIMAL TIME AND MINIMAL NORM CONTROL PROBLEMS FOR THE HEAT EUATION SHULIN IN AND GENGSHENG WANG arxiv:1609.06860v1 [math.oc] 22 Sep 2016 Abstract. This paper presents the equivalence
More informationEXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC SCHRÖDINGER EQUATIONS
Electronic Journal of Differential Equations, Vol. 017 (017), No. 15, pp. 1 7. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationMinimal periods of semilinear evolution equations with Lipschitz nonlinearity
Minimal periods of semilinear evolution equations with Lipschitz nonlinearity James C. Robinson a Alejandro Vidal-López b a Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. b Departamento
More informationTitle of your thesis here
Title of your thesis here Your name here A Thesis presented for the degree of Doctor of Philosophy Your Research Group Here Department of Mathematical Sciences University of Durham England Month and Year
More informationWEAK CONVERGENCE OF RESOLVENTS OF MAXIMAL MONOTONE OPERATORS AND MOSCO CONVERGENCE
Fixed Point Theory, Volume 6, No. 1, 2005, 59-69 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.htm WEAK CONVERGENCE OF RESOLVENTS OF MAXIMAL MONOTONE OPERATORS AND MOSCO CONVERGENCE YASUNORI KIMURA Department
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationFractional Cahn-Hilliard equation
Fractional Cahn-Hilliard equation Goro Akagi Mathematical Institute, Tohoku University Abstract In this note, we review a recent joint work with G. Schimperna and A. Segatti (University of Pavia, Italy)
More informationTake-Home Final Examination of Math 55a (January 17 to January 23, 2004)
Take-Home Final Eamination of Math 55a January 17 to January 3, 004) N.B. For problems which are similar to those on the homework assignments, complete self-contained solutions are required and homework
More informationNonexistence of solutions for quasilinear elliptic equations with p-growth in the gradient
Electronic Journal of Differential Equations, Vol. 2002(2002), No. 54, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Nonexistence
More informationGlobal weak solutions to a sixth order Cahn-Hilliard type equation
Global weak solutions to a sixth order Cahn-Hilliard type equation M. D. Korzec, P. Nayar, P. Rybka December 8, 21 Abstract In this paper we study a sixth order Cahn-Hilliard type equation that arises
More informationExistence and Uniqueness Results for Nonlinear Implicit Fractional Differential Equations with Boundary Conditions
Existence and Uniqueness Results for Nonlinear Implicit Fractional Differential Equations with Boundary Conditions Mouffak Benchohra a,b 1 and Jamal E. Lazreg a, a Laboratory of Mathematics, University
More informationINFINITELY MANY SOLUTIONS FOR A CLASS OF SUBLINEAR SCHRÖDINGER EQUATIONS
Journal of Applied Analysis and Computation Volume 8, Number 5, October 018, 1475 1493 Website:http://jaac-online.com/ DOI:10.11948/018.1475 INFINITELY MANY SOLUTIONS FOR A CLASS OF SUBLINEAR SCHRÖDINGER
More informationRecent Trends in Differential Inclusions
Recent Trends in Alberto Bressan Department of Mathematics, Penn State University (Aveiro, June 2016) (Aveiro, June 2016) 1 / Two main topics ẋ F (x) differential inclusions with upper semicontinuous,
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied athematics http://jipam.vu.edu.au/ Volume 4, Issue 5, Article 98, 2003 ASYPTOTIC BEHAVIOUR OF SOE EQUATIONS IN ORLICZ SPACES D. ESKINE AND A. ELAHI DÉPARTEENT
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 133 Higher order variational problems on two-dimensional domains Michael Bildhauer and Martin
More informationarxiv: v2 [math.ap] 28 Nov 2016
ONE-DIMENSIONAL SAIONARY MEAN-FIELD GAMES WIH LOCAL COUPLING DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES arxiv:1611.8161v [math.ap] 8 Nov 16 Abstract. A standard assumption in mean-field game
More informationOptimal Control Algorithm for Complex Heat Transfer Model
Optimal Control Algorithm for Complex Heat Transfer Model Alexander Chebotarev, Gleb Grenkin, and Andrey Kovtanyuk Far Eastern Federal University, Sukhanova st. 8, 690950 Vladivostok, Russia, 2 Institute
More informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE XIAOBING FENG AND ANDREAS PROHL Abstract. In this
More informationBLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED
BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED TAOUFIK HMIDI AND SAHBI KERAANI Abstract. In this note we prove a refined version of compactness lemma adapted to the blowup analysis
More information