Unsteady-state temperature distribution in a convecting fin of constant area

A solution for the unsteady-state temperature distribution in a fin of constant area dissipating heat only by convection to an environment of constant temperature, is obtained. The partial differential equation is separated into an ordinary differential equation with position as the independent variable, and a partial differential equation with position and time as the independent variables. The problem is solved for either a step function in temperature or a step function in heat flow rate, for zero time, at one boundary while the other boundary is insulated. The initial condition is taken as an arbitrary constant. The unspecified boundary values (temperature or heat flow rate) are presented for both cases by utilizing dimensionless plots. Experimental verification is presented for the case of constant heat flow rate boundary condition.     

Exact solutions for nonlinear transient flow model including a quadratic gradient term

The models of the nonlinear radial flow for the infinite and finite reservoirs including a quadratic gradient term were presented. The exact solution was given in real space for flow equation including quadratic gradiet term for both constant-rate and constant pressure production cases in an infinite system by using generalized Weber transform.Analytical solutions for flow equation including quadratic gradient term were also obtained by using the Hankel transform for a finite circular reservoir case. Both closed and constant pressure outer boundary conditions are considered. Moreover, both constant rate and constant pressure inner boundary conditions are considered. The difference between the nonlinear pressure solution and linear pressure solution is analyzed. The difference may be reached about 8% in the long time. The effect of the quadratic gradient term in the large time well test is considered.     

考虑二次梯度项及动边界的双重介质低渗透油藏流动分析

在传统试井模型的非线性偏微分方程中根据弱可压缩流体的假设,忽略了二次梯度项,对于低渗透油藏这种方法是有疑问的.低渗透问题一个显著的特点就是流体的流动边界随着时间不断向外扩展.为了更好地研究双重介质低渗透油藏中流体的流动问题,考虑了二次梯度项及活动边界的影响,同时考虑了低渗透油藏的非达西渗流特征,建立了双重介质低渗透油藏流动模型.采用Douglas-Jones预估-校正差分方法获得了无限大地层定产量生产时模型的数值解,分别讨论了不同参数变化时压力的变化规律及活动边界随时间的传播规律,还分析了考虑和忽略二次梯度项影响时模型数值解之间的差异随时间的变化规律,做出了典型压力曲线图版,这些结果可用于实际试井分析.     

Solution of the incompressible Navier–Stokes equations by the method of lines

The numerical method of lines (NUMOL) is a numerical technique used to solve efficiently partial differential equations. In this paper, the NUMOL is applied to the solution of the two‐dimensional unsteady Navier–Stokes equations for incompressible laminar flows in Cartesian coordinates. The Navier–Stokes equations are first discretized (in space) on a staggered grid as in the Marker and Cell scheme. The discretized Navier–Stokes equations form an index 2 system of differential algebraic equations, which are afterwards reduced to a system of ordinary differential equations (ODEs), using the discretized form of the continuity equation. The pressure field is computed solving a discrete pressure Poisson equation. Finally, the resulting ODEs are solved using the backward differentiation formulas. The proposed method is illustrated with Dirichlet boundary conditions through applications to the driven cavity flow and to the backward facing step flow. Copyright © 2015 John Wiley & Sons, Ltd.     

A Quasilinear Based Procedure for Saturated/Unsaturated Water Movement in Soils

The quasilinear form of Richards equation for one-dimensional unsaturated flow in soils can be readily solved for a wide variety of conditions. However, it cannot explain saturated/unsaturated flow and the constant diffusivity assumption, used to linearise the transient quasilinear equation, can introduce significant error. This paper presents a quasi-analytical solution to transient saturated/unsaturated flow based on the quasilinear equation, with saturated flow explained by a transformed Darcy's equation. The procedure presented is based on the modified finite analytic method. With this approach, the problem domain is divided into elements, with the element equations being solutions to a constant coefficient form of the governing partial differential equation. While the element equations are based on a constant diffusivity assumption, transient diffusivity behaviour is incorporated by time stepping. Profile heterogeneity can be incorporated into the procedure by allowing flow properties to vary from element to element. Two procedures are presented for the temporal solution; a Laplace transform procedure and a finite difference scheme. An advantage of the Laplace transform procedure is the ability to incorporate transient boundary condition behaviour directly into the analytical solutions. The scheme is shown to work well for two different flow problems, for three soil types. The technique presented can yield results of high accuracy if the spatial discretisation is sufficient, or alternatively can produce approximate solutions with low computational overheads by using large sized elements. Error was shown to be stable, linearly related to element size.     

On the comparison of the solutions obtained by using two different transform methods for the second problem of Stokes for Newtonian fluids

The unsteady flow over a plane wall which is initially at rest and the plate begins suddenly to oscillate in own plane is considered. The solution subject to the boundary and initial conditions is obtained by applying to the governing equation the Laplace transform method or Fourier transform method. A comparison of the solutions obtained by two transform methods for flow considered is given. It is shown that the solution obtained by the Laplace transform method or Fourier transform method is the sum of the steady-state and the transient parts. The transient parts are found in terms of definite integrands whose integrals are oscillatory functions. Therefore, the transient parts are expressed in terms of the tabulated functions.     

A computer model of pulse wave and input impedance in human arteries

To predict the propagation of pressure and flow pulses in arterial system and the variation of vascular input impedance, a branched and tapered tube model is studied through one-dimensional transient flow analysis. Coupling the continuity and momentum equations yields a group of quasilinear hyperbolic partial differential equations which can be solved numerically by using the method of characteristics. Several boundary conditions of the arterial system are also simplified suitably. The propagation of the pulses of the arterial system and the vascular input impedance is calculated on computer by using the dimensions and the physiological data of the arterial system. The results point out that the pressure and flow pulses of the arterial system and the vascular input impedance produced by this theoretical model is consistent quite well with the experimental results published.     

管内上随体Maxwell流体非定常流动

本文研究了上随体Maxwell流体在圆管内非定常流动规律,对于上随体Maxwell流体模型,导出了特殊的运动方程,分别应用隐式差分格式和Kantorovich变分法,求得数值解,对两类方法的结果进行比较,揭示了粘弹流效应对管内非定常流动规津的影响,根据上述研究认为,以上的特殊的变分方法适应于研究非定常流动。     

Rotating flow of a second-order fluid on a porous plate

A study is made of the unsteady flow engendered in a second-order incompressible, rotating fluid by an infinite porous plate exhibiting non-torsional oscillation of a given frequency. The porous character of the plate and the non-Newtonian effect of the fluid increase the order of the partial differential equation (it increases up to third order). The solution of the initial value problem is obtained by the method of Laplace transform. The effect of material parameters on the flow is given explicitly and several limiting cases are deduced. It is found that a non-Newtonian effect is present in the velocity field for both the unsteady and steady-state cases. Once again for a second-order fluid, it is also found that except for the resonant case the asymptotic steady solution exists for blowing. Furthermore, the structure of the associated boundary layers is determined.     

A VARIATIONAL INEQUALITY APPROACH FOR NON-STEADY SEEPAGE FLOW THROUGH THREE-DIMENSIONAL FRACTURE NETWORK

为了求解裂隙岩体有自由面非稳定渗流问题,将Darcy定律延拓至整个研究区域,使得潜在溢出边界条件满足Signorini型边界条件,建立了三维裂隙网络非稳定渗流问题的抛物型变分不等式（parabolic variational inequality,PVI）提法,并证明其与偏微分方程（partial differential equation,PDE）提法的等价性,从而将自由面上的流量条件以及潜在溢出边界上的互补条件转化成自然边界条件,降低该问题求解难度。同时给出了基于PVI提法的有限元数值求解方法,通过与交叉裂隙模型理论解的对比分析,证明了该方法的正确性。最后将该方法对含复杂三维裂隙网络的边坡进行非稳定渗流分析,计算结果表明该方法对于复杂裂隙网络求解具有较强的可靠性和适应性。     

Algorithm for analysis of flows in ribbed annuli

Spectral methods for analyses of steady flows in annuli bounded by walls with either axi‐symmetric or longitudinal ribs are developed. The physical boundary conditions are enforced using the immersed boundary conditions concept. In the former case, the Stokes stream function is used to eliminate pressure and to reduce system of field equations to a single fourth‐order partial differential equation. The ribs are assumed to be periodic in the axial direction and this permits representation of the solution in terms of the Fourier expansion. In the latter case, the problem is reduced to the Laplace equations for the flow modifications that can be expressed in terms of the Fourier expansions. The modal functions, which are functions of the radial coordinate, are discretized using Chebyshev polynomials. The problem formulations are closed using either the fixed volume flow rate constraint or the fixed pressure gradient constraint. Various tests have been carried out in order to demonstrate the spectral accuracy of the discretizations, as well as the spectral accuracy of the enforcement of the flow boundary conditions at the ribbed walls using the immersed boundary conditions concept. Special linear solver that takes advantage of the matrix structure has been implemented in order to reduce computational time and memory requirements. It is shown that the algorithm has superior performance when one is interested in the analysis of a large number of geometries, as part of the coefficient matrix that corresponds to the field equation is always the same and one needs to change only the part of the matrix that corresponds to the boundary relations when changing geometry of the flow domain. Copyright © 2011 John Wiley & Sons, Ltd.     

Uniqueness of the self-similar solutions of three-dimensional laminar boundary-layer equations

The equations of the three-dimensional laminar boundary layer on lines of flow outflow and inflow are studied for conical outer flow under the assumption that the Prandtl number and the productρμ are constant. It is shown that in the case of a positive velocity gradient of the secondary flow (α₁>0) the additional conditions which result from the physical flow pattern determine a unique solution of the system of boundary-layer equations. For a negative velocity gradient of the secondary flow (α₁≤0) these conditions are satisfied by two solutions. An approximate solution is obtained for the boundary layer equations which is in rather good agreement with the numerical integration results. Compressible gas flow in a three-dimensional laminar boundary layer is described by a system of nonlinear differential equations whose solution is not unique for given boundary conditions. Therefore additional conditions resulting from the physical pattern of the gas flow are imposed on the resulting solution. In the solution of problems with a negative pressure gradient these additional conditions are sufficient for a unique selection of the solution of the boundary-layer equations. However, in the case of a positive pressure gradient the solution of the boundary-layer equations satisfying the boundary and additional conditions may not be unique. In particular, in [1] in a study of a three-dimensional laminar boundary layer in the vicinity of the stagnation point it was shown that for $$c = {{\frac{{\partial v_e }}{{\partial y}}} \mathord{\left/ {\vphantom {{\frac{{\partial v_e }}{{\partial y}}} {\frac{{\partial u_e }}{{\partial x}}}}} \right. \kern-\nulldelimiterspace} {\frac{{\partial u_e }}{{\partial x}}}} > 0$$ the solution is unique, while for c<0 there are two solutions. In the present paper we study the question of the uniqueness of the self-similar solution of the three-dimensional laminar boundary-layer equations on lines of flow outflow and inflow for a conical outer flow.     

An exact analytical solution of the unsteady magnetohydrodynamics nonlinear dynamics of laminar boundary layer due to an impulsively linear stretching sheet

In this paper, we investigate the theoretical analysis for the unsteady magnetohydrodynamic laminar boundary layer flow due to impulsively stretching sheet. The third-order highly nonlinear partial differential equation modeling the unsteady boundary layer flow brought on by an impulsively stretching flat sheet was solved by applying Adomian decomposition method and Pade approximants. The exact analytical solution so obtained is in terms of rapidly converging power series and each of the variants are easily computable. Variations in parameters such as mass transfer (suction/injection) and Chandrasekhar number on the velocity are observed by plotting the graphs. This particular problem is technically sound and has got applications in expulsion process and related process in fluid dynamics problems.     

三维势流场的比例边界有限元求解方法

比例边界有限元法(SBFEM)是线性偏微分方程的一种新的数值求解方法。该方法只对计算域边界利用Galerkin方法进行数值离散,相对于有限元方法(FEM)减少了一个空间坐标的维数,而在减少的空间坐标方向利用解析方法进行求解;相对于边界元法(BEM),比例边界有限元方法不需要基本解,避免了奇异积分的计算,所以它结合了有限元和边界元方法的优点。本文建立了利用比例边界有限元法求解三维Laplace方程的数值模型并用于计算三维物体周围的水流场,将计算结果与解析解和边界元方法进行了对比,结果表明此方法可以很好地模拟水流场,且具有较高的计算精度。     

Analytic solution for axisymmetric flow over a nonlinearly stretching sheet

An analysis is performed for the boundary-layer flow of a viscous fluid over a nonlinear axisymmetric stretching sheet. By introducing new nonlinear similarity transformations, the partial differential equations governing the flow are reduced to an ordinary differential equation. The resulting ordinary differential equation is solved using the homotopy analysis method (HAM). Analytic solution is given in the form of an infinite series. Convergence of the obtained series solution is explicitly established. The solution for an axisymmetric linear stretching sheet is obtained as a special case.     

MHD boundary-layer flow of an upper-convected Maxwell fluid in a porous channel

Two-dimensional magnetohydrodynamic (MHD) boundary layer flow of an upper-convected Maxwell fluid is investigated in a channel. The walls of the channel are taken as porous. Using the similarity transformations and boundary layer approximations, the nonlinear partial differential equations are reduced to an ordinary differential equation. The developed nonlinear equation is solved analytically using the homotopy analysis method. An expression for the analytic solution is derived in the form of a series. The convergence of the obtained series is shown. The effects of the Reynolds number Re, Deborah number De and Hartman number M are shown through graphs and discussed for both the suction and injection cases.     

薄板弯曲问题的神经网络方法

为发展神经网络方法在求解薄板弯曲问题中的应用,基于Kirchhoff板理论,提出一种采用全连接层求解薄板弯曲四阶偏微分控制方程的神经网络方法。首先在求解域、边界中随机生成数据点作为神经网络输入层的参数,由前向传播系统求出预测解;其次计算预测解在域内及边界处的误差,利用反向传播系统优化神经网络系统的计算参数;最后,不断训练神经网络使误差收敛,从而得到薄板弯曲的挠度精确解。以不同边界、荷载条件的三角形、椭圆形、矩形薄板为例,利用所提方法求解其偏微分方程,与理论解或有限元解对比,讨论了影响神经网络方法收敛的因素。研究表明,本文方法对求解薄板弯曲问题的四阶偏微分控制方程具有一定的适应性,其收敛性受多种条件影响。相比有限元,该方法收敛速度较慢,在复杂的边界条件下收敛性不佳,然而其不基于变分原理,无需计算刚度矩阵,便可获得较高的计算精度。     

Effect of quadratic pressure gradient term on a one-dimensional moving boundary problem based on modified Darcy’s law

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 equations with moving boundary conditions into a closed system with fixed boundary conditions; then the solution can be stably 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 conclusion that such moving boundary problems must incorporate the quadratic pressure gradient term in their governing equations; the sensitive effects of the quadratic pressure gradient term tend to diminish, with the dimensionless threshold pressure gradient increasing for the one-dimensional problem.     

一种改进格林元方法及在渗流问题中的应用

采用格林公式和基本解推导出直接边界积分方程来求解渗流问题.边界积分方程数值离散基于格林元方法(Green element methond),改进了原方法中压力和压力导数的求解方法,命名为混合边界元方法(Mixed boundary element method).相较于格林元类方法,该方法显式考虑了求解节点的外法向流量值和压力值,并使求得的数值解在求解区域上能够连续,符合实际的物理过程,在不增加额外未知数的情况下提高了计算精度.分析了不同网格类型对模拟计算结果的影响,并对稳定渗流问题、非稳定(瞬态)渗流问题和非稳态问题进行了实例计算,结果显示改进方法提高了计算精度,并对各类渗流问题有较好的适应性.     

Finite element solution of multi-dimensional two-phase flow through porous media with arbitrary heating conditions

Mechanistic models for flow regime transitions and drag forces proposed in an earlier work are employed to predict two-phase flow characteristics in multi-dimensional porous layers. The numerical scheme calls for elimination of velocities in favor of pressure and void fraction. The momentum equations for vapor and liquid then can be reduced to a system of two partial differential equations (PDEs) which must be solved simultaneously for pressure and void fraction.

Solutions are obtained both in two-dimensional cartesian and in axi-symmetric coordinate systems. The porous layers in both cases are composed of regions with different permeabilities. The finite element method is employed by casting the PDEs in their equivalent variational forms. Two classes of boundary conditions (specified pressure and specified fluid fluxes) can be incorporated in the solution. Volumetric heating can be included as a source term. The numerical procedure is thus suitable for a wide variety of geometry and heating conditions. Numerical solutions are also compared with available experimental data.