Advances in Mathematical Physics
Volume 2015 (2015), Article ID 657620, 9 pages
http://dx.doi.org/10.1155/2015/657620

Research Article

Anisotropic Characteristics of Turbulence Dissipation in Swirling Flow: A Direct Numerical Simulation Study

Xingtuan Yang,¹ Nan Gui,1,2 Gongnan Xie,³ Jie Yan,⁴ Jiyuan Tu,1,5 and Shengyao Jiang¹

¹Institute of Nuclear and New Energy Technology and Key Laboratory of Advanced Reactor Engineering and Safety, Ministry of Education, Beijing 100084, China
²College of Mechanical and Transportation Engineering, China University of Petroleum, Beijing 102249, China
³School of Mechanical Engineering, Northwestern Polytechnical University, Xi’an 710072, China
⁴China Academy of Space Technology, Beijing 100094, China
⁵School of Aerospace, Mechanical & Manufacturing Engineering, RMIT University, Melbourne, VIC 3083, Australia

Received 7 August 2014; Revised 13 December 2014; Accepted 15 December 2014

Academic Editor: Alina Adriana Minea

Copyright © 2015 Xingtuan Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This study investigates the anisotropic characteristics of turbulent energy dissipation rate in a rotating jet flow via direct numerical simulation. The turbulent energy dissipation tensor, including its eigenvalues in the swirling flows with different rotating velocities, is analyzed to investigate the anisotropic characteristics of turbulence and dissipation. In addition, the probability density function of the eigenvalues of turbulence dissipation tensor is presented. The isotropic subrange of PDF always exists in swirling flows relevant to small-scale vortex structure. Thus, with remarkable large-scale vortex breakdown, the isotropic subrange of PDF is reduced in strongly swirling flows, and anisotropic energy dissipation is proven to exist in the core region of the vortex breakdown. More specifically, strong anisotropic turbulence dissipation occurs concentratively in the vortex breakdown region, whereas nearly isotropic turbulence dissipation occurs dispersively in the peripheral region of the strong swirling flows.

1. Introduction

Vortex breakdown is an intriguing and important phenomenon that occurs in a variety of natural and technological swirling flows, for example, swirling combustor, cyclone, bathtub vortex, hurricanes and tornadoes, and spiral galaxies. In general, the appearance of swirl is caused by the impartment of rotating motion upon the jet which makes the flow more complicated and has been widely studied in the literature [1–10]. The scope of this study does not include summarizing existing vast literature of swirling flow investigation; thus, only a few investigations on this topic were covered.

For example, Chen and Sun [1] addressed the nonlinear 3D instability of a specific type of viscous swirling flow, the Ossen vortex, by using direct numerical simulation at Re = 5000. They considered the global optimal perturbation as the initial perturbation and characterized different flow regimes in axisymmetric cases. Wang and Chen [2] studied vortex breakdown by solving 3D unsteady Navier-Stokes equations for swirling pipe flows, including the flow structures in the bubble domain and the tails behind the breakdown vortex. Moreover, Shtern et al. [4] studied symmetry breaking in a meridional steady motion of viscous incompressible fluids using the laminar axisymmetric “vortex dynamo.” They demonstrated the feasibility of a supercritical pitchfork bifurcation from an initially nonswirling flow to a steady swirling regime, as Re exceeds a critical value. In addition, with regard to model study, di Pierro and Abid [5] investigated modeled inviscid swirling flows, addressing the weak variation of axial and azimuthal velocities. They checked the asymptotic results using numerically computed growth rates of linearized Euler equations for a family of variable-density Batchelor-like vortices as base flows and so forth.

Turbulent flows are well known to have scientific and practical importance. However, the full spectrums of the length and time scales of turbulent flows are impossible to solve via computer simulation. Thus, some types of modeling for Reynolds stresses are needed to simulate high Reynolds turbulence. For example, the turbulent kinetic energy and dissipation are involved in the standard form of the two-equation Reynolds stress turbulence model based on the Boussinesq-type approximation [11]:

where is the turbulent kinetic energy, is the time-averaged velocity, and is the isotropic turbulent viscosity. However, this turbulence model is fairly ineffective in simulating anisotropic turbulence due to the theoretical deficiency for anisotropic nature of turbulent motions, although it is widely applied in engineering. Thus, numerous studies have been carried out to further understand the physical aspects of anisotropic turbulent flows and vortex dynamics [12–14].

In conclusion, swirling flow is present in anisotropic turbulent flows and is a striking and intriguing case of the generation and breakdown of strong vortices. However, the characteristics of the anisotropic nature of turbulent motion and energy dissipation are not clear, especially in strong swirling turbulent jet flows. Thus, the present study carried out a numerical study on the characteristics of energy dissipation tensor in a rotating jet flow. Three swirling numbers are used, corresponding to weak, intermediate, and strong swirling levels. The relations of the components of energy dissipation corresponding to normal and shearing turbulent fluctuations are explored, including their probability density function. The specific locations in the swirling flows, corresponding to the regions of extremely anisotropic and nearly isotropic turbulent dissipation, are also presented. The locations indicate the correlation of anisotropic turbulence dissipation to the large-scale structure of vortex breakdown and the correlation of isotropic turbulence dissipation to small-scale vortices.

2. Numerical Description

2.1. Governing Equations

For incompressible Newtonian fluids, the governing equations can be expressed in dimensionless form as follows:

where and are the velocity and pressure, respectively. Re = is the Reynolds number, in which is the inflow velocity, is the jet diameter at the inlet, and is the kinematic viscosity.

Equations (2) and (3) deal with a three-dimensional time-dependent flow problem without body force. To solve them, the finite difference method is applied, where the convection term is discretized by upwind compact schemes [15], and the space derivatives and pressure-gradient terms are discretized by fourth-order compact difference schemes [16], respectively. The third-order explicit schemes are used to deal with the boundary points and to maintain the global fourth-order spatial accuracy. The time stepping process is integrated by using the fourth-order Runge-Kutta schemes [17]. The pressure-Poisson equation is solved to obtain pressure via the fourth-order finite difference method [18]. The validation of the codes has been done in a recently published work on swirling flows [19].

2.2. Boundary Conditions

As shown in Figure 1(a), the flow configuration contains a rectangular flow domain of , where is the diameter of the jet inlet and the jets are injected from the inlet with a mean velocity . The entire flow domain is discretized by Cartesian mesh grids, and the mesh scale is not larger than . The Kolmogorov length scale is estimated in the same order of the finest mesh scale, namely, [19, 20]. According to the suggestion by Moin and Mahesh [21], the mesh scale in the same order of the Kolmogorov scale is fine enough to capture the smallest scale of turbulence. Thus, the grid testing procedure is omitted here. The time step is , and 20000 time steps are simulated for each case. The nonreflecting boundary condition is utilized for the outlet condition [22], and the sidewalls are set as nonslipping wall boundaries.

(a)

(b)

Figure 1: Sketch of simulation setup: typical vortex streets (a) and inflow velocity profiles (b).

For axial inflow velocity, a hyperbolic tangent profile is used within the range of :

where is the center position of the jet. The reference system is centered at the inlet with . is the inflow momentum thickness (Table 1) and the ratio of the jet width to the inflow momentum thickness is . For azimuthal inflow velocity, a polynomial expression is used; that is,

where the coefficients are listed in Table 2. The combined profiles of axial and azimuthal velocities are shown in Figure 1(b). In addition, no initial turbulence is introduced to show the intrinsic full evolution of coherent vortex structures and interactions.

Table 1: Typical real values of variables used in the numerical simulation.

Scales of the flow domain (mm)
20 * 10 * 10

Scale of jet diameter (mm)
1.0

Grids number
512 * 256 * 256

Scale of resolution (m)
39

Mean inlet axial velocity (m/s)
71.71

Fluid density (kg/m³)
1.29

Fluid viscosity (Pa⋅s)
1.85 * 10⁻⁵

Inflow momentum thickness (/d)
1/20

Reynolds number Re
5000

Swirl number S
Case 1: 0.28;

Case 2: 0.45;

Case 3: 0.59;

Time step
0.001

Number of simulation step
20000

Table 2: Coefficients of the expressions for inflow velocity .

0.28
0
0.72
0.24
16.72
−1.40
−62.38

0.45
0
1.14
0.38
26.56
−2.22
−99.07

0.59
0
1.20
0.50
34.92
−2.92
−130.26

In addition, swirl number is defined as the ratio of the axial flux of angular momentum to the axial momentum; that is,

Based on (4) and (5), three swirl numbers are used (Table 2). From the engineering viewpoint, is considered as low swirl jets, whereas is considered as strong swirling flows and is intermediate.

2.3. Reynolds Stress Transport Equation

The Reynolds-averaged Navier-Stokes equation [23] is usually used to solve the turbulent flows in industrial scales:

where is the time-averaged operator. is the Reynolds stress tensor, which can be solved through the Reynolds stress transport equation and is given by

where is the turbulent energy dissipation tensor. The Reynolds stress tensor

is relevant to the dissipation of turbulent kinetic energy and is a symmetric positive definite tensor. Thus, it corresponds to a diagonalizable matrix, and the anisotropic nature of the turbulent kinetic energy dissipation is indicated from the eigenvalues of the diagonalizable matrix. The anisotropic tensor has three unequal eigenvalues (i.e., at least two of them are unequal). Thus, the discrepancy between the eigenvalues, or equivalently, their relation and distribution, can indicate the anisotropic nature of turbulent kinetic energy dissipation in swirling flows.

3. Results and Discussions

3.1. Coherent Vortex Structures

The typical vortex structures are shown in Figure 1(a) when and . For , the Kelvin-Helmholtz instability dominates the vortex evolution due to the existence of a shear layer between the jets and ambient fluids. In contrast, the strong rotating effect dominates the evolution of the vortex structures for . A cone-type vortex breakdown is established, and the K-H instability causes the breakdown of the “cone” feature from its terminal into the turbulent vortex streets.

Figure 2 shows that the typical cross-sectional visualizations of the vortex structures are visualized for (referring to the jet diameter for the region width). The early evolution of the vortex for (Figure 2(a)) has a full ring-type structure, with a positive stream-wise velocity () in the inner side of the vortex ring (flood by red color in Figure 2(a) for 0.2 < < 0.8) and a negative stream-wise velocity () in the outer side of the ring (flood by blue color in Figure 2(a) for ). Central recirculation then occurs in the center of the jet, where the velocities are weakly negative. The ring structure of vortex for (Figure 2(b)) evolves into more complicated structures. The diameter of the vortex ring tube becomes larger and starts to break down, and braid vortices are then formed. The complex vortex structures indicate the existence of turbulent vortex motions and dissipation of turbulent energy.

Figure 2: Cross-sectional visualization of isovortex cores (, flood by axial velocity ) at at (a) and (b) in the central subregion.

Based on the observations of Figure 2, the intrinsic characteristics of the vortex can be reflected by turbulent energy dissipations. However, these characteristics need to be explored to show the anisotropic/isotropic turbulent energy dissipation tensor and its correlation to vortex structures, among others.

3.2. Dissipation Tensor

3.2.1. Relation of Components of

As previously mentioned, dissipates turbulent kinetic energy . is a symmetric positive definite tensor, which can be diagonalized. In general, three eigenvalues and three eigenvectors can be obtained in the diagonalization process. The eigenvectors are orthogonal and denote the three characteristic directions. The corresponding eigenvalues indicate the magnitudes of dissipation in each direction. Thus, an anisotropic tensor should have three (at least two) different eigenvalues, which correspond to different levels of dissipation in the three characteristic directions. In this way, the characteristics of anisotropic turbulent motion can be reflected by the characteristics of the eigenvalues and eigenvectors of dissipation tensor . was used because Re = 5000 is a constant.

First, Figure 3 illustrates the relation of the components of the tensor . Each point in Figure 3 corresponds to the energy dissipation characteristics at each point location in the flow, with part in the -axis and part in the -axis, because the dissipation tensor has different values for different locations in the flow domain. The magnitude of fluctuation is almost the same as , despite the swirling levels. Thus the energy dissipation of is almost the same as that of . It is possibly caused by the transport of turbulent kinetic energy between and , which results in dynamic balances of turbulent kinetic energy and dissipation between these types of turbulent fluctuation.

Figure 3: Relation of the components of the dissipation tensor ( in the -axis and in the -axis, with and , normalized by the root-mean-square values and of them, resp.).

3.2.2. Distribution of the Eigenvalues

Second, tensor was diagonalized, where , , and are the three eigenvalues of the tensor, which were sorted as for simplicity. The distribution and relation of are illustrated in Figure 4 using and for the case study. Figure 4(a) () shows that the distribution of has a widespread area below the line of , indicating that it has in most regions. Similar results were obtained in Figure 4(c) for , and extreme discrepancies between and were observed, even as large as more than two or three orders. This observation indicates the extreme anisotropic characteristics of turbulent kinetic energy dissipation when the swirling level is sufficiently strong. Moreover, for (Figures 4(b) and 4(d)), although extreme events seem to be prohibited due to the possible existence of superior limit line, the distribution area is open and enlarged as increases. In conclusion, Figure 4 shows the validation of the anisotropic characteristics of turbulent energy dissipation tensor, especially for strong swirling flows where extremely anisotropic events may occur.

(a)

(b)

(c)

(d)

Figure 4: Relation of eigenvalues. (a) , ; (b) , ; (c) , ; (d) , .

3.2.3. Probability Density Functions

Based on the combined distributions (e.g., and ), the probability density function (PDF) of one eigenvalue can be obtained by integration of the other eigenvalue throughout the flow domain; that is, . In general, PDF is defined as the concentration or density of located within the range of ; that is,

where is the total number of points within the flow.

Figure 5 shows the PDFs for different swirling flows. For all the cases, that is, (Figure 5(a)), (Figure 5(b)), and (Figure 5(c)), the following similar trends have been observed: In most ranges, > > for the same eigenvalue of and for the same value of PDF , which means is true in most ranges; that is, the anisotropic turbulent dissipation occurs in most ranges. A subrange of , which has , always exists, except that for most strong swirling flows . This phenomenon indicates the existence of subrange isotropic turbulent dissipation when the swirling level is not strong. On the other hand, turbulent energy dissipation in the major characteristic direction corresponding to the largest eigenvalues is always larger than the other characteristic directions for strong swirling flows (), which indicates highly anisotropic turbulent energy dissipation characteristics in the main flow region.

(a)

(b)

(c)

(d)

Figure 5: Probability density function for (a), (b), (c), and the linear subrange of them (d).

Moreover, Figures 5(a) to 5(c) show that the isotropic subrange can be characterized by a regression line, designated as , , and for , , and , respectively. The data in these subranges are illustrated in Figure 5(d) for clarity. Figure 5(d) shows a perfect linearity between and for the isotropic subrange, which obtained the following equation using linear regression:

The regression coefficients are listed in Table 3, which shows that the power-law exits in the isotropic subrange of turbulent energy dissipation and that the power exponent is approximately −1 to −0.8. Moreover, the power-law occurs mainly for low turbulent energy dissipations, and the width of strong swirling flows is reduced. In other words, weak swirling flow has wider subranges of isotropic turbulent energy dissipation than the strong swirling flow and vice versa. This phenomenon indicates the significant role of swirling motion in augmenting the anisotropic characteristics of turbulence. To speak specifically, it is possible to assume that turbulent kinetic energy dissipation is more intensive and more nonlinear in strongly swirling flows than in weakly swirling flows. Therefore, the subranges of linear relations between probability distribution functions and eigenvalues of dissipation rates are reduced in strongly swirling flows compared to those in weakly swirling flows.

Table 3: Coefficients of the regression lines of , , and .

= 0.28
−0.823
−1.735
0.018

= 0.45
−0.976
−1.257
0.055

= 0.59
−0.883
−1.245
0.057

3.2.4. Correlation to Vortex Structures

Previously mentioned results are based on statistical analysis; thus visualizing the locations of extremely anisotropic and nearly isotropic turbulence dissipation should be helpful to further understand anisotropic turbulent swirling flows. For this reason, the data of > 100 (see Figures 4(a) and 4(c)) and −0.8 < , , < 0.4 (see Figure 5(d) or equivalently 0.158 < , , < 2.5) was extracted and the locations for these data in the background of vorticities of and 0.59 were shown.

Figure 6 shows that the locations of extremely anisotropic turbulence dissipation (e.g., > 100, designated by colored spheres) are concentrated in the central region of the jet (Figure 6(a)), corresponding to strong twisted vortices when the swirl level is low. On the other hand, the locations of nearly isotropic turbulence dissipation (0.158 < , , < 2.5) are dispersed in the peripheral region of the jet (Figure 6(b)), corresponding to the region with small-scale and weak turbulence. Moreover, the locations of extremely anisotropic turbulence dissipation are concentrated in the central region of vortex breakdown (Figure 6(c)), corresponding to strong swirling large-scale vortex structure when the swirl level is large. The locations of nearly isotropic turbulence dissipation are dispersed in the peripheral region of the vortex breakdown, such as strong small-scale vortices (Figure 6(d)).

(a)

(b)

(c)

(d)

Figure 6: Locations of anisotropic turbulence dissipation ( > 100) for (a) and (c) and locations of nearly isotropic turbulence dissipation (, , or 0.158 < , , < 2.5) for (b) and (d), respectively, in the background of vorticities.

Figure 6 conclusively shows that strong anisotropic turbulence dissipation occurs concentratively in the vortex breakdown region or is closely related to the large-scale vortex structure. On the other hand, nearly isotropic turbulence dissipation occurs dispersively in the peripheral region of the strong swirling flows, that is, closely related to small-scale vortices.

4. Conclusion

This work was carried out to investigate the physical aspects of anisotropic turbulent motions and dissipations in swirling flows. Based on the observation and analysis of the DNS results, the following results were found.

(1)Turbulent swirling jet flows have evident coherent structures relevant to the remarkable VB phenomenon and recirculation of the flows, which dominates anisotropic turbulent motions and dissipations.(2)The evidently nonzero components of dissipation tensor relevant to validates the occurrence of anisotropic turbulent fluctuations and motions in swirling flows.(3)With diagonalized dissipation tensor, the relation and distribution of the three eigenvalues show strong anisotropic characteristics of turbulent energy dissipation, especially in strong swirling flows.(4)Based on the probability density functions, turbulent dissipation of swirling flows is anisotropic in most regions, but an isotropic subrange of turbulent dissipation still exists, especially for low swirling flows. On the contrary, the isotropic subrange of PDF is reduced for strong swirling flows. The power-law form of PDF for the isotropic subrange was obtained through linear regression.(5)More importantly, strong anisotropic turbulence dissipation occurs concentratively in the vortex breakdown region or is closely related to the large-scale vortex structure, whereas the nearly isotropic turbulence dissipation occurs dispersively in the peripheral region of the strong swirling flows, that is, closely related to small-scale vortices.

Nomenclatures

Scalars

, , : | Coefficients
: | The jet diameter at the inlet
, : | Low and high limiting value
: | Turbulent kinetic energy
: | Number density
: | Fluid pressure
, : | Radius, radial center
Re: | Reynolds number
: | Swirl number
: | Time
, : | Axial and azimuthal velocity of fluids
, : | Velocity of fluids
, : | Mean velocities of fluids
: | The inflow velocity
, : | Spatial variables
, , : | Regression lines
: | Kronecker function
: | Time step
: | Mesh spacing
, : | Turbulent energy dissipation tensor
, , : | Three eigenvalues of
: | The Kolmogorov length scale
: | Inflow momentum thickness
: | Kinematic viscosity of fluids
: | Turbulent viscosity of fluids
: | Reynolds stress tensor.

Operators

: | Hamiltonian operator,
: | Laplace operator
: | Partial derivative
: | Assemble averaging process
: | Probability density function
: | Infinitesimal increment
: | Summation operator.

Subscripts

: | Center
: | High
: | Index
: | Low
0: | Initial
: | Turbulence
: | Axial direction
: | Azimuthal direction
: | Fluctuation value
: | Eigenvalue.

Abbreviations

DNS: | Direct numerical simulation
K-H: | Kelvin-Helmholtz
PDF: | Probability density function
VB: | Vortex breakdown.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful for the support of this study by the National Natural Science Foundation of China (51106180), the Foundation for the Author of National Excellent Doctoral Dissertation of PR China (FANEDD, Grant no. 201438), and the China Postdoctoral Science Foundation (2013M540964).

References

C. Chen and D.-J. Sun, “Numerical investigation of a swirling flow under the optimal perturbation,” Journal of Hydrodynamics Series B, vol. 22, no. 5, supplement 1, pp. 237–241, 2010.

Z. Wang and S. Q. Chen, “Structures of confined vortex breakdown in constant diameter pipe flow,” Journal of Hydrodynamics, vol. 21, no. 3, pp. 341–346, 2009.

S.-Y. Jin, Y.-Z. Liu, W.-Z. Wang, Z.-M. Cao, and H. S. Koyama, “Numerical evaluation of two-fluid mixing in a swirl micro-mixer,” Journal of Hydrodynamics, vol. 18, no. 5, pp. 542–546, 2006.

V. Shtern, M. Goldshtik, and F. Hussain, “Generation of swirl due to symmetry breaking,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 49, no. 4, pp. 2881–2886, 1994.

B. di Pierro and M. Abid, “Instabilities of variable-density swirling flows,” Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, vol. 82, no. 4, Article ID 046312, 2010.

M. R. Ruith, P. Chen, E. Meiburg, and T. Maxworthy, “Three-dimensional vortex breakdown in swirling jets and wakes: direct numerical simulation,” Journal of Fluid Mechanics, vol. 486, pp. 331–378, 2003.

L. Li, X. Y. Qiu, S. Jin, J. Xiao, and S.-Y. Gong, “Weakly swirling turbulent flow in turbid water hydraulic separation device,” Journal of Hydrodynamics B, vol. 20, no. 3, pp. 347–355, 2008.

N. Gui, J. Fan, K. Cen, and S. Chen, “A direct numerical simulation study of coherent oscillation effects of swirling flows,” Fuel, vol. 89, no. 12, pp. 3926–3933, 2010.

N. Gui, J. Fan, and S. Chen, “Numerical study of particle-vortex interaction and turbulence modulation in swirling jets,” Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, vol. 82, no. 5, Article ID 056323, 2010.

N. Gui, J. Fan, and S. Chen, “The effects of flow structure and particle mass loading on particle dispersion in particle-laden swirling jets,” Physics Letters, Section A: General, Atomic and Solid State Physics, vol. 375, no. 4, pp. 839–844, 2011.

B. E. Launder and D. B. Spalding, “The numerical computation of turbulent flows,” Computer Methods in Applied Mechanics and Engineering, vol. 3, no. 2, pp. 269–289, 1974.

L. N. Jones, P. H. Gaskell, H. M. Thompson, X. J. Gu, and D. R. Emerson, “Anisotropic, isothermal, turbulent swirling flow in a complex combustor geometry,” International Journal for Numerical Methods in Fluids, vol. 47, no. 10-11, pp. 1053–1059, 2005.

D. R. Radenković, J. M. Burazer, and Đ. M. Novković, “Anisotropy analysis of turbulent swirl flow,” FME Transactions, vol. 42, no. 1, pp. 19–25, 2014.

A. Escue and J. Cui, “Comparison of turbulence models in simulating swirling pipe flows,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 34, no. 10, pp. 2840–2849, 2010.

D. X. Fu and Y. W. Ma, “A high order accurate difference scheme for complex flow fields,” Journal of Computational Physics, vol. 134, no. 1, pp. 1–15, 1997.

S. K. Lele, “Compact finite difference schemes with spectral-like resolution,” Journal of Computational Physics, vol. 103, no. 1, pp. 16–42, 1992.

A. Jameson and W. Schmidt, “Some recent developments in numerical methods for transonic flows,” Computer Methods in Applied Mechanics and Engineering, vol. 51, no. 1–3, pp. 467–493, 1985.

U. A. Krishnaiah, R. Manohar, and J. W. Stephenson, “Fourth-order finite difference methods for three-dimensional general linear elliptic problems with variable coefficients,” Numerical Methods for Partial Differential Equations, vol. 3, no. 3, pp. 229–240, 1987.

N. Gui, J. Yan, Z. Li, and J. Fan, “Direct numerical simulation of confined swirling jets,” International Journal of Computational Fluid Dynamics, vol. 28, no. 1-2, pp. 76–88, 2014.

N. Gui, J. R. Fan, and S. Chen, “Numerical study of particle-particle collision in swirling jets: a DEM-DNS coupling simulation,” Chemical Engineering Science, vol. 65, no. 10, pp. 3268–3278, 2010.

P. Moin and K. Mahesh, “Direct numerical simulation: a tool in turbulence research,” Annual Review of Fluid Mechanics, vol. 30, pp. 539–578, 1998.

I. Orlanski, “A simple boundary condition for unbounded hyperbolic flows,” Journal of Computational Physics, vol. 21, no. 3, pp. 251–269, 1976.

R. W. Johnson, The Handbook of Fluid Dynamics, CRC Press, New York, NY, USA, 1998.