共查询到20条相似文献,搜索用时 0 毫秒
1.
Infrared thermography has increasingly gained importance because of environmental and technological advancements of this method and is applied in a variety of disciplines related to non-isothermal flow. However, it has not been used so far for quantitative thermal analysis in saturated porous media. This article suggests infrared thermographic approach to obtain the entire surface temperature distribution(s) in water-saturated porous media. For this purpose, infrared thermal analysis is applied with in situ calibration for a better understanding of the heat transfer processes in porous media. Calibration is achieved with a combination of invasive sensors which are inserted into the medium and non-invasive thermal sensors in which sensors are not inserted to measure temperatures but it works through the detection of infrared radiation emitted from the surface. Thermocouples of relatively thin diameter are used to minimize the disturbance for flow. Thermocouples give the temperature values at specified positions inside the porous medium, and these values are compared with the values suggested by the infrared thermographic device at the same positions, in the calibration exercise. The calibration process was repeated for different temperatures and flow rates to get the temperature distributions of the whole material inside the system. This technique enables us to measure accurate two-dimensional temperature distributions, which is not possible by using thermocouples only. Continuous point heat sources at different flow rates and temperatures are studied experimentally. Additionally, it offers numerical simulations of the experiments utilizing a finite element-based model. A two-dimensional density and viscosity-dependent flow and transport model accounting for thermal dispersion is utilized to simulate the experimental results. Possible small heat losses from the surface are incorporated in the model according to the properties and thickness of the Plexiglass material used for the construction of the experiment tank. The numerical results agree well with the experimental observations. 相似文献
2.
Pore-fluid flow associated with seepage instabilities can play an important role in controlling large mineralisation patterns within the upper crust of the Earth. To demonstrate this process, two kinds of seepage instability problems in fluid-saturated porous rocks are considered in this paper. The first kind of seepage instability problem is caused by the temperature-induced buoyancy of pore fluid, so that it can be called the buoyancy-driven seepage instability problem, while the second kind of seepage instability problem is caused by chemical dissolution reactions that are commonly encountered in the upper crust of the Earth, so that it can be called the chemical-dissolution-driven seepage instability problem. After the mathematical governing equations of and computational methods for these two kinds of seepage instability problems are introduced, two numerical examples are used to elucidate how and why these two kinds of seepage instabilities can provide favorable places for the formation of large mineralisation patterns within the upper crust of the Earth. The related computational simulation results have demonstrated that: (1) the convective pore-fluid flow caused by the buoyancy-driven seepage instability not only can dissolve minerals at the lower part of the upper crust, but also can transport the dissolved minerals from the lower part to the upper part of the upper crust, resulting in large mineralisation patterns near the surface of the Earth's upper crust. (2) The chemical-dissolution-driven seepage instability in fluid-saturated porous rock can provide some favorable places, such as finger-like channels created by porosity enhancement in the porous rock, for the formation of large mineralisation patterns within the upper crust of the Earth. 相似文献
3.
Chongbin Zhao Lynn B. Reid Klaus Regenauer-Lieb Thomas Poulet 《Computational Geosciences》2012,16(3):735-755
In dealing with chemical-dissolution-front propagation problems in fluid-saturated porous media, the chemical dissolution front represented by the porosity of the medium may have a very steep slope (i.e., a very large porosity gradient) at the dissolution front, depending on the mineral dissolution ratio that is defined as the equilibrium concentration of the dissolved minerals in the pore-fluid to the solid molar density of the dissolvable minerals in the solid matrix. When the mineral dissolution ratio approaches zero, the theoretical value of the porosity gradient tends to infinity at the chemical dissolution front. Even for a very small value of the mineral dissolution ratio, which is very common in geochemical systems, the porosity gradient can be large enough to cause the solution hard to converge when the conventional finite element method is used to solve a chemical dissolution problem in a fluid-saturated porous medium where the pore-fluid is compressible. To improve the convergent speed of solution, a porosity-gradient replacement approach, in which the term involving porosity-gradient computation is replaced by a new term consisting of pore-fluid density-gradient and pressure-gradient computation, is first proposed and then incorporated into the finite element method in this study. Through comparing the numerical results obtained from the proposed approach with the theoretical solutions for a benchmark problem, it has been demonstrated that not only can the solution divergence be avoid, but also the accurate simulation results can be obtained when the proposed porosity-gradient replacement approach is used to solve chemical-dissolution-front propagation problems in fluid-saturated porous media including pore-fluid compressibility. 相似文献
4.
5.
6.
This paper deals with the theoretical aspects of nonaqueous phase liquid (NAPL)‐dissolution‐induced instability in two‐dimensional fluid‐saturated porous media including solute dispersion effects.After some weaknesses associated with the previous work are analyzed and overcome, a comprehensive dimensionless number, known as the Zhao number, is proposed to represent the main driving force and three controlling mechanisms of an NAPL‐dissolution system that has a finite domain. The linear stability analysis is carried out to derive the critical value of the comprehensive dimensionless number of the NAPL‐dissolution system in a limit case as the ratio of the equilibrium concentration to the density of the NAPL approaches zero. As a result, a theoretical criterion that can be used to assess the instability of planar NAPL‐dissolution fronts in two‐dimensional fluid‐saturated porous media of finite domains has been established. Not only can the present theoretical results be used for the theoretical understanding of the effect of solute dispersion on the instability of an NAPL‐dissolution front in the fluid‐saturated porous medium of either a finite domain or an infinite domain, but also they can be used as benchmark solutions for verifying numerical methods employed to simulate detailed morphological evolution processes of NAPL‐dissolution fronts in two‐dimensional fluid‐saturated porous media. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
7.
Ettore Vidotto Rainer Helmig Martin Schneider Barbara Wohlmuth 《Computational Geosciences》2018,22(6):1487-1502
In this paper, we present a fast streamline-based numerical method for the two-phase flow equations in high-rate flooding scenarios for incompressible fluids in heterogeneous and anisotropic porous media. A fractional flow formulation is adopted and a discontinuous Galerkin method (DG) is employed to solve the pressure equation. Capillary effects can be neglected in high-rate flooding scenarios. This allows us to present an improved streamline approach in combination with the one-dimensional front tracking method to solve the transport equation. To handle the high computational costs of the DG approximation, domain decomposition is applied combined with an algebraic multigrid preconditioner to solve the linear system. Special care at the interior interfaces is required and the streamline tracer has to include a dynamic communication strategy. The method is validated in various two- and three-dimensional tests, where comparisons of the solutions in terms of approximation of flow front propagation with standard fully implicit finite-volume methods are provided. 相似文献
8.
Practical civil engineering problems are usually formulated in an infinite half-space domain, and a selected finite domain is required to analyze the dynamic responses of a fluid-saturated porous medium by the finite element method (FEM). Devising a method to deal with the boundaries of the finite domain is the key issue for this open system. In this paper, a two-dimensional spring–dashpot artificial boundary (SDAB) for transient analysis in a fluid-saturated porous media is developed. Based on Biot’s dynamic theory of fluid-saturated porous media, the normal and tangential boundary stress formulae are deduced for out-going cylindrical body waves. The boundary stress is proportional to displacement and velocity, thus continuously distributed dashpots and springs can be placed on the artificial boundaries in the normal and tangential directions to simulate the energy absorption of the infinite media outside of the finite domain for the interior distributed source problems. In this paper, the input seismic motion can be realized by applying an equivalent load on the SDAB for the seismic scattering problems of exterior distributed sources. Numerical examples are given and the analyzed results show that the SDAB and the method of wave motion input have good stability and acceptable accuracy. 相似文献
9.
A mixed finite element–boundary element solution for the analysis of two-dimensional flow in porous media composed of rock blocks and discrete fractures is described. The rock blocks are modelled implicitly by using boundary elements whereas finite elements are adopted to model the discrete fractures. The computational procedure has been implemented in a hybrid code which has been validated first by comparing the numerical results with the closed-form solution for flow in a porous aquifer intercepted by a vertical fracture only. Then, a more complex problem has been solved where a pervious, homogeneous and isotropic matrix containing a net of fractures is considered. The results obtained are shown to describe satisfactorily the main features of the flow problem under study. © 1997 by John Wiley & Sons, Ltd. 相似文献
10.
饱和土半空间中圆柱形孔洞对平面P波的散射 总被引:1,自引:0,他引:1
在 Biot饱和多孔介质动力学理论的基础上,首次建立了求解饱和土半空间中圆柱形孔洞对平面 P 波散射问题的波函数展开法。首先。分析了具有圆柱形孔洞的饱和土半空间场地在平面 P 波入射下产生散射波系,并将入射波和散射波的波函数在圆柱坐标下展开。然后,引入边界条件,求出散射波函数的待定系数,从而,得到饱和土半空间中圆柱形孔洞对平面P波的散射问题的解析解。根据所得的波函数的解,可求解区域内的位移、应力的值,同时,分析了入射波频率、入射角对柱面上的应力集中因子的影响。 相似文献
11.
12.
单层不可压缩饱和多孔介质一维瞬态响应精确解 总被引:1,自引:0,他引:1
基于Biot理论,考虑惯性、黏滞和机械耦合作用,假定固体颗粒和流体均不可压缩,得到了表面任意竖向荷载作用下单层饱和多孔介质一维瞬态响应的精确解。导出了以固体骨架位移表示的无量纲控制方程,并将边界条件齐次化。求解对应无黏滞耦合作用的特征值问题,得到一组满足齐次边界条件、关于空间坐标的正交函数基。利用变异系数法和基函数的正交性,得到一系列相互解耦的、关于时间的二阶常微分方程及相应的初始条件,并采用状态空间法求解常微分方程,得到位移分量。对整体平衡方程关于空间坐标积分,根据边界条件可确定总应力,并进而求得孔隙压力。通过算例验证所得解法的正确性 相似文献
13.
In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton–Rogers–Lapwood problem in a fluid-saturated porous medium. The Horton–Rogers–Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton–Rogers–Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton–Rogers–Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton–Rogers–Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. © 1997 by John Wiley & Sons, Ltd. 相似文献
14.
Under certain physically reasonable assumptions, three-phase flow of immiscible, incompressible fluids can be described by a 2×2 nongenuinely nonlinear, hyperbolic system. We combine analytical solutions to the corresponding Riemann problem with an efficient front-tracking method to study Cauchy and initial-boundary value problems. Unlike finite difference methods, the front-tracking method treats all waves as discontinuities by evolving shocks exactly and approximating rarefactions by small entropy-violating discontinuities. This way, the method can track individual waves and give very accurate (or even exact) resolution of discontinuities. We demonstrate the applicability of the method through several numerical examples, including a streamline simulation of a water-alternating-gas (WAG) injection process in a three-dimensional, heterogeneous, shallow-marine formation. 相似文献
15.
重非水相污染物(DNAPL)在地下介质中运移和分布受多种因素控制,包括DNAPL本身的物理化学性质,土的性质,泄漏条件等等。由于介质的非均质性,使得多相流运移行为更为复杂。基于地下水随机理论构建渗透率随机场,采用蒙特卡罗方法探讨泄漏速率对非均质饱和介质中DNAPL运移的影响。数值结果表明,在泄漏总量一定的情况下,泄漏速率越低,介质非均质性对DNAPL运移的影响程度越高。反之,DNAPL的渗漏速率越高,小尺度地层的非均质性影响越低。由于DNAPL运移过程中在垂直方向受重力的影响,污染羽在空间上的质心位置(一阶矩)以及展布范围(二阶矩)在垂直方向上的变异程度要高于水平方向。 相似文献
16.
17.
Conservation equations for mass, momentum, energy, and entropy are formulated for the phases and interfaces of a three‐phase system consisting of a solid and two immiscible fluids. The microscale equations are averaged to the macroscale by integration over a representative elementary volume. Thermodynamic statements for each of the phases and interface entities are also formulated at the microscale and then averaged to the macroscale. This departure from most uses of thermodynamics in macroscale analysis ensures consistency between models and parameters at the two scales. The expressions for the macroscale rates of change of internal energy are obtained by differentiating the derived forms for energy and making use of averaging theorems. These thermodynamic expressions, along with the conservation equations, serve as constraints on the entropy inequality. A linearization of the resulting equations is employed to investigate the theoretical origins of the Biot coefficient that relates the hydrostatic part of the total stress tensor to the normal force applied at the solid surface by the pore fluids. The results here are placed in the context of other formulations and expressions that appear in the literature. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
18.
Variationally consistent homogenization is exploited for the analysis of transient uncoupled consolidation in micro‐heterogeneous porous solids, whereby the classical approach of first‐order homogenization for stationary problems is extended to transient problems. Homogenization is then carried out in the spatial domain on representative volume elements (RVE), which are introduced in quadrature points in standard fashion. Along with the classical averages, a higher‐order conservation quantity is obtained. An iterative FE2‐algorithm is devised for the case of nonlinear permeability and storage coefficients, and it is applied to pore pressure changes in asphalt‐concrete (particle composite). Various parametric studies are carried out, in particular, with respect to the influence of the ‘substructure length scale’ that is represented by the size of the RVE's. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
19.
Three node-centered finite volume discretizations for multiphase porous media flow are presented and compared. By combination
of these methods two additional discretization methods are generated. The ability of these schemes to describe flows at textural
interfaces of different geologic formations is investigated. It was found that models with nonzero-entry pressures for the
capillary pressure-saturation relationship in conjunction with the Box discretization may give rise to spurious oscillations
for flows around low permeable lenses. Furthermore, the applicability and sensitivity of the discretization methods with regard
to the used computational grids is discussed. The schemes are used for the numerical study of two-phase flow in porous media
with zones of different material properties.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
20.
Quantification of the saturation degree of nonaqueous phase liquid in the vadose zone using hydrocarbon partitioning tracers 总被引:1,自引:1,他引:0
Partitioning tracer tests, as an alternative to the core sampling method, were conducted to quantify the degree of saturation
of water and nonaqueous phase liquids (NAPL) in the vadose zone. Hydrocarbon gases, which have less effect on global warming
than conventional tracers, were used as partitioning tracers. Column tests using CH4, C3H5, and C4H10 as non-partitioning and partitioning tracers were performed to determine the retardation factor and partition coefficient
of the tracer into water and NAPL. The retardation factors of these tracers were estimated to be in the range of 1.0–7.0 based
on breakthrough curves of the tracers. The partition coefficient of C3H5 to water and diesel phase was calculated to be 0.57 and 8.45, respectively. For a heavier tracer, C4H10, the partition coefficient to the water and diesel phases was 1.2 and 40.5, respectively. The average value of water and
diesel saturation estimated from column tests agreed well with known values in unsaturated soil. A residence time longer than
7.5 h within soil pores was found to provide local equilibrium partitioning of the tracer to the diesel phase. The concentration
of tracer had no effect on the partitioning process. 相似文献