The distribution of remaining oil is often described qualitatively. The remaining oil distributed in the whole reservoir is calculated according to the characteristics of the space distribution of the saturation of re...The distribution of remaining oil is often described qualitatively. The remaining oil distributed in the whole reservoir is calculated according to the characteristics of the space distribution of the saturation of remaining oil. Logging data are required to accomplish this. However, many such projects cannot be completed. Since the old study of remaining oil distribution could not be quantified efficiently, the ＂dynamic two-step method＂ is presented. Firstly, the water cut of every flow unit in one well at one time is calculated according to the comprehensive water cut of a single well at one time. Secondly, the remaining oil saturation of the flow unit of the well at one time is calculated based on the water cut of the flow unit at a given time. The results show that ＂dynamic two-step method＂ has characteristics of simplicity and convenience, and is especially suitable for the study of remaining oil distribution at high water-cut stage. The distribution of remaining oil presented banding and potato form, remaining oil was relatively concentrated in faultage neighborhood and imperfect well netting position, and the net thickness of the place was great. This proposal can provide an effective way to forecast remaining oil distribution and enhance oil recovery, especially applied at the high water-cut stage.展开更多
A two-dimensional numerical model based on the Navier-Stokes equations and computational Lagrangian-Eulerian advection remap-volume of fluid （CLEAR-VOF） method was developed to simulate wave and flow problems. The N...A two-dimensional numerical model based on the Navier-Stokes equations and computational Lagrangian-Eulerian advection remap-volume of fluid （CLEAR-VOF） method was developed to simulate wave and flow problems. The Navier-Stokes equations were discretized with a three-step finite element method that has a third-order accuracy. In the CLEAR-VOF method, the VOF function F was calculated in the Lagrangian manner and allowed the complicated free surface to be accurately captured. The propagation of regular waves and solitary waves over a flat bottom, and shoaling and breaking of solitary waves on two different slopes were simulated with this model, and the numerical results agreed with experimental data and theoretical solutions. A benchmark test of dam-collapse flow was also simulated with an unstructured mesh, and the capability of the present model for wave and flow simulations with unstructured meshes, was verified. The results show that the model is effective for numerical simulation of wave and flow problems with both structured and unstructured meshes.展开更多
The present contributions deal with the development and implementation of the direct and adjoint incompressible Navie-Stokes equations in a matrix-free time-stepping context making use of the open-source OpenFOAM~ （o...The present contributions deal with the development and implementation of the direct and adjoint incompressible Navie-Stokes equations in a matrix-free time-stepping context making use of the open-source OpenFOAM~ （open source field operation and manipulation） C ＋ ＋ tool- box. It is shown that a few minor modifications in the source code can enable these multi-purpose computational fluid dynamics （CFD） solvers to carry out global instability analysis of three- dimensional flows of engineering interest, exploiting all additional capabilities of such codes. The strength of this methodology is demonstrated through the analysis of an interesting selection of ooen and wall-bounded flows.展开更多
In practical problems, the variation of the free surface around a circular pier is severe. For the Fluent CFD, the multiphase flow models cannot be used together with the LES model to calculate the free surface. Prese...In practical problems, the variation of the free surface around a circular pier is severe. For the Fluent CFD, the multiphase flow models cannot be used together with the LES model to calculate the free surface. Present paper provides a two-step method which joints the 2D compressible ideal-gas equations and the LES model to calculate the 3D flow field with free surface around the pier. The effects of the free surface on the flow structures are studied in detail.展开更多
The simulation of injection molding process requires a stable algorithm to model the molten polymer with non-isothermal non-Newtonian property.In this paper,a staggered and iterative scheme is particularly designed to...The simulation of injection molding process requires a stable algorithm to model the molten polymer with non-isothermal non-Newtonian property.In this paper,a staggered and iterative scheme is particularly designed to solve the velocity-pressure-temperature variables.In consideration of the polymer characteristic of high viscosity and low thermal conductivity,the non-Newtonian momentum-mass conservation equations are solved by the Crank-Nicolson method based split (CNBS) scheme,and the energy conservation equation with convective character is discretized by the characteristic Galerkin (CG) method.In addition,an arbitrary Lagrangian Eulerian (ALE) free surface tracking and mesh generation method is introduced to catch the front of the fluid flow.The efficiency of the proposed scheme is demonstrated by numerical experiments including a lid-driven cavity flow problem and an injection molding problem.展开更多
In the present work, a numerical solution is described for turbulent forced convection flow of an absorbing, emitting, scattering and gray fluid over a two-dimensional backward facing step in a horizontal duct. The AK...In the present work, a numerical solution is described for turbulent forced convection flow of an absorbing, emitting, scattering and gray fluid over a two-dimensional backward facing step in a horizontal duct. The AKN low-Reynolds-number model is employed to predict turbulent flows with separation and heat transfer, while the radiation part of the problem is modeled by the discrete ordinate method (DOM). Discretized forms of the governing equations for fluid flow are obtained by finite volume approach and solved using SIMPLE algorithm. Results are presented for the distributions of Nusselt numbers as a function of the controlling parameters like radiation-conduction parameter (RC) and optical thickness.展开更多
针对油浸式电力变压器瞬态温升计算效率过低的问题,该文提出本征正交分解-αATS(proper orthogonal decomposition-adaptive time stepping based onαfactor,POD-αATS)降阶自适应变步长瞬态计算方法。首先,推导变压器绕组瞬态温升计...针对油浸式电力变压器瞬态温升计算效率过低的问题,该文提出本征正交分解-αATS(proper orthogonal decomposition-adaptive time stepping based onαfactor,POD-αATS)降阶自适应变步长瞬态计算方法。首先,推导变压器绕组瞬态温升计算的有限元离散方程;其次,采用POD降阶算法改善传统瞬态计算中存在的条件数过大及方程阶数过高的问题;同时对于瞬态计算中的时间步长选择问题,提出适用于非线性问题的αATS变步长策略;然后,为验证方法的有效性,基于110 kV油浸式电力变压器绕组的基本结构建立二维八分区数值计算模型,同时将计算结果与基于110 kV绕组的温升实验结果进行对比。数值计算及实验结果表明,所提算法与全阶定步长算法在流场和温度场中的精度几乎相同,且流场计算效率提升约45倍,温度场计算效率提升约38倍,计算速度得到显著提高。这一点在温升实验中同样得到验证,说明该文所提算法的准确性、高效性及一定的工程实用性。展开更多
文摘The distribution of remaining oil is often described qualitatively. The remaining oil distributed in the whole reservoir is calculated according to the characteristics of the space distribution of the saturation of remaining oil. Logging data are required to accomplish this. However, many such projects cannot be completed. Since the old study of remaining oil distribution could not be quantified efficiently, the ＂dynamic two-step method＂ is presented. Firstly, the water cut of every flow unit in one well at one time is calculated according to the comprehensive water cut of a single well at one time. Secondly, the remaining oil saturation of the flow unit of the well at one time is calculated based on the water cut of the flow unit at a given time. The results show that ＂dynamic two-step method＂ has characteristics of simplicity and convenience, and is especially suitable for the study of remaining oil distribution at high water-cut stage. The distribution of remaining oil presented banding and potato form, remaining oil was relatively concentrated in faultage neighborhood and imperfect well netting position, and the net thickness of the place was great. This proposal can provide an effective way to forecast remaining oil distribution and enhance oil recovery, especially applied at the high water-cut stage.
基金supported by the National Natural Science Foundation of China (Grant No. 50679008)
文摘A two-dimensional numerical model based on the Navier-Stokes equations and computational Lagrangian-Eulerian advection remap-volume of fluid （CLEAR-VOF） method was developed to simulate wave and flow problems. The Navier-Stokes equations were discretized with a three-step finite element method that has a third-order accuracy. In the CLEAR-VOF method, the VOF function F was calculated in the Lagrangian manner and allowed the complicated free surface to be accurately captured. The propagation of regular waves and solitary waves over a flat bottom, and shoaling and breaking of solitary waves on two different slopes were simulated with this model, and the numerical results agreed with experimental data and theoretical solutions. A benchmark test of dam-collapse flow was also simulated with an unstructured mesh, and the capability of the present model for wave and flow simulations with unstructured meshes, was verified. The results show that the model is effective for numerical simulation of wave and flow problems with both structured and unstructured meshes.
基金Support of the Marie Curie Grant PIRSES-GA-2009-247651‘‘FP7-PEOPLE-IRSES:ICOMASEF-Instability and Control of Massively Separated Flows”
文摘The present contributions deal with the development and implementation of the direct and adjoint incompressible Navie-Stokes equations in a matrix-free time-stepping context making use of the open-source OpenFOAM~ （open source field operation and manipulation） C ＋ ＋ tool- box. It is shown that a few minor modifications in the source code can enable these multi-purpose computational fluid dynamics （CFD） solvers to carry out global instability analysis of three- dimensional flows of engineering interest, exploiting all additional capabilities of such codes. The strength of this methodology is demonstrated through the analysis of an interesting selection of ooen and wall-bounded flows.
文摘In practical problems, the variation of the free surface around a circular pier is severe. For the Fluent CFD, the multiphase flow models cannot be used together with the LES model to calculate the free surface. Present paper provides a two-step method which joints the 2D compressible ideal-gas equations and the LES model to calculate the 3D flow field with free surface around the pier. The effects of the free surface on the flow structures are studied in detail.
基金the National Natural Science Foundation of China(No. 50873060)the YuYao Technology Division Grand Science and Technology Special Project
文摘The simulation of injection molding process requires a stable algorithm to model the molten polymer with non-isothermal non-Newtonian property.In this paper,a staggered and iterative scheme is particularly designed to solve the velocity-pressure-temperature variables.In consideration of the polymer characteristic of high viscosity and low thermal conductivity,the non-Newtonian momentum-mass conservation equations are solved by the Crank-Nicolson method based split (CNBS) scheme,and the energy conservation equation with convective character is discretized by the characteristic Galerkin (CG) method.In addition,an arbitrary Lagrangian Eulerian (ALE) free surface tracking and mesh generation method is introduced to catch the front of the fluid flow.The efficiency of the proposed scheme is demonstrated by numerical experiments including a lid-driven cavity flow problem and an injection molding problem.
文摘In the present work, a numerical solution is described for turbulent forced convection flow of an absorbing, emitting, scattering and gray fluid over a two-dimensional backward facing step in a horizontal duct. The AKN low-Reynolds-number model is employed to predict turbulent flows with separation and heat transfer, while the radiation part of the problem is modeled by the discrete ordinate method (DOM). Discretized forms of the governing equations for fluid flow are obtained by finite volume approach and solved using SIMPLE algorithm. Results are presented for the distributions of Nusselt numbers as a function of the controlling parameters like radiation-conduction parameter (RC) and optical thickness.
文摘针对油浸式电力变压器瞬态温升计算效率过低的问题,该文提出本征正交分解-αATS(proper orthogonal decomposition-adaptive time stepping based onαfactor,POD-αATS)降阶自适应变步长瞬态计算方法。首先,推导变压器绕组瞬态温升计算的有限元离散方程;其次,采用POD降阶算法改善传统瞬态计算中存在的条件数过大及方程阶数过高的问题;同时对于瞬态计算中的时间步长选择问题,提出适用于非线性问题的αATS变步长策略;然后,为验证方法的有效性,基于110 kV油浸式电力变压器绕组的基本结构建立二维八分区数值计算模型,同时将计算结果与基于110 kV绕组的温升实验结果进行对比。数值计算及实验结果表明,所提算法与全阶定步长算法在流场和温度场中的精度几乎相同,且流场计算效率提升约45倍,温度场计算效率提升约38倍,计算速度得到显著提高。这一点在温升实验中同样得到验证,说明该文所提算法的准确性、高效性及一定的工程实用性。
