Uniqueness of a vector decomposition











up vote
0
down vote

favorite












I want to express a divergence-free vector $vec{w}(x,y)$ with non-zero curl in the domain $x,y geq 0$ as
$$ vec{w} = b nabla c ; .$$
The vector is normal to the boundary $partial Omega$. Furthermore, the value of $c$ on the boundary is zero ($c=0$ on $x,y=0$). Based on a helpful answer I received to this question in another post, I have made some progress. However, I am now unsure about the uniqueness of the problem.



To illustrate the problem, I will assume the model
$$ vec{w} = frac{1}{sqrt{x^2+y^2}} left( begin{array}{c} y \ -x end{array} right) ; .$$
This vector has the properties that it is divergence-free, has non-zero curl, and is normal to the boundary as mentioned above. From the theory of first-order linear PDE's it follows that $c$ has the general solution
$$ c = phi left( frac{x}{y} right) $$
where $phi$ is an arbitrary function of the variable $eta equiv x/y$. From this follows the solution for $b$ as
$$ b = frac{y^2}{sqrt{x^2+y^2}} frac{1}{phi^{prime}} $$
where the prime indicates a derivative. Taking the gradient results in
$$ nabla b = left( begin{array}{c}
-frac{eta}{(1+eta^2)^{3/2}} frac{1}{phi^{prime}} - frac{phi^{primeprime}}{{phi^{prime}}^2} frac{1}{sqrt{1+eta^2}} \
frac{2}{sqrt{1+eta^2}} frac{1}{phi^{prime}} - frac{1}{(1+eta^2)^{3/2}} frac{1}{phi^{prime}} + frac{phi^{primeprime}}{{phi^{prime}}^2} frac{eta}{sqrt{1+eta^2}}
end{array} right) ; . $$

These solutions satisfy the rotational and the divergence-free condition and the equation $vec{w} = b nabla c$. In order to impose the boundary conditions for $c$, one needs $phi(0) = 0$ and $phi(eta) = 0$ as $eta rightarrow infty$, as well as that $phi^{primeprime} / phi^{prime}$ remains finite at $eta=0$ and $eta=infty$.



Question:



As can be seen from the existence of an arbitrary function, the conditions so far do not suffice to determine $b$ and $c$. What further conditions can I reasonably impose to uniquely determine $b$ and $c$?
The term $b dot c$ can be regarded as a physical source term where $dot{()}$ indicates a derivative. Ideally, this source term would indicate no spurious sources. I have considered setting $nabla^2 b = 0$ which gives a 3rd order ODE but I see no formal justification for doing so.










share|cite|improve this question




























    up vote
    0
    down vote

    favorite












    I want to express a divergence-free vector $vec{w}(x,y)$ with non-zero curl in the domain $x,y geq 0$ as
    $$ vec{w} = b nabla c ; .$$
    The vector is normal to the boundary $partial Omega$. Furthermore, the value of $c$ on the boundary is zero ($c=0$ on $x,y=0$). Based on a helpful answer I received to this question in another post, I have made some progress. However, I am now unsure about the uniqueness of the problem.



    To illustrate the problem, I will assume the model
    $$ vec{w} = frac{1}{sqrt{x^2+y^2}} left( begin{array}{c} y \ -x end{array} right) ; .$$
    This vector has the properties that it is divergence-free, has non-zero curl, and is normal to the boundary as mentioned above. From the theory of first-order linear PDE's it follows that $c$ has the general solution
    $$ c = phi left( frac{x}{y} right) $$
    where $phi$ is an arbitrary function of the variable $eta equiv x/y$. From this follows the solution for $b$ as
    $$ b = frac{y^2}{sqrt{x^2+y^2}} frac{1}{phi^{prime}} $$
    where the prime indicates a derivative. Taking the gradient results in
    $$ nabla b = left( begin{array}{c}
    -frac{eta}{(1+eta^2)^{3/2}} frac{1}{phi^{prime}} - frac{phi^{primeprime}}{{phi^{prime}}^2} frac{1}{sqrt{1+eta^2}} \
    frac{2}{sqrt{1+eta^2}} frac{1}{phi^{prime}} - frac{1}{(1+eta^2)^{3/2}} frac{1}{phi^{prime}} + frac{phi^{primeprime}}{{phi^{prime}}^2} frac{eta}{sqrt{1+eta^2}}
    end{array} right) ; . $$

    These solutions satisfy the rotational and the divergence-free condition and the equation $vec{w} = b nabla c$. In order to impose the boundary conditions for $c$, one needs $phi(0) = 0$ and $phi(eta) = 0$ as $eta rightarrow infty$, as well as that $phi^{primeprime} / phi^{prime}$ remains finite at $eta=0$ and $eta=infty$.



    Question:



    As can be seen from the existence of an arbitrary function, the conditions so far do not suffice to determine $b$ and $c$. What further conditions can I reasonably impose to uniquely determine $b$ and $c$?
    The term $b dot c$ can be regarded as a physical source term where $dot{()}$ indicates a derivative. Ideally, this source term would indicate no spurious sources. I have considered setting $nabla^2 b = 0$ which gives a 3rd order ODE but I see no formal justification for doing so.










    share|cite|improve this question


























      up vote
      0
      down vote

      favorite









      up vote
      0
      down vote

      favorite











      I want to express a divergence-free vector $vec{w}(x,y)$ with non-zero curl in the domain $x,y geq 0$ as
      $$ vec{w} = b nabla c ; .$$
      The vector is normal to the boundary $partial Omega$. Furthermore, the value of $c$ on the boundary is zero ($c=0$ on $x,y=0$). Based on a helpful answer I received to this question in another post, I have made some progress. However, I am now unsure about the uniqueness of the problem.



      To illustrate the problem, I will assume the model
      $$ vec{w} = frac{1}{sqrt{x^2+y^2}} left( begin{array}{c} y \ -x end{array} right) ; .$$
      This vector has the properties that it is divergence-free, has non-zero curl, and is normal to the boundary as mentioned above. From the theory of first-order linear PDE's it follows that $c$ has the general solution
      $$ c = phi left( frac{x}{y} right) $$
      where $phi$ is an arbitrary function of the variable $eta equiv x/y$. From this follows the solution for $b$ as
      $$ b = frac{y^2}{sqrt{x^2+y^2}} frac{1}{phi^{prime}} $$
      where the prime indicates a derivative. Taking the gradient results in
      $$ nabla b = left( begin{array}{c}
      -frac{eta}{(1+eta^2)^{3/2}} frac{1}{phi^{prime}} - frac{phi^{primeprime}}{{phi^{prime}}^2} frac{1}{sqrt{1+eta^2}} \
      frac{2}{sqrt{1+eta^2}} frac{1}{phi^{prime}} - frac{1}{(1+eta^2)^{3/2}} frac{1}{phi^{prime}} + frac{phi^{primeprime}}{{phi^{prime}}^2} frac{eta}{sqrt{1+eta^2}}
      end{array} right) ; . $$

      These solutions satisfy the rotational and the divergence-free condition and the equation $vec{w} = b nabla c$. In order to impose the boundary conditions for $c$, one needs $phi(0) = 0$ and $phi(eta) = 0$ as $eta rightarrow infty$, as well as that $phi^{primeprime} / phi^{prime}$ remains finite at $eta=0$ and $eta=infty$.



      Question:



      As can be seen from the existence of an arbitrary function, the conditions so far do not suffice to determine $b$ and $c$. What further conditions can I reasonably impose to uniquely determine $b$ and $c$?
      The term $b dot c$ can be regarded as a physical source term where $dot{()}$ indicates a derivative. Ideally, this source term would indicate no spurious sources. I have considered setting $nabla^2 b = 0$ which gives a 3rd order ODE but I see no formal justification for doing so.










      share|cite|improve this question















      I want to express a divergence-free vector $vec{w}(x,y)$ with non-zero curl in the domain $x,y geq 0$ as
      $$ vec{w} = b nabla c ; .$$
      The vector is normal to the boundary $partial Omega$. Furthermore, the value of $c$ on the boundary is zero ($c=0$ on $x,y=0$). Based on a helpful answer I received to this question in another post, I have made some progress. However, I am now unsure about the uniqueness of the problem.



      To illustrate the problem, I will assume the model
      $$ vec{w} = frac{1}{sqrt{x^2+y^2}} left( begin{array}{c} y \ -x end{array} right) ; .$$
      This vector has the properties that it is divergence-free, has non-zero curl, and is normal to the boundary as mentioned above. From the theory of first-order linear PDE's it follows that $c$ has the general solution
      $$ c = phi left( frac{x}{y} right) $$
      where $phi$ is an arbitrary function of the variable $eta equiv x/y$. From this follows the solution for $b$ as
      $$ b = frac{y^2}{sqrt{x^2+y^2}} frac{1}{phi^{prime}} $$
      where the prime indicates a derivative. Taking the gradient results in
      $$ nabla b = left( begin{array}{c}
      -frac{eta}{(1+eta^2)^{3/2}} frac{1}{phi^{prime}} - frac{phi^{primeprime}}{{phi^{prime}}^2} frac{1}{sqrt{1+eta^2}} \
      frac{2}{sqrt{1+eta^2}} frac{1}{phi^{prime}} - frac{1}{(1+eta^2)^{3/2}} frac{1}{phi^{prime}} + frac{phi^{primeprime}}{{phi^{prime}}^2} frac{eta}{sqrt{1+eta^2}}
      end{array} right) ; . $$

      These solutions satisfy the rotational and the divergence-free condition and the equation $vec{w} = b nabla c$. In order to impose the boundary conditions for $c$, one needs $phi(0) = 0$ and $phi(eta) = 0$ as $eta rightarrow infty$, as well as that $phi^{primeprime} / phi^{prime}$ remains finite at $eta=0$ and $eta=infty$.



      Question:



      As can be seen from the existence of an arbitrary function, the conditions so far do not suffice to determine $b$ and $c$. What further conditions can I reasonably impose to uniquely determine $b$ and $c$?
      The term $b dot c$ can be regarded as a physical source term where $dot{()}$ indicates a derivative. Ideally, this source term would indicate no spurious sources. I have considered setting $nabla^2 b = 0$ which gives a 3rd order ODE but I see no formal justification for doing so.







      pde vectors vector-analysis boundary-value-problem






      share|cite|improve this question















      share|cite|improve this question













      share|cite|improve this question




      share|cite|improve this question








      edited 2 days ago

























      asked Nov 22 at 21:06









      Crenguta

      406




      406



























          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: 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%2f3009661%2funiqueness-of-a-vector-decomposition%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown






























          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%2f3009661%2funiqueness-of-a-vector-decomposition%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

          Sphinx de Gizeh

          Dijon

          Guerrita