There is an increasing interest in low voltage direct current (LVDC) distribution grids due to advancements in power electronics enabling efficient and economical electrical networks in the DC paradigm. Power flow equations in LVDC grids are non-linear and non-convex due to the presence of constant power nodes. Depending on the implementation, power flow equations may lead to more than one solution and unrealistic solutions; therefore, the uniqueness of the solution should not be taken for granted. This paper proposes a new power flow solver based on a graph theory for LVDC grids having radial or meshed configurations. The solver provides a unique solution. Two test feeders composed of 33 nodes and 69 nodes are considered to validate the effectiveness of the proposed method. The proposed method is compared with a fixed-point methodology called direct load flow (DLF) having a mathematical formulation equivalent to a backward forward sweep (BFS) class of solvers in the case of radial distribution networks but that can handle meshed networks more easily thanks to the use of connectivity matrices. In addition, the convergence and uniqueness of the solution is demonstrated using a Banach fixed-point theorem. The performance of the proposed method is tested for different loading conditions. The results show that the proposed method is robust and has fast convergence characteristics even with high loading conditions. All simulations are carried out in MATLAB 2020b software.