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.     

Asymptotic Solution for a Kind of Boundary Layer Problem

This paper studies the partial differential equation with a small coefficient in the highest-order item. This kind of equation is also named as boundary layer problem. The Burgers equation and modified Burgers equation are analyzed in this approach. First, these equations are transferred into the strong nonlinear ones, and then the corresponding strong nonlinear equations are solved based on the perturbation method. The results from the asymptotic method are comparable with those obtained from numerical computation. An erratum to this article is available at .     

Nonlinear Disturbances and Weak Discontinuities in a Supersonic Boundary Layer

The processes of wave disturbance propagation in a supersonic boundary layer with self-induced pressure [1–4] are analyzed. The application of a new mathematical apparatus, namely, the theory of characteristics for systems of differential equations with operator coefficients [5–8], makes it possible to obtain generalized characteristics of the discrete and continuous spectra of the governing system of equations. It is shown that the discontinuities in the derivatives of the solution of the boundary layer equations are concentrated on the generalized characteristics. It is established that in the process of flow evolution the amplitude of the weak discontinuity in the derivatives may increase without bound, which indicates the possibility of breaking of nonlinear waves traveling in the boundary layer.     

New Solution of Axisymmetric Contact Problem of Elasticity

We consider two problems of elasticity, namely, the Boussinesq problem about the action of a lumped force on a half-space and the related problem about the interaction of the half-space with a cylindrical rigid punch with plane base. In the classical statement, these problems have singular solutions. In the Boussinesq problem, the displacement under the action of the force is infinitely large, and in the punch problem, the infinitely large variable is the pressure on the punch boundary. In the present paper, these problems are solved with the use of relations of generalized elasticity derived regarding a medium element of small but finite dimensions rather than a traditional infinitesimal element. The structure parameter of the medium contained in the solutions can be determined experimentally. The obtained generalized solutions of the problems under study are regular.     

Solution of the Inverse Boundary Value Problem of Aerohydrodynamics for a Two-Element Airfoil

The basic inverse boundary value problem of aerohydrodynamics for a two-element airfoil is analytically and numerically solved in the complete formulation. The problems of designing biplane airfoils and an airfoil with a trailing-edge flap or leading edge flap (slat) are solved for a given distribution over the unknown contours of the velocity or the pressure as a function of the contour arc abscissa of the airfoils which depends on a finite number of the parameters.     

An Elastic-viscoplastic Quasistatic Contact Problem: Existence and Uniqueness of a Weak Solution

In this paper, the contact between an elastic-viscoplastic body and a deformable obstacle is studied. The effect of the damage, due to internal tension or compression and caused by the opening and growth of micro-cracks and micro-cavities, is also considered. The variational formulation leads to a coupled system of evolutionary equations. An existence and uniqueness result is established by using approximate problems, the pseudomonotone operators theory, Schauder fixed-point theorem and a comparison result.     

Delta Shock Wave Solution of the Riemann Problem for the Non-homogeneous Modified Chaplygin Gasdynamics

The motivation of the present study is to derive the solution of the Riemann problem for modified Chaplygin gas equations in the presence of constant external force. The analysis leads to the fact that in some special circumstances delta shock appears in the solution of the Riemann problem. Also, the Rankine–Hugoniot relations for delta shock wave which are utilized to determine the strength, position and propagation speed of the delta shocks have been derived. Delta shock wave solution to the Riemann problem for the modified Chaplygin gas equation is obtained. It is found that the external force term, appearing in the governing equations, influences the Riemann solution for the modified Chaplygin gas equation.

     

Solution of a Riemann problem for elasticity

This paper describes a numerical algorithm for the Riemann solution for nonlinear elasticity. We assume that the material is hyperelastic, which means that the stress-strain relations are given by the specific internal energy. Our results become more explicit under further assumptions: that the material is isotropic and that the Riemann problem is uniaxial. We assume that any umbilical points lie outside the region of physical relevance. Our main conclusion is that the Riemann solution can be obtained by the iterative solution of functional equations (Godunov iterations) each defined in one- or two-dimensional spaces.Supported in part by AFOSR-88-0025.     

Local and Global Injective Solution of the Rotationally Symmetric Sphere Problem

There are problems in the classical linear theory of elasticity whose closed form solutions, while satisfying the governing equations of equilibrium together with well-posed boundary conditions, predict the existence of regions, often quite small, inside the body where material overlaps. Of course, material overlapping is not physically realistic, and one possible way to prevent it combines linear theory with the requirement that the deformation field be injective. A formulation of minimization problems in classical linear elasticity proposed by Fosdick and Royer [3] imposes this requirement through a Lagrange multiplier technique. An existence theorem for minimizers of plane problems is also presented. In general, however, it is not certain that such minimizers exist. Here, the Euler–Lagrange equations corresponding to a family of three-dimensional problems is investigated. In classical linear elasticity, these problems do not have bounded solutions inside a body of anisotropic material for a range of material parameters. For another range of parameters, bounded solutions do exist but yield stresses that are infinite at a point inside the body. In addition, these solutions are not injective in a region surrounding this point, yielding unrealistic behavior such as overlapping of material. Applying the formulation of Fosdick and Royer on this family of problems, it is shown that both the displacements and the constitutive part of the stresses are bounded for all values of the material parameters and that the injectivity constraint is preserved. In addition, a penalty functional formulation of the constrained elastic problems is proposed, which allows to devise a numerical approach to compute the solutions of these problems. The approach consists of finding the displacement field that minimizes an augmented potential energy functional. This augmented functional is composed of the potential energy of linear elasticity theory and of a penalty functional divided by a penalty parameter. A sequence of solutions is then constructed, parameterized by the penalty parameter, that converges to a function that satisfies the first variation conditions for a minimizer of the constrained minimization problem when this parameter tends to infinity. This approach has the advantages of being mathematically appealling and computationally simple to implement.     

二维弹性接触问题中随机边界元法与可靠性分析

发展了二维弹性接触问题中的随机边界元法，推导并建立了相应的随机边界元基本方程，并将所发展的方法用于静强度的可靠性分析，讨论了其数值解技术。通过算例分析表明，本文发展的方法是可行的。     

摩擦约束塑性力学变分不等原理的半反推法

带摩擦约束的弹塑性接触问题,由于摩擦约束条件是一种判别性的条件,它的变分问题的逆问题的研究比较困难。本文对弹塑性接触力学中的变分不等问题的逆问题进行了研究,改进了半反推法并将其应用到弹塑性变分不等原理的研究中,导出了摩擦约束弹塑性增量广义变分不等原理中的能量泛函,消除了用拉氏乘子法可能产生的临界变分现象,在证明中,巧妙地处理了增量表示的接触摩擦边界条件,避免了使用非线性泛函分析和凸分析,简化了证明。     

Approximate Solution of an Axisymmetric Contact Problem with Allowance for Tangential Displacements on the Contact Surface

A structurally nonlinear contact problem of a punch shaped like a paraboloid of revolution is studied. An equation for the contactpressure density is derived with allowance for the radial tangential displacements of the boundary points of an elastic halfspace. A method for constructing a closedform approximate solution is proposed. The effect of the tangential displacements on the main contact parameters is discussed.     

Problem of Breakdown of Solitary Waves and Discontinuities

The evolution of initial data of the solitary-wave type with time is investigated numerically. The solitary wave amplitude decreases due to the generation of short-wave radiation. This solution is interpreted as the solution with a discontinuity qualitatively analogous to the solution of the problem of the breakdown of an arbitrary discontinuity in dissipationless systems. The solitary wave amplitude reduction rate is estimated, first for a generalized Korteweg-de Vries equation and then for plasma waves. Features of the investigation are analyzed for cold and hot-electron plasmas.     

Analysis of Constitutive Assumptions for the Strain Energy of a Generalized Elastic Membrane in a Nonlinear Contact Problem

For very thin shell-like structures it is common to ignore bending effects and model the structure using simple membrane theory. However, since the thickness of the membrane is not modeled explicitly in simple membrane theory it is not possible to use the three-dimensional strain energy function directly. Approximations must be introduced like the assumptions of: no thickness changes, generalized plane stress or incompressibility. In contrast, the theory of a Cosserat generalized membrane uses the three-dimensional strain energy function directly, it includes both thickness changes and shear deformation and it allows contact conditions to be formulated on the interface of the membrane with another body instead of on the middle surface of the membrane. A specific nonlinear contact problem is used to study these effects and comparison is made with solutions of a hierarchy of theories which include different levels of deformation through the thickness of the membrane and different formulations of the contact conditions. The results indicate that within the context of a simple membrane the assumption of generalized plane stress is best for this problem and that a generalized contact condition extends the range validity of the simple membrane solution to thicker membranes.     

Asymptotic Solution of a Harmonic Contact Problem for a Permeable Stamp on a Liquid Saturated Base

A plane harmonic problem of vertical vibrations of a rigid permeable stamp on a liquid saturated poroelastic base is considered. The equations of two-phase Biot media, which take into account inertial and viscous interactions of phases, are used. The asymptotic properties of the contact stress at low vibration frequencies are studied.     

Propagation of a Stationary Shock Layer in the Presence of a Boundary

We are interested in the large-time behavior of a stationary shock layer in the presence of a boundary effect. It is shown that the stationary shock propagates away from the boundary with speed inversely proportional to the time . Consequently, the location of the wave front is of order . We show therefore that the zero-viscosity and time-asymptotic limits do not commute. Our analysis involves detailed point-wise estimates, which allow us to trace the shock. (Accepted April 16, 1996)     

Solution Construction to a Class of Riemann Problems of Multiphase Flow in Porous Media

Transport in Porous Media - Fluid displacement in porous media can usually be formulated as a Riemann problem. Finding the solution to such a problem helps shed light on the dynamics of flow and...     

3D Passive Walkers: Finding Periodic Gaits in the Presence of Discontinuities

This paper studies repetitive gaits found in a 3D passivewalking mechanism descending an inclined plane. By using directnumerical integration and implementing a semi-analytical scheme forstability analysis and root finding, we follow the corresponding limitcycles under parameter variations. The 3D walking model, which is fullydescribed in the paper, contains both force discontinuities andimpact-like instantaneous changes of state variables. As a result, thestandard use of the variational equations is suitably modified. Theproblem of finding initial conditions for the 3D walker is solved bystarting in an almost planar configuration, making it possible to useparameters and initial conditions found for planar walkers. The walkeris gradually transformed into a 3D walker having dynamics in all spatialdirections. We present such a parameter variation showing the stabilityand the amplitude of the hip sway motion. We also show the dependence ofgait cycle measurements, such as stride time, stride length, averagevelocity, and power consumption, on the plane inclination. The paperconcludes with a discussion of some ideas on how to extend the present3D walker using the tools derived in this paper.     

Stability of Contact Discontinuities for the 1-D Compressible Navier-Stokes Equations

In this paper, we study the large-time asymptotic behavior of solutions of the one-dimensional compressible Navier-Stokes system toward a contact discontinuity, which is one of the basic wave patterns for the compressible Euler equations. It is proved that such a weak contact discontinuity is a metastable wave pattern, in the sense introduced in [24], for the 1-D compressible Navier-Stokes system for polytropic fluid by showing that a viscous contact wave, which approximates the contact discontinuity on any finite-time interval for small heat conduction and then runs away from it for large time, is nonlinearly stable with a uniform convergence rate provided that the initial excess mass is zero. This result is proved by an elaborate combination of elementary energy estimates with a weighted characteristic energy estimate, which makes full use of the underlying structure of the viscous contact wave.