Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Formal description of the PhysicalOpticsIntegral function #1

Open
terekib opened this issue Apr 19, 2023 · 5 comments
Open

Formal description of the PhysicalOpticsIntegral function #1

terekib opened this issue Apr 19, 2023 · 5 comments

Comments

@terekib
Copy link

terekib commented Apr 19, 2023

Hi!
I'm trying to develop a radar simulation (FMCW radar), so I am reading a lot of papers. Still I can't figure out the PhysicalOpticsIntegral function. I usually find formulas similar to your code, but always different in may ways. Could you help me understand why you are doing it that way or could you give me some links to papers on which your code is based on? Thank you very much :)

The ones I looked at:
https://ntrs.nasa.gov/api/citations/19950004523/downloads/19950004523.pdf
https://www.jpier.org/issues/volume.html?paper=10061807
https://dokumen.tips/documents/ray-tube-integration-in-shooting-and-bouncing-ray-method.html?page=3

@pingpongballz
Copy link

pingpongballz commented Dec 13, 2023

ok i'm not too sure, but here goes:

from [1] and [3], We can see that $\begin{bmatrix} E_\theta \\\ E_\phi \end{bmatrix} = \sum_{rays} \begin{bmatrix} B_\theta \\\ B_\phi \end{bmatrix} \frac{jk}{2\pi} (\delta A) S(\theta , \phi) e^{jkr}$ (eq 5) . Here, $ S(\theta , \phi)$ is assumed to be unity, if the ray tube area is sufficiently small. For clarity, also define $E_s = \begin{bmatrix} E_\theta \\\ E_\phi \end{bmatrix} \frac{e^{-jkr}}{r}$ . Also note that $j = -i$

$B_\theta =0.5( E_{ap} \times -\phi +Z_0 H_{ap} \times \theta ) \cdot \hat{k'}$, k' is the exit ray unit vector. (eq 6)[1]
$B_\phi =0.5( E_{ap} \times \theta +Z_0 H_{ap} \times \phi ) \cdot \hat{k'}$, k' is the exit ray unit vector. (eq 7)[1]

So this indeed leads to lines 225 to 234 in SbrSolver.hpp. Only difference is that the $Z_0$ term is missing. This can be explained as the author of the code defined $H = \hat{k} \times E$ and not $H = \frac{1}{Z_0} \hat{k} \times E$, so no need to multiply the $Z_0$ term, as $Z_0(\frac{1}{Z_0} \hat{k} \times E) = \hat{k} \times E$

However, in [2], $\hat{k'} $ is defined as the normal of the patch where the bounce happens, which is contradictory to [1] and [3]. this kinda sus if u ask me

[1] https://ntrs.nasa.gov/api/citations/19950004523/downloads/19950004523.pdf
[2] https://www.jpier.org/issues/volume.html?paper=10061807
[3] https://dokumen.tips/documents/ray-tube-integration-in-shooting-and-bouncing-ray-method.html?page=3

@terekib
Copy link
Author

terekib commented Dec 15, 2023

Thank you very much!

@xiongyuup
Copy link

Where is the material information of the media set in the RaytrAMP project

@pingpongballz
Copy link

Hi,

From what i remember, the only material supported is PEC. If you want to use other materials, the nature the code is written prohibits this.

A different formulation is required, one that accounts for the displacement current density that is induced within the dielectric, and calls for different boundary conditions other than PEC.

@xiongyuup
Copy link

Oh my God, thank you very much. Thank you very much. I thought this project was the same as those projects of particle light scattering. Recently, I am looking for a method that can quickly calculate the electromagnetic scattering characteristics of naval ships on the sea surface. I wonder if you have come across such codes. In addition, may I ask for your contact information?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants