Finite Difference Discretization of a Fourth Order Partial Differential Term

Hypatio
Messages
147
Reaction score
1
What is a finite-difference discretization for the partial differential term:

[tex]\frac{\partial^4\phi}{\partial x^2\partial y^2}[/tex]

Thanks in advance.
 
Last edited:
Physics news on Phys.org
There are as usual an infinite number, the choice of which depends upon the situation. A common choice would be
u(x-1,y+1) -2u(x+0,y+1) +u(x+1,y+1)
-2u(x-1,y+0) +4u(x+0,y+0) -2u(x+1,y+0)
+u(x-1,y-1) -2u(x+0,y-1) +u(x+1,y-1)
 
Last edited:
lurflurf said:
There are as usual an infinite number, the choice of which depends upon the situation. A common choice would be
u(x-1,y+1) -2u(x+0,y+1) u(x+1,y+1)
-2u(x-1,y+0) +4u(x+0,y+0) -2u(x+1,y+0)
u(x-1,y-1) -2u(x+0,y-1) u(x+1,y-1)

Thanks, but is there a formulation which does not use bivariate (?) terms? eg. 2u(x+0,y+1)*u(x+1,y+1)

I do not understand how it is possible to create a solvable matrix with a mix of bivariate terms in it, which is what I am trying to to do (stress function solution using gaussian elimination).
 
Last edited:
^Sorry that was implied addition not implied multiplication. Writen out in full we have
16uxxyy~u(x-1,y+1)-2u(x+0,y+1)+u(x+1,y+1)-2u(x-1,y+0)+4u(x+0,y+0)-2u(x+1,y+0)+u(x-1,y-1) -2u(x+0,y-1) +u(x+1,y-1)

or unscaled

(4st)2 uxxyy~u(x-s,y+t)-2u(x+0,y+t)+u(x+s,y+t)-2u(x-s,y+0)+4u(x+0,y+t)-2u(x+s,y+0)+u(x-s,y-t) -2u(x+0,y-t) +u(x+s,y-t)

A higher order or biased expansion could be used if needed.
 
Ah, thank you for alleviating my fears, and for the help.
 

Similar threads

  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 12 ·
Replies
12
Views
8K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 5 ·
Replies
5
Views
5K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K