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.     

The buckling of elastic spherical caps

The John equations are used to model the buckling of a simply-supported elastic spherical cap that is subjected to a constant uniform external load . The Liapunov-Schmidt method is used to solve these equations. We show that solutions possessing circular, pear-shaped, elliptical, triangular, square-shaped, pentagonal and a variety of other symmetries branch from the unbuckled state of the shell. The stability of these solutions is discussed. Some numerics that complement the analytical results are also included.     

A novel simulation model for ring type ultrasonic motor

In this paper, a novel mathematical model for a traveling wave ultrasonic motor, developed by Alenia Spazio, now Alcatel Alenia Space Italia S.p.A. (Roma, Italy), within an Italian Space Agency (ASI) program, is described. The dynamic equations for the stator and the rotors of the ultrasonic motor are assembled into a differential system, whose equations are coupled by terms which represent interface generalized forces. Neglecting transient conditions, the complete mathematic model of the system is solved and an iterative process is developed, in order to obtain the motor’s running curves for different operation parameters, geometric dimensions and physical features of the system. The algorithm is implemented in Matlab environment and a graphical user interface is constructed for user-friendly managing. The model, also validated by means of experimental tests, can be used for parametric analyses with respect to different parameters, in order to optimize motor’s configuration. It represents a simple but powerful aid to determine final motor design that can satisfy specifications or to predict motor’s behavior under different working conditions.     

The determination of a general time creep compliance relation of linear viscoelastic materials under constant load and its extension to nonlinear viscoelastic behavior for the Burger model

Linear viscoelastic materials yield a creep function which only depends on time if creep experiments are performed under constant stress ₀. In practice, this condition is very difficult to realize, and as a consequence, the experiments are performed under constant force. For small strains the difference between the conditions of constant stress and constant force is negligible. Otherwise, the decrease in cross-section has to be taken into account and leads to increasing stress in the course of time for creep experiments under constant load. The Boltzmann superposition principle is solved under the condition of constant load and for strains . The creep complicance C(t; ₀) defined by the ratio becomes, in principle, dependent on the initial stress ₀. As a consequence, a set of creep compliance curves cannot be approximated with a simple parameter fit. Already the application of the solution on the Burger model yields a creep compliance curve with all three creep ranges. Furthermore, the mathematical structure of the time creep compliance relation of the Burger model allows nonlinear viscoelastic extension via the introduction of the yield strength _max and a nonlinearity parameter n _l . The creep behavior of PBT and PC can be described in the range of long times up to initial stresses ₀, being 75% for PBT and 60% for PC of the yield stress _max with only two or one free fit parameter, respectively.     

On a New Route to Chaos in Railway Dynamics

Cooperrider's mathematical model of a railway bogie running on a straight track has been thoroughly investigated due to its interesting nonlinear dynamics (see True [1] for a survey). In this article a detailed numerical investigation is made of the dynamics in a speed range, where many solutions exist, but only a couple of which are stable. One of them is a chaotic attractor.Cooperrider's bogie model is described in Section 2, and in Section 3 we explain the method of numerical investigation. In Section 4 the results are shown. The main result is that the chaotic attractor is created through a period-doubling cascade of the secondary period in an asymptotically stable quasiperiodic oscillation at decreasing speed. Several quasiperiodic windows were found in the chaotic motion.This route to chaos was first described by Franceschini [9], who discovered it in a seven-mode truncation of the plane incompressible Navier–Stokes equations. The problem investigated by Franceschini is a smooth dynamical system in contrast to the dynamics of the Cooperrider truck model. The forcing in the Cooperrider model includes a component, which has the form of a very stiff linear spring with a dead band simulating an elastic impact. The dynamics of the Cooperrider truck is therefore non-smooth.The quasiperiodic oscillation is created in a supercritical Neimark bifurcation at higher speeds from an asymmetric unstable periodic oscillation, which gains stability in the bifurcation. The bifurcating quasiperiodic solution is initially unstable, but it gains stability in a saddle-node bifurcation when the branch turns back toward lower speeds.The chaotic attractor disappears abruptly in what is conjectured to be a blue sky catastrophe, when the speed decreases further.     

Transient behaviour of the mechanical loss factor of high density polyethylene during creep and relaxation

The mechanical loss factor of high density polyethylene has been measured during creep and stress relaxation experiments. The application of a constant stress (creep) or strain (relaxation) resulted in an instantaneous increase in tan above the value obtained in absence of the creep or relaxation process. During the flow the tan-value decreased and approached an apparent equilibrium value slightly above the value obtained without the static load or deformation. This behaviour was observed both at 1 and 10 Hz and at 23 and 60 °C. It is suggested that this time dependence of the mechanical loss factor is associated with the basic mechanisms responsible for the creep or relaxation process itself.     

On the linear instability of the Hall-Stewartson vortex

It is demonstrated that the Hall-Stewartson leading-edge vortex is linearly unstable to viscous perturbations of the center-mode type. Center modes are found to occur in two reigons of Reynolds-number-wave-number space, in limits in which the axial wave number is large. The appropriate center-mode equations in these neighborhoods are established, and it emerges that the two sets are identical. The single system of equations, which depends on the azimuthal wave number m and a distance parameter only, is solved numerically for various values of m and . Highly unstable modes are found for large positive , and the results are shown to be in good agreement with proposed asymptotic expansions when >1. To lowest order, unstable modes have phase surfaces that rotate with the fluid: in addition constant phase surfaces propagate upstream but the group velocity is directed downstream. The growth rate of the instability increases faster than Reynolds number to the quarter power. This, together with the finding that the length scale of the unstable modes found goes to zero as the Reynolds number tends to infinity, makes this instability an unusual one.This work was supported by the Air Force Office of Scientific Research under contract AFOSR-89-0346 monitored by Dr. L. Sakell, and by the U.S. Army Research Office at the Mathematical Sciences Institute of Cornell University.     

Nonlinear dynamic buckling of a viscoelastic model under impact loading

A simple nonlinear buckling analysis is applied to a one-degree-of-freedom arch under impact loading in which viscous damping may also be included. Such a loading consists of a falling body striking centrally the joint mass of the arch in such a way that a completely plastic impact can be postulated. When there is no damping the exact dynamic buckling load for such a kind of loading-associated with an unbounded motion can be established by using a static criterion (approach). More specifically, it was shown that the dynamic buckling load corresponds to that unstable equilibrium state where the total potential energy of the system is zero. Furthermore, it was proved that the second variation of the total potential energy at the foregoing unstable equilibrium state is negative definite. This implies that the curve loading versus displacement resulting by the vanishing of the total potential energy has always a maximum on the afore mentioned unstable state. It was also found that the system may become sensitive to initial conditions. If damping is included the foregoing static criterion yields lower bound buckling estimates. These findings were verified by employing a highly efficient approximate technique as well as the numerical scheme of Runge-Kutta for solving any nonlinear initial-value problem.     

Mathematical Model for the Numerical Solution of Nonstationary Problems in Solid Mechanics by a Modified Godunov Method

A mathematical model of substance behavior under developed elastoplastic strains is worked out for solving onedimensional problems of solid mechanics. The model is based on the fundamental laws of conservation of mass, momentum, and total energy, Wilkins model, kinetic model of substance destruction, and modified Godunov method for the numerical solution of problems in mathematical physics. A hybrid difference scheme is constructed, which approximates acoustics equations with constant coefficients in smooth flows for the case of plane symmetry with the second order in time and space.     

Damage evolution of ductile material under impact loading

Nucleation, growth and coalescence of micro-voids result in the fracture of materials. Most mathematical models neglect nucleation and introduce initial damage, assuming it as a material constant. However, the original damage, which is formed during material working, is a material constant. The initial damage is a model parameter and depends on the load. Apparently, the predictability of such a model is poor.This paper made comparison and analysis of the four classical void growth models and showed their similarities. At the beginning of damage evolution, all the models follow a linear relationship in the form , where c is the size of micro voids and k is a parameter which relates the material and loading condition. With the concept of statistical micro-damage and the assumption of uniform void radius for new voids, a damage evolution equation was deduced based on the above void growth model. With this equation the effects of nucleation and growth at the beginning of the damage stage on the whole process of damage evolution can be calculated. The transition time from the nucleation dominant phase to the growth dominant phase can be determined. When the transition time is applied to the damage failure model of ductile material proposed by Johnson, the initial damage (f₀), a model parameter in the original model, can also be determined. The results of the derived damage evolution equation agree well with the previous research results.     

Effect of pulsations on the stability of a gas column

The effect of radial pulsations on the stability of a compressible cylindrical gas column surrounded by an ambient liquid is discussed. In the absence of pulsations, the stationary interface is susceptible to the Rayleigh capillary instability, which promotes the growth of longitudinal waves whose wave length is larger than 2 times the column radius, irrespective of the Reynolds number. A Floquet stability analysis for potential flow shows that the pulsations further destabilize the interface by extending the range of unstable wave numbers to a sequence of islands. A similar stability analysis for Stokes flow shows that the pulsations also have a destabilizing influence, though the presence of an insoluble surfactant has a competing stabilizing influence that may cause an overall reduction in the range of unstable wave numbers.     

Spontaneous rotation in the exact solution of magnetohydrodynamic equations for flow between two stationary impermeable disks

Magnetohydrodynamic (MHD) flow of a viscous electrically conducting incompressible fluid between two stationary impermeable disks is considered. A homogeneous electric current density vector normal to the surface is specified on the upper disk, and the lower disk is nonconducting. The exact von Karman solution of the complete system of MHD equations is studied in which the axial velocity and the magnetic field depend only on the axial coordinate. The problem contains two dimensionless parameters: the electric current density on the upper plate Y and the Batchelor number (magnetic Prandtl number). It is assumed that there is no external source that produces an axial magnetic field. The problem is solved for a Batchelor number of 0–2. Fluid flow is caused by the electric current. It is shown that for small values of Y, the fluid velocity vector has only axial and radial components. The velocity of motion increases with increasing Y, and at a critical value of Y, there is a bifurcation of the new steady flow regime with fluid rotation, while the flow without rotation becomes unstable. A feature of the obtained new exact solution is the absence of an axial magnetic field necessary for the occurrence of an azimuthal component of the ponderomotive force, as is the case in the MHD dynamo. A new mechanism for the bifurcation of rotation in MHD flow is found.     

Bifurcation points of the elastic-plastic equilibrium in a simple frame model

Summary This is a study of the characteristics of equilibrium bifurcation in a geometrically symmetric simple frame model in which the beam is loaded symmetrically and gradually. The deformability of the model is concentrated in two Shanley cells for whose elements a nonholonomic bilinear force-lengthening (or force-contraction) relation is adopted.The equilibrium bifurcation characteristics of this model are compared with those of Shanley's uniformly stressed column model. The characteristics of the two models coincide when instantaneous plasticisation of the deformable elements for the frame model are assumed, as occurs in the column model. But with a succession of different plastic states as the load increases, differences emerge. The most important being: in the frame model the borderline between the interval of unstable bifurcation points and the interval of stable bifurcation points is not necessarily, as it is in the column model, marked by a neutral (stationary load) bifurcation point.

---

Sommario S'indagano i caratteri della biforcazione dell'equilibrio di un modello di portale, simmetrico nella geometria e simmetricamente caricato sul traverso in modo graduale. La deformabilità del modello è concentrata in due celle alla Shanley, per i cui elementi si adotta una relazione bilineare anolonoma tra forze ed allungamenti (od accorciamenti).Si confrontano tali caratteri con quelli della biforcazione dell'equilibrio del ben noto modello di asta caricata di punta studiato da Shanley.Quando si ipotizzi, per il modello di portale, una plasticizzazione istantanea degli elementi deformabili (proprio come avviene nel modello di asta), si ha perfetta coincidenza dei caratteri in esame.All'opposto, quando si ipotizzi un succedersi di diversi stati plastici al crescere del carico, si constatano delle differenze; in particolare, per il modello di portale esiste la possibilità, a differenza di quanto avviene nel modello di asta, che non sia un punto di biforcazione neutra (a carico stazionario) a segnare il confine tra l'intervallo dei punti di biforcazione instabile e quello dei punti di biforcazione stabile.

---

---

Study supported by the Plasticity Group of the National Research Council.     

Formation of the shear-induced state in dilute cationic surfactant solutions

In cationic surfactant solutions a change of state occurs due to mechanical stresses. In the dilute regime of rodlike micelles the formation of a so-called Shear-Induced State (SIS) occurs above a critical shear rate. In this context dilute means that there is no sterical interaction between rodlike micelles, the solution is below the overlap concentration. Employing a mathematical model, it is shown that aggregation forces are weak compared to hydrodynamic forces. The mathematical formulation is based on a model of Israelachvili which describes the chemical potential of micelles. Hydrodynamic forces are calculated with a rigid-dumbbell model. SIS formation can be explained by the destruction of rodlike micelles.     

Memory effects in a non-uniform flow: A study of the behaviour of a tubular film of viscoelastic fluid

Summary Formal use of constitutive equations such as that ofOldroyd in the mathematical model of a flow leads, in general, to a higher order differential equation than is obtained for a purely viscous fluid, and so we expect to need more boundary conditions in order to specify the problem completely. (These extra boundary conditions may be thought of as arising from the need to specify what the fluid remembers of the flow outside the region of interest.) In flows which are uniform spatially, or uniform with time for a material element, the uniformity will provide the extra information and so no extra conditions are needed. Similarly for confined flows, where no new fluid enters the region of interest, no information about flow outside this region is needed.Here the steady flow of a tubular film of a viscoelastic fluid is studied with the particular aim of examining the effect of these extra boundary conditions in a situation where they may be expected to have some significant influence on the flow as a whole. The flow, while being geometrically complex, is essentially an elongational free-surface flow involving the biaxial stretching of a thin axisymmetric tubular film. Features of the constitutive equations studied are the presence of a non-zero relaxation time and the possibility of a variable viscosity. One effect of the non-zero relaxation time is that a tube of constant radius (possible but unstable for aNewtonian fluid) is not dynamically possible. Preliminary computational results suggest that the effect of the extra upstream boundary conditions is not large, and also have failed to show any major difference between the two generalisations of theMaxwell model which have been used.With 1 figure     

Some Results on Dynamics of Special Double-Loop ΣΔ-Modulators

The dynamic behaviour of a specific two-dimensional state space model with discontinuity is studied. This model arises from the study of double-loop -modulators with constant input. Using mathematical tools we explain certain simulation results, and some properties are derived. Simulations based on time-varying input are also provided.     

A Peierls criterion for the onset of deformation twinning at a crack tip

A criterion for the onset of deformation twinning (DT) is derived within the Peierls framework for dislocation emission from a crack tip due to Rice (J. Mech. Phys. Solids 40(2) (1992) 239). The critical stress intensity factor (SIF) is obtained for nucleation of a two-layer microtwin, which is taken to be a precursor to DT. The nucleation of the microtwin is controlled by the unstable twinning energyγ_ut, a new material parameter identified in the analysis. γ_ut plays the same role for DT as γ_us, the unstable stacking energy introduced by Rice, plays for dislocation emission. The competition between dislocation emission and DT at the crack tip is quantified by the twinning tendencyT defined as the ratio of the critical SIFs for dislocation nucleation and microtwin formation. DT is predicted when T>1 and dislocation emission when T<1. For the case where the external loading is proportional to a single load parameter, T is proportional to . The predictions of the criterion are compared with atomistic simulations for aluminum of Hai and Tadmor (Acta Mater. 51 (2003) 117) for a number of different crack configurations and loading modes. The criterion is found to be qualitatively exact for all cases, predicting the correct deformation mode and activated slip system. Quantitatively, the accuracy of the predicted nucleation loads varies from 5% to 56%. The sources of error are known and may be reduced by appropriate extensions to the model.     

Changes in the Seam-Filtration Mode with Stress Redistribution in Country Rock

A mathematical model for the mechanism of increasing oil output of productive seams is developed. The model involves a deliberate conversion of segments of the fault zone of the country rock to a supercritical state, which leads to a local redistribution of stresses in the block massif of rocks and an increase in contour and seam pressures. Based on solving the problem of restricted filtration, it is shown that the use of the proposed mechanism can ensure a relative increase in well production of 5–8%.     

Analytical Solution to the Riemann Problem of Three-Phase Flow in Porous Media

In this paper we study one-dimensional three-phase flow through porous media of immiscible, incompressible fluids. The model uses the common multiphase flow extension of Darcys equation, and does not include gravity and capillarity effects. Under these conditions, the mathematical problem reduces to a 2 × 2 system of conservation laws whose essential features are: (1) the system is strictly hyperbolic; (2) both characteristic fields are nongenuinely nonlinear, with single, connected inflection loci. These properties, which are natural extensions of the two-phase flow model, ensure that the solution is physically sensible. We present the complete analytical solution to the Riemann problem (constant initial and injected states) in detail, and describe the characteristic waves that may arise, concluding that only nine combinations of rarefactions, shocks and rarefaction-shocks are possible. We demonstrate that assuming the saturation paths of the solution are straight lines may result in inaccurate predictions for some realistic systems. Efficient algorithms for computing the exact solution are also given, making the analytical developments presented here readily applicable to interpretation of lab displacement experiments, and implementation of streamline simulators.     

Analysis of the Dynamic Stability of Electrical Graded Piezoelectric Circular Cylindrical Shells

I. INTRODUCTION In recent years, piezoelectric materials have been widely used as sensors and actuators in the areaof smart structures, and shell-type piezoelectric smart structures have become the focus of research.Various analytical and ?nite element models for piezoelectric shells have been developed[1??3]. Inaddition, a series of exact studies on the bending and buckling of piezoelectric circular cylindrical shellhave been presented by Chen et al.[4??7]. Recently, the use of functi…