Search results

The purpose of the paper is to disclose the effect of the relative position (d) between the impeller and non-vane cavity on the hydraulic performance and unsteady characteristics of vortex pump.

Three groups of vortex pump models with different impeller installation positions were analyzed and studied by combining experimental and CFD (Computational Fluid Dynamics) numerical calculations.

The steady numerical results show that as the width (d) of the impeller moves into the non-vane cavity increases, the proportion of circulation flow in the non-vane cavity is reduced and both the pump head and efficiency are on the rise. The unsteady numerical results and the Enstrophy analysis show that the dynamic and static interference between the circulation flow and the volute tongue is the main reason for the pressure pulsation with a frequency of 2fn in the vortex pump. With the increase of the d value, the dynamic and static interference between the circulation flow and the volute tongue is enhanced. The pulsation amplitude at the volute tongue of the d = 16.5 mm model increases about six times compared with the d = 0 mm model; the distribution of the vortex core in the non-vane cavity is closely related to the position of the impeller, and the peak of the Enstrophy of the circulation flow vortex belt always appears at the top of the impeller.

The research results provide a theoretical foundation for the optimization and improvement of the vortex pump.

The purpose of this study is to investigate the accuracy and applicability of the Flowfield Dependent Variation (FDV) method for large-eddy simulations (LES) of decaying isotropic turbulence.

In an earlier paper, the FDV method was successfully demonstrated for simulations of laminar flows with speeds varying from low subsonic to high supersonic Mach numbers. In the current study, the FDV method, implemented in a finite element framework, is used to perform LESs of decaying isotropic turbulence. The FDV method is fundamentally derived from the Lax–Wendroff Scheme (LWS) by replacing the explicit time derivatives in LWS with a weighted combination of explicit and implicit time derivatives. The increased implicitness and the inherent numerical dissipation of FDV contribute to the scheme’s numerical stability and monotonicity. Understanding the role of numerical dissipation that is inherent to the FDV method is essential for the maturation of FDV into a robust scheme for LES of turbulent flows. Accordingly, three types of LES of decaying isotropic turbulence were performed. The first two types of LES utilized explicit subgrid scale (SGS) models, namely, the constant-coefficient Smagorinsky and dynamic Smagorinsky models. In the third, no explicit SGS model was employed; instead, the numerical dissipation inherent to FDV was used to emulate the role played by explicit SGS models. Such an approach is commonly known as Implicit LES (ILES). A new formulation was also developed for quantifying the FDV numerical viscosity that principally arises from the convective terms of the filtered Navier–Stokes equations.

The temporal variation of the turbulent kinetic energy and enstrophy and the energy spectra are presented and analyzed. At all grid resolutions, the temporal profiles of kinetic energy showed good agreement with t^(−1.43) theoretical scaling in the fully developed turbulent flow regime, where t represents time. The energy spectra also showed reasonable agreement with the Kolmogorov’s k^(−5/3) power law in the inertial subrange, with the spectra moving closer to the Kolmogorov scaling at higher-grid resolutions. The intrinsic numerical viscosity and the dissipation rate of the FDV scheme are quantified, both in physical and spectral spaces, and compared with those of the two SGS LES runs. Furthermore, at a finite number of flow realizations, the numerical viscosities of FDV and of the Streamline Upwind/Petrov–Galerkin (SUPG) finite element method are compared. In the initial stages of turbulence development, all three LES cases have similar viscosities. But, once the turbulence is fully developed, implicit LES is less dissipative compared to the two SGS LES runs. It was also observed that the SUPG method is significantly more dissipative than the three LES approaches.

Just as any computational method, the limitations are based on the available computational resources.

Solving problems involving turbulent flows is by far the biggest challenge facing engineers and scientists in the twenty-first century, this is the road that the authors have embarked upon in this paper and the road ahead of is very long.

Understanding turbulence is a very lofty goal and a challenging one as well; however, if the authors succeed, the rewards are limitless.

The derivation of an explicit expression for the numerical viscosity tensor of FDV is an important contribution of this study, and is a crucial step forward in elucidating the fundamental properties of the FDV method. The comparison of viscosities for the three LES cases and the SUPG method has important implications for the application of ILES approach for turbulent flow simulations.

The purpose of this paper is to suggest a new approach to the numerical simulation of shallow‐water flows both in plane domains and on the sphere.

The approach involves the technique of splitting of the model operator by geometric coordinates and by physical processes. Specially chosen temporal and spatial approximations result in one‐dimensional finite difference schemes that conserve the mass and the total energy. Therefore, the mass and the total energy of the whole two‐dimensional split scheme are kept constant too.

Explicit expressions for the schemes of arbitrary approximation orders in space are given. The schemes are shown to be mass‐ and energy‐conserving, and hence absolutely stable because the square root of the total energy is the norm of the solution. The schemes of the first four approximation orders are then tested by simulating nonlinear solitary waves generated by a model topography. In the analysis, the primary attention is given to the study of the time‐space structure of the numerical solutions.

The approach can be used for the numerical simulation of shallow‐water flows in domains of both Cartesian and spherical geometries, providing the solution adequate from the physical and mathematical standpoints in the sense of keeping its mass and total energy constant even when fully discrete shallow‐water models are applied.

This study is concerned with the direct numerical simulation (DNS) of a turbulent channel flow by an improved vortex in cell (VIC) method. The paper aims to discuss these issues.

First, two improvements for VIC method are proposed to heighten the numerical accuracy and efficiency. A discretization method employing a staggered grid is presented to ensure the consistency among the discretized equations as well as to prevent the numerical oscillation of the solution. A correction method for vorticity is also proposed to compute the vorticity field satisfying the solenoidal condition. Second, the DNS for a turbulent channel flow is conducted by the improved VIC method. The Reynolds number based on the friction velocity and the channel half width is 180.

It is highlighted that the simulated turbulence statistics, such as the mean velocity, the Reynolds shear stress and the budget of the mean enstrophy, agree well with the existing DNS results. It is also shown that the organized flow structures in the near-wall region, such as the streaks and the streamwise vortices, are favourably captured. These demonstrate the high applicability of the improved VIC method to the DNS for wall turbulent flows.

This study enables the VIC method to perform the DNS for wall turbulent flows.

Numerical simulations are performed for studying the vorticity dynamics of a dipole colliding with the wall in a bounded flow and the wake structure and separated flow properties past a circular cylinder at the values of Reynolds numbers.

The near wake statistics of separated fluid flows are investigated by using the lattice Boltzmann method (LBM) in a two-dimensional framework. A multi-block technique is applied to accurately resolve the flow characteristics by the grid refinement near the wall and preserve the stability of the numerical solution at relatively high Reynolds numbers.

The results show that the rolling-up of the boundary layer occurs due to the shear-layer instabilities near the surface which causes a boundary layer detachment from the wall and consequently leads to the formation of small-scale vortices. These shear-layer vortices shed at higher frequencies than the large-scale Strouhal vortices which result in small-scale high-frequency fluctuations in the velocity field in the very near wake. The present study also demonstrates that the efficiency of the multi-block LBM used for predicting the statistical features of flow problems is comparable with the solvers based on the Navier-Stokes equations.

Studying the separated flow characteristics in aerospace applications.

Applying a multi-block lattice Boltzmann method (LBM) for simulation of separated fluid flows at high-Reynolds numbers. Studying of the near wake statistics of unsteady separated fluid flows using the multi-block LBM. Comparison of flow characteristics obtained based on the LBM with those of reported based on the Navier-Stokes equations.

Vorticity formulations for the incompressible Navier‐Stokes equations have certain advantages over primitive‐variable formulations including the fact that the number of equations to be solved is reduced through the elimination of the pressure variable, identical satisfaction of the incompressibility constraint and the continuity equation, and an implicitly higher‐order approximation of the velocity components. For the most part, vorticity methods have been used to solve exterior isothermal problems. In this research, a vorticity formulation is used to study the natural convection flows in differentially‐heated enclosures. The numerical algorithm is divided into three steps: two kinematic steps and one kinetic step. The kinematics are governed by the generalized Helmholtz decomposition (GHD) which is solved using a boundary element method (BEM) whereas the kinetics are governed by the vorticity equation which is solved using a finite element method (FEM). In the first kinematic step, vortex sheet strengths are determined from a novel Galerkin implementation of the GHD. These vortex sheet strengths are used to determine Neumann boundary conditions for the vorticity equation. (The thermal boundary conditions are already known.) In the second kinematic step, the interior velocity field is determined using the regular (non‐Galerkin) form of the GHD. This step, in a sense, linearizes the convective acceleration terms in both the vorticity and energy equations. In the third kinetic step, the coupled vorticity and energy equations are solved using a Galerkin FEM to determine the updated values of the vorticity and thermal fields. Two benchmark problems are considered to show the robustness and versatility of this formulation including natural convection in an 8×1 differentially‐heated enclosure at a near critical Rayleigh number.

The purpose of this paper is to review the available reduced order modeling approaches in the literature for predicting the flow and specially temperature fields inside data centers in terms of the involved design parameters.

This paper begins with a motivation for flow/thermal modeling needs for designing an energy‐efficient thermal management system in data centers. Recent studies on air velocity and temperature field simulations in data centers through computational fluid dynamics/heat transfer (CFD/HT) are reviewed. Meta‐modeling and reduced order modeling are tools to generate accurate and rapid surrogate models for a complex system. These tools, with a focus on low‐dimensional models of turbulent flows are reviewed. Reduced order modeling techniques based on turbulent coherent structures identification, in particular the proper orthogonal decomposition (POD) are explained and reviewed in more details. Then, the available approaches for rapid thermal modeling of data centers are reviewed. Finally, recent studies on generating POD‐based reduced order thermal models of data centers are reviewed and representative results are presented and compared for a case study.

It is concluded that low‐dimensional models are needed in order to predict the multi‐parameter dependent thermal behavior of data centers accurately and rapidly for design and control purposes. POD‐based techniques have shown great approximation for multi‐parameter thermal modeling of data centers. It is believed that wavelet‐based techniques due to the their ability to separate between coherent and incoherent structures – something that POD cannot do – can be considered as new promising tools for reduced order thermal modeling of complex electronic systems such as data centers

The paper reviews different numerical methods and provides the reader with some insight for reduced order thermal modeling of complex convective systems such as data centers.

A computational procedure based on a hybrid Lagrangian‐Eulerian discrete‐vortical element formulation and conformal transformation schemes are employed in this study to simulate the interaction of an air jet with swirling air flow inside a two‐dimensional cylinder. Such an investigation is of importance to many flow‐related industrial and environmental problems, such as mixing, cooling, combustion and dispersion of air‐borne or water‐borne contaminants because of the role of vortices in the global transport of matter and heat. The basis for the simulation is discussed and numerical results compared with theoretical results for the velocity field and streamfunction obtained by the method of images. The swirling air motion and the features of a real jet are well simulated and numerical results are validated by predictions of theory to within 20 per cent. To illustrate the merging and interaction processes of vortices and the formation of large eddies, velocity vectors, particle trajectories and streamline contours are presented.