Difference between revisions of "User:Tohline/Appendix/CGH/ParallelApertures2D"

From VistrailsWiki
Jump to navigation Jump to search
(Created page with '<!-- __FORCETOC__ will force the creation of a Table of Contents --> <!-- __NOTOC__ will force TOC off --> =CGH: 2D Rectangular Appertures that are Parallel to the Image Screen=…')
 
 
(11 intermediate revisions by the same user not shown)
Line 1: Line 1:
<!-- __FORCETOC__ will force the creation of a Table of Contents -->
__FORCETOC__<!-- will force the creation of a Table of Contents -->
<!-- __NOTOC__ will force TOC off -->
<!-- __NOTOC__ will force TOC off -->
=CGH:  2D Rectangular Appertures that are Parallel to the Image Screen=
=CGH:  2D Rectangular Apertures that are Parallel to the Image Screen=
This chapter is intended primarily to replicate [http://www.phys.lsu.edu/faculty/tohline/phys4412/howto/slit1d.html &sect;I.B from the online class notes] &#8212; see also an updated [[User:Tohline/Appendix/Ramblings#Computer-Generated_Holography|Table of Contents]] &#8212; that I developed in conjunction with a course that I taught in 1999 on the topic of ''Computer Generated Holography (CGH)'' for a subset of LSU physics majors who were interested in computational science.
This chapter is intended primarily to replicate [http://www.phys.lsu.edu/faculty/tohline/phys4412/howto/slit2d.html &sect;I.B from the online class notes] &#8212; see also an updated [[User:Tohline/Appendix/Ramblings#Computer-Generated_Holography|Table of Contents]] &#8212; that I developed in conjunction with a course that I taught in 1999 on the topic of ''Computer Generated Holography (CGH)'' for a subset of LSU physics majors who were interested in computational science.  This discussion parallels the somewhat more detailed one presented in [[User:Tohline/Appendix/CGH/ParallelApertures|&sect;I.A]] on the one-dimensional aperture oriented parallel to the image screen.




{{LSU_HBook_header}}
{{LSU_HBook_header}}


==Propagation of Light==
Here we reference heavily the traveling, wave-like nature of light and, as is customary, use <math>~c</math> to represent its speed of propagation through space.  Consider a single ray of light of wavelength,


==Utility of FFT Techniques==
Consider the amplitude (and phase) of light that is incident at a location (x<sub>1</sub>, y<sub>1</sub>) on an image screen that is located a distance Z from a rectangular aperture of width ''w'' and height ''h''.  By analogy with [[User:Tohline/Appendix/CGH/ParallelApertures#Utility_of_FFT_Techniques|our accompanying discussion in the context of 1D apertures]], the complex number, A, representing the light amplitude and phase at (x<sub>1</sub>, y<sub>1</sub>) will be,
<table border="0" cellpadding="5" align="center">


<tr>
  <td align="right">
<math>~A(x_1, y_1)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \sum_j \sum_k
a_{jk} e^{i(2\pi D_{jk} /\lambda + \phi_{jk})}
\, ,
</math>
  </td>
</tr>
</table>
where, here, the summations are taken over all "j,k" elements of light across the entire 2D aperture, and now the distance D<sub>jk</sub> is given by the expression,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~D^2_{jk}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
(X_j - x_1)^2 + (Y_k - y_1)^2 + Z^2 
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
Z^2 + y_1^2 - 2y_1 Y_k + Y_k^2 + x_1^2 - 2x_1 X_j + X_j^2
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
L^2 \biggl[1 - \frac{2(x_1 X_j + y_1 Y_k ) }{L^2}  + \frac{X_j^2 + Y_k^2}{L^2} \biggr] \, ,
</math>
  </td>
</tr>
</table>
and,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~L</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
[Z^2 + y_1^2  + x_1^2]^{1 / 2} \, .
</math>
  </td>
</tr>
</table>
If <math>~|X_j/L| \ll 1</math> and <math>~|Y_k/L| \ll 1</math> we can drop the quadratic terms in favor of the linear ones in the expression for  D<sub>jk</sub> and deduce that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~D_{jk}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~
L \biggl[1 - \frac{2(x_1 X_j + y_1 Y_k ) }{L^2}  \biggr]^{1 / 2}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~
L \biggl[1 - \frac{(x_1 X_j + y_1 Y_k ) }{L^2}  \biggr] \, .
</math>
  </td>
</tr>
</table>
Hence, the double-summation expression for the amplitude at screen location (x<sub>1</sub>, y<sub>1</sub>) becomes,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~A(x_1, y_1)</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ A_0 \sum_j \sum_k
a_{jk} e^{i\phi_{jk} } \cdot \exp\biggl\{ -i \biggl[ \frac{2\pi(x_1 X_j + y_1Y_k}{\lambda L} \biggr] \biggr\}
\, ,
</math>
  </td>
</tr>
</table>
where,
<div align="center">
<math>~A_0 \equiv e^{i(2\pi L/\lambda)}</math>.
</div>
When written in this form, it should be apparent why discrete Fourier transform techniques &#8212; specifically, 2D-FFT techniques &#8212; are useful tools for evaluation of the complex amplitude, A(x<sub>1</sub>, y<sub>1</sub>).
<font color="red"><b>NOTE:</b></font>  If x<sub>1</sub> and/or y<sub>1</sub> are ever comparable in size to Z &#8212; which may be the case for large apertures or for apertures tilted by nearly 90&deg; to the image screen &#8212; then the variation of L with image screen position cannot be ignored and, accordingly, the coefficient A<sub>0</sub> cannot be moved outside of the double summation.
==Analytic Results==
Now, as in [[User:Tohline/Appendix/CGH/ParallelApertures#Parallels_With_Example_.232|&sect;I.A]], we would like to consider the case where the rectangular aperture is divided into an infinite number of divisions in both the X and Y dimensions and convert the summations in this last expression into integrals whose limits in both directions are given by the edges of the aperture.  If we specifically consider the case where the aperture is assumed to be uniformly bright (''i.e.,'' a<sub>jk</sub>  = a<sub>0</sub> dX DY, and a<sub>0</sub> is the brightness per unit area), and the phase <math>~\phi_{jk} = 0</math> at all locations on the aperture, we can write,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~A(x_1, y_1)</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ A_0 a_0 \iint \exp\biggl\{ -i \biggl[ \frac{2\pi(x_1 X_j + y_1Y_k}{\lambda L} \biggr] \biggr\} dX dY
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ A_0 a_0
\int \exp\biggl\{ -i \biggl[ \frac{2\pi x_1 X_j }{\lambda L} \biggr] \biggr\} dX
\cdot \int \exp\biggl\{ -i \biggl[ \frac{2\pi y_1Y_k}{\lambda L} \biggr] \biggr\} dY
\, .
</math>
  </td>
</tr>
</table>
Both of these integrals can be completed in the same fashion as described in [[User:Tohline/Appendix/CGH/ParallelApertures#Parallels_With_Example_.232|&sect;I.A]] for the 1D slit, giving,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~A(x_1, y_1)</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ A_0 a_0 (wh) e^{-i\Theta_1} ~ e^{-i\vartheta_1} \mathrm{sinc}(\alpha_1) \mathrm{sinc}(\beta_1)
\, ,
</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\alpha_1</math>
  </td>
  <td align="center">
<math>~\equiv </math>
  </td>
  <td align="left">
<math>~\frac{\pi x_1 (X_1 - X_2)}{\lambda L}= \frac{\pi x_1 h}{\lambda L} \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Theta_1</math>
  </td>
  <td align="center">
<math>~\equiv </math>
  </td>
  <td align="left">
<math>~\frac{\pi x_1 (X_1 + X_2)}{\lambda L} \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\beta_1</math>
  </td>
  <td align="center">
<math>~\equiv </math>
  </td>
  <td align="left">
<math>~\frac{\pi y_1 (Y_1 - Y_2)}{\lambda L}= \frac{\pi y_1 w}{\lambda L} \, ,</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\vartheta_1</math>
  </td>
  <td align="center">
<math>~\equiv </math>
  </td>
  <td align="left">
<math>~\frac{\pi y_1 (Y_1 + Y_2)}{\lambda L} \, .</math>
  </td>
</tr>
</table>
It is worth noting that this derivation closely parallels the one presented in &sect;8.5.1 (p. 393) of Born &amp; Wolf (1980) (see the [[#See_Also|reference below]]).  Specifically, our double-integral expression is identical to the "Fraunhofer diffraction integral" written down by Born &amp; Wolf at the beginning of their &sect;8.5.1; and if we follow Born &amp; Wolf's lead and position the origin of our coordinate system at the center of the rectangle, then <math>~\Theta_1 = \vartheta_1 = 0</math> and the intensity <math>~I(P)</math> at point <math>~P</math> that is given by their equation (1) precisely matches the expression for the square of the amplitude, <math>~A^*A(x_1, y_1)</math>, that is obtained from our last expression.


=See Also=
=See Also=
* Updated [[User:Tohline/Appendix/Ramblings#Computer-Generated_Holography|Table of Contents]]
* Updated [[User:Tohline/Appendix/Ramblings#Computer-Generated_Holography|Table of Contents]]
* Born, M. and Wolf, E. (1980) ''Principles of Optics,'' 3<sup>rd</sup> Edition.  New York:  Pergamon Press.  [See especially their &sect;&sect;8.5 and 8.10.]  A link to the 6<sup>th</sup> edition of this book can be found [https://www.sciencedirect.com/book/9780080264820/principles-of-optics here].




{{LSU_HBook_footer}}
{{LSU_HBook_footer}}

Latest revision as of 17:42, 10 February 2020

CGH: 2D Rectangular Apertures that are Parallel to the Image Screen

This chapter is intended primarily to replicate §I.B from the online class notes — see also an updated Table of Contents — that I developed in conjunction with a course that I taught in 1999 on the topic of Computer Generated Holography (CGH) for a subset of LSU physics majors who were interested in computational science. This discussion parallels the somewhat more detailed one presented in §I.A on the one-dimensional aperture oriented parallel to the image screen.


Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |


Utility of FFT Techniques

Consider the amplitude (and phase) of light that is incident at a location (x1, y1) on an image screen that is located a distance Z from a rectangular aperture of width w and height h. By analogy with our accompanying discussion in the context of 1D apertures, the complex number, A, representing the light amplitude and phase at (x1, y1) will be,

<math>~A(x_1, y_1)</math>

<math>~=</math>

<math>~ \sum_j \sum_k a_{jk} e^{i(2\pi D_{jk} /\lambda + \phi_{jk})} \, , </math>

where, here, the summations are taken over all "j,k" elements of light across the entire 2D aperture, and now the distance Djk is given by the expression,

<math>~D^2_{jk}</math>

<math>~\equiv</math>

<math>~ (X_j - x_1)^2 + (Y_k - y_1)^2 + Z^2 </math>

 

<math>~=</math>

<math>~ Z^2 + y_1^2 - 2y_1 Y_k + Y_k^2 + x_1^2 - 2x_1 X_j + X_j^2 </math>

 

<math>~=</math>

<math>~ L^2 \biggl[1 - \frac{2(x_1 X_j + y_1 Y_k ) }{L^2} + \frac{X_j^2 + Y_k^2}{L^2} \biggr] \, , </math>

and,

<math>~L</math>

<math>~\equiv</math>

<math>~ [Z^2 + y_1^2 + x_1^2]^{1 / 2} \, . </math>

If <math>~|X_j/L| \ll 1</math> and <math>~|Y_k/L| \ll 1</math> we can drop the quadratic terms in favor of the linear ones in the expression for Djk and deduce that,

<math>~D_{jk}</math>

<math>~\approx</math>

<math>~ L \biggl[1 - \frac{2(x_1 X_j + y_1 Y_k ) }{L^2} \biggr]^{1 / 2} </math>

 

<math>~\approx</math>

<math>~ L \biggl[1 - \frac{(x_1 X_j + y_1 Y_k ) }{L^2} \biggr] \, . </math>

Hence, the double-summation expression for the amplitude at screen location (x1, y1) becomes,

<math>~A(x_1, y_1)</math>

<math>~\approx</math>

<math>~ A_0 \sum_j \sum_k a_{jk} e^{i\phi_{jk} } \cdot \exp\biggl\{ -i \biggl[ \frac{2\pi(x_1 X_j + y_1Y_k}{\lambda L} \biggr] \biggr\} \, , </math>

where,

<math>~A_0 \equiv e^{i(2\pi L/\lambda)}</math>.

When written in this form, it should be apparent why discrete Fourier transform techniques — specifically, 2D-FFT techniques — are useful tools for evaluation of the complex amplitude, A(x1, y1).

NOTE: If x1 and/or y1 are ever comparable in size to Z — which may be the case for large apertures or for apertures tilted by nearly 90° to the image screen — then the variation of L with image screen position cannot be ignored and, accordingly, the coefficient A0 cannot be moved outside of the double summation.

Analytic Results

Now, as in §I.A, we would like to consider the case where the rectangular aperture is divided into an infinite number of divisions in both the X and Y dimensions and convert the summations in this last expression into integrals whose limits in both directions are given by the edges of the aperture. If we specifically consider the case where the aperture is assumed to be uniformly bright (i.e., ajk = a0 dX DY, and a0 is the brightness per unit area), and the phase <math>~\phi_{jk} = 0</math> at all locations on the aperture, we can write,

<math>~A(x_1, y_1)</math>

<math>~\approx</math>

<math>~ A_0 a_0 \iint \exp\biggl\{ -i \biggl[ \frac{2\pi(x_1 X_j + y_1Y_k}{\lambda L} \biggr] \biggr\} dX dY </math>

 

<math>~=</math>

<math>~ A_0 a_0 \int \exp\biggl\{ -i \biggl[ \frac{2\pi x_1 X_j }{\lambda L} \biggr] \biggr\} dX \cdot \int \exp\biggl\{ -i \biggl[ \frac{2\pi y_1Y_k}{\lambda L} \biggr] \biggr\} dY \, . </math>

Both of these integrals can be completed in the same fashion as described in §I.A for the 1D slit, giving,

<math>~A(x_1, y_1)</math>

<math>~\approx</math>

<math>~ A_0 a_0 (wh) e^{-i\Theta_1} ~ e^{-i\vartheta_1} \mathrm{sinc}(\alpha_1) \mathrm{sinc}(\beta_1) \, , </math>

where,

<math>~\alpha_1</math>

<math>~\equiv </math>

<math>~\frac{\pi x_1 (X_1 - X_2)}{\lambda L}= \frac{\pi x_1 h}{\lambda L} \, ,</math>

<math>~\Theta_1</math>

<math>~\equiv </math>

<math>~\frac{\pi x_1 (X_1 + X_2)}{\lambda L} \, ,</math>

<math>~\beta_1</math>

<math>~\equiv </math>

<math>~\frac{\pi y_1 (Y_1 - Y_2)}{\lambda L}= \frac{\pi y_1 w}{\lambda L} \, ,</math>

<math>~\vartheta_1</math>

<math>~\equiv </math>

<math>~\frac{\pi y_1 (Y_1 + Y_2)}{\lambda L} \, .</math>

It is worth noting that this derivation closely parallels the one presented in §8.5.1 (p. 393) of Born & Wolf (1980) (see the reference below). Specifically, our double-integral expression is identical to the "Fraunhofer diffraction integral" written down by Born & Wolf at the beginning of their §8.5.1; and if we follow Born & Wolf's lead and position the origin of our coordinate system at the center of the rectangle, then <math>~\Theta_1 = \vartheta_1 = 0</math> and the intensity <math>~I(P)</math> at point <math>~P</math> that is given by their equation (1) precisely matches the expression for the square of the amplitude, <math>~A^*A(x_1, y_1)</math>, that is obtained from our last expression.

See Also

  • Updated Table of Contents
  • Born, M. and Wolf, E. (1980) Principles of Optics, 3rd Edition. New York: Pergamon Press. [See especially their §§8.5 and 8.10.] A link to the 6th edition of this book can be found here.


Whitworth's (1981) Isothermal Free-Energy Surface

© 2014 - 2021 by Joel E. Tohline
|   H_Book Home   |   YouTube   |
Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS |
Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation