Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability简

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.     

弹性地基上多孔功能梯度材料矩形板的自由振动与临界屈曲载荷分析

多孔功能梯度材料(FGM)构件的特性与孔隙率和孔隙分布形式有密切关系。本文基于经典板理论,考虑不同孔隙分布形式时修正的混合率模型,研究Winkler弹性地基上四边受压多孔FGM矩形板的自由振动与临界屈曲载荷特性。首先利用Hamilton原理和物理中面的定义推导Winkler弹性地基上四边受压多孔FGM矩形板自由振动的控制微分方程并进行无量纲化,然后应用微分变换法(DTM)对无量纲控制微分方程和边界条件进行变换,得到计算无量纲固有频率和临界屈曲载荷的代数特征方程。将问题退化为孔隙率为零时的FGM矩形板并与已有文献进行对比以验证其有效性。最后计算并分析了梯度指数、孔隙率、地基刚度系数、长宽比、四边受压载荷及边界条件对多孔FGM矩形板无量纲固有频率的影响以及各参数对无量纲临界屈曲载荷的影响。     

The effect of radial diffusion on the performance of a liquid-liquid displacement process

The effect of radial diffusion on the performance of a liquid-liquid displacement process is considered in fluid flow between porous parallel plates and through a porous tube, as examples of a two-zone problem in unsteady-state mass transfer. The double Laplace transformation is applied to the system equations. In obtaining the inversion of the Laplace transformed equations the first inversion (with respect to the transformed dimensionless axial distance) is performed by use of the residue method, and then the second inversion (with respect to the transformed dimensionless time) is performed by use of the numerical Laplace transform technique advanced by Bellman et al. A numerical example is shown and discussed.     

Influence of diffusion and kinetic characteristics on the combustion of a carbon particle

The spherically symmetric steady-state problem of heterogeneous combustion for a carbon particle is discussed. The exact solution of this problem is obtained. The dimensionless mass flux through the particle surface is graphically represented as a function of oxygen concentration in the external domain. The rate ratio for the diffusion process and the kinetic process is introduced as a dimensionless characteristic of the combustion process. The dimensionless mass fluxes through the particle surface are graphically represented as functions of external constitutive parameters. The resulting solution is approximated. The diffusion and kinetic modes of combustion are considered as limiting cases.     

Unsteady convective mass transfer in a three-zone system

A three zone unsteady state mass transfer problem is considered for flow between parallel plates. Solute is transferred from the moving zone to the other two, in one of which a first order reaction consuming the solute occurs. (An industrially important example of this situation would be the absorption of a gas from a lean mixture into a liquid-saturated ion exchange membrane.) The double Laplace transformation is applied to the system equations. In obtaining the inversion of the transformed equations the first inversion (with respect to the transformed dimensionless axial distance) is performed by use of the residue method, and then the second inversion (with respect to the transformed dimensionless time) is performed by use of the numerical technique developed by Bellman et al. Some numerical results are presented and discussed.     

Fountain flow of pseudoplastic and viscoplastic fluids

Numerical simulations have been undertaken for the benchmark problem of fountain flow present in injection-mold filling. The Finite Element Method (FEM) is used to provide numerical results for both cases of planar and axisymmetric domains under steady-state conditions. The Herschel–Bulkley model of viscoplasticity is used, which reduces with appropriate modifications to the Bingham, power-law and Newtonian models. The present results extend previous ones regarding the shape of the front, which is essential in correctly capturing the flow field. In particular the centreline front position is found as a function of the dimensionless power-law index (in the case of pseudoplasticity) and the dimensionless yield stress (in the case of viscoplasticity). The pressures from the simulations have been used to compute the excess pressure losses in the system (front pressure correction or exit correction). Both shear-thinning and shear-thickening lead to more extended front positions relative to the Newtonian values, which are 0.895 for the planar case and 0.835 for the axisymmetric one. Viscoplasticity leads also to more extended front positions as the dimensionless yield stress goes from zero (Newtonian behaviour) to higher values of the yield stress. In both cases of non-Newtonian behaviour, the front tends to follow the development of the fully developed Poiseuille velocity profile, which tends towards a plug-like profile at the extreme cases of non-Newtonianness. The front pressure (exit) correction increases monotonically with the decrease in the power-law index and the increase in the dimensionless yield stress.     

Solution of a Plane Hydrofracture Problem with Stress Contrast

A plane hydrofracture problem for the Khristianovich–Geertsma–de Klerk model is extended and solved in the case where a confining stress closing a fracture is not constant in the direction of its propagation. A method is developed for solving the problem with an arbitrary stress contrast. It is stated that the transition through a contact with positive (negative) contrast occurs with fracture arresting (acceleration), whose intensity is controlled by a dimensionless parameter derived from theoretical considerations and numerical results.     

Fluid flow simulation in a random elliptical porous medium by using the lattice Boltzmann method

A fluid flow through an isotropic porous medium with randomly arranged elliptical particles is simulated by the lattice Boltzmann method. The dimensionless pressure drop and the dimensionless permeability are evaluated as functions of the Reynolds number. The effect of the aspect ratio of the major to minor semi-axis of the ellipse on the dimensionless permeability is considered for different values of porosity. The pressure drop is thoroughly investigated as a function of fluid viscosity for different values of the aspect ratio and porosity. The influence of various parameters of the problem on the mean tortuosity of the medium is considered.     

直接有限元法求解广义磁热弹二维旋转问题

为了验证直接有限元法求解广义磁热弹耦合旋转问题的有效性及准确性,该文基于Lord和Shulman(L-S)广义热弹性理论,采用直接有限元方法,求解了置于磁场中的旋转半无限大体受热冲击作用的动态响应问题.文中给出了L-S型广义磁热弹耦合旋转问题的控制方程,建立了L-S型广义磁热弹旋转问题的虚位移原理,推导得到了相应的有限...     

Annular extrudate swell of pseudoplastic and viscoplastic fluids

Numerical simulations have been undertaken for the benchmark problem of annular extrudate swell present in pipe extrusion and parison formation in blow molding. The finite element method (FEM) is used to provide numerical results for different inner/outer diameter ratios κ under steady-state conditions. The Herschel-Bulkley model of viscoplasticity is used with the Papanastasiou regularization, which reduces with appropriate parameter choices to the Bingham–Papanastasiou, power-law and Newtonian models. The present results provide the shape of the extrudate, and in particular the thickness and diameter swells, as a function of the dimensionless power-law index (in the case of pseudoplasticity) and the dimensionless yield stress (in the case of viscoplasticity). The pressures from the simulations have been used to compute the excess pressure losses in the system (exit correction). While shear-thinning leads to reduced swelling relative to the Newtonian values for all κ-values, the opposite is true for shear-thickening fluids, which exhibit considerable swelling. Viscoplasticity leads to decreased extrudate swell as the dimensionless yield stress goes from zero (Newtonian behaviour) to an asymptotic value of fully plastic behaviour. The exit correction decreases to zero with a decrease in the power-law index to zero and an increase in the dimensionless yield stress to its asymptotic limit. However, the decrease is not monotonic: for power-law fluids it has maxima in the range of power-law indices between 0.8 and 0.6, while for viscoplastic fluids it has maxima around Bingham number values of 5.     

一维六方准晶中纳米尺度开裂孔洞的III型断裂力学

研究了一维六方准晶中纳米尺度开裂孔洞的Ⅲ型断裂力学问题。基于复变弹性理论和表面弹性理论获得了考虑表面效应时椭圆孔边裂纹的应力场、应力强度因子和能量释放率的解析表达;讨论了缺陷尺寸、裂纹/孔洞比、耦合系数和施加载荷对应力强度因子和能量释放率的影响。研究表明：考虑表面效应且缺陷的尺寸在纳米尺度时,声子场和相位子场的无量纲应力强度因子以及无量纲能量释放率具有明显的尺寸依赖;裂纹相对尺寸较小时,表面效应对声子场和相位子场的无量纲应力强度因子影响较小;纳米尺度时无量纲能量释放率随耦合系数的增加而增大;耦合系数一定时,无量纲能量释放率受到椭圆孔尺寸影响;随着声子场载荷的增大,无量纲能量释放率先减小后增加,最后趋于稳定;无量纲能量释放率随相位子场载荷的增大单调减小,非常小和非常大的声子场载荷(或相位子场载荷)屏蔽了相位子场载荷(或声子场载荷)的影响。     

Resistance in the core flow of a fiber suspension in a tube

Using a more simple formulation than that of [1], the problem of core flow of a fiber suspension in a straight tube of circular cross section is considered. The solution is sought by the small parameter method. The first approximation thus found is compared in detail with experiment. Formulas are obtained for the tube resistance coefficient and the dimensionless thickness of the near-wall layer. The laminar flow in the first approximation is characterized by three dimensionless complexes. Both the resistance coefficient and the dimensionless thickness of the near-wall layer depend on only two of the complexes, and so partial simulation of the flow is possible.     

一维六方准晶中纳米尺度开裂孔洞的III型断裂力学

研究了一维六方准晶中纳米尺度开裂孔洞的Ⅲ型断裂力学问题。基于复变弹性理论和表面弹性理论获得了考虑表面效应时椭圆孔边裂纹的应力场、应力强度因子和能量释放率的解析表达;讨论了缺陷尺寸、裂纹/孔洞比、耦合系数和施加载荷对应力强度因子和能量释放率的影响。研究表明：考虑表面效应且缺陷的尺寸在纳米尺度时,声子场和相位子场的无量纲应力强度因子以及无量纲能量释放率具有明显的尺寸依赖;裂纹相对尺寸较小时,表面效应对声子场和相位子场的无量纲应力强度因子影响较小;纳米尺度时无量纲能量释放率随耦合系数的增加而增大;耦合系数一定时,无量纲能量释放率受到椭圆孔尺寸影响;随着声子场载荷的增大,无量纲能量释放率先减小后增加,最后趋于稳定;无量纲能量释放率随相位子场载荷的增大单调减小,非常小和非常大的声子场载荷(或相位子场载荷)屏蔽了相位子场载荷(或声子场载荷)的影响。     

A profile tracking system for investigating the behaviour of discharges in moving environments

A profile tracking system (PTS) is developed to make measurements at previously inaccessible dimensionless distances from a discharge in a moving environment. With this system the measuring components remain stationary as the source is towed away from them. Image processing techniques are then employed to track concentration profiles through a sequence of images, which are a constant distance from the flow source. Data obtained from the PTS are shown to be both reliable and accurate. The benefits and limitations of the PTS when compared with normally configured laser-induced fluorescence systems are discussed. PTS is employed to gather data at dimensionless distances that are an order of magnitude greater than those previously obtained for the co-flowing jet problem. Received: 30 June 1999/Accepted: 21 March 2001     

Properties and appropriate conditions of stress reduction factor and thermal shock resistance parameters for ceramics

Through introducing the analytical solution of the transient heat conduction problem of the plate with convection into the thermal stress field model of the elastic plate, the stress reduction factor is presented explicitly in its dimensionless form. A new stress reduction factor is introduced for the purpose of comparison. The properties and appropriate conditions of the stress reduction factor, the first and second thermal shock resistance (TSR) parameters for the high and low Biot numbers, respectively, and the approximation formulas for the intermediate Biot number-interval are discussed. To investigate the TSR of ceramics more accurately, it is recommended to combine the heat transfer theory with the theory of thermoelasticity or fracture mechanics or use a numerical method. The critical rupture temperature difference and the critical rup- ture dimensionless time can be used to characterize the TSR of ceramics intuitively and legibly.     

Casson Fluid Flow Due to Non-Coaxial Rotation of a Porous Disk and the Fluid at Infinity Through a Porous Medium

The flow of the Casson fluid due to non-coaxial rotation of a disk and the fluid at infinity is investigated. Partial differential equations are made dimensionless and coupled. The exact solution of the resultant nonlinear initial-boundary-value problem is solved by applying the Laplace transform. The shear stresses at the disk surface and the steady state stresses are computed. The effects of dimensionless parameters on the dimensionless primary and secondary velocities are analyzed.     

Functionally graded piezoelectric strip with a penny-shaped crack under electromechanical loadings

In this paper, the mixed-mode penny-shaped crack problem for a functionally graded piezoelectric material (FGPM) strip is considered. It is assumed that the electroelastic properties of the strip vary continuously along the thickness of the strip, and that the strip is under in-plane electromechanical loadings. The problem is formulated in terms of a system of singular integral equations. The stress and electric displacement intensity factors are presented for various values of dimensionless parameters representing the crack size, the crack location, and the material nonhomogeneity.     

Effect of Reactive Surface Areas Associated with Different Particle Shapes on Chemical-Dissolution Front Instability in Fluid-Saturated Porous Rocks

In this article, the effect of reactive surface areas associated with different particle shapes on the reactive infiltration instability in a fluid-saturated porous medium is investigated through analytically deriving the dimensionless pore-fluid pressure-gradient of a coupled system between porosity, pore-fluid flow and reactive chemical-species transport within two idealized porous media consisting of spherical and cubic grains respectively. Compared with the critical dimensionless pore-fluid pressure-gradient of the coupled system, the derived dimensionless pore-fluid pressure-gradient can be used to assess the instability of a chemical dissolution front within the fluid-saturated porous medium. The related theoretical analysis has demonstrated that (1) since the shape coefficient of spherical grains is greater than that of cubic grains, the chemical system consisting of spherical grains is more unstable than that consisting of cubic grains, and (2) the instability likelihood of a natural porous medium, which is comprised of irregular grains, is smaller than that of an idealized porous medium, which is comprised of regular spherical grains. To simulate the complicated morphological evolution of a chemical dissolution front in the case of the chemical dissolution system becoming supercritical, a numerical procedure is proposed for solving this kind of problem. The related numerical results have demonstrated that the reactive surface area associated with different particle shapes can have a significant influence on the morphological evolution of an unstable chemical-dissolution front within fluid-saturated porous rocks.     

Dynamic stability of a pipe conveying fluid with an uncertain computational model

This paper deals with the problem of a pipe conveying fluid of interest in several engineering applications, such as micro-systems or drill-string dynamics. The deterministic stability analysis developed by Paidoussis and Issid (1974) is extended to the case for which there are model uncertainties induced by modeling errors in the computational model. The aim of this work is twofold: (1) to propose a probabilistic model for the fluid–structure interaction considering modeling errors and (2) to analyze the stability and reliability of the stochastic system. The Euler–Bernoulli beam model is used to model the pipe and the plug flow model is used to take into account the internal flow in the pipe. The resulting differential equation is discretized by means of the finite element method and a reduced-order model is constructed from some eigenmodes of the beam. A probabilistic approach is used to model uncertainties in the fluid–structure interaction. The proposed strategy takes into account global uncertainties related to the noninertial coupled fluid forces (related to damping and stiffness). The resulting random eigenvalue problem is used to analyze flutter and divergence unstable modes of the system for different values of the dimensionless flow speed. The numerical results show the random response of the system for different levels of uncertainty, and the reliability of the system for different dimensionless speeds and levels of uncertainty.     

A finite crack in a semi-infinite strip of a graded piezoelectric material under electric loading

In this paper, the mixed-mode crack problem for a functionally graded piezoelectric material (FGPM) strip is considered. It is assumed that the electroelastic properties of the strip vary continuously along the thickness of the strip, and that the strip is under in-plane electric loading. The problem is formulated in terms of a system of singular integral equations. The stress and electric displacement intensity factors are presented for various values of dimensionless parameters representing the crack size, the crack location, and the material nonhomogeneity.