<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2012.312258</article-id><article-id pub-id-type="publisher-id">AM-25444</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Pontryagin’s Maximum Principle for a Advection-Diffusion-Reaction Equation
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>oujun</surname><given-names>Xu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Cuie</surname><given-names>Xiao</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hui</surname><given-names>Zhu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Department of Mathematics and Computation Sciences, Hunan City University, Yiyang, China</addr-line></aff><aff id="aff1"><addr-line>School of Mathematics and Physics, University of South China, Hengyang, China</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>youjunxu@163.com(OX)</email>;<email>xiaocuie@163.com(CX)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>12</day><month>12</month><year>2012</year></pub-date><volume>03</volume><issue>12</issue><fpage>1888</fpage><lpage>1891</lpage><history><date date-type="received"><day>July</day>	<month>2,</month>	<year>2012</year></date><date date-type="rev-recd"><day>November</day>	<month>19,</month>	<year>2012</year>	</date><date date-type="accepted"><day>November</day>	<month>26,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  In this paper we investigate optimal control problems governed by a advection-diffusion-reaction equation. We present a method for deriving conditions in the form of Pontryagin’s principle. The main tools used are the Ekeland’s variational principle combined with penalization and spike variation techniques.
 
</p></abstract><kwd-group><kwd>Optimal Control; Pontryagin’s Maximum Principle; State Constraint</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Consider the following controlled advection convection diffusion equations:</p><disp-formula id="scirp.25444-formula142267"><label>(1.1)</label><graphic position="anchor" xlink:href="8-20870\33ab45c1-0bbb-479c-80f9-410c28594c1d.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="8-20870\f707a5a2-bc9a-4a18-bc91-9362099fef23.jpg" /> is a convex bounded domain with a smooth boundary<img src="8-20870\93de5408-b2cd-4448-8709-82490c94e0a2.jpg" />, the diffusity <img src="8-20870\40f784d4-047e-4b73-b35a-871ea819af90.jpg" /> with</p><p><img src="8-20870\1ef89b0d-db7f-49a0-8a69-8625c999fbbe.jpg" />, the reaction <img src="8-20870\c4aa892b-9102-4ead-946c-50f6b35c79a8.jpg" /> with</p><p><img src="8-20870\0d77043b-4548-4ab1-8cb8-7c8890d9802c.jpg" />, and the advective field<img src="8-20870\bf63daf7-e0a4-45cd-934d-855a4b83d892.jpg" />with <img src="8-20870\2fa51a5e-c5e1-4adc-868f-7e58ab41a71c.jpg" /> and <img src="8-20870\e35efaeb-361a-44f4-83a3-a3382cbd9a6f.jpg" /> are assigned functions. Here<img src="8-20870\ea04e015-1fa1-4870-ac20-ccaa94a983db.jpg" />, with <img src="8-20870\e9283282-af42-42e0-94be-0bcbe8940e55.jpg" /> being a separable metric space. Function<img src="8-20870\0dcc0d8b-5461-43a7-8231-8ad194046051.jpg" />, called a control, is taken from the set</p><p><img src="8-20870\9893680c-2f73-4666-aba4-34fb6260d8c3.jpg" /></p><p>Under some mild conditions, for any<img src="8-20870\984e2e37-76fc-4d52-a541-8d6ed435d0cd.jpg" />, (1.1) admits a unique weak solution <img src="8-20870\bf880ac3-8369-4a96-96de-00c73ebce051.jpg" /> which is called the state(corresponding to the control<img src="8-20870\0ae4ff8b-34ad-47dc-b8d1-9b32b6031f7f.jpg" />). The performance of the control is measured by the cost functional</p><disp-formula id="scirp.25444-formula142268"><label>(1.2)</label><graphic position="anchor" xlink:href="8-20870\63237836-bbd7-4b79-8688-a4e503ff4f7c.jpg"  xlink:type="simple"/></disp-formula><p>for some given map<img src="8-20870\75f8cab8-3c0f-4d18-95f7-bbac35c72169.jpg" />. Our optimal control problem can be stated as follows.</p><p>Problem (C). Find a <img src="8-20870\8a5c8af7-73a4-4f46-80ae-ba2c96d24ffd.jpg" /> such that</p><disp-formula id="scirp.25444-formula142269"><label>(1.3)</label><graphic position="anchor" xlink:href="8-20870\6ce8a3f8-9145-47cb-b56f-472262d10829.jpg"  xlink:type="simple"/></disp-formula><p>And the state constraint of form:</p><disp-formula id="scirp.25444-formula142270"><label>(1.4)</label><graphic position="anchor" xlink:href="8-20870\e9295c3c-782a-40ef-a2d8-822f1f05d8fb.jpg"  xlink:type="simple"/></disp-formula><p>In this paper, we make the following assumptions.</p><p>(H1) Set <img src="8-20870\be14a55d-2666-4cd7-868e-4fb28d9067fa.jpg" /> is a convex bounded domain with a smooth boundary<img src="8-20870\e7ffd08d-931d-49a3-bbb2-509b12ce5f04.jpg" />.</p><p>(H2) Set <img src="8-20870\e8757fc0-3ee8-4ccc-8a77-752ac1c5e38e.jpg" /> is a separable metric space.</p><p>(H3) The function <img src="8-20870\b644f45b-4f39-4350-8887-c02bfa144284.jpg" /> has the following properties: <img src="8-20870\2e4b42cf-d1d7-43df-b478-06fa23a73dbb.jpg" />is measurable on<img src="8-20870\bdf4a29c-ea55-433f-a5ac-7f952b0f2dc1.jpg" />, and <img src="8-20870\7053f7bb-d4f1-4f1c-a446-52e023524daa.jpg" /> continuous on <img src="8-20870\1e05d93f-b658-4b6e-b00a-09d300cb7276.jpg" /> and for any<img src="8-20870\c0cb1a27-2f98-4992-954c-7c6a56e0aa9b.jpg" />, a constant</p><p><img src="8-20870\56bdba48-753c-4b68-b16e-e0ff88a7a3d6.jpg" />, such that <img src="8-20870\19a91e0c-8c46-472a-94ac-a83e0a43821f.jpg" /></p><p>(H4) Function<img src="8-20870\2ae2cd9c-52f2-4657-9c94-b1e671b00599.jpg" /> is measurable in <img src="8-20870\b8cfcf44-1968-429d-9e13-5d5192315675.jpg" />and continuous in<img src="8-20870\9873e3b1-5a5e-475e-bb16-4f01542abfd2.jpg" /> for almost all<img src="8-20870\1ff1bc58-9a0e-49ae-bcc0-992096e1ecd3.jpg" />. Moreover, for any<img src="8-20870\406509cb-9ebb-4803-84d5-0f8174fc9030.jpg" />, there exists a <img src="8-20870\00c045f7-ac85-4317-8736-33d01c69e89b.jpg" /> such that</p><disp-formula id="scirp.25444-formula142271"><label>(1.5)</label><graphic position="anchor" xlink:href="8-20870\6447b6b6-b77a-4d8e-b3d3-38fb49fb560f.jpg"  xlink:type="simple"/></disp-formula><p>(H5) <img src="8-20870\77070c39-39ad-45c5-8c21-f682af321631.jpg" />is a Banach space with strictly convex dual<img src="8-20870\2c4e7811-dfd2-455a-817f-afc4a3b30dac.jpg" />, <img src="8-20870\9d7b53ac-b8a7-4533-97b3-445cc13ac402.jpg" />is continuously Fr&#233;chet differentiable, and <img src="8-20870\c227f518-1fb1-4ed0-b72b-9eb8a3d36224.jpg" /> is closed and convex set.</p><p>(H6) <img src="8-20870\e02d4c09-26d9-448a-ac03-3586116c1d94.jpg" />has finite condimensionality in <img src="8-20870\bbe2c228-7674-4744-aa1f-9fab70246997.jpg" /> for some<img src="8-20870\b785ca95-8d3b-4375-b1ba-8e61c1a35e45.jpg" />, where<img src="8-20870\344ace12-fe78-4b3c-847e-93c33688dbc4.jpg" />.</p><p>Definition 1.1 (see [<xref ref-type="bibr" rid="scirp.25444-ref1">1</xref>]) Let <img src="8-20870\5d13eea2-1759-44e7-9dd5-5da10b3f00c3.jpg" /> is a Banach space and <img src="8-20870\26615773-31c1-49c2-b0f0-b50f23d74c33.jpg" /> is a subspace of<img src="8-20870\5d5e748b-2947-4160-b108-88720e2fed0a.jpg" />. We say that <img src="8-20870\cbefe66f-a7c5-4334-b1b6-a587690b47ce.jpg" /> is finite codimensional in <img src="8-20870\e027c4da-bd02-4380-9218-45c12292473a.jpg" /> if there exists <img src="8-20870\7780d15b-0cc8-4907-8d37-4a9c107b676d.jpg" /> such that</p><p><img src="8-20870\3ded39af-cafc-48db-aa82-0dfb77939627.jpg" /></p><p>A subset <img src="8-20870\2cd6022d-08b6-4787-8bda-d9e8ba478859.jpg" /> of <img src="8-20870\0260bd52-2539-4f53-8b47-c9934e254c4b.jpg" /> is said to be finite codimensional in <img src="8-20870\d9358cbf-eafd-4507-ace4-34cb10f7cb20.jpg" /> if for some<img src="8-20870\fea850d8-064a-4b42-b46d-82a61cb837a4.jpg" />, <img src="8-20870\947c807e-7086-461d-a3e1-53d5ce1e1d79.jpg" />the closed subspace spanned by <img src="8-20870\350f7469-2320-49d7-8b1a-c5fe13e9ecf5.jpg" /> is a finite codimensional subspace of <img src="8-20870\8a257645-1d98-49b9-a3f2-fc857a2b8272.jpg" /> and <img src="8-20870\55220322-0210-4b81-a8f1-00889b194667.jpg" /> the closed convex hull of <img src="8-20870\e3afafb8-dd60-408c-a459-9530249d3fd7.jpg" /> has a nonempty interior in this subspace.</p><p>Lemma 1.2. Let (H1) - (H3) hold. Then, for any<img src="8-20870\3ece5200-10ca-43b5-99d7-33fad33d25fa.jpg" />, (1.1) admits a unique weak solution</p><p><img src="8-20870\779209f4-f84d-4296-971a-500d3452b323.jpg" />.</p><p>Furthermore, there exists a constant<img src="8-20870\e7b2e442-4a38-4fb6-bfd5-17796b48ca54.jpg" />, independent of</p><disp-formula id="scirp.25444-formula142272"><label>(1.6)</label><graphic position="anchor" xlink:href="8-20870\ac8c8690-d4ed-44af-8223-f1550a81d5c4.jpg"  xlink:type="simple"/></disp-formula><p>The weak solution <img src="8-20870\80da646e-5353-4b4c-879c-0965d5270833.jpg" /> of the state Equation (1.1) is determined by</p><p><img src="8-20870\8392f7d9-fb1f-49bb-a1c2-a4d4c3cf348d.jpg" /></p><p>using the bilinear form <img src="8-20870\bb0d353c-031a-4f71-aece-e2d9b6c33c00.jpg" /> given by</p><p><img src="8-20870\6a0f7096-c2dd-4c31-aa35-3c0b5dba1998.jpg" /></p><p>Existence and uniqueness of the solution to (1.1) follow from the above hypotheses on the problem data (see [<xref ref-type="bibr" rid="scirp.25444-ref2">2</xref>]). Let <img src="8-20870\f25c7edf-21f2-4a6d-be80-bb8dad3360cb.jpg" /> be the set of all pairs <img src="8-20870\be2f3277-7d3d-4f18-9fd5-638bcae32d41.jpg" /> satisfying (1.1) and (1.4) is called an admissible set. Any <img src="8-20870\a75591c0-d37b-4686-878b-333403eaf0ee.jpg" /> is called an admissible pair. The pair</p><p><img src="8-20870\adde4b71-6646-45e0-bb9c-c8ec7a8d5741.jpg" />, moveover satisfies<img src="8-20870\df643b32-20f0-4ba2-b4fa-f5a3013dc5e2.jpg" /></p><p>for all <img src="8-20870\6109c821-d881-4748-99ef-b0ebbb7f6a32.jpg" /> is called an optimal pair. If it exists, refer to <img src="8-20870\0391e699-e1ba-4d77-bb76-7efb579d8cd2.jpg" /> and <img src="8-20870\40412a55-2648-48fa-b960-05fc19768657.jpg" /> as an optimal state and control, respectively.</p><p>Now, let <img src="8-20870\7bb23ba8-ac53-4b74-bdaa-4349f7aa5d51.jpg" /> be an optimal pair of Problem (C).</p><p>Let <img src="8-20870\997ddd12-6084-4ee6-ba21-37dfe042e621.jpg" />be the unique solution of the following problem:</p><disp-formula id="scirp.25444-formula142273"><label>(1.7)</label><graphic position="anchor" xlink:href="8-20870\d6cc2014-8328-48db-9c76-9f4da31ac149.jpg"  xlink:type="simple"/></disp-formula><p>And define the reachable set of variational system (1.7)</p><disp-formula id="scirp.25444-formula142274"><label>(1.8)</label><graphic position="anchor" xlink:href="8-20870\00123330-54f9-4778-9f00-5211b3b7e92d.jpg"  xlink:type="simple"/></disp-formula><p>Now, let us state the first order necessary conditions of an optimal control to Problem (C) as follows.</p><p>Theorem 1.3. (Pontryagin’s maximum principle) Let (H1) - (H6) hold. Let <img src="8-20870\fbda6ad1-de8a-4820-95e6-d0ad998073e9.jpg" /> be an optimal pair of Problem (C). Then there exists a triplet</p><p><img src="8-20870\17a7bb37-36cb-4ae8-8719-ff890f4f4710.jpg" /></p><p>such that</p><disp-formula id="scirp.25444-formula142275"><label>(19)</label><graphic position="anchor" xlink:href="8-20870\d1e7956d-313e-44c2-8a63-649850fb657b.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.25444-formula142276"><label>(1.10)</label><graphic position="anchor" xlink:href="8-20870\1ce66585-a8c1-428e-810c-b861dc6cc7e7.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.25444-formula142277"><label>(1.11)</label><graphic position="anchor" xlink:href="8-20870\068a3dc6-8094-4509-a711-ee74992b6d1e.jpg"  xlink:type="simple"/></disp-formula><p>where</p><p><img src="8-20870\8cd65d82-44d2-4dfe-97c4-bf5ca7ae6ed2.jpg" /></p><p>(1.9), (1.10), and (1.11) are called the transversality condition, the adjoint system(along the given optimal pair), and the maximum condition, respectively.</p><p>Many authors (Dede [<xref ref-type="bibr" rid="scirp.25444-ref3">3</xref>], Yan [<xref ref-type="bibr" rid="scirp.25444-ref4">4</xref>], Becker [<xref ref-type="bibr" rid="scirp.25444-ref5">5</xref>], Stefano [<xref ref-type="bibr" rid="scirp.25444-ref6">6</xref>], Collis [<xref ref-type="bibr" rid="scirp.25444-ref7">7</xref>]) have already considered control problems for convection-diffusion equations from theoretical or numerical point of view. In the work mentioned above, the control set is convex. However, in many practical cases, the control set can not convex. This stimulates us to study Problem (C). To get Pontryagin’s Principle, we use a method based on penalization of state constraints, and Ekeland’s principle combined with diffuse perturbations [<xref ref-type="bibr" rid="scirp.25444-ref8">8</xref>].</p><p>In the next section, we will prove Pontryagin’s maximum principle of optimal control of Problem (C).</p></sec><sec id="s2"><title>2. Proof of the Maximum Principle</title><p>This section is devoted to the proof of the maximum principle.</p><p>Proof of Theorem 1.3. Firstly, let</p><p><img src="8-20870\1c391a82-d4ee-4e54-afcc-5357ba423e16.jpg" />where <img src="8-20870\3e66ab6a-628b-49d5-86c9-bf815e250872.jpg" /> is the Lebesgue measure of<img src="8-20870\cfbad99f-cdfa-4b1e-bafb-abc0843cf8d3.jpg" />. We can easily prove that <img src="8-20870\39bd18cf-fa43-466b-9fd8-5f100f2250a8.jpg" /> is a complete metric space. Let <img src="8-20870\91bb1435-bcb2-4ba4-9366-04bd50a20e1c.jpg" /> be anoptimal pair of Problem (C). For any <img src="8-20870\9f5cff11-c4aa-4a5e-b379-a086f3e335fb.jpg" /> be the corresponding state, emphasizing the dependence on the control. Without loss of generality, we may assume that<img src="8-20870\7d56f2ba-143d-4fab-95e5-ff80b9eef90f.jpg" />. For any <img src="8-20870\b0990437-f659-4bfa-8817-53d08787e7a8.jpg" /> define</p><disp-formula id="scirp.25444-formula142278"><label>(2.1)</label><graphic position="anchor" xlink:href="8-20870\53df72e9-a6e1-467d-af63-1b5bf086654b.jpg"  xlink:type="simple"/></disp-formula><p>where<img src="8-20870\8ec263a5-daee-4bc0-bf2a-e1d67a599e94.jpg" />, and <img src="8-20870\20267ec3-74cd-4519-8411-e458746ac74d.jpg" /> is an optimal control.</p><p>Clearly, this function is continuous on the (complete) metric space<img src="8-20870\c8358731-ab7d-4a89-9e8f-d0e6b817317d.jpg" />. Also, we have</p><disp-formula id="scirp.25444-formula142279"><label>(2.2)</label><graphic position="anchor" xlink:href="8-20870\c867e057-27d2-48f4-aca3-dcd5a4d1d19c.jpg"  xlink:type="simple"/></disp-formula><p>Hence, by Ekeland’s variational principle, we can find a<img src="8-20870\e38c36fd-dee8-453a-aa76-5a1f7c905f79.jpg" />, such that</p><disp-formula id="scirp.25444-formula142280"><label>(2.3)</label><graphic position="anchor" xlink:href="8-20870\b39a032d-2d3a-4c38-b962-e8ac6bd62d5a.jpg"  xlink:type="simple"/></disp-formula><p>Let <img src="8-20870\42e98322-a056-4d5b-8226-4958ca800902.jpg" /> and <img src="8-20870\98f75e4b-12b9-441a-aecc-12d81be976d7.jpg" /> be fixed and let<img src="8-20870\049ff372-7afd-4607-82aa-465ba2dec632.jpg" />, we know that for any<img src="8-20870\440014db-eca4-4766-8657-f1539d390b4b.jpg" />, there exists a measurable set <img src="8-20870\175295e0-a4ba-4243-a52d-37d182bd78a0.jpg" /> with the property <img src="8-20870\0d7340c0-5e20-453a-9701-c67e4af29edb.jpg" /> such that if we define</p><p><img src="8-20870\2dd01b8c-0c8c-4a09-9b33-77191956e176.jpg" /></p><p>and let <img src="8-20870\48d0386e-7e87-4015-9071-ce06c688da85.jpg" /> be the corresponding state, then</p><disp-formula id="scirp.25444-formula142281"><label>(2.4)</label><graphic position="anchor" xlink:href="8-20870\5f0034b1-b7d4-4b7f-82f8-f0fc83b391db.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="8-20870\8d6846b3-7954-4009-abe3-4dc063d04947.jpg" /> and <img src="8-20870\c5756a1b-55a0-4367-a721-9ce4c97d32e1.jpg" /> satisfying the following</p><disp-formula id="scirp.25444-formula142282"><label>(2.5)</label><graphic position="anchor" xlink:href="8-20870\65d2aa31-098e-493e-baa4-33de266b72e6.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.25444-formula142283"><label>(2.6)</label><graphic position="anchor" xlink:href="8-20870\5128324d-ec09-4dc5-9265-0e034b0d48c4.jpg"  xlink:type="simple"/></disp-formula><p>with</p><disp-formula id="scirp.25444-formula142284"><label>(2.7)</label><graphic position="anchor" xlink:href="8-20870\42bb3266-3ae2-48ec-a571-d88b5e519f27.jpg"  xlink:type="simple"/></disp-formula><p>We take<img src="8-20870\7736c708-f3e6-4b67-b71d-6e6b1be5a369.jpg" />. It follows that</p><disp-formula id="scirp.25444-formula142285"><label>(2.8)</label><graphic position="anchor" xlink:href="8-20870\9f8b47d4-4054-4148-8a9c-a467bf4148cc.jpg"  xlink:type="simple"/></disp-formula><p>where</p><p><img src="8-20870\13e98468-32c2-4fdd-a362-b03a72039d4b.jpg" /></p><p><img src="8-20870\e762856a-4f06-487b-9811-60e22dc12e0b.jpg" />denotes the subdifferential of<img src="8-20870\aea0d8f8-6d36-4cd2-841d-aff9ad16e6eb.jpg" />.</p><p>Next, we define <img src="8-20870\a0d341a3-1a13-4a3c-a5d7-3e696b42799c.jpg" /> as follows:</p><disp-formula id="scirp.25444-formula142286"><label>(2.9)</label><graphic position="anchor" xlink:href="8-20870\c88657d7-f5bc-4877-b309-094dfc3ac1af.jpg"  xlink:type="simple"/></disp-formula><p>By (2.1) and chapter 4 of [<xref ref-type="bibr" rid="scirp.25444-ref8">8</xref>], (2.8) becomes</p><disp-formula id="scirp.25444-formula142287"><label>(2.10)</label><graphic position="anchor" xlink:href="8-20870\c9f26089-d7e5-44ef-b028-76b964e2216a.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.25444-formula142288"><label>(2.11)</label><graphic position="anchor" xlink:href="8-20870\57b969f4-167e-4347-b25f-db974a772402.jpg"  xlink:type="simple"/></disp-formula><p>On the other hand, by the definition of the subdifferential, we have</p><disp-formula id="scirp.25444-formula142289"><label>(2.12)</label><graphic position="anchor" xlink:href="8-20870\0fe0ef50-1f10-4c23-94db-8b7c45f982fe.jpg"  xlink:type="simple"/></disp-formula><p>Next, from the first relation in (2.3) and by some calculations, we have</p><disp-formula id="scirp.25444-formula142290"><label>(2.13)</label><graphic position="anchor" xlink:href="8-20870\91d1e2bc-d3ff-4803-8e2e-fd01fc03f576.jpg"  xlink:type="simple"/></disp-formula><p>Consequently,</p><disp-formula id="scirp.25444-formula142291"><label>(2.14)</label><graphic position="anchor" xlink:href="8-20870\90e4c32e-cded-42af-b04f-70a6bd19ad0f.jpg"  xlink:type="simple"/></disp-formula><p>From (2.5) and (2.6), we have</p><disp-formula id="scirp.25444-formula142292"><label>(2.15)</label><graphic position="anchor" xlink:href="8-20870\b15add32-8620-4d0a-959a-2955ad6be9d2.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="8-20870\72596e92-5cdd-4a65-a715-69856f854836.jpg" /> is the solution of system (1.7) and</p><disp-formula id="scirp.25444-formula142293"><label>(2.16)</label><graphic position="anchor" xlink:href="8-20870\7460cb02-9749-404f-9263-32c6247bba39.jpg"  xlink:type="simple"/></disp-formula><p>From (2.10), (2.12) and (2.15), we have</p><disp-formula id="scirp.25444-formula142294"><label>(2.17)</label><graphic position="anchor" xlink:href="8-20870\31f88f6b-9827-4103-8a67-5bcd40c47ba9.jpg"  xlink:type="simple"/></disp-formula><p>with <img src="8-20870\8069f3b3-a625-48e2-819d-cf6661414017.jpg" /> Because <img src="8-20870\403c4827-d562-45ab-a294-b4383f8de772.jpg" /> has finite condimensionality in<img src="8-20870\a962d08a-596d-432b-915a-f17947803f5c.jpg" />, we can extract some subsequence, still denoted by itself, such that</p><p><img src="8-20870\3a98d11b-4e53-4fa9-9a0f-ac8ae810d988.jpg" /></p><p>From (2.17), we have</p><disp-formula id="scirp.25444-formula142295"><label>(2.18)</label><graphic position="anchor" xlink:href="8-20870\39b6ad7d-1c6b-41c2-9c9f-b8590a3c879b.jpg"  xlink:type="simple"/></disp-formula><p>Now, let</p><p><img src="8-20870\0429faaf-418e-4e23-a1f1-bff14e6f9a49.jpg" />.</p><p>Then</p><p><img src="8-20870\1b5592d2-f35b-493a-8705-427d1078d34c.jpg" />.</p><p>Then we have</p><disp-formula id="scirp.25444-formula142296"><label>(2.19)</label><graphic position="anchor" xlink:href="8-20870\13051e64-4082-4693-8546-f9d9a8c498cc.jpg"  xlink:type="simple"/></disp-formula><p>Take<img src="8-20870\814606f5-5b1a-41d3-833c-07503b5d9d62.jpg" />, we obtain (1.9).</p><p>Next, we let <img src="8-20870\c1730771-6fa9-4d63-8982-f88526ef77d9.jpg" /> to get</p><disp-formula id="scirp.25444-formula142297"><label>(2.20)</label><graphic position="anchor" xlink:href="8-20870\1ae1af80-9f8d-41b2-9d38-1515c4ae11a7.jpg"  xlink:type="simple"/></disp-formula><p>Because<img src="8-20870\32361672-9c91-4167-b6d0-042c0d89f67f.jpg" />, for the given<img src="8-20870\655e25b4-d798-414d-bedf-d6f2727352d8.jpg" />, there exists a unique solution <img src="8-20870\1e113322-47b0-4cdd-9156-f425b5d71fae.jpg" /> of the adjoint Equation (1.10). Then, from (1.6), (2.16), and (2.2), we have</p><disp-formula id="scirp.25444-formula142298"><label>(2.21)</label><graphic position="anchor" xlink:href="8-20870\8a5f1dd9-ce50-4ecf-877f-f98daeae5559.jpg"  xlink:type="simple"/></disp-formula><p>There, (1.11) follows. Finally, by (1.10), if <img src="8-20870\0080a761-b51d-46a8-b045-0d3075efc522.jpg" />, then<img src="8-20870\f0f252a7-7c76-4494-ae25-fcd09b1617d6.jpg" />. Thus, in the case where</p><p><img src="8-20870\6858e7bd-bffc-41c8-ab45-97bfbb16e727.jpg" /></p><p>we must have<img src="8-20870\a82709a3-0e2c-46c7-9039-e0c5ea321478.jpg" />, because<img src="8-20870\cff3142e-3864-4083-a68a-d556846d820a.jpg" />.</p></sec><sec id="s3"><title>3. Conclusion</title><p>We have already attained Pontryagin’s Maximum Principle for the advection-diffusion-reaction equation. It seems to us that this method can be used in treating many other relevant problems.</p></sec><sec id="s4"><title>REFERENCES</title></sec><sec id="s5"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.25444-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">H. W. Lou and J. M. Yong, “Optimal Controls for Semilinear Elliptic Equations with Leaing Term Containing Controls,” SIAM Journal on Control and Optimization, Vol. 48, No. 4, 2009, pp. 2366-2387.  
doi:10.1137/080740301</mixed-citation></ref><ref id="scirp.25444-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">J. T. Oden and J. N. Reddy, “Variational Methods in Theoretical Mechanics,” Springer, Berlin and Heidelberg, 1983. doi:10.1007/978-3-642-68811-9</mixed-citation></ref><ref id="scirp.25444-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">L. Dede and A. Quarteroni, “Optimal Control and Numerical Adaptivity for Advection-Diffusion Equations,” Mathematical Modelling and Numerical Analysis, Vol. 39, No. 2, 2005, pp. 1019-1040.</mixed-citation></ref><ref id="scirp.25444-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">N. N. Yan and Z. J. Zhou, “A Priori and a Posteriori Error Analysis of Edge Stabilization Galerkin Method for the Optimal Control Problem Governed by Convection-Dominated Diffusion Equation,” Journal of Computational and Applied Mathematics, Vol. 223, No. 1, 2009, pp. 198-217. doi:10.1016/j.cam.2008.01.006</mixed-citation></ref><ref id="scirp.25444-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">R. Becker and B. Vexler, “Optimal Control of the Convection-Diffusion Equation Using Stabilized Finite Element Methods,” Numerische Mathematik, Vol. 106, No. 3, 2007, pp. 349-367. doi:10.1007/s00211-007-0067-0</mixed-citation></ref><ref id="scirp.25444-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">S. Micheletti and S. Perotto, “An Anisotropic Mesh Adaptation Procedure for an Optimal Control Problem of the Advection-Diffusion-Reaction Equation,” MOX-Report No. 15, 2008.</mixed-citation></ref><ref id="scirp.25444-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">S. S. Collis and M. Heinkenschloss, “Analysis of the Streamline Upwind/Petrov Galerkin Method Applied to the Solution of Optimal Control Problems,” Technical Report 02-01, Department of Computational and Applied Mathematics, Rice University, Houston, 2002.</mixed-citation></ref><ref id="scirp.25444-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">X. Li and J. Yong, “Optimal Control Theory for InfiniteDimensional Systems,” Birkh?user, Boston, 1995.</mixed-citation></ref></ref-list></back></article>