  We introduce a new Eulerian simulation framework for liquid animation that
leverages both finite element and finite volume methods. In contrast to
previous methods where the whole simulation domain is discretized either using
the finite volume method or finite element method, our method spatially merges
them together using two types of discretization being tightly coupled on its
seams while enforcing second order accurate boundary conditions at free
surfaces. We achieve our formulation via a variational form using new shape
functions specifically designed for this purpose. By enabling a mixture of the
two methods, we can take advantage of the best of two worlds. For example,
finite volume method (FVM) result in sparse linear systems; however, complexity
is encountered when unstructured grids such as tetrahedral or Voronoi elements
are used. Finite element method (FEM), on the other hand, result in comparably
denser linear systems, but the complexity remains the same even if unstructured
elements are chosen; thereby facilitating spatial adaptivity. In this paper, we
propose to use FVM for the majority parts to retain the sparsity of linear
systems and FEM for parts where the grid elements are allowed to be freely
deformed. An example of this application is locally moving grids. We show that
by adapting the local grid movement to an underlying nearly rigid motion,
numerical diffusion is noticeably reduced; leading to better preservation of
structural details such as sharp edges, thin sheets and spindles of liquids.
