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.     

A moving finite element method for the population balance equation

A moving finite element algorithm has been compared against the upwind-differencing and Smolarkiewicz methods for the population balance equation of multicomponent particle growth processes. Analytical solutions and an error function have been used to test the numerical methods. The moving finite elements technique is much more accurate than other methods for a wide range of parameters. Since this method uses moving grids, it is able to model very narrow particle size distributions. It is also shown that the method can be extended to solve condensational growth problems which include particle curvature and non-continuum mass transfer effects.     

极坐标平面问题的四阶差分及应用

推导出了极坐标系下双调和方程的差分公式,用逐次超松弛迭代法求出圆板平面应力问题的差分解,并和解析解作比较,验证了差分公式的正确性,为解决圆域及其类似区域的平面问题提供了新方法．     

Global solution of the inverse problem for a class of nonlinear evolution equations of dispersive type

The inverse problem for a class of nonlinear evolution equations of dispersive type was reduced to Cauchy problem of nonlinear evolution equation in an abstract space. By means of the semigroup method and equipping equivalent norm technique, the existence and uniqueness theorem of global solution was obtained for this class of abstract evolution equations, and was applied to the inverse problem discussed here. The existence and uniqueness theorem of global solution was given for this class of nonlinear evolution equations of dispersive type. The results extend and generalize essentially the related results of the existence and uniqueness of local solution presented by YUAN Zhong-xin. Contributed by Chen Yu-shu Foundation item: the National Natural Science Foundation of China (Significance 199990510); the National Key Basic Research Special Foundation of China (G1998020316); Liuhui Center for Applied Mathematics, Nankai University & Tianjin University Biography: Chen Fang-qi (1963-)     

Transversal vibrations of the axially moving sandwich belts

Inspired by literature on free transversal vibrations of one axially moving belt, we derive and solve analytically coupled partial differential equations of the transversal vibrations of an axially moving sandwich double-belt system. A numerical experiment and visualization are carried out.     

Numerical perturbation method for approximate solution of poisson's equation on a moderately deforming grid

In problems such as the computation of incompressible flows with moving boundaries, it may be necessary to solve Poisson's equation on a large sequence of related grids. In this paper the LU decomposition of the matrix A ₀ representing Poisson's equation discretized on one grid is used to efficiently obtain an approximate solution on a perturbation of that grid. Instead of doing an LU decomposition of the new matrix A , the RHS is perturbed by a Taylor expansion of A ^?1 about A ₀. Each term in the resulting series requires one ‘backsolve’ using the original LU . Tests using Laplace's equation on a square/rectangle deformation look promising; three and seven correction terms for deformations of 20% and 40% respectively yielded better than 1% accuracy. As another test, Poisson's equation was solved in an ellipse (fully developed flow in a duct) of aspect ratio 2/3 by perturbing about a circle; one correction term yielded better than 1% accuracy. Envisioned applications other than the computation of unsteady incompressible flow include: three-dimensional parabolic problems in tubes of varying cross-section, use of ‘elimination’ techniques other than LU decomposition, and the solution of PDEs other than Poisson's equation.     

A general solution of the axisymmetric contact problem for biphasic cartilage layers

A unilateral axisymmetric contact problem for articular cartilage layers is considered. The articular cartilages bonded to subchondral bones are modeled as biphasic materials consisting of a solid phase and a fluid phase. It is assumed that the subchondral bones are rigid and shaped like bodies of revolution with arbitrary convex profiles. The obtained closed-form analytical solution is valid over time periods compared with the typical diffusion time and can be used for increasing loading.     

Forced nonlinear oscillations of semi-infinite cables and beams resting on a unilateral elastic substrate

In this work, we study the nonlinear oscillations of mechanical systems resting on a (unilateral) elastic substrate reacting in compression only. We consider both semi-infinite cables and semi-infinite beams, subject to a constant distributed load and to a harmonic displacement applied to the finite boundary. Due to the nonlinearity of the substrate, the problem falls in the realm of free-boundary problems, because the position of the points where the system detaches from the substrate, called Touch Down Points (TDP), is not known in advance. By an appropriate change of variables, the problem is transformed into a fixed-boundary problem, which is successively approached by a perturbative expansion method. In order to detect the main mechanical phenomenon, terms up to the second order have to be considered. Two different regimes have been identified in the behaviour of the system, one below (called subcritical) and one above (called supercritical) a certain critical excitation frequency. In the latter, energy is lost by radiation at infinity, while in the former this phenomenon does not occur and various resonances are observed instead; their number depends on the statical configuration around which the system performs nonlinear oscillations.     

一类非线性周期系统响应的精细积分法

对于一类非线性周期／变系数微分方程,提出基于精细积分法的数值解法,处理非线性周期／变系数微分方程系统的响应问题,其积分策略是：采用精细积分格式处理常系数部分;采用线性插值格式处理非线性周期／变系数部分,既继承精细积分方程高度准确的特点,又保证足够的精度与较小的计算量。通过数值算例,与以往与用的微分方程直接数值积分法（如预估－校正哈明法）求得的解加以比较表明,对于给定的精度要求,精细积分法更经济有效     

The study of quadruple integral equations involving inverse Mellin transforms and its application to a contact problem for a wedge shaped elastic solid in antiplane shear distribution

Closed form solution of quadruple integral equations involving inverse Mellin transforms has been obtained. The solution of quadruple integral equations is used in solving a two dimensional four-part mixed boundary value contact problem for an elastic wedge-shaped region as an application. Closed form expression for shear stress has been obtained. Finally, numerical results for shear stress are obtained and shown graphically.     

Exact solution of external tangential contact problem for a transversely isotropic elastic half-space

Summary  The following mixed boundary-value problem for a transversely isotropic elastic half-space is considered. Arbitrary tangential displacements are prescribed at the exterior of a circle, while the interior of the circle is free of tangential stress, and the normal stress vanishes all over the boundary. The governing integral equation is solved exactly, in closed form, and in terms of elementary functions. The method of continuation of solutions previously published by the author has been used here. Several examples are considered. No similar results has been reported before, even in the case of an isotropic body. Received 8 May 2000; accepted for publication 20 July 2000     

A self- similar solution of a shock propagation in a mixture of a non-ideal gas and small solid particles

Similarity solutions are obtained for unsteady, one-dimensional self-similar flow behind a strong shock wave, driven by a moving piston, in a dusty gas. The dusty gas is assumed to consist of a mixture of small solid particles and a non-ideal gas, in which solid particles are continuously distributed. It is assumed that the equilibrium flow-condition is maintained and variable energy input is continuously supplied by the piston. Solutions are obtained under both the isothermal and adiabatic conditions of the flow-field. The spherical case is worked out in detail to investigate to what extent the flow-field behind the shock is influenced by the non-idealness of the gas in the mixture as well as by the mass concentration of the solid particles, by the ratio of density of the solid particles to the initial density of the mixture and by the energy input due to moving piston. A comparison is also made between isothermal and adiabatic cases.     

Green's function of the displacement boundary value problem for a heat source in an infinite plane with an arbitrary shaped rigid inclusion

Summary Green's functions of the displacement boundary value problem are derived within two-dimensional thermoelasticity for a heat source in an infinite plane with an arbitrary shaped rigid inclusion. The following two cases are considered: either rigid-body displacements and rigid-body rotations of the inclusion are allowed or no rigid-body displacements and no rigid-body rotations of the inclusion are possible. To solve these problems, fundamental solutions are developed for a point heat source, for rigid-body rotations of the inclusion, and for concentrated loads acting on the inclusion. Complex stress functions, temperature function, a rational mapping function and the thermal dislocation method are used for the analysis. In analytical examples, distributions of stresses are developed for an infinite plane with a rectangular rigid inclusion. Received 5 August 1998; accepted for publication 1 December 1998     

Impact of a long thin body on a cylindrical cavity in liquid: A plane problem

An approach is developed to the investigation of the shock interaction between a long thin cylindrical body and a cylindrical cavity in an infinite compressible perfect liquid. This process accompanies the supercavitation of the body. Three typical cases of cross-sectional dimensions of the body and the cavity are examined. For each case, a mixed nonstationary boundary-value problem with an unknown moving boundary is formulated. The unknown quantities are expanded into Fourier series. An auxiliary problem is solved using the Laplace transform to establish the relationship between the pressure and the velocity on the cavity surface. As a result, the problem is reduced to an infinite system of Volterra equations of the second kind solved simultaneously with the equation of transverse motion and the equation of the contact boundary. An asymptotic solution valid at the initial stage of interaction is obtained for all the three cases, and a numerical solution is found for the most typical case __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 6, pp. 32–53, June 2006.     

Comment on “Analysing flow and heat transfer of a viscoelastic fluid over a semi-infinite horizontal moving flat plate, IJNLM, 43 (2008) 772”

In this report we make corrections to a paper by Rafael Cortell (IJNLM, 43 (2008) 772). In this paper (IJNLM, 43 (2008) 772) the author considered an incorrect equation (by including a wrong term instead of correct term) for his considered flow and heat transfer problem and reported unrealistic results as a consequence. In this communication we have considered the correct equation (replacing the incorrect term by correct one) and solved it analytically. Our analysis shows that the unrealistic behavior as observed by Cortell (IJNLM, 43 (2008) 772) persists no more if one considers the correct term.     

The 3-D stress intensity factor for a half plane crack in a transversely isotropic solid due to the motion of loads on the crack faces

Three-dimensional analysis of a half plane crack in a transversely isotropic solid is performed. The crack is subjected to a pair of normal point loads moving in a direction perpendicular to the crack edge on its faces. Transform methods are used to reduce the boundary value problem to a single integral equation that can be solved by the Wiener-Hopf technique. The Cagniard-de Hoop method is employed to invert the transforms. An exact expression is derived for the mode I stress intensity factor as a function of time and position along the crack edge. Some features of the solution are discussed through numerical results. The project supported by the Guangdong Provincial Natural Science Foundation and the Science Foundation of Shantou University     

A new formulation for solving 3-D time dependent rolling contact problems of a rigid body on a viscoelastic half-space

The paper deals with a new formulation for solving the rolling contact problem without friction of a rigid body on a viscoelastic half-space in three dimensions. Assuming that the material behavior is independent of time for a sufficiently short time duration, the viscoelastic contact problem is transformed into elastic like problems. Then the contact problem is solved using a direct numerical method at each time step. The convergence of the method in time and space is good for a spherical indenter. The dissymmetry of the contact patch due to hysteresis was found in three dimensions for the spherical indenter and two cylinders of different width. Finally the method was tested for a sinusoidal varying speed and shows a good efficiency.     

A note on internal damping induced self-excited vibration in a rotor by considering source loading of a DC motor drive

It is well known that rotors become unstable beyond a certain threshold spinning speed due to the non-conservative circulatory forces, which arise out of rotating internal damping. In this note, it is shown that if the source loading of the non-ideal drive is considered then this instability manifests itself as a constant rotor spinning speed and a constant amplitude whirl orbit about an unstable equilibrium. A DC motor drive is considered and the corresponding steady-state spinning frequency and whirl orbit amplitude are analytically derived as functions of the drive and the rotor system parameters.     

A numerical investigation of the effects of the spanwise length on the 3-D wake of a circular cylinder

The numerical prediction of vortex-induced vibrations has been the focus of numerous investigations to date using tools such as computational fluid dynamics. In particular, the flow around a circular cylinder has raised much attention as it is present in critical engineering problems such as marine cables or risers. Limitations due to the computational cost imposed by the solution of a large number of equations have resulted in the study of mostly 2-D flows with only a few exceptions. The discrepancies found between experimental data and 2-D numerical simulations suggested that 3-D instabilities occurred in the wake of the cylinder that affect substantially the characteristics of the flow. The few 3-D numerical solutions available in the literature confirmed such a hypothesis. In the present investigation the effect of the spanwise extension of the solution domain on the 3-D wake of a circular cylinder is investigated for various Reynolds numbers between 40 and 1000. By assessing the minimum spanwise extension required to predict accurately the flow around a circular cylinder, the infinitely long cylinder is reduced to a finite length cylinder, thus making numerical solution an effective way of investigating flows around circular cylinders. Results are presented for three different spanwise extensions, namely πD/2, πD and 2πD. The analysis of the force coefficients obtained for the various Reynolds numbers together with a visualization of the three-dimensionalities in the wake of the cylinder allowed for a comparison between the effects of the three spanwise extensions. Furthermore, by showing the different modes of vortex shedding present in the wake and by analysing the streamwise components of the vorticity, it was possible to estimate the spanwise wavelengths at the various Reynolds numbers and to demonstrate that a finite spanwise extension is sufficient to accurately predict the flow past an infinitely long circular cylinder.