  In order to understand the formation of the multipolar structures of the
pre-planetary nebula (PPN) CRL 618, we perform 3D simulations using a
multi-directional bullet model. The optical lobes of CRL 618 and fast molecular
outflows at the tips of the lobes have been found to have similar expansion
ages of ~ 100 yr. Additional fast molecular outflows were found near the source
along the outflow axes with ages of ~ 45 yr, suggesting a second episode of
bullet ejections. Thus, in our simulations, two episodes of bullet ejections
are assumed. The shaping process is simulated using the ZEUS-3D hydrodynamics
code that includes molecular and atomic cooling. In addition, molecular
chemistry is also included to calculate the CO intensity maps. Our results show
the following: (1) Multi-epoch bullets interacting with the toroidal dense core
can produce the collimated multiple lobes as seen in CRL 618. The total mass of
the bullets is ~ 0.034 solar mass, consistent with the observed high-velocity
CO emission in fast molecular outflows. (2) The simulated CO J=3-2 intensity
maps show that the low-velocity cavity wall and the high-velocity outflows
along the lobes are reasonably consistent with the observations. The
position-velocity diagram of the outflows along the outflow axes shows a linear
increase of velocity with distance, similar to the observations. The ejections
of these bullets could be due to magneto-rotational explosions or nova-like
explosions around a binary companion.
