This project has moved. For the latest updates, please go here.

Heterogeneous Fluence

Apr 5, 2013 at 11:55 PM
Can a layered system be used for the "FluenceofRhoAndZ" function in MATLAB? I have the basic cases working but would like to try layers.

I would also like to know what options (if any) are around for fluence calculations in gradients - say optical properties that vary according to some determined function in space.

Apr 6, 2013 at 7:16 PM
Hi Albert,

Yes, you can use a layered system with FluenceOfRhoAndZ. Any axi-symmetric tissue definition is fine.

No, we don't have the capability to define optical properties that vary according to some function in space. Sorry.

Apr 6, 2013 at 7:21 PM
great - how do you program it? I noticed that in the comments you can enter an array for the different optical properties but how to fix their position?

Im not that clever so Ill need it spelled out ...

Apr 6, 2013 at 8:30 PM
You specify a MultiLayerTissue with 2 or more tissue layers in the infile. There are a couple of sample infiles provided with the Monte Carlo CommandLine download that you can start with. For example infile_one_layer_FluenceOfRhoAndZ_RadianceOfRhoAndZAndAngle.xml specifies a 1 layer tissue model with detectors FluenceOfRhoAndZ and RadianceOfRhoAndZAndAngle (you can delete the Radiance detector). Note that the MultiLayerTissue for a single layer of tissue is comprised of layers: air-tissue-air. You want to specify air-tissueLayer1-tissueLayer2...-air. You can refer to this page towards the bottom describes how to add a tissue layer.
Apr 7, 2013 at 4:48 AM
Just to clarify further:
  • the FluenceOfRhoAndZ function in VtsSolvers.m (located in the /vts folder of the Matlab package here) mirrors the functionality in the Fluence/Interrogation Solver Panel in the online tool, and only supports the diffusion solvers at this point.
  • there are examples of how to call this function in the downloadable Matlab tools in the file vts_solver_demo.m, (located here in our source tree)
  • the FluenceOfRhoAndZ detector is an option in the classic Monte Carlo tools, and supports all of the heterogeneous capabilities - one or more non-overlapping inclusions (ellipsoids/spheres, cyllinders, voxels), and layers.
  • this latter tool is available in two ways via Matlab: 1) run the command line tool and use the provided Matlab scripts to load/visualize, or 2) run Monte Carlo simulations directly from VtsMonteCarlo.m - there are examples in the downloadable Matlab tools in the files vts_mc_demo.m and /vts_tests/mc_tests.m (and can be viewed here and here. respectively in our source tree)
Apr 7, 2013 at 2:02 PM
thanks - ill give a try