Plane Wave Expansion: Solving Issues with Bohren and Huffman

AI Thread Summary
The discussion centers on challenges faced while implementing a plane wave expansion in spherical harmonics using MATLAB, based on the work of Bohren and Huffman. The user has encountered issues when defining certain components as zero, leading to nonsensical results due to terms becoming infinite as they are divided by sin(m*phi), which approaches zero. Despite the code being largely complete, the user suspects a deeper issue related to the equations referenced in Bohren and Huffman. They are seeking assistance or examples of similar expansions within a MATLAB framework to resolve these problems. Any guidance from others familiar with this field would be greatly appreciated.
falcao
Messages
9
Reaction score
0
I've been working on a plane wave expansion in spherical harmonics.. as demonstrated in Bohren and Huffman (if anyone's familiar with it_)... I have more or less finished my code in Matlab and I'm putting in values to expand and compare with the exact solution...


in short, I'm having trouble when defining some of my components as zero... getting results going crazy all over the place.! now I appreciate this would seem like a problem with the code but looking at the equations in Bohren and Huffman(4.19) (and in my code) there seems to be a different problem!..

one of the terms is divisible by a sin(m*phi) value and of course this value is going to go to zero with my phi and then the term goes to infinity! so my result ends up being nonsensical..>

if anyone is familiar with the field t'would be nice to hear!:redface:
 
Physics news on Phys.org
or even if anyone has any examples of such an expansion within a MATLAB framework>?!
 
honestly... any help would be massively appreciated!

please!:cry::confused::cry:
 
Thread 'Gauss' law seems to imply instantaneous electric field'
Imagine a charged sphere at the origin connected through an open switch to a vertical grounded wire. We wish to find an expression for the horizontal component of the electric field at a distance ##\mathbf{r}## from the sphere as it discharges. By using the Lorenz gauge condition: $$\nabla \cdot \mathbf{A} + \frac{1}{c^2}\frac{\partial \phi}{\partial t}=0\tag{1}$$ we find the following retarded solutions to the Maxwell equations If we assume that...
Maxwell’s equations imply the following wave equation for the electric field $$\nabla^2\mathbf{E}-\frac{1}{c^2}\frac{\partial^2\mathbf{E}}{\partial t^2} = \frac{1}{\varepsilon_0}\nabla\rho+\mu_0\frac{\partial\mathbf J}{\partial t}.\tag{1}$$ I wonder if eqn.##(1)## can be split into the following transverse part $$\nabla^2\mathbf{E}_T-\frac{1}{c^2}\frac{\partial^2\mathbf{E}_T}{\partial t^2} = \mu_0\frac{\partial\mathbf{J}_T}{\partial t}\tag{2}$$ and longitudinal part...
Thread 'Recovering Hamilton's Equations from Poisson brackets'
The issue : Let me start by copying and pasting the relevant passage from the text, thanks to modern day methods of computing. The trouble is, in equation (4.79), it completely ignores the partial derivative of ##q_i## with respect to time, i.e. it puts ##\partial q_i/\partial t=0##. But ##q_i## is a dynamical variable of ##t##, or ##q_i(t)##. In the derivation of Hamilton's equations from the Hamiltonian, viz. ##H = p_i \dot q_i-L##, nowhere did we assume that ##\partial q_i/\partial...
Back
Top