Finite difference for non-uniform unstructured mesh/stencilHow to obtain (prove) 5-stencil formula for 2nd derivative?Finite difference implicit schema for wave equation 1d not unconditionally stable?Calculating gradient from finite difference resultsDiscretization of the Anisotropic Diffusion Operator for Finite Difference MethodNumerically Solving a Poisson Equation with Neumann Boundary Conditionsfinite difference method with non-uniform meshFinite-Difference Approximation of Mixed Derivative in Spherical Coordinates?Stiffness Matrix for Galerkin Method (Finite Element Approx)Finite difference method for non-uniform gridApproximating finite differences by higher order derivatives of continuous functions

Why was Sir Cadogan fired?

Is there a hemisphere-neutral way of specifying a season?

How obscure is the use of 令 in 令和?

What exactly is ineptocracy?

Do creatures with a listed speed of "0 ft., fly 30 ft. (hover)" ever touch the ground?

How to find if SQL server backup is encrypted with TDE without restoring the backup

Ambiguity in the definition of entropy

How could indestructible materials be used in power generation?

How to travel to Japan while expressing milk?

What is the most common color to indicate the input-field is disabled?

How to prevent "they're falling in love" trope

How to coordinate airplane tickets?

Finitely generated matrix groups whose eigenvalues are all algebraic

What is the fastest integer factorization to break RSA?

GFCI outlets - can they be repaired? Are they really needed at the end of a circuit?

Why are UK visa biometrics appointments suspended at USCIS Application Support Centers?

How to enclose theorems and definition in rectangles?

Why do I get negative height?

How to install cross-compiler on Ubuntu 18.04?

How can saying a song's name be a copyright violation?

What is the opposite of "eschatology"?

Was the old ablative pronoun "med" or "mēd"?

How to stretch the corners of this image so that it looks like a perfect rectangle?

Rotate ASCII Art by 45 Degrees



Finite difference for non-uniform unstructured mesh/stencil


How to obtain (prove) 5-stencil formula for 2nd derivative?Finite difference implicit schema for wave equation 1d not unconditionally stable?Calculating gradient from finite difference resultsDiscretization of the Anisotropic Diffusion Operator for Finite Difference MethodNumerically Solving a Poisson Equation with Neumann Boundary Conditionsfinite difference method with non-uniform meshFinite-Difference Approximation of Mixed Derivative in Spherical Coordinates?Stiffness Matrix for Galerkin Method (Finite Element Approx)Finite difference method for non-uniform gridApproximating finite differences by higher order derivatives of continuous functions













0












$begingroup$


Below I have shown my non-uniform unstructured mesh (as in there is no pattern between the relative size of $h_i$ and $h_i+1$ etc.



NU-US-M



I've made the following equations using nodes $i-1$, $i$, and $i+1$.



  1. $phi_i-1=phi_i+h_i-1fracdphidx|_i+frach_i-1^22!fracd^2phidx^2|_i+frach_i-1^33!fracd^3phidx^3|_i+...$

  2. $phi_i+1=phi_i+h_ifracdphidx|_i+frach_i^22!fracd^2phidx^2|_i+frach_i^33!fracd^3phidx^3|_i+...$

Combining these in a linear combination, then solving for the first derivative
$$fracdphidx|_i=fracAphi_i-1+Bphi_i+1-(A+B)phi_iAh_i-1+Bh_i-frac12left(fracAh_i-1^2+Bh_i^2Ah_i-1+Bh_iright)fracd^2phidx^2|_i-frac16left(fracAh_i-1^3+Bh_i^3Ah_i-1+Bh_iright)fracd^3phidx^3|_i+/-...$$



According to my notes from class, I should be getting a final answer of
$$fracdphidx|_i=fracphi_i+1-phi_i-1x_i+1-x_i-1-fracleft(h_iright)^2-left(h_i-1right)^22left(x_i+1-x_i-1right)fracd^2phidx^2|_i-fracleft(h_iright)^3-left(h_i-1right)^36left(x_i+1-x_i-1right)fracd^3phidx^3|_i+O(h^3).$$



Based on what I am expecting to get, and the fact that I want a second order accurate final equation I can say the following.



  1. $Ah_i-1^2+Bh_i^2=0$

  2. $A+B=0$

However, I can't get the expected answer. Can I get some help with the error I'm making?










share|cite|improve this question









$endgroup$
















    0












    $begingroup$


    Below I have shown my non-uniform unstructured mesh (as in there is no pattern between the relative size of $h_i$ and $h_i+1$ etc.



    NU-US-M



    I've made the following equations using nodes $i-1$, $i$, and $i+1$.



    1. $phi_i-1=phi_i+h_i-1fracdphidx|_i+frach_i-1^22!fracd^2phidx^2|_i+frach_i-1^33!fracd^3phidx^3|_i+...$

    2. $phi_i+1=phi_i+h_ifracdphidx|_i+frach_i^22!fracd^2phidx^2|_i+frach_i^33!fracd^3phidx^3|_i+...$

    Combining these in a linear combination, then solving for the first derivative
    $$fracdphidx|_i=fracAphi_i-1+Bphi_i+1-(A+B)phi_iAh_i-1+Bh_i-frac12left(fracAh_i-1^2+Bh_i^2Ah_i-1+Bh_iright)fracd^2phidx^2|_i-frac16left(fracAh_i-1^3+Bh_i^3Ah_i-1+Bh_iright)fracd^3phidx^3|_i+/-...$$



    According to my notes from class, I should be getting a final answer of
    $$fracdphidx|_i=fracphi_i+1-phi_i-1x_i+1-x_i-1-fracleft(h_iright)^2-left(h_i-1right)^22left(x_i+1-x_i-1right)fracd^2phidx^2|_i-fracleft(h_iright)^3-left(h_i-1right)^36left(x_i+1-x_i-1right)fracd^3phidx^3|_i+O(h^3).$$



    Based on what I am expecting to get, and the fact that I want a second order accurate final equation I can say the following.



    1. $Ah_i-1^2+Bh_i^2=0$

    2. $A+B=0$

    However, I can't get the expected answer. Can I get some help with the error I'm making?










    share|cite|improve this question









    $endgroup$














      0












      0








      0





      $begingroup$


      Below I have shown my non-uniform unstructured mesh (as in there is no pattern between the relative size of $h_i$ and $h_i+1$ etc.



      NU-US-M



      I've made the following equations using nodes $i-1$, $i$, and $i+1$.



      1. $phi_i-1=phi_i+h_i-1fracdphidx|_i+frach_i-1^22!fracd^2phidx^2|_i+frach_i-1^33!fracd^3phidx^3|_i+...$

      2. $phi_i+1=phi_i+h_ifracdphidx|_i+frach_i^22!fracd^2phidx^2|_i+frach_i^33!fracd^3phidx^3|_i+...$

      Combining these in a linear combination, then solving for the first derivative
      $$fracdphidx|_i=fracAphi_i-1+Bphi_i+1-(A+B)phi_iAh_i-1+Bh_i-frac12left(fracAh_i-1^2+Bh_i^2Ah_i-1+Bh_iright)fracd^2phidx^2|_i-frac16left(fracAh_i-1^3+Bh_i^3Ah_i-1+Bh_iright)fracd^3phidx^3|_i+/-...$$



      According to my notes from class, I should be getting a final answer of
      $$fracdphidx|_i=fracphi_i+1-phi_i-1x_i+1-x_i-1-fracleft(h_iright)^2-left(h_i-1right)^22left(x_i+1-x_i-1right)fracd^2phidx^2|_i-fracleft(h_iright)^3-left(h_i-1right)^36left(x_i+1-x_i-1right)fracd^3phidx^3|_i+O(h^3).$$



      Based on what I am expecting to get, and the fact that I want a second order accurate final equation I can say the following.



      1. $Ah_i-1^2+Bh_i^2=0$

      2. $A+B=0$

      However, I can't get the expected answer. Can I get some help with the error I'm making?










      share|cite|improve this question









      $endgroup$




      Below I have shown my non-uniform unstructured mesh (as in there is no pattern between the relative size of $h_i$ and $h_i+1$ etc.



      NU-US-M



      I've made the following equations using nodes $i-1$, $i$, and $i+1$.



      1. $phi_i-1=phi_i+h_i-1fracdphidx|_i+frach_i-1^22!fracd^2phidx^2|_i+frach_i-1^33!fracd^3phidx^3|_i+...$

      2. $phi_i+1=phi_i+h_ifracdphidx|_i+frach_i^22!fracd^2phidx^2|_i+frach_i^33!fracd^3phidx^3|_i+...$

      Combining these in a linear combination, then solving for the first derivative
      $$fracdphidx|_i=fracAphi_i-1+Bphi_i+1-(A+B)phi_iAh_i-1+Bh_i-frac12left(fracAh_i-1^2+Bh_i^2Ah_i-1+Bh_iright)fracd^2phidx^2|_i-frac16left(fracAh_i-1^3+Bh_i^3Ah_i-1+Bh_iright)fracd^3phidx^3|_i+/-...$$



      According to my notes from class, I should be getting a final answer of
      $$fracdphidx|_i=fracphi_i+1-phi_i-1x_i+1-x_i-1-fracleft(h_iright)^2-left(h_i-1right)^22left(x_i+1-x_i-1right)fracd^2phidx^2|_i-fracleft(h_iright)^3-left(h_i-1right)^36left(x_i+1-x_i-1right)fracd^3phidx^3|_i+O(h^3).$$



      Based on what I am expecting to get, and the fact that I want a second order accurate final equation I can say the following.



      1. $Ah_i-1^2+Bh_i^2=0$

      2. $A+B=0$

      However, I can't get the expected answer. Can I get some help with the error I'm making?







      numerical-methods finite-differences finite-difference-methods






      share|cite|improve this question













      share|cite|improve this question











      share|cite|improve this question




      share|cite|improve this question










      asked Mar 20 at 18:37









      WnGatRC456WnGatRC456

      10811




      10811




















          0






          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',
          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%2f3155853%2ffinite-difference-for-non-uniform-unstructured-mesh-stencil%23new-answer', 'question_page');

          );

          Post as a guest















          Required, but never shown

























          0






          active

          oldest

          votes








          0






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes















          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%2f3155853%2ffinite-difference-for-non-uniform-unstructured-mesh-stencil%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

          Lowndes Grove History Architecture References Navigation menu32°48′6″N 79°57′58″W / 32.80167°N 79.96611°W / 32.80167; -79.9661132°48′6″N 79°57′58″W / 32.80167°N 79.96611°W / 32.80167; -79.9661178002500"National Register Information System"Historic houses of South Carolina"Lowndes Grove""+32° 48' 6.00", −79° 57' 58.00""Lowndes Grove, Charleston County (260 St. Margaret St., Charleston)""Lowndes Grove"The Charleston ExpositionIt Happened in South Carolina"Lowndes Grove (House), Saint Margaret Street & Sixth Avenue, Charleston, Charleston County, SC(Photographs)"Plantations of the Carolina Low Countrye

          random experiment with two different functions on unit interval Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 23, 2019 at 00:00UTC (8:00pm US/Eastern)Random variable and probability space notionsRandom Walk with EdgesFinding functions where the increase over a random interval is Poisson distributedNumber of days until dayCan an observed event in fact be of zero probability?Unit random processmodels of coins and uniform distributionHow to get the number of successes given $n$ trials , probability $P$ and a random variable $X$Absorbing Markov chain in a computer. Is “almost every” turned into always convergence in computer executions?Stopped random walk is not uniformly integrable

          How should I support this large drywall patch? Planned maintenance scheduled April 23, 2019 at 00:00UTC (8:00pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?How do I cover large gaps in drywall?How do I keep drywall around a patch from crumbling?Can I glue a second layer of drywall?How to patch long strip on drywall?Large drywall patch: how to avoid bulging seams?Drywall Mesh Patch vs. Bulge? To remove or not to remove?How to fix this drywall job?Prep drywall before backsplashWhat's the best way to fix this horrible drywall patch job?Drywall patching using 3M Patch Plus Primer