For the safety of railroad operations,sand barriers are utilized to mitigate wind-sand disaster effects.These disasters,characterized by multi-directional wind patterns,result in diverse angles among the barriers.In t...For the safety of railroad operations,sand barriers are utilized to mitigate wind-sand disaster effects.These disasters,characterized by multi-directional wind patterns,result in diverse angles among the barriers.In this study,using numerical simulations,we examined the behavior of High Density Polyethylene(HDPE)sheet sand barriers under different wind angles,focusing on flow field distribution,windproof efficiency,and sedimentation erosion dynamics.This study discovered that at a steady wind speed,airflow velocity varies as the angle between the airflow and the HDPE barrier changes.Specifically,a 90°angle results in the widest low-speed airflow area on the barrier’s downwind side.If the airflow is not perpendicular to the barrier,it prompts a lateral airflow movement which decreases as the angle expands.The windproof efficiency correlates directly with this angle but inversely with the wind’s speed.Notably,with a wind angle of 90°,wind speed drops by 81%.The minimum wind speed is found at 5.1H(the sand barrier height)on the barrier’s downwind side.As the angle grows,the barrier’s windproof efficiency improves,extending its protective reach.Sedimentation is most prominent on the barrier’s downwind side,as the wind angle shifts from 30°to 90°,the sand sedimentation area on the barrier’s downwind side enlarges by 14.8H.As the angle grows,sedimentation intensifies,eventually overtakes the forward erosion and enlarges the sedimentation area.展开更多
Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pr...Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.展开更多
Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed bounda...Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.展开更多
This paper reviews the authors' recent studies on compressible turbulence by using direct numerical simulation (DNS),including DNS of isotropic(decaying) turbulence, turbulent mixing-layer,turbulent boundary-laye...This paper reviews the authors' recent studies on compressible turbulence by using direct numerical simulation (DNS),including DNS of isotropic(decaying) turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.展开更多
To improve the mechanical properties of the parts fabricated by Laser Direct Metal Shaping(LDMS),it is of great significance to understand the distribution regularities of transient temperature field during LDMS proce...To improve the mechanical properties of the parts fabricated by Laser Direct Metal Shaping(LDMS),it is of great significance to understand the distribution regularities of transient temperature field during LDMS process.Based on the“el- ement birth and death”technique of finite element method,a three-dimensional multi-track and multi-layer model for the transient temperature field analysis of LDMS is developed by ANSYS Parametric Design Language(APDL)for the first time.In the fab- ricated modal,X-direction parallel reciprocating scanning paths is introduced.Using the same process parameters,the simulation results show good agreement with the microstructure features of samples which fabricated by LDMS.展开更多
Direct numerical simulation of a spatially developing turbulent boundary layer over a compliant wall with anisotropic wall material properties is performed. The Reynolds number varies from 300 to approximately 860 alo...Direct numerical simulation of a spatially developing turbulent boundary layer over a compliant wall with anisotropic wall material properties is performed. The Reynolds number varies from 300 to approximately 860 along the streamwise direction, based on the external flow velocity and the momentum thickness. Eight typical cases are selected for numerical investigation under the guidance of the monoharmonic analysis. The instantaneous flow fields exhibit the traveling wavy motion of the compliant wall, and the frequency-wavenumber power spectrum of wall pressure fluctuation is computed to quantify the mutual influence of the wall compliance and the turbulent flow at different wave numbers. It is shown that the Reynolds shear stress and the pressure fluctuation are generally enhanced by the wall compliance with the parameters considered in the present study. A dynamical decomposition of the skin-friction coefficient is derived, and a new term (CW) appears due to the wall-induced Reynolds shear stress. The influence of the anisotropic compliant wall motion on the turbulent boundary layer through the wall-induced negative Reynolds shear stress is discussed. To elucidate the underlying mechanism, the budget analysis of the Reynolds stresses transportation is further carried out. The impact of the wall compliance on the turbulent flow is disclosed by examining the variations of the diffusion and velocity-pressure correlation terms. It is shown that increase of the Reynolds stresses inside the flow domain is caused by enhancement of the velocity-pressure correlation term, possibly through the long-range influence of the wall compliance on the pressure field, rather than diffusion of the wall-induced Reynolds shear stress into the fluid flow.展开更多
Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics h...Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.展开更多
Direct numerical simulation of decaying homogeneous isotropic turbulence (DHIT) of a polymer solution is performed. In order to understand the polymer effect on turbulence or additive-turbulence interaction, we dire...Direct numerical simulation of decaying homogeneous isotropic turbulence (DHIT) of a polymer solution is performed. In order to understand the polymer effect on turbulence or additive-turbulence interaction, we directly investigate the influence of polymers on velocity gradient tensor including vorticity and strain. By visualizing vortex tubes and sheets, we observe a remarkable inhibition of vortex structures in an intermediate-scale field and a small-scale field but not for a large scale field in DHIT with polymers. The geometric study indicates a strong relevance among the vorticity vector, rate-of-strain tensor, and polymer conformation tensor. Joint probability density functions show that the polymer effect can increase "strain generation resistance" and "vorticity generation resistance", i.e., inhibit the generation of vortex sheets and tubes, ultimately leading to turbulence inhibition effects.展开更多
A numerical method was developed to directly simulate the compressible, particle-laden turbulent jets.The fourth order compact finite difference schemes were used to discretize the space derivatives. The Lagrangian me...A numerical method was developed to directly simulate the compressible, particle-laden turbulent jets.The fourth order compact finite difference schemes were used to discretize the space derivatives. The Lagrangian method was adopted to simulate the particle motion based on one-way coupling. It is found that the turbulent intensity profiles attain self-similar status in the jet downstream regions. At the Stokes number of 1, particles are concentrated largely in the outer boundaries of the large-scale vortex structures with the most uneven distribution and the widest dispersion in the lateral direction. Particles at the much smaller Stokes numbers are distributed evenly in the flow field, and the lateral dispersion is also considerable. Distribution of particles at much larger Stokes numbers is more uniform and the lateral dispersion becomes small. In addition, the inflow conditions have different effects on the particle dispersion. The direct numerical simulation (DNS) results accord with the previous experiments and numerical studies.展开更多
Our previous experimental studies have confirmed that viscoelastic-fluid-based nanofluid(VFBN) prepared by suspending nanoparticles in a viscoelastic base fluid(VBF, behaves drag reduction at turbulent flow state) can...Our previous experimental studies have confirmed that viscoelastic-fluid-based nanofluid(VFBN) prepared by suspending nanoparticles in a viscoelastic base fluid(VBF, behaves drag reduction at turbulent flow state) can reduce turbulent flow resistance as compared with water and enhance heat transfer as compared with VBF. Direct numerical simulation(DNS) is performed in this study to explore the mechanisms of heat transfer enhancement(HTE) and flow drag reduction(DR) for the VFBN turbulent flow. The Giesekus model is used as the constitutive equation for VFBN. Our previously proposed thermal dispersion model is adopted to take into account the thermal dispersion effects of nanoparticles in the VFBN turbulent flow. The DNS results show similar behaviors for flow resistance and heat transfer to those obtained in our previous experiments. Detailed analyses are conducted for the turbulent velocity, temperature, and conformation fields obtained by DNSs for different fluid cases, and for the friction factor with viscous, turbulent, and elastic contributions and heat transfer rate with conductive, turbulent and thermal dispersion contributions of nanoparticles, respectively. The mechanisms of HTE and DR of VFBN turbulent flows are then discussed. Based on analogy theory, the ratios of Chilton–Colburn factor to friction factor for different fluid flow cases are investigated, which from another aspect show the significant enhancement in heat transfer performance for some cases of water-based nanofluid and VFBN turbulent flows.展开更多
Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS...Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS) methods for different Reynolds numbers. A formula is derived to express the relation between fluctuating velocities and the friction drag coefficient. With the application of electromagnetic force, the in-depth relations among the fluctuating velocities near the wall, Reynolds stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The fluctuating velocities along the streamwise and normal directions are suppressed significantly,while the fluctuating velocity along the spanwise direction is enhanced dramatically due to the spanwise electromagnetic force. However, the values of Reynolds stress depend on the fluctuating velocities along the streamwise and normal directions rather than that along the spanwise direction. Therefore, the significant effect of drag reduction is obtained. Moreover,the maximum drag reduction is weakened due to the decay of control effect for fluctuating velocities as the Reynolds number increases.展开更多
Direct numerical simulation(DNS) of gas–solid flow at high resolution has been carried out by coupling the lattice Boltzmann method(LBM) for gas flow and the discrete element method(DEM) for solid particles. However,...Direct numerical simulation(DNS) of gas–solid flow at high resolution has been carried out by coupling the lattice Boltzmann method(LBM) for gas flow and the discrete element method(DEM) for solid particles. However,the body force periodic boundary condition(FPBC) commonly used to cut down the huge computational cost of such simulation has faced accuracy concerns. In this study, a novel two-region periodic boundary condition(TPBC) is presented to remedy this problem, with the flow driven in the region with body force and freely evolving in the other region. With simulation cases for simple circulating fluidized bed risers, the validity and advantages of TPBC are demonstrated with more reasonable heterogeneity of the particle distribution as compared to the corresponding case with FPBC.展开更多
Accurate descriptions of matrix diffusion across the fracture/matrix interface are critical to assessing contaminant migration in fractured media. The classical transfer probability method is only applicable for relat...Accurate descriptions of matrix diffusion across the fracture/matrix interface are critical to assessing contaminant migration in fractured media. The classical transfer probability method is only applicable for relatively large diffusion coefficients and small fracture spacings, due to an intrinsic assumption of an equilibrium concentration profile in the matrix blocks. Motivated and required by practical applications, we propose a direct numerical simulation (DNS) approach without any empirical assumptions. A three-step Lagrangian algorithm was developed and validated to directly track the particle dynamics across the fracture/matrix interface, where particle's diffusive displacement across the discontinuity is controlled by an analytical, one-side reflection probability. Numerical experiments show that the DNS approach is especially efficient for small diffusion coefficients and large fracture spacings, alleviating limitations of the classical modeling approach.展开更多
Solute diffusion controlled solidification model was applied to simulate the columnar to equiaxed transition (CET) during directional solidification of Ti-44Al alloy. The simulation results show that the solutal int...Solute diffusion controlled solidification model was applied to simulate the columnar to equiaxed transition (CET) during directional solidification of Ti-44Al alloy. The simulation results show that the solutal interactions from growing equiaxed grains play an important role on CET. The effects of the applied thermal gradient and pulling velocity, the equiaxed seed spacing and nucleation undercooling on the CET are investigated in the present simulation. The simulated results indicated that the columnar branch spacing depends not only on the thermal gradient and the pulling velocity, but also on number of the seeds. A spacing adjustment can occur through initiation of seeds that develop into new columnar grains. The dependence of the CET on the thermal gradient and pulling velocity, qualitatively agrees with the analytical CET model of Hunt,展开更多
Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data ...Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data (such as the distribution of fractures, low permeability zones, faults and major stress) and simulation experiences to design the grid direction of simulation model qualitatively. The research of the paper indicates that the key to determine the grid direction is to determine the principal permeability direction. Under the circumstances of few static materials, a new grid direction determination method has been developed by using field data (well location map and inter-well permeability) on the bases of Darcy’s law and tensor analysis theory. The grid direction of WZ11-7 Oilfield simulation model has been determined using four production wells and two production zones (L1 and L3) in WZ11-7-2 well group, the results are in conformity with the geological studied major stress. Therefore, this method can give insights into the numerical simulation study.展开更多
Direct numerical simulations(DNSs) of purely elastic turbulence in rectilinear shear flows in a three-dimensional(3D) parallel plate channel were carried out,by which numerical databases were established.Based on ...Direct numerical simulations(DNSs) of purely elastic turbulence in rectilinear shear flows in a three-dimensional(3D) parallel plate channel were carried out,by which numerical databases were established.Based on the numerical databases,the present paper analyzed the structural and statistical characteristics of the elastic turbulence including flow patterns,the wall effect on the turbulent kinetic energy spectrum,and the local relationship between the flow motion and the microstructures' behavior.Moreover,to address the underlying physical mechanism of elastic turbulence,its generation was presented in terms of the global energy budget.The results showed that the flow structures in elastic turbulence were 3D with spatial scales on the order of the geometrical characteristic length,and vortex tubes were more likely to be embedded in the regions where the polymers were strongly stretched.In addition,the patterns of microstructures' elongation behave like a filament.From the results of the turbulent kinetic energy budget,it was found that the continuous energy releasing from the polymers into the main flow was the main source of the generation and maintenance of the elastic turbulent status.展开更多
A kind of direct numerical simulation method suitable for supercritical carbon dioxide jet flow has been discussed in this paper. The form of dimensionless nonconservative compressible Navier-Stokes equations in a two...A kind of direct numerical simulation method suitable for supercritical carbon dioxide jet flow has been discussed in this paper. The form of dimensionless nonconservative compressible Navier-Stokes equations in a two-dimensional cartesian coordinate system is derived in detail. High accurate finite difference compact schemes based on non-uniform grid system are introduced to solve the equations. The simulation results of the three vortex pairing phenomenon of plane mixing layer and a compressible axisymmetric jet flow field show that the discussed numerical simulation method is feasible to calculate the supercritical carbon dioxide jet fluid. And it is found that the difficulties of splitting the convective terms in conservation Navier-Stokes equations, which are brought by the supercritical carbon dioxide fluid pressure state equation, can be avoided by solving the nonconservative compressible Navier-Stokes equations.展开更多
基金financially supported by the Natural Science Foundation of Gansu Province,China(22JR5RA050,20JR10RA231)the fellowship of the China Postdoctoral Science Foundation(2021M703466)the Basic Research Innovation Group Project of Gansu Province,China(21JR7RA347).
文摘For the safety of railroad operations,sand barriers are utilized to mitigate wind-sand disaster effects.These disasters,characterized by multi-directional wind patterns,result in diverse angles among the barriers.In this study,using numerical simulations,we examined the behavior of High Density Polyethylene(HDPE)sheet sand barriers under different wind angles,focusing on flow field distribution,windproof efficiency,and sedimentation erosion dynamics.This study discovered that at a steady wind speed,airflow velocity varies as the angle between the airflow and the HDPE barrier changes.Specifically,a 90°angle results in the widest low-speed airflow area on the barrier’s downwind side.If the airflow is not perpendicular to the barrier,it prompts a lateral airflow movement which decreases as the angle expands.The windproof efficiency correlates directly with this angle but inversely with the wind’s speed.Notably,with a wind angle of 90°,wind speed drops by 81%.The minimum wind speed is found at 5.1H(the sand barrier height)on the barrier’s downwind side.As the angle grows,the barrier’s windproof efficiency improves,extending its protective reach.Sedimentation is most prominent on the barrier’s downwind side,as the wind angle shifts from 30°to 90°,the sand sedimentation area on the barrier’s downwind side enlarges by 14.8H.As the angle grows,sedimentation intensifies,eventually overtakes the forward erosion and enlarges the sedimentation area.
基金supported by the National Natural Science Foundation of China(Project Nos.12272270,11972261).
文摘Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.
基金financially supported by the National Key R&D Program of China (Grant No.2022YFB2603000)the National Natural Science Foundation of China (Grant Nos.51779172 and 52179076)。
文摘Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.
基金supported by the National Basic Research Program of China(2009CB724100)the National Natural Science Foundation of China(10632050,10872205,11072248).
文摘This paper reviews the authors' recent studies on compressible turbulence by using direct numerical simulation (DNS),including DNS of isotropic(decaying) turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.
文摘To improve the mechanical properties of the parts fabricated by Laser Direct Metal Shaping(LDMS),it is of great significance to understand the distribution regularities of transient temperature field during LDMS process.Based on the“el- ement birth and death”technique of finite element method,a three-dimensional multi-track and multi-layer model for the transient temperature field analysis of LDMS is developed by ANSYS Parametric Design Language(APDL)for the first time.In the fab- ricated modal,X-direction parallel reciprocating scanning paths is introduced.Using the same process parameters,the simulation results show good agreement with the microstructure features of samples which fabricated by LDMS.
基金the National Natural Science Foundation of China (Grants 11772172 and 11490551).
文摘Direct numerical simulation of a spatially developing turbulent boundary layer over a compliant wall with anisotropic wall material properties is performed. The Reynolds number varies from 300 to approximately 860 along the streamwise direction, based on the external flow velocity and the momentum thickness. Eight typical cases are selected for numerical investigation under the guidance of the monoharmonic analysis. The instantaneous flow fields exhibit the traveling wavy motion of the compliant wall, and the frequency-wavenumber power spectrum of wall pressure fluctuation is computed to quantify the mutual influence of the wall compliance and the turbulent flow at different wave numbers. It is shown that the Reynolds shear stress and the pressure fluctuation are generally enhanced by the wall compliance with the parameters considered in the present study. A dynamical decomposition of the skin-friction coefficient is derived, and a new term (CW) appears due to the wall-induced Reynolds shear stress. The influence of the anisotropic compliant wall motion on the turbulent boundary layer through the wall-induced negative Reynolds shear stress is discussed. To elucidate the underlying mechanism, the budget analysis of the Reynolds stresses transportation is further carried out. The impact of the wall compliance on the turbulent flow is disclosed by examining the variations of the diffusion and velocity-pressure correlation terms. It is shown that increase of the Reynolds stresses inside the flow domain is caused by enhancement of the velocity-pressure correlation term, possibly through the long-range influence of the wall compliance on the pressure field, rather than diffusion of the wall-induced Reynolds shear stress into the fluid flow.
文摘Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.
基金supported by the National Natural Science Foundation of China (Grant No. 10872060)the Fundamental Research Funds for the Central Universities (Grant No. HIT.BRET2.2010008)
文摘Direct numerical simulation of decaying homogeneous isotropic turbulence (DHIT) of a polymer solution is performed. In order to understand the polymer effect on turbulence or additive-turbulence interaction, we directly investigate the influence of polymers on velocity gradient tensor including vorticity and strain. By visualizing vortex tubes and sheets, we observe a remarkable inhibition of vortex structures in an intermediate-scale field and a small-scale field but not for a large scale field in DHIT with polymers. The geometric study indicates a strong relevance among the vorticity vector, rate-of-strain tensor, and polymer conformation tensor. Joint probability density functions show that the polymer effect can increase "strain generation resistance" and "vorticity generation resistance", i.e., inhibit the generation of vortex sheets and tubes, ultimately leading to turbulence inhibition effects.
基金Supported by the Natural Science Foundation of Zhejiang Province (No. 502047 and No. M503094)National Basic Research Program of China (No. 2003CB214500).
文摘A numerical method was developed to directly simulate the compressible, particle-laden turbulent jets.The fourth order compact finite difference schemes were used to discretize the space derivatives. The Lagrangian method was adopted to simulate the particle motion based on one-way coupling. It is found that the turbulent intensity profiles attain self-similar status in the jet downstream regions. At the Stokes number of 1, particles are concentrated largely in the outer boundaries of the large-scale vortex structures with the most uneven distribution and the widest dispersion in the lateral direction. Particles at the much smaller Stokes numbers are distributed evenly in the flow field, and the lateral dispersion is also considerable. Distribution of particles at much larger Stokes numbers is more uniform and the lateral dispersion becomes small. In addition, the inflow conditions have different effects on the particle dispersion. The direct numerical simulation (DNS) results accord with the previous experiments and numerical studies.
基金supported by the National Natural Science Foundation of China(Grant No.51276046)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20112302110020)+1 种基金the China Postdoctoral Science Foundation(Grant No.2014M561037)the President Fund of University of Chinese Academy of Sciences,China(Grant No.Y3510213N00)
文摘Our previous experimental studies have confirmed that viscoelastic-fluid-based nanofluid(VFBN) prepared by suspending nanoparticles in a viscoelastic base fluid(VBF, behaves drag reduction at turbulent flow state) can reduce turbulent flow resistance as compared with water and enhance heat transfer as compared with VBF. Direct numerical simulation(DNS) is performed in this study to explore the mechanisms of heat transfer enhancement(HTE) and flow drag reduction(DR) for the VFBN turbulent flow. The Giesekus model is used as the constitutive equation for VFBN. Our previously proposed thermal dispersion model is adopted to take into account the thermal dispersion effects of nanoparticles in the VFBN turbulent flow. The DNS results show similar behaviors for flow resistance and heat transfer to those obtained in our previous experiments. Detailed analyses are conducted for the turbulent velocity, temperature, and conformation fields obtained by DNSs for different fluid cases, and for the friction factor with viscous, turbulent, and elastic contributions and heat transfer rate with conductive, turbulent and thermal dispersion contributions of nanoparticles, respectively. The mechanisms of HTE and DR of VFBN turbulent flows are then discussed. Based on analogy theory, the ratios of Chilton–Colburn factor to friction factor for different fluid flow cases are investigated, which from another aspect show the significant enhancement in heat transfer performance for some cases of water-based nanofluid and VFBN turbulent flows.
基金Project supported by the National Natural Science Foundation of China(Grant No.11672135)a Foundation for the Author of National Excellent Doctoral Dissertation of China(Grant No.201461)
文摘Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS) methods for different Reynolds numbers. A formula is derived to express the relation between fluctuating velocities and the friction drag coefficient. With the application of electromagnetic force, the in-depth relations among the fluctuating velocities near the wall, Reynolds stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The fluctuating velocities along the streamwise and normal directions are suppressed significantly,while the fluctuating velocity along the spanwise direction is enhanced dramatically due to the spanwise electromagnetic force. However, the values of Reynolds stress depend on the fluctuating velocities along the streamwise and normal directions rather than that along the spanwise direction. Therefore, the significant effect of drag reduction is obtained. Moreover,the maximum drag reduction is weakened due to the decay of control effect for fluctuating velocities as the Reynolds number increases.
基金Supported by the National Natural Science Foundation of China(21821005,91834303)Science Challenge Project(TZ2016001)+1 种基金the Key Research Program of Frontier Science of the Chinese Academy of Sciences(QYZDJ-SSW-JSC029)the Strategic Priority Research Program of the CAS(XDA21030700).
文摘Direct numerical simulation(DNS) of gas–solid flow at high resolution has been carried out by coupling the lattice Boltzmann method(LBM) for gas flow and the discrete element method(DEM) for solid particles. However,the body force periodic boundary condition(FPBC) commonly used to cut down the huge computational cost of such simulation has faced accuracy concerns. In this study, a novel two-region periodic boundary condition(TPBC) is presented to remedy this problem, with the flow driven in the region with body force and freely evolving in the other region. With simulation cases for simple circulating fluidized bed risers, the validity and advantages of TPBC are demonstrated with more reasonable heterogeneity of the particle distribution as compared to the corresponding case with FPBC.
基金supported by the United States Department of Energythe Desert Research Institute IR&D Funds
文摘Accurate descriptions of matrix diffusion across the fracture/matrix interface are critical to assessing contaminant migration in fractured media. The classical transfer probability method is only applicable for relatively large diffusion coefficients and small fracture spacings, due to an intrinsic assumption of an equilibrium concentration profile in the matrix blocks. Motivated and required by practical applications, we propose a direct numerical simulation (DNS) approach without any empirical assumptions. A three-step Lagrangian algorithm was developed and validated to directly track the particle dynamics across the fracture/matrix interface, where particle's diffusive displacement across the discontinuity is controlled by an analytical, one-side reflection probability. Numerical experiments show that the DNS approach is especially efficient for small diffusion coefficients and large fracture spacings, alleviating limitations of the classical modeling approach.
文摘Solute diffusion controlled solidification model was applied to simulate the columnar to equiaxed transition (CET) during directional solidification of Ti-44Al alloy. The simulation results show that the solutal interactions from growing equiaxed grains play an important role on CET. The effects of the applied thermal gradient and pulling velocity, the equiaxed seed spacing and nucleation undercooling on the CET are investigated in the present simulation. The simulated results indicated that the columnar branch spacing depends not only on the thermal gradient and the pulling velocity, but also on number of the seeds. A spacing adjustment can occur through initiation of seeds that develop into new columnar grains. The dependence of the CET on the thermal gradient and pulling velocity, qualitatively agrees with the analytical CET model of Hunt,
文摘Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data (such as the distribution of fractures, low permeability zones, faults and major stress) and simulation experiences to design the grid direction of simulation model qualitatively. The research of the paper indicates that the key to determine the grid direction is to determine the principal permeability direction. Under the circumstances of few static materials, a new grid direction determination method has been developed by using field data (well location map and inter-well permeability) on the bases of Darcy’s law and tensor analysis theory. The grid direction of WZ11-7 Oilfield simulation model has been determined using four production wells and two production zones (L1 and L3) in WZ11-7-2 well group, the results are in conformity with the geological studied major stress. Therefore, this method can give insights into the numerical simulation study.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51276046 and 51506037)the Foundation for Innovative Research Groups of the National Natural Science Foundation of China(Grant No.51421063)+2 种基金the China Postdoctoral Science Foundation(Grant No.2016M591526)the Heilongjiang Postdoctoral Fund,China(Grant No.LBH-Z15063)the China Postdoctoral International Exchange Program
文摘Direct numerical simulations(DNSs) of purely elastic turbulence in rectilinear shear flows in a three-dimensional(3D) parallel plate channel were carried out,by which numerical databases were established.Based on the numerical databases,the present paper analyzed the structural and statistical characteristics of the elastic turbulence including flow patterns,the wall effect on the turbulent kinetic energy spectrum,and the local relationship between the flow motion and the microstructures' behavior.Moreover,to address the underlying physical mechanism of elastic turbulence,its generation was presented in terms of the global energy budget.The results showed that the flow structures in elastic turbulence were 3D with spatial scales on the order of the geometrical characteristic length,and vortex tubes were more likely to be embedded in the regions where the polymers were strongly stretched.In addition,the patterns of microstructures' elongation behave like a filament.From the results of the turbulent kinetic energy budget,it was found that the continuous energy releasing from the polymers into the main flow was the main source of the generation and maintenance of the elastic turbulent status.
文摘A kind of direct numerical simulation method suitable for supercritical carbon dioxide jet flow has been discussed in this paper. The form of dimensionless nonconservative compressible Navier-Stokes equations in a two-dimensional cartesian coordinate system is derived in detail. High accurate finite difference compact schemes based on non-uniform grid system are introduced to solve the equations. The simulation results of the three vortex pairing phenomenon of plane mixing layer and a compressible axisymmetric jet flow field show that the discussed numerical simulation method is feasible to calculate the supercritical carbon dioxide jet fluid. And it is found that the difficulties of splitting the convective terms in conservation Navier-Stokes equations, which are brought by the supercritical carbon dioxide fluid pressure state equation, can be avoided by solving the nonconservative compressible Navier-Stokes equations.