A relatively high formation pressure gradient can exist in seepage flow in low-permeable porous media with a threshold pressure gradient, and a significant error may then be caused in the model computation by neglecti...A relatively high formation pressure gradient can exist in seepage flow in low-permeable porous media with a threshold pressure gradient, and a significant error may then be caused in the model computation by neglecting the quadratic pressure gradient term in the governing equations. Based on these concerns, in consideration of the quadratic pressure gradient term, a basic moving boundary model is constructed for a one-dimensional seepage flow problem with a threshold pressure gradient. Owing to a strong nonlinearity and the existing moving boundary in the mathematical model, a corresponding numerical solution method is presented. First, a spatial coordinate transformation method is adopted in order to transform the system of partial differential equa- tions with moving boundary conditions into a closed system with fixed boundary conditions; then the solution can be sta- bly numerically obtained by a fully implicit finite-difference method. The validity of the numerical method is verified by a published exact analytical solution. Furthermore, to compare with Darcy's flow problem, the exact analytical solution for the case of Darcy's flow considering the quadratic pressure gradient term is also derived by an inverse Laplace transform. A comparison of these model solutions leads to the conclu- sion that such moving boundary problems must incorporate the quadratic pressure gradient term in their governing equa- tions; the sensitive effects of the quadratic pressure gradient term tend to diminish, with the dimensionless threshold pres- sure gradient increasing for the one-dimensional problem.展开更多
基金funding by the project (Grant 51404232) sponsored by the National Natural Science Foundation of Chinathe National Science and Technology Major Project (Grant 2011ZX05038003)+1 种基金the China Postdoctoral Science Foundation project (Grant 2014M561074)the China Scholarship Council for its generous financial support of the research
文摘A relatively high formation pressure gradient can exist in seepage flow in low-permeable porous media with a threshold pressure gradient, and a significant error may then be caused in the model computation by neglecting the quadratic pressure gradient term in the governing equations. Based on these concerns, in consideration of the quadratic pressure gradient term, a basic moving boundary model is constructed for a one-dimensional seepage flow problem with a threshold pressure gradient. Owing to a strong nonlinearity and the existing moving boundary in the mathematical model, a corresponding numerical solution method is presented. First, a spatial coordinate transformation method is adopted in order to transform the system of partial differential equa- tions with moving boundary conditions into a closed system with fixed boundary conditions; then the solution can be sta- bly numerically obtained by a fully implicit finite-difference method. The validity of the numerical method is verified by a published exact analytical solution. Furthermore, to compare with Darcy's flow problem, the exact analytical solution for the case of Darcy's flow considering the quadratic pressure gradient term is also derived by an inverse Laplace transform. A comparison of these model solutions leads to the conclu- sion that such moving boundary problems must incorporate the quadratic pressure gradient term in their governing equa- tions; the sensitive effects of the quadratic pressure gradient term tend to diminish, with the dimensionless threshold pres- sure gradient increasing for the one-dimensional problem.