Slowly chirped two-dimensional photonic crystal cavities are promising devices for creating photonic Bose-Einstein condensates. Before experimentally achieving such a condensate, one first has to thoroughly investigate the electromagnetic eigenmodes in such crystals. However, slowly chirped photonic crystals leading to cavities for light will easily have sizes in the order of tens of micrometers. Therefore simulating the behaviour of light in these crystals is very time consuming. In this thesis we demonstrate a novel and intuitive approach to obtain the envelopes of the electromagnetic modes in these crystals. An enormous advantage of this approach is that it can calculate the energies and the envelopes of these eigenmodes to a high accuracy in a few seconds. We model a chirped photonic crystal using a local density approach; we assign a potential energy for light, extracted from photonic bandstructure calculations, to each unit cell of the crystal. We also obtain an effective mass from the curvature of the photonic band at this energy. With these ingredients we are left with the task of solving the corresponding Schrödinger equation, which is an elegant and far less time-consuming exercise than calculating the envelopes of the electromagnetic modes using finite-difference time-domain simulations. In this thesis it is shown that for one-and two-dimensional quadratically chirped photonic crystals the agreement between the envelopes obtained by these simulations and the analytical ones resulting from this model is larger than 90% for the lowest energy eigenmodes. Moreover, even with small distortions on this quadratic behaviour of the widths of a onedimensional chirped photonic crystal, the corresponding numerically obtained solutions of the Schrödinger equation clearly have an excellent agreement with the simulated mode envelopes.