Josephson traveling-wave parametric amplifiers (JTWPAs) are critical for superconducting quantum computers in enabling simultaneous high-fidelity measurements of many qubits. Despite their importance, translating designs to working hardware remains a significant obstacle, with numerical modeling increasingly being used to help close this gap. Unfortunately, existing time-domain (TD) circuit modeling methods struggle to close this gap due to their inability to model the physical layout of a device. In principle, finite-element time-domain (FETD) methods can model the JTWPA physical layout, but these simulations are complicated by the wide bandwidths and multiscale features involved. Here, we develop an unconditionally stable 1-D FETD model for JTWPAs that can model dispersive losses over wide bandwidths to begin overcoming these challenges. We validate our method against harmonic balance (HB) simulations of a Floquet-mode JTWPA before simulating a realistic multiplexed qubit readout scenario where HB is impractical.