  We introduce a computational framework which avoids solving explicitly
hydrodynamic equations and is suitable to study the pre-merger evolution of
black hole-neutron star binary systems. The essence of the method consists of
constructing a neutron star model with a black hole companion and freezing the
internal degrees of freedom of the neutron star during the course of the
evolution of the space-time geometry. We present the main ingredients of the
framework, from the formulation of the problem to the appropriate computational
techniques to study these binary systems. In addition, we present numerical
results of the construction of initial data sets and evolutions that
demonstrate the feasibility of this approach.
