The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direc...The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direction of the basic flows.By defining an energy functional,it is proven that plane parallel shear flows are unconditionally nonlinearly exponentially stable for tilted streamwise perturbation when the Reynolds number is below a certain critical value and the boundary conditions are either rigid or stress-free.In the case of stress-free boundaries,by taking advantage of the poloidal-toroidal decomposition of a solenoidal field to define energy functionals,it can be even shown that plane parallel shear flows are unconditionally nonlinearly exponentially stable for all Reynolds numbers,where the tilted perturbation can be either spanwise or streamwise.展开更多
The geometric properties of fracture surfaces significantly influence shear-seepage in rock fractures,introducing complexities to fracture modelling.The present study focuses on the hydro-mechanical behaviours of roug...The geometric properties of fracture surfaces significantly influence shear-seepage in rock fractures,introducing complexities to fracture modelling.The present study focuses on the hydro-mechanical behaviours of rough rock fractures during shear-seepage processes to reveal how dilatancy and fracture asperities affect these phenomena.To achieve this,an improved shear-flow model(SFM)is proposed with the incorporation of dilatancy effect and asperities.In particular,shear dilatancy is accounted for in both the elastic and plastic stages,in contrast to some existing models that only consider it in the elastic stage.Depending on the computation approaches for the peak dilatancy angle,three different versions of the SFM are derived based on Mohr-Coulomb,joint roughness coefficient-joint compressive strength(JRC-JCS),and Grasselli’s theories.Notably,this is a new attempt that utilizes Grasselli’s model in shearseepage analysis.An advanced parameter optimization method is introduced to accurately determine model parameters,addressing the issue of local optima inherent in some conventional methods.Then,model performance is evaluated against existing experimental results.The findings demonstrate that the SFM effectively reproduces the shear-seepage characteristics of rock fracture across a wide range of stress levels.Further sensitivity analysis reveals how dilatancy and asperity affect hydraulic properties.The relation between hydro-mechanical properties(dilatancy displacement and hydraulic conductivity)and asperity parameters is analysed.Several profound understandings of the shear-seepage process are obtained by exploring the phenomenon under various conditions.展开更多
Idealized numerical simulations have been carried out to reveal the complexity in the development of asymmetric convection in a tropical cyclone(TC)under the influence of an environment with either uniform flow,vertic...Idealized numerical simulations have been carried out to reveal the complexity in the development of asymmetric convection in a tropical cyclone(TC)under the influence of an environment with either uniform flow,vertical wind shear(VWS),or both.Results show that rainwater is enhanced to the right of the motion in the outer rainband,but such enhancement occurs in the upshear-left area of the inner-core region.Additionally,due to the asymmetries introduced by environmental flow,wavenumber-1 temperature and height anomalies develop at a radius of~1000 km in the upper levels.A sub-vortex aside from the TC center encompassing the wavenumber-1 warm center appears,and asymmetric horizontal winds emerge,which,in turn,changes the storm-scale(within 400 km)VWS.Deep convection in the inner core closely follows the changing storm-scale VWS when its magnitude is larger than 2 m s^(-1) and is located downshear of the storm-scale VWS in all the experiments with environmental flow.In the outer rainbands,the maximum boundary layer convergence is mainly controlled by the direction of motion and is located in the rear-right quadrant.These results extend upon the findings of previous studies in three aspects:(1)The discovery of the roughly linear combination effect from the uniform flow and large-scale VWS;(2)The development of upper-level asymmetric winds on a 1000-km scale through the interaction between the TC vortex and environmental flow,resulting in changes in the storm-scale VWS pattern within the TC area;(3)The revelation that TC asymmetric convection closely aligns with the direction-varying storm-scale VWS instead of the initially designated VWS.展开更多
A series of fully three-dimensional(3 D) numerical simulations of flow past a free-to-oscillate curved flexible riser in shear flow were conducted at Reynolds number of 185–1015. The numerical results obtained by the...A series of fully three-dimensional(3 D) numerical simulations of flow past a free-to-oscillate curved flexible riser in shear flow were conducted at Reynolds number of 185–1015. The numerical results obtained by the two-way fluid–structure interaction(FSI) simulations are in good agreement with the experimental results reported in the earlier study. It is further found that the frequency transition is out of phase not only in the inline(IL) and crossflow(CF) directions but also along the span direction. The mode competition leads to the non-zero nodes of the rootmean-square(RMS) amplitude and the relatively chaotic trajectories. The fluid–structure interaction is to some extent reflected by the transverse velocity of the ambient fluid, which reaches the maximum value when the riser reaches the equilibrium position. Moreover, the local maximum transverse velocities occur at the peak CF amplitudes, and the values are relatively large when the vibration is in the resonance regions. The 3 D vortex columns are shed nearly parallel to the axis of the curved flexible riser. As the local Reynolds number increases from 0 at the bottom of the riser to the maximum value at the top, the wake undergoes a transition from a two-dimensional structure to a 3 D one. More irregular small-scale vortices appeared at the wake region of the riser, undergoing large amplitude responses.展开更多
The particle dynamics in an annular shear granular flow is studied using the discrete element method, and the influences of packing fraction, shear rate and friction coefficient are analyzed. We demonstrate the existe...The particle dynamics in an annular shear granular flow is studied using the discrete element method, and the influences of packing fraction, shear rate and friction coefficient are analyzed. We demonstrate the existence of a critical packing fraction exists in the shear granular flow. When the packing fraction is lower than this critical value, the mean tangential velocity profile exhibits a rate-independent feature. However, when the packing fraction exceeds this critical value, the tangential velocity profile becomes rate-dependent and varies gradually from linear to nonlinear with increasing shear rate. Furthermore, we find a continuous transition from the unjammed state to the jammed state in a shear granular flow as the packing fraction increases. In this transforming process, the force distribution varies distinctly and the contact force network also exhibits different features.展开更多
A time domain model is presented to study the vibrations of long slender cylinders placed in shear flow. Long slender cylinders such as risers and tension legs are widely used in the field of ocean engineering. They a...A time domain model is presented to study the vibrations of long slender cylinders placed in shear flow. Long slender cylinders such as risers and tension legs are widely used in the field of ocean engineering. They are subjected to vortex-induced vibrations(VIV) when placed within a transverse incident flow. A three dimensional model coupled with wake oscillators is formulated to describe the response of the slender cylinder in cross-flow and in-line directions. The wake oscillators are distributed along the cylinder and the vortex-shedding frequency is derived from the local current velocity. A non-linear fiuid force model is accounted for the coupled effect between cross-flow and in-line vibrations. The comparisons with the published experimental data show that the dynamic features of VIV of long slender cylinder placed in shear flow can be obtained by the proposed model,such as the spanwise average displacement,vibration frequency,dominant mode and the combination of standing and traveling waves. The simulation in a uniform flow is also conducted and the result is compared with the case of nonuniform flow. It is concluded that the flow shear characteristic has significantly changed the cylinder vibration behavior.展开更多
An analytical solution of the governing equations of the interacting shear flows for unsteady oblique stagnation point flow is obtained. It has the same form as that of the exact solution obtained from the complete NS...An analytical solution of the governing equations of the interacting shear flows for unsteady oblique stagnation point flow is obtained. It has the same form as that of the exact solution obtained from the complete NS equations and physical analysis and relevant discussions are then presented.展开更多
We compare the space-time correlations calculated from direct numerical simulation (DNS) and large-eddy simulation (LES) of turbulent channel flows. It is found from the comparisons that the LES with an eddy-visco...We compare the space-time correlations calculated from direct numerical simulation (DNS) and large-eddy simulation (LES) of turbulent channel flows. It is found from the comparisons that the LES with an eddy-viscosity subgrid scale (SGS) model over-predicts the space-time corre- lations than the DNS. The overpredictions are further quantified by the integral scales of directional correlations and convection velocities. A physical argument for the overpre- diction is provided that the eddy-viscosity SGS model alone does not includes the backscatter effects although it correctly represents the energy dissipations of SGS motions. This argument is confirmed by the recently developed elliptic model for space-time correlations in turbulent shear flows. It suggests that enstrophy is crucial to the LES prediction of spacetime correlations. The random forcing models and stochastic SGS models are proposed to overcome the overpredictions on space-time correlations.展开更多
A numerical study based on a wake oscillator model was conducted to determine the response performance of vortex-induced vibration(VIV) on a long flexible cylinder with pinned-pinned boundary conditions subjected to l...A numerical study based on a wake oscillator model was conducted to determine the response performance of vortex-induced vibration(VIV) on a long flexible cylinder with pinned-pinned boundary conditions subjected to linear and exponential shear flows. The coupling equations of a structural vibration model and wake oscillator model were solved using a standard central finite difference method of the second order. The VIV response characteristics including the structural displacement, structural frequency, structural wavenumber, standing wave behavior,travelling wave behavior, structural velocity, lift force coefficient and transferred energy from the fluid to the structure with different flow profiles were compared. The numerical results show that the VIV displacement is a combination of standing waves and travelling waves. For linear shear flow, standing waves and travelling waves dominate the VIV response within the low-velocity and high-velocity zones, respectively. The negative values of the transferred energy only occur within the low-velocity zone. However, for exponential shear flow, travelling waves dominate the VIV response and the negative energy occurs along the entire length of the cylinder.展开更多
This paper studies some interesting features of two-dimensional granular shearing flow by using molecular dynamic approach for a specific granular system. The obtained results show that the probability distribution fu...This paper studies some interesting features of two-dimensional granular shearing flow by using molecular dynamic approach for a specific granular system. The obtained results show that the probability distribution function of velocities of particles is Gaussian at the central part, but diverts from Gaussian distribution nearby the wall. The macroscopic stress along the vertical direction has large fluctuation around a constant value, the non-zero average velocity occurs mainly near the moving wall, which forms a shearing zone.. In the shearing movement, the volume of the granular material behaves in a random manner. The equivalent fl'iction coefficient between moving slab and granular material correlates with the moving speed at low velocity, and approaches constant as the velocity is large enough.展开更多
This work presents a numerical investigation on steady internal, external and surface flows of a liquid sphere immersed in a simple shear flow at low and intermediate Reynolds numbers. The control volume formulation i...This work presents a numerical investigation on steady internal, external and surface flows of a liquid sphere immersed in a simple shear flow at low and intermediate Reynolds numbers. The control volume formulation is adopted to solve the governing equations of two-phase flow in a 3-D spherical coordinate system. Numerical results show that the streamlines for Re = 0 are closed Jeffery orbits on the surface of the liquid sphere, and also closed curves outside and inside the liquid sphere. However, the streamlines have intricate and non-closed structures for Re ≠ 0. The flow structure is dependent on the values of Reynolds number and interior-to-exterior viscosity ratio.展开更多
The experimental results of the deformation and breakup of a single drop immersed in a Newtonian liq-uid and subjected to a constant shear rate which generated by counter rotating Couette apparatus were presented in t...The experimental results of the deformation and breakup of a single drop immersed in a Newtonian liq-uid and subjected to a constant shear rate which generated by counter rotating Couette apparatus were presented in this paper. From experimental observations, the breakup occurred by three mechanisms, namely, necking, end pinching, and capillary instability. Quantitative results for the deformation and breakup of drop are presented. The maximum diameter and Sauter mean diameter of daughter drops and capillary thread radius are linearly related to the inverse shear rate and independent of the initial drop size, the dimensionless wavelength which is the wave-length divided by the thread width at breakup is independent of the shear rate and initial drop size, and the deforma-tion of threads follows a pseudo-affine deformation for Cai/Cac larger than 2.展开更多
The coherent structures and the chaotic phenomena in the transition of the axisymmetric countercurrent mixing shear flow were investigated experimentally. Two kinds of self-excited oscillation modes could exist in the...The coherent structures and the chaotic phenomena in the transition of the axisymmetric countercurrent mixing shear flow were investigated experimentally. Two kinds of self-excited oscillation modes could exist in the axisymmetric countercurrent mixing shear flow. One is the shear layer self-excited oscillation mode corresponding to the high Reynolds number regime and the other is the jet column self-excited oscillation mode corresponding to the low Reynolds number regime in the case of the velocity ratio ranging from I to 1.5. Analyzing the auto-power spectrum, self-correlation-function and three dimensional reconstructed phase trajectory, the route to chaos through three Hopf bifurcations intercepted by an intermittence of the dynamical system corresponding to the axisymmetric countercurrent mixing shear flow was discovered when the velocity ratio is equal to 1.32.展开更多
The flow of the weak electrolyte solution can be controlled by Lorentz force achieved with the suitable magnetic and electric fields, and it has the advantages of vortex street suppression, drag reduction, lift enhanc...The flow of the weak electrolyte solution can be controlled by Lorentz force achieved with the suitable magnetic and electric fields, and it has the advantages of vortex street suppression, drag reduction, lift enhancement and oscillatory suppression for the flow over a bluff body. The electro-magnetic control of vortex-induced vibration (VIV) of a circular cylinder in the shear flow was investigated numerically in the exponential-polar coordinates attached on the moving cylinder for Re=150. With the effect of background vorticity, the vortex street of VIV cylinder was composed of two parallel rows with an opposite sign of the vortices which inclines toward the lower side and the strength of upper vortex is larger than that of lower vortex. The lift force vibrated periodically with the effect of vortex shedding and the mean value was negative due to the background vorticity. The Lorentz force for controlling the VIV cylinder was classified into the field Lorentz force and the wall Lorentz force. The field Lorentz force suppresses the lift oscillation, and in turn, suppresses the VIV, whereas the wall Lorentz force increases the lift.展开更多
The structure and dynamics of confined single polymer chain in a dilute solution, either in equilibrium or at different shear rates in the uniform shear flow fields, were investigated by means of dissipative particle ...The structure and dynamics of confined single polymer chain in a dilute solution, either in equilibrium or at different shear rates in the uniform shear flow fields, were investigated by means of dissipative particle dynamics simulations. The no-slip boundary condition without density fluctuation near the wall was taken into account to mimic the environment of a nanochannel. The dependences of the radius of gyration, especially in three different di- rections, and the density profile of the chain mass center on the strength of the confinement and the Weissenberg number(Wn) was studied. The effect of the interaction between polymer and solvent on the density profile was also investigated in the cases of moderate and strong Wn. In the high shear flow, the polymer migrates to the center of the channel with increasing Wn. There is only one density profile peak in the channel center in the uniform shear flow, which is in agreement with the results of the experiments and theory.展开更多
Within the magnetohydrodynamics (MHD) frame, we analyse the effect of viscosity on magneto-Rayleigh Taylor (MRT) instability in a Z-pinch configuration by using an exact method and an approximate method separately...Within the magnetohydrodynamics (MHD) frame, we analyse the effect of viscosity on magneto-Rayleigh Taylor (MRT) instability in a Z-pinch configuration by using an exact method and an approximate method separately. It is demonstrated that the plasma viscosity indeed has a stabilization effect on the MRT mode in the whole wavenumber region, and its influence increases with the perturbation wavenumber increasing. After the characteristics and feasibility of the approximate method have been investigated, we apply it to the stability analysis of viscous plasma where a sheared axial flow (SAF) is involved, and we attain an analytical dispersion relation. It is suggested that the viscosity and the SAF are complemental with each other, and a wide wavenumber range of perturbation is possible to be restrained if the SAF and the viscosity are large enough. Finally, we calculate the possible value of viscosity parameter according to the current experimental conditions, and the results show that since the value of viscosity is much less than the threshold value, its mitigation effect is small enough to be neglected. The role of the viscosity in the stabilization becomes considerable only if special techniques are so developed that the Z-pinch plasma viscosity can be increased greatly.展开更多
This study experimentally analyzes the nonlinear flow characteristics and channelization of fluid through rough-walled fractures during the shear process using a shear-flow-visualization apparatus.A series of fluid fl...This study experimentally analyzes the nonlinear flow characteristics and channelization of fluid through rough-walled fractures during the shear process using a shear-flow-visualization apparatus.A series of fluid flow and visualization tests is performed on four transparent fracture specimens with various shear displacements of 1 mm,3 mm,5 mm,7 mm and 10 mm under a normal stress of 0.5 MPa.Four granite fractures with different roughnesses are selected and quantified using variogram fractal dimensions.The obtained results show that the critical Reynolds number tends to increase with increasing shear displacement but decrease with increasing roughness of fracture surface.The flow paths are more tortuous at the beginning of shear because of the wide distribution of small contact spots.As the shear displacement continues to increase,preferential flow paths are more distinctly observed due to the decrease in the number of contact spots caused by shear dilation;yet the area of single contacts in-creases.Based on the experimental results,an empirical mathematical equation is proposed to quantify the critical Reynolds number using the contact area ratio and fractal dimension.展开更多
Considering that the fluid is inviscid and incompressible and the flow is irrotational in a fixed frame of reference and using the multiple scale analysis method, we derive a nonlinear Schrodinger equation(NLSE) descr...Considering that the fluid is inviscid and incompressible and the flow is irrotational in a fixed frame of reference and using the multiple scale analysis method, we derive a nonlinear Schrodinger equation(NLSE) describing the evolution dynamics of gravity-capillary wavetrains in arbitrary constant depth. The gravity-capillary waves(GCWs) are influenced by a linear shear flow(LSF) which consists of a uniform flow and a shear flow with constant vorticity. The modulational instability(MI) of GCWs with the LSF is analyzed using the NLSE. The MI is effectively modified by the LSF. In infinite depth, there are four asymptotes which are the boundaries between MI and modulational stability(MS) in the instability diagram. In addition, the dimensionless free surface elevation as a function of time for different dimensionless water depth,surface tension, uniform flow and vorticity is exhibited. It is found that the decay of free surface elevation and the steepness of free surface amplitude change over time, which are greatly affected by the water depth, surface tension, uniform flow and vorticity.展开更多
基金supported by the National Natural Science Foundation of China(21627813)。
文摘The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direction of the basic flows.By defining an energy functional,it is proven that plane parallel shear flows are unconditionally nonlinearly exponentially stable for tilted streamwise perturbation when the Reynolds number is below a certain critical value and the boundary conditions are either rigid or stress-free.In the case of stress-free boundaries,by taking advantage of the poloidal-toroidal decomposition of a solenoidal field to define energy functionals,it can be even shown that plane parallel shear flows are unconditionally nonlinearly exponentially stable for all Reynolds numbers,where the tilted perturbation can be either spanwise or streamwise.
基金support from the National Natural Science Foundation of China(Grant Nos.51991392 and 42293355).
文摘The geometric properties of fracture surfaces significantly influence shear-seepage in rock fractures,introducing complexities to fracture modelling.The present study focuses on the hydro-mechanical behaviours of rough rock fractures during shear-seepage processes to reveal how dilatancy and fracture asperities affect these phenomena.To achieve this,an improved shear-flow model(SFM)is proposed with the incorporation of dilatancy effect and asperities.In particular,shear dilatancy is accounted for in both the elastic and plastic stages,in contrast to some existing models that only consider it in the elastic stage.Depending on the computation approaches for the peak dilatancy angle,three different versions of the SFM are derived based on Mohr-Coulomb,joint roughness coefficient-joint compressive strength(JRC-JCS),and Grasselli’s theories.Notably,this is a new attempt that utilizes Grasselli’s model in shearseepage analysis.An advanced parameter optimization method is introduced to accurately determine model parameters,addressing the issue of local optima inherent in some conventional methods.Then,model performance is evaluated against existing experimental results.The findings demonstrate that the SFM effectively reproduces the shear-seepage characteristics of rock fracture across a wide range of stress levels.Further sensitivity analysis reveals how dilatancy and asperity affect hydraulic properties.The relation between hydro-mechanical properties(dilatancy displacement and hydraulic conductivity)and asperity parameters is analysed.Several profound understandings of the shear-seepage process are obtained by exploring the phenomenon under various conditions.
基金supported by the National Natural Science Foundation of China(Grant number 42075072)support from The Startup Foundation for Introducing Talent of the Nanjing University of Information Science and Technology.
文摘Idealized numerical simulations have been carried out to reveal the complexity in the development of asymmetric convection in a tropical cyclone(TC)under the influence of an environment with either uniform flow,vertical wind shear(VWS),or both.Results show that rainwater is enhanced to the right of the motion in the outer rainband,but such enhancement occurs in the upshear-left area of the inner-core region.Additionally,due to the asymmetries introduced by environmental flow,wavenumber-1 temperature and height anomalies develop at a radius of~1000 km in the upper levels.A sub-vortex aside from the TC center encompassing the wavenumber-1 warm center appears,and asymmetric horizontal winds emerge,which,in turn,changes the storm-scale(within 400 km)VWS.Deep convection in the inner core closely follows the changing storm-scale VWS when its magnitude is larger than 2 m s^(-1) and is located downshear of the storm-scale VWS in all the experiments with environmental flow.In the outer rainbands,the maximum boundary layer convergence is mainly controlled by the direction of motion and is located in the rear-right quadrant.These results extend upon the findings of previous studies in three aspects:(1)The discovery of the roughly linear combination effect from the uniform flow and large-scale VWS;(2)The development of upper-level asymmetric winds on a 1000-km scale through the interaction between the TC vortex and environmental flow,resulting in changes in the storm-scale VWS pattern within the TC area;(3)The revelation that TC asymmetric convection closely aligns with the direction-varying storm-scale VWS instead of the initially designated VWS.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.11502220 and51479126)the Youth Science and Technology Foundation of Sichuan Province(Grant No.2017JQ0055)the Youth Scientific and Technological Innovation Team of the Safety of Deep-Water Pipe Strings of Southwest Petroleum University(Grant No.2017CXTD06)
文摘A series of fully three-dimensional(3 D) numerical simulations of flow past a free-to-oscillate curved flexible riser in shear flow were conducted at Reynolds number of 185–1015. The numerical results obtained by the two-way fluid–structure interaction(FSI) simulations are in good agreement with the experimental results reported in the earlier study. It is further found that the frequency transition is out of phase not only in the inline(IL) and crossflow(CF) directions but also along the span direction. The mode competition leads to the non-zero nodes of the rootmean-square(RMS) amplitude and the relatively chaotic trajectories. The fluid–structure interaction is to some extent reflected by the transverse velocity of the ambient fluid, which reaches the maximum value when the riser reaches the equilibrium position. Moreover, the local maximum transverse velocities occur at the peak CF amplitudes, and the values are relatively large when the vibration is in the resonance regions. The 3 D vortex columns are shed nearly parallel to the axis of the curved flexible riser. As the local Reynolds number increases from 0 at the bottom of the riser to the maximum value at the top, the wake undergoes a transition from a two-dimensional structure to a 3 D one. More irregular small-scale vortices appeared at the wake region of the riser, undergoing large amplitude responses.
基金supported by the Key Project of the National Natural Science Foundation of China(10532040)the Programme of Changjiang Scholars and Innovative Research Team in University (IRT 0628)
文摘The particle dynamics in an annular shear granular flow is studied using the discrete element method, and the influences of packing fraction, shear rate and friction coefficient are analyzed. We demonstrate the existence of a critical packing fraction exists in the shear granular flow. When the packing fraction is lower than this critical value, the mean tangential velocity profile exhibits a rate-independent feature. However, when the packing fraction exceeds this critical value, the tangential velocity profile becomes rate-dependent and varies gradually from linear to nonlinear with increasing shear rate. Furthermore, we find a continuous transition from the unjammed state to the jammed state in a shear granular flow as the packing fraction increases. In this transforming process, the force distribution varies distinctly and the contact force network also exhibits different features.
基金supported by the National Natural Science Foundation of China (10532070)the Knowledge Innovation Program of Chinese Academy of Sciences (KJCX2-YW-L07)the LNM Initial Funding for Young Investigators
文摘A time domain model is presented to study the vibrations of long slender cylinders placed in shear flow. Long slender cylinders such as risers and tension legs are widely used in the field of ocean engineering. They are subjected to vortex-induced vibrations(VIV) when placed within a transverse incident flow. A three dimensional model coupled with wake oscillators is formulated to describe the response of the slender cylinder in cross-flow and in-line directions. The wake oscillators are distributed along the cylinder and the vortex-shedding frequency is derived from the local current velocity. A non-linear fiuid force model is accounted for the coupled effect between cross-flow and in-line vibrations. The comparisons with the published experimental data show that the dynamic features of VIV of long slender cylinder placed in shear flow can be obtained by the proposed model,such as the spanwise average displacement,vibration frequency,dominant mode and the combination of standing and traveling waves. The simulation in a uniform flow is also conducted and the result is compared with the case of nonuniform flow. It is concluded that the flow shear characteristic has significantly changed the cylinder vibration behavior.
文摘An analytical solution of the governing equations of the interacting shear flows for unsteady oblique stagnation point flow is obtained. It has the same form as that of the exact solution obtained from the complete NS equations and physical analysis and relevant discussions are then presented.
基金supported by the National Basic Research Program of China (973 Program) (2007CB814800)the National Natural Science Foundation of China (10325211 and 10628206)
文摘We compare the space-time correlations calculated from direct numerical simulation (DNS) and large-eddy simulation (LES) of turbulent channel flows. It is found from the comparisons that the LES with an eddy-viscosity subgrid scale (SGS) model over-predicts the space-time corre- lations than the DNS. The overpredictions are further quantified by the integral scales of directional correlations and convection velocities. A physical argument for the overpre- diction is provided that the eddy-viscosity SGS model alone does not includes the backscatter effects although it correctly represents the energy dissipations of SGS motions. This argument is confirmed by the recently developed elliptic model for space-time correlations in turbulent shear flows. It suggests that enstrophy is crucial to the LES prediction of spacetime correlations. The random forcing models and stochastic SGS models are proposed to overcome the overpredictions on space-time correlations.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51609206,51522902 and 51579040)
文摘A numerical study based on a wake oscillator model was conducted to determine the response performance of vortex-induced vibration(VIV) on a long flexible cylinder with pinned-pinned boundary conditions subjected to linear and exponential shear flows. The coupling equations of a structural vibration model and wake oscillator model were solved using a standard central finite difference method of the second order. The VIV response characteristics including the structural displacement, structural frequency, structural wavenumber, standing wave behavior,travelling wave behavior, structural velocity, lift force coefficient and transferred energy from the fluid to the structure with different flow profiles were compared. The numerical results show that the VIV displacement is a combination of standing waves and travelling waves. For linear shear flow, standing waves and travelling waves dominate the VIV response within the low-velocity and high-velocity zones, respectively. The negative values of the transferred energy only occur within the low-velocity zone. However, for exponential shear flow, travelling waves dominate the VIV response and the negative energy occurs along the entire length of the cylinder.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10872005 and 10128204)
文摘This paper studies some interesting features of two-dimensional granular shearing flow by using molecular dynamic approach for a specific granular system. The obtained results show that the probability distribution function of velocities of particles is Gaussian at the central part, but diverts from Gaussian distribution nearby the wall. The macroscopic stress along the vertical direction has large fluctuation around a constant value, the non-zero average velocity occurs mainly near the moving wall, which forms a shearing zone.. In the shearing movement, the volume of the granular material behaves in a random manner. The equivalent fl'iction coefficient between moving slab and granular material correlates with the moving speed at low velocity, and approaches constant as the velocity is large enough.
基金Supported by the National Basic Research Program of China(2013CB632601)the National Science Fund for Distinguished Young Scholars(21025627)+1 种基金the National Natural Science Foundation of China(21276256,21106150)the National High Technology Research and Development Program of China(2012AA03A606)
文摘This work presents a numerical investigation on steady internal, external and surface flows of a liquid sphere immersed in a simple shear flow at low and intermediate Reynolds numbers. The control volume formulation is adopted to solve the governing equations of two-phase flow in a 3-D spherical coordinate system. Numerical results show that the streamlines for Re = 0 are closed Jeffery orbits on the surface of the liquid sphere, and also closed curves outside and inside the liquid sphere. However, the streamlines have intricate and non-closed structures for Re ≠ 0. The flow structure is dependent on the values of Reynolds number and interior-to-exterior viscosity ratio.
基金Supported by the National Natural Science Foundation of China (No.50536020 and No.10172069).
文摘The experimental results of the deformation and breakup of a single drop immersed in a Newtonian liq-uid and subjected to a constant shear rate which generated by counter rotating Couette apparatus were presented in this paper. From experimental observations, the breakup occurred by three mechanisms, namely, necking, end pinching, and capillary instability. Quantitative results for the deformation and breakup of drop are presented. The maximum diameter and Sauter mean diameter of daughter drops and capillary thread radius are linearly related to the inverse shear rate and independent of the initial drop size, the dimensionless wavelength which is the wave-length divided by the thread width at breakup is independent of the shear rate and initial drop size, and the deforma-tion of threads follows a pseudo-affine deformation for Cai/Cac larger than 2.
文摘The coherent structures and the chaotic phenomena in the transition of the axisymmetric countercurrent mixing shear flow were investigated experimentally. Two kinds of self-excited oscillation modes could exist in the axisymmetric countercurrent mixing shear flow. One is the shear layer self-excited oscillation mode corresponding to the high Reynolds number regime and the other is the jet column self-excited oscillation mode corresponding to the low Reynolds number regime in the case of the velocity ratio ranging from I to 1.5. Analyzing the auto-power spectrum, self-correlation-function and three dimensional reconstructed phase trajectory, the route to chaos through three Hopf bifurcations intercepted by an intermittence of the dynamical system corresponding to the axisymmetric countercurrent mixing shear flow was discovered when the velocity ratio is equal to 1.32.
基金Sponsored by the National Nature Science Foundation of China (11202102,11172140)
文摘The flow of the weak electrolyte solution can be controlled by Lorentz force achieved with the suitable magnetic and electric fields, and it has the advantages of vortex street suppression, drag reduction, lift enhancement and oscillatory suppression for the flow over a bluff body. The electro-magnetic control of vortex-induced vibration (VIV) of a circular cylinder in the shear flow was investigated numerically in the exponential-polar coordinates attached on the moving cylinder for Re=150. With the effect of background vorticity, the vortex street of VIV cylinder was composed of two parallel rows with an opposite sign of the vortices which inclines toward the lower side and the strength of upper vortex is larger than that of lower vortex. The lift force vibrated periodically with the effect of vortex shedding and the mean value was negative due to the background vorticity. The Lorentz force for controlling the VIV cylinder was classified into the field Lorentz force and the wall Lorentz force. The field Lorentz force suppresses the lift oscillation, and in turn, suppresses the VIV, whereas the wall Lorentz force increases the lift.
基金Supported by the National Natural Science Foundation of China(No.20774036)Fok Ying Tung Education Foundation (No.114008)
文摘The structure and dynamics of confined single polymer chain in a dilute solution, either in equilibrium or at different shear rates in the uniform shear flow fields, were investigated by means of dissipative particle dynamics simulations. The no-slip boundary condition without density fluctuation near the wall was taken into account to mimic the environment of a nanochannel. The dependences of the radius of gyration, especially in three different di- rections, and the density profile of the chain mass center on the strength of the confinement and the Weissenberg number(Wn) was studied. The effect of the interaction between polymer and solvent on the density profile was also investigated in the cases of moderate and strong Wn. In the high shear flow, the polymer migrates to the center of the channel with increasing Wn. There is only one density profile peak in the channel center in the uniform shear flow, which is in agreement with the results of the experiments and theory.
基金Project supported by the National Science Foundation of China (Grant Nos 10575014 and 10635050)
文摘Within the magnetohydrodynamics (MHD) frame, we analyse the effect of viscosity on magneto-Rayleigh Taylor (MRT) instability in a Z-pinch configuration by using an exact method and an approximate method separately. It is demonstrated that the plasma viscosity indeed has a stabilization effect on the MRT mode in the whole wavenumber region, and its influence increases with the perturbation wavenumber increasing. After the characteristics and feasibility of the approximate method have been investigated, we apply it to the stability analysis of viscous plasma where a sheared axial flow (SAF) is involved, and we attain an analytical dispersion relation. It is suggested that the viscosity and the SAF are complemental with each other, and a wide wavenumber range of perturbation is possible to be restrained if the SAF and the viscosity are large enough. Finally, we calculate the possible value of viscosity parameter according to the current experimental conditions, and the results show that since the value of viscosity is much less than the threshold value, its mitigation effect is small enough to be neglected. The role of the viscosity in the stabilization becomes considerable only if special techniques are so developed that the Z-pinch plasma viscosity can be increased greatly.
基金This study has been partially funded by National Key Research and Development Program of China(Grant No.2020YFA0711800)the National Natural Science Foundation of China(Grant No.51979272)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2021QE069).
文摘This study experimentally analyzes the nonlinear flow characteristics and channelization of fluid through rough-walled fractures during the shear process using a shear-flow-visualization apparatus.A series of fluid flow and visualization tests is performed on four transparent fracture specimens with various shear displacements of 1 mm,3 mm,5 mm,7 mm and 10 mm under a normal stress of 0.5 MPa.Four granite fractures with different roughnesses are selected and quantified using variogram fractal dimensions.The obtained results show that the critical Reynolds number tends to increase with increasing shear displacement but decrease with increasing roughness of fracture surface.The flow paths are more tortuous at the beginning of shear because of the wide distribution of small contact spots.As the shear displacement continues to increase,preferential flow paths are more distinctly observed due to the decrease in the number of contact spots caused by shear dilation;yet the area of single contacts in-creases.Based on the experimental results,an empirical mathematical equation is proposed to quantify the critical Reynolds number using the contact area ratio and fractal dimension.
基金Project supported by the National Natural Science Foundation of China(Grant No.41830533)the National Key Research and Development Program of China(Grant Nos.2016YFC1401404 and 2017YFA0604102).
文摘Considering that the fluid is inviscid and incompressible and the flow is irrotational in a fixed frame of reference and using the multiple scale analysis method, we derive a nonlinear Schrodinger equation(NLSE) describing the evolution dynamics of gravity-capillary wavetrains in arbitrary constant depth. The gravity-capillary waves(GCWs) are influenced by a linear shear flow(LSF) which consists of a uniform flow and a shear flow with constant vorticity. The modulational instability(MI) of GCWs with the LSF is analyzed using the NLSE. The MI is effectively modified by the LSF. In infinite depth, there are four asymptotes which are the boundaries between MI and modulational stability(MS) in the instability diagram. In addition, the dimensionless free surface elevation as a function of time for different dimensionless water depth,surface tension, uniform flow and vorticity is exhibited. It is found that the decay of free surface elevation and the steepness of free surface amplitude change over time, which are greatly affected by the water depth, surface tension, uniform flow and vorticity.