We demonstrate in this paper how to use direct modeling of heat transfer and circuital equations based on "packaged" simulation tools to produce a model suitable for the study of quench propagation in a system of magnetically coupled solenoids protected without an external dump resistor. Our application example is a system of three, layer-wound and bath-cooled coils, all using a NbTi conductor. The same model can be easily adapted to study coils with different topologies, geometry and conductors. We demonstrate the use of the model by parametric analysis of physical parameters in the system, showing under which conditions a quench-back due to inductive coupling is triggered. The results, and in particular longitudinal and transverse propagation velocities, are consistent with expected analytical scalings. The advantage of the model is that it provides self-consistent results, i.e. based on first principles rather than assumptions on quench propagation, and it gives access to local details that are not intuitive, nor easily measurable. (C) 2012 Elsevier Ltd. All rights reserved.