Abstract: We present the first spatially adaptive Eulerian fluid animation method to support challenging viscous liquid effects such as folding, coiling, and variable viscosity. We propose a tetrahedral node-based embedded finite volume method for fluid viscosity, adapted from popular techniques for Lagrangian deformable objects. Applied in an Eulerian fashion with implicit integration, this scheme stably and efficiently supports high viscosity fluids while yielding symmetric positive definite linear systems. To integrate this scheme into standard tetrahedral mesh-based fluid simulators, which store normal velocities on faces rather than velocity vectors at nodes, we offer two methods to reconcile these representations. The first incorporates a mapping between different degrees of freedom into the viscosity solve itself. The second uses a FLIP-like approach to transfer velocity data between nodes and faces before and after the linear solve. The former offers tighter coupling by enabling the linear solver to act directly on the face velocities of the staggered mesh, while the latter provides a sparser linear system and a simpler implementation. We demonstrate the effectiveness of our approach with animations of spatially varying viscosity, realistic rotational motion, and viscous liquid buckling and coiling.