The strong suppression of bottomonia production in ultra-relativistic heavy-ion collisions is a smoking gun for the creation of a deconfined quarkgluon plasma (QGP). In this proceedings contribution, I review recent work that aims to provide a more comprehensive and systematic understanding of bottomonium dynamics in the QGP through the use of pNRQCD and an open quantum systems approach. This approach allows one to evolve the heavyquarkonium reduced density matrix, taking into account non-unitary effective Hamiltonian evolution of the wave-function and quantum jumps between different angular momentum and color states. In the case of a strong coupled QGP in which Ebind ≪ T, mD ≪ 1=a0, the corresponding evolution equation is Markovian and can therefore be mapped to a Lindblad evolution equation. To solve the resulting Lindblad equation, we make use of a stochastic unraveling called the quantum trajectories algorithm and couple the non-abelian quantum evolution to a realistic 3+1D viscous hydrodynamical background. Using a large number of Monte-Carlo sampled bottomonium trajectories, we make predictions for bottomonium RAA and elliptic flow as a function of centrality and transverse momentum and compare to data collected by the ALICE, ATLAS, and CMS collaborations.