This paper describes the development and testing of a 3D finite differ
ence code written specifically to model turbulence in an open channel
with a moving free surface. The code has been developed so that either
a full simulation or a large eddy simulation (LES) of the turbulence
may be performed. The free surface may undergo arbitrarily large defor
mations but the slope may not exceed a limit related to the aspect rat
io of the mesh and so the possibility of breaking waves is excluded. T
he LES application demands numerical approximations which conserve mas
s, momentum, and total energy with high precision, and it permits wave
motion with very little numerical dispersion or dissipation. We descr
ibe a novel numerical method for tracking the free surface using a spl
it-merge technique which combines the volume of fluid and height funct
ion methods in a way that is conservative. (C) 1995 Academic Press, In
c.