Numerical analysis of a higher order time relaxation model of fluids
Vincent J. Ervin
1, William J. Layton
2, Monika Neda
2 1 Department of Mathematical Sciences, Clemson University, Clemson, SC, 29634-0975, USA
2 Department of Mathematics, University of Pittsburgh, Pittsburgh, PA, 15260, USA
Received by the editors February 20, 2006
Abstract
We study the numerical errors in finite elements
discretizations of a time relaxation model of fluid motion:
u(t) + u center dot del u + del(p) - nu del u + chi u* = f and del
center dot u = 0
In this model, introduced by Stolz, Adams, and Kleisser, u* is a
generalized fluctuation and chi the time relaxation parameter. The goal
of inclusion of the chi u* is to drive unresolved fluctuations to sero
exponentially. We study convergence of discretization of model to the
model's solution as h, Delta t -> 0. Next we complement this with an
experimental study of the effect the time relaxation term (and a
nonlinear extension of it) has on the large scales of a flow near a
transitional point. We close by showing that the time relaxation term
does not alter shock speeds in the inviscid, compressible case, giving
analytical confirmation of a result of Stolz, Adams, and Kleiser.
AMS subject classifications: 65N30
Key words: time relaxation; deconvolution; turbulence
Email: vjervin@clemson.edu (V. J. Ervin), wjl@pitt.edu (W. J. Layton), mon5@pitt.edu (M. Neda)
|