Skip to main content
Title capitalization
Link
DanielSank
  • 25k
  • 7
  • 81
  • 115

Validity of Finite Difference Methodfinite difference method

added tags, minor latex fixes
Source Link
Kyle Kanos
  • 28.8k
  • 41
  • 69
  • 135

I am using the Finite Difference Method to solve Poisson's equation

$\frac{\partial \phi}{\partial z^2} = \frac{\rho}{\epsilon}$$$\frac{\partial \phi}{\partial z^2} = \frac{\rho}{\epsilon}$$

To do it is discretized according to the Finite Difference Approximation of the second order derivative yielding the following set of equations for each grid point

$\frac{1}{2\Delta z^2}(\phi_{i+1}+\phi_{i-1}-2\phi_{i}) = \frac{\rho_{i}}{\epsilon} $$$\frac{1}{2\Delta z^2}(\phi_{i+1}+\phi_{i-1}-2\phi_{i}) = \frac{\rho_{i}}{\epsilon} $$

My question is the following: Is there any simple way to see how many points and thus what spacing $\Delta z$ one should choose to produce reliable results? Obviously you should choose more than 5 but should I choose 500, 5000 or 50000? Asked another way: Over what length scale will the electrostatic potential vary significantly? I suppose if this is known one should make the spacing smaller than this length scale.

I once did a similar thing for the Schrödinger Equation. In this case it was straightforward to see what the characteristic wavelength of the solutions were and then to choose the spacing much smaller than this.

I am using the Finite Difference Method to solve Poisson's equation

$\frac{\partial \phi}{\partial z^2} = \frac{\rho}{\epsilon}$

To do it is discretized according to the Finite Difference Approximation of the second order derivative yielding the following set of equations for each grid point

$\frac{1}{2\Delta z^2}(\phi_{i+1}+\phi_{i-1}-2\phi_{i}) = \frac{\rho_{i}}{\epsilon} $

My question is the following: Is there any simple way to see how many points and thus what spacing $\Delta z$ one should choose to produce reliable results? Obviously you should choose more than 5 but should I choose 500, 5000 or 50000? Asked another way: Over what length scale will the electrostatic potential vary significantly? I suppose if this is known one should make the spacing smaller than this length scale.

I once did a similar thing for the Schrödinger Equation. In this case it was straightforward to see what the characteristic wavelength of the solutions were and then to choose the spacing much smaller than this.

I am using the Finite Difference Method to solve Poisson's equation

$$\frac{\partial \phi}{\partial z^2} = \frac{\rho}{\epsilon}$$

To do it is discretized according to the Finite Difference Approximation of the second order derivative yielding the following set of equations for each grid point

$$\frac{1}{2\Delta z^2}(\phi_{i+1}+\phi_{i-1}-2\phi_{i}) = \frac{\rho_{i}}{\epsilon} $$

My question is the following: Is there any simple way to see how many points and thus what spacing $\Delta z$ one should choose to produce reliable results? Obviously you should choose more than 5 but should I choose 500, 5000 or 50000? Asked another way: Over what length scale will the electrostatic potential vary significantly? I suppose if this is known one should make the spacing smaller than this length scale.

I once did a similar thing for the Schrödinger Equation. In this case it was straightforward to see what the characteristic wavelength of the solutions were and then to choose the spacing much smaller than this.

Source Link
user13514
  • 681
  • 6
  • 10

Validity of Finite Difference Method

I am using the Finite Difference Method to solve Poisson's equation

$\frac{\partial \phi}{\partial z^2} = \frac{\rho}{\epsilon}$

To do it is discretized according to the Finite Difference Approximation of the second order derivative yielding the following set of equations for each grid point

$\frac{1}{2\Delta z^2}(\phi_{i+1}+\phi_{i-1}-2\phi_{i}) = \frac{\rho_{i}}{\epsilon} $

My question is the following: Is there any simple way to see how many points and thus what spacing $\Delta z$ one should choose to produce reliable results? Obviously you should choose more than 5 but should I choose 500, 5000 or 50000? Asked another way: Over what length scale will the electrostatic potential vary significantly? I suppose if this is known one should make the spacing smaller than this length scale.

I once did a similar thing for the Schrödinger Equation. In this case it was straightforward to see what the characteristic wavelength of the solutions were and then to choose the spacing much smaller than this.