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
$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.
I've made the following equations using nodes $i-1$, $i$, and $i+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+...$
- $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.
- $Ah_i-1^2+Bh_i^2=0$
- $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
$endgroup$
add a comment |
$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.
I've made the following equations using nodes $i-1$, $i$, and $i+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+...$
- $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.
- $Ah_i-1^2+Bh_i^2=0$
- $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
$endgroup$
add a comment |
$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.
I've made the following equations using nodes $i-1$, $i$, and $i+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+...$
- $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.
- $Ah_i-1^2+Bh_i^2=0$
- $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
$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.
I've made the following equations using nodes $i-1$, $i$, and $i+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+...$
- $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.
- $Ah_i-1^2+Bh_i^2=0$
- $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
numerical-methods finite-differences finite-difference-methods
asked Mar 20 at 18:37
WnGatRC456WnGatRC456
10811
10811
add a comment |
add a comment |
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
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
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.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
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