Finite Difference Discretization of a Fourth Order Partial Differential Term

Click For Summary

Discussion Overview

The discussion centers around the finite difference discretization of a specific fourth order partial differential term, \(\frac{\partial^4\phi}{\partial x^2\partial y^2}\). Participants explore various formulations and approaches to discretize this term, considering both theoretical and practical implications.

Discussion Character

  • Technical explanation
  • Mathematical reasoning
  • Debate/contested

Main Points Raised

  • One participant asks for a finite-difference discretization of the term \(\frac{\partial^4\phi}{\partial x^2\partial y^2}\).
  • Another participant provides a common discretization approach, emphasizing that there are many possible choices depending on the situation.
  • A different participant presents a similar discretization but expresses a desire for a formulation that avoids bivariate terms, questioning how to construct a solvable matrix with mixed terms.
  • A subsequent post clarifies a previous misunderstanding regarding notation, providing a detailed expansion of the discretization.
  • One participant expresses gratitude for the assistance received in the discussion.

Areas of Agreement / Disagreement

Participants do not reach a consensus on a specific discretization method, as multiple approaches are presented and some participants express uncertainty about the use of bivariate terms.

Contextual Notes

There are unresolved questions regarding the formulation of discretizations that avoid bivariate terms and the implications for constructing solvable matrices.

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
2K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 12 ·
Replies
12
Views
7K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 5 ·
Replies
5
Views
4K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K