A Cartesian grid generation methodology is developed for unsteady control volume computational fluid dynamic (CFD) solvers. Arbitrary combinations and numbers of moving and stationary boundaries are allowed to exist in the two-dimensional Cartesian grid template. Specific definitions for the possible cases, CFD solver requirements and moving geometry handling algorithms of Cartesian grids are described. Applications are selected from bio-fluid dynamics and aerospace propulsion to demonstrate the capability of the method.