Evolution of the Diffusion-Induced Flow over a Disk Submerged in a Stratified Viscous Fluid
Autor: | P. V. Matyushin |
---|---|
Rok vydání: | 2019 |
Předmět: |
Convection
Physics Buoyancy 010102 general mathematics Mechanics engineering.material Viscous liquid 01 natural sciences Instability 010305 fluids & plasmas Vortex Physics::Fluid Dynamics Computational Mathematics Neutral buoyancy Modeling and Simulation 0103 physical sciences engineering 0101 mathematics Boussinesq approximation (water waves) Convection cell |
Zdroj: | Mathematical Models and Computer Simulations. 11:479-487 |
ISSN: | 2070-0490 2070-0482 |
DOI: | 10.1134/s2070048219030141 |
Popis: | The paper presents the results of the mathematically simulated evolution of a 3D diffusion-induced flow over a disk (with diameter d and thickness H = 0.76 d) immersed in a linearly density-stratified incompressible viscous fluid (described by a system of the Navier–Stokes equations in the Boussinesq approximation). The disk rests at the level of the neutral buoyancy (coinciding with its symmetry axis z) and disturbs the homogeneity of the background diffusion flux in the fluid forming a complex system of slow currents (internal gravitational waves). Over time, two thin horizontal convection cells are formed at the upper and lower parts of the disk stretching parallel to the z axis and adjacent to the base cell with thickness d/2. This work is the first to analyze in detail the fundamental mechanism for the formation of each new half-wave near the vertical axis x (passing through the center of the disk) during half the buoyancy period of the fluid Tb. This mechanism is based on gravitational instability. The emergence of this instability is first detected at 0.473Tb at a height of 3.9 d above the center of the disk. The same mechanism is also implemented over the place where the body moves in the horizontal direction. The 3D vortex structure of the flow is visualized by constructing the isosurfaces of the imaginary part of the complex-conjugate eigenvalues of the velocity gradient tensor. In mathematical simulation, we employed a numerical method SMIF that has proved itself over three decades with an explicit hybrid finite-difference scheme for the approximation of the convective terms of the equations (second-order approximation, monotonicity). |
Databáze: | OpenAIRE |
Externí odkaz: |