Analysis of pseudo-spectral methods used for numerical simulation of turbulence

09/01/2021
by   Tapan K. Sengupta, et al.
0

Global spectral analysis (GSA) is used as a tool to test the accuracy of numerical methods with the help of canonical problems of convection and convection-diffusion equation which admit exact solutions. Similarly, events in turbulent flows computed by direct numerical simulation (DNS) are often calibrated with theoretical results of homogeneous isotropic turbulence due to Kolmogorov, as given in Turbulence -U. Frisch, Cambridge Univ. Press, UK (1995). However, numerical methods for the simulation of this problem are not calibrated, as by using GSA of convection and/or convection-diffusion equation. This is with the exception in "A critical assessment of simulations for transitional and turbulence flows-Sengupta, T.K., In Proc. of IUTAM Symp. on Advances in Computation, Modeling and Control of Transitional and Turbulent Flows, pp 491-532, World Sci. Publ. Co. Pte. Ltd., Singapore (2016)", where such a calibration has been advocated with the help of convection equation. For turbulent flows, an extreme event is characterized by the presence of length scales smaller than the Kolmogorov length scale, a heuristic limit for the largest wavenumber present without being converted to heat. With growing computer power, recently many simulations have been reported using a pseudo-spectral method, with spatial discretization performed in Fourier spectral space and a two-stage, Runge-Kutta (RK2) method for time discretization. But no analyses are reported to ensure high accuracy of such simulations. Here, an analysis is reported for few multi-stage Runge-Kutta methods in the Fourier spectral framework for convection and convection-diffusion equations. We identify the major source of error for the RK2-Fourier spectral method using GSA and also show how to avoid this error and specify numerical parameters for achieving highest accuracy possible to capture extreme events in turbulent flows.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
08/21/2022

Determining Parameter Ranges for High Accuracy Large Eddy Simulation by Lax-Wendroff Method

The analysis of Lax-Wendroff (LW) method is performed by the generic mod...
research
08/19/2022

Quantifying Resolutions for DNS and LES for Lax-Wendroff Method: Application to Uniform/Non-Uniform Compact Schemes

The global spectral analysis (GSA) of numerical methods ensures that the...
research
03/07/2020

Quantum turbulence simulations using the Gross-Pitaevskii equation: high-performance computing and new numerical benchmarks

This paper is concerned with the numerical investigation of Quantum Turb...
research
06/10/2022

A Fast Spectral Solver for the Heat Equation, with Applications to Navier-Stokes

We develop a spectral method to solve the heat equation in a closed cyli...
research
05/19/2023

Stabilized finite element methods for the time-spectral convection-diffusion equation

Discretizing a solution in the Fourier domain rather than the time domai...
research
07/28/2019

Fourier spectral methods for nonlocal models

Efficient and accurate spectral solvers for nonlocal models in any spati...
research
07/20/2021

Effective maximum principles for spectral methods

Many physical problems such as Allen-Cahn flows have natural maximum pri...

Please sign up or login with your details

Forgot password? Click here to reset