Method of characteristics example: using different change of variables to textbook, characteristic speed, and example clarification

The name of the pictureThe name of the pictureThe name of the pictureClash Royale CLAN TAG#URR8PPP











up vote
1
down vote

favorite
1












The following method of characteristics example is from chapter 4 of Essential Partial Differential Equations by Griffiths, Dold, and Silvester:




Example 4.2



(Half-plane problem) Solve the PDE $pu_x + qu_y = f$ (with $p$ and $q$ constant) in the domain $alpha x + beta y > 0$ given that $u = g(x)$ on the line $mathscrl: alpha x + beta y = 0$, where $q > 0$ and $beta > 0$.



The characteristic equations are readily solved to give



$$x = pt + C_1, y = qt + C_2 (4.8)$$



where $C_1$ and $C_2$ are constant along any characteristic. The assumption $q > 0$ means that $y$ increases along a characteristic as $t$ increases -- thus the parameter $t$ may be viewed as a time-like variable.



In order to determine the solution at a point $P(x, y)$, we trace the characteristic through this point backwards in $t$ until it intersects the line $mathscrl$. We shall suppose that this occurs at the point $Q$ having coordinates $x = beta s$, $y = - alpha s$ when $t = 0$. Then, from (4.8), we find $C_1 = beta s$, $C_2 = -alpha s$. $s$ varies along the line $mathscrl$ and each choice $s$ selects a different characteristic:



$$x = pt + beta s, y = qt - alpha s (4.9)$$



Note that the characteristics have inherited their parameterization from that of the line $mathscrl$. Next, we integrate the ODE $dfracdudt = f$:



$$int_0^t fracdudt' dt' = int_0^t f(x(s, t'), y(s, t')) dt'$$



to give



$$u(s, t) = u(s, 0) + F(t), F(t) = int_0^t f(x(s, t'), y(s, t')) dt'$$



Here, $A(s) = u(s, 0)$ is constant on any characteristic but varies , with $s$, from one characteristic to another. The initial condition gives $u = g(x) = g(beta s)$ at $t = 0$ from which we deduce (using (4.9)) that $u(s, 0) = g(beta s)$ so that



$$u = g(beta s) + int_0^t f(x(s, t'), y(s, t')) dt' (4.10)$$



Finally, to express $u$ as a function of $x$ and $y$, the equation (4.9) have to be solved to give $s$ and $t$ in terms of $x$ and $y$:



$$t = fracalpha x + beta yalpha p + beta q, s = fracqx - pyalpha p + beta q$$



which are valid so long as $alpha p + beta q not= 0$. This condition has a geometric interpretation. Since $alpha p + beta q$ is the scalar product of the vector $(p, q)$ (which is parallel to the characteristics), and the vector $(alpha, beta)$ (which is orthogonal to the line $mathscrl$), the change of variables (4.9) is one-to-one if, and only if, (p, q) is not parallel to $(alpha, beta)$. That is the initial condition should not be specified on any line parallel to the characteristics. This is an important conclusion that applied more widely and, when it holds, we have a unique solution at a point $P$ that depends only on the initial value $g$ given at the foot $Q$ of the characteristic through $P$ and the values of $f$ along $PQ$.



Finally, note that in the simple homogeneous case the coefficients of the PDE $pu_x + qu_y = 0$ form a vector $(p, q)$, whereas the solution (4.10) depends only on $s$, that is, $qx - py$ whose coefficients form a vector $(q, -p)$ that is orthogonal to $(p, q)$. This means that boundary information is simply conveyed along characteristics. Since information travels a horizontal distance $p / q$ for each unit of vertical distance, $p/q$ is known as the characteristic speed of the equation.




  1. When I did the change of variables for this example, I set $x(0) = s$ and $y(0) = frac-alpha sbeta$ so that we get $x(t) = pt + s$ and $y(t) = qt - fracalpha sbeta$. However, when finding the change of variables so that we get $t$ and $s$ in terms of $x$ and $y$, I get the same equations $t = fracalpha x + beta yalpha p + beta q, s = fracqx - pyalpha p + beta q$ as the example. So, despite doing the change of variables differently, is my work still correct?

And, at the end, I don't understand what is meant by this:




Finally, note that in the simple homogeneous case the coefficients of the PDE $pu_x + qu_y = 0$ form a vector $(p, q)$, whereas the solution (4.10) depends only on $s$, that is, $qx - py$ whose coefficients form a vector $(q, -p)$ that is orthogonal to $(p, q)$. This means that boundary information is simply conveyed along characteristics.




  1. First of all, since $u = u(s, t)$ isn't (4.10) dependent on both $s$ and $t$, not just $s$?


  2. Also, wouldn't $s$ be dependent on $fracqx - pyalpha p + beta q$, not just $qx - py$, since $s = fracqx - pyalpha p + beta q$?


  3. And why does the last part mean that boundary information is simply conveyed along characteristics? This isn't clear to me.



Since information travels a horizontal distance $p / q$ for each unit of vertical distance, $p/q$ is known as the characteristic speed of the equation.




  1. Why does information travels a horizontal distance $p / q$ for each unit of vertical distance? I don't see how the author got this.

Thank you for any help you could kindly provide.







share|cite|improve this question

























    up vote
    1
    down vote

    favorite
    1












    The following method of characteristics example is from chapter 4 of Essential Partial Differential Equations by Griffiths, Dold, and Silvester:




    Example 4.2



    (Half-plane problem) Solve the PDE $pu_x + qu_y = f$ (with $p$ and $q$ constant) in the domain $alpha x + beta y > 0$ given that $u = g(x)$ on the line $mathscrl: alpha x + beta y = 0$, where $q > 0$ and $beta > 0$.



    The characteristic equations are readily solved to give



    $$x = pt + C_1, y = qt + C_2 (4.8)$$



    where $C_1$ and $C_2$ are constant along any characteristic. The assumption $q > 0$ means that $y$ increases along a characteristic as $t$ increases -- thus the parameter $t$ may be viewed as a time-like variable.



    In order to determine the solution at a point $P(x, y)$, we trace the characteristic through this point backwards in $t$ until it intersects the line $mathscrl$. We shall suppose that this occurs at the point $Q$ having coordinates $x = beta s$, $y = - alpha s$ when $t = 0$. Then, from (4.8), we find $C_1 = beta s$, $C_2 = -alpha s$. $s$ varies along the line $mathscrl$ and each choice $s$ selects a different characteristic:



    $$x = pt + beta s, y = qt - alpha s (4.9)$$



    Note that the characteristics have inherited their parameterization from that of the line $mathscrl$. Next, we integrate the ODE $dfracdudt = f$:



    $$int_0^t fracdudt' dt' = int_0^t f(x(s, t'), y(s, t')) dt'$$



    to give



    $$u(s, t) = u(s, 0) + F(t), F(t) = int_0^t f(x(s, t'), y(s, t')) dt'$$



    Here, $A(s) = u(s, 0)$ is constant on any characteristic but varies , with $s$, from one characteristic to another. The initial condition gives $u = g(x) = g(beta s)$ at $t = 0$ from which we deduce (using (4.9)) that $u(s, 0) = g(beta s)$ so that



    $$u = g(beta s) + int_0^t f(x(s, t'), y(s, t')) dt' (4.10)$$



    Finally, to express $u$ as a function of $x$ and $y$, the equation (4.9) have to be solved to give $s$ and $t$ in terms of $x$ and $y$:



    $$t = fracalpha x + beta yalpha p + beta q, s = fracqx - pyalpha p + beta q$$



    which are valid so long as $alpha p + beta q not= 0$. This condition has a geometric interpretation. Since $alpha p + beta q$ is the scalar product of the vector $(p, q)$ (which is parallel to the characteristics), and the vector $(alpha, beta)$ (which is orthogonal to the line $mathscrl$), the change of variables (4.9) is one-to-one if, and only if, (p, q) is not parallel to $(alpha, beta)$. That is the initial condition should not be specified on any line parallel to the characteristics. This is an important conclusion that applied more widely and, when it holds, we have a unique solution at a point $P$ that depends only on the initial value $g$ given at the foot $Q$ of the characteristic through $P$ and the values of $f$ along $PQ$.



    Finally, note that in the simple homogeneous case the coefficients of the PDE $pu_x + qu_y = 0$ form a vector $(p, q)$, whereas the solution (4.10) depends only on $s$, that is, $qx - py$ whose coefficients form a vector $(q, -p)$ that is orthogonal to $(p, q)$. This means that boundary information is simply conveyed along characteristics. Since information travels a horizontal distance $p / q$ for each unit of vertical distance, $p/q$ is known as the characteristic speed of the equation.




    1. When I did the change of variables for this example, I set $x(0) = s$ and $y(0) = frac-alpha sbeta$ so that we get $x(t) = pt + s$ and $y(t) = qt - fracalpha sbeta$. However, when finding the change of variables so that we get $t$ and $s$ in terms of $x$ and $y$, I get the same equations $t = fracalpha x + beta yalpha p + beta q, s = fracqx - pyalpha p + beta q$ as the example. So, despite doing the change of variables differently, is my work still correct?

    And, at the end, I don't understand what is meant by this:




    Finally, note that in the simple homogeneous case the coefficients of the PDE $pu_x + qu_y = 0$ form a vector $(p, q)$, whereas the solution (4.10) depends only on $s$, that is, $qx - py$ whose coefficients form a vector $(q, -p)$ that is orthogonal to $(p, q)$. This means that boundary information is simply conveyed along characteristics.




    1. First of all, since $u = u(s, t)$ isn't (4.10) dependent on both $s$ and $t$, not just $s$?


    2. Also, wouldn't $s$ be dependent on $fracqx - pyalpha p + beta q$, not just $qx - py$, since $s = fracqx - pyalpha p + beta q$?


    3. And why does the last part mean that boundary information is simply conveyed along characteristics? This isn't clear to me.



    Since information travels a horizontal distance $p / q$ for each unit of vertical distance, $p/q$ is known as the characteristic speed of the equation.




    1. Why does information travels a horizontal distance $p / q$ for each unit of vertical distance? I don't see how the author got this.

    Thank you for any help you could kindly provide.







    share|cite|improve this question























      up vote
      1
      down vote

      favorite
      1









      up vote
      1
      down vote

      favorite
      1






      1





      The following method of characteristics example is from chapter 4 of Essential Partial Differential Equations by Griffiths, Dold, and Silvester:




      Example 4.2



      (Half-plane problem) Solve the PDE $pu_x + qu_y = f$ (with $p$ and $q$ constant) in the domain $alpha x + beta y > 0$ given that $u = g(x)$ on the line $mathscrl: alpha x + beta y = 0$, where $q > 0$ and $beta > 0$.



      The characteristic equations are readily solved to give



      $$x = pt + C_1, y = qt + C_2 (4.8)$$



      where $C_1$ and $C_2$ are constant along any characteristic. The assumption $q > 0$ means that $y$ increases along a characteristic as $t$ increases -- thus the parameter $t$ may be viewed as a time-like variable.



      In order to determine the solution at a point $P(x, y)$, we trace the characteristic through this point backwards in $t$ until it intersects the line $mathscrl$. We shall suppose that this occurs at the point $Q$ having coordinates $x = beta s$, $y = - alpha s$ when $t = 0$. Then, from (4.8), we find $C_1 = beta s$, $C_2 = -alpha s$. $s$ varies along the line $mathscrl$ and each choice $s$ selects a different characteristic:



      $$x = pt + beta s, y = qt - alpha s (4.9)$$



      Note that the characteristics have inherited their parameterization from that of the line $mathscrl$. Next, we integrate the ODE $dfracdudt = f$:



      $$int_0^t fracdudt' dt' = int_0^t f(x(s, t'), y(s, t')) dt'$$



      to give



      $$u(s, t) = u(s, 0) + F(t), F(t) = int_0^t f(x(s, t'), y(s, t')) dt'$$



      Here, $A(s) = u(s, 0)$ is constant on any characteristic but varies , with $s$, from one characteristic to another. The initial condition gives $u = g(x) = g(beta s)$ at $t = 0$ from which we deduce (using (4.9)) that $u(s, 0) = g(beta s)$ so that



      $$u = g(beta s) + int_0^t f(x(s, t'), y(s, t')) dt' (4.10)$$



      Finally, to express $u$ as a function of $x$ and $y$, the equation (4.9) have to be solved to give $s$ and $t$ in terms of $x$ and $y$:



      $$t = fracalpha x + beta yalpha p + beta q, s = fracqx - pyalpha p + beta q$$



      which are valid so long as $alpha p + beta q not= 0$. This condition has a geometric interpretation. Since $alpha p + beta q$ is the scalar product of the vector $(p, q)$ (which is parallel to the characteristics), and the vector $(alpha, beta)$ (which is orthogonal to the line $mathscrl$), the change of variables (4.9) is one-to-one if, and only if, (p, q) is not parallel to $(alpha, beta)$. That is the initial condition should not be specified on any line parallel to the characteristics. This is an important conclusion that applied more widely and, when it holds, we have a unique solution at a point $P$ that depends only on the initial value $g$ given at the foot $Q$ of the characteristic through $P$ and the values of $f$ along $PQ$.



      Finally, note that in the simple homogeneous case the coefficients of the PDE $pu_x + qu_y = 0$ form a vector $(p, q)$, whereas the solution (4.10) depends only on $s$, that is, $qx - py$ whose coefficients form a vector $(q, -p)$ that is orthogonal to $(p, q)$. This means that boundary information is simply conveyed along characteristics. Since information travels a horizontal distance $p / q$ for each unit of vertical distance, $p/q$ is known as the characteristic speed of the equation.




      1. When I did the change of variables for this example, I set $x(0) = s$ and $y(0) = frac-alpha sbeta$ so that we get $x(t) = pt + s$ and $y(t) = qt - fracalpha sbeta$. However, when finding the change of variables so that we get $t$ and $s$ in terms of $x$ and $y$, I get the same equations $t = fracalpha x + beta yalpha p + beta q, s = fracqx - pyalpha p + beta q$ as the example. So, despite doing the change of variables differently, is my work still correct?

      And, at the end, I don't understand what is meant by this:




      Finally, note that in the simple homogeneous case the coefficients of the PDE $pu_x + qu_y = 0$ form a vector $(p, q)$, whereas the solution (4.10) depends only on $s$, that is, $qx - py$ whose coefficients form a vector $(q, -p)$ that is orthogonal to $(p, q)$. This means that boundary information is simply conveyed along characteristics.




      1. First of all, since $u = u(s, t)$ isn't (4.10) dependent on both $s$ and $t$, not just $s$?


      2. Also, wouldn't $s$ be dependent on $fracqx - pyalpha p + beta q$, not just $qx - py$, since $s = fracqx - pyalpha p + beta q$?


      3. And why does the last part mean that boundary information is simply conveyed along characteristics? This isn't clear to me.



      Since information travels a horizontal distance $p / q$ for each unit of vertical distance, $p/q$ is known as the characteristic speed of the equation.




      1. Why does information travels a horizontal distance $p / q$ for each unit of vertical distance? I don't see how the author got this.

      Thank you for any help you could kindly provide.







      share|cite|improve this question













      The following method of characteristics example is from chapter 4 of Essential Partial Differential Equations by Griffiths, Dold, and Silvester:




      Example 4.2



      (Half-plane problem) Solve the PDE $pu_x + qu_y = f$ (with $p$ and $q$ constant) in the domain $alpha x + beta y > 0$ given that $u = g(x)$ on the line $mathscrl: alpha x + beta y = 0$, where $q > 0$ and $beta > 0$.



      The characteristic equations are readily solved to give



      $$x = pt + C_1, y = qt + C_2 (4.8)$$



      where $C_1$ and $C_2$ are constant along any characteristic. The assumption $q > 0$ means that $y$ increases along a characteristic as $t$ increases -- thus the parameter $t$ may be viewed as a time-like variable.



      In order to determine the solution at a point $P(x, y)$, we trace the characteristic through this point backwards in $t$ until it intersects the line $mathscrl$. We shall suppose that this occurs at the point $Q$ having coordinates $x = beta s$, $y = - alpha s$ when $t = 0$. Then, from (4.8), we find $C_1 = beta s$, $C_2 = -alpha s$. $s$ varies along the line $mathscrl$ and each choice $s$ selects a different characteristic:



      $$x = pt + beta s, y = qt - alpha s (4.9)$$



      Note that the characteristics have inherited their parameterization from that of the line $mathscrl$. Next, we integrate the ODE $dfracdudt = f$:



      $$int_0^t fracdudt' dt' = int_0^t f(x(s, t'), y(s, t')) dt'$$



      to give



      $$u(s, t) = u(s, 0) + F(t), F(t) = int_0^t f(x(s, t'), y(s, t')) dt'$$



      Here, $A(s) = u(s, 0)$ is constant on any characteristic but varies , with $s$, from one characteristic to another. The initial condition gives $u = g(x) = g(beta s)$ at $t = 0$ from which we deduce (using (4.9)) that $u(s, 0) = g(beta s)$ so that



      $$u = g(beta s) + int_0^t f(x(s, t'), y(s, t')) dt' (4.10)$$



      Finally, to express $u$ as a function of $x$ and $y$, the equation (4.9) have to be solved to give $s$ and $t$ in terms of $x$ and $y$:



      $$t = fracalpha x + beta yalpha p + beta q, s = fracqx - pyalpha p + beta q$$



      which are valid so long as $alpha p + beta q not= 0$. This condition has a geometric interpretation. Since $alpha p + beta q$ is the scalar product of the vector $(p, q)$ (which is parallel to the characteristics), and the vector $(alpha, beta)$ (which is orthogonal to the line $mathscrl$), the change of variables (4.9) is one-to-one if, and only if, (p, q) is not parallel to $(alpha, beta)$. That is the initial condition should not be specified on any line parallel to the characteristics. This is an important conclusion that applied more widely and, when it holds, we have a unique solution at a point $P$ that depends only on the initial value $g$ given at the foot $Q$ of the characteristic through $P$ and the values of $f$ along $PQ$.



      Finally, note that in the simple homogeneous case the coefficients of the PDE $pu_x + qu_y = 0$ form a vector $(p, q)$, whereas the solution (4.10) depends only on $s$, that is, $qx - py$ whose coefficients form a vector $(q, -p)$ that is orthogonal to $(p, q)$. This means that boundary information is simply conveyed along characteristics. Since information travels a horizontal distance $p / q$ for each unit of vertical distance, $p/q$ is known as the characteristic speed of the equation.




      1. When I did the change of variables for this example, I set $x(0) = s$ and $y(0) = frac-alpha sbeta$ so that we get $x(t) = pt + s$ and $y(t) = qt - fracalpha sbeta$. However, when finding the change of variables so that we get $t$ and $s$ in terms of $x$ and $y$, I get the same equations $t = fracalpha x + beta yalpha p + beta q, s = fracqx - pyalpha p + beta q$ as the example. So, despite doing the change of variables differently, is my work still correct?

      And, at the end, I don't understand what is meant by this:




      Finally, note that in the simple homogeneous case the coefficients of the PDE $pu_x + qu_y = 0$ form a vector $(p, q)$, whereas the solution (4.10) depends only on $s$, that is, $qx - py$ whose coefficients form a vector $(q, -p)$ that is orthogonal to $(p, q)$. This means that boundary information is simply conveyed along characteristics.




      1. First of all, since $u = u(s, t)$ isn't (4.10) dependent on both $s$ and $t$, not just $s$?


      2. Also, wouldn't $s$ be dependent on $fracqx - pyalpha p + beta q$, not just $qx - py$, since $s = fracqx - pyalpha p + beta q$?


      3. And why does the last part mean that boundary information is simply conveyed along characteristics? This isn't clear to me.



      Since information travels a horizontal distance $p / q$ for each unit of vertical distance, $p/q$ is known as the characteristic speed of the equation.




      1. Why does information travels a horizontal distance $p / q$ for each unit of vertical distance? I don't see how the author got this.

      Thank you for any help you could kindly provide.









      share|cite|improve this question












      share|cite|improve this question




      share|cite|improve this question








      edited Aug 7 at 10:22
























      asked Aug 5 at 23:03









      handler's handle

      1508




      1508

























          active

          oldest

          votes











          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: false,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          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%2f2873412%2fmethod-of-characteristics-example-using-different-change-of-variables-to-textbo%23new-answer', 'question_page');

          );

          Post as a guest



































          active

          oldest

          votes













          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes










           

          draft saved


          draft discarded


























           


          draft saved


          draft discarded














          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2873412%2fmethod-of-characteristics-example-using-different-change-of-variables-to-textbo%23new-answer', 'question_page');

          );

          Post as a guest













































































          Comments

          Popular posts from this blog

          What is the equation of a 3D cone with generalised tilt?

          Color the edges and diagonals of a regular polygon

          Relationship between determinant of matrix and determinant of adjoint?