7. Conclusion and discussion
In this paper, we have systematically presented a novel and versatile method to address a suite of research challenges encountered in modal reduction based real-time deformation and arbitrary cutting simulation of heterogeneous objects (with multiple sub-domains and large variations of material distribution). The most critical idea of our novel approach is to conduct a dynamic interchange between adaptively integrating material-aware and/or geometry-structure-aware simulation with full-physics capability and performing deformation reconstruction based on sub-domain-specific local modal’ reuse, and all of the above numerical procedures have been implemented in a CUDA-centric parallel computation framework. The novel technical components within our new framework include: the space–time-varying local modal generation from previous-cycles’ fully-physical simulation, the adaptive alternation scheme between sub-domain physical simulation and modal reuse, the cross-domain coupling of spatially-varying numerous deformations (i.e., some sub-domains’ deformation is from physical simulation, and others are from modal-subspace reconstruction), and the sub-domain-level parallel implicit integration solvers supporting CUDA-enabled numerical computation, which collectively equip our method with remarkable advantages in terms of realtime efficiency, high-accuracy simulation, unconditional stableness, and practical versatility. Currently, our prototype system is still a proof-of-concept only at the experimental stage, hence, it is not yet of practical use due to certain limitations. For example, we should introduce some collision detection function (Li et al., 2014b) into our current system to make it possible to support more complex interactions far beyond the simple cutting operations. And we should also incorporate certain mature techniques (Barbicˇ and James, 2010; Teng et al., 2014) to accommodate self-collision caused by elastic deformation. In addition, our ongoing research efforts are concentrated on seeking an efficient optimization method to guarantee the physical accuracy in an absolute sense. Meanwhile, it also deserves our efforts to extend our method to handle more sophisticated fluid-elasticity coupling phenomena. In terms of limitations, it should be noted that, our domain coupling method’s requirement for construction consistency of sub-domain interface may not be appropriate in practical applications and our current method is lack of flexibility in supporting more comprehensive adaptivity when handling much more complex boundary interfaces.