  We develop a special phase field/diffusive interface method to model the
nuclear architecture reorganization process. In particular, we use a Lagrange
multiplier approach in the phase field model to preserve the specific physical
and geometrical constraints for the biological events. We develop several
efficient and robust linear and weakly nonlinear schemes for this new model. To
validate the model and numerical methods, we present ample numerical
simulations which in particular reproduce several processes of nuclear
architecture reorganization from the experiment literature.
