Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinni...Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.展开更多
To ensure the control of the precision of air-fuel ratio(AFR)of port fuel injection(PFI)spark ignition(SI)engines,a chaos radial basis function(RBF)neural network is used to predict the air intake flow of the engine.T...To ensure the control of the precision of air-fuel ratio(AFR)of port fuel injection(PFI)spark ignition(SI)engines,a chaos radial basis function(RBF)neural network is used to predict the air intake flow of the engine.The data of air intake flow is proved to be multidimensionally nonlinear and chaotic.The RBF neural network is used to train the reconstructed phase space of the data.The chaos algorithm is employed to optimize the weights of output layer connection and the radial basis center of Gaussian function in hidden layer.The simulation results obtained from Matlab/Simulink illustrate that the model has higher accuracy compared to the conventional RBF model.The mean absolute error and the mean relative error of the chaos RBF model can reach 0.0017 and 0.48,respectively.展开更多
This paper presents a KIVA-3 code based numerical model for three-dimensional transient intake flow in the intake port-valve-cylinder system of internal combustion engine using body-fitted technique, which can be used...This paper presents a KIVA-3 code based numerical model for three-dimensional transient intake flow in the intake port-valve-cylinder system of internal combustion engine using body-fitted technique, which can be used in numerical study on internal combustion engine with vertical and inclined valves, and has higher calculation precision. A numerical simulation (on the intake process of a two-valve engine with a semi-sphere combustion chamber and a radial intake port) is provided for analysis of the velocity field and pressure field of different plane at different crank angles. The results revealed the formation of the tumble motion, the evolution of flow field parameters and the variation of tumble ratios as important information for the design of engine in-take system.展开更多
In this paper, the effect of non-uniform heat flux on heat transfer in boundary layer stagnation-point flow over a shrinking sheet is studied. The variable boundary heat fluxes are considered of two types: direct pow...In this paper, the effect of non-uniform heat flux on heat transfer in boundary layer stagnation-point flow over a shrinking sheet is studied. The variable boundary heat fluxes are considered of two types: direct power-law variation with the distance along the sheet and inverse power-law variation with the distance. The governing partial differential equations (PDEs) are transformed into non linear self-similar ordinary differential equations (ODEs) by similarity transformations, and then those are solved using very efficient shooting method. The direct variation and inverse variation of heat flux along the sheet have completely different effects on the temperature distribution. Moreover, the heat transfer characteristics in the presence of non-uniform heat flux for several values of physical parameters are also found to be interesting.展开更多
Currently, the flow field of annular seals disturbed by the circular whirl motion of rotors is usually solved using computational fluid dynamics(CFD) to evaluate the five rotordynamic coefficients. The simulations are...Currently, the flow field of annular seals disturbed by the circular whirl motion of rotors is usually solved using computational fluid dynamics(CFD) to evaluate the five rotordynamic coefficients. The simulations are based on the traditional quasi-steady method. In this work, an improved quasi-steady method along with the transient method was presented to compute the rotordynamic coefficients of a long seal. By comparisons with experimental data, the shortcomings of quasi-steady methods have been identified. Then, the effects of non-uniform incoming flow on seal dynamic coefficients were studied by transient simulations. Results indicate that the long seal has large cross stiffness k and direct mass M which are not good for rotor stability, while the transient method is more suitable for the long seal for its excellent performance in predicting M. When the incoming flow is non-uniform, the stiffness coefficients vary with the eccentric directions. Based on the rotordynamic coefficients under uniform incoming flow, the linearized fluid force formulas, which can consider the effects of non-uniform incoming flow, have been presented and can well explain the varying-stiffness phenomenon.展开更多
In hydraulics,when we deal with the problem of sand particles moving relative to the surrounding water,Stokes'formula of resistance has usually been used to render the velocity of sedimentation of the particles.Bu...In hydraulics,when we deal with the problem of sand particles moving relative to the surrounding water,Stokes'formula of resistance has usually been used to render the velocity of sedimentation of the particles.But such an approach has not been proved rigorously,and its accuracy must be carefully considered.In this paper,we discuss the problem of a sphere moving in a non-uniform flow field,on the basis of the fundamental theory of hydrodynamics.We introduce two assumptions:i)the diameter of the sphere is much smaller than the linear dimension of the flow field,and ii)the velocity of the sphere relative to the surrounding water is very small.Using these two assumptions,we solve the linearized Navier-Stokes equations and equations of continuity by the method of Laplace transform,and finally we obtain a formula for the resistance acting on a sphere moving in a non-uniform flow field.展开更多
Based on finite volume method, the pressure drop and heat transfer characteristics of one smooth tube and ten different axisymmetric corrugated tubes, including two with uniform corrugation and eight with non-uniform ...Based on finite volume method, the pressure drop and heat transfer characteristics of one smooth tube and ten different axisymmetric corrugated tubes, including two with uniform corrugation and eight with non-uniform corrugation, have been studied. A physical model of the corrugated tube was built, then the numerical simulation of the model was carried out and the numerical simulation results were compared with the empirical formula.The results show that: the friction factor decreases with the increase of Reynolds number ranging from 6000 to 57000, the value of which in the corrugated tubes with non-uniform corrugation(tube 03–10) are smaller than those with uniform corrugation(tube 01–02). The geometry parameters of tube(01) have advantages on the heat transfer enhancement in low Reynolds number flow region(from 6000 to 13000) and tube(07–08)have advantages on the heat transfer enhancement in high Reynolds number flow region(from 13000 to 57000). The vortex, existed in each area between two adjacent corrugations called second flow region, is the root of the enhancement on heat transfer in the corrugated tubes. The effectiveness factor decreases with the increasing of Reynolds number and the performances of the corrugated tubes with pitch of 12.5 mm have advantages than these of 10 mm under the same corrugation geometric parameter.展开更多
An analysis has been carried out to study the effect of hydrodynamic laminar boundary layer flow and heat transfer of a dusty fluid over an unsteady stretching surface in the presence of non-uniform heat source/sink. ...An analysis has been carried out to study the effect of hydrodynamic laminar boundary layer flow and heat transfer of a dusty fluid over an unsteady stretching surface in the presence of non-uniform heat source/sink. Heat transfer characteristics are examined for two different kinds of boundary conditions, namely 1) variable wall temperature and 2) variable heat flux. The governing partial differential equations are transformed to system of ordinary differential equations. These equations are solved numerically by applying RKF-45 method. The effects of various physical parameters such as magnetic parameter, dust interaction parameter, number density, Prandtl number, Eckert number, heat source/sink parameter and unsteadiness parameter on velocity and temperature profiles are studied.展开更多
Exact simulation of the acoustic performance is essential to the engineering application for a vehicle intake system. The rectangular-pulse method based on the computational fluid dynamics approach was employed for ca...Exact simulation of the acoustic performance is essential to the engineering application for a vehicle intake system. The rectangular-pulse method based on the computational fluid dynamics approach was employed for calculating the transmission loss. Firstly, the transmission loss of the single-cavity element was simulated without any airflow, and the effects of different structural parameters on the acoustic performance were investigated comprehensively. Secondly, the static transmission loss of the perforated intake pipe was obtained by the rectangular-pulse method, which is proved to be accurate enough compared with the result by finite element method. Thirdly, under the different conditions of the mean airflow and the operating temperature, the specific transmission loss was acquired respectively. In general, the peaks of the transmission loss are shifted to the lower frequency range because of the reverse airflow, but the amplitudes are irregularly changed. Besides, when the operating temperature increases, the peaks are shifted to the higher frequencies. Finally, with the designed perforated pipe installed to the intake system, the road tests were proceeded to evaluate the actual acoustic performance, and the result indicates that the intake sound pressure level is greatly attenuated. Typically in the range of 600–1500 Hz, the insertion loss of the intake noise at the decelerating moment is almost 20 d B(A), and the overall noise is reduced more than 14.2 d B(A). In conclusion, the perforated intake pipe has been proved excellent in improving the acoustic performance of intake system and could provide the guidance for the automotive engineering application.展开更多
Most of water flow in open channel or in river belongs to steady non-uniform flow. The surface profiles are caused by changes of channel section. It is very important to analyze its computation. According to the regu...Most of water flow in open channel or in river belongs to steady non-uniform flow. The surface profiles are caused by changes of channel section. It is very important to analyze its computation. According to the regularity of its surface change, the suitable sectional dimensions of open channel or flood control work can be designed. Commonly, computation of non-uniform flow adopts the traditional methods by hand or by graphic method. The speed and precision of computation are restricted. In this paper, a software to calculate water surface profile is introduced. The software is put forward by using C++ .By means of interpolate method and dialogue between user and computer, we can calculate the water surface profile much more quickly and exactly.展开更多
In the metallurgical industries, it is very important to characterize the flow of molten metals in open channels given that they are transported through these devices to different plant sections. Howeve...In the metallurgical industries, it is very important to characterize the flow of molten metals in open channels given that they are transported through these devices to different plant sections. However, unlike the flow of water which has been studied since ancient times, the flow of molten metals in open channels has received little attention. The unsteady non-uniform flow of blast furnace molten pig iron in a rectangular open channel is analyzed in this work by numerical solution of the Saint-Venant equations. The influence of mesh size on the convergence of molten metal height is studied to determine the proper mesh and time step sizes. A sinusoidal inflow pulse is imposed at the entrance of the channel in order to analyze the propagation of the resulting wave. The influence of the angle of inclination of the channel and the roughness coefficient of the walls on the amplitude and the dynamic behavior of the height of the molten metal are analyzed. Phase portraits of the channel state variables are constructed and interpreted. Numerical simulations show that as the angle of inclination of the channel increases, the amplitude of the formed wave decreases. From 10 degrees onwards, the peak of the wave descends even below the initial height. On the other hand, the roughness coefficient affects the molten pig iron height profiles in an inverse way than the angle of inclination. The amplitude of the formed wave increases as the roughness coefficient increases.展开更多
基金This work is supported by the National Natural Science Foundation of China(No.52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(No.BYESS2023262)Science Foundation of China University of Petroleum,Beijing(No.2462022BJRC004).
文摘Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.
基金Project(51176014)supported by the National Natural Science Foundation of ChinaProject(2016JJ2003)supported by Natural Scienceof Hunan Province,ChinaProject(KF1605)supported by Key Laboratory of Safety Design and Reliability Technology of Engineering Vehicle in Hunan Province,China。
文摘To ensure the control of the precision of air-fuel ratio(AFR)of port fuel injection(PFI)spark ignition(SI)engines,a chaos radial basis function(RBF)neural network is used to predict the air intake flow of the engine.The data of air intake flow is proved to be multidimensionally nonlinear and chaotic.The RBF neural network is used to train the reconstructed phase space of the data.The chaos algorithm is employed to optimize the weights of output layer connection and the radial basis center of Gaussian function in hidden layer.The simulation results obtained from Matlab/Simulink illustrate that the model has higher accuracy compared to the conventional RBF model.The mean absolute error and the mean relative error of the chaos RBF model can reach 0.0017 and 0.48,respectively.
文摘This paper presents a KIVA-3 code based numerical model for three-dimensional transient intake flow in the intake port-valve-cylinder system of internal combustion engine using body-fitted technique, which can be used in numerical study on internal combustion engine with vertical and inclined valves, and has higher calculation precision. A numerical simulation (on the intake process of a two-valve engine with a semi-sphere combustion chamber and a radial intake port) is provided for analysis of the velocity field and pressure field of different plane at different crank angles. The results revealed the formation of the tumble motion, the evolution of flow field parameters and the variation of tumble ratios as important information for the design of engine in-take system.
基金the National Board for Higher Mathematics(NBHM),DAE,Mumbai,India
文摘In this paper, the effect of non-uniform heat flux on heat transfer in boundary layer stagnation-point flow over a shrinking sheet is studied. The variable boundary heat fluxes are considered of two types: direct power-law variation with the distance along the sheet and inverse power-law variation with the distance. The governing partial differential equations (PDEs) are transformed into non linear self-similar ordinary differential equations (ODEs) by similarity transformations, and then those are solved using very efficient shooting method. The direct variation and inverse variation of heat flux along the sheet have completely different effects on the temperature distribution. Moreover, the heat transfer characteristics in the presence of non-uniform heat flux for several values of physical parameters are also found to be interesting.
基金Project(51276213)supported by the National Natural Science Foundation of ChinaProject(2013BAF01B00)supported by the National Science and Technology Support Program of China
文摘Currently, the flow field of annular seals disturbed by the circular whirl motion of rotors is usually solved using computational fluid dynamics(CFD) to evaluate the five rotordynamic coefficients. The simulations are based on the traditional quasi-steady method. In this work, an improved quasi-steady method along with the transient method was presented to compute the rotordynamic coefficients of a long seal. By comparisons with experimental data, the shortcomings of quasi-steady methods have been identified. Then, the effects of non-uniform incoming flow on seal dynamic coefficients were studied by transient simulations. Results indicate that the long seal has large cross stiffness k and direct mass M which are not good for rotor stability, while the transient method is more suitable for the long seal for its excellent performance in predicting M. When the incoming flow is non-uniform, the stiffness coefficients vary with the eccentric directions. Based on the rotordynamic coefficients under uniform incoming flow, the linearized fluid force formulas, which can consider the effects of non-uniform incoming flow, have been presented and can well explain the varying-stiffness phenomenon.
文摘In hydraulics,when we deal with the problem of sand particles moving relative to the surrounding water,Stokes'formula of resistance has usually been used to render the velocity of sedimentation of the particles.But such an approach has not been proved rigorously,and its accuracy must be carefully considered.In this paper,we discuss the problem of a sphere moving in a non-uniform flow field,on the basis of the fundamental theory of hydrodynamics.We introduce two assumptions:i)the diameter of the sphere is much smaller than the linear dimension of the flow field,and ii)the velocity of the sphere relative to the surrounding water is very small.Using these two assumptions,we solve the linearized Navier-Stokes equations and equations of continuity by the method of Laplace transform,and finally we obtain a formula for the resistance acting on a sphere moving in a non-uniform flow field.
文摘Based on finite volume method, the pressure drop and heat transfer characteristics of one smooth tube and ten different axisymmetric corrugated tubes, including two with uniform corrugation and eight with non-uniform corrugation, have been studied. A physical model of the corrugated tube was built, then the numerical simulation of the model was carried out and the numerical simulation results were compared with the empirical formula.The results show that: the friction factor decreases with the increase of Reynolds number ranging from 6000 to 57000, the value of which in the corrugated tubes with non-uniform corrugation(tube 03–10) are smaller than those with uniform corrugation(tube 01–02). The geometry parameters of tube(01) have advantages on the heat transfer enhancement in low Reynolds number flow region(from 6000 to 13000) and tube(07–08)have advantages on the heat transfer enhancement in high Reynolds number flow region(from 13000 to 57000). The vortex, existed in each area between two adjacent corrugations called second flow region, is the root of the enhancement on heat transfer in the corrugated tubes. The effectiveness factor decreases with the increasing of Reynolds number and the performances of the corrugated tubes with pitch of 12.5 mm have advantages than these of 10 mm under the same corrugation geometric parameter.
文摘An analysis has been carried out to study the effect of hydrodynamic laminar boundary layer flow and heat transfer of a dusty fluid over an unsteady stretching surface in the presence of non-uniform heat source/sink. Heat transfer characteristics are examined for two different kinds of boundary conditions, namely 1) variable wall temperature and 2) variable heat flux. The governing partial differential equations are transformed to system of ordinary differential equations. These equations are solved numerically by applying RKF-45 method. The effects of various physical parameters such as magnetic parameter, dust interaction parameter, number density, Prandtl number, Eckert number, heat source/sink parameter and unsteadiness parameter on velocity and temperature profiles are studied.
基金Project(51705454)supported by the National Natural Science Foundation of China
文摘Exact simulation of the acoustic performance is essential to the engineering application for a vehicle intake system. The rectangular-pulse method based on the computational fluid dynamics approach was employed for calculating the transmission loss. Firstly, the transmission loss of the single-cavity element was simulated without any airflow, and the effects of different structural parameters on the acoustic performance were investigated comprehensively. Secondly, the static transmission loss of the perforated intake pipe was obtained by the rectangular-pulse method, which is proved to be accurate enough compared with the result by finite element method. Thirdly, under the different conditions of the mean airflow and the operating temperature, the specific transmission loss was acquired respectively. In general, the peaks of the transmission loss are shifted to the lower frequency range because of the reverse airflow, but the amplitudes are irregularly changed. Besides, when the operating temperature increases, the peaks are shifted to the higher frequencies. Finally, with the designed perforated pipe installed to the intake system, the road tests were proceeded to evaluate the actual acoustic performance, and the result indicates that the intake sound pressure level is greatly attenuated. Typically in the range of 600–1500 Hz, the insertion loss of the intake noise at the decelerating moment is almost 20 d B(A), and the overall noise is reduced more than 14.2 d B(A). In conclusion, the perforated intake pipe has been proved excellent in improving the acoustic performance of intake system and could provide the guidance for the automotive engineering application.
文摘Most of water flow in open channel or in river belongs to steady non-uniform flow. The surface profiles are caused by changes of channel section. It is very important to analyze its computation. According to the regularity of its surface change, the suitable sectional dimensions of open channel or flood control work can be designed. Commonly, computation of non-uniform flow adopts the traditional methods by hand or by graphic method. The speed and precision of computation are restricted. In this paper, a software to calculate water surface profile is introduced. The software is put forward by using C++ .By means of interpolate method and dialogue between user and computer, we can calculate the water surface profile much more quickly and exactly.
文摘In the metallurgical industries, it is very important to characterize the flow of molten metals in open channels given that they are transported through these devices to different plant sections. However, unlike the flow of water which has been studied since ancient times, the flow of molten metals in open channels has received little attention. The unsteady non-uniform flow of blast furnace molten pig iron in a rectangular open channel is analyzed in this work by numerical solution of the Saint-Venant equations. The influence of mesh size on the convergence of molten metal height is studied to determine the proper mesh and time step sizes. A sinusoidal inflow pulse is imposed at the entrance of the channel in order to analyze the propagation of the resulting wave. The influence of the angle of inclination of the channel and the roughness coefficient of the walls on the amplitude and the dynamic behavior of the height of the molten metal are analyzed. Phase portraits of the channel state variables are constructed and interpreted. Numerical simulations show that as the angle of inclination of the channel increases, the amplitude of the formed wave decreases. From 10 degrees onwards, the peak of the wave descends even below the initial height. On the other hand, the roughness coefficient affects the molten pig iron height profiles in an inverse way than the angle of inclination. The amplitude of the formed wave increases as the roughness coefficient increases.