The growing use of composite materials for aerospace applications has resulted in a need for quantitative nondestructive evaluation (NDE) methods appropriate for characterizing damage in composite components. NDE simulation tools, such as ultrasound models, can aid in enabling optimized inspection methods and establishing confidence in inspection capabilities. In this paper, a mathematical approach using the Lebedev Finite Difference (LFD) method is presented for ultrasonic wave simulation in composites. Boundary condition equations for implementing stress-free boundaries (necessary for simulation of NDE scenarios) are also presented. Quantitative comparisons between LFD guided wave ultrasound simulation results, experimental guided wave data, and dispersion curves are described. Additionally, stability tests are performed to establish the LFD code behavior in the presence of stress-free boundaries and low-symmetry anisotropy. Results show that LFD is an appropriate approach for simulating ultrasound in anisotropic composite materials and that the method is stable in the presence of low-symmetry anisotropy and stress-free boundaries. Studies presented in this paper include guided wave simulation in hexagonal, monoclinic, triclinic and layered composite laminates.