2006/10/09 by Zhaohua Yin, Z. Yin, Tao Tang +2
Engineering · Mathematics · Physics and Astronomy · #Computational Fluid Dynamics and Aerodynamics #Computational Physics (physics.comp-ph) #FOS: Physical sciences #Fluid Dynamics (physics.flu-dyn) #Fluid Dynamics and Turbulent Flows #Navier-Stokes equation solutions #Plasma Physics (physics.plasm-ph) #physics.comp-ph #physics.flu-dyn #physics.plasm-ph
paper · pdf · doi:10.48550/arxiv.physics/0610053
25 pages, 16 figures
arxiv created 2006/10/09 · openalex publication_date 2006/10/09 · arxiv updated 2009/12/01 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
To investigate the finite time singularity in three-dimensional (3D) Euler flows, the simplified model of 3D axisymmetric incompressible fluids (i.e., two-dimensional Boussinesq approximation equations) is studied numerically. The system describes a cap-like hot zone of fluid rising from the bottom, while the edges of the cap lag behind, forming eye-like vortices. The hot liquid is driven by the buoyancy and meanwhile attracted by the vortices, which leads to the singularity-forming mechanism in our simulation. In the previous 2D Boussinesq simulations, the symmetricial initial data is used. However, it is observed that the adoption of symmetry leads to coordinate singularity. Moreover, as demonstrated in this work that the locations of peak values for the vorticity and the temperature gradient becomes far apart as t approaches the predicted blow-up time. This suggests that the symmetry assumption may be unreasonable for searching solution blow-ups. One of the main contributions of this work is to propose an appropriate asymmetric initial condition, which avoids coordinate singularity and also makes the blow-up to occur much earlier than that given by the previously simulations. The shorter simulation time suppresses the development of the round-off error. On the numerical side, the pseudo-spectral method with filtering technique is adopted. The resolutions adopted in this study vary from 10242, 20482, 40962 to 61442. With our proposed asymmetric initial condition, it is shown that the 40962 and 61442 runs yield convergent results when t is fairly close to the predicted blow-up time. Moreover, as expected the locations of peak values for the vorticity and the temperature gradient are very close to each other as t approaches the predicted blow-up time.