This section shows how a multilayer 3D OLED with PC patterning can be simulated with FDTD using the methodology described in the 2D example The process is very similar to the 2D example from the previous section, and we will omit a lot of the detail here. Therefore, we highly recommend going through the following sections in detail before proceeding on to one of the 3D OLED examples
Even though the OLED structure in this example clearly has some symmetry, we will disregard this for now and treat this as an OLED structure with no symmetry. The files provided here are only intended to demonstrate this general approach, and should only be used when the actual OLED device has no symmetry! In the next section, we will describe how to take advantage of the lattice symmetry to reduce the number of simulations required.

We will consider a structure similar to the one propose by Chutinan et al. Initially, we assume that all materials, except the cathode, are lossless dielectrics, but this can be easily changed. Also, we will start with a simulation span of only 8x8 um2 in the x and y plane. The simulation file OLED_simple.fsp contains a parametrized OLED structure. Please note that all the children of this structure are fully deleted each time a parameter is changed, so you should not add or remove objects from this group. 
Here, we might want to consider 4x4 dipole locations across the unit cell (more may be required, but we have found this to work reasonably well).
For each dipole location, we need 3 simulations for the 3 dipole orientations (x/y/z). In principle, this means that we need 4x4x3=48 simulations (please see the next page for a discussion on how to reduce the total number of simulations). As in the simple 2D example, we use the parameter sweeping tool to automate the process of sweeping through the different dipole locations and orientations.
Note: to add more dipole locations, simply increase the "Number of points" in the pattern_dipole_position_x and pattern_dipole_position_y sweeps.
The process of calculating the radiative decay rate is identical to the approach used in Simple 2D OLED. The script OLED_3D_sweepresults.lsf can be used to run the sweeps and obtain the results.
Even though the extraction efficiency analysis can be carried out with the parameter sweep project as shown in the simple 2D OLED example (ie. by averaging the far field projections for each dipole orientation/location). Here, we will carry out this analysis with a separate script OLED_3D_farfield.lsf.
Please see Using symmetry to reduce the number of simulations and 3D OLED with square symmetry, where symmetry is applied and one can obtain reasonably similar results in a fraction of the time that it takes to run the 48 simulations required here.