Are you ready to take your data science career to the next step, or break into data science? With Springboard’s Data Science Career Track, you’ll master data science topics, have personalized career guidance, weekly calls with a data science expert, and a job guarantee.
Solved
3D refraction of ray / vector operations
Posted on 2014-02-27
Dear experts,
I've been working at a problem for a while. If anyone can help, I'd really appreciate it. I'm working on refraction of an ultrasound beam as it crosses a planar interface between a Rexolite wedge and a steel sample block. I've got it mostly worked out, but am struggling with a (the) key parameter: How to define the (3D) vector after refraction across the plane. I can get the angle with relation to the plane-normal by Snell's Law, and I'm assuming the vector will be in the plane defined by the normal and the incident wave, and the line will contain the point at which it intersects the plane. This seems like it should be easy given the information I've already calculated, but I just can't seem to get the keystone in the top of the arch. 20140227-BeamPathWork-CloseUp.pdf
Can anyone help?
Thank You,
-Matt
I've been working at a problem for a while. If anyone can help, I'd really appreciate it. I'm working on refraction of an ultrasound beam as it crosses a planar interface between a Rexolite wedge and a steel sample block. I've got it mostly worked out, but am struggling with a (the) key parameter: How to define the (3D) vector after refraction across the plane. I can get the angle with relation to the plane-normal by Snell's Law, and I'm assuming the vector will be in the plane defined by the normal and the incident wave, and the line will contain the point at which it intersects the plane. This seems like it should be easy given the information I've already calculated, but I just can't seem to get the keystone in the top of the arch. 20140227-BeamPathWork-CloseUp.pdf
Can anyone help?
Thank You,
-Matt
6-
6
12 Comments
Active today
If your unit incident ray is i and your unit normal is n,
and the ratio of wavelength between the two media is r,
then your refracted ray will be
ri+n(rn·i-sqrt(1-r^2+(rn·i)^2)))
and the ratio of wavelength between the two media is r,
then your refracted ray will be
ri+n(rn·i-sqrt(1-r^2+(rn·i)^2)))
Sorry, I have a follow-up: Does this equation have a name? I wanted to look more into it to see if there were any assumptions (small-angle, etc.). I'm getting angles that differ slightly from direct Snell results.
Active today
If the number in the square root is negative, you'd have total internal reflection.
What are the numbers you are using when you get slightly different angles?
What are the numbers you are using when you get slightly different angles?
If I try the simple case:
i=<0,0,1>
n=<1/Sqrt[2],0,1/Sqrt[2]>
r=0.73
I get <-0.28445, 0., 1.01465> which is 29.3 degrees from -n
If I use, Snell's Law:
Theta2 = 45 Degree;
Theta1 = ArcSin[Sin[Theta2]*0.73] = 31.1
I'm sure it's user error, but I'm searching for it.
Thanks again!
-Matt
i=<0,0,1>
n=<1/Sqrt[2],0,1/Sqrt[2]>
r=0.73
I get <-0.28445, 0., 1.01465> which is 29.3 degrees from -n
If I use, Snell's Law:
Theta2 = 45 Degree;
Theta1 = ArcSin[Sin[Theta2]*0.73] = 31.1
I'm sure it's user error, but I'm searching for it.
Thanks again!
-Matt
Hmm... That gives the right angle, but looks to be the mirror image across the refraction plane (i.e. reflected). BeamPathWork-3-copy.pdf
Do you have a source, or are you deriving?
Thanks for the help!
Do you have a source, or are you deriving?
Thanks for the help!
I'm sorry, can you spell it out for me? I tried changing all of the i's to (-i)'s, and a few other permutations, but I failed.
Active today
i=<0,0,1>
n=<0.707106781186548,0,0.7
r=0.73
ir-n(ir·n-sqrt(1-r^2+(ir·n
-ir-n(-ir·n-sqrt(1-r^2+(-i
Is this what you got, and does it make sense with your conventions?
n=<0.707106781186548,0,0.7
r=0.73
ir-n(ir·n-sqrt(1-r^2+(ir·n
-ir-n(-ir·n-sqrt(1-r^2+(-i
Is this what you got, and does it make sense with your conventions?
Active today
I'm not sure what convention to use, in http:#a39894769, the norm was pointing with the incident ray, but in http:#a39896563, the norm is pointing against the incident ray, so the sign convention changes.
With
i={0,0,1}
n={0.70710678118655,0,-0.7
r=0.73
the ray that makes sense to me when I draw it is
ri - n(ri·n + sqrt(1 - r^2 + (ri·n)^2)) = {-0.240619517519044,0,0.97
With
i={0,0,1}
n={0.70710678118655,0,-0.7
r=0.73
the ray that makes sense to me when I draw it is
ri - n(ri·n + sqrt(1 - r^2 + (ri·n)^2)) = {-0.240619517519044,0,0.97
Featured Post
There are many ways to learn to code these days. From coding bootcamps like Flatiron School to online courses to totally free beginner resources. The best way to learn to code depends on many factors, but the most important one is you. See what course is best for you.
Question has a verified solution.
Are you are experiencing a similar issue? Get a personalized answer when you ask a related question.
Have a better answer? Share it in a comment.
Join & Write a Comment Already a member? Login.
Complex Numbers are funny things. Many people have a basic understanding of them, some a more advanced. The confusion usually arises when that pesky i (or j for Electrical Engineers) appears and understanding the meaning of a square root of a nega…
Have you ever thought of installing a power system that generates solar electricity to power your house? Some may say yes, while others may tell me no. But have you noticed that people around you are now considering installing such systems in their …
This is a video describing the growing solar energy use in Utah. This is a topic that greatly interests me and so I decided to produce a video about it.
Although Jacob Bernoulli (1654-1705) has been credited as the creator of "Binomial Distribution Table", Gottfried Leibniz (1646-1716) did his dissertation on the subject in 1666; Leibniz you may recall is the co-inventor of "Calculus" and beat Isaac…
Suggested Courses
Course of the Month9 days, 3 hours left to enroll
595 members asked questions and received personalized solutions in the past 7 days.
Join the community of 500,000 technology professionals and ask your questions.