DC smart grids are a promising solution for the efficient integration of renewable energy sources and loads. Still, their widespread adoption is hindered by significant challenges related to fault response, identification, and clearance. The traditional DC fault analysis method is a useful tool for straightforwardly understanding the behaviour of fault current contributions from DC converters in LVDC networks during a fault. However, when a system with multiple converters and non-negligible fault impedance need to be considered, its accuracy is severely limited due to the assumptions included in the problem solution, thus leading to the following: (a) the dependency of the results’ reliability on fault impedance values and/or other converter fault current contributions; (b) the inaccuracy of the diode current estimation; and (c) the inaccuracy of the conductor joule integral. Thus, these results’ data may be unreliable for designing protection systems for one converter or for an entire network. In order to overcome these issues, this paper proposes an innovative, simple numerical approach to DC fault current evaluation, which can be adopted when the number of converters become significant, or the network is complex. This method arises from the primary interest in solving the circuit to extract the indicators (current peak value and time, joule integral, etc.) necessary for designing circuit protections. This approach proved to grant two main advantages over traditional methods: (a) it provides accurate results, with no need to introduce any specific assumption; (b) it can be structured to manage an arbitrary number of converters; and (c) it reduces the computational processing times and resources necessary to simulate an entire DC network in comparison to other circuit solution software.