Analytical solution for PDE-system's IBVP to validate method of lines












0












$begingroup$


I am in search of initial-boundary value problems which are posed in the form of a system of coupled PDE's (not a single PDE) and for which an analytical solution exists at least for some specific initial-boundary conditions.



In particular I would be interested in a PDE's system of one dynamical ($partial_r$ on lhs) equation coupled to one constraint ($partial_t$ on lhs) equation.



The generic character of the question is due to the fact that l have not yet come acrross any PDE's system with the above characteristics. Being more specific would probably make the question overdetermined.



The ultimate purpose of my search is to use the analytical solutions as a numerical test for Mathematica Method of Lines code.



I would appreciate any recommendations.






pde numerical-methods boundary-value-problem







share|cite|improve this question











$endgroup$












  • $begingroup$
    The heat equation has various benchmark problems with known solutions.
    $endgroup$
    – User123456789
    Dec 20 '18 at 10:48










  • $begingroup$
    This question is an enhanced version of a question I deleted. The deleted question did not lead to any usefull results except the comment by Skip linking to this PDE's system.
    $endgroup$
    – dkstack
    Dec 20 '18 at 10:49










  • $begingroup$
    @User123456789 these benchmark problems concern systems of coupled PDE's or single equations?
    $endgroup$
    – dkstack
    Dec 20 '18 at 10:52










  • $begingroup$
    You are right, I misread. You can maybe think of Stokes' equation then
    $endgroup$
    – User123456789
    Dec 20 '18 at 10:55










  • $begingroup$
    @User123456789 I searched in wikipedia about Stoke's equation but found nothing about PDE's system with an exact solution. Is there any link about the Stoke's equation you refer to?
    $endgroup$
    – dkstack
    Dec 20 '18 at 11:12


















0












$begingroup$


I am in search of initial-boundary value problems which are posed in the form of a system of coupled PDE's (not a single PDE) and for which an analytical solution exists at least for some specific initial-boundary conditions.



In particular I would be interested in a PDE's system of one dynamical ($partial_r$ on lhs) equation coupled to one constraint ($partial_t$ on lhs) equation.



The generic character of the question is due to the fact that l have not yet come acrross any PDE's system with the above characteristics. Being more specific would probably make the question overdetermined.



The ultimate purpose of my search is to use the analytical solutions as a numerical test for Mathematica Method of Lines code.



I would appreciate any recommendations.






pde numerical-methods boundary-value-problem







share|cite|improve this question











$endgroup$












  • $begingroup$
    The heat equation has various benchmark problems with known solutions.
    $endgroup$
    – User123456789
    Dec 20 '18 at 10:48










  • $begingroup$
    This question is an enhanced version of a question I deleted. The deleted question did not lead to any usefull results except the comment by Skip linking to this PDE's system.
    $endgroup$
    – dkstack
    Dec 20 '18 at 10:49










  • $begingroup$
    @User123456789 these benchmark problems concern systems of coupled PDE's or single equations?
    $endgroup$
    – dkstack
    Dec 20 '18 at 10:52










  • $begingroup$
    You are right, I misread. You can maybe think of Stokes' equation then
    $endgroup$
    – User123456789
    Dec 20 '18 at 10:55










  • $begingroup$
    @User123456789 I searched in wikipedia about Stoke's equation but found nothing about PDE's system with an exact solution. Is there any link about the Stoke's equation you refer to?
    $endgroup$
    – dkstack
    Dec 20 '18 at 11:12
















0












0








0


0



$begingroup$


I am in search of initial-boundary value problems which are posed in the form of a system of coupled PDE's (not a single PDE) and for which an analytical solution exists at least for some specific initial-boundary conditions.



In particular I would be interested in a PDE's system of one dynamical ($partial_r$ on lhs) equation coupled to one constraint ($partial_t$ on lhs) equation.



The generic character of the question is due to the fact that l have not yet come acrross any PDE's system with the above characteristics. Being more specific would probably make the question overdetermined.



The ultimate purpose of my search is to use the analytical solutions as a numerical test for Mathematica Method of Lines code.



I would appreciate any recommendations.






pde numerical-methods boundary-value-problem







share|cite|improve this question











$endgroup$




I am in search of initial-boundary value problems which are posed in the form of a system of coupled PDE's (not a single PDE) and for which an analytical solution exists at least for some specific initial-boundary conditions.



In particular I would be interested in a PDE's system of one dynamical ($partial_r$ on lhs) equation coupled to one constraint ($partial_t$ on lhs) equation.



The generic character of the question is due to the fact that l have not yet come acrross any PDE's system with the above characteristics. Being more specific would probably make the question overdetermined.



The ultimate purpose of my search is to use the analytical solutions as a numerical test for Mathematica Method of Lines code.



I would appreciate any recommendations.






pde numerical-methods boundary-value-problem




pde numerical-methods boundary-value-problem






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Dec 21 '18 at 11:21







dkstack

















asked Dec 20 '18 at 10:46









dkstackdkstack

204




204












  • $begingroup$
    The heat equation has various benchmark problems with known solutions.
    $endgroup$
    – User123456789
    Dec 20 '18 at 10:48










  • $begingroup$
    This question is an enhanced version of a question I deleted. The deleted question did not lead to any usefull results except the comment by Skip linking to this PDE's system.
    $endgroup$
    – dkstack
    Dec 20 '18 at 10:49










  • $begingroup$
    @User123456789 these benchmark problems concern systems of coupled PDE's or single equations?
    $endgroup$
    – dkstack
    Dec 20 '18 at 10:52










  • $begingroup$
    You are right, I misread. You can maybe think of Stokes' equation then
    $endgroup$
    – User123456789
    Dec 20 '18 at 10:55










  • $begingroup$
    @User123456789 I searched in wikipedia about Stoke's equation but found nothing about PDE's system with an exact solution. Is there any link about the Stoke's equation you refer to?
    $endgroup$
    – dkstack
    Dec 20 '18 at 11:12




















  • $begingroup$
    The heat equation has various benchmark problems with known solutions.
    $endgroup$
    – User123456789
    Dec 20 '18 at 10:48










  • $begingroup$
    This question is an enhanced version of a question I deleted. The deleted question did not lead to any usefull results except the comment by Skip linking to this PDE's system.
    $endgroup$
    – dkstack
    Dec 20 '18 at 10:49










  • $begingroup$
    @User123456789 these benchmark problems concern systems of coupled PDE's or single equations?
    $endgroup$
    – dkstack
    Dec 20 '18 at 10:52










  • $begingroup$
    You are right, I misread. You can maybe think of Stokes' equation then
    $endgroup$
    – User123456789
    Dec 20 '18 at 10:55










  • $begingroup$
    @User123456789 I searched in wikipedia about Stoke's equation but found nothing about PDE's system with an exact solution. Is there any link about the Stoke's equation you refer to?
    $endgroup$
    – dkstack
    Dec 20 '18 at 11:12


















$begingroup$
The heat equation has various benchmark problems with known solutions.
$endgroup$
– User123456789
Dec 20 '18 at 10:48




$begingroup$
The heat equation has various benchmark problems with known solutions.
$endgroup$
– User123456789
Dec 20 '18 at 10:48












$begingroup$
This question is an enhanced version of a question I deleted. The deleted question did not lead to any usefull results except the comment by Skip linking to this PDE's system.
$endgroup$
– dkstack
Dec 20 '18 at 10:49




$begingroup$
This question is an enhanced version of a question I deleted. The deleted question did not lead to any usefull results except the comment by Skip linking to this PDE's system.
$endgroup$
– dkstack
Dec 20 '18 at 10:49












$begingroup$
@User123456789 these benchmark problems concern systems of coupled PDE's or single equations?
$endgroup$
– dkstack
Dec 20 '18 at 10:52




$begingroup$
@User123456789 these benchmark problems concern systems of coupled PDE's or single equations?
$endgroup$
– dkstack
Dec 20 '18 at 10:52












$begingroup$
You are right, I misread. You can maybe think of Stokes' equation then
$endgroup$
– User123456789
Dec 20 '18 at 10:55




$begingroup$
You are right, I misread. You can maybe think of Stokes' equation then
$endgroup$
– User123456789
Dec 20 '18 at 10:55












$begingroup$
@User123456789 I searched in wikipedia about Stoke's equation but found nothing about PDE's system with an exact solution. Is there any link about the Stoke's equation you refer to?
$endgroup$
– dkstack
Dec 20 '18 at 11:12






$begingroup$
@User123456789 I searched in wikipedia about Stoke's equation but found nothing about PDE's system with an exact solution. Is there any link about the Stoke's equation you refer to?
$endgroup$
– dkstack
Dec 20 '18 at 11:12












2 Answers
2






active

oldest

votes


















1












$begingroup$

You could use the "method of manufactured solutions (MMS)", which is nicely described in



https://www.comsol.de/blogs/verify-simulations-with-the-method-of-manufactured-solutions/



Essentially, you take some system of PDE, summarized as $$F[u]=r_F,$$
where $[u]$ includes the partial derivatives of $u$ that are used, and right side $r_B$ which is to be chosen later. Denote the boundary conditions similarly as $$B[u]=r_B.$$ Now select some function $p$ that is in the class of admissible functions and perhaps contains qualities expected in a solution.



Then the test system is $$F[u]=F[p],~~ B[u]=B[p],$$ that is, the right sides are computed as $r_F=F[p]$ and $r_B=B[p]$, and $p$ is obviously the exact solution.




  • discussion of MMS in a very simple context, with accordingly simple example https://scicomp.stackexchange.com/q/30562/6839






share|cite|improve this answer











$endgroup$





















    1












    $begingroup$

    The simplest may be a 1D linear hyperbolic system of conservation laws $q_t + A q_x = 0$, such as the system of linear acoustics
    begin{aligned}
    p_t + u_0 p_x + K_0 u_x &= 0,\
    u_t + p_x/rho_0 + u_0 u_x &= 0,
    end{aligned}
    the system of linear elasticity, or the system of electromagnetism. In all cases, $q in Bbb R^2$ and $A$ is a matrix with real eigenvalues. Almost every initial-and-boundary-value problem for such systems can be solved via the method of characteristics.






    share|cite|improve this answer











    $endgroup$













      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',
      autoActivateHeartbeat: false,
      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%2f3047405%2fanalytical-solution-for-pde-systems-ibvp-to-validate-method-of-lines%23new-answer', 'question_page');
      }
      );

      Post as a guest















      Required, but never shown

























      2 Answers
      2






      active

      oldest

      votes








      2 Answers
      2






      active

      oldest

      votes









      active

      oldest

      votes






      active

      oldest

      votes









      1












      $begingroup$

      You could use the "method of manufactured solutions (MMS)", which is nicely described in



      https://www.comsol.de/blogs/verify-simulations-with-the-method-of-manufactured-solutions/



      Essentially, you take some system of PDE, summarized as $$F[u]=r_F,$$
      where $[u]$ includes the partial derivatives of $u$ that are used, and right side $r_B$ which is to be chosen later. Denote the boundary conditions similarly as $$B[u]=r_B.$$ Now select some function $p$ that is in the class of admissible functions and perhaps contains qualities expected in a solution.



      Then the test system is $$F[u]=F[p],~~ B[u]=B[p],$$ that is, the right sides are computed as $r_F=F[p]$ and $r_B=B[p]$, and $p$ is obviously the exact solution.




      • discussion of MMS in a very simple context, with accordingly simple example https://scicomp.stackexchange.com/q/30562/6839






      share|cite|improve this answer











      $endgroup$


















        1












        $begingroup$

        You could use the "method of manufactured solutions (MMS)", which is nicely described in



        https://www.comsol.de/blogs/verify-simulations-with-the-method-of-manufactured-solutions/



        Essentially, you take some system of PDE, summarized as $$F[u]=r_F,$$
        where $[u]$ includes the partial derivatives of $u$ that are used, and right side $r_B$ which is to be chosen later. Denote the boundary conditions similarly as $$B[u]=r_B.$$ Now select some function $p$ that is in the class of admissible functions and perhaps contains qualities expected in a solution.



        Then the test system is $$F[u]=F[p],~~ B[u]=B[p],$$ that is, the right sides are computed as $r_F=F[p]$ and $r_B=B[p]$, and $p$ is obviously the exact solution.




        • discussion of MMS in a very simple context, with accordingly simple example https://scicomp.stackexchange.com/q/30562/6839






        share|cite|improve this answer











        $endgroup$
















          1












          1








          1





          $begingroup$

          You could use the "method of manufactured solutions (MMS)", which is nicely described in



          https://www.comsol.de/blogs/verify-simulations-with-the-method-of-manufactured-solutions/



          Essentially, you take some system of PDE, summarized as $$F[u]=r_F,$$
          where $[u]$ includes the partial derivatives of $u$ that are used, and right side $r_B$ which is to be chosen later. Denote the boundary conditions similarly as $$B[u]=r_B.$$ Now select some function $p$ that is in the class of admissible functions and perhaps contains qualities expected in a solution.



          Then the test system is $$F[u]=F[p],~~ B[u]=B[p],$$ that is, the right sides are computed as $r_F=F[p]$ and $r_B=B[p]$, and $p$ is obviously the exact solution.




          • discussion of MMS in a very simple context, with accordingly simple example https://scicomp.stackexchange.com/q/30562/6839






          share|cite|improve this answer











          $endgroup$



          You could use the "method of manufactured solutions (MMS)", which is nicely described in



          https://www.comsol.de/blogs/verify-simulations-with-the-method-of-manufactured-solutions/



          Essentially, you take some system of PDE, summarized as $$F[u]=r_F,$$
          where $[u]$ includes the partial derivatives of $u$ that are used, and right side $r_B$ which is to be chosen later. Denote the boundary conditions similarly as $$B[u]=r_B.$$ Now select some function $p$ that is in the class of admissible functions and perhaps contains qualities expected in a solution.



          Then the test system is $$F[u]=F[p],~~ B[u]=B[p],$$ that is, the right sides are computed as $r_F=F[p]$ and $r_B=B[p]$, and $p$ is obviously the exact solution.




          • discussion of MMS in a very simple context, with accordingly simple example https://scicomp.stackexchange.com/q/30562/6839







          share|cite|improve this answer














          share|cite|improve this answer



          share|cite|improve this answer








          edited Dec 20 '18 at 21:54

























          answered Dec 20 '18 at 21:24









          LutzLLutzL

          58.7k42055




          58.7k42055























              1












              $begingroup$

              The simplest may be a 1D linear hyperbolic system of conservation laws $q_t + A q_x = 0$, such as the system of linear acoustics
              begin{aligned}
              p_t + u_0 p_x + K_0 u_x &= 0,\
              u_t + p_x/rho_0 + u_0 u_x &= 0,
              end{aligned}
              the system of linear elasticity, or the system of electromagnetism. In all cases, $q in Bbb R^2$ and $A$ is a matrix with real eigenvalues. Almost every initial-and-boundary-value problem for such systems can be solved via the method of characteristics.






              share|cite|improve this answer











              $endgroup$


















                1












                $begingroup$

                The simplest may be a 1D linear hyperbolic system of conservation laws $q_t + A q_x = 0$, such as the system of linear acoustics
                begin{aligned}
                p_t + u_0 p_x + K_0 u_x &= 0,\
                u_t + p_x/rho_0 + u_0 u_x &= 0,
                end{aligned}
                the system of linear elasticity, or the system of electromagnetism. In all cases, $q in Bbb R^2$ and $A$ is a matrix with real eigenvalues. Almost every initial-and-boundary-value problem for such systems can be solved via the method of characteristics.






                share|cite|improve this answer











                $endgroup$
















                  1












                  1








                  1





                  $begingroup$

                  The simplest may be a 1D linear hyperbolic system of conservation laws $q_t + A q_x = 0$, such as the system of linear acoustics
                  begin{aligned}
                  p_t + u_0 p_x + K_0 u_x &= 0,\
                  u_t + p_x/rho_0 + u_0 u_x &= 0,
                  end{aligned}
                  the system of linear elasticity, or the system of electromagnetism. In all cases, $q in Bbb R^2$ and $A$ is a matrix with real eigenvalues. Almost every initial-and-boundary-value problem for such systems can be solved via the method of characteristics.






                  share|cite|improve this answer











                  $endgroup$



                  The simplest may be a 1D linear hyperbolic system of conservation laws $q_t + A q_x = 0$, such as the system of linear acoustics
                  begin{aligned}
                  p_t + u_0 p_x + K_0 u_x &= 0,\
                  u_t + p_x/rho_0 + u_0 u_x &= 0,
                  end{aligned}
                  the system of linear elasticity, or the system of electromagnetism. In all cases, $q in Bbb R^2$ and $A$ is a matrix with real eigenvalues. Almost every initial-and-boundary-value problem for such systems can be solved via the method of characteristics.







                  share|cite|improve this answer














                  share|cite|improve this answer



                  share|cite|improve this answer








                  edited Dec 20 '18 at 13:58

























                  answered Dec 20 '18 at 13:48









                  Harry49Harry49

                  6,97631238




                  6,97631238






























                      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.




                      draft saved


                      draft discarded














                      StackExchange.ready(
                      function () {
                      StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3047405%2fanalytical-solution-for-pde-systems-ibvp-to-validate-method-of-lines%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

                      Ellipse (mathématiques)

                      Image
                      Pour les articles homonymes, voir Ellipse. L'ellipse est le lieu des points dont la somme des distances à deux points fixes, dits foyers, est constante. Section du cône ou projection du cercle. En géométrie, une ellipse est une courbe plane fermée obtenue par l’intersection d’un cône de révolution avec un plan, à condition que celui-ci coupe l'axe de rotation du cône ou du cylindre : c'est une conique d'excentricité strictement comprise entre 0 et 1. On peut également la définir comme le lieu des points dont la somme des distances à deux points fixes, dits foyers, est constante (sa construction par la méthode du jardinier est très simple). Dans la vie courante, l’ellipse est la forme qu'on perçoit en regardant un cercle en perspective, ou la figure formée par l’ombre d'un disque sur une surface plane. On retrouve aussi, en première approximation [ 1 ] , des ellipses dans les trajectoires des corps célestes (planètes, comètes ou satellites arti...
                      Read more

                      Quarter-circle Tiles

                      Image
                      up vote 0 down vote favorite I have $99$ identical square tiles, each with a quarter-circle drawn on it like this: [asy] size(1.5cm); draw(Arc((2,0),1,90,180),red+1); draw((0,0)--(2,0)--(2,2)--(0,2)--(0,0)); [/asy] When I arrange the tiles in a $9times 11$ rectangular grid, each with a random orientation, what is the expected value of the number of full circles I form? I think this problem has to do with finding the chance any given 2x2 square has a circle, but I can't find it. expected-value share | cite | improve this question asked Nov 20 at 15:03 6minecraftninja 1 2 ...
                      Read more

                      Mont Emei

                      Image
                      Mont Emei.mw-parser-output .entete.map{background-image:url("//upload.wikimedia.org/wikipedia/commons/7/7a/Picto_infobox_map.png")} Vue du sommet Wanfo du mont Emei Géographie Altitude 3 099  m , Wanfo Coordonnées 29° 31′ 11″ nord, 103° 19′ 57″ est Administration Pays   Chine Province Sichuan Ville-préfecture Leshan Géolocalisation sur la carte : Sichuan Mont Emei Géolocalisation sur la carte : Chine Mont Emei modifier   Paysage panoramique du mont Emei, incluant le paysage panoramique du grand Bouddha de Leshan *   Patrimoine mondial de l'UNESCO Passerelle en bois sur le versant ouest Pays   Chine Subdivision Sichuan Type Mixte Critères (iv) (vi) (x) Superficie 15 400 ha Numéro d’identification 779 Zone géographique Asie et Pacifique  ** Année d’inscription 1996 (20 e session) * Descriptif officiel UNESCO ** Classification géog...
                      Read more