Sunday, 21 September 2014

acoustics - Nonsensical dispersion relations for elastic wave propagation


In an earlier question about Einstein notation, a link was provided to a medical paper which used acoustic propagation to noninvasively detect the orientation of muscle fibers. In short, muscle fibers are roughly a transverse isotropic material, with Hooke's law given by σ=(C11C112C66C13000C112C66C11C13000C13C13C33000000C44000000C44000000C66)(ϵ11ϵ22ϵ332ϵ232ϵ312ϵ12)

where ϵjk=12(uj,k+uk,j) where u is the material displacement from equilibrium. Solving the equation of motion σ=ρ¨u using plane wave solutions of the form u=Aei(ωtkr) yields the dispersion relation ΓA=0
where Γ=ρω2I(C11k21+C66k22+C44k23(C66C11)k1k2(C13+C44)k1k3(C66C11)k1k2C66k21+C11k22+C44k23(C13+C44)k2k3(C13+C44)k1k3(C13+C44)k2k3C44k21+C44k22+C33k23)
where I is the identity matrix. Since the k2 and k3 directions are identical, WLOG you can assume k2=0 to give Γ=(ρω2C11k21C44k230(C13+C44)k1k30ρω2C66k21C44k230(C13+C44)k1k30ρω2C44k21C33k23).
By choosing a particular material and frequency ω, plotting the regions where the dispersion relation holds as a function of wavevector k=(k1,0,k3), one gets beautiful images, such as the following picture of the allowed regions using data values for graphite-epoxy fiber lamina (see Table 3) using {C33 -> 23.4, C11 -> 2.1, C13 -> 0.943, C66 -> 0.527, C44 -> 1.03, \[Rho] -> 1, \[Omega] -> 4} and 6.8k1,k36.8:


enter image description here


The horizontal axis is the k3-axis, and the vertical axis is the k1-axis. The red mode corresponds to a pure shear-wave, and the green and blue modes correspond to quasi-transverse and quasi-longitudinal modes.


If you just make up some random unphysical numbers for the Cij, you occasionally get something which looks like this, using {C11 -> 0.192954, C33 -> 0.318008, C44 -> 0.269456, C13 -> 0.782727, C66 -> 0.195521, \[Rho] -> 1, \[Omega] -> 1} and 6.8k1,k36.8:


enter image description here



In the above, there is one propagation mode (in blue) which extends out to infinite radius in k-space, which seems to physically makes no sense, since ω is fixed, so you'd expect the dispersion relation to hold only within a bounded annulus in k-space (infinite spatial frequency makes no sense).


So this is probably a simple question, but is there any type of material which can have this sort of dispersion relation, or is there some reason why real materials can never behave in this manner (ie, they need to have bounded spatial wavenumber for a given acoustic frequency)?


For reference, here is the code I used to derive the above equations and make the images. The images above are essentially defined by cj(k1,k3)=ea|λj(k1,0,k3)|

where cj is the jth RGB color channel, λj(k1,0,k3) is the jth eigenvalue of Γ evaulated at (k1,0,k3), and a is a constant which affects the visibility of the curves.




No comments:

Post a Comment

Understanding Stagnation point in pitot fluid

What is stagnation point in fluid mechanics. At the open end of the pitot tube the velocity of the fluid becomes zero.But that should result...