Mehrstellenverfahren for different grid spacing along the three space directions

  • Thread starter Thread starter kolmog
  • Start date Start date
  • Tags Tags
    Grid Space
kolmog
Messages
10
Reaction score
0
Hi everybody,

I need to discretize a differential equation. The grid that I am considering has a different grid spacing along the three space directions. That is, hx different than hy and hz. For that purpose, I would like to discretize the differential equation through the Mehrstellenverfahren scheme. The problem is that in all the references that I have found, the equation with the corresponding weights assumes that the grid spacing is the same along each of the three space directions (hx=hx=hx=h). I show you the corresponding equation in the attached file.

I would like to use this scheme in my particular situation. Does someone know how this equation would look like?

Thank you very much in advance!

Miquel
 

Attachments

  • eq1.jpg
    eq1.jpg
    17.6 KB · Views: 404
Physics news on Phys.org
Thanks for the post! Sorry you aren't generating responses at the moment. Do you have any further information, come to any new conclusions or is it possible to reword the post?
 
Hi,

As I see that maybe not many people use this scheme to discretize differential equations, I would reword this post. I have discretized a differential equation (generalized Poisson equation) in a real space grid. If you take a look to the attached file, this is the first equation. The left hand side of this equation can be rewritten like in the second equation. Then, I have used Fornberg weights to discretize first and second derivatives. Do you know if, a part from Bernt Fornbverg weights, there exist a better option?
 

Attachments

  • discr.jpg
    discr.jpg
    11.7 KB · Views: 455
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...
Back
Top