Comparison of numerical methods for the calculation of two-dimensional turbulence.

*(English)*Zbl 0678.76048Summary: The estimate derived by W. D. Henshaw, the second author and L. G. Reyna for the smallest scale present in solutions of the two- dimensional incompressible Navier-Stokes equations [ICASE Rep. No.88-8, 49 pp. (1988)] is employed to obtain convergent pseudospectral approximations. These solutions are then compared with those obtained by a number of commonly used numerical methods.

If the viscosity term is deleted and the energy in high wave numbers removed by setting the amplitudes of all wave numbers above a certain point in the spectrum to zero, the “chopped” solution differs considerably from the convergent solution, even at early times. In the case that the regular viscosity is replaced by a hyperviscosity term, i.e., the square of the Laplacian, we also derive an estimate for the smallest scale present. If the coefficient of hyperviscosity is chosen so that the spectrum of the hyperviscosity solution disappears at the same point as for the regular viscosity solution, the hyperviscosity solution is also completely different from the convergent solution. If we “tune” the hyperviscosity coefficient, then the solutions are similar in amplitude or phase, but not both.

The solution obtained by a second-order difference method with twice the number of points as the pseudospectral model, or a fourth-order difference method with the same number of points as the pseudospectral model, is essentially identical to the convergent solution. This is reasonable since most of the energy of the solution is contained in the lower part of the spectrum.

If the viscosity term is deleted and the energy in high wave numbers removed by setting the amplitudes of all wave numbers above a certain point in the spectrum to zero, the “chopped” solution differs considerably from the convergent solution, even at early times. In the case that the regular viscosity is replaced by a hyperviscosity term, i.e., the square of the Laplacian, we also derive an estimate for the smallest scale present. If the coefficient of hyperviscosity is chosen so that the spectrum of the hyperviscosity solution disappears at the same point as for the regular viscosity solution, the hyperviscosity solution is also completely different from the convergent solution. If we “tune” the hyperviscosity coefficient, then the solutions are similar in amplitude or phase, but not both.

The solution obtained by a second-order difference method with twice the number of points as the pseudospectral model, or a fourth-order difference method with the same number of points as the pseudospectral model, is essentially identical to the convergent solution. This is reasonable since most of the energy of the solution is contained in the lower part of the spectrum.

##### MSC:

76F99 | Turbulence |

35B30 | Dependence of solutions to PDEs on initial and/or boundary data and/or on parameters of PDEs |

65M99 | Numerical methods for partial differential equations, initial value and time-dependent initial-boundary value problems |

65N99 | Numerical methods for partial differential equations, boundary value problems |

##### Keywords:

two-dimensional incompressible Navier-Stokes equations; convergent pseudospectral approximations; second-order difference method; pseudospectral model; fourth-order difference method
PDF
BibTeX
XML
Cite

\textit{G. L. Browning} and \textit{H. O. Kreiss}, Math. Comput. 52, No. 186, 369--388 (1989; Zbl 0678.76048)

Full Text:
DOI

##### References:

[1] | D. L. Book, NRL Plasma Formulary, Naval Research Laboratory, 1980, 60 pp. |

[2] | M. E. Brachet, M. Meneguzzi & P. L. Sulem, ”Small-scale dynamics of high-Reynolds-number two-dimensional turbulence,” Phys. Rev. Lett., v. 57, 1986, pp. 683-686. |

[3] | G. S. Deem & N. J. Zabusky, ”Vortex waves: stationary V-states, interactions, recurrence and breaking,” Phys. Rev. Lett., v. 40, 1978, pp. 859-862. |

[4] | B. Fornberg, ”A numerical study of 2-d turbulence,” J. Comput. Phys., v. 25, 1977, pp. 1-31. · Zbl 0461.76040 |

[5] | D. G. Fox & S. A. Orszag, ”Pseudospectral approximation to two-dimensional turbulence,” J. Comput. Phys., v. 11, 1973, pp. 612-619. · Zbl 0261.76041 |

[6] | W. D. Henshaw, H.-O. Kreiss & L. G. Reyna, On the Smallest Scale for the Incompressible Navier-Stokes Equations, ICASE Report No. 88-8, 1988, 49 pp. · Zbl 0708.76043 |

[7] | J. R. Herring, S. A. Orszag, R. H. Kraichnan & D. G. Fox, ”Decay of two-dimensional turbulence,” J. Fluid Mech., v. 66, 1974, pp. 417-444. · Zbl 0296.76026 |

[8] | Heinz-Otto Kreiss and Joseph Oliger, Comparison of accurate methods for the integration of hyperbolic equations, Tellus 24 (1972), 199 – 215 (English, with Russian summary). |

[9] | D. K. Lilly, ”Numerical simulation of developing and decaying of two-dimensional turbulence,” J. Fluid Mech., v. 45, 1971, pp. 395-415. · Zbl 0242.76027 |

[10] | J. C. McWilliams, ”The emergence of isolated coherent vortices in turbulent flow,” J. Fluid Mech., v. 146, 1984, pp. 21-43. · Zbl 0561.76059 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.