The simplest method is indeed based on plotting of zero isolines (or zero iso-surfaces) of a level set function. Some further details and extensive discussion of the problems related to implementation of the level set method is presented at
I would like to draw your attention to the reinitialization procedure discussed in section 3.2. This procedure is crucial to obtain a high fidelity rendering of the boundary.
I hope that this is useful in resolving your question.