A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function w...A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function with three parameters, a mass conservation law and a concept of turbulent ellipses. A solution to the model was obtained by using a finite difference method and an extrapolation method. Formulas of calculating development index not only before but also after water breaks through an oil well in the condition of two-phase fluid nonlinear flow in the media were derived. An example was discussed. Water saturation distribution was presented. The moving law of drainage front was found. Laws of change of pressure difference with time were recognized. Results show that there is much difference of water saturation distribution between nonlinear flow and linear flow; that drainage front by water moves faster, water breaks through sooner and the index gets worse because of the nonlinear flow; and that dimensionless pressure difference gets larger at the same dimensionless time and difficulty of oil development becomes bigger by the nonlinear flow. Thus, it is necessary that influence of nonlinear flow on development indexes of the oil fields be taken into account. The results provide water-flooding development of the oilfields with scientific basis.展开更多
Non-Newtonan fluid is a kind of fluid whose components of stresstensor aren’t theliear funtions of compoents of the strain rate tensor.Non-Newtonianfluid is beingprocessed in many kinds of modern industry,Stability ...Non-Newtonan fluid is a kind of fluid whose components of stresstensor aren’t theliear funtions of compoents of the strain rate tensor.Non-Newtonianfluid is beingprocessed in many kinds of modern industry,Stability of flows for Non- Newtonianfluid is of important applicatuib,In this article we calculate subcritical thrdshold of flow which oecurs in polymer-processing when the melting substance is driven throughtwo parallel fixed boundaries.展开更多
A Jeffery-Hamel （J-H） flow model of the non-Newtonian fluid type inside a convergent wedge （inclined walls） with a wall friction is derived by a nonlinear ordinary differential equation with appropriate boundary c...A Jeffery-Hamel （J-H） flow model of the non-Newtonian fluid type inside a convergent wedge （inclined walls） with a wall friction is derived by a nonlinear ordinary differential equation with appropriate boundary conditions based on similarity relationships. Unlike the usual power law model, this paper develops nonlinear viscosity based only on a tangential coordinate function due to the radial geometry shape. Two kinds of solutions are developed, i.e., analytical and semi-analytical （numerical） solutions with suitable assumptions. As a result of the parametric examination, it has been found that the Newtonian normalized velocity gradually decreases with the tangential direction progress. Also, an increase in the friction coefficient leads to a decrease in the normalized Newtonian velocity profile values. However, an increase in the Reynolds number causes an increase in the normalized velocity function values. Additionally, for the small values of wedge semi-angle, the present solutions are in good agreement with the previous results in the literature.展开更多
A fluid–structure interaction method combining a nonlinear finite element algorithm with a preconditioning finite volume method is proposed in this paper to simulate parachute transient dynamics. This method uses a t...A fluid–structure interaction method combining a nonlinear finite element algorithm with a preconditioning finite volume method is proposed in this paper to simulate parachute transient dynamics. This method uses a three-dimensional membrane–cable fabric model to represent a parachute system at a highly folded configuration. The large shape change during parachute inflation is computed by the nonlinear Newton–Raphson iteration and the linear system equation is solved by the generalized minimal residual（GMRES） method. A membrane wrinkling algorithm is also utilized to evaluate the special uniaxial tension state of membrane elements on the parachute canopy. In order to avoid large time expenses during structural nonlinear iteration, the implicit Hilber–Hughes–Taylor（HHT） time integration method is employed. For the fluid dynamic simulations, the Roe and HLLC（Harten–Lax–van Leer contact） scheme has been modified and extended to compute flow problems at all speeds. The lower–upper symmetric Gauss–Seidel（LUSGS） approximate factorization is applied to accelerate the numerical convergence speed. Finally,the test model of a highly folded C-9 parachute is simulated at a prescribed speed and the results show similar characteristics compared with experimental results and previous literature.展开更多
A nonlinear boundary slip model consisting of an initial slip length and a critical shear rate was used to study the nonlinear boundary slip of squeeze fluid film confined between two approaching spheres. It is found ...A nonlinear boundary slip model consisting of an initial slip length and a critical shear rate was used to study the nonlinear boundary slip of squeeze fluid film confined between two approaching spheres. It is found that the initial slip length controls the slip behavior at small shear rate, but the critical shear rate controls the boundary slip at high shear rate. The boundary slip at the squeeze fluid film of spherical surfaces is a strongly nonlinear function of the radius coordinate. At the center or far from the center of the squeeze film, the slip length equals the initial slip length due to the small shear rate. However, in the high shear rate regime the slip length increases very much. The hydrodynamic force of the spherical squeeze film decreases with increasing the initial slip length and decreasing the critical shear rate. The effect of initial slip length on the hydrodynamic force seems less than that of the critical shear rate. When the critical shear rate is very small the hydrodynamic force increases very slowly with a decrease in minimum film thickness. The theoretical predictions agree well with the experiment measurements.展开更多
This work is focused on the effects of heat source/sink, viscous dissipation, radiation and work done by deformation on flow and heat transfer of a viscoelastic fluid over a nonlinear stretching sheet. The similarity ...This work is focused on the effects of heat source/sink, viscous dissipation, radiation and work done by deformation on flow and heat transfer of a viscoelastic fluid over a nonlinear stretching sheet. The similarity transformations have been used to convert the governing partial differential equations into a set of nonlinear ordinary differential equations. These equations are then solved numerically using a very efficient implicit finite difference method. Favorable comparison with previously published work is performed and it is found to be in excellent agreement. The results of this parametric study are shown in several plots and tables and the physical aspects of the problem are highlighted and discussed.展开更多
In this paper, He’s variational iteration method is successfully employed to solve a nonlinear boundary value problem arising in the study of thin film flow of a third grade fluid down an inclined plane. For comparis...In this paper, He’s variational iteration method is successfully employed to solve a nonlinear boundary value problem arising in the study of thin film flow of a third grade fluid down an inclined plane. For comparison, the same problem is solved by the Adomian decomposition method. The results show that the difference between the two solutions is negligible. The conclusion is that this technique may be considered an alternative and efficient method for finding approximate solutions of both linear and nonlinear boundary value problems. Furthermore, the variational iteration method has an advantage over the decomposition method in that it solves the nonlinear problems without using the Adomian polynomials.展开更多
A numerical method for simulating nonlinear fluid-rigid structure interaction problems is developed. The structure is assumed to undergo large rigid body motions and the fluid flow is governed by nonlinear, viscous or...A numerical method for simulating nonlinear fluid-rigid structure interaction problems is developed. The structure is assumed to undergo large rigid body motions and the fluid flow is governed by nonlinear, viscous or non-viscous, field equations with nonlinear boundary conditions applied to the free surface and fluid-solid interaction interfaces. An Arbitrary-LagrangianEulerian (ALE) mesh system is used to construct the numerical model. A multi-block numerical scheme of study is adopted allowing for the relative motion between moving overset grids, which are independent of one another. This provides a convenient method to overcome the difficulties in matching fluid meshes with large solid motions. Nonlinear numerical equations describing nonlinear fluid-solid interaction dynamics are derived through a numerical discretization scheme of study. A coupling iteration process is used to solve these numerical equations. Numerical examples are presented to demonstrate applications of the model developed.展开更多
Nonlinearly dynamic stability of flexible liquid-conveying pipe in fluid structure interaction was analyzed by using modal disassembling technique. The effects of Poisson, Junction and Friction couplings in the wave-f...Nonlinearly dynamic stability of flexible liquid-conveying pipe in fluid structure interaction was analyzed by using modal disassembling technique. The effects of Poisson, Junction and Friction couplings in the wave-flowing-vibration system on the pipe dynamic stability were included in the analytical model constituted by four nonlinear differential equations. An analyzing example of cantilevered pipe was done to illustrate the dynamic stability,characteristics on the pipe in the full coupling mechanisms, and the phase curves related to the first four modal motions were drawn. The results show that the dynamic stable characteristics of the pipe are very complicated in the complete coupling mechanisms, and the kinds of the singularity points corresponding to the various modal motions are different.展开更多
In order to simulate and analyze the dynamic characteristics of the parachute from advanced tactical parachute system(ATPS),a nonlinear finite element algorithm and a preconditioning finite volume method are employed ...In order to simulate and analyze the dynamic characteristics of the parachute from advanced tactical parachute system(ATPS),a nonlinear finite element algorithm and a preconditioning finite volume method are employed and developed to construct three dimensional parachute fluid-structure interaction(FSI)model.Parachute fabric material is represented by membrane-cable elements,and geometrical nonlinear algorithm is employed with wrinkling technique embedded to simulate the large deformations of parachute structure by applying the NewtonRaphson iteration method.On the other hand,the time-dependent flow surrounding parachute canopy is simulated using preconditioned lower-upper symmetric Gauss-Seidel(LU-SGS)method.The pseudo solid dynamic mesh algorithm is employed to update the flow-field mesh based on the complex and arbitrary motion of parachute canopy.Due to the large amount of computation during the FSI simulation,massage passing interface(MPI)parallel computation technique is used for all those three modules to improve the performance of the FSI code.The FSI method is tested to simulate one kind of ATPS parachutes to predict the parachute configuration and anticipate the parachute descent speeds.The comparison of results between the proposed method and those in literatures demonstrates the method to be a useful tool for parachute designers.展开更多
Waves of finite amplitude on a thin layer of non-Newtonian fluid modelled as a power-law fluid are considered. In the long wave approximation, the system of equations taking into account the viscous and nonlinear effe...Waves of finite amplitude on a thin layer of non-Newtonian fluid modelled as a power-law fluid are considered. In the long wave approximation, the system of equations taking into account the viscous and nonlinear effects has the hyper- bolic type. For the two-parameter family of periodic waves in the film flow on a vertical wall the modulation equations for nonlinear wave trains are derived and investigated. The stability criterium for roll waves based on the hyperbolicity of the modulation equations is suggested. It is shown that the evolution of stable roll waves can be described by self-similar solutions of the modulation equations.展开更多
文摘A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function with three parameters, a mass conservation law and a concept of turbulent ellipses. A solution to the model was obtained by using a finite difference method and an extrapolation method. Formulas of calculating development index not only before but also after water breaks through an oil well in the condition of two-phase fluid nonlinear flow in the media were derived. An example was discussed. Water saturation distribution was presented. The moving law of drainage front was found. Laws of change of pressure difference with time were recognized. Results show that there is much difference of water saturation distribution between nonlinear flow and linear flow; that drainage front by water moves faster, water breaks through sooner and the index gets worse because of the nonlinear flow; and that dimensionless pressure difference gets larger at the same dimensionless time and difficulty of oil development becomes bigger by the nonlinear flow. Thus, it is necessary that influence of nonlinear flow on development indexes of the oil fields be taken into account. The results provide water-flooding development of the oilfields with scientific basis.
文摘Non-Newtonan fluid is a kind of fluid whose components of stresstensor aren’t theliear funtions of compoents of the strain rate tensor.Non-Newtonianfluid is beingprocessed in many kinds of modern industry,Stability of flows for Non- Newtonianfluid is of important applicatuib,In this article we calculate subcritical thrdshold of flow which oecurs in polymer-processing when the melting substance is driven throughtwo parallel fixed boundaries.
文摘A Jeffery-Hamel （J-H） flow model of the non-Newtonian fluid type inside a convergent wedge （inclined walls） with a wall friction is derived by a nonlinear ordinary differential equation with appropriate boundary conditions based on similarity relationships. Unlike the usual power law model, this paper develops nonlinear viscosity based only on a tangential coordinate function due to the radial geometry shape. Two kinds of solutions are developed, i.e., analytical and semi-analytical （numerical） solutions with suitable assumptions. As a result of the parametric examination, it has been found that the Newtonian normalized velocity gradually decreases with the tangential direction progress. Also, an increase in the friction coefficient leads to a decrease in the normalized Newtonian velocity profile values. However, an increase in the Reynolds number causes an increase in the normalized velocity function values. Additionally, for the small values of wedge semi-angle, the present solutions are in good agreement with the previous results in the literature.
文摘A fluid–structure interaction method combining a nonlinear finite element algorithm with a preconditioning finite volume method is proposed in this paper to simulate parachute transient dynamics. This method uses a three-dimensional membrane–cable fabric model to represent a parachute system at a highly folded configuration. The large shape change during parachute inflation is computed by the nonlinear Newton–Raphson iteration and the linear system equation is solved by the generalized minimal residual（GMRES） method. A membrane wrinkling algorithm is also utilized to evaluate the special uniaxial tension state of membrane elements on the parachute canopy. In order to avoid large time expenses during structural nonlinear iteration, the implicit Hilber–Hughes–Taylor（HHT） time integration method is employed. For the fluid dynamic simulations, the Roe and HLLC（Harten–Lax–van Leer contact） scheme has been modified and extended to compute flow problems at all speeds. The lower–upper symmetric Gauss–Seidel（LUSGS） approximate factorization is applied to accelerate the numerical convergence speed. Finally,the test model of a highly folded C-9 parachute is simulated at a prescribed speed and the results show similar characteristics compared with experimental results and previous literature.
基金Project supported by the National Natural Science Foundation of China (Nos.10332010, 10272028 and 10421002)the Ph. D. Programs Foundation of Ministry of Education of China (No.2003141013)
文摘A nonlinear boundary slip model consisting of an initial slip length and a critical shear rate was used to study the nonlinear boundary slip of squeeze fluid film confined between two approaching spheres. It is found that the initial slip length controls the slip behavior at small shear rate, but the critical shear rate controls the boundary slip at high shear rate. The boundary slip at the squeeze fluid film of spherical surfaces is a strongly nonlinear function of the radius coordinate. At the center or far from the center of the squeeze film, the slip length equals the initial slip length due to the small shear rate. However, in the high shear rate regime the slip length increases very much. The hydrodynamic force of the spherical squeeze film decreases with increasing the initial slip length and decreasing the critical shear rate. The effect of initial slip length on the hydrodynamic force seems less than that of the critical shear rate. When the critical shear rate is very small the hydrodynamic force increases very slowly with a decrease in minimum film thickness. The theoretical predictions agree well with the experiment measurements.
文摘This work is focused on the effects of heat source/sink, viscous dissipation, radiation and work done by deformation on flow and heat transfer of a viscoelastic fluid over a nonlinear stretching sheet. The similarity transformations have been used to convert the governing partial differential equations into a set of nonlinear ordinary differential equations. These equations are then solved numerically using a very efficient implicit finite difference method. Favorable comparison with previously published work is performed and it is found to be in excellent agreement. The results of this parametric study are shown in several plots and tables and the physical aspects of the problem are highlighted and discussed.
文摘In this paper, He’s variational iteration method is successfully employed to solve a nonlinear boundary value problem arising in the study of thin film flow of a third grade fluid down an inclined plane. For comparison, the same problem is solved by the Adomian decomposition method. The results show that the difference between the two solutions is negligible. The conclusion is that this technique may be considered an alternative and efficient method for finding approximate solutions of both linear and nonlinear boundary value problems. Furthermore, the variational iteration method has an advantage over the decomposition method in that it solves the nonlinear problems without using the Adomian polynomials.
文摘A numerical method for simulating nonlinear fluid-rigid structure interaction problems is developed. The structure is assumed to undergo large rigid body motions and the fluid flow is governed by nonlinear, viscous or non-viscous, field equations with nonlinear boundary conditions applied to the free surface and fluid-solid interaction interfaces. An Arbitrary-LagrangianEulerian (ALE) mesh system is used to construct the numerical model. A multi-block numerical scheme of study is adopted allowing for the relative motion between moving overset grids, which are independent of one another. This provides a convenient method to overcome the difficulties in matching fluid meshes with large solid motions. Nonlinear numerical equations describing nonlinear fluid-solid interaction dynamics are derived through a numerical discretization scheme of study. A coupling iteration process is used to solve these numerical equations. Numerical examples are presented to demonstrate applications of the model developed.
基金Foundation items:the National Natural Science Foundation of China(50079007)the Hydraulic Science Foundation of China Hydraulic Ministry(SZ9830)the Natural Science Foundation of Yunnan Province(98E003G)
文摘Nonlinearly dynamic stability of flexible liquid-conveying pipe in fluid structure interaction was analyzed by using modal disassembling technique. The effects of Poisson, Junction and Friction couplings in the wave-flowing-vibration system on the pipe dynamic stability were included in the analytical model constituted by four nonlinear differential equations. An analyzing example of cantilevered pipe was done to illustrate the dynamic stability,characteristics on the pipe in the full coupling mechanisms, and the phase curves related to the first four modal motions were drawn. The results show that the dynamic stable characteristics of the pipe are very complicated in the complete coupling mechanisms, and the kinds of the singularity points corresponding to the various modal motions are different.
文摘In order to simulate and analyze the dynamic characteristics of the parachute from advanced tactical parachute system(ATPS),a nonlinear finite element algorithm and a preconditioning finite volume method are employed and developed to construct three dimensional parachute fluid-structure interaction(FSI)model.Parachute fabric material is represented by membrane-cable elements,and geometrical nonlinear algorithm is employed with wrinkling technique embedded to simulate the large deformations of parachute structure by applying the NewtonRaphson iteration method.On the other hand,the time-dependent flow surrounding parachute canopy is simulated using preconditioned lower-upper symmetric Gauss-Seidel(LU-SGS)method.The pseudo solid dynamic mesh algorithm is employed to update the flow-field mesh based on the complex and arbitrary motion of parachute canopy.Due to the large amount of computation during the FSI simulation,massage passing interface(MPI)parallel computation technique is used for all those three modules to improve the performance of the FSI code.The FSI method is tested to simulate one kind of ATPS parachutes to predict the parachute configuration and anticipate the parachute descent speeds.The comparison of results between the proposed method and those in literatures demonstrates the method to be a useful tool for parachute designers.
文摘Waves of finite amplitude on a thin layer of non-Newtonian fluid modelled as a power-law fluid are considered. In the long wave approximation, the system of equations taking into account the viscous and nonlinear effects has the hyper- bolic type. For the two-parameter family of periodic waves in the film flow on a vertical wall the modulation equations for nonlinear wave trains are derived and investigated. The stability criterium for roll waves based on the hyperbolicity of the modulation equations is suggested. It is shown that the evolution of stable roll waves can be described by self-similar solutions of the modulation equations.
