My research in lattice gauge theory during the past year is described. Several projects were completed dealing with QCD simulations including dynamical fermions. Under the DOE Grand Challenge program, a large scale calculation of the QCD spectrum with two light flavors of dynamical staggered quarks was carried out. This calculation is one of the most significant efforts to data to take into account the effects of dynamical fermions. Smaller lattice spacing and lighter quark masses were used than in previous attempts. QCD thermodynamics was studied on the ST-100 array processor and on an ETA supercomputer at the John von Neumann Supercomputer Center. On the ST-100, a study with two flavors of dynamical staggered quarks with am{sub q} = 0.025 and 0.0125 was carried out on a 12{sup 3} {times} 8 lattice. These results give a rough estimate of the crossover couplings where we see the restoration of chiral symmetry. A study of QCD with dynamical Wilson fermions was carried out with N{sub t} = 4 to try to bring the study of QCD with dynamical Wilson fermions to the level that has been attained with staggered fermions over the past two years. We have calculated screening lengths to elucidate the properties of the high temperature phase. In the pure gluon theory, claims that the finite temperature deconfinement transition is second order, rather than first order, were investigated using a finite size scaling analysis. Our results support a first order transition. Finally, work was done to port computer code to new environments involving parallelism in order to pursue more ambitious calculations on more powerful hardware than the ST-100 and ETA10 used for the calculations reported here.