Multi-relaxation-time lattice Boltzmann simulations of lid driven flows using graphics processing unit

Large eddy simulation(LES) using the Smagorinsky eddy viscosity model is added to the two-dimensional nine velocity components(D2Q9) lattice Boltzmann equation(LBE) with multi-relaxation-time(MRT) to simulate incompressible turbulent cavity flows with the Reynolds numbers up to 1×10~7. To improve the computation efficiency of LBM on the numerical simulations of turbulent flows, the massively parallel computing power from a graphic processing unit(GPU) with a computing unified device architecture(CUDA) is introduced into the MRT-LBE-LES model. The model performs well, compared with the results from others, with an increase of 76 times in computation efficiency.It appears that the higher the Reynolds numbers is, the smaller the Smagorinsky constant should be, if the lattice number is fixed. Also, for a selected high Reynolds number and a selected proper Smagorinsky constant, there is a minimum requirement for the lattice number so that the Smagorinsky eddy viscosity will not be excessively large.     

Investigation of the turbulent model for pressure fluctuations with spectral theory

The pressure fluctuations in turbulent shear flows are investigated with the theory of spectral analysis.An expression for pressure spectra is analytically derived in terms of velocity spectra.This derivation is based on a formal solution of the Navier-Stokes equation and quasi-normal assumption to express the third and fourth order velocity correlations in terms of double velocity correlation.Then,a turbulent model for the computation of pressure fluctuation intensity with Renolds stress and mean flow velocity gradients is established.The turbulent constants in the model are calculated from the assumptions about the general behaviour of velocity spectra in high Renolds number flows.Comparison with direct simulation of turbulent boundary layer is made.It is found that the turbulent-turbulent,cross correlation,and turbulent-shear source terms for mean square value of pressure fluctuation are about the same magnitude.     

Effect of particles on turbulent thermal field of channel flow with different Prandtl numbers

The direct numerical simulation(DNS) of heat transfer in a fully developed non-isothermal particle-laden turbulent channel flow is performed.The focus of this paper is on the modulation of the particles on turbulent thermal statistics in the particle-laden flow with three Prandtl numbers(P r = 0.71,1.5,and 3.0) and a shear Reynolds number(Reτ = 180).Some typical thermal statistics,including normalized mean temperature and their fluctuations,turbulent heat fluxes,Nusselt number and so on,are analyzed.The results show that the particles have less effects on turbulent thermal fields with the increase of Prandtl number.Two reasons can explain this.First,the correlation between fluid thermal field and velocity field decreases as the Prandtl number increases,and the modulation of turbulent velocity field induced by the particles has less influence on the turbulent thermal field.Second,the heat exchange between turbulence and particles decreases for the particle-laden flow with the larger Prandtl number,and the thermal feedback of the particles to turbulence becomes weak.     

Large eddy simulation of hot and cold fluids mixing in a T-junction for predicting thermal fluctuations

Temperature fluctuations in a mixing T-junction have been simulated on the FLUENT platform using the large eddy simulation （LES） turbulent flow model and a sub-grid scale Smagorinsky-Lilly model. The normalized mean and root mean square temperatures for describing time-averaged temperature and temperature fluctuation intensity, and the velocity are obtained. The power spectrum densities of temperature fluctuations, which are key parameters for thermal fatigue analysis and lifetime evaluation, are analyzed. Simulation results are consistent with experimental data published in the literature, showing that the LES is reliable. Several mixing processes under different conditions are simulated in order to analyze the effects of varying Reynolds number and Richardson number on the mixing course and thermal fluctuations.     

Effects of viscous dissipation on thermally developing forced convection in a porous saturated circular tube with an isoflux wall

The viscous dissipation effect on forced convection in a porous saturated circular tube with an isoflux wall is investigated on the basis of the Brinkman flow model. For the thermally developing region, a numerical study is reported while a perturbation analysis is presented to find expressions for the temperature profile and the Nusselt number for the fully developed region. The fully developed Nusselt number found by numerical solution for the developing region is compared with that of asymptotic analysis and a good degree of agreement is observed.     

Similarity solutions of vertical plane wall plume based on finite analytic method

The turbulent flow of vertical plane wall plume with concentration variation was studied with the finite analytical method. The k-epsilon model with the effect of buoyancy on turbulent kinetic energy and its dissipation rate was adopted. There were similarity solutions in the uniform environment for the system of equations including the equation of continuity, the equation of momentum along the flow direction and concentration, and equations of k, epsilon. The finite analytic method was applied to obtain the similarity solution. The calculated data of velocity, relative density difference, the kinetic energy of turbulence and its dissipation rate distribution for vertical plane plumes are in good agreement with the experimental data at the turbulent Schmidt number equal to 1.0. The variations of their maximum value along the direction of main flow were also given. It shows that the present model is good, i.e., the effect of buoyancy on turbulent kinetic energy and its dissipation rate should be taken into account, and the finite analytic method is effective.     

Propagation of wave at the boundary surface of transversely isotropic thermoelastic material with voids and isotropic elastic half-space

The purpose of this research is to study the effect of voids on the surface wave propagation in a layer of a transversely isotropic thermoelastic material with voids lying over an isotropic elastic half-space. The frequency equation is derived after developing a mathematical model for welded and smooth contact boundary conditions. The dispersion curves giving the phase velocity and attenuation coefficient via wave number are plotted graphically to depict the effects of voids and anisotropy for welded contact boundary conditions. The specific loss and amplitudes of the volume fraction field, the normal stress, and the temperature change for welded contact are obtained and shown graphically for a particular model to depict the voids and anisotropy effects. Some special cases are also deduced from the present investigation.     

A Boussinesq model with alleviated nonlinearity and dispersion

The classical Boussinesq equation is a weakly nonlinear and weakly dispersive equation, which has been widely applied to simulate wave propagation in off-coast shallow waters. A new form of the Boussinesq model for an uneven bottoms is derived in this paper. In the new model, nonlinearity is reduced without increasing the order of the highest derivative in the differential equations. Dispersion relationship of the model is improved to the order of Pade （2,2） by adjusting a parameter in the model based on the long wave approximation. Analysis of the linear dispersion, linear shoaling and nonlinearity of the present model shows that the performances in terms of nonlinearity, dispersion and shoaling of this model are improved. Numerical results obtained with the present model are in agreement with experimental data.     

Distribution of non-spherical nanoparticles in turbulent flow of ventilation chamber considering fluctuating particle number density

The Reynolds-averaged general dynamic equation(RAGDE)for the nanoparticle size distribution function is derived,including the contribution to particle coagulation resulting from the fluctuating concentration.The equation together with that of a turbulent gas flow is solved numerically in the turbulent flow of a ventilation chamber with a jet on the wall based on the proposed model relating the fluctuating coagulation to the gradient of mean concentration.Some results are compared with the experimental data.The results show that the proposed model relating the fluctuating coagulation to the gradient of mean concentration is reasonable,and it is necessary to consider the contribution to coagulation resulting from the fluctuating concentration in such a flow.The changes of the particle number concentration M₀ and the geometric mean diameter dg are more obvious in the core area of the jet,but less obvious in other areas.With the increase in the initial particle number concentration m00,the values of M₀ and the standard deviation of the particle sizeσdecrease,but the value of d_g increases.The decrease in the initial particle diameter leads to the reduction of M₀ andσand the increase in d_g.With the increase in the Reynolds number,particles have few chances of collision,and hence the coagulation rate is reduced,leading to the increase in M₀ andσand the decrease in d_g.     

Experimental and statistical assessments of the mechanical strength reliability of gamma alumina catalyst supports

The mechanical strength of solid catalysts is considered an important factor in terms of ensuring the reliable performance of industrial reactors.In this work,a pelletizing method was used to form gamma alumina support for catalysts.Response surface methodology(RSM) was employed to analyze and model the effects of various manufacturing parameters on the crushing strength of the supports.These parameters were binder concentration,compaction pressure,calcination temperature,and drying mode.The suggested model was verified by applying an analysis of variance to assess its validity with regard to crushing strength.The mechanical reliability of various supports was also determined by calculating their Weibull modulus values through linear regression of the Weibull equation.The material with the highest mechanical strength reliability will have both a high mean crushing strength and a high Weibull modulus,and the best values obtained for a support in this work were 70.7 MPa and 6.63,respectively.The conditions used to form this sample were:20mass%binder concentration,861 MPa compaction pressure,466 ℃ calcination temperature,and gentle drying.     

Effects of temperature-dependent viscosity variation on entropy generation,heat and fluid flow through a porous-saturated duct of rectangular cross-section

Effect of temperature-dependent viscosity on fully developed forced convection in a duct of rectangular cross-section occupied by a fluid-saturated porous medium is investigated analytically. The Darcy flow model is applied and the viscosity-temperature relation is assumed to be an inverse-linear one. The case of uniform heat flux on the walls, i.e. the H boundary condition in the terminology of Kays and Crawford [12], is treated. For the case of a fluid whose viscosity decreases with temperature, it is found that the effect of the variation is to increase the Nusselt number for heated walls. Having found the velocity and the temperature distribution, the second law of thermodynamics is invoked to find the local and average entropy generation rate. Expressions for the entropy generation rate, the Bejan number, the heat transfer irreversibility, and the fluid flow irreversibility are presented in terms of the Brinkman number, the Peclet number, the viscosity variation number, the dimensionless wall heat flux, and the aspect ratio (width to height ratio). These expressions let a parametric study of the problem based on which it is observed that the entropy generated due to flow in a duct of square cross-section is more than those of rectangular counterparts while increasing the aspect ratio decreases the entropy generation rate similar to what previously reported for the clear flow case by Ratts and Raut [14].     

MHD flow and heat transfer from continuous surface in uniform free stream of non-Newtonian fluid

An analysis is carried out to study the steady flow and heat transfer charac- teristics from a continuous flat surface moving in a parallel free stream of an electrically conducting non-Newtonian viscoelastic fluid.The flow is subjected to a transverse uni- form magnetic field.The constitutive equation of the fluid is modeled by that for a second grade fluid.Numerical results are obtained for the distribution of velocity and temperature profiles.The effects of various physical parameters like viscoelastic param- eter,magnetic parameter and Prandtl number on various momentum and heat transfer characteristics are discussed in detail and shown graphically.     

Numerical simulation of pedestrian flow past a circular obstruction

In this paper, a revisiting Hughes' dynamic continuum model is used to investigate and predict the essential macroscopic characteristics of pedestrian flow, such as flow, density and average speed, in a two dimensional continuous walking facility scattered with a circular obstruction. It is assumed that pedestrians prefer to walk a path with the lowest instantaneous travel cost from origin to destination, under the consideration of the current traffic conditions and the tendency to avoid a high-density region and an obstruction. An algorithm for the pedestrian flow model is based on a cellcentered finite volume method for a scalar conservation law equation, a fast sweeping method for an Eikonal-type equation and a second-order TVD Runge-Kutta method for the time integration on unstructured meshes. Numerical results demonstrate the effectiveness of the algorithm. It is verified that density distribution of pedestrian flow is influenced by the position of the obstruction and the path-choice behavior of pedestrians.     

Green’s function solution for transient heat conduction in annular fin during solidification of phase change material

The Green’s function method is applied for the transient temperature of an annular fin when a phase change material (PCM) solidifies on it. The solidification of the PCMs takes place in a cylindrical shell storage. The thickness of the solid PCM on the fin varies with time and is obtained by the Megerlin method. The models are found with the Bessel equation to form an analytical solution. Three different kinds of boundary conditions are investigated. The comparison between analytical and numerical solutions is given. The results demonstrate that the significant accuracy is obtained for the temperature distribution for the fin in all cases.     

Nonlinear convection in a non-Darcy porous medium

In this paper, the natural convection in a non-Darcy porous medium is studied using a temperature-concentration-dependent density relation. The effect of the two parameters responsible for the nonlinear convection is analyzed for different values of the inertial parameter, dispersion parameters, Rayleigh number, Lewis number, Soret number, and Dufour number. In the aiding buoyancy, the tangential velocity increases steeply with an increase in the nonlinear temperature parameter and the nonlinear concentration parameter when the inertial effect is zero. However, when the inertial effect is non-zero, the effect of the nonlinear temperature parameter and the nonlinear concentration parameter on the tangential velocity is marginal. The concentration distribution varies appreciably and spreads in different ranges for different values of the double dispersion parameters, the inertial effect parameter, and also for the parameters which control the nonlinear temperature and the nonlinear concentration. Heat and mass transfer varies extensively with an increase in the nonlinear temperature parameter and the nonlinear concentration parameter depending on Dacry and non-Darcy porous media. The variation in heat and mass transfer when all the effects, i.e., the inertial effect, double dispersion ef- fects, and Soret and Dufour effects, are simultaneously zero and non-zero. The combined effects of the nonlinear temperature parameter, the nonlinear concentration parameter and buoyancy are analyzed. The effect of the nonlinear temperature parameter and the nonlinear concentration parameter and also the cross diffusion effects on heat and mass transfer are observed to be more in Darcy porous media compared with those in non- Darcy porous media. In the opposing buoyancy, the effect of the temperature parameter is to increase the heat and mass transfer rate, whereas that of the concentration parameter is to decrease.     

Large eddy simulation of vertical turbulent jets under JONSWAP waves

The effect of random waves on vertical plane turbulent jets is studied numerically and the mechanism behind the interaction of the jet and waves is analyzed. The large eddy simulation method is used and the σ-coordinate system is adopted. Turbulence is modeled by a dynamic coherent eddy model. The σ-coordinate transformation is introduced to map the irregular physical domain with a wavy free surface and an uneven bottom onto a regular computational domain. The fractional step method is used to solve the fil...     

TEMPERATURE PROFILES OF LOCAL THERMAL NONEQUILIBRIUM FOR THERMAL DEVELOPING FORCED CONVECTION IN POROUS MEDIUM PARALLEL PLATE CHANNEL

Based on the two-energy equation model, taking into account viscous dissipation due to the interaction between solid skeleton and pore fluid flow, temperature expressions of the solid skeleton and pore fluid flow are obtained analytically for the thermally developing forced convection in a saturated porous medium parallel plate channel, with walls being at constant temperature. It is proved that the temperatures of the two phases for the local thermal nonequilibrium approach to the temperature derived from the one-energy equation model for the local thermal equilibrium when the heat exchange coefficient goes to infinite. The temperature profiles are shown in figures for different dimensionless parameters and the effects of the parameters on the local thermal nonequilibrium are revealed by parameter study.     

Lattice Boltzmann simulations of turbulent shear flow between parallel porous walls

The effects of two parallel porous walls are investigated, consisting of the Darcy number and the porosity of a porous medium, on the behavior of turbulent shear flows as well as skin-friction drag. The turbulent channel flow with a porous surface is directly simulated by the lattice Boltzmann method （LBM）. The Darcy-Brinkman- Forcheimer （DBF） acting force term is added in the lattice Boltzmann equation to simu- late the turbulent flow bounded by porous walls. It is found that there are two opposite trends （enhancement or reduction） for the porous medium to modify the intensities of the velocity fluctuations and the Reynolds stresses in the near wall region. The parametric study shows that flow modification depends on the Darcy number and the porosity of the porous medium. The results show that, with respect to the conventional impermeable wall, the degree of turbulence modification does not depend on any simple set of param- eters obviously. Moreover, the drag in porous wall-bounded turbulent flow decreases if the Darcy number is smaller than the order of O（10-4） and the porosity of porous walls is up to 0.4.     

Stability analysis of a noise control system in a duct by using delay differential equation

The paper deals with the criteria for the closed- loop stability of a noise control system in a duct. To study the stability of the system, the model of delay differential equation is derived from the propagation of acoustic wave governed by a partial differential equation of hyperbolic type. Then, a simple feedback controller is designed, and its closed- loop stability is analyzed on the basis of the derived model of delay differential equation. The obtained criteria reveal the influence of the controller gain and the positions of a sensor and an actuator on the closed-loop stability. Finally, numerical simulations are presented to support the theoretical results.     

EXPERIMENTAL STUDY OF MEASUREMENT FOR DISSIPATION RATE SCALING EXPONENT IN HEATED WALL TURBULENCE

Experimental investigations have been devoted to the study of scaling law of coarse-grained dissipation rate structure function for velocity and temperature fluctuation of non-isotropic and inhomogeneous turbulent flows at moderate Reynolds number. Much attention has been paid to the case of turbulent boundary layer, which is typically the non-istropic and inhomogeneous trubulence because of the dynamically important existence of organized coherent structure burst process in the near wall region . Longitudinal velocity and temperature have been measured at different vertical positions in turbulent boundary layer over a heated and unheated flat plate in a wind tunnel using hot wire anemometer. The influence of non-isotropy and inhomogeneity and heating the wall on the scaling law of the dissipation rate structure function is studied because of the existence of organized coherent structure burst process in the near wall region . The scaling law of coarse-grained dissipation rate structure function is foun