This work aims to develop a highly flexible platform for 3D SNR hydro and cosmic ray (CR) simulation, which is easily expandable to accept new physics, and adaptable to the complex environments of individual supernova remnants. As a preliminary work in progress, the simulation is tested under toy model configurations, in a cubic mesh with a relatively low resolution (21x21x21 binning). Several assumptions and simplifications are made in the physics of cosmic ray acceleration and diffusion. Once the framework becomes mature, however, we can immediately proceed to perform realistic modeling of any given SNR, match in broad-band with observations, and extract useful information.