Underhand cut-and-fill mining has been widely used in underground mining operations,especially when the rock mass or orebody is of poor quality or prone to rockburst due to high stress.In such cases,mining workers sho...Underhand cut-and-fill mining has been widely used in underground mining operations,especially when the rock mass or orebody is of poor quality or prone to rockburst due to high stress.In such cases,mining workers should carry out all production activities under the cemented backfill roof or sill mat instead of a highly fractured and unstable rock roof or a strong rock roof with a high potential of rockburst.Therefore,the stability and required strength of the sill mat are critical issues for mining engineers.In 1991,Mitchell considered that sill mat could fail by caving,sliding,rotation,and flexure.Mitchell also proposed an analytical solution to determine the minimum required strength of the sill mat for each type of failure based on two stiff or immobile rock walls.However,recent publications using numerical modeling and field measurements indicate that the compressive stresses in the sill mat induced by rock wall closure due to a stope excavation beneath the sill mat can be significant.It is thus highly necessary to investigate the required strength of the sill mat by considering rock wall closure.In this study,the crushing failure of sill mat due to rock wall closure generated by underground excavation and a new failure mode called"crushing and caving”is revealed by numerical modeling.An analytical solution corresponding to each failure mode is then developed to estimate the minimum required cohesion(cmin)of the sill mat.A criterion is also proposed to determine if the sill mat fails by crushing or crushing-and-caving failure.The proposed analytical solution does not involve any correction coefficients.The validity of the proposed analytical solution is demonstrated by numerical modeling.The proposed analytical solution can thus be employed to predict the cmin of sill mat subjected to wall closure generated by underlying stope excavation.展开更多
A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equ...A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.展开更多
The spatial-temporal evolution of coherent structures （CS） is significant for turbulence control and drag re- duction. Among the CS, low and high speed streak structures show typical burst phenomena. The analysis wa...The spatial-temporal evolution of coherent structures （CS） is significant for turbulence control and drag re- duction. Among the CS, low and high speed streak structures show typical burst phenomena. The analysis was based on a time series of three-dimensional and three-component （3D-3C） velocity fields of the flat plate turbulent boundary layer （TBL） measured by a Tomographic and Time-resolved PIV （Tomo TRPIV） system. Using multi-resolution wavelet transform and conditional sampling method, we extracted the intrinsic topologies and found that the streak structures appear in bar-like patterns. Furthermore, we seized locations and velocity information of transient CS, and then calculated the propagation velocity of CS based on spatial-temporal cross-correlation scanning. This laid a foundation for further studies on relevant dynamics properties.展开更多
Adopting Yoshizawa＇s two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at...Adopting Yoshizawa＇s two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.展开更多
With the increased penetration of wind energy in our nation’s energy portfolio, wind farms are placed in a way close to each other. Thus, their wakes have to be fully considered in the design and operation of a wind ...With the increased penetration of wind energy in our nation’s energy portfolio, wind farms are placed in a way close to each other. Thus, their wakes have to be fully considered in the design and operation of a wind farm. In this study, we investigate the wake of a wind farm using large-eddy simulation with wind turbine rotor modelled by the actuator disk model. The simulated results show that the wake of a wind farm can persist for a long distance in its downstream. For the considered wind farm layout, the velocity in the wake recovers 95% of that of the undisturbed inflow at 55 rotor diameters downstream from its last row, suggesting that the wake of a wind farm should be fully considered in the optimal design and operation for its downstream wind farms.展开更多
The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance in...The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Kairmain vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.展开更多
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.展开更多
Nonlinear energy transfer is represented through eddy viscosity and stochastic forcing within the framework of resolvent analysis.Previous investigations estimate the contribution of eddy-viscosity-enhanced resolvent ...Nonlinear energy transfer is represented through eddy viscosity and stochastic forcing within the framework of resolvent analysis.Previous investigations estimate the contribution of eddy-viscosity-enhanced resolvent opera-tor to nonlinear energy transfer.The present article estimates the contribution of stochastic forcing to nonlinear energy transfer and demonstrates that the contribution of stochastic forcing cannot be ignored.These results are achieved by numerically comparing the eddy-viscosity-enhanced resolvent operator and stochastic forcing with nonlinear energy transfer in turbulent channel flows.Furthermore,the numerical results indicate that composite resolvent operators can improve the prediction of nonlinear energy transfer.展开更多
In the underhand cut-and-fill mining method,a sill mat(i.e.an artificial horizontal pillar)constructed by cemented backfill is essential to prevent mine workers from being directly exposed under problematic rock roofs...In the underhand cut-and-fill mining method,a sill mat(i.e.an artificial horizontal pillar)constructed by cemented backfill is essential to prevent mine workers from being directly exposed under problematic rock roofs.A critical issue is to determine the minimum required strength of the sill mat to ensure a safe and cost-effective design.Until now,Mitchell’s analytical solution is the only available option,considering two stiff and immobile rock walls.Unavoidable rock wall closure associated with stope excavation below the sill mat was neglected.This,along with other undefined parameters,explains why Mitchell’s solution is rarely used in sill mat design.A new analytical solution for determining the minimum required strength of the sill mat accounting for wall closure is necessary.In this study,a closed-form analytical solution for estimating rock wall closure generated by stope excavation below a sill mat is developed by using Salamon’s and Flamant’s models.The proposed analytical solution does not contain any coefficients of correction or calibration.Despite several assumptions(or somewhat of oversimplifications)necessary to render a simple analytical solution possible,good agreements are obtained between the rock wall closures predicted by applying the proposed analytical solution and those obtained numerically with FLAC3D for many cases with arbitrarily chosen geometrical and material parameters.The proposed analytical solution is therefore validated and can be used to evaluate the rock wall closure generated by stope excavation below a sill mat.展开更多
In this paper, ultrasonic （20 kHz） fatigue tests were performed on specimens of a high-strength steel in very high cycle fatigue （VHCF） regime. Experimental results showed that for most tested specimens failed in ...In this paper, ultrasonic （20 kHz） fatigue tests were performed on specimens of a high-strength steel in very high cycle fatigue （VHCF） regime. Experimental results showed that for most tested specimens failed in a VHCF regime, a fatigue crack originated from the interior of specimen with a fish-eye pattern, which contained a fine granular area （FGA） centered by an inclusion as the crack origin. Then, a two-parameter model is proposed to predict the fatigue life of high-strength steels with fish-eye mode failure in a VHCF regime, which takes into account the inclusion size and the FGA size. The model was verified by the data of present experiments and those in the literature. Furthermore, an analytic formula was obtained for estimating the equivalent crack growth rate within the FGA. The results also indicated that the stress intensity factor range at the front of the FGA varies within a small range, which is irrespective of stress amplitude and fatigue life.展开更多
Flapping-powered propulsion is used by many animals to locomote through air or water. Here we review recent experimental and numerical studies on self-propelled mechanical systems powered by a flapping motion. These s...Flapping-powered propulsion is used by many animals to locomote through air or water. Here we review recent experimental and numerical studies on self-propelled mechanical systems powered by a flapping motion. These studies improve our understanding of the mutual interaction between actively flapping bodies and surrounding fluids. The results obtained in these works provide not only new insights into biolocomotion but also useful information for the biomimetic design of artificial flyers and swimmers.展开更多
A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the ...A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.展开更多
A new shedding phenomenon of ventilated partial cavitations is observed around an axisymmetric projectile in a horizontal launching experiment.The experiment system is established based on SHPB launching and high spee...A new shedding phenomenon of ventilated partial cavitations is observed around an axisymmetric projectile in a horizontal launching experiment.The experiment system is established based on SHPB launching and high speed photography.A numerical simulation is carried out based on the homogeneous mixture approach,and its predicted evolutions of cavities are compared with the experimental results.The cavity breaks off by the interaction between the gas injection and the re-entry jet at the middle location of the projectile,which is obviously different from natural cavitation.The mechanism of cavity breaking and shedding is investigated,and the influences of important factors are also discussed.展开更多
Eddy-damping quasinormal Markovian （EDQNM） theory is employed to calculate the resolved-scale spectrum and transfer spectrum, based on which we investigate the resolved-scale scaling law. Results show that the scali...Eddy-damping quasinormal Markovian （EDQNM） theory is employed to calculate the resolved-scale spectrum and transfer spectrum, based on which we investigate the resolved-scale scaling law. Results show that the scaling law of the resolved-scale turbulence, which is affected by several factors, is far from that of the full-scale turbulence and should be corrected. These results are then applied to an existing subgrid model to improve its performance. A series of simulations are performed to verify the necessity of a fixed scaling law in the subgrid modeling.展开更多
When using a miniature single sensor boundary layer probe, the time sequences of the stream-wise velocity in the turbulent boundary layer （TBL） are measured by using a hot wire anemometer. Beneath the fully develope...When using a miniature single sensor boundary layer probe, the time sequences of the stream-wise velocity in the turbulent boundary layer （TBL） are measured by using a hot wire anemometer. Beneath the fully developed TBL, the wall pressure fluctuations are attained by a microphone mechanism with high spatial resolution. Analysis on the statistic and spectrum properties of velocity and wall pressure reveals the relationship between the wall pressure fluctuation and the energy-containing structure in the buffer layer of the TBL. Wavelet transform shows the multi-scale natures of coherent structures contained in both signals of velocity and pressure. The most intermittent wall pressure scale is associated with the coherent structure in the buffer layer. Meanwhile the most energetic scale of velocity fluctuation at y＋ = 14 provides a specific frequency f9 ≈ 147 Hz for wall actuating control with Ret = 996.展开更多
The present experimental work focuses on a new model for space-time correlation and the scale-dependencies of convection velocity and sweep velocity in turbulent boundary layer over a flat wail. A turbulent boundary l...The present experimental work focuses on a new model for space-time correlation and the scale-dependencies of convection velocity and sweep velocity in turbulent boundary layer over a flat wail. A turbulent boundary layer flow at Reo = 2460 is measured by tomographic particle image velocimetry （tomographic PIV）. It is demonstrated that arch, cane, and hairpin vortices are dominant in the logarithmic layer. Hairpins and hairpin packets are responsible for the elongated low-momentum zones observed in the instantaneous flow field. The conditionally-averaged coherent structures systemically illustrate the key roles of hairpin vortice in the turbulence dynamic events, such as ejection and sweep events and energy transport. The space-time correlations of instantaneous streamwise fluctuation velocity are calculated and confirm the new elliptic model for the space-time correlation instead of Taylor hypothesis. The convection velocities derived from the space-time correlation and conditionally-averaged method both suggest the scaling with the local mean velocity in the logarithmic layer. Convection velocity result based on Fourier decomposition （FD） shows stronger scale- dependency in the spanwise direction than in streamwise direction. Compared with FD, the proper orthogonal decomposition （POD） has a distinct distribution of convection velocity for the large- and small-scales which are separated in light of their contributions of turbulent kinetic energy.展开更多
Turbulent flows over rough surfaces widely exist in nature and industry.Investigating its mechanism is of theoretical and practical significance.In this work we simulate the turbulent channel flow with rough walls usi...Turbulent flows over rough surfaces widely exist in nature and industry.Investigating its mechanism is of theoretical and practical significance.In this work we simulate the turbulent channel flow with rough walls using large-eddy simulation with rough elements resolved using the curvilinear immersed boundary method and compare the results obtained in this work with those in the paper by Yuan and Piomelli(J.Fluid Mech.,vol.760,pp.R1,2014),where the volume of fluid method was employed for modeling rough elements.The mean streamwise velocity profiles predicted by the two methods agree well with each other.Differences in Reynolds stresses and dispersive stresses are observed,which are attributed to the different approaches in dealing with the complex geometry of the rough surface.展开更多
An experimental measurement was performed us- ing time-resolved particle image velocimetry （TRPIV） to in- vestigate the spatial topological character of coherent struc- tures in wall-bounded turbulence of polymer ad...An experimental measurement was performed us- ing time-resolved particle image velocimetry （TRPIV） to in- vestigate the spatial topological character of coherent struc- tures in wall-bounded turbulence of polymer additive solu- tion. The fully developed near-wall turbulent flow fields with and without polymer additives at the same Reynolds number were measured by TRPIV in a water channel. The compar- isons of turbulent statistics confirm that due to viscoelastic structure of long-chain polymers, the wall-normal velocity fluctuation and Reynolds shear stress in the near-wall region are suppressed significantly. Furthermore, it is noted that such a behavior of polymers is closely related to the decease of the motion of the second and forth quadrants, i.e., the ejection and sweep events, in the near-wall region. The spa- tial topological mode of coherent structures during bursts has been extracted by the new mu-level criteria based on locally averaged velocity structure function. Although the general shapes of coherent structures are unchanged by polymer additives, the fluctuating velocity, velocity gradient, velocity strain rate and vorticity of coherent structures during burst events are suppressed in the polymer additive solution com- pared with that in water. The results show that due to the polymer additives the occurrence and intensity of coherent structures are suppressed, leading to drag reduction.展开更多
Flow around a real-life underwater vehicle often happens at a high Reynolds number with flow structures at different scales from the boundary layer around a blade to that around the hull. This poses a great challenge ...Flow around a real-life underwater vehicle often happens at a high Reynolds number with flow structures at different scales from the boundary layer around a blade to that around the hull. This poses a great challenge for large-eddy simulation of an underwater vehicle aiming at resolving all relevant flow scales. In this work, we propose to model the hull with appendages using the immersed boundary method, and model the propeller using the actuator disk model without resolving the geometry of the blade. The proposed method is then applied to simulate the flow around Defense Advanced Research Projects Agency(DARPA) suboff. An overall acceptable agreement is obtained for the pressure and friction coefficients. Complex flow features are observed in the near wake of suboff. In the far wake, the core region is featured by a jet because of the actuator disk, surrounded by an annular region with velocity deficit due to the body of suboff.展开更多
The mechanical properties of laser welded joints under impact loadings such as explosion and car crash etc. are critical for the engineering designs. The hardness, static and dynamic mechanical properties of AISI304 a...The mechanical properties of laser welded joints under impact loadings such as explosion and car crash etc. are critical for the engineering designs. The hardness, static and dynamic mechanical properties of AISI304 and AISI316L dissimilar stainless steel welded joints by CO2 laser were experimentally studied. The dynamic strainstress curves at the strain rate around 103 s-1 were obtained by the split Hopkinson tensile bar （SHTB）. The static mechanical properties of the welded joints have little changes with the laser power density and all fracture occurs at 316 L side. However, the strain rate sensitivity has a strong depen- dence on laser power density. The value of strain rate factor decreases with the increase of laser power density. The welded joint which may be applied for the impact loading can be obtained by reducing the laser power density in the case of welding quality assurance.展开更多
基金financial support from the Young Scientist Project of the National Key Research and Development Program of China(Grant No.2021YFC2900600)Beijing Nova Program(Grant No.20220484057)+1 种基金The authors acknowledge the financial support from the Natural Sciences and Engineering Research Council of Canada(Grant No.RGPIN-2018-06902)industrial partners of the Research Institute on Mines and the Environment(RIME UQAT-Polytechnique:https://irme.ca/en/).
文摘Underhand cut-and-fill mining has been widely used in underground mining operations,especially when the rock mass or orebody is of poor quality or prone to rockburst due to high stress.In such cases,mining workers should carry out all production activities under the cemented backfill roof or sill mat instead of a highly fractured and unstable rock roof or a strong rock roof with a high potential of rockburst.Therefore,the stability and required strength of the sill mat are critical issues for mining engineers.In 1991,Mitchell considered that sill mat could fail by caving,sliding,rotation,and flexure.Mitchell also proposed an analytical solution to determine the minimum required strength of the sill mat for each type of failure based on two stiff or immobile rock walls.However,recent publications using numerical modeling and field measurements indicate that the compressive stresses in the sill mat induced by rock wall closure due to a stope excavation beneath the sill mat can be significant.It is thus highly necessary to investigate the required strength of the sill mat by considering rock wall closure.In this study,the crushing failure of sill mat due to rock wall closure generated by underground excavation and a new failure mode called"crushing and caving”is revealed by numerical modeling.An analytical solution corresponding to each failure mode is then developed to estimate the minimum required cohesion(cmin)of the sill mat.A criterion is also proposed to determine if the sill mat fails by crushing or crushing-and-caving failure.The proposed analytical solution does not involve any correction coefficients.The validity of the proposed analytical solution is demonstrated by numerical modeling.The proposed analytical solution can thus be employed to predict the cmin of sill mat subjected to wall closure generated by underlying stope excavation.
基金Project supported by the National Natural Science Foundation of China(Nos.91752118,11672305,11232011,and 11572331)the Strategic Priority Research Program(No.XDB22040104)the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences(No.QYZDJ-SSWSYS002)
文摘A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.
基金supported by the National Natural Science Foundation of China(11332006,11272233,and 11411130150)the National Basic Research Programm of China(2012CB720101)
文摘The spatial-temporal evolution of coherent structures （CS） is significant for turbulence control and drag re- duction. Among the CS, low and high speed streak structures show typical burst phenomena. The analysis was based on a time series of three-dimensional and three-component （3D-3C） velocity fields of the flat plate turbulent boundary layer （TBL） measured by a Tomographic and Time-resolved PIV （Tomo TRPIV） system. Using multi-resolution wavelet transform and conditional sampling method, we extracted the intrinsic topologies and found that the streak structures appear in bar-like patterns. Furthermore, we seized locations and velocity information of transient CS, and then calculated the propagation velocity of CS based on spatial-temporal cross-correlation scanning. This laid a foundation for further studies on relevant dynamics properties.
基金Supported by the National Natural Science Foundation of China under Grant No 10472115, the Programme for New Century Excellent Talents in University of China, and the Opening Project of the State Key Laboratory of Nonlinear Mechanics.
文摘Adopting Yoshizawa＇s two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.
基金supported by the National Natural Science Foundation of China(Nos.11988102,12172360)Institute of Mechanics and Chinese Academy of Sciences。
文摘With the increased penetration of wind energy in our nation’s energy portfolio, wind farms are placed in a way close to each other. Thus, their wakes have to be fully considered in the design and operation of a wind farm. In this study, we investigate the wake of a wind farm using large-eddy simulation with wind turbine rotor modelled by the actuator disk model. The simulated results show that the wake of a wind farm can persist for a long distance in its downstream. For the considered wind farm layout, the velocity in the wake recovers 95% of that of the undisturbed inflow at 55 rotor diameters downstream from its last row, suggesting that the wake of a wind farm should be fully considered in the optimal design and operation for its downstream wind farms.
基金the National High-Tech Research and Development Program of China (863 Program)(2006AA09Z350)the Knowledge Innovation Program of the Chinese Academy of Sciences (KJCX2-YW-L02)
文摘The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Kairmain vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.
基金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.
基金supported by the National Natural Science Foundation of China(NSFC)Basic Science Center Program for Multiscale Problems in Nonlinear Mechanics(Grant No.11988102).
文摘Nonlinear energy transfer is represented through eddy viscosity and stochastic forcing within the framework of resolvent analysis.Previous investigations estimate the contribution of eddy-viscosity-enhanced resolvent opera-tor to nonlinear energy transfer.The present article estimates the contribution of stochastic forcing to nonlinear energy transfer and demonstrates that the contribution of stochastic forcing cannot be ignored.These results are achieved by numerically comparing the eddy-viscosity-enhanced resolvent operator and stochastic forcing with nonlinear energy transfer in turbulent channel flows.Furthermore,the numerical results indicate that composite resolvent operators can improve the prediction of nonlinear energy transfer.
基金financial support from the Young Scientist Project of the National Key Research and Development Program of China(Grant No.2021YFC2900600)the Beijing Nova Program(Grant No.20220484057)+1 种基金The authors acknowledge the financial support from the Natural Sciences and Engineering Research Council of Canada(Grant No.RGPIN-2018-06902)industrial partners of the Research Institute on Mines and the Environment(RIME UQAT-Polytechnique:https://irme.ca/en/).
文摘In the underhand cut-and-fill mining method,a sill mat(i.e.an artificial horizontal pillar)constructed by cemented backfill is essential to prevent mine workers from being directly exposed under problematic rock roofs.A critical issue is to determine the minimum required strength of the sill mat to ensure a safe and cost-effective design.Until now,Mitchell’s analytical solution is the only available option,considering two stiff and immobile rock walls.Unavoidable rock wall closure associated with stope excavation below the sill mat was neglected.This,along with other undefined parameters,explains why Mitchell’s solution is rarely used in sill mat design.A new analytical solution for determining the minimum required strength of the sill mat accounting for wall closure is necessary.In this study,a closed-form analytical solution for estimating rock wall closure generated by stope excavation below a sill mat is developed by using Salamon’s and Flamant’s models.The proposed analytical solution does not contain any coefficients of correction or calibration.Despite several assumptions(or somewhat of oversimplifications)necessary to render a simple analytical solution possible,good agreements are obtained between the rock wall closures predicted by applying the proposed analytical solution and those obtained numerically with FLAC3D for many cases with arbitrarily chosen geometrical and material parameters.The proposed analytical solution is therefore validated and can be used to evaluate the rock wall closure generated by stope excavation below a sill mat.
基金supported by the National Basic Research Program of China (Grant 2012CB937500)the National Natural Science Foundations of China (Grants 11172304 and 11202210)
文摘In this paper, ultrasonic （20 kHz） fatigue tests were performed on specimens of a high-strength steel in very high cycle fatigue （VHCF） regime. Experimental results showed that for most tested specimens failed in a VHCF regime, a fatigue crack originated from the interior of specimen with a fish-eye pattern, which contained a fine granular area （FGA） centered by an inclusion as the crack origin. Then, a two-parameter model is proposed to predict the fatigue life of high-strength steels with fish-eye mode failure in a VHCF regime, which takes into account the inclusion size and the FGA size. The model was verified by the data of present experiments and those in the literature. Furthermore, an analytic formula was obtained for estimating the equivalent crack growth rate within the FGA. The results also indicated that the stress intensity factor range at the front of the FGA varies within a small range, which is irrespective of stress amplitude and fatigue life.
基金supported by the Chinese Academy of Sciences (Grants KJCX-SW-L08, KJCX3-SYW-S01)the National Natural Science Foundation of China (Grants 11021262, 11023001, 11232011, 11372331)
文摘Flapping-powered propulsion is used by many animals to locomote through air or water. Here we review recent experimental and numerical studies on self-propelled mechanical systems powered by a flapping motion. These studies improve our understanding of the mutual interaction between actively flapping bodies and surrounding fluids. The results obtained in these works provide not only new insights into biolocomotion but also useful information for the biomimetic design of artificial flyers and swimmers.
基金supported by the National Natural Science Foundation of China (11302238, 11232011. and 11572331)support from the Strategic Priority Research Program (XDB22040104)+1 种基金the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences (QYZDJ-SSW-SYS002)the National Basic Research Program of China (973 Program 2013CB834100: Nonlinear science)
文摘A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.
文摘A new shedding phenomenon of ventilated partial cavitations is observed around an axisymmetric projectile in a horizontal launching experiment.The experiment system is established based on SHPB launching and high speed photography.A numerical simulation is carried out based on the homogeneous mixture approach,and its predicted evolutions of cavities are compared with the experimental results.The cavity breaks off by the interaction between the gas injection and the re-entry jet at the middle location of the projectile,which is obviously different from natural cavitation.The mechanism of cavity breaking and shedding is investigated,and the influences of important factors are also discussed.
基金supported by the National Natural Science Foundation of China(11202013 and 51136003)the National Basic Research Program of China(2012CB720200)the Opening fundof State Key Laboratory of Nonlinear Mechanics
文摘Eddy-damping quasinormal Markovian （EDQNM） theory is employed to calculate the resolved-scale spectrum and transfer spectrum, based on which we investigate the resolved-scale scaling law. Results show that the scaling law of the resolved-scale turbulence, which is affected by several factors, is far from that of the full-scale turbulence and should be corrected. These results are then applied to an existing subgrid model to improve its performance. A series of simulations are performed to verify the necessity of a fixed scaling law in the subgrid modeling.
基金Project supported by the National Basic Research Program of China(Grant Nos.2012CB720101 and 2012CB720103)the National Natural Science Foundation of China(Grant Nos.11272233,11332006,and 11411130150)
文摘When using a miniature single sensor boundary layer probe, the time sequences of the stream-wise velocity in the turbulent boundary layer （TBL） are measured by using a hot wire anemometer. Beneath the fully developed TBL, the wall pressure fluctuations are attained by a microphone mechanism with high spatial resolution. Analysis on the statistic and spectrum properties of velocity and wall pressure reveals the relationship between the wall pressure fluctuation and the energy-containing structure in the buffer layer of the TBL. Wavelet transform shows the multi-scale natures of coherent structures contained in both signals of velocity and pressure. The most intermittent wall pressure scale is associated with the coherent structure in the buffer layer. Meanwhile the most energetic scale of velocity fluctuation at y＋ = 14 provides a specific frequency f9 ≈ 147 Hz for wall actuating control with Ret = 996.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11332006 and 11272233)the National Key Basic Research and Development Program of China(Grant No.2012CB720101)
文摘The present experimental work focuses on a new model for space-time correlation and the scale-dependencies of convection velocity and sweep velocity in turbulent boundary layer over a flat wail. A turbulent boundary layer flow at Reo = 2460 is measured by tomographic particle image velocimetry （tomographic PIV）. It is demonstrated that arch, cane, and hairpin vortices are dominant in the logarithmic layer. Hairpins and hairpin packets are responsible for the elongated low-momentum zones observed in the instantaneous flow field. The conditionally-averaged coherent structures systemically illustrate the key roles of hairpin vortice in the turbulence dynamic events, such as ejection and sweep events and energy transport. The space-time correlations of instantaneous streamwise fluctuation velocity are calculated and confirm the new elliptic model for the space-time correlation instead of Taylor hypothesis. The convection velocities derived from the space-time correlation and conditionally-averaged method both suggest the scaling with the local mean velocity in the logarithmic layer. Convection velocity result based on Fourier decomposition （FD） shows stronger scale- dependency in the spanwise direction than in streamwise direction. Compared with FD, the proper orthogonal decomposition （POD） has a distinct distribution of convection velocity for the large- and small-scales which are separated in light of their contributions of turbulent kinetic energy.
基金supported by the National Natural Science Foundation of China (NSFC) Basic Science Center Program for “Multiscale Problems in Nonlinear Mechanics” (Grant No. 11988102)the NSFC Program (Grant No. 11772337)+3 种基金the Science Challenge Program (Grant No. TZ2016001)the Strategic Priority Research Program,Chinese Academy of Sciences (CAS) (Grant No. XDB22040104)the Key Research Program of Frontier Sciences, CAS (Grant No. QYZDJ-SSW-SYS002)the CAS Center for Excellence in Complex System Mechanics
文摘Turbulent flows over rough surfaces widely exist in nature and industry.Investigating its mechanism is of theoretical and practical significance.In this work we simulate the turbulent channel flow with rough walls using large-eddy simulation with rough elements resolved using the curvilinear immersed boundary method and compare the results obtained in this work with those in the paper by Yuan and Piomelli(J.Fluid Mech.,vol.760,pp.R1,2014),where the volume of fluid method was employed for modeling rough elements.The mean streamwise velocity profiles predicted by the two methods agree well with each other.Differences in Reynolds stresses and dispersive stresses are observed,which are attributed to the different approaches in dealing with the complex geometry of the rough surface.
基金supported by the National Natural Science Foundation of China(11272233)National Key Basic Research and Development Program(2012CB720101)2012 opening subjects of The State Key Laboratory of Nonlinear Mechanics(LNM),Institute of Mechanics,Chinese Academy of Sciences
文摘An experimental measurement was performed us- ing time-resolved particle image velocimetry （TRPIV） to in- vestigate the spatial topological character of coherent struc- tures in wall-bounded turbulence of polymer additive solu- tion. The fully developed near-wall turbulent flow fields with and without polymer additives at the same Reynolds number were measured by TRPIV in a water channel. The compar- isons of turbulent statistics confirm that due to viscoelastic structure of long-chain polymers, the wall-normal velocity fluctuation and Reynolds shear stress in the near-wall region are suppressed significantly. Furthermore, it is noted that such a behavior of polymers is closely related to the decease of the motion of the second and forth quadrants, i.e., the ejection and sweep events, in the near-wall region. The spa- tial topological mode of coherent structures during bursts has been extracted by the new mu-level criteria based on locally averaged velocity structure function. Although the general shapes of coherent structures are unchanged by polymer additives, the fluctuating velocity, velocity gradient, velocity strain rate and vorticity of coherent structures during burst events are suppressed in the polymer additive solution com- pared with that in water. The results show that due to the polymer additives the occurrence and intensity of coherent structures are suppressed, leading to drag reduction.
基金supported by the National Natural Science Foundation of China(NSFC)Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”(No.11988102)NSFC(No.12002345)China Postdoctoral Science Foundation(No.2020M680027)。
文摘Flow around a real-life underwater vehicle often happens at a high Reynolds number with flow structures at different scales from the boundary layer around a blade to that around the hull. This poses a great challenge for large-eddy simulation of an underwater vehicle aiming at resolving all relevant flow scales. In this work, we propose to model the hull with appendages using the immersed boundary method, and model the propeller using the actuator disk model without resolving the geometry of the blade. The proposed method is then applied to simulate the flow around Defense Advanced Research Projects Agency(DARPA) suboff. An overall acceptable agreement is obtained for the pressure and friction coefficients. Complex flow features are observed in the near wake of suboff. In the far wake, the core region is featured by a jet because of the actuator disk, surrounded by an annular region with velocity deficit due to the body of suboff.
基金supported by the National Natural Science Foundation of China (10832011)
文摘The mechanical properties of laser welded joints under impact loadings such as explosion and car crash etc. are critical for the engineering designs. The hardness, static and dynamic mechanical properties of AISI304 and AISI316L dissimilar stainless steel welded joints by CO2 laser were experimentally studied. The dynamic strainstress curves at the strain rate around 103 s-1 were obtained by the split Hopkinson tensile bar （SHTB）. The static mechanical properties of the welded joints have little changes with the laser power density and all fracture occurs at 316 L side. However, the strain rate sensitivity has a strong depen- dence on laser power density. The value of strain rate factor decreases with the increase of laser power density. The welded joint which may be applied for the impact loading can be obtained by reducing the laser power density in the case of welding quality assurance.
