ترغب بنشر مسار تعليمي؟ اضغط هنا

Numerical investigations on the finite time singularity in two-dimensional Boussinesq equations

86   0   0.0 ( 0 )
 نشر من قبل Zhaohua Yin
 تاريخ النشر 2006
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

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 $1024^2$, $2048^2$, $4096^2$ to $6144^2$. With our proposed asymmetric initial condition, it is shown that the $4096^2$ and $6144^2$ 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.



قيم البحث

اقرأ أيضاً

Most fluid flow problems that are vital in engineering applications involve at least one of the following features: turbulence, shocks, and/or material interfaces. While seemingly different phenomena, these flows all share continuous generation of hi gh wavenumber modes, which we term the $k_infty$ irregularity. In this work, an inviscid regularization technique called observable regularization is proposed for the simulation of two-phase compressible flows. The proposed approach regularizes the equations at the level of the partial differential equation and as a result, any numerical method can be used to solve the system of equations. The regularization is accomplished by introducing an observability limit that represents the length scale below which one cannot properly model or continue to resolve flow structures. An observable volume fraction equation is derived for capturing the material interface, which satisfies the pressure equilibrium at the interface. The efficacy of the observable regularization method is demonstrated using several test cases, including a one-dimensional material interface tracking, one-dimensional shock-tube and shock-bubble problems, and two-dimensional simulations of a shock interacting with a cylindrical bubble. The results show favorable agreement, both qualitatively and quantitatively, with available exact solutions or numerical and experimental data from the literature. The computational saving by using the current method is estimated to be about one order of magnitude in two-dimensional computations and significantly higher in three-dimensional computations. Lastly, the effect of the observability limit and best practices to choose its value are discussed.
351 - Yifei Guan , James Riley , 2019
The study focuses on the 3D electro-hydrodynamic (EHD) instability for flow between to parallel electrodes with unipolar charge injection with and without cross-flow. Lattice Boltzmann Method (LBM) with two-relaxation time (TRT) model is used to stud y flow pattern. In the absence of cross-flow, the base-state solution is hydrostatic, and the electric field is one-dimensional. With strong charge injection and high electrical Rayleigh number, the system exhibits electro-convective vortices. Disturbed by different perturbation patterns, such as rolling pattern, square pattern, and hexagon pattern, the flow patterns develop according to the most unstable modes. The growth rate and the unstable modes are examined using dynamic mode decomposition (DMD) of the transient numerical solutions. The interactions between the applied Couette and Poiseuille cross-flows and electroconvective vortices lead to the flow patterns change. When the cross-flow velocity is greater than a threshold value, the spanwise structures are suppressed; however, the cross-flow does not affect the streamwise patterns. The dynamics of the transition is analyzed by DMD. Hysteresis in the 3D to 2D transition is characterized by the non-dimensional parameter Y, a ratio of the coulombic force to viscous term in the momentum equation. The change from 3D to 2D structures enhances the convection marked by a significant increase in the electric Nusselt number.
We numerically analyze Non-Oberbeck-Boussinesq (NOB) effects in two-dimensional Rayleigh-Benard flow in glycerol, which shows a dramatic change in the viscosity with temperature. The results are presented both as functions of the Rayleigh number (Ra) up to $10^8$ (for fixed temperature difference between the top and bottom plates) and as functions of non-Oberbeck-Boussinesqness or NOBness ($Delta$) up to 50 K (for fixed Ra). For this large NOBness the center temperature $T_c$ is more than 5 K larger than the arithmetic mean temperature $T_m$ between top and bottom plate and only weakly depends on Ra. To physically account for the NOB deviations of the Nusselt numbers from its Oberbeck-Boussinesq values, we apply the decomposition of $Nu_{NOB}/Nu_{OB}$ into the product of two effects, namely first the change in the sum of the top and bottom thermal BL thicknesses, and second the shift of the center temperature $T_c$ as compared to $T_m$. While for water the origin of the $Nu$ deviation is totally dominated by the second effect (cf. Ahlers et al., J. Fluid Mech. 569, pp. 409 (2006)) for glycerol the first effect is dominating, in spite of the large increase of $T_c$ as compared to $T_m$.
In this numerical study, an original approach to simulate non-isothermal viscoelastic fluid flows at high Weissenberg numbers is presented. Stable computations over a wide range of Weissenberg numbers are assured by using the root conformation approa ch in a finite volume framework on general unstructured meshes. The numerical stabilization framework is extended to consider thermo-rheological properties in Oldroyd-B type viscoelastic fluids. The temperature dependence of the viscoelastic fluid is modeled with the time-temperature superposition principle. Both Arrhenius and WLF shift factors can be chosen, depending on the flow characteristics. The internal energy balance takes into account both energy and entropy elasticity. Partitioning is achieved by a constant split factor. An analytical solution of the balance equations in planar channel flow is derived to verify the results of the main field variables and to estimate the numerical error. The more complex entry flow of a polyisobutylene-based polymer solution in an axisymmetric 4:1 contraction is studied and compared to experimental data from the literature. We demonstrate the stability of the method in the experimentally relevant range of high Weissenberg numbers. The results at different imposed wall temperatures, as well as Weissenberg numbers, are found to be in good agreement with experimental data. Furthermore, the division between energy and entropy elasticity is investigated in detail with regard to the experimental setup.
Using complementary numerical approaches at high resolution, we study the late-time behaviour of an inviscid, incompressible two-dimensional flow on the surface of a sphere. Starting from a random initial vorticity field comprised of a small set of i ntermediate wavenumber spherical harmonics, we find that -- contrary to the predictions of equilibrium statistical mechanics -- the flow does not evolve into a large-scale steady state. Instead, significant unsteadiness persists, characterised by a population of persistent small-scale vortices interacting with a large-scale oscillating quadrupolar vorticity field. Moreover, the vorticity develops a stepped, staircase distribution, consisting of nearly homogeneous regions separated by sharp gradients. The persistence of unsteadiness is explained by a simple point vortex model characterising the interactions between the four main vortices which emerge.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا