Earth's deep mantle dynamics are closely linked to the evolution of Large Low Shear Velocity Provinces (LLSVPs) at the 2900 km deep Core-Mantle Boundary (CMB). These dynamics drive large convective mantle cells that underpin Earth's plate tectonics. Our study explores the formation and deformation of LLSVPs through geodynamic modeling and seismic anisotropy analysis. Seismic anisotropy, indicating wave propagation varying with direction, occurs in the lowermost few hundred kilometers of the mantle, especially near LLSVP edges. High anisotropy in these regions may suggest significant deformation from mantle flow and the generation of mantle plumes. Using compressible 3-D global mantle convection models in ASPECT, we examine the flow behavior of the lower mantle and simulate seismic anisotropy with ECOMAN. Our results indicate that the lower mantle is generally isotropic, except in regions near LLSVP margins, where vertically polarized shear waves (Vsv) are faster due to high finite strain, while other areas are characterized by faster horizontally polarized shear waves (Vsh). Simulations are running to test the robustness of these findings and assess the impact of various factors, including different LLSVP densities, viscosities, and temperatures. Comparing simulated and observed shear wave radial anisotropy will enhance our understanding of the lower mantle's rheology and dynamics and how large scale structures in the lower mantle control the Earth's thermo-chemical evolution.