Solving Poisson-Boltzmann equation in Cylindrical and Spherical Coordinates

Click For Summary
SUMMARY

The discussion focuses on solving the Poisson-Boltzmann equation, represented as ∇²Φ = k²Φ, using both cylindrical and spherical coordinates. The method involves assuming a separable solution of the form Φ(r) = A(ρ)B(φ)C(z) for cylindrical coordinates or Φ(r) = R(r)Θ(θ)G(φ) for spherical coordinates. Participants emphasize the importance of boundary conditions to determine the symmetry of the solution and the necessity of eliminating variables to achieve a constant on one side of the equation. The Laplacian operator's definition in the respective coordinate systems is crucial for the solution process.

PREREQUISITES
  • Understanding of the Laplacian operator in cylindrical and spherical coordinates
  • Familiarity with separation of variables technique in partial differential equations
  • Knowledge of boundary conditions and their impact on solutions
  • Basic concepts of the Poisson-Boltzmann equation
NEXT STEPS
  • Study the derivation of the Laplacian operator in cylindrical coordinates
  • Learn about boundary value problems in partial differential equations
  • Explore the method of separation of variables in spherical coordinates
  • Investigate the physical applications of the Poisson-Boltzmann equation in electrostatics
USEFUL FOR

Students and researchers in applied mathematics, physics, and engineering, particularly those working on electrostatics and fluid dynamics involving the Poisson-Boltzmann equation.

thepopasmurf
Messages
73
Reaction score
0

Homework Statement



I don't have a specific problem in mind, it's more that I forgot how to solve the particular equation from first principles.

\nabla^{2} \Phi = k^{2}\Phi

Places I've looked so far have just quoted the results but I would like the complete method or the appropriate substitution.


Homework Equations



The relevant equations would be the definition of the Laplacian operator in cylindrical or spherical coordinates.



The Attempt at a Solution



My attempt was to sub in

\Phi = r^a \times exp(br)

and solve for a and b, but I got both a and b equal to plus and minus 1 without a way to eliminate the postitive a solution.
 
Physics news on Phys.org
Well, the usual 'first principles' way would be to assume a solution of the form ##\Phi(\mathbf{r}) = A(\rho)B(\phi)C(z)## in cylindrical coordinates or ##\Phi(\mathbf{r}) = R(r)\Theta(\theta)G(\phi)## in spherical coordinates. If you then divide the entire equation by ##\Phi## again, you can rearrange terms until one side of the equation depends only on one of the three variables, implying it must be a constant, which the other side of the equation is then also equal to.

For a very simple example, say you had the 2d problem ##\nabla^2 \Phi(x,y) = 0##, and you let ##\Phi(x,y) = X(x)Y(y)##. Plugging this in and dividing by ##\Phi## again gives ##\partial_x X(x)/X(x) + \partial_y Y(y)/Y(y) = 0##, or ##\partial_x X(x)/X(x) = -\partial_y Y(y)/Y(y)##, and hence it must be the case that ##\partial_x X(x)/X(x) = k## and ##-\partial_y Y(y)/Y(y) = k##, where k is some constant (not your equation's k, in this example!).

Also, keep in mind you will need some boundary conditions for the problem. Based on your boundary conditions, you may be able to argue that the solution should be spherically or cylindrically symmetric, so you can through out the angular or height dependences, but if you want an absolutely generally solution, you have to keep those terms as you could have boundaries conditions which vary with angle or height.
 
Last edited:

Similar threads

  • · Replies 7 ·
Replies
7
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
3
Views
2K
  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 17 ·
Replies
17
Views
4K
  • · Replies 9 ·
Replies
9
Views
2K