CABARET scheme implementation for free shear layer modeling
In present paper we reexamine the properties of CABARET numerical scheme formulated for a weakly compressible fluid flow basing the results of free shear layer modeling. Kelvin-Helmholtz instability and successive generation of two-dimensional turbulence provide a wide field for a scheme analysis in...
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | Russian |
| Published: |
Institute of Computer Science
2017-12-01
|
| Series: | Компьютерные исследования и моделирование |
| Subjects: | |
| Online Access: | http://crm.ics.org.ru/uploads/crmissues/crm_2017_6/2017_06_03.pdf |
| _version_ | 1831690261534605312 |
|---|---|
| author | Yury Matveevich Kulikov Eduard Evgenievich Son |
| author_facet | Yury Matveevich Kulikov Eduard Evgenievich Son |
| author_sort | Yury Matveevich Kulikov |
| collection | DOAJ |
| description | In present paper we reexamine the properties of CABARET numerical scheme formulated for a weakly compressible fluid flow basing the results of free shear layer modeling. Kelvin-Helmholtz instability and successive generation of two-dimensional turbulence provide a wide field for a scheme analysis including temporal evolution of the integral energy and enstrophy curves, the vorticity patterns and energy spectra, as well as the dispersion relation for the instability increment. The most part of calculations is performed for Reynolds number $\text{Re} = 4 \times 10^5$ for square grids sequentially refined in the range of $128^2-2048^2$ nodes. An attention is paid to the problem of underresolved layers generating a spurious vortex during the vorticity layers roll-up. This phenomenon takes place only on a coarse grid with $128^2$ nodes, while the fully regularized evolution pattern of vorticity appears only when approaching $1024^2$-node grid. We also discuss the vorticity resolution properties of grids used with respect to dimensional estimates for the eddies at the borders of the inertial interval, showing that the available range of grids appears to be sufficient for a good resolution of small-scale vorticity patches. Nevertheless, we claim for the convergence achieved for the domains occupied by large-scale structures.
The generated turbulence evolution is consistent with theoretical concepts imposing the emergence of large vortices, which collect all the kinetic energy of motion, and solitary small-scale eddies. The latter resemble the coherent structures surviving in the filamentation process and almost noninteracting with other scales. The dissipative characteristics of numerical method employed are discussed in terms of kinetic energy dissipation rate calculated directly and basing theoretical laws for incompressible (via enstrophy curves) and compressible (with respect to the strain rate tensor and dilatation) fluid models. The asymptotic behavior of the kinetic energy and enstrophy cascades comply with two-dimensional turbulence laws $E(k) \propto k^{-3}, \omega^2(k) \propto k^{-1}$. Considering the instability increment as a function of dimensionless wave number shows a good agreement with other papers, however, commonly used method of instability growth rate calculation is not always accurate, so some modification is proposed. Thus, the implemented CABARET scheme possessing remarkably small numerical dissipation and good vorticity resolution is quite competitive approach compared to other high-order accuracy methods |
| first_indexed | 2024-12-20T11:06:44Z |
| format | Article |
| id | doaj.art-09d24f7925c7439cbd378bc51ccd52a4 |
| institution | Directory Open Access Journal |
| issn | 2076-7633 2077-6853 |
| language | Russian |
| last_indexed | 2024-12-20T11:06:44Z |
| publishDate | 2017-12-01 |
| publisher | Institute of Computer Science |
| record_format | Article |
| series | Компьютерные исследования и моделирование |
| spelling | doaj.art-09d24f7925c7439cbd378bc51ccd52a42022-12-21T19:42:51ZrusInstitute of Computer ScienceКомпьютерные исследования и моделирование2076-76332077-68532017-12-019688190310.20537/2076-7633-2017-9-6-881-9032631CABARET scheme implementation for free shear layer modelingYury Matveevich KulikovEduard Evgenievich SonIn present paper we reexamine the properties of CABARET numerical scheme formulated for a weakly compressible fluid flow basing the results of free shear layer modeling. Kelvin-Helmholtz instability and successive generation of two-dimensional turbulence provide a wide field for a scheme analysis including temporal evolution of the integral energy and enstrophy curves, the vorticity patterns and energy spectra, as well as the dispersion relation for the instability increment. The most part of calculations is performed for Reynolds number $\text{Re} = 4 \times 10^5$ for square grids sequentially refined in the range of $128^2-2048^2$ nodes. An attention is paid to the problem of underresolved layers generating a spurious vortex during the vorticity layers roll-up. This phenomenon takes place only on a coarse grid with $128^2$ nodes, while the fully regularized evolution pattern of vorticity appears only when approaching $1024^2$-node grid. We also discuss the vorticity resolution properties of grids used with respect to dimensional estimates for the eddies at the borders of the inertial interval, showing that the available range of grids appears to be sufficient for a good resolution of small-scale vorticity patches. Nevertheless, we claim for the convergence achieved for the domains occupied by large-scale structures. The generated turbulence evolution is consistent with theoretical concepts imposing the emergence of large vortices, which collect all the kinetic energy of motion, and solitary small-scale eddies. The latter resemble the coherent structures surviving in the filamentation process and almost noninteracting with other scales. The dissipative characteristics of numerical method employed are discussed in terms of kinetic energy dissipation rate calculated directly and basing theoretical laws for incompressible (via enstrophy curves) and compressible (with respect to the strain rate tensor and dilatation) fluid models. The asymptotic behavior of the kinetic energy and enstrophy cascades comply with two-dimensional turbulence laws $E(k) \propto k^{-3}, \omega^2(k) \propto k^{-1}$. Considering the instability increment as a function of dimensionless wave number shows a good agreement with other papers, however, commonly used method of instability growth rate calculation is not always accurate, so some modification is proposed. Thus, the implemented CABARET scheme possessing remarkably small numerical dissipation and good vorticity resolution is quite competitive approach compared to other high-order accuracy methodshttp://crm.ics.org.ru/uploads/crmissues/crm_2017_6/2017_06_03.pdfCABARET numerical schemeweakly compressible fluidKelvin – Helmholtz instabilityvorticityenstrophyinstability incrementunderresolved layersspurious vortexrollupinertial intervalcoherent structuresfilamentationdissipation ratedilatation |
| spellingShingle | Yury Matveevich Kulikov Eduard Evgenievich Son CABARET scheme implementation for free shear layer modeling Компьютерные исследования и моделирование CABARET numerical scheme weakly compressible fluid Kelvin – Helmholtz instability vorticity enstrophy instability increment underresolved layers spurious vortex rollup inertial interval coherent structures filamentation dissipation rate dilatation |
| title | CABARET scheme implementation for free shear layer modeling |
| title_full | CABARET scheme implementation for free shear layer modeling |
| title_fullStr | CABARET scheme implementation for free shear layer modeling |
| title_full_unstemmed | CABARET scheme implementation for free shear layer modeling |
| title_short | CABARET scheme implementation for free shear layer modeling |
| title_sort | cabaret scheme implementation for free shear layer modeling |
| topic | CABARET numerical scheme weakly compressible fluid Kelvin – Helmholtz instability vorticity enstrophy instability increment underresolved layers spurious vortex rollup inertial interval coherent structures filamentation dissipation rate dilatation |
| url | http://crm.ics.org.ru/uploads/crmissues/crm_2017_6/2017_06_03.pdf |
| work_keys_str_mv | AT yurymatveevichkulikov cabaretschemeimplementationforfreeshearlayermodeling AT eduardevgenievichson cabaretschemeimplementationforfreeshearlayermodeling |