kh2d-solver: A Python Library for Idealized Two-Dimensional Incompressible Kelvin-Helmholtz Instability
Fuente:
arXiv
Guardado en:
| Autores principales: | , , , , , , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2025
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866909797890129920 |
|---|---|
| author | Herho, Sandy H. S. Trilaksono, Nurjanna J. Fajary, Faiz R. Napitupulu, Gandhi Anwar, Iwan P. Khadami, Faruq Irawan, Dasapta E. |
| author_facet | Herho, Sandy H. S. Trilaksono, Nurjanna J. Fajary, Faiz R. Napitupulu, Gandhi Anwar, Iwan P. Khadami, Faruq Irawan, Dasapta E. |
| contents | We present an open-source Python library for simulating two-dimensional incompressible Kelvin-Helmholtz instabilities in stratified shear flows. The solver employs a fractional-step projection method with spectral Poisson solution via Fast Sine Transform, achieving second-order spatial accuracy. Implementation leverages NumPy, SciPy, and Numba JIT compilation for efficient computation. Four canonical test cases explore Reynolds numbers 1000--5000 and Richardson numbers 0.1--0.3: classical shear layer, double shear configuration, rotating flow, and forced turbulence. Statistical analysis using Shannon entropy and complexity indices reveals that double shear layers achieve 2.8$\times$ higher mixing rates than forced turbulence despite lower Reynolds numbers. The solver runs efficiently on standard desktop hardware, with 384$\times$192 grid simulations completing in approximately 31 minutes. Results demonstrate that mixing efficiency depends on instability generation pathways rather than intensity measures alone, challenging Richardson number-based parameterizations and suggesting refinements for subgrid-scale representation in climate models. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2509_16080 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | kh2d-solver: A Python Library for Idealized Two-Dimensional Incompressible Kelvin-Helmholtz Instability Herho, Sandy H. S. Trilaksono, Nurjanna J. Fajary, Faiz R. Napitupulu, Gandhi Anwar, Iwan P. Khadami, Faruq Irawan, Dasapta E. Fluid Dynamics Atmospheric and Oceanic Physics Computational Physics Physics Education Primary: 76E05, 76M22, 65M70, Secondary: 76D05, 65M12, 35Q30, 76F06, 65T50 We present an open-source Python library for simulating two-dimensional incompressible Kelvin-Helmholtz instabilities in stratified shear flows. The solver employs a fractional-step projection method with spectral Poisson solution via Fast Sine Transform, achieving second-order spatial accuracy. Implementation leverages NumPy, SciPy, and Numba JIT compilation for efficient computation. Four canonical test cases explore Reynolds numbers 1000--5000 and Richardson numbers 0.1--0.3: classical shear layer, double shear configuration, rotating flow, and forced turbulence. Statistical analysis using Shannon entropy and complexity indices reveals that double shear layers achieve 2.8$\times$ higher mixing rates than forced turbulence despite lower Reynolds numbers. The solver runs efficiently on standard desktop hardware, with 384$\times$192 grid simulations completing in approximately 31 minutes. Results demonstrate that mixing efficiency depends on instability generation pathways rather than intensity measures alone, challenging Richardson number-based parameterizations and suggesting refinements for subgrid-scale representation in climate models. |
| title | kh2d-solver: A Python Library for Idealized Two-Dimensional Incompressible Kelvin-Helmholtz Instability |
| topic | Fluid Dynamics Atmospheric and Oceanic Physics Computational Physics Physics Education Primary: 76E05, 76M22, 65M70, Secondary: 76D05, 65M12, 35Q30, 76F06, 65T50 |
| url | https://arxiv.org/abs/2509.16080 |