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

Numerical simulation of compressible flows using two-phase observable Euler equations

72   0   0.0 ( 0 )
 نشر من قبل Kamran Mohseni
 تاريخ النشر 2018
  مجال البحث فيزياء
والبحث باللغة English




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

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 high 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.



قيم البحث

اقرأ أيضاً

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.
A two-fluid Discrete Boltzmann Model(DBM) for compressible flows based on Ellipsoidal Statistical Bhatnagar-Gross-Krook(ES-BGK) is presented. The model has flexible Prandtl number or specific heat ratio. Mathematically, the model is composed of two c oupled Discrete Boltzmann Equations(DBE). Each DBE describes one component of the fluid. Physically, the model is equivalent to a macroscopic fluid model based on Navier-Stokes(NS) equations, and supplemented by a coarse-grained model for thermodynamic non-equilibrium behaviors. To obtain a flexible Prandtl number, a coefficient is introduced in the ellipsoidal statistical distribution function to control the viscosity. To obtain a flexible specific heat ratio, a parameter is introduced in the energy kinetic moments to control the extra degree of freedom. For binary mixture, the correspondence between the macroscopic fluid model and the DBM may be several-to-one. Five typical benchmark tests are used to verify and validate the model. Some interesting non-equilibrium results, which are not available in the NS model or the single-fluid DBM, are presented.
We present a high-order implicit large-eddy simulation (ILES) approach for simulating transitional turbulent flows. The approach consists of an Interior Embedded Discontinuous Galerkin (IEDG) method for the discretization of the compressible Navier-S tokes equations and a parallel preconditioned Newton-GMRES solver for the resulting nonlinear system of equations. The IEDG method arises from the marriage of the Embedded Discontinuous Galerkin (EDG) method and the Hybridizable Discontinuous Galerkin (HDG) method. As such, the IEDG method inherits the advantages of both the EDG method and the HDG method to make itself well-suited for turbulence simulations. We propose a minimal residual Newton algorithm for solving the nonlinear system arising from the IEDG discretization of the Navier-Stokes equations. The preconditioned GMRES algorithm is based on a restricted additive Schwarz (RAS) preconditioner in conjunction with a block incomplete LU factorization at the subdomain level. The proposed approach is applied to the ILES of transitional turbulent flows over a NACA 65-(18)10 compressor cascade at Reynolds number 250,000 in both design and off-design conditions. The high-order ILES results show good agreement with a subgrid-scale LES model discretized with a second-order finite volume code while using significantly less degrees of freedom. This work shows that high-order accuracy is key for predicting transitional turbulent flows without a SGS model.
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.
Sand traps are used to measure Aeolian flux. Since they modify the surrounding wind velocity field their gauging represents an important challenge. We use numerical simulations under the assumption of homogeneous turbulence based on FLUENT to systema tically study the flow field and trapping efficiency of one of the most common devices based on a hollow cylinder with two slits. In particular, we investigate the dependence on the wind speed, the Stokes number, the permeability of the membrane on the slit and the saltation height.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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