GOALS: Gain familiarity with the use of Monte Carlo Simulations to provide fluence and radiance predictions
Bring up the GUI
I. Monte Carlo Simulations of Narrow Collimated Beam Irradiation
Goal: This portion of the lab is to examine the impact of optical properties on the amplitude and axial/lateral dispersion of the light in turbid tissues.
- Select the Monte Carlo Solver Panel.
- Click Download Prototype Input Files button. This download sample infiles in a zip file. Extract the contents of the zip file.
- In the Input File Specification, click the Load Input File button. Select the file: infile_one_layer_ROfRho_FluenceOfRhoAndZ.txt. This is a Monte Carlo simulation input file that specifies a Discrete Absorption Weighting (DAW) simulation to produce reflectance as a function of source-detector (ρ) separation, R(ρ) and fluence in cylindrical coordinates (ρ,z).
- Set the optical properties of the tissue layer to be μa = 0.01mm-1, μ's = 1mm-1, g=0.8.
- Set Number of Photons to 1000.
- Click the Run Simulation button.
- On the Map View plot, note the Mean Depth magnitude, axial and lateral penetration depths of the fluence distribution.
- Increase the absorption coefficient to μa = 1mm-1, keeping all other properties constant.
- Click the Run Monte Carlo Simulation button.
- Note the magnitude of the Mean Depth, axial and lateral penetration depths of the fluence distribution.
- Now repeat this exercise keeping μa constant at 0.1mm-1 and μ's values of 0 and 1mm-1, noting the mean depth of penetration.
- How does increasing absorption change the magnitude, axial and lateral penetration depths?
- How does increasing scattering change the magnitude, axial and lateral penetration depths?
II. Visualizing Radiance versus Fluence distributions using the Monte Carlo CommandLine Application
Goal: This portion of the lab exercise is to examine the directional aspect of radiance in turbid tissues.
- On the Desktop, holding down the shift key, right click the MCCL folder and select Open command window here.
- Execute the Monte Carlo Command Line Application using the command is mc infile=infile_one_layer_FluenceOfRhoAndZ_RadianceOfRhoAndZAndAngle.txt. This infile specifies a perpendicular point source impinging on the tissue at (0,0,0). The tissue is homogeneous with optical properties μa = 0.01mm-1, μ's = 1mm-1, g=0.8 and n = 1.4. N=10,000 photons are launched. Executing this command creates a folder called "one_layer_FluenceOfRhoAndZ_RadianceOfRhoAndZAndAngle" in the MCCL directory with fluence and radiance detector results.
- Go into the MCCL folder.
- Click on load_results_script. This will bring up MATLAB and also a window to edit load_results_script.m. Change datanames to be one_layer_FluenceOfRhoAndZ_RadianceOfRhoAndZAndAngle. Enter load_results_script at the MATLAB prompt.
- Three figures should appear, one plot of Fluence, and two plots of radiance in the lower and upper hemispheres. Note that the lower hemisphere is notated as [0-pi/2] and the upper hemisphere as [pi/2-pi]. These angles are with reference to the positive z-axis which is into the tissue.
- In what way can the angular information shown in the radiance plots be useful?
- Which plot of radiance looks most like Fluence and why?
Additional Questions: Consider the situation where you using a laser-based therapy to treat an embedded tumor. For this application, it is critical that you maximize the axial penetration of the light field. However during the treatment, the tissue absorption may increase due to increased blood flow and scattering may decrease due to morphological changes in the tissue. In this context comment on the following:
- How does an increase in absorption impact (a) the axial penetration of the light field and (b) the lateral dispersion of the light field?
- How does a decrease in scattering impact (a) the axial penetration of the light field and (b) the lateral dispersion of the light field?
- Would increasing the beam diameter improve (a) the axial penetration of the light field and/or (b) the lateral dispersion of the light field? Is this impact the same regardless of the optical properties?