Neutron-induced reactions with charged particle emission play an important role in a variety of research fields ranging from fundamental nuclear physics and nuclear astrophysics, to applications of nuclear technologies, to energy production, and material science. Recently, the capability to study reactions with radioactive targets has become important to significantly advance research in explosive nucleosynthesis and nuclear applications. To achieve the relevant research goals and study (n, x) reactions over a broad neutron beam energy range, the Low Energy Neutron-induced charged-particle (Z) chamber (LENZ) at Los Alamos Neutron Science Center (LANSCE) was developed, along with varied ancillary instrumentation, to enable the aforementioned research program. For the (n, x) reactions of interest at low energies, a precise simulation of the discrete spectrum of emitted charged particles is essential. In addition, since LANSCE is a user facility, a simulation of the LENZ setup easily accessible by users has a high value. With these goals in mind, we have developed a detailed simulation using the Geant4 toolkit. In this work, we present the implementation and the validation of the simulation using experimental data from recent campaigns with the LENZ apparatus. Specifically, we benchmark the simulation against a similar MCNP-based application, and determine the realistic range of applicability for the utilized probability biasing technique. We describe the implementation of an evaluated library with angular distribution and partial cross-section data, and we perform a validation of the simulation based on comparisons of simulated spectra with experimental ones for a number of targets used in previous experimental campaigns. Lastly, we discuss the limitations, caveats, and assets of the simulation code and techniques used.