The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investig...The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investigated. Results show that the disturbance velocity presents the sine vibration features with the change of time, and the vibration period is identical to the vibration of the local wall. The disturbance velocity presents the fluctuation feature downstream, and the streamwise wave length approximates to the results from the Orr-Sommerfeld equation （OSE）. The growth rate from direct numerical simulation（DNS） is a little greater than that from OSE, and their trends are almost consistent. Under the condition of Re= 2 800, the disturbance amplitude gradually grows in the given computational region with the period T=30. However, it firstly increases and then decreases with the period T= 20. The disturbance harmonic of the former is obviously larger than that of the latter. The maximum streamwise and vertical disturbance velocities from DNS do not fully coincide with those from OSE at the vicinity of the local vibration wall, but coincide well with the former when they travel downstream. The 2-D disturbance induced by the local micro-vibration represents the form of Tollmien-Schlichting （T-S） wave on the boundary layer.展开更多
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.展开更多
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.展开更多
We report the results of protein folding （219M, C34, N36, 2KES, 2KHK） by the method of accelerated molecular dynamics （aMD） at room temperature with the implicit solvent model. Starting from the linear structures,...We report the results of protein folding （219M, C34, N36, 2KES, 2KHK） by the method of accelerated molecular dynamics （aMD） at room temperature with the implicit solvent model. Starting from the linear structures, these proteins successfully fold to the native structure in a lO0-ns aMD simulation. In contrast, they are failed under the traditional MD simulation in the same simulation time. Then we find that the lowest root mean square deviations of helix structures from the native structures are 0.36 A, 0.63 A, 0.52 A, 1.1 A and 0.78 A. What is more, native contacts, cluster and free energy analyses show that the results of the aMD method are in accordance with the experiment very well. All analyses show that the aMD can accelerate the simulation process, thus we may apply it to the field of computer aided drug designs.展开更多
In this paper,we present a direct numerical simulation（DNS） of elastic turbulence of viscoelastic fluid at vanishingly low Reynolds number（Re = 1） in a three-dimensional straight channel flow for the first time,us...In this paper,we present a direct numerical simulation（DNS） of elastic turbulence of viscoelastic fluid at vanishingly low Reynolds number（Re = 1） in a three-dimensional straight channel flow for the first time,using the Giesekus constitutive model for the fluid.In order to generate and maintain the turbulent fluid motion in the straight channel,a sinusoidal force term is added to the momentum equation,and then the elastic turbulence is numerically realized with an initialized chaotic velocity field and a stretched conformation field.Statistical and structural characteristics of the elastic turbulence therein are analyzed based on the detailed information obtained from the DNS.The fluid mixing enhancement effect of elastic turbulence is also demonstrated for the potential applications of this phenomenon.展开更多
The direct implicit particle-in-cell is a powerful kinetic method for researching plasma characteristics.However,it is time-consuming to obtain the future electromagnetic field in such a method since the field equatio...The direct implicit particle-in-cell is a powerful kinetic method for researching plasma characteristics.However,it is time-consuming to obtain the future electromagnetic field in such a method since the field equations contain time-dependent matrix coefficients.In this work,we propose to explicitly push particles and obtain the future electromagnetic field based on the information about the particles in the future.The new method retains the form of implicit particle pusher,but the future field is obtained by solving the traditional explicit equation.Several numerical experiments,including the motion of charged particle in electromagnetic field,plasma sheath,and free diffusion of plasma into vacuum,are implemented to evaluate the performance of the method.The results demonstrate that the proposed method can suppress finite-grid-instability resulting from the coarse spatial resolution in electron Debye length through the strong damping of high-frequency plasma oscillation,while accurately describe low-frequency plasma phenomena,with the price of losing the numerical stability at large time-step.We believe that this work is helpful for people to research the bounded plasma by using particle-in-cell simulations.展开更多
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.展开更多
The swimming performance of rod-shaped microswimmers in a channel was numerically investigated using the two-dimensional lattice Boltzmann method(LBM).We considered variable-length squirmer rods,assembled from circula...The swimming performance of rod-shaped microswimmers in a channel was numerically investigated using the two-dimensional lattice Boltzmann method(LBM).We considered variable-length squirmer rods,assembled from circular squirmer models with self-propulsion mechanisms,and analyzed the effects of the Reynolds number(Re),aspect ratio(ε),squirmer-type factor(β)and blockage ratio(κ)on swimming efficiency(η)and power expenditure(P).The results show no significant difference in power expenditure between pushers(microswimmers propelled from the tail)and pullers(microswimmers propelled from the head)at the low Reynolds numbers adopted in this study.However,the swimming efficiency of pushers surpasses that of pullers.Moreover,as the degree of channel blockage increases(i.e.,κincreases),the squirmer rod consumes more energy while swimming,and its swimming efficiency also increases,which is clearly reflected whenε≤3.Notably,squirmer rods with a larger aspect ratioεand aβvalue approaching 0 can achieve high swimming efficiency with lower power expenditure.The advantages of self-propelled microswimmers are manifested whenε>4 andβ=±1,where the squirmer rod consumes less energy than a passive rod driven by an external field.These findings underscore the potential for designing more efficient microswimmers by carefully considering the interactions between the microswimmer geometry,propulsion mechanism and fluid dynamic environment.展开更多
Turbulent spots play a key role in the formation of the turbulence and the transition. The generation and evolution of turbulent spots using the wall impulse model in the plane Couette flow are studied by direct numer...Turbulent spots play a key role in the formation of the turbulence and the transition. The generation and evolution of turbulent spots using the wall impulse model in the plane Couette flow are studied by direct numerical simulation of Navier-Stokes equations. A group of three-dimensional coupling compact difference schemes with high accuracy and high resolution is used in the numerical calculation. The important characteristics of turbulent spots based on the results of examples are analyzed, including the formation of random pulse, the generation of Reynolds stress, the growth of disturbance amplitude, and the continuous change of spot shape, especially the complex evolution process of the streamwise vortices. Computational results confirm that basic properties of turbulent spots in the laminar flow are similar to those in the turbulent flow.展开更多
Numerical simulations are presented about the effects of gas rarefaction on hypersonic flow field.Due to the extremely difficult experiment,limited wind-tunnel conditions and high cost,most problems in rarefied flow r...Numerical simulations are presented about the effects of gas rarefaction on hypersonic flow field.Due to the extremely difficult experiment,limited wind-tunnel conditions and high cost,most problems in rarefied flow regime are investigated through numerical methods,in which the direct simulation Monte-Carlo(DSMC)method is widely adopted.And the unstructured DSMC method is employed here.Flows around a vertical plate at a given velocity 7 500 m/s are simulated.For gas rarefaction is judged by the free-stream Knudsen number(Kn),two vital factors are considered:molecular number density and the plate′s length.Cases in which Kn varies from 0.035 to13.36 are simulated.Flow characters in the whole rarefied regime are described,and flow-field structure affected by Knis analyzed.Then,the dimensionless position D*of a certain velocity in the stagnation line is chosen as the marker of flow field to measure its variation.Through flow-field tracing and least-square numerical method analyzing,it is proved that hypersonic rarefied flow field expands outward linearly with the increase of 1/2Kn.An empirical method is proposed,which can be used for the prediction of the hypersonic flow-field structure at a given inflow velocity,especially the shock wave position.展开更多
A direct numerical simulation（DNS） is performed to investigate the control effect and mechanism of turbulent channel flow with the distribution of spanwise Lorentz force. A sinusoidal distribution of constant spanwi...A direct numerical simulation（DNS） is performed to investigate the control effect and mechanism of turbulent channel flow with the distribution of spanwise Lorentz force. A sinusoidal distribution of constant spanwise Lorentz force is selected, of which the control effects, such as flow characters, mean Reynolds stress, and drag reductions, at different parameters of amplitude A and wave number k_x are discussed. The results indicate that the control effects vary with the parameter A and k_x. With the increase of A, the drag reduction rate D_r first increases and then decreases rapidly at low k_x,and slowly at high k_x. The low drag reduction（or even drag increase） is due to a weak suppression or even the enhancements of the random velocity fluctuation and mean Reynolds stress. The efficient drag reduction is due to the quasi-streamwise vortex structure induced by Lorentz force, which contributes to suppressing the random velocity fluctuation and mean Reynolds stress, and the negative vorticity improves the distribution of streamwise velocity. Therefore, the optimal control effect with a drag reduction of up to 58% can be obtained.展开更多
We report a numerical study of the Prandtl-number(Pr)effects in two-dimensional turbulent Rayleigh-Bénard convection.The simulations were conducted in a square box over the Pr range from 0.25 to 100 and over the ...We report a numerical study of the Prandtl-number(Pr)effects in two-dimensional turbulent Rayleigh-Bénard convection.The simulations were conducted in a square box over the Pr range from 0.25 to 100 and over the Rayleigh number(Ra)range from 10^(7) to 10^(10).We find that both the strength and the stability of the large-scale flow decrease with the increasing of Pr,and the flow pattern becomes plume-dominated at high Pr.The evolution in flow pattern is quantified by the Reynolds number(Re),with the Ra and the Pr scaling exponents varying from 0.54 to 0.67 and-0.87 to-0.93,respectively.It is further found that the non-dimensional heat flux at small Ra diverges strongly for different Pr,but their difference becomes marginal as Ra increases.For the thermal boundary layer,the spatially averaged thicknesses for all the Pr numbers can be described byδθ~Ra^(-0.30) approximately,but the local values vary a lot for different Pr,which become more uniform with Pr increasing.展开更多
This work is motivated by previous experimental and numerical studies which reveal that the hairpin vortex could be formed by the interaction between spanwise adjacent low-speed streaks. To prove that such an interact...This work is motivated by previous experimental and numerical studies which reveal that the hairpin vortex could be formed by the interaction between spanwise adjacent low-speed streaks. To prove that such an interaction mechanism is still applicable in the normal direction, two sinuous low-speed streaks with the same streamwise phase are set to be in the upper half and bottom half of a small size channel, respectively, and their evolution and interaction are investigated by direct numerical simulation. A new kind of hairpin-like vortical structure, distributed in the normal direction and straddled across both halves of the channel, is found during the cross-interaction process of the low-speed streaks. The influence of such a normal-distributed hairpin-like vortex （NHV） on the turbulent statistical regularity is also revealed. It is observed that the NHV can lead to a sudden surge of wall skin friction, but the value of the normal velocity as well as the streamwise and spanwise vorticity sharply decrease to zero in the center of the channel.展开更多
Direct numerical simulations （DNS） have now become a well established tool to examine complex multiphase flows. Such flows typically exhibit a large range of scales and it is generally necessary to use different des...Direct numerical simulations （DNS） have now become a well established tool to examine complex multiphase flows. Such flows typically exhibit a large range of scales and it is generally necessary to use different descriptions of the flow depending on the scale that we are examining. Here we discuss multiphase flows from a multiscale perspective. Those include both how DNS are providing insight and understanding for modeling of scales much larger than the ＂dominant scale＂ （defined where surface tension, viscous forces or inertia are important）, as well as how DNS are often limited by the need to resolve processes taking place on much smaller scales. Both problems can be cast into a language introduced for general classes of multiscale problems and reveal that while the classification may be new, the issues are not.展开更多
Wall temperature significantly affects stability and receptivity of the boundary layer. Changing the wall temperature locally may therefore be an effective laminar flow control technique. However, the situation is com...Wall temperature significantly affects stability and receptivity of the boundary layer. Changing the wall temperature locally may therefore be an effective laminar flow control technique. However, the situation is complicated when the wall temperature distribution is nonuniform, and researchers have experimentally found that local wall cooling may delay the onset of transition. We attempt to clarify the physical mechanisms whereby the local wall temperature affects the transition and the stability of a hypersonic boundary layer. A numerical investigation of the disturbance evolution in a Mach-6 sharp cone boundary layer with local wall heating or cooling is conducted. Direct numerical simulation(DNS) is performed for the single-frequency and broadband disturbance evolution caused by random forcing. We vary the local wall temperature and the location of heating/cooling, and then use the eNmethod to estimate the transition onset. Our results show that local wall cooling amplifies high-frequency unstable waves while stabilizing low-frequency unstable waves, with local heating amplifying all unstable waves locally. The disturbance amplitude and second-mode peak frequency obtained by DNS agree well with the previous experimental results. Local cooling/heating has a dual effect on the stability of the hypersonic boundary layer. For local cooling, while it effectively inhibits the growth of the low-frequency unstable waves that dominate the transition downstream, it also further destabilizes the downstream flow. In addition, while upstream cooling can delay the transition, excessive cooling may promote it;local heating always slightly promotes the transition.Finally, recommendations are given for practical engineering applications based on the present results.展开更多
基金Supported by the National Natural Science Foundation of China(10672052)the Advanced TalentStart-Up Foundation of Jiangsu University(08JDG018)~~
文摘The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investigated. Results show that the disturbance velocity presents the sine vibration features with the change of time, and the vibration period is identical to the vibration of the local wall. The disturbance velocity presents the fluctuation feature downstream, and the streamwise wave length approximates to the results from the Orr-Sommerfeld equation （OSE）. The growth rate from direct numerical simulation（DNS） is a little greater than that from OSE, and their trends are almost consistent. Under the condition of Re= 2 800, the disturbance amplitude gradually grows in the given computational region with the period T=30. However, it firstly increases and then decreases with the period T= 20. The disturbance harmonic of the former is obviously larger than that of the latter. The maximum streamwise and vertical disturbance velocities from DNS do not fully coincide with those from OSE at the vicinity of the local vibration wall, but coincide well with the former when they travel downstream. The 2-D disturbance induced by the local micro-vibration represents the form of Tollmien-Schlichting （T-S） wave on the boundary layer.
基金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.
基金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 National Natural Science Foundation of China under Grant Nos 31200545,11274206 and 11574184
文摘We report the results of protein folding （219M, C34, N36, 2KES, 2KHK） by the method of accelerated molecular dynamics （aMD） at room temperature with the implicit solvent model. Starting from the linear structures, these proteins successfully fold to the native structure in a lO0-ns aMD simulation. In contrast, they are failed under the traditional MD simulation in the same simulation time. Then we find that the lowest root mean square deviations of helix structures from the native structures are 0.36 A, 0.63 A, 0.52 A, 1.1 A and 0.78 A. What is more, native contacts, cluster and free energy analyses show that the results of the aMD method are in accordance with the experiment very well. All analyses show that the aMD can accelerate the simulation process, thus we may apply it to the field of computer aided drug designs.
基金Project supported by the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (Grant No. 51121004)the National Natural Science Foundation of China (Grant No. 51076036)+2 种基金the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20112302110020)the Fundamental Research Funds for the Central Universities,China (Grant No. HIT.BRET1.2010008)the Japan Society for the Promotion of Science (JSPS) Research Fellowship
文摘In this paper,we present a direct numerical simulation（DNS） of elastic turbulence of viscoelastic fluid at vanishingly low Reynolds number（Re = 1） in a three-dimensional straight channel flow for the first time,using the Giesekus constitutive model for the fluid.In order to generate and maintain the turbulent fluid motion in the straight channel,a sinusoidal force term is added to the momentum equation,and then the elastic turbulence is numerically realized with an initialized chaotic velocity field and a stretched conformation field.Statistical and structural characteristics of the elastic turbulence therein are analyzed based on the detailed information obtained from the DNS.The fluid mixing enhancement effect of elastic turbulence is also demonstrated for the potential applications of this phenomenon.
基金Project supported by the National Key Research and Development Program of China (Grant No.2022YFE03050001)partly by the National Natural Science Foundation of China (Grant No.12175160)the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD)。
文摘The direct implicit particle-in-cell is a powerful kinetic method for researching plasma characteristics.However,it is time-consuming to obtain the future electromagnetic field in such a method since the field equations contain time-dependent matrix coefficients.In this work,we propose to explicitly push particles and obtain the future electromagnetic field based on the information about the particles in the future.The new method retains the form of implicit particle pusher,but the future field is obtained by solving the traditional explicit equation.Several numerical experiments,including the motion of charged particle in electromagnetic field,plasma sheath,and free diffusion of plasma into vacuum,are implemented to evaluate the performance of the method.The results demonstrate that the proposed method can suppress finite-grid-instability resulting from the coarse spatial resolution in electron Debye length through the strong damping of high-frequency plasma oscillation,while accurately describe low-frequency plasma phenomena,with the price of losing the numerical stability at large time-step.We believe that this work is helpful for people to research the bounded plasma by using particle-in-cell simulations.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12372251 and 12132015)the Fundamental Research Funds for the Provincial Universities of Zhejiang(Grant No.2023YW69)。
文摘The swimming performance of rod-shaped microswimmers in a channel was numerically investigated using the two-dimensional lattice Boltzmann method(LBM).We considered variable-length squirmer rods,assembled from circular squirmer models with self-propulsion mechanisms,and analyzed the effects of the Reynolds number(Re),aspect ratio(ε),squirmer-type factor(β)and blockage ratio(κ)on swimming efficiency(η)and power expenditure(P).The results show no significant difference in power expenditure between pushers(microswimmers propelled from the tail)and pullers(microswimmers propelled from the head)at the low Reynolds numbers adopted in this study.However,the swimming efficiency of pushers surpasses that of pullers.Moreover,as the degree of channel blockage increases(i.e.,κincreases),the squirmer rod consumes more energy while swimming,and its swimming efficiency also increases,which is clearly reflected whenε≤3.Notably,squirmer rods with a larger aspect ratioεand aβvalue approaching 0 can achieve high swimming efficiency with lower power expenditure.The advantages of self-propelled microswimmers are manifested whenε>4 andβ=±1,where the squirmer rod consumes less energy than a passive rod driven by an external field.These findings underscore the potential for designing more efficient microswimmers by carefully considering the interactions between the microswimmer geometry,propulsion mechanism and fluid dynamic environment.
文摘Turbulent spots play a key role in the formation of the turbulence and the transition. The generation and evolution of turbulent spots using the wall impulse model in the plane Couette flow are studied by direct numerical simulation of Navier-Stokes equations. A group of three-dimensional coupling compact difference schemes with high accuracy and high resolution is used in the numerical calculation. The important characteristics of turbulent spots based on the results of examples are analyzed, including the formation of random pulse, the generation of Reynolds stress, the growth of disturbance amplitude, and the continuous change of spot shape, especially the complex evolution process of the streamwise vortices. Computational results confirm that basic properties of turbulent spots in the laminar flow are similar to those in the turbulent flow.
文摘Numerical simulations are presented about the effects of gas rarefaction on hypersonic flow field.Due to the extremely difficult experiment,limited wind-tunnel conditions and high cost,most problems in rarefied flow regime are investigated through numerical methods,in which the direct simulation Monte-Carlo(DSMC)method is widely adopted.And the unstructured DSMC method is employed here.Flows around a vertical plate at a given velocity 7 500 m/s are simulated.For gas rarefaction is judged by the free-stream Knudsen number(Kn),two vital factors are considered:molecular number density and the plate′s length.Cases in which Kn varies from 0.035 to13.36 are simulated.Flow characters in the whole rarefied regime are described,and flow-field structure affected by Knis analyzed.Then,the dimensionless position D*of a certain velocity in the stagnation line is chosen as the marker of flow field to measure its variation.Through flow-field tracing and least-square numerical method analyzing,it is proved that hypersonic rarefied flow field expands outward linearly with the increase of 1/2Kn.An empirical method is proposed,which can be used for the prediction of the hypersonic flow-field structure at a given inflow velocity,especially the shock wave position.
基金supported by the National Natural Science Foundation of China(Grant Nos.11672135 and 11202102)the Fundamental Research Funds for the Central Universities,China(Grant No.30916011347)a Foundation for the Author of National Excellent Doctoral Dissertation,China(Grant No.201461)
文摘A direct numerical simulation（DNS） is performed to investigate the control effect and mechanism of turbulent channel flow with the distribution of spanwise Lorentz force. A sinusoidal distribution of constant spanwise Lorentz force is selected, of which the control effects, such as flow characters, mean Reynolds stress, and drag reductions, at different parameters of amplitude A and wave number k_x are discussed. The results indicate that the control effects vary with the parameter A and k_x. With the increase of A, the drag reduction rate D_r first increases and then decreases rapidly at low k_x,and slowly at high k_x. The low drag reduction（or even drag increase） is due to a weak suppression or even the enhancements of the random velocity fluctuation and mean Reynolds stress. The efficient drag reduction is due to the quasi-streamwise vortex structure induced by Lorentz force, which contributes to suppressing the random velocity fluctuation and mean Reynolds stress, and the negative vorticity improves the distribution of streamwise velocity. Therefore, the optimal control effect with a drag reduction of up to 58% can be obtained.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11961160719,11702128,91752201,and 11772362)the Shenzhen Fundamental Research Program(Grant No.JCYJ20190807160413162)+1 种基金the Fundamental Research Funds for the Central Universities(Sun Yat-sen University under Grant No.19lgzd15)the Department of Science and Technology of Guangdong Province,China(Grant No.2019B21203001).
文摘We report a numerical study of the Prandtl-number(Pr)effects in two-dimensional turbulent Rayleigh-Bénard convection.The simulations were conducted in a square box over the Pr range from 0.25 to 100 and over the Rayleigh number(Ra)range from 10^(7) to 10^(10).We find that both the strength and the stability of the large-scale flow decrease with the increasing of Pr,and the flow pattern becomes plume-dominated at high Pr.The evolution in flow pattern is quantified by the Reynolds number(Re),with the Ra and the Pr scaling exponents varying from 0.54 to 0.67 and-0.87 to-0.93,respectively.It is further found that the non-dimensional heat flux at small Ra diverges strongly for different Pr,but their difference becomes marginal as Ra increases.For the thermal boundary layer,the spatially averaged thicknesses for all the Pr numbers can be described byδθ~Ra^(-0.30) approximately,but the local values vary a lot for different Pr,which become more uniform with Pr increasing.
基金Project supported by the Natural Science Foundation of Zhejiang Province,China(Grant No.LQ18A020002)the National Natural Science Foundation of China(Grant No.41776105)the Foundation of the Educational Committee of Zhejiang Province,China(Grant No.Y201737053)
文摘This work is motivated by previous experimental and numerical studies which reveal that the hairpin vortex could be formed by the interaction between spanwise adjacent low-speed streaks. To prove that such an interaction mechanism is still applicable in the normal direction, two sinuous low-speed streaks with the same streamwise phase are set to be in the upper half and bottom half of a small size channel, respectively, and their evolution and interaction are investigated by direct numerical simulation. A new kind of hairpin-like vortical structure, distributed in the normal direction and straddled across both halves of the channel, is found during the cross-interaction process of the low-speed streaks. The influence of such a normal-distributed hairpin-like vortex （NHV） on the turbulent statistical regularity is also revealed. It is observed that the NHV can lead to a sudden surge of wall skin friction, but the value of the normal velocity as well as the streamwise and spanwise vorticity sharply decrease to zero in the center of the channel.
文摘Direct numerical simulations （DNS） have now become a well established tool to examine complex multiphase flows. Such flows typically exhibit a large range of scales and it is generally necessary to use different descriptions of the flow depending on the scale that we are examining. Here we discuss multiphase flows from a multiscale perspective. Those include both how DNS are providing insight and understanding for modeling of scales much larger than the ＂dominant scale＂ （defined where surface tension, viscous forces or inertia are important）, as well as how DNS are often limited by the need to resolve processes taking place on much smaller scales. Both problems can be cast into a language introduced for general classes of multiscale problems and reveal that while the classification may be new, the issues are not.
基金supported the National Natural Science Foundation of China (Grant No. 92271102)。
文摘Wall temperature significantly affects stability and receptivity of the boundary layer. Changing the wall temperature locally may therefore be an effective laminar flow control technique. However, the situation is complicated when the wall temperature distribution is nonuniform, and researchers have experimentally found that local wall cooling may delay the onset of transition. We attempt to clarify the physical mechanisms whereby the local wall temperature affects the transition and the stability of a hypersonic boundary layer. A numerical investigation of the disturbance evolution in a Mach-6 sharp cone boundary layer with local wall heating or cooling is conducted. Direct numerical simulation(DNS) is performed for the single-frequency and broadband disturbance evolution caused by random forcing. We vary the local wall temperature and the location of heating/cooling, and then use the eNmethod to estimate the transition onset. Our results show that local wall cooling amplifies high-frequency unstable waves while stabilizing low-frequency unstable waves, with local heating amplifying all unstable waves locally. The disturbance amplitude and second-mode peak frequency obtained by DNS agree well with the previous experimental results. Local cooling/heating has a dual effect on the stability of the hypersonic boundary layer. For local cooling, while it effectively inhibits the growth of the low-frequency unstable waves that dominate the transition downstream, it also further destabilizes the downstream flow. In addition, while upstream cooling can delay the transition, excessive cooling may promote it;local heating always slightly promotes the transition.Finally, recommendations are given for practical engineering applications based on the present results.
