We present the Quantum Computer Aided Design (QCAD) simulator that targets modeling multi-dimensional quantum devices, particularly silicon multi-quantum dots (QDs) developed for quantum bits (qubits). This finite-element simulator has three differentiating features: (i) its core contains nonlinear Poisson, effective mass Schrodinger, and Configuration Interaction solvers that have massively parallel capability for high simulation throughput, and can be run individually or combined self-consistently for 1D/2D/3D quantum devices; (ii) the core solvers show superior convergence even at near-zero-Kelvin temperatures, which is critical for modeling quantum computing devices; and (iii) it interfaces directly with the full-featured optimization engine Dakota. In this work, we describe the capabilities and implementation of the QCAD simulation tool, and show how it can be used to both analyze existing experimental QD devices through capacitance calculations, and aid in the design of few-electron multi-QDs. In particular, we observe that computed capacitances are in rough agreement with experiment, and that quantum confinement increases capacitance when the number of electrons is fixed in a quantum dot. Coupling of QCAD with the optimizer Dakota allows for rapid identification and improvement of device layouts that are likely to exhibit few-electron quantum dot characteristics. FIG. 1. (Color online) Schematic diagram of the QCAD code structure.