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 ''splain this hydrostatic paradox in tiny words'
This is (ostensibly) not a trick shot or video*. The scale was balanced before any blue water was added. 550mL of blue water was added to the left side. only 60mL of water needed to be added to the right side to re-balance the scale. Apparently, the scale will balance when the height of the two columns is equal. The left side of the scale only feels the weight of the column above the lower "tail" of the funnel (i.e. 60mL). So where does the weight of the remaining (550-60=) 490mL go...
Consider an extremely long and perfectly calibrated scale. A car with a mass of 1000 kg is placed on it, and the scale registers this weight accurately. Now, suppose the car begins to move, reaching very high speeds. Neglecting air resistance and rolling friction, if the car attains, for example, a velocity of 500 km/h, will the scale still indicate a weight corresponding to 1000 kg, or will the measured value decrease as a result of the motion? In a second scenario, imagine a person with a...
Scalar and vector potentials in Coulomb gauge Assume Coulomb gauge so that $$\nabla \cdot \mathbf{A}=0.\tag{1}$$ The scalar potential ##\phi## is described by Poisson's equation $$\nabla^2 \phi = -\frac{\rho}{\varepsilon_0}\tag{2}$$ which has the instantaneous general solution given by $$\phi(\mathbf{r},t)=\frac{1}{4\pi\varepsilon_0}\int \frac{\rho(\mathbf{r}',t)}{|\mathbf{r}-\mathbf{r}'|}d^3r'.\tag{3}$$ In Coulomb gauge the vector potential ##\mathbf{A}## is given by...
Back
Top