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. Copyright © 2004 John Wiley & Sons, Ltd.

[1]  Meinhard E. Mayer,et al.  Navier-Stokes Equations and Turbulence , 2008 .

[2]  Roger Temam,et al.  Navier-Stokes Equations and Turbulence by C. Foias , 2001 .

[3]  Assad A. Oberai,et al.  Variational formulation of the Germano identity for the Navier–Stokes equations , 2005 .

[4]  Nagi N. Mansour,et al.  Decay of Isotropic Turbulence at Low Reynolds Number , 1994 .

[5]  Endre Süli,et al.  Modeling subgrid viscosity for advection–diffusion problems , 2000 .

[6]  D. Lilly,et al.  A proposed modification of the Germano subgrid‐scale closure method , 1992 .

[7]  Eugenio Oñate,et al.  Derivation of stabilized equations for numerical solution of advective-diffusive transport and fluid flow problems , 1998 .

[8]  Franco Brezzi,et al.  $b=\int g$ , 1997 .

[9]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[10]  P. Moin,et al.  A dynamic subgrid‐scale eddy viscosity model , 1990 .

[11]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[12]  Paul E. Barbone,et al.  Nearly H1-optimal finite element methods , 2001 .

[13]  T. Hughes,et al.  The variational multiscale method—a paradigm for computational mechanics , 1998 .

[14]  K. Lilly.  On the application of the eddy viscosity concept in the Inertial sub-range of turbulence , 1966 .

[15]  P. Moin,et al.  A dynamic localization model for large-eddy simulation of turbulent flows , 1995, Journal of Fluid Mechanics.

[16]  T. Hughes.  Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods , 1995 .

[17]  R. Rogallo.  Numerical experiments in homogeneous turbulence , 1981 .