Boundary layer type with initial value problem











up vote
0
down vote

favorite












Consider the initial value problem
$sqrt{epsilon} , u'' + u' - u = e^{2t}$ , with $u(0)=1$, and $u'(0)=1/sqrt{epsilon}$.



I am trying to use a matched asymptotic expansion to find the leading term of the solution.



First by setting $epsilon=0$, and using the initial condition $u(0)=1$, we get the outer solution $u_0(t)=c_1 e^t + e^{2t}$.



I set the change of variables $tau=t/sqrt{epsilon}$.



Writing the original equation in terms of the new variable I got:



$U'' + U' -sqrt{epsilon} , U = sqrt{epsilon} , e^{2tau epsilon^{1/2}}$



with $U(0)=0$ ; $U'(0)=1$.



By setting $epsilon = 0 $, we get the boundary layer equation :
$U'' + U' = 0$ with same initial conditions.



The solution of this equation is then $U(tau)= 2-e^{-tau}$.



Finally matching terms:



$displaystyle lim_{tau to infty} U(tau)=displaystyle lim_{t to 0} u_0(t)$



So that $2=c_1+1$ then $c_1=1$.



Therefore the solution will be $u(t)=U(t/sqrt{epsilon}) +u_0(t) - 2= - e^{-t/sqrt{epsilon}} + e^{2t} + e^t $, which corresponds to adding the two solutions found earlier and subtracting the common term.



Is my work correct?
If so, how come my solution satisfies the first initial condition but not the second?



Any suggestions will be greatly appreciate it.
Thanks!










share|cite|improve this question




























    up vote
    0
    down vote

    favorite












    Consider the initial value problem
    $sqrt{epsilon} , u'' + u' - u = e^{2t}$ , with $u(0)=1$, and $u'(0)=1/sqrt{epsilon}$.



    I am trying to use a matched asymptotic expansion to find the leading term of the solution.



    First by setting $epsilon=0$, and using the initial condition $u(0)=1$, we get the outer solution $u_0(t)=c_1 e^t + e^{2t}$.



    I set the change of variables $tau=t/sqrt{epsilon}$.



    Writing the original equation in terms of the new variable I got:



    $U'' + U' -sqrt{epsilon} , U = sqrt{epsilon} , e^{2tau epsilon^{1/2}}$



    with $U(0)=0$ ; $U'(0)=1$.



    By setting $epsilon = 0 $, we get the boundary layer equation :
    $U'' + U' = 0$ with same initial conditions.



    The solution of this equation is then $U(tau)= 2-e^{-tau}$.



    Finally matching terms:



    $displaystyle lim_{tau to infty} U(tau)=displaystyle lim_{t to 0} u_0(t)$



    So that $2=c_1+1$ then $c_1=1$.



    Therefore the solution will be $u(t)=U(t/sqrt{epsilon}) +u_0(t) - 2= - e^{-t/sqrt{epsilon}} + e^{2t} + e^t $, which corresponds to adding the two solutions found earlier and subtracting the common term.



    Is my work correct?
    If so, how come my solution satisfies the first initial condition but not the second?



    Any suggestions will be greatly appreciate it.
    Thanks!










    share|cite|improve this question


























      up vote
      0
      down vote

      favorite









      up vote
      0
      down vote

      favorite











      Consider the initial value problem
      $sqrt{epsilon} , u'' + u' - u = e^{2t}$ , with $u(0)=1$, and $u'(0)=1/sqrt{epsilon}$.



      I am trying to use a matched asymptotic expansion to find the leading term of the solution.



      First by setting $epsilon=0$, and using the initial condition $u(0)=1$, we get the outer solution $u_0(t)=c_1 e^t + e^{2t}$.



      I set the change of variables $tau=t/sqrt{epsilon}$.



      Writing the original equation in terms of the new variable I got:



      $U'' + U' -sqrt{epsilon} , U = sqrt{epsilon} , e^{2tau epsilon^{1/2}}$



      with $U(0)=0$ ; $U'(0)=1$.



      By setting $epsilon = 0 $, we get the boundary layer equation :
      $U'' + U' = 0$ with same initial conditions.



      The solution of this equation is then $U(tau)= 2-e^{-tau}$.



      Finally matching terms:



      $displaystyle lim_{tau to infty} U(tau)=displaystyle lim_{t to 0} u_0(t)$



      So that $2=c_1+1$ then $c_1=1$.



      Therefore the solution will be $u(t)=U(t/sqrt{epsilon}) +u_0(t) - 2= - e^{-t/sqrt{epsilon}} + e^{2t} + e^t $, which corresponds to adding the two solutions found earlier and subtracting the common term.



      Is my work correct?
      If so, how come my solution satisfies the first initial condition but not the second?



      Any suggestions will be greatly appreciate it.
      Thanks!










      share|cite|improve this question















      Consider the initial value problem
      $sqrt{epsilon} , u'' + u' - u = e^{2t}$ , with $u(0)=1$, and $u'(0)=1/sqrt{epsilon}$.



      I am trying to use a matched asymptotic expansion to find the leading term of the solution.



      First by setting $epsilon=0$, and using the initial condition $u(0)=1$, we get the outer solution $u_0(t)=c_1 e^t + e^{2t}$.



      I set the change of variables $tau=t/sqrt{epsilon}$.



      Writing the original equation in terms of the new variable I got:



      $U'' + U' -sqrt{epsilon} , U = sqrt{epsilon} , e^{2tau epsilon^{1/2}}$



      with $U(0)=0$ ; $U'(0)=1$.



      By setting $epsilon = 0 $, we get the boundary layer equation :
      $U'' + U' = 0$ with same initial conditions.



      The solution of this equation is then $U(tau)= 2-e^{-tau}$.



      Finally matching terms:



      $displaystyle lim_{tau to infty} U(tau)=displaystyle lim_{t to 0} u_0(t)$



      So that $2=c_1+1$ then $c_1=1$.



      Therefore the solution will be $u(t)=U(t/sqrt{epsilon}) +u_0(t) - 2= - e^{-t/sqrt{epsilon}} + e^{2t} + e^t $, which corresponds to adding the two solutions found earlier and subtracting the common term.



      Is my work correct?
      If so, how come my solution satisfies the first initial condition but not the second?



      Any suggestions will be greatly appreciate it.
      Thanks!







      differential-equations perturbation-theory






      share|cite|improve this question















      share|cite|improve this question













      share|cite|improve this question




      share|cite|improve this question








      edited Nov 23 at 16:54

























      asked Nov 23 at 3:27









      user569959

      937




      937






















          1 Answer
          1






          active

          oldest

          votes

















          up vote
          0
          down vote













          As the approximation is to the first order in $sqrtϵ$, you can also only expect to satisfy the initial conditions to first order. And indeed for the combined solution $y'(0)=frac1{sqrtϵ}+3=frac1{sqrtϵ}(1+3sqrtϵ)$ has only a first order (relative) error.






          share|cite|improve this answer























          • Thank you for your answer LutzL. I just realized there is a typo in the initial condition for $U$. It should read $U(0)=1$ not $U(0)=0$. In that case I think $U(tau)=2-e^{-tau}$ would be correct.
            – user569959
            Nov 23 at 17:05













          Your Answer





          StackExchange.ifUsing("editor", function () {
          return StackExchange.using("mathjaxEditing", function () {
          StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
          StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
          });
          });
          }, "mathjax-editing");

          StackExchange.ready(function() {
          var channelOptions = {
          tags: "".split(" "),
          id: "69"
          };
          initTagRenderer("".split(" "), "".split(" "), channelOptions);

          StackExchange.using("externalEditor", function() {
          // Have to fire editor after snippets, if snippets enabled
          if (StackExchange.settings.snippets.snippetsEnabled) {
          StackExchange.using("snippets", function() {
          createEditor();
          });
          }
          else {
          createEditor();
          }
          });

          function createEditor() {
          StackExchange.prepareEditor({
          heartbeatType: 'answer',
          convertImagesToLinks: true,
          noModals: true,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          imageUploader: {
          brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
          contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
          allowUrls: true
          },
          noCode: true, onDemand: true,
          discardSelector: ".discard-answer"
          ,immediatelyShowMarkdownHelp:true
          });


          }
          });














          draft saved

          draft discarded


















          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3009942%2fboundary-layer-type-with-initial-value-problem%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown

























          1 Answer
          1






          active

          oldest

          votes








          1 Answer
          1






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes








          up vote
          0
          down vote













          As the approximation is to the first order in $sqrtϵ$, you can also only expect to satisfy the initial conditions to first order. And indeed for the combined solution $y'(0)=frac1{sqrtϵ}+3=frac1{sqrtϵ}(1+3sqrtϵ)$ has only a first order (relative) error.






          share|cite|improve this answer























          • Thank you for your answer LutzL. I just realized there is a typo in the initial condition for $U$. It should read $U(0)=1$ not $U(0)=0$. In that case I think $U(tau)=2-e^{-tau}$ would be correct.
            – user569959
            Nov 23 at 17:05

















          up vote
          0
          down vote













          As the approximation is to the first order in $sqrtϵ$, you can also only expect to satisfy the initial conditions to first order. And indeed for the combined solution $y'(0)=frac1{sqrtϵ}+3=frac1{sqrtϵ}(1+3sqrtϵ)$ has only a first order (relative) error.






          share|cite|improve this answer























          • Thank you for your answer LutzL. I just realized there is a typo in the initial condition for $U$. It should read $U(0)=1$ not $U(0)=0$. In that case I think $U(tau)=2-e^{-tau}$ would be correct.
            – user569959
            Nov 23 at 17:05















          up vote
          0
          down vote










          up vote
          0
          down vote









          As the approximation is to the first order in $sqrtϵ$, you can also only expect to satisfy the initial conditions to first order. And indeed for the combined solution $y'(0)=frac1{sqrtϵ}+3=frac1{sqrtϵ}(1+3sqrtϵ)$ has only a first order (relative) error.






          share|cite|improve this answer














          As the approximation is to the first order in $sqrtϵ$, you can also only expect to satisfy the initial conditions to first order. And indeed for the combined solution $y'(0)=frac1{sqrtϵ}+3=frac1{sqrtϵ}(1+3sqrtϵ)$ has only a first order (relative) error.







          share|cite|improve this answer














          share|cite|improve this answer



          share|cite|improve this answer








          edited Nov 23 at 17:18

























          answered Nov 23 at 8:38









          LutzL

          54.3k41953




          54.3k41953












          • Thank you for your answer LutzL. I just realized there is a typo in the initial condition for $U$. It should read $U(0)=1$ not $U(0)=0$. In that case I think $U(tau)=2-e^{-tau}$ would be correct.
            – user569959
            Nov 23 at 17:05




















          • Thank you for your answer LutzL. I just realized there is a typo in the initial condition for $U$. It should read $U(0)=1$ not $U(0)=0$. In that case I think $U(tau)=2-e^{-tau}$ would be correct.
            – user569959
            Nov 23 at 17:05


















          Thank you for your answer LutzL. I just realized there is a typo in the initial condition for $U$. It should read $U(0)=1$ not $U(0)=0$. In that case I think $U(tau)=2-e^{-tau}$ would be correct.
          – user569959
          Nov 23 at 17:05






          Thank you for your answer LutzL. I just realized there is a typo in the initial condition for $U$. It should read $U(0)=1$ not $U(0)=0$. In that case I think $U(tau)=2-e^{-tau}$ would be correct.
          – user569959
          Nov 23 at 17:05




















          draft saved

          draft discarded




















































          Thanks for contributing an answer to Mathematics Stack Exchange!


          • Please be sure to answer the question. Provide details and share your research!

          But avoid



          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.


          Use MathJax to format equations. MathJax reference.


          To learn more, see our tips on writing great answers.





          Some of your past answers have not been well-received, and you're in danger of being blocked from answering.


          Please pay close attention to the following guidance:


          • Please be sure to answer the question. Provide details and share your research!

          But avoid



          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.


          To learn more, see our tips on writing great answers.




          draft saved


          draft discarded














          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3009942%2fboundary-layer-type-with-initial-value-problem%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown





















































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown

































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown







          Popular posts from this blog

          Wiesbaden

          Marschland

          Dieringhausen