The unsteady magnetohydrodynamical(MHD)free convection flow of an incompressible,electrically conducting hybrid nanofluid within a vertical cylindrical geometry is investigated,incorporating the effects of thermal rad...The unsteady magnetohydrodynamical(MHD)free convection flow of an incompressible,electrically conducting hybrid nanofluid within a vertical cylindrical geometry is investigated,incorporating the effects of thermal radiation,viscous dissipation,and internal heat generation.The system is subjected to a time-periodic boundary temperature condition.The Laplace and finite Hankel transforms are used to derive the exact solutions for the velocity and temperature distributions.The effects of various key physical parameters,including the Richardson number,the Eckert number,the radiation parameter,the heat source parameter,and the nanoparticle volume fraction,are considered.The numerical results reveal that increasing the volume fraction significantly enhances the thermal conductivity and temperature,while the magnetic field intensity and viscous dissipation strongly influence the fluid motion and heat transport.Additionally,the pulsating boundary conditions produce distinct oscillatory behaviors in both the velocity and temperature fields.These findings provide important insights into optimizing the heat transfer performance in cylindrical systems such as electronic cooling modules and energy storage devices operating under dynamic thermal conditions.展开更多
The motion of an elliptical rigid particle in a lid-driven cavity flow was numerically simulated using the immersed boundary lattice Boltzmann method(IB-LBM).The effects of the particle's initial orientation angle...The motion of an elliptical rigid particle in a lid-driven cavity flow was numerically simulated using the immersed boundary lattice Boltzmann method(IB-LBM).The effects of the particle's initial orientation angle,initial position,aspect ratio,and size on the motion characteristics were investigated.The computational results indicate that the particle's motion undergoes two distinct stages:a starting stage that involves moving from the release position to a limit cycle,and a periodic stage that involves moving on the limit cycle.The initial orientation of the particle has a minimal impact on both stages of motion.In contrast,the time it takes for the particle to reach the limit cycle may vary depending on the release position.Furthermore,particles with a larger aspect ratio exhibit a greater maximum velocity magnitude;an increase in particle size causes the particle trajectory to contract more toward the center of the cavity,decreasing the centrifugal force experienced by the particle.展开更多
By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improv...By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number （Re） flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.展开更多
Steady,laminar mixed convection inside a lid-driven square cavity filled with nanofluid is investigated numerically.We consider the case where the right and left walls are moving downwards and upwards respectively and...Steady,laminar mixed convection inside a lid-driven square cavity filled with nanofluid is investigated numerically.We consider the case where the right and left walls are moving downwards and upwards respectively and maintained at different temperatures while the other two horizontal ones are kept adiabatic and impermeable.The set of nonlinear coupled governing mass,momentum,and energy equations are solved using an extensively validated and a highly accurate finite difference method of fourth-order.Comparisons with previously conducted investigations on special configurations are performed and show an excellent agreement.Meanwhile,attention is focused on the heat transfer enhancement when different nano-particles:Cu,Ag,Al2O3,TiO2 and Fe3O4 are incorporated separately in different base fluids such as:Water,Ethylene-glycol,Methanol and Kerosene oil.In this framework,the numerical results related to several mixtures are presented and concern flow pattern and heat transfer curves for various values of Richardson number[Ri=0.1,1 and 10].It turns out that the choice of the efficient binary mixture for an optimal heat transfer depends not only on the thermophysical properties of the nanofluids but also on the range of the Richardson number.Special attention is devoted to shedding light on the effect of the shape of the nanoparticles on the heat transfer in the case of Water-Ag nanofluid.It is concluded that the spherical shape is more suitable for a better heat transfer enhancement in comparison to the cylindrical ones.展开更多
To date, there are very few studies on the second Hopf bifurcation in a driven square cavity, although there are intensive investigations focused on the first Hopf bifurcation in literature, due to the difficulties of...To date, there are very few studies on the second Hopf bifurcation in a driven square cavity, although there are intensive investigations focused on the first Hopf bifurcation in literature, due to the difficulties of theoretical analyses and numerical simulations. In this paper, we study the characteristics of the second Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme recently developed by us. We numerically identify the critical Reynolds number of the second Hopf bifurcation located in the interval of(11093.75, 11094.3604) by bisection. In addition, we find that there are two dominant frequencies in its spectral diagram when the flow is in the status of the second Hopf bifurcation, while only one dominant frequency is identified if the flow is in the first Hopf bifurcation via the Fourier analysis. More interestingly, the flow phase portrait of velocity components is found to make transition from a regular elliptical closed form for the first Hopf bifurcation to a non-elliptical closed form with self-intersection for the second Hopf bifurcation. Such characteristics disclose flow in a quasi-periodic state when the second Hopf bifurcation occurs.展开更多
A streamline upwind/Petrov-Galerkin （SUPG） finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for-...A streamline upwind/Petrov-Galerkin （SUPG） finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.展开更多
To date,there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity,due to the difficulties in theoretical analysis and numerical simulations.In this paper,we study the ch...To date,there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity,due to the difficulties in theoretical analysis and numerical simulations.In this paper,we study the characteristics of the third Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme rectently developed by us.We numerically identify the critical Reynolds number of the third Hopf bifurcation located in the interval of(13944.7021,13946.5333)by the method of bisection.Through Fourier analysis,it is discovered that the flow becomes chaotic with a characteristic of period-doubling bifurcation when the Reynolds number is beyond the third bifurcation critical interval.Nonlinear time series analysis further ascertains the flow chaotic behaviors via the phase diagram,Kolmogorov entropy and maximal Lyapunov exponent.The phase diagram changes interestingly from a closed curve with self-intersection to an unclosed curve and the attractor eventually becomes strange when the flow becomes chaotic.展开更多
This research work numerically analyzes 2D,steady state,mixed convective heat transfer for Newtonian fluids in lid driven square enclosure with centered triangular block（blockage—10%or 30%）maintained either at the ...This research work numerically analyzes 2D,steady state,mixed convective heat transfer for Newtonian fluids in lid driven square enclosure with centered triangular block（blockage—10%or 30%）maintained either at the constant wall temperature or constant heat flux thermal conditions.The fluid flow in the enclosure is initiated by top moving wall in＋x-direction,while all other walls are stationary.The top and bottom walls are thermally insulated.In particular,the governing field equations are solved for range of governing parameters such as,Reynolds number（1–1000）,Prandtl number（1–100）,and Grashof number展开更多
Hydrodynamic mixed convection in a lid-driven hexagonal cavity with corner heater is numerically simulated in this paper by employing finite element method. The working fluid is assigned as air with a Prandtl num-ber ...Hydrodynamic mixed convection in a lid-driven hexagonal cavity with corner heater is numerically simulated in this paper by employing finite element method. The working fluid is assigned as air with a Prandtl num-ber of 0.71 throughout the simulation. The left and right walls of the hex-agonal cavity are kept thermally insulated and the lid moves top to bottom at a constant speed U0. The top left and right walls of the enclosure are maintained at cold temperature Tc. The bottom right wall is considered with a corner heater whereas the bottom remaining part is adiabatic and inside the cavity a square shape heated block Th. The focus of the work is to investigate the effect of Hartmann number, Richardson number, Grashof number and Reynolds number on the fluid flow and heat transfer characteristics inside the enclosure. A set of graphical results is presented in terms of streamlines, isotherms, local Nusselt number, velocity profiles, temperature profiles and average Nusselt numbers. The results reveal that heat transfer rate increases with increasing Richardson number and Hartmann number. It is also observed that, Hartmann number is a good control parameter for heat transfer in fluid flow in hexagonal cavity.展开更多
In the case of lid-driven deep cavity flow, the effects of different resolutions of local grid refinement have been studied in the frame of multiple relaxation times (MRT) lattice Boltzmann method (LBM). In all the ca...In the case of lid-driven deep cavity flow, the effects of different resolutions of local grid refinement have been studied in the frame of multiple relaxation times (MRT) lattice Boltzmann method (LBM). In all the cases, the aspect ratio and Reynolds number are set as 1.5 and 3.200, respectively. First, the applied method is validated by comparing it with two reported works, with which agreements are reached. Then, six separate degrees of local grid refinement at the upper left corner, i.e. purely coarse grid, including 1/64, 1/32, 1/16, 1/8, 1/4 refinements of the lattice number in the width direction have been studied in detail. The results give the following indications:① The refinement degrees lower than 1/8 produce similar results;② For single corner refinement, 1/4 refinement is adequate for clearing the noises in the singularity zone to a large extent;③ New noise around the interface between coarse and fine zones are introduced by local grid refinement. Finally, refinement of entire subzone neighboring the lid is examined to avoid introducing new noises and it has been found effective.展开更多
This study aims to analyze mixed convection in a square cavity with two moving vertical walls by finite volume method.The cavity filled with Non-Newtonian fluid of Bingham model is heated from below and cooled by the ...This study aims to analyze mixed convection in a square cavity with two moving vertical walls by finite volume method.The cavity filled with Non-Newtonian fluid of Bingham model is heated from below and cooled by the other walls.This study has been conducted for certain parameters of Reynolds number(Re=1-100),Richardson number(Ri=1-20),Prandtl number(Pr=1-500),and Bingham number has been studied from 0 to 10.The results indicate that the increase in yield stress drops the heat transfer and the flow become flatter,while increasing Reynolds number augments it.The convective transport is dominant when increasing Richardson number which leads to enhance heat transfer in the cavity for both Newtonian and Non-Newtonian fluid.A correlation of Nusselt number is given in function of different parameters.展开更多
The concentration and velocity fields of two refractive index matched miscible shear-thinning fluids in a lid-driven cavity were investigated by using planar laser-induced fluorescence and particle image velocimetry,a...The concentration and velocity fields of two refractive index matched miscible shear-thinning fluids in a lid-driven cavity were investigated by using planar laser-induced fluorescence and particle image velocimetry,as well by computational fluid dynamics.Quantitative analyses show that the results obtained by flow simulations with the species transport model are in good agreement with the experimental results.The effects of different parameters were studied by using the intensity of segregation.For two fluids with the same rheological parameters,the relative amounts of liquids H_(1)/H and the power-law index n dominate the mixing process while the Reynolds number Re plays a marginal role.As for two fluids with density difference,buoyancy has significant influence on the mixing process.The dimensionless group Ar/Re(redefined such as to include shear thinning behavior)is proposed for assessing the effect of buoyancy and rheological properties on the mixing of miscible shear-thinning fluids.展开更多
Numerical simulation of the bifurcation of Bingham fluid streamline topologies in rectangular double-lid-driven cavity, with varying aspect (height to width) ratio A, is presented. The lids on the top and bottom move ...Numerical simulation of the bifurcation of Bingham fluid streamline topologies in rectangular double-lid-driven cavity, with varying aspect (height to width) ratio A, is presented. The lids on the top and bottom move at the same speed but in opposite directions so that symmetric flow patterns are generated. Similar to the Newtonian case, bifurcations occur as the aspect ratio decreases. Special to Bingham fluids, the non-Newtonian indicator, Bingham number B, also governs the bifurcation besides the bifurcation parameter A.展开更多
An investigation on flow and heat transfer due to mixed convection, in a lid-driven rectangular cavity filled with Cu- water nanofluids and submitted to uniform heat flux along with its vertical short sides, has been ...An investigation on flow and heat transfer due to mixed convection, in a lid-driven rectangular cavity filled with Cu- water nanofluids and submitted to uniform heat flux along with its vertical short sides, has been conducted numerically by solving the full governing equations with the finite volume method and the SIMPLER algorithm. In the case of a slender enclosure, these equations are considerably reduced by using the parallel flow concept. Solutions, for the flow and temperature fields, and the heat transfer rate, have been obtained depending on the governing parameters, which are the Reynolds, the Richardson numbers and the solid volume fraction of nanoparticles. A perfect agreement has been found between the results of the two approaches for a wide range of the abovementioned parameters. It has been shown that at low and high Richardson numbers, the convection is ensured by lid and buoyancy-driven effects, respectively, whereas between these extremes, both mechanisms compete. Moreover, the addition of Cu-nanoparticles, into the pure water, has been seen enhancing and degrading heat transfer by lid and buoyancy-driven effects, respectively.展开更多
This study is devoted to the computational fluid dynamics (CFD) modeling of steady laminar mixed convection flow and heat transfer in lid driven cavity (10 ≤ Re ≤ 1000). The ratio of the height to the width of the c...This study is devoted to the computational fluid dynamics (CFD) modeling of steady laminar mixed convection flow and heat transfer in lid driven cavity (10 ≤ Re ≤ 1000). The ratio of the height to the width of the cavity is ranged over H/L = 0.5 to 1.5. The governing equations are solved using commercial finite volume package FLUENT to visualize the nature of the flow and estimate the heat transfer inside the cavity for different aspect ratio. The simulation results are presented in terms of average Nusselt number of the hot wall, velocity profile, and temperature contours. It was found that the average Nusselt number inside the cavity is strongly governed by the aspect ratio as well as the Reynolds number. A parametric study is conducted to demonstrate the effect of aspect ratio on the flow and heat transfer characteristics. It is found that heat transfer enhancement was obtained by decreasing the aspect ratio and/or increasing the Reynolds number.展开更多
This paper investigates a polygonal finite element(PFE)to solve a two-dimensional(2D)incompressible steady fluid problem in a cavity square.It is a well-known standard benchmark(i.e.,lid-driven cavity flow)-to evaluat...This paper investigates a polygonal finite element(PFE)to solve a two-dimensional(2D)incompressible steady fluid problem in a cavity square.It is a well-known standard benchmark(i.e.,lid-driven cavity flow)-to evaluate the numerical methods in solving fluid problems controlled by the Navier-Stokes(N-S)equation system.The approximation solutions provided in this research are based on our developed equal-order mixed PFE,called Pe1Pe1.It is an exciting development based on constructing the mixed scheme method of two equal-order discretisation spaces for both fluid pressure and velocity fields of flows and our proposed stabilisation technique.In this research,to handle the nonlinear problem of N-S,the Picard iteration scheme is applied.Our proposed method’s performance and convergence are validated by several simulations coded by commercial software,i.e.,MATLAB.For this research,the benchmark is executed with variousReynolds numbers up to the maximum Re=1000.All results then numerously compared to available sources in the literature.展开更多
Lid-driven mixed convection has been given immense importance due to its wide range of applications. A T-shaped cavity is introduced and pertinent parameters controlling mixed convection phenomenon are analyzed in thi...Lid-driven mixed convection has been given immense importance due to its wide range of applications. A T-shaped cavity is introduced and pertinent parameters controlling mixed convection phenomenon are analyzed in this paper. Water-Al2O3 nanofluid is considered inside the cavity to augment heat transfer rate. Galerkin weighted residual method of finite element analysis is applied for the numerical simulations. Numerical solution is obtained for different solid volume fractions of nanofluid（？= 0- 0.15）, Grashof numbers（Gr= 0.1 ？5 000） and Reynolds numbers（Re= 0.31 ？1 000） in laminar flow regime. Special attention is given on the analysis of flow at the pure mixed convection regime. It is found that Grashof, Reynolds and Richardson numbers along with solid volume fraction of nanofluid have significant effect on heat transfer characteristics inside the cavity. Results are presented using streamline and isotherm contours along with related variation of average Nusselt numbers of the heated wall and average fluid temperature inside the cavity.展开更多
In this study,numerical analysis and optimization in a single and multiple walled carbon nanotube-water nanofluid filled lid driven cavity having an inner elliptic obstacle were performed by using finite element metho...In this study,numerical analysis and optimization in a single and multiple walled carbon nanotube-water nanofluid filled lid driven cavity having an inner elliptic obstacle were performed by using finite element method and COBYLA optimization solver.The top wall is moving with constant speed and vertical walls are kept at constant temperatures.An optimal size of the inner elliptic obstacle was determined by using an optimization study to maximize the average heat transfer along the hot wall of the cavity.Numerical simulation was performed by for various values of Richardson numbers(between 0.05 and 50)and various solid particle volume fraction(between 0 and 0.06)for single and multiple-walled carbon nanotubes water nanofluid.A larger obstacle(higher values of radii in the major and minor axis)with lower values of Richardson number results in higher heat transfer rates.The average Nusselt number versus solid particle volume fraction shows a linear trend and the discrepancy between the average Nusselt number for the cavity with the optimized obstacle and other obstacles becomes higher with higher particle volume fraction.The average heat transfer enhances significantly which is about 120.20%for single wall carbon nanotube-water nanofluid at solid volume fraction of 0.06 when compared to pure water.The discrepancy between the average Nusselt number for single and multiple walled carbon nanotubes becomes higher for higher values of Richardson number and solid particle volume fraction.A polynomial type correlation was proposed for the average Nusselt number along the hot wall which is fifth order for Richardson number an first order for nanoparticle volume fraction.展开更多
A lattice Boltzmann model combined with curvilinear coordinate is proposed for lid-driven cavity three-dimensional (3D) flows. For particle velocity distribution, the particle collision process is performed in physica...A lattice Boltzmann model combined with curvilinear coordinate is proposed for lid-driven cavity three-dimensional (3D) flows. For particle velocity distribution, the particle collision process is performed in physical domain, and the particle streaming process is carried out in the corresponding computational domain, which is transferred from the physical domain using interpolation method. For the interpolation calculation, a second-order upwind interpolation method is adopted on internal lattice nodes in flow fields while a second-order central interpolation algorithm is employed at neighbor-boundary lattice nodes. Then the above-mentioned model and algorithms are used to numerically simulate the 3D flows in the lid-driven cavity at Reynolds numbers of 100, 400 and 1000 on non-uniform meshes. Various vortices on the x-y, y-z and x-z symmetrical planes are successfully predicted, and their changes in position with the Reynolds number increasing are obtained. The velocity profiles of u component along the vertical centerline and w component along the horizontal centerline are both in good agreement with the data in literature and the calculated results on uniform meshes. Besides, the velocity vector distributions on various cross sections in lid-driven cavity predicted on non-uniform meshes are compared with those simulated on uniform meshes and those in the literature. All the comparisons and validations show that the 3D lattice Boltzmann model and all the numerical algorithms on non-uniform meshes are accurate and reliable to predict effectively flow fields.展开更多
基金Project supported by the National Natural Science Foundation of China(No.12250410244)the Jiangsu Funding Program for Excellent Postdoctoral Talent of China(No.2023ZB884)+2 种基金the Foreign Expert Project funding of China(No.WGXZ2023017L)the Shuang-Chuang(SC)Doctor Program of Jiangsu Provincethe Longshan Scholar Program of Nanjing University of Information Science&Technology。
文摘The unsteady magnetohydrodynamical(MHD)free convection flow of an incompressible,electrically conducting hybrid nanofluid within a vertical cylindrical geometry is investigated,incorporating the effects of thermal radiation,viscous dissipation,and internal heat generation.The system is subjected to a time-periodic boundary temperature condition.The Laplace and finite Hankel transforms are used to derive the exact solutions for the velocity and temperature distributions.The effects of various key physical parameters,including the Richardson number,the Eckert number,the radiation parameter,the heat source parameter,and the nanoparticle volume fraction,are considered.The numerical results reveal that increasing the volume fraction significantly enhances the thermal conductivity and temperature,while the magnetic field intensity and viscous dissipation strongly influence the fluid motion and heat transport.Additionally,the pulsating boundary conditions produce distinct oscillatory behaviors in both the velocity and temperature fields.These findings provide important insights into optimizing the heat transfer performance in cylindrical systems such as electronic cooling modules and energy storage devices operating under dynamic thermal conditions.
文摘The motion of an elliptical rigid particle in a lid-driven cavity flow was numerically simulated using the immersed boundary lattice Boltzmann method(IB-LBM).The effects of the particle's initial orientation angle,initial position,aspect ratio,and size on the motion characteristics were investigated.The computational results indicate that the particle's motion undergoes two distinct stages:a starting stage that involves moving from the release position to a limit cycle,and a periodic stage that involves moving on the limit cycle.The initial orientation of the particle has a minimal impact on both stages of motion.In contrast,the time it takes for the particle to reach the limit cycle may vary depending on the release position.Furthermore,particles with a larger aspect ratio exhibit a greater maximum velocity magnitude;an increase in particle size causes the particle trajectory to contract more toward the center of the cavity,decreasing the centrifugal force experienced by the particle.
基金Project supported by the National Natural Science Foundation of China （Grant No 70271069）.
文摘By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number （Re） flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.
文摘Steady,laminar mixed convection inside a lid-driven square cavity filled with nanofluid is investigated numerically.We consider the case where the right and left walls are moving downwards and upwards respectively and maintained at different temperatures while the other two horizontal ones are kept adiabatic and impermeable.The set of nonlinear coupled governing mass,momentum,and energy equations are solved using an extensively validated and a highly accurate finite difference method of fourth-order.Comparisons with previously conducted investigations on special configurations are performed and show an excellent agreement.Meanwhile,attention is focused on the heat transfer enhancement when different nano-particles:Cu,Ag,Al2O3,TiO2 and Fe3O4 are incorporated separately in different base fluids such as:Water,Ethylene-glycol,Methanol and Kerosene oil.In this framework,the numerical results related to several mixtures are presented and concern flow pattern and heat transfer curves for various values of Richardson number[Ri=0.1,1 and 10].It turns out that the choice of the efficient binary mixture for an optimal heat transfer depends not only on the thermophysical properties of the nanofluids but also on the range of the Richardson number.Special attention is devoted to shedding light on the effect of the shape of the nanoparticles on the heat transfer in the case of Water-Ag nanofluid.It is concluded that the spherical shape is more suitable for a better heat transfer enhancement in comparison to the cylindrical ones.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11601013 and 91530325)。
文摘To date, there are very few studies on the second Hopf bifurcation in a driven square cavity, although there are intensive investigations focused on the first Hopf bifurcation in literature, due to the difficulties of theoretical analyses and numerical simulations. In this paper, we study the characteristics of the second Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme recently developed by us. We numerically identify the critical Reynolds number of the second Hopf bifurcation located in the interval of(11093.75, 11094.3604) by bisection. In addition, we find that there are two dominant frequencies in its spectral diagram when the flow is in the status of the second Hopf bifurcation, while only one dominant frequency is identified if the flow is in the first Hopf bifurcation via the Fourier analysis. More interestingly, the flow phase portrait of velocity components is found to make transition from a regular elliptical closed form for the first Hopf bifurcation to a non-elliptical closed form with self-intersection for the second Hopf bifurcation. Such characteristics disclose flow in a quasi-periodic state when the second Hopf bifurcation occurs.
基金the National Natural Science Foundation of China (Grants 41372301 and 51349011)the Preeminent Youth Talent Project of Southwest University of Science and Technology (Grant 13zx9109)
文摘A streamline upwind/Petrov-Galerkin （SUPG） finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.
基金Project supported by the National Natural Science Foundation of China(Grant No.12162001)the Natural Science Foundation of Ningxia(Grant No.2019AAC03129)the Construction Project of First-Class Disciplines in Ningxia Higher Education(Grant No.NXYLXK2017B09)。
文摘To date,there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity,due to the difficulties in theoretical analysis and numerical simulations.In this paper,we study the characteristics of the third Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme rectently developed by us.We numerically identify the critical Reynolds number of the third Hopf bifurcation located in the interval of(13944.7021,13946.5333)by the method of bisection.Through Fourier analysis,it is discovered that the flow becomes chaotic with a characteristic of period-doubling bifurcation when the Reynolds number is beyond the third bifurcation critical interval.Nonlinear time series analysis further ascertains the flow chaotic behaviors via the phase diagram,Kolmogorov entropy and maximal Lyapunov exponent.The phase diagram changes interestingly from a closed curve with self-intersection to an unclosed curve and the attractor eventually becomes strange when the flow becomes chaotic.
文摘This research work numerically analyzes 2D,steady state,mixed convective heat transfer for Newtonian fluids in lid driven square enclosure with centered triangular block（blockage—10%or 30%）maintained either at the constant wall temperature or constant heat flux thermal conditions.The fluid flow in the enclosure is initiated by top moving wall in＋x-direction,while all other walls are stationary.The top and bottom walls are thermally insulated.In particular,the governing field equations are solved for range of governing parameters such as,Reynolds number（1–1000）,Prandtl number（1–100）,and Grashof number
文摘Hydrodynamic mixed convection in a lid-driven hexagonal cavity with corner heater is numerically simulated in this paper by employing finite element method. The working fluid is assigned as air with a Prandtl num-ber of 0.71 throughout the simulation. The left and right walls of the hex-agonal cavity are kept thermally insulated and the lid moves top to bottom at a constant speed U0. The top left and right walls of the enclosure are maintained at cold temperature Tc. The bottom right wall is considered with a corner heater whereas the bottom remaining part is adiabatic and inside the cavity a square shape heated block Th. The focus of the work is to investigate the effect of Hartmann number, Richardson number, Grashof number and Reynolds number on the fluid flow and heat transfer characteristics inside the enclosure. A set of graphical results is presented in terms of streamlines, isotherms, local Nusselt number, velocity profiles, temperature profiles and average Nusselt numbers. The results reveal that heat transfer rate increases with increasing Richardson number and Hartmann number. It is also observed that, Hartmann number is a good control parameter for heat transfer in fluid flow in hexagonal cavity.
基金Supported by Science and Technology Development Planning of Shandong Province,P.R.China(2016GGX104018)
文摘In the case of lid-driven deep cavity flow, the effects of different resolutions of local grid refinement have been studied in the frame of multiple relaxation times (MRT) lattice Boltzmann method (LBM). In all the cases, the aspect ratio and Reynolds number are set as 1.5 and 3.200, respectively. First, the applied method is validated by comparing it with two reported works, with which agreements are reached. Then, six separate degrees of local grid refinement at the upper left corner, i.e. purely coarse grid, including 1/64, 1/32, 1/16, 1/8, 1/4 refinements of the lattice number in the width direction have been studied in detail. The results give the following indications:① The refinement degrees lower than 1/8 produce similar results;② For single corner refinement, 1/4 refinement is adequate for clearing the noises in the singularity zone to a large extent;③ New noise around the interface between coarse and fine zones are introduced by local grid refinement. Finally, refinement of entire subzone neighboring the lid is examined to avoid introducing new noises and it has been found effective.
文摘This study aims to analyze mixed convection in a square cavity with two moving vertical walls by finite volume method.The cavity filled with Non-Newtonian fluid of Bingham model is heated from below and cooled by the other walls.This study has been conducted for certain parameters of Reynolds number(Re=1-100),Richardson number(Ri=1-20),Prandtl number(Pr=1-500),and Bingham number has been studied from 0 to 10.The results indicate that the increase in yield stress drops the heat transfer and the flow become flatter,while increasing Reynolds number augments it.The convective transport is dominant when increasing Richardson number which leads to enhance heat transfer in the cavity for both Newtonian and Non-Newtonian fluid.A correlation of Nusselt number is given in function of different parameters.
基金The financial supports from the National Natural Science Foundation of China(22178014)。
文摘The concentration and velocity fields of two refractive index matched miscible shear-thinning fluids in a lid-driven cavity were investigated by using planar laser-induced fluorescence and particle image velocimetry,as well by computational fluid dynamics.Quantitative analyses show that the results obtained by flow simulations with the species transport model are in good agreement with the experimental results.The effects of different parameters were studied by using the intensity of segregation.For two fluids with the same rheological parameters,the relative amounts of liquids H_(1)/H and the power-law index n dominate the mixing process while the Reynolds number Re plays a marginal role.As for two fluids with density difference,buoyancy has significant influence on the mixing process.The dimensionless group Ar/Re(redefined such as to include shear thinning behavior)is proposed for assessing the effect of buoyancy and rheological properties on the mixing of miscible shear-thinning fluids.
文摘Numerical simulation of the bifurcation of Bingham fluid streamline topologies in rectangular double-lid-driven cavity, with varying aspect (height to width) ratio A, is presented. The lids on the top and bottom move at the same speed but in opposite directions so that symmetric flow patterns are generated. Similar to the Newtonian case, bifurcations occur as the aspect ratio decreases. Special to Bingham fluids, the non-Newtonian indicator, Bingham number B, also governs the bifurcation besides the bifurcation parameter A.
文摘An investigation on flow and heat transfer due to mixed convection, in a lid-driven rectangular cavity filled with Cu- water nanofluids and submitted to uniform heat flux along with its vertical short sides, has been conducted numerically by solving the full governing equations with the finite volume method and the SIMPLER algorithm. In the case of a slender enclosure, these equations are considerably reduced by using the parallel flow concept. Solutions, for the flow and temperature fields, and the heat transfer rate, have been obtained depending on the governing parameters, which are the Reynolds, the Richardson numbers and the solid volume fraction of nanoparticles. A perfect agreement has been found between the results of the two approaches for a wide range of the abovementioned parameters. It has been shown that at low and high Richardson numbers, the convection is ensured by lid and buoyancy-driven effects, respectively, whereas between these extremes, both mechanisms compete. Moreover, the addition of Cu-nanoparticles, into the pure water, has been seen enhancing and degrading heat transfer by lid and buoyancy-driven effects, respectively.
文摘This study is devoted to the computational fluid dynamics (CFD) modeling of steady laminar mixed convection flow and heat transfer in lid driven cavity (10 ≤ Re ≤ 1000). The ratio of the height to the width of the cavity is ranged over H/L = 0.5 to 1.5. The governing equations are solved using commercial finite volume package FLUENT to visualize the nature of the flow and estimate the heat transfer inside the cavity for different aspect ratio. The simulation results are presented in terms of average Nusselt number of the hot wall, velocity profile, and temperature contours. It was found that the average Nusselt number inside the cavity is strongly governed by the aspect ratio as well as the Reynolds number. A parametric study is conducted to demonstrate the effect of aspect ratio on the flow and heat transfer characteristics. It is found that heat transfer enhancement was obtained by decreasing the aspect ratio and/or increasing the Reynolds number.
基金This work was supported by the VLIR-UOS TEAM Project,VN2017TEA454A 103,‘An innovative solution to protect Vietnamese coastal riverbanks from floods and erosion’funded by the Flemish Government.
文摘This paper investigates a polygonal finite element(PFE)to solve a two-dimensional(2D)incompressible steady fluid problem in a cavity square.It is a well-known standard benchmark(i.e.,lid-driven cavity flow)-to evaluate the numerical methods in solving fluid problems controlled by the Navier-Stokes(N-S)equation system.The approximation solutions provided in this research are based on our developed equal-order mixed PFE,called Pe1Pe1.It is an exciting development based on constructing the mixed scheme method of two equal-order discretisation spaces for both fluid pressure and velocity fields of flows and our proposed stabilisation technique.In this research,to handle the nonlinear problem of N-S,the Picard iteration scheme is applied.Our proposed method’s performance and convergence are validated by several simulations coded by commercial software,i.e.,MATLAB.For this research,the benchmark is executed with variousReynolds numbers up to the maximum Re=1000.All results then numerously compared to available sources in the literature.
基金the support provided by the Department of Mechanical Engineering, Bangladesh University of Engineering and Technology (BUET), Dhaka, Bangladesh during this research work
文摘Lid-driven mixed convection has been given immense importance due to its wide range of applications. A T-shaped cavity is introduced and pertinent parameters controlling mixed convection phenomenon are analyzed in this paper. Water-Al2O3 nanofluid is considered inside the cavity to augment heat transfer rate. Galerkin weighted residual method of finite element analysis is applied for the numerical simulations. Numerical solution is obtained for different solid volume fractions of nanofluid（？= 0- 0.15）, Grashof numbers（Gr= 0.1 ？5 000） and Reynolds numbers（Re= 0.31 ？1 000） in laminar flow regime. Special attention is given on the analysis of flow at the pure mixed convection regime. It is found that Grashof, Reynolds and Richardson numbers along with solid volume fraction of nanofluid have significant effect on heat transfer characteristics inside the cavity. Results are presented using streamline and isotherm contours along with related variation of average Nusselt numbers of the heated wall and average fluid temperature inside the cavity.
文摘In this study,numerical analysis and optimization in a single and multiple walled carbon nanotube-water nanofluid filled lid driven cavity having an inner elliptic obstacle were performed by using finite element method and COBYLA optimization solver.The top wall is moving with constant speed and vertical walls are kept at constant temperatures.An optimal size of the inner elliptic obstacle was determined by using an optimization study to maximize the average heat transfer along the hot wall of the cavity.Numerical simulation was performed by for various values of Richardson numbers(between 0.05 and 50)and various solid particle volume fraction(between 0 and 0.06)for single and multiple-walled carbon nanotubes water nanofluid.A larger obstacle(higher values of radii in the major and minor axis)with lower values of Richardson number results in higher heat transfer rates.The average Nusselt number versus solid particle volume fraction shows a linear trend and the discrepancy between the average Nusselt number for the cavity with the optimized obstacle and other obstacles becomes higher with higher particle volume fraction.The average heat transfer enhances significantly which is about 120.20%for single wall carbon nanotube-water nanofluid at solid volume fraction of 0.06 when compared to pure water.The discrepancy between the average Nusselt number for single and multiple walled carbon nanotubes becomes higher for higher values of Richardson number and solid particle volume fraction.A polynomial type correlation was proposed for the average Nusselt number along the hot wall which is fifth order for Richardson number an first order for nanoparticle volume fraction.
基金supported by the National Natural Science Foundation of China (Grant Nos. 51179192, 50779069, 51139007)the Program for New Century Excellent Talents in University (NCET) (Grant No. NETC-10-0784)+1 种基金the National Hi-Tech Research and Development Program of China ("863" Project) (Grant No. 2011AA100505)the Chinese Universities Scientific Fund (Grant No. 2013RC045)
文摘A lattice Boltzmann model combined with curvilinear coordinate is proposed for lid-driven cavity three-dimensional (3D) flows. For particle velocity distribution, the particle collision process is performed in physical domain, and the particle streaming process is carried out in the corresponding computational domain, which is transferred from the physical domain using interpolation method. For the interpolation calculation, a second-order upwind interpolation method is adopted on internal lattice nodes in flow fields while a second-order central interpolation algorithm is employed at neighbor-boundary lattice nodes. Then the above-mentioned model and algorithms are used to numerically simulate the 3D flows in the lid-driven cavity at Reynolds numbers of 100, 400 and 1000 on non-uniform meshes. Various vortices on the x-y, y-z and x-z symmetrical planes are successfully predicted, and their changes in position with the Reynolds number increasing are obtained. The velocity profiles of u component along the vertical centerline and w component along the horizontal centerline are both in good agreement with the data in literature and the calculated results on uniform meshes. Besides, the velocity vector distributions on various cross sections in lid-driven cavity predicted on non-uniform meshes are compared with those simulated on uniform meshes and those in the literature. All the comparisons and validations show that the 3D lattice Boltzmann model and all the numerical algorithms on non-uniform meshes are accurate and reliable to predict effectively flow fields.
