When plasmas come in contact with the boundaries that confine them, various complex processes occur between the plasma and the materials in the boundary. These processes, called plasma-material interactions (PMI) lead to physical and chemical modifications in the materials and in the plasma. In the case of a tokamak, a magnetic confinement fusion reactor, the interactions between the plasma and the material in the bounding walls can negatively impact the performance and service life of the reactor. Furthermore, PMI are also found in other areas of significant engineering interest, such as plasma-based spacecraft propulsion engines, where interactions affect the transport properties of the plasma and consequently the performance of the engine. Therefore, gaining a fundamental understanding of the various plasma-material interactions is necessary for the development and improvement of these devices. PMI are dictated by the plasma sheath, a layer of net positive charge that forms at the plasma-boundary interface. The sheath regulates the energy and particle fluxes to the boundary, mediating the interactions. Sheaths, however, are only stable and well-developed when the ions enter the sheath with a speed equal to or greater than the `Bohm speed'. The Bohm speed is a landmark result in sheath theory and various mathematical expressions for it have been derived from fluid and kinetic treatment of plasmas. Although these models are widely used, they are only accurate in cases where the thickness of the sheath is negligible when compared to the scale length of the plasma in consideration. These cases are said to satisfy the `asymptotic limit'. To resolve this, a new Bohm speed model that considers the effects of transport terms such as the electron heat flux, thermal force, and temperature isotropization has been recently proposed [Y. Li et al., Physical Review Letters (2022)]. The model is verified using particle-in-cell (PIC) kinetic simulations and is shown to accurately predict the Bohm speed in cases away from the asymptotic limit. This thesis investigates the new model using the continuum kinetic approach on the Gkeyll software framework. The continuum kinetic approach numerically solves the Vlasov-Maxwell equations using the discontinuous Galerkin method and captures the kinetic phenomena of the plasma without needing to track individual particles. Multiple collisional cases ranging from a Knudsen number of 20 to 5000 are considered in a 1X3V simulation domain using the Lenard-Bernstein collisional operator. The results of the continuum kinetic simulations are benchmarked to the PIC simulation results. It is concluded that across a wide range of collisionalities, the continuum kinetic method captures much of the same physics as the PIC method while offering noise-free results. However, there is a discrepancy between the Bohm speed prediction and the simulation results in the continuum kinetic case. This discrepancy is explored and significant error in the collisional integral derived transport terms between the continuum kinetic method and PIC method is found, suggesting that the difference in collisional operator may be the source of the discrepancy. Nevertheless, the sheath profiles developed in the PIC simulations and the continuum kinetic simulations are in reasonable agreement. / Master of Science / Nuclear fusion is a process in which two light atomic nuclei (like hydrogen) fuse to form a heavier nucleus (like helium) and release tremendous amounts of energy. The resultant energy from these reactions powers the sun and has the potential to provide clean energy for our terrestrial needs. Harnessing fusion energy has been one of the biggest scientific and engineering challenges of our time due to various reasons. One of these reasons is the interaction of plasma, which is the fuel for the fusion reaction, and the materials of the walls that bound the plasma. These interactions are called plasma-material interactions (PMI) and can affect the longevity and performance of fusion reactors. The main governing phenomenon behind these interactions is the plasma sheath, a layer of plasma that is formed when the plasma encounters a boundary. For a sheath to form it is also necessary that ions in the plasma possess a speed greater than the so-called `Bohm speed'. While many expressions have been derived for the Bohm speed, these expressions are only valid when there is a clear sheath entrance that divides the bulk plasma and the sheath. This condition is not satisfied in many cases of interest. Instead, a sheath-transition region is found to exist between the bulk plasma and the sheath. A recently proposed Bohm speed model [Y. Li et al., Physical Review Letters (2022)] resolves this. This model is accurate in cases where the sheath-transition region exists and is derived by considering previously overlooked transport physics. In this work, this new model is studied using a different computational approach known as `continuum kinetics' using an open-source solver called Gkeyll. The results of the continuum kinetic simulations are compared to the results used to verify the model.