A dynamic approach for evaluating parameters in a numerical method

A new methodology for evaluating unknown parameters in a numerical method for solving a partial differential equation is developed. The main result is the identification of a functional form for the parameters which is derived by requiring the numerical method to yield ‘optimal’ solutions over a set of finite-dimensional function spaces. The functional depends upon the numerical solution, the forcing function, the set of function spaces, and the definition of the optimal solution. It does not require exact or approximate analytical solutions of the continuous problem, and is derived from an extension of the variational Germano identity. This methodology is applied to the one-dimensional, linear advection–diffusion problem to yield a non-linear dynamic diffusivity method. It is found that this method yields results that are commensurate to the SUPG method. The same methodology is then used to evaluate the Smagorinsky eddy viscosity for the large eddy simulation of the decay of homogeneous isotropic turbulence in three dimensions. In this case the resulting method is found to be more accurate than the constant-coefficient and the traditional dynamic versions of the Smagorinsky model