As emerging two-dimensional(2D)materials,carbides and nitrides(MXenes)could be solid solutions or organized structures made up of multi-atomic layers.With remarkable and adjustable electrical,optical,mechanical,and el...As emerging two-dimensional(2D)materials,carbides and nitrides(MXenes)could be solid solutions or organized structures made up of multi-atomic layers.With remarkable and adjustable electrical,optical,mechanical,and electrochemical characteristics,MXenes have shown great potential in brain-inspired neuromorphic computing electronics,including neuromorphic gas sensors,pressure sensors and photodetectors.This paper provides a forward-looking review of the research progress regarding MXenes in the neuromorphic sensing domain and discussed the critical challenges that need to be resolved.Key bottlenecks such as insufficient long-term stability under environmental exposure,high costs,scalability limitations in large-scale production,and mechanical mismatch in wearable integration hinder their practical deployment.Furthermore,unresolved issues like interfacial compatibility in heterostructures and energy inefficiency in neu-romorphic signal conversion demand urgent attention.The review offers insights into future research directions enhance the fundamental understanding of MXene properties and promote further integration into neuromorphic computing applications through the convergence with various emerging technologies.展开更多
The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller （HALLER, G. Exact theor...The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller （HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257-311 （2004））. By analyzing the distribution of the finite-time Lyapunov exponent （FTLE） along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE. This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.展开更多
Owing to their rolling friction,two-dimensional piston pumps are highly suitable as power components for electro-hydrostatic actuators(EHAs).These pumps are particularly advantageous for applications requiring high ef...Owing to their rolling friction,two-dimensional piston pumps are highly suitable as power components for electro-hydrostatic actuators(EHAs).These pumps are particularly advantageous for applications requiring high efficiency and reliability.However,the ambiguity surrounding the output flow characteristics of individual two-dimensional pumps poses a significant challenge in achieving precise closed-loop control of the EHA positions.To address this issue,this study established a comprehensive numerical model that included gap leakage to analyze the impact of leakage on the output flow characteristics of a two-dimensional piston pump.The validity of the numerical analysis was indirectly confirmed through meticulous measurements of the leakage and volumetric efficiency,ensuring robust results.The research findings indicated that,at lower pump speeds,leakage significantly affected the output flow rate,leading to potential inefficiencies in the system.Conversely,at higher rotational speeds,the impact of leakage was less pronounced,implying that the influence of leakage on the pump outlet flow must be carefully considered and managed for EHAs to perform position servo control.Additionally,the research demonstrates that two-dimensional motion does not have a unique or additional effect on pump leakage,thus simplifying the design considerations.Finally,the study concluded that maintaining an oil-filled leakage environment is beneficial because it helps reduce the impact of leakage and enhances the overall volumetric efficiency of the pump system.展开更多
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.展开更多
Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional （2D） multi-fluid flows, as well as ...Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional （2D） multi-fluid flows, as well as construct the continuous Roe matrix for the whole flow field. The interface capture equations and fluid dynamic conservative equations are coupled together and solved by using any high-resolution schemes that usually suit for the single-fluid flows. Some numerical examples are given to illustrate the solution of 1D and 2D multi-fluid Riemann problems.展开更多
If the vertically-mixing time is comparable with the period of oscillatory current, the contaminant contraction may occur, then the coefficient of shear dispersion will be negative (singularity). According to the two-...If the vertically-mixing time is comparable with the period of oscillatory current, the contaminant contraction may occur, then the coefficient of shear dispersion will be negative (singularity). According to the two-dimensional delay-diffusion equation derived by Liu in another paper, The equations for X(t), Y(t), central displacement as well as dispersion tensor were derived. Generally speaking, δDij/δτ is positive when τ is small; but when τ is large, the memory functions may be negative. The analytical expressions for Dij and X,Y were also given.展开更多
Debris flows are recurrent natural hazards in many mountainous regions.This paper presents a numerical study on the propagation of debris flows in natural erodible open channels,in which the bed erosion and sedimentat...Debris flows are recurrent natural hazards in many mountainous regions.This paper presents a numerical study on the propagation of debris flows in natural erodible open channels,in which the bed erosion and sedimentation processes are important.Based on the Bingham fluid theory,a mathematical model of the two-dimensional non-constant debris flow is developed.The governing equations include the continuity and momentum conservation equations of debris flow,the sediment convection-diffusion equation,the bed erosion-deposition equation and the bed-sediment size gradation adjustment equation.The yield stress and shear stress components are included to describe the dynamic rheological properties.The upwind control-volume Finite Volume Method (FVM) is applied to discretize the convection terms.The improved SIMPLE algorithm with velocity-free-surface coupled correction is developed to solve the equations on non-orthogonal,quadrilateral grids.The model is applied to simulate a debris flow event in Jiangjia Gully,Yunnan Province and to predict the flow pattern and bed erosion-deposition processes.The results show the effectiveness of the proposed numercial model in debris flow simulation and potential hazard analysis.展开更多
A graphics processing unit(GPU)-accelerated discontinuous Galerkin(DG)method is presented for solving two-dimensional laminar flows.The DG method is ported from central processing unit to GPU in a way of achieving GPU...A graphics processing unit(GPU)-accelerated discontinuous Galerkin(DG)method is presented for solving two-dimensional laminar flows.The DG method is ported from central processing unit to GPU in a way of achieving GPU speedup through programming under the compute unified device architecture(CUDA)model.The CUDA kernel subroutines are designed to meet with the requirement of high order computing of DG method.The corresponding data structures are constructed in component-wised manners and the thread hierarchy is manipulated in cell-wised or edge-wised manners associated with related integrals involved in solving laminar Navier-Stokes equations,in which the inviscid and viscous flux terms are computed by the local lax-Friedrichs scheme and the second scheme of Bassi&Rebay,respectively.A strong stability preserving Runge-Kutta scheme is then used for time marching of numerical solutions.The resulting GPU-accelerated DG method is first validated by the traditional Couette flow problems with different mesh sizes associated with different orders of approximation,which shows that the orders of convergence,as expected,can be achieved.The numerical simulations of the typical flows over a circular cylinder or a NACA 0012 airfoil are then carried out,and the results are further compared with the analytical solutions or available experimental and numerical values reported in the literature,as well as with a performance analysis of the developed code in terms of GPU speedups.This shows that the costs of computing time of the presented test cases are significantly reduced without losing accuracy,while impressive speedups up to 69.7 times are achieved by the present method in comparison to its CPU counterpart.展开更多
The regions with shear stress and mean velocity gradient of opposite sign often exist in complex turbulent shear flows.In these cases,the eddy viscosity hypothesis breaks down.Hinze regards the,departure from eddy vis...The regions with shear stress and mean velocity gradient of opposite sign often exist in complex turbulent shear flows.In these cases,the eddy viscosity hypothesis breaks down.Hinze regards the,departure from eddy viscosity hypothesis as a result from transportation of mean momentum over distance by the large structures and arrives at a shear stress expression including the second order derivatives of the mean velocity.However,his expression greatly overestimates the shear stress.This implies that the flow particles are unlikely to have enough memory of the mean momentum over distance.By assuming the departure from eddy viscosity hypothesis as a result from transportation of the shear stress contained in smaller eddies over distance by the large structures,the present author has arrived at a new shear stress expression.The shear stress estimated so far is in good agreement with the experiments.展开更多
Stably stratified flows over a two-dimensional hill are investigated in a channel of finite depth using a three-dimensional direct numerical simulation (DNS). The present study follows onto our previous two-dimensiona...Stably stratified flows over a two-dimensional hill are investigated in a channel of finite depth using a three-dimensional direct numerical simulation (DNS). The present study follows onto our previous two-dimensional DNS studies of stably stratified flows over a hill in a channel of finite depth and provides a more realistic simulation of atmospheric flows than our previous studies. A hill with a constant cross-section in the spanwise (y) direction is placed in a 3-D computational domain. As in the previous 2-D simulations, to avoid the effect of the ground boundary layer that develops upstream of the hill, no-slip conditions are imposed only on the hill surface and the surface downstream of the hill;slip conditions are imposed on the surface upstream of the hill. The simulated 3-D flows are discussed by comparing them to the simulated 2-D flows with a focus on the effect of the stable stratification on the non-periodic separation and reattachment of the flow behind the hill. In neutral (K = 0, where K is a non-dimensional stability parameter) and weakly stable (K = 0.8) conditions, 3-D flows over a hill differ clearly from 2-D flows over a hill mainly because of the three-dimensionality of the flow, that is the development of a spanwise flow component in the 3-D flows. In highly stable conditions (K = 1, 1.3), long-wavelength lee waves develop downstream of the hill in both 2-D and 3-D flows, and the behaviors of the 2-D and 3-D flows are similar in the vicinity of the hill. In other words, the spanwise component of the 3-D flows is strongly suppressed in highly stable conditions, and the flow in the vicinity of the hill becomes approximately two-dimensional in the x and z directions.展开更多
Boundary conditions （BCs） play an essential role in lattice Boltzmann （LB） simulations. This paper investigates several most commonly applied BCs by evaluating the relative L2-norm errors of the LB simulations for...Boundary conditions （BCs） play an essential role in lattice Boltzmann （LB） simulations. This paper investigates several most commonly applied BCs by evaluating the relative L2-norm errors of the LB simulations for two-dimensional （2-D） Poiseuille flow. It is found that the relative L2-norm error resulting from FHML＇s BC is smaller than that from other BCs as a whole. Then, based on the FHML＇s BC, it formulates an LB model for simulating fluid flows in 2-D channel with complex geometries. Afterwards, the flows between two inclined plates, in a pulmonary blood vessel and in a blood vessel with local expansion region, are simulated. The numerical results are in good agreement with the analytical predictions and clearly show that the model is effective. It is expected that the model can be extended to simulate some real biologic flows, such as blood flows in arteries, vessels with stenosises, aneurysms and bifurcations,展开更多
Several kinds of explicit and implicit finite-difference schemes directly solving the discretized velocity distribution functions are designed with precision of different orders by analyzing the inner characteristics ...Several kinds of explicit and implicit finite-difference schemes directly solving the discretized velocity distribution functions are designed with precision of different orders by analyzing the inner characteristics of the gas-kinetic numerical algorithm for Boltzmann model equation. The peculiar flow phenomena and mechanism from various flow regimes are revealed in the numerical simulations of the unsteady Sod shock-tube problems and the two-dimensional channel flows with different Knudsen numbers. The numerical remainder-effects of the difference schemes are investigated aad analyzed based on the computed results. The ways of improving the computational efficiency of the gaskinetic numerical method and the computing principles of difference discretization are discussed.展开更多
Geometric inversion is applied to two-dimensional Stokes flow in view to find new Stokes flow solutions. The principle of this method and the relations between the reference and inverse fluid velocity fields are prese...Geometric inversion is applied to two-dimensional Stokes flow in view to find new Stokes flow solutions. The principle of this method and the relations between the reference and inverse fluid velocity fields are presented. They are followed by applications to the flow between two parallel plates induced by a rotating or a translating cylinder. Thus hydrodynamic characteristics of flow around circular bodies obtained by inversion of the plates are thus deduced. Typically fluid flow patterns around two circular cylinders in contact placed in the centre of a rotating or a translating circular cylinder are illustrated.展开更多
Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contai...Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contains both structural and aerodynamic nonlinearities. Hopf bifurcation theory is used to analyze the flutter speed of the system. The effects of system parameters on the flutter speed are studied. The 4th order Runge-Kutta method is used to calculate the stable limit cycle responses and chaotic motions of the aeroelastic system. Results show that the number and the stability of equilibrium points of the system vary with the increase of flow speed. Besides the simple limit cycle response of period 1, there are also period-doubling responses and chaotic motions in the flutter system. The route leading to chaos in the aeroelastic model used here is the period-doubling bifurcation. The chaotic motions in the system occur only when the flow speed is higher than the linear divergent speed and the initial condition is very small. Moreover, the flow speed regions in which the system behaves chaos axe very narrow.展开更多
To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitr...To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation（BIE） representation solves the two-dimensional convected Helmholtz equation（CHE） and its fundamental solution, which must satisfy a new Sommerfeld radiation condition（SRC） in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green＇s kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.展开更多
An experimental investigation was performed to investigate two-dimensional axial velocity field at downstream of the 90°double bend pipe with and without inlet swirling condition. The main objectives are to fi...An experimental investigation was performed to investigate two-dimensional axial velocity field at downstream of the 90°double bend pipe with and without inlet swirling condition. The main objectives are to find separation region and observe the influence of inlet swirling flow on the velocity fluctuation using ultrasound technique. The experiments were carried out in the pipe at Reynolds number Re = 1 × 104. In case of inlet swirling flow condition, a rotary swirler was used as swirling generator, and the swirl number was setup S = 1. The ultrasonic measurements were taken at four downstream locations of the second bend pipe. Phased Array Ultrasonic Velocity Profiler (Phased Array UVP) technique was applied to obtain the two-dimensional velocity of the fluid and the axial and tangential velocity fluctuation. It was found that the secondary reverse flow became smaller at the downstream from the bend when the inlet condition on the first bend was swirling flow. In addition, inlet swirling condition influenced mainly on the tangential velocity fluctuation, and its maximum turbulence intensity was 40%.展开更多
By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improv...By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number （Re） flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.展开更多
The proliferation of wearable biodevices has boosted the development of soft,innovative,and multifunctional materials for human health monitoring.The integration of wearable sensors with intelligent systems is an over...The proliferation of wearable biodevices has boosted the development of soft,innovative,and multifunctional materials for human health monitoring.The integration of wearable sensors with intelligent systems is an overwhelming tendency,providing powerful tools for remote health monitoring and personal health management.Among many candidates,two-dimensional(2D)materials stand out due to several exotic mechanical,electrical,optical,and chemical properties that can be efficiently integrated into atomic-thin films.While previous reviews on 2D materials for biodevices primarily focus on conventional configurations and materials like graphene,the rapid development of new 2D materials with exotic properties has opened up novel applications,particularly in smart interaction and integrated functionalities.This review aims to consolidate recent progress,highlight the unique advantages of 2D materials,and guide future research by discussing existing challenges and opportunities in applying 2D materials for smart wearable biodevices.We begin with an in-depth analysis of the advantages,sensing mechanisms,and potential applications of 2D materials in wearable biodevice fabrication.Following this,we systematically discuss state-of-the-art biodevices based on 2D materials for monitoring various physiological signals within the human body.Special attention is given to showcasing the integration of multi-functionality in 2D smart devices,mainly including self-power supply,integrated diagnosis/treatment,and human–machine interaction.Finally,the review concludes with a concise summary of existing challenges and prospective solutions concerning the utilization of2D materials for advanced biodevices.展开更多
This paper aims to numerically explore the characteristics of unsteady cavitating flow around a NACA0015 hydrofoil,with a focus on vorticity attributes.The simulation utilizes a homogeneous mixture model coupled with ...This paper aims to numerically explore the characteristics of unsteady cavitating flow around a NACA0015 hydrofoil,with a focus on vorticity attributes.The simulation utilizes a homogeneous mixture model coupled with a filter-based density correction turbulence model and a modified Zwart cavitation model.The study investigates the dynamic cavitation features of the thermal fluid around the hydrofoil at various incoming flow velocities.It systematically elucidates the evolution of cavitation and vortex dynamics corresponding to each velocity condition.The results indicate that with increasing incoming flow velocity,distinct cavitation processes take place in the flow field.展开更多
Large size titanium alloy parts are widely used in aerospace.However,they are difficult to manufacture using mechanical cutting technology because of severe tool wear.Electrochemical jet machining is a promising techn...Large size titanium alloy parts are widely used in aerospace.However,they are difficult to manufacture using mechanical cutting technology because of severe tool wear.Electrochemical jet machining is a promising technology to achieve high efficiency,because it has high machining flexibility and no machining tool wear.However,reports on the macro electrochemical jet machining of large size titanium alloy parts are very scarce,because it is difficult to achieve effective constraint of the flow field in macro electrochemical jet machining.In addition,titanium alloy is very sensitive to fluctuation of the flow field,and a turbulent flow field would lead to serious stray corrosion.This paper reports a series of investigations of the electrochemical jet machining of titanium alloy parts.Based on the flow analysis and experiments,the machining flow field was effectively constrained.TB6 titanium alloy part with a perimeter of one meter was machined.The machined surface was smooth with no obvious machining defects.The machining process was particularly stable with no obvious spark discharge.The research provides a reference for the application of electrochemical jet machining technology to achieve large allowance material removal in the machining of large titanium alloy parts.展开更多
基金supported by the NSFC(12474071)Natural Science Foundation of Shandong Province(ZR2024YQ051,ZR2025QB50)+6 种基金Guangdong Basic and Applied Basic Research Foundation(2025A1515011191)the Shanghai Sailing Program(23YF1402200,23YF1402400)funded by Basic Research Program of Jiangsu(BK20240424)Open Research Fund of State Key Laboratory of Crystal Materials(KF2406)Taishan Scholar Foundation of Shandong Province(tsqn202408006,tsqn202507058)Young Talent of Lifting engineering for Science and Technology in Shandong,China(SDAST2024QTB002)the Qilu Young Scholar Program of Shandong University。
文摘As emerging two-dimensional(2D)materials,carbides and nitrides(MXenes)could be solid solutions or organized structures made up of multi-atomic layers.With remarkable and adjustable electrical,optical,mechanical,and electrochemical characteristics,MXenes have shown great potential in brain-inspired neuromorphic computing electronics,including neuromorphic gas sensors,pressure sensors and photodetectors.This paper provides a forward-looking review of the research progress regarding MXenes in the neuromorphic sensing domain and discussed the critical challenges that need to be resolved.Key bottlenecks such as insufficient long-term stability under environmental exposure,high costs,scalability limitations in large-scale production,and mechanical mismatch in wearable integration hinder their practical deployment.Furthermore,unresolved issues like interfacial compatibility in heterostructures and energy inefficiency in neu-romorphic signal conversion demand urgent attention.The review offers insights into future research directions enhance the fundamental understanding of MXene properties and promote further integration into neuromorphic computing applications through the convergence with various emerging technologies.
基金supported by the National Natural Science Foundation of China(Nos.11372340 and 11732016)
文摘The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller （HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257-311 （2004））. By analyzing the distribution of the finite-time Lyapunov exponent （FTLE） along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE. This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.
基金Supported by National Natural Science Foundation of China(Grant No.52205072).
文摘Owing to their rolling friction,two-dimensional piston pumps are highly suitable as power components for electro-hydrostatic actuators(EHAs).These pumps are particularly advantageous for applications requiring high efficiency and reliability.However,the ambiguity surrounding the output flow characteristics of individual two-dimensional pumps poses a significant challenge in achieving precise closed-loop control of the EHA positions.To address this issue,this study established a comprehensive numerical model that included gap leakage to analyze the impact of leakage on the output flow characteristics of a two-dimensional piston pump.The validity of the numerical analysis was indirectly confirmed through meticulous measurements of the leakage and volumetric efficiency,ensuring robust results.The research findings indicated that,at lower pump speeds,leakage significantly affected the output flow rate,leading to potential inefficiencies in the system.Conversely,at higher rotational speeds,the impact of leakage was less pronounced,implying that the influence of leakage on the pump outlet flow must be carefully considered and managed for EHAs to perform position servo control.Additionally,the research demonstrates that two-dimensional motion does not have a unique or additional effect on pump leakage,thus simplifying the design considerations.Finally,the study concluded that maintaining an oil-filled leakage environment is beneficial because it helps reduce the impact of leakage and enhances the overall volumetric efficiency of the pump system.
基金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.
文摘Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional （2D） multi-fluid flows, as well as construct the continuous Roe matrix for the whole flow field. The interface capture equations and fluid dynamic conservative equations are coupled together and solved by using any high-resolution schemes that usually suit for the single-fluid flows. Some numerical examples are given to illustrate the solution of 1D and 2D multi-fluid Riemann problems.
基金Project supported by National Natural Science Fundation of China
文摘If the vertically-mixing time is comparable with the period of oscillatory current, the contaminant contraction may occur, then the coefficient of shear dispersion will be negative (singularity). According to the two-dimensional delay-diffusion equation derived by Liu in another paper, The equations for X(t), Y(t), central displacement as well as dispersion tensor were derived. Generally speaking, δDij/δτ is positive when τ is small; but when τ is large, the memory functions may be negative. The analytical expressions for Dij and X,Y were also given.
基金supported by the National Basic Research Program of China (973 Program)(Grant No.2011CB409902)the Knowledge Innovation Project of the Chinese Academy of Sciences (No.KZCX2-YW-302)
文摘Debris flows are recurrent natural hazards in many mountainous regions.This paper presents a numerical study on the propagation of debris flows in natural erodible open channels,in which the bed erosion and sedimentation processes are important.Based on the Bingham fluid theory,a mathematical model of the two-dimensional non-constant debris flow is developed.The governing equations include the continuity and momentum conservation equations of debris flow,the sediment convection-diffusion equation,the bed erosion-deposition equation and the bed-sediment size gradation adjustment equation.The yield stress and shear stress components are included to describe the dynamic rheological properties.The upwind control-volume Finite Volume Method (FVM) is applied to discretize the convection terms.The improved SIMPLE algorithm with velocity-free-surface coupled correction is developed to solve the equations on non-orthogonal,quadrilateral grids.The model is applied to simulate a debris flow event in Jiangjia Gully,Yunnan Province and to predict the flow pattern and bed erosion-deposition processes.The results show the effectiveness of the proposed numercial model in debris flow simulation and potential hazard analysis.
基金partially supported by the National Natural Science Foundation of China(No.11972189)the Natural Science Foundation of Jiangsu Province(No.BK20190391)+1 种基金the Natural Science Foundation of Anhui Province(No.1908085QF260)the Priority Academic Program Development of Jiangsu Higher Education Institutions。
文摘A graphics processing unit(GPU)-accelerated discontinuous Galerkin(DG)method is presented for solving two-dimensional laminar flows.The DG method is ported from central processing unit to GPU in a way of achieving GPU speedup through programming under the compute unified device architecture(CUDA)model.The CUDA kernel subroutines are designed to meet with the requirement of high order computing of DG method.The corresponding data structures are constructed in component-wised manners and the thread hierarchy is manipulated in cell-wised or edge-wised manners associated with related integrals involved in solving laminar Navier-Stokes equations,in which the inviscid and viscous flux terms are computed by the local lax-Friedrichs scheme and the second scheme of Bassi&Rebay,respectively.A strong stability preserving Runge-Kutta scheme is then used for time marching of numerical solutions.The resulting GPU-accelerated DG method is first validated by the traditional Couette flow problems with different mesh sizes associated with different orders of approximation,which shows that the orders of convergence,as expected,can be achieved.The numerical simulations of the typical flows over a circular cylinder or a NACA 0012 airfoil are then carried out,and the results are further compared with the analytical solutions or available experimental and numerical values reported in the literature,as well as with a performance analysis of the developed code in terms of GPU speedups.This shows that the costs of computing time of the presented test cases are significantly reduced without losing accuracy,while impressive speedups up to 69.7 times are achieved by the present method in comparison to its CPU counterpart.
文摘The regions with shear stress and mean velocity gradient of opposite sign often exist in complex turbulent shear flows.In these cases,the eddy viscosity hypothesis breaks down.Hinze regards the,departure from eddy viscosity hypothesis as a result from transportation of mean momentum over distance by the large structures and arrives at a shear stress expression including the second order derivatives of the mean velocity.However,his expression greatly overestimates the shear stress.This implies that the flow particles are unlikely to have enough memory of the mean momentum over distance.By assuming the departure from eddy viscosity hypothesis as a result from transportation of the shear stress contained in smaller eddies over distance by the large structures,the present author has arrived at a new shear stress expression.The shear stress estimated so far is in good agreement with the experiments.
文摘Stably stratified flows over a two-dimensional hill are investigated in a channel of finite depth using a three-dimensional direct numerical simulation (DNS). The present study follows onto our previous two-dimensional DNS studies of stably stratified flows over a hill in a channel of finite depth and provides a more realistic simulation of atmospheric flows than our previous studies. A hill with a constant cross-section in the spanwise (y) direction is placed in a 3-D computational domain. As in the previous 2-D simulations, to avoid the effect of the ground boundary layer that develops upstream of the hill, no-slip conditions are imposed only on the hill surface and the surface downstream of the hill;slip conditions are imposed on the surface upstream of the hill. The simulated 3-D flows are discussed by comparing them to the simulated 2-D flows with a focus on the effect of the stable stratification on the non-periodic separation and reattachment of the flow behind the hill. In neutral (K = 0, where K is a non-dimensional stability parameter) and weakly stable (K = 0.8) conditions, 3-D flows over a hill differ clearly from 2-D flows over a hill mainly because of the three-dimensionality of the flow, that is the development of a spanwise flow component in the 3-D flows. In highly stable conditions (K = 1, 1.3), long-wavelength lee waves develop downstream of the hill in both 2-D and 3-D flows, and the behaviors of the 2-D and 3-D flows are similar in the vicinity of the hill. In other words, the spanwise component of the 3-D flows is strongly suppressed in highly stable conditions, and the flow in the vicinity of the hill becomes approximately two-dimensional in the x and z directions.
基金Project supported by the National Natural Science Foundation of China (Grant No 10765002)Guangxi Natural Science Foundation (Grant No 0542045)
文摘Boundary conditions （BCs） play an essential role in lattice Boltzmann （LB） simulations. This paper investigates several most commonly applied BCs by evaluating the relative L2-norm errors of the LB simulations for two-dimensional （2-D） Poiseuille flow. It is found that the relative L2-norm error resulting from FHML＇s BC is smaller than that from other BCs as a whole. Then, based on the FHML＇s BC, it formulates an LB model for simulating fluid flows in 2-D channel with complex geometries. Afterwards, the flows between two inclined plates, in a pulmonary blood vessel and in a blood vessel with local expansion region, are simulated. The numerical results are in good agreement with the analytical predictions and clearly show that the model is effective. It is expected that the model can be extended to simulate some real biologic flows, such as blood flows in arteries, vessels with stenosises, aneurysms and bifurcations,
基金supported by the National Natural Science Foundation of China (No.10621062)the Research Fund for Next Generation of General Armament Department (No.9140A13050207KG29)
文摘Several kinds of explicit and implicit finite-difference schemes directly solving the discretized velocity distribution functions are designed with precision of different orders by analyzing the inner characteristics of the gas-kinetic numerical algorithm for Boltzmann model equation. The peculiar flow phenomena and mechanism from various flow regimes are revealed in the numerical simulations of the unsteady Sod shock-tube problems and the two-dimensional channel flows with different Knudsen numbers. The numerical remainder-effects of the difference schemes are investigated aad analyzed based on the computed results. The ways of improving the computational efficiency of the gaskinetic numerical method and the computing principles of difference discretization are discussed.
文摘Geometric inversion is applied to two-dimensional Stokes flow in view to find new Stokes flow solutions. The principle of this method and the relations between the reference and inverse fluid velocity fields are presented. They are followed by applications to the flow between two parallel plates induced by a rotating or a translating cylinder. Thus hydrodynamic characteristics of flow around circular bodies obtained by inversion of the plates are thus deduced. Typically fluid flow patterns around two circular cylinders in contact placed in the centre of a rotating or a translating circular cylinder are illustrated.
基金supported by the National Natural Science Foundation of China and China Academy of Engineering Physics(No. 10576024).
文摘Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contains both structural and aerodynamic nonlinearities. Hopf bifurcation theory is used to analyze the flutter speed of the system. The effects of system parameters on the flutter speed are studied. The 4th order Runge-Kutta method is used to calculate the stable limit cycle responses and chaotic motions of the aeroelastic system. Results show that the number and the stability of equilibrium points of the system vary with the increase of flow speed. Besides the simple limit cycle response of period 1, there are also period-doubling responses and chaotic motions in the flutter system. The route leading to chaos in the aeroelastic model used here is the period-doubling bifurcation. The chaotic motions in the system occur only when the flow speed is higher than the linear divergent speed and the initial condition is very small. Moreover, the flow speed regions in which the system behaves chaos axe very narrow.
基金supported by National Engineering School of Tunis (No.13039.1)
文摘To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation（BIE） representation solves the two-dimensional convected Helmholtz equation（CHE） and its fundamental solution, which must satisfy a new Sommerfeld radiation condition（SRC） in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green＇s kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.
文摘An experimental investigation was performed to investigate two-dimensional axial velocity field at downstream of the 90°double bend pipe with and without inlet swirling condition. The main objectives are to find separation region and observe the influence of inlet swirling flow on the velocity fluctuation using ultrasound technique. The experiments were carried out in the pipe at Reynolds number Re = 1 × 104. In case of inlet swirling flow condition, a rotary swirler was used as swirling generator, and the swirl number was setup S = 1. The ultrasonic measurements were taken at four downstream locations of the second bend pipe. Phased Array Ultrasonic Velocity Profiler (Phased Array UVP) technique was applied to obtain the two-dimensional velocity of the fluid and the axial and tangential velocity fluctuation. It was found that the secondary reverse flow became smaller at the downstream from the bend when the inlet condition on the first bend was swirling flow. In addition, inlet swirling condition influenced mainly on the tangential velocity fluctuation, and its maximum turbulence intensity was 40%.
基金Project supported by the National Natural Science Foundation of China （Grant No 70271069）.
文摘By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number （Re） flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.
基金the support from the National Natural Science Foundation of China(22272004,62272041)the Fundamental Research Funds for the Central Universities(YWF-22-L-1256)+1 种基金the National Key R&D Program of China(2023YFC3402600)the Beijing Institute of Technology Research Fund Program for Young Scholars(No.1870011182126)。
文摘The proliferation of wearable biodevices has boosted the development of soft,innovative,and multifunctional materials for human health monitoring.The integration of wearable sensors with intelligent systems is an overwhelming tendency,providing powerful tools for remote health monitoring and personal health management.Among many candidates,two-dimensional(2D)materials stand out due to several exotic mechanical,electrical,optical,and chemical properties that can be efficiently integrated into atomic-thin films.While previous reviews on 2D materials for biodevices primarily focus on conventional configurations and materials like graphene,the rapid development of new 2D materials with exotic properties has opened up novel applications,particularly in smart interaction and integrated functionalities.This review aims to consolidate recent progress,highlight the unique advantages of 2D materials,and guide future research by discussing existing challenges and opportunities in applying 2D materials for smart wearable biodevices.We begin with an in-depth analysis of the advantages,sensing mechanisms,and potential applications of 2D materials in wearable biodevice fabrication.Following this,we systematically discuss state-of-the-art biodevices based on 2D materials for monitoring various physiological signals within the human body.Special attention is given to showcasing the integration of multi-functionality in 2D smart devices,mainly including self-power supply,integrated diagnosis/treatment,and human–machine interaction.Finally,the review concludes with a concise summary of existing challenges and prospective solutions concerning the utilization of2D materials for advanced biodevices.
文摘This paper aims to numerically explore the characteristics of unsteady cavitating flow around a NACA0015 hydrofoil,with a focus on vorticity attributes.The simulation utilizes a homogeneous mixture model coupled with a filter-based density correction turbulence model and a modified Zwart cavitation model.The study investigates the dynamic cavitation features of the thermal fluid around the hydrofoil at various incoming flow velocities.It systematically elucidates the evolution of cavitation and vortex dynamics corresponding to each velocity condition.The results indicate that with increasing incoming flow velocity,distinct cavitation processes take place in the flow field.
基金the National Natural Science Foundation of China(No.52205468)China Postdoctoral Science Foundation(No.2022M710061 and No.2023T160277)Natural Science Foundation of Jiangsu Province(No.BK20210755)。
文摘Large size titanium alloy parts are widely used in aerospace.However,they are difficult to manufacture using mechanical cutting technology because of severe tool wear.Electrochemical jet machining is a promising technology to achieve high efficiency,because it has high machining flexibility and no machining tool wear.However,reports on the macro electrochemical jet machining of large size titanium alloy parts are very scarce,because it is difficult to achieve effective constraint of the flow field in macro electrochemical jet machining.In addition,titanium alloy is very sensitive to fluctuation of the flow field,and a turbulent flow field would lead to serious stray corrosion.This paper reports a series of investigations of the electrochemical jet machining of titanium alloy parts.Based on the flow analysis and experiments,the machining flow field was effectively constrained.TB6 titanium alloy part with a perimeter of one meter was machined.The machined surface was smooth with no obvious machining defects.The machining process was particularly stable with no obvious spark discharge.The research provides a reference for the application of electrochemical jet machining technology to achieve large allowance material removal in the machining of large titanium alloy parts.
