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.展开更多
Multi-dimensional direct numerical simulations （DNS） of astrophysical detonations in degenerate matter have revealed that the nuclear burning is typically characterized by cellular structure caused by transverse ins...Multi-dimensional direct numerical simulations （DNS） of astrophysical detonations in degenerate matter have revealed that the nuclear burning is typically characterized by cellular structure caused by transverse instabilities in the detonation front. Type Ia supernova modelers often use one- dimensional DNS of detonations as inputs or constraints for their whole star simulations. While these one-dimensional studies are useful tools, the true nature of the detonation is multi-dimensional. The multi-dimensional structure of the burning influences the speed, stability, and the composition of the detonation and its burning products, and therefore, could have an impact on the spectra of Type Ia supernovae. Considerable effort has been expended modeling Type Ia supernovae at densities above 1 ~ 107 g.cm-3 where the complexities of turbulent burning dominate the flame propagation. How- ever, most full star models turn the nuclear burning schemes off when the density falls below 1 ~ 107 g.cm-3 and distributed burning begins. The deflagration to detonation transition （DDT） is believed to occur at just these densities and consequently they are the densities important for studying the properties of the subsequent detonation. This work will review the status of DNS studies of detona- tions and their possible implications for Type Ia supernova models. It will cover the development of Detonation theory from the first simple Chapman-Jouguet （C J） detonation models to the current models based on the time-dependent, compressible, reactive flow Euler equations of fluid dynamics.展开更多
Turbulent channel flows with consideration of the buoyancy effect of the bubble phase is investigated by means of the Direct Numerical Simulation （DNS）. This two-phase system is solved by a two-way coupling Lagrangi...Turbulent channel flows with consideration of the buoyancy effect of the bubble phase is investigated by means of the Direct Numerical Simulation （DNS）. This two-phase system is solved by a two-way coupling Lagrangian-Eulerian approach. The Reynolds number based on the friction velocity and the half-width of the channel is 194, and the gravitational acceleration varies from -0.5 to 0.5, ranging from the upflow to the downflow cases. This study aims to reveal the influence of buoyancy on the turbulence behavior and the bubble motion. Some typical statistical quantities, including the averaged velocities and velocity fluctuations for the fluid and bubble phases, as well as the flow structures of the turbulence fluctuations, are analyzed.展开更多
Shock wave/boundary layer interaction in a 24°turning angle of the compression ramp at Mach number 2.9 controlled by steady microjet is investigated using direct numerical simulation.Three different jet spacings ...Shock wave/boundary layer interaction in a 24°turning angle of the compression ramp at Mach number 2.9 controlled by steady microjet is investigated using direct numerical simulation.Three different jet spacings which are termed as sparse,moderate and dense are considered,and the induced vortex system and shock structures are compared.A moderate jet spacing configuration is found to generate counter-rotating vortex pairs that transport high-momentum fluid towards the vicinity of wall and strengthen the boundary layer to resist separation,reducing the separation region.The dense jet spacing configuration creates a larger momentum deficit region,reducing the friction downstream of the corner.Analysis of pressure and pressure gradient reveals that dense jet spacing configuration reduces the intensity of separation shock.The impact of varying jet spacings on the turbulent kinetic energy transport mechanism is also investigated by decomposing the budget terms in the transport equation.Furthermore,the spectral characteristics of the separation region are studied using power spectral density and dynamic mode decomposition methods,revealing that moderate jet spacing configuration suppresses low-frequency fluctuations in the separation region.展开更多
We successfully perform the three-dimensional tracking in a turbulent fluid flow of small axisymmetrical particles that are neutrally-buoyant and bottom-heavy,i.e.,they have a non-homogeneous mass distribu-tion along ...We successfully perform the three-dimensional tracking in a turbulent fluid flow of small axisymmetrical particles that are neutrally-buoyant and bottom-heavy,i.e.,they have a non-homogeneous mass distribu-tion along their symmetry axis.We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate.The experiment is complemented by a series of simulations based on realistic Navier-Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravi-tational torque stability number and by the particle aspect ratio.We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential ori-entation and tumbling rate of the particles.We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.展开更多
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.展开更多
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.展开更多
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.展开更多
The interaction of an impinging oblique shock wave with an angle of 30°and a supersonic turbulent boundary layer at Ma_(∞)=2.9 and Re_(θ)=2400 over a wavy-wall is investigated through direct numerical simulatio...The interaction of an impinging oblique shock wave with an angle of 30°and a supersonic turbulent boundary layer at Ma_(∞)=2.9 and Re_(θ)=2400 over a wavy-wall is investigated through direct numerical simulation and compared with the interaction on a flat-plate under the same flow conditions.A sinusoidal wave with amplitude to wavelength ratio of 0.26 moves in the streamwise direction and is uniformly distributed across the spanwise direction.The influences of the wavy-wall on the interaction,including the characterization of the flow field,the skin-friction,pressure and the budget of turbulence kinetic energy,are systematically studied.The region of separation grows slightly and decomposes into four bubbles.Local peaks of skin-friction are observed at the rear part of the interaction region.The low-frequency shock motion can be seen in the wall pressure spectra.Analyses of the turbulence kinetic energy budget indicate that both diffusion and transport significantly increase near the crests,balanced by an amplified dissipation in the near-wall region.Proper orthogonal decomposition analyses show that the most energetic structures are associated with the separated shock and the shear layer over the bubbles.Only the bubbles in the first two troughs are dominated by a low-frequency enlargement or shrinkage.展开更多
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.展开更多
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.展开更多
The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the n...The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boun...Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.展开更多
In this article dedicated to the modeling of vertical mass transfers between the biofilm and the bulk flow, we have, in the first instance, presented the methodology used, followed by the presentation of various resul...In this article dedicated to the modeling of vertical mass transfers between the biofilm and the bulk flow, we have, in the first instance, presented the methodology used, followed by the presentation of various results obtained through analyses conducted on velocity fields, different fluxes, and overall transfer coefficients. Due to numerical constraints (resolution of relevant spatial scales), we have restricted the analysis to low Schmidt numbers (S<sub>c</sub><sub></sub>=0.1, S<sub>c</sub></sub>=1, and S<sub>c</sub></sub>=10) and a single roughness Reynolds number (Re<sub>*</sub>=150). The analysis of instantaneous concentration fields from various simulations revealed logarithmic concentration profiles above the canopy. In this zone, the concentration is relatively homogeneous for longer times. The analysis of results also showed that the contribution of molecular diffusion to the total flux depends on the Schmidt number. This contribution is negligible for Schmidt numbers S<sub>c</sub></sub>≥0.1, but nearly balances the turbulent flux for S<sub>c</sub></sub>=0.1. In the canopy, the local Sherwood number, given by the ratio of the total flux (within or above the canopy) to the molecular diffusion flux at the wall, also depends on the Schmidt number and varies significantly between the canopy and the region above. The exchange velocity, a purely hydrodynamic parameter, is independent of the Schmidt number and is on the order of 10% of in the present case. This study also reveals that nutrient absorption by organisms near the wall depends on the Schmidt number. Such absorption is facilitated by lower Schmidt numbers.展开更多
Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models(GCMs)but also in various higher-resoluti...Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models(GCMs)but also in various higher-resolution limited-area models such as cloud-resolving models(CRMs)and large-eddy simulation(LES)models.Instead of giving a comprehensive review of existing microphysical parameterizations that have been developed over the years,this study concentrates purposely on several topics that we believe are understudied but hold great potential for further advancing bulk microphysics parameterizations:multi-moment bulk microphysics parameterizations and the role of the spectral shape of hydrometeor size distributions;discrete vs“continuous”representation of hydrometeor types;turbulence-microphysics interactions including turbulent entrainment-mixing processes and stochastic condensation;theoretical foundations for the mathematical expressions used to describe hydrometeor size distributions and hydrometeor morphology;and approaches for developing bulk microphysics parameterizations.Also presented are the spectral bin scheme and particle-based scheme(especially,super-droplet method)for representing explicit microphysics.Their advantages and disadvantages are elucidated for constructing cloud models with detailed microphysics that are essential to developing processes understanding and bulk microphysics parameterizations.Particle-resolved direct numerical simulation(DNS)models are described as an emerging technique to investigate turbulence-microphysics interactions at the most fundamental level by tracking individual particles and resolving the smallest turbulent eddies in turbulent clouds.Outstanding challenges and future research directions are explored as well.展开更多
文摘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.
文摘Multi-dimensional direct numerical simulations （DNS） of astrophysical detonations in degenerate matter have revealed that the nuclear burning is typically characterized by cellular structure caused by transverse instabilities in the detonation front. Type Ia supernova modelers often use one- dimensional DNS of detonations as inputs or constraints for their whole star simulations. While these one-dimensional studies are useful tools, the true nature of the detonation is multi-dimensional. The multi-dimensional structure of the burning influences the speed, stability, and the composition of the detonation and its burning products, and therefore, could have an impact on the spectra of Type Ia supernovae. Considerable effort has been expended modeling Type Ia supernovae at densities above 1 ~ 107 g.cm-3 where the complexities of turbulent burning dominate the flame propagation. How- ever, most full star models turn the nuclear burning schemes off when the density falls below 1 ~ 107 g.cm-3 and distributed burning begins. The deflagration to detonation transition （DDT） is believed to occur at just these densities and consequently they are the densities important for studying the properties of the subsequent detonation. This work will review the status of DNS studies of detona- tions and their possible implications for Type Ia supernova models. It will cover the development of Detonation theory from the first simple Chapman-Jouguet （C J） detonation models to the current models based on the time-dependent, compressible, reactive flow Euler equations of fluid dynamics.
基金supported by the National Natural Science Foundation of China(Grant Nos.10772173,10972211 and 11072236)the Fundamental Research Funds for the Central Universitiesthe Science and Technology Innovation Foundation of the Chinese Academy of Sciences(Grant No.CXJJ-11-M69)
文摘Turbulent channel flows with consideration of the buoyancy effect of the bubble phase is investigated by means of the Direct Numerical Simulation （DNS）. This two-phase system is solved by a two-way coupling Lagrangian-Eulerian approach. The Reynolds number based on the friction velocity and the half-width of the channel is 194, and the gravitational acceleration varies from -0.5 to 0.5, ranging from the upflow to the downflow cases. This study aims to reveal the influence of buoyancy on the turbulence behavior and the bubble motion. Some typical statistical quantities, including the averaged velocities and velocity fluctuations for the fluid and bubble phases, as well as the flow structures of the turbulence fluctuations, are analyzed.
基金Supported by the National Natural Science Foundation of China(Nos.11972064,92052104)the Key Laboratory of Computational Aerodynamics,AVIC Aerodynamics Research Institute,China(No.YL2022XFX0405)the Fundamental Research Funds for the Central Universities,China.
文摘Shock wave/boundary layer interaction in a 24°turning angle of the compression ramp at Mach number 2.9 controlled by steady microjet is investigated using direct numerical simulation.Three different jet spacings which are termed as sparse,moderate and dense are considered,and the induced vortex system and shock structures are compared.A moderate jet spacing configuration is found to generate counter-rotating vortex pairs that transport high-momentum fluid towards the vicinity of wall and strengthen the boundary layer to resist separation,reducing the separation region.The dense jet spacing configuration creates a larger momentum deficit region,reducing the friction downstream of the corner.Analysis of pressure and pressure gradient reveals that dense jet spacing configuration reduces the intensity of separation shock.The impact of varying jet spacings on the turbulent kinetic energy transport mechanism is also investigated by decomposing the budget terms in the transport equation.Furthermore,the spectral characteristics of the separation region are studied using power spectral density and dynamic mode decomposition methods,revealing that moderate jet spacing configuration suppresses low-frequency fluctuations in the separation region.
基金supported by the National Natural Science Foundation of China (Grant 11988102)
文摘We successfully perform the three-dimensional tracking in a turbulent fluid flow of small axisymmetrical particles that are neutrally-buoyant and bottom-heavy,i.e.,they have a non-homogeneous mass distribu-tion along their symmetry axis.We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate.The experiment is complemented by a series of simulations based on realistic Navier-Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravi-tational torque stability number and by the particle aspect ratio.We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential ori-entation and tumbling rate of the particles.We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.
基金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.
基金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.
基金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.
基金co-supported by the National Key Research and Development Program of China(Nos.2019YFA0405300 and 2016YFA0401200)the National Natural Science Foundation of China(Nos.11972356 and 91852203)+1 种基金LHD Youth Innovation Fund(No.LHD2019CX04)National Numerical Wind Tunnel Project。
文摘The interaction of an impinging oblique shock wave with an angle of 30°and a supersonic turbulent boundary layer at Ma_(∞)=2.9 and Re_(θ)=2400 over a wavy-wall is investigated through direct numerical simulation and compared with the interaction on a flat-plate under the same flow conditions.A sinusoidal wave with amplitude to wavelength ratio of 0.26 moves in the streamwise direction and is uniformly distributed across the spanwise direction.The influences of the wavy-wall on the interaction,including the characterization of the flow field,the skin-friction,pressure and the budget of turbulence kinetic energy,are systematically studied.The region of separation grows slightly and decomposes into four bubbles.Local peaks of skin-friction are observed at the rear part of the interaction region.The low-frequency shock motion can be seen in the wall pressure spectra.Analyses of the turbulence kinetic energy budget indicate that both diffusion and transport significantly increase near the crests,balanced by an amplified dissipation in the near-wall region.Proper orthogonal decomposition analyses show that the most energetic structures are associated with the separated shock and the shear layer over the bubbles.Only the bubbles in the first two troughs are dominated by a low-frequency enlargement or shrinkage.
基金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.
基金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 （50606026, 50736006）.
文摘The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.
基金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 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.
基金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 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(Nos.12202456 and12172360)the Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”of the National Natural Science Foundation of China(No.11988102)the China Postdoctoral Science Foundation(No.2021M693241)。
文摘Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.
文摘In this article dedicated to the modeling of vertical mass transfers between the biofilm and the bulk flow, we have, in the first instance, presented the methodology used, followed by the presentation of various results obtained through analyses conducted on velocity fields, different fluxes, and overall transfer coefficients. Due to numerical constraints (resolution of relevant spatial scales), we have restricted the analysis to low Schmidt numbers (S<sub>c</sub><sub></sub>=0.1, S<sub>c</sub></sub>=1, and S<sub>c</sub></sub>=10) and a single roughness Reynolds number (Re<sub>*</sub>=150). The analysis of instantaneous concentration fields from various simulations revealed logarithmic concentration profiles above the canopy. In this zone, the concentration is relatively homogeneous for longer times. The analysis of results also showed that the contribution of molecular diffusion to the total flux depends on the Schmidt number. This contribution is negligible for Schmidt numbers S<sub>c</sub></sub>≥0.1, but nearly balances the turbulent flux for S<sub>c</sub></sub>=0.1. In the canopy, the local Sherwood number, given by the ratio of the total flux (within or above the canopy) to the molecular diffusion flux at the wall, also depends on the Schmidt number and varies significantly between the canopy and the region above. The exchange velocity, a purely hydrodynamic parameter, is independent of the Schmidt number and is on the order of 10% of in the present case. This study also reveals that nutrient absorption by organisms near the wall depends on the Schmidt number. Such absorption is facilitated by lower Schmidt numbers.
基金supported by the US Department of Energy(DOE)’s Office of Science Atmospheric Systems Research(ASR)Programthe Office of Energy Efficiency and Renewable Energy(EERE)Solar Energy Technologies Office(SETO)award(33504)+3 种基金the Brookhaven National Laboratory(BNL)’s Laboratory Directed Research&Development Program(LDRD)(22-065)The Brookhaven National Laboratory is operated by the Brookhaven Science Associates,LLC(BSA),for the US Department of Energy under Contract No.DESC0012704supported by JSPS KAKENHI Grant No.26286089MEXT KAKENHI Grant No.18H04448。
文摘Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models(GCMs)but also in various higher-resolution limited-area models such as cloud-resolving models(CRMs)and large-eddy simulation(LES)models.Instead of giving a comprehensive review of existing microphysical parameterizations that have been developed over the years,this study concentrates purposely on several topics that we believe are understudied but hold great potential for further advancing bulk microphysics parameterizations:multi-moment bulk microphysics parameterizations and the role of the spectral shape of hydrometeor size distributions;discrete vs“continuous”representation of hydrometeor types;turbulence-microphysics interactions including turbulent entrainment-mixing processes and stochastic condensation;theoretical foundations for the mathematical expressions used to describe hydrometeor size distributions and hydrometeor morphology;and approaches for developing bulk microphysics parameterizations.Also presented are the spectral bin scheme and particle-based scheme(especially,super-droplet method)for representing explicit microphysics.Their advantages and disadvantages are elucidated for constructing cloud models with detailed microphysics that are essential to developing processes understanding and bulk microphysics parameterizations.Particle-resolved direct numerical simulation(DNS)models are described as an emerging technique to investigate turbulence-microphysics interactions at the most fundamental level by tracking individual particles and resolving the smallest turbulent eddies in turbulent clouds.Outstanding challenges and future research directions are explored as well.
