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.展开更多
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.展开更多
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 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.展开更多
Because the real atmosphere is non-uniformly saturated, the generalized potential temperature is introduced. The convective vorticity vector, which can depict the occurrence and development of mesoscale deep convectiv...Because the real atmosphere is non-uniformly saturated, the generalized potential temperature is introduced. The convective vorticity vector, which can depict the occurrence and development of mesoscale deep convective systems, is modified and re-derived in a nonuniformly saturated moist atmosphere (C*). Then, a case study is performed for a frontal rainfall event which occurred near the middle and lower reaches of the Yangtze River in China. The diagnostic results of C* show that, in the lower troposphere, the vertical component of C* (Cz*) can diagnose developments and movements of precipitation and convection better than those of Cm (Cmz, in saturated moist flow) and C (Cz, in dry flow). Cz* is a good predictor for precipitation analyses as well.展开更多
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.展开更多
A numerical study on the acoustic radiation of a propeller interacting with non-uniform inflow has been conducted. Real geometry of a marine propeller DTMB 4118 is used in the calculation, and sliding mesh technique i...A numerical study on the acoustic radiation of a propeller interacting with non-uniform inflow has been conducted. Real geometry of a marine propeller DTMB 4118 is used in the calculation, and sliding mesh technique is adopted to deal with the rotational motion of the propeller. The performance of the DES (Detached Eddy Simulation) approach at capturing the unsteady forces and moments on the propeller is compared with experiment. Far-field sound radiation is predicted by the formation 1A developed by Farassat, an integral solution of FW-H (Ffowcs Williams-Hawkings) equation in time domain. The sound pressure and directivity patterns of the propeller operating in two specific velocity distributions are discussed.展开更多
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.展开更多
This study presents the effect of non-uniform heat source on the magneto-hydrodynamic flow of nanofluid across an expanding plate with consideration of the homogeneous-heterogeneous reactions and thermal radiation eff...This study presents the effect of non-uniform heat source on the magneto-hydrodynamic flow of nanofluid across an expanding plate with consideration of the homogeneous-heterogeneous reactions and thermal radiation effects.A nanofluid’s dynamic viscosity and effective thermal conductivity are specified with Corcione correlation.According to this correlation,the thermal conductivity is carried out by the Brownian motion.Similarity transformations reduce the governing equations concerned with energy,momentum,and concentration of nanofluid and then numerically solved.The influences of the effective parameters,e.g.,the internal heat source parameters,the volume fraction of nanofluid,the radiation parameter,the homogeneous reaction parameter,the magnetic parameter,the heterogeneous parameter and the Schmidt number are studied on the heat and flow transfer features.Further,regarding the effective parameters of the present work,the correlation for the Nusselt number has been developed.The outcomes illustrate that with the raising of the heterogeneous parameter and the homogeneous reaction parameter,the concentration profile diminishes.In addition,the outcomes point to a reverse relationship between the Nusselt number and the internal heat source parameters.展开更多
The present article aims to investigate the Graetz-Nusselt problem for blood as a non-Newtonian fluid obeying the power-law constitutive equation and flowing inside the axisymmetric tube subjected to nonuniform surfac...The present article aims to investigate the Graetz-Nusselt problem for blood as a non-Newtonian fluid obeying the power-law constitutive equation and flowing inside the axisymmetric tube subjected to nonuniform surface heat flux.After the flow field is determined by solving the continuity and the momentum equations,the energy equation is handled by employing the separation of variables method.The resulting Eigen functions and Eigen values are numerically calculated using MATLAB built-in solver BVP4C.The analysis is first conducted for the situation of constant heat flux and subsequently generalized to apply to the case of sinusoidal variation of wall heat flux along the tube length,using Duhamel’s Theorem.Furthermore,an approximate analytic solution is determined,employing an integral approach to solve the boundary layer equations.With respect to the comparison,the results of approximate solution display acceptable congruence with those of exact solution with an average error of 7.4%.Interestingly,with decreasing the power-law index,the discrepancy between the two presented methods significantly reduces.Eventually,the influences of the controlling parameters such as surface heat flux and power-law index on the non-Newtonian fluid flow’s thermal characteristics and structure are elaborately discussed.It is found that switching from constant wall heat flux to non-uniform wall heat flux that sinusoidally varies along the tube length significantly improves the simulation’s accuracy due to the better characterization of the heat transport phenomenon in non-Newtonian fluid flow through the tube.In the presence of sinusoidally varying wall heat flux with an amplitude of 200 W/m 2 and when the power-law index is 0.25,the maximum arterial wall temperature is found to be about 311.56 K.展开更多
The aim of this work is to study the effect of non-uniform single and double slot suction/injection into a steady mixed convection boundary layer flow over a vertical cone, while the axis of the cone is inline with th...The aim of this work is to study the effect of non-uniform single and double slot suction/injection into a steady mixed convection boundary layer flow over a vertical cone, while the axis of the cone is inline with the flow. The governing boundary layer equations are transformed into a non-dimensional form by a group of non-similar trans- formations. The resulting coupled non-linear partial differential equations are solved nu- merically by employing the quasi-linearization technique and an implicit finite-difference scheme. Numerical computations are performed for different values of the dimensionless parameters to display the velocity and temperature profiles graphically. Also, numerical results are presented for the skin friction and heat transfer coefficients. Results indicate that the skin friction and heat transfer coefficients increase with non-uniform slot suction, but the effect of non-uniform slot injection is just opposite.展开更多
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.
文摘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.
基金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.
基金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.
基金supported by the National Natural Science Foundation of China under grant Nos. 40805001 and 40433016
文摘Because the real atmosphere is non-uniformly saturated, the generalized potential temperature is introduced. The convective vorticity vector, which can depict the occurrence and development of mesoscale deep convective systems, is modified and re-derived in a nonuniformly saturated moist atmosphere (C*). Then, a case study is performed for a frontal rainfall event which occurred near the middle and lower reaches of the Yangtze River in China. The diagnostic results of C* show that, in the lower troposphere, the vertical component of C* (Cz*) can diagnose developments and movements of precipitation and convection better than those of Cm (Cmz, in saturated moist flow) and C (Cz, in dry flow). Cz* is a good predictor for precipitation analyses as well.
文摘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.
基金supported by the National Natural Science Foundation of China (Grant No. 11272213)
文摘A numerical study on the acoustic radiation of a propeller interacting with non-uniform inflow has been conducted. Real geometry of a marine propeller DTMB 4118 is used in the calculation, and sliding mesh technique is adopted to deal with the rotational motion of the propeller. The performance of the DES (Detached Eddy Simulation) approach at capturing the unsteady forces and moments on the propeller is compared with experiment. Far-field sound radiation is predicted by the formation 1A developed by Farassat, an integral solution of FW-H (Ffowcs Williams-Hawkings) equation in time domain. The sound pressure and directivity patterns of the propeller operating in two specific velocity distributions are discussed.
文摘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.
文摘This study presents the effect of non-uniform heat source on the magneto-hydrodynamic flow of nanofluid across an expanding plate with consideration of the homogeneous-heterogeneous reactions and thermal radiation effects.A nanofluid’s dynamic viscosity and effective thermal conductivity are specified with Corcione correlation.According to this correlation,the thermal conductivity is carried out by the Brownian motion.Similarity transformations reduce the governing equations concerned with energy,momentum,and concentration of nanofluid and then numerically solved.The influences of the effective parameters,e.g.,the internal heat source parameters,the volume fraction of nanofluid,the radiation parameter,the homogeneous reaction parameter,the magnetic parameter,the heterogeneous parameter and the Schmidt number are studied on the heat and flow transfer features.Further,regarding the effective parameters of the present work,the correlation for the Nusselt number has been developed.The outcomes illustrate that with the raising of the heterogeneous parameter and the homogeneous reaction parameter,the concentration profile diminishes.In addition,the outcomes point to a reverse relationship between the Nusselt number and the internal heat source parameters.
文摘The present article aims to investigate the Graetz-Nusselt problem for blood as a non-Newtonian fluid obeying the power-law constitutive equation and flowing inside the axisymmetric tube subjected to nonuniform surface heat flux.After the flow field is determined by solving the continuity and the momentum equations,the energy equation is handled by employing the separation of variables method.The resulting Eigen functions and Eigen values are numerically calculated using MATLAB built-in solver BVP4C.The analysis is first conducted for the situation of constant heat flux and subsequently generalized to apply to the case of sinusoidal variation of wall heat flux along the tube length,using Duhamel’s Theorem.Furthermore,an approximate analytic solution is determined,employing an integral approach to solve the boundary layer equations.With respect to the comparison,the results of approximate solution display acceptable congruence with those of exact solution with an average error of 7.4%.Interestingly,with decreasing the power-law index,the discrepancy between the two presented methods significantly reduces.Eventually,the influences of the controlling parameters such as surface heat flux and power-law index on the non-Newtonian fluid flow’s thermal characteristics and structure are elaborately discussed.It is found that switching from constant wall heat flux to non-uniform wall heat flux that sinusoidally varies along the tube length significantly improves the simulation’s accuracy due to the better characterization of the heat transport phenomenon in non-Newtonian fluid flow through the tube.In the presence of sinusoidally varying wall heat flux with an amplitude of 200 W/m 2 and when the power-law index is 0.25,the maximum arterial wall temperature is found to be about 311.56 K.
文摘The aim of this work is to study the effect of non-uniform single and double slot suction/injection into a steady mixed convection boundary layer flow over a vertical cone, while the axis of the cone is inline with the flow. The governing boundary layer equations are transformed into a non-dimensional form by a group of non-similar trans- formations. The resulting coupled non-linear partial differential equations are solved nu- merically by employing the quasi-linearization technique and an implicit finite-difference scheme. Numerical computations are performed for different values of the dimensionless parameters to display the velocity and temperature profiles graphically. Also, numerical results are presented for the skin friction and heat transfer coefficients. Results indicate that the skin friction and heat transfer coefficients increase with non-uniform slot suction, but the effect of non-uniform slot injection is just opposite.
文摘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.