We created a simulation model to investigate the characteristics of fluorescence in two-photon-excited samples. In the model, the sample is a diffusible solution of fluorophore molecules, which is divided into cubic cells and illuminated by a train of focused laser pulses described as a Gaussian beam. Simulating the state transitions according to a multilevel photodynamic model (also including photobleaching and intersystem crossing), the simulator provides the expected number and the spatial distribution of emitted photons over time. Our simulations demonstrated how the illumination laser power, diffusion, and the photodynamic parameters of the fluorophore affect fluorescence. We revealed the unusual fluorescent profile that evolves as photobleaching progresses: the most photons are not emitted from the focus (where a "dark hole" appears) but from an ellipsoid around the focus. The model could be adapted to several fluorescent techniques (such as two-photon microscopy and fluorescence recovery after photobleaching). Furthermore, it might help to optimize the operating parameters of the measurement devices (e.g., in order to reach higher image quality and lower photobleaching).