Harmonic Maps and Ideal Fluid Flows

Using harmonic maps we provide an approach towards obtaining explicit solutions to the incompressible two-dimensional Euler equations. More precisely, the problem of finding all solutions which in Lagrangian variables (describing the particle paths of the flow) present a labelling by harmonic functions is reduced to solving an explicit nonlinear differential system in \mathbb Cⁿ{\mathbb {C^n}} with n = 3 or n = 4. While the general solution is not available in explicit form, structural properties of the system permit us to identify several families of explicit solutions.     

On the solutions in the linear theory of micropolar viscoelasticity

In the present paper the linear theory of micropolar viscoelasticity is considered. The explicit expression of fundamental solution of the system of equations of steady vibrations is constructed by means of elementary functions and its basic properties are established. The Green's formulas in the considered theory are obtained. The formulas of integral representations of Somigliana-type of regular vector and regular (classical) solution are presented. The representation formulas of Galerkin-type solution of the system of nonhomogeneous equations and of the general solution of the system of homogeneous equations by means of eight metaharmonic functions are presented. The completeness of these solutions is proved.     

NEW THEORY FOR EQUATIONS OF NON-FUCHSIAN TYPE——REPRESENTATION THEOREM OF TREE SERIES SOLUTION (Ⅱ)

Our main result consists in proving the representation theorem. Irregular integral is a new type of analytic function, represented by a compound Taylor-Fourier tree series, in which each coefficient of the Fourier series is a Taylor series, while the Taylor coefficients are tree series in terms of equations parameters, higher order correction terms to each coefficient having tree structure with inexhaustible proliferation.The solution obtained is proved to be convergent absolutely and uniformly in the region defined by coefficient functions of the original equation, provided the structure parameter is less than unity. Direct substitution shows that our tree series solution satisfies the equation explicity generation by generation.As compared with classical theory our method not only furnishes explicit expression of irregular integral, leading to the solution of Poincare problem, but also provides possibility of extending the scope of investigation for analytic theory to equations with various kinds of singularities in a unifying way.Exact explicit analytic expression for irregular integrals can be obtained by means of correspondence principle.It is not difficult to prove the convergence of the tree series solution obtained. Direct substitution shows it satisfies the equation.The tree series is automorphic, which agrees completely with Poincaré’s conjecture.     

Plane Strain Bending of Cylindrical Sectors of Admissible Compressible Hyperelastic Materials

In the theory of nonlinear elasticity of rubber-like materials, if a homogeneous isotropic compressible material is described by a strain–energy function that is a homogeneous function of the principal stretches, then the equations of equilibrium for axisymmetric deformations reduce to a separable first-order ordinary differential equation. For a particular class of such strain–energy functions, this property is used to obtain a general parametric solution to the equilibrium equation for plane strain bending of cylindrical sectors. Specification of the arbitrary function that appears in such strain–energy functions yields some parametric solutions. In some cases, the parameter can be eliminated to yield closed-form solutions in implicit or explicit form. Other possible forms for the arbitrary constitutive function that are likely to yield such solutions are also indicated.     

New development in Poincaré's problem of irregular integrals

In connection with non-Fuchsian equations Poincaré has made an important conclusion: It is impossible to obtain explicit expressions of irregular integrals .To elucidate the essence of Poincaré's problem, we establish correspondence theorem. Irregular integrals are analytic functions of new kind, possessing tree structure; part of which can be represented by conventional recursive series, while its remaining part is expressed by the so-called tree series, not subjecting to any recursive relation at all.In contrast to the numerical solution calculated by infinite determinant of classical theory (Hill-Poincaré-von Koch), our method yields naturally exact analytic solution in explicit form. The method proposed may be used to construct a unifying theory for general equations with variable coefficients, having various kinds of singularities as singular lines.The significance of Poincaré conjecture is discussed, the tree series obtained belong to higher automorphic functions.     

On solutions of the problem in stresses with the use of Maxwell stress functions

The spatial problems of elasticity are mainly solved in displacements [1, 2], i.e., the Lamé equations are taken as the initial equations. This is related to the lack of general solutions for the system of basic equations of elasticity expressed in stresses. In this connection, a new variational statement of the problem in stresses was developed in [3, 4]; this statement consists in solving six generalized equations of compatibility for six independent components of the stress tensor, while the three equilibrium equations are transferred to the set of boundary conditions. This method is more convenient for the numerical solution of problems in stresses and has been tested when solving various boundary value problems. In the present paper, analyzing the completeness of the Saint-Venant identities and using the Maxwell stress functions, we obtain a new resolving system of three differential equations of strain compatibility for the three desired stress functions φ, ξ, and ψ. This system is an alternative to the three Lamé equilibrium equations for three desired displacement components u, v, w and is simpler in structure. Moreover, both of these systems of resolving equations can be solved by the new recursive-operator method [5, 6]. In contrast to well-known methods for constructing general solutions of linear differential equations and their systems, the solutions obtained by the recursive-operator method are constructed as operator-power series acting on arbitrary analytic functions of real variables (not necessarily harmonic), and the series coefficients are determined from recursive relations (matrix in the case of systems of equations). The arbitrary functions contained in the general solution can be determined directly either from the boundary conditions (the obtained system of inhomogeneous equations with a right-hand side can also be solved by the recursive-operator method [6]) or by choosing them from various classes of analytic functions (elementary, special); a complete set of particular solutions can be obtained in the same function classes, and the coefficients of linear combinations of particular solutions can be determined by the Trefftz method, the least-squares method, and the collocation method.     

Young Measures Generated by Ideal Incompressible Fluid Flows

In their seminal paper, DiPerna and Majda (Commun Math Phys 108(4):667–689, 1987) introduced the notion of a measure-valued solution for the incompressible Euler equations in order to capture complex phenomena present in limits of approximate solutions, such as persistence of oscillation and development of concentrations. Furthermore, they gave several explicit examples exhibiting such phenomena. In this paper we show that any measure-valued solution can be generated by a sequence of exact weak solutions. In particular this gives rise to a very large, arguably too large, set of weak solutions of the incompressible Euler equations.     

Removable Singularities for Fully Nonlinear Elliptic Equations

We obtain a removability result for the fully nonlinear uniformly elliptic equations F(D ² u)+f(u)=0. The main theorem states that every solution to the equation in a punctured ball (without any restrictions on the behaviour near the centre of the ball) is extendable to the solution in the entire ball provided the function f satisfies certain sharp conditions depending on F. Previously such results were known for linear and quasilinear operators F. In comparison with the semi- or quasilinear theory the techniques for the fully nonlinear equations are new and based on the use of the viscosity notion of generalised solution rather than the distributional or the weak solutions. Accepted May 3, 2000?Published online November 16, 2000     

Asymptotic Axially Symmetric Deformations for Perfectly Elastic Neo-Hookean and Mooney Materials

For axially symmetric deformations of the perfectly elastic neo-Hookean and Mooney materials, formal series solutions are determined in terms of expansions in appropriate powers of 1/R, where R is the cylindrical polar coordinate for the material coordinates. Remarkably, for both the neo-Hookean and Mooney materials, the first three terms of such expansions can be completely determined analytically in terms of elementary integrals. From the incompressibility condition and the equilibrium equations, the six unknown deformation functions, appearing in the first three terms can be reduced to five formal integrations involving in total seven arbitrary constants A, B, C, D, E, H and k ², and a further five integration constants, making a total of 12 integration constants for the deformation field. The solutions obtained for the neo-Hookean material are applied to the problem of the axial compression of a cylindrical rubber tube which has bonded metal end-plates. The solution so determined is approximate in two senses; namely as an approximate solution of the governing equations and for which the stress free and displacement boundary conditions are satisfied in an average manner only. The resulting load-deflection relation is shown graphically. The solution so determined, although approximate, attempts to solve a problem not previously tackled in the literature.      

A σ‐coordinate three‐dimensional numerical model for surface wave propagation

A three‐dimensional numerical model based on the full Navier–Stokes equations (NSE) in σ‐coordinate is developed in this study. The σ‐coordinate transformation is first introduced to map the irregular physical domain with the wavy free surface and uneven bottom to the regular computational domain with the shape of a rectangular prism. Using the chain rule of partial differentiation, a new set of governing equations is derived in the σ‐coordinate from the original NSE defined in the Cartesian coordinate. The operator splitting method (Li and Yu, Int. J. Num. Meth. Fluids 1996; 23 : 485–501), which splits the solution procedure into the advection, diffusion, and propagation steps, is used to solve the modified NSE. The model is first tested for mass and energy conservation as well as mesh convergence by using an example of water sloshing in a confined tank. Excellent agreements between numerical results and analytical solutions are obtained. The model is then used to simulate two‐ and three‐dimensional solitary waves propagating in constant depth. Very good agreements between numerical results and analytical solutions are obtained for both free surface displacements and velocities. Finally, a more realistic case of periodic wave train passing through a submerged breakwater is simulated. Comparisons between numerical results and experimental data are promising. The model is proven to be an accurate tool for consequent studies of wave‐structure interaction. Copyright © 2002 John Wiley & Sons, Ltd.     

Three-dimensional Green's functions in anisotropic trimaterials

Three-dimensional elastostatic Green's functions in anisotropic trimaterials are derived, for the first time, by applying the generalized Stroh's formalism and Fourier transforms. The Green's functions are expressed as a series summation with the first term corresponding to the full-space solution and other terms to the image solutions due to the interfaces. The most remarkable feature of the present solution is that the image solutions can be expressed by a simple line integral over a finite interval [0,2π]. By partitioning the trimaterial Green's function into a full-space solution and a complementary part, the line integral involves only regular functions if the singularity is within one of the three materials, being treated analytically owning to the explicit expression of the full-space solution. When the singularity is on the interface, which occurs if the field and source points are both on the same interface, the involved singularity is handled with the interfacial Green's functions.A numerical example is presented for a trimaterial system made of two anisotropic half spaces bonded perfectly by an isotropic adhesive layer, showing clearly the effect of material layering on the Green's displacements and stresses. Furthermore, by comparing the present Green's solution to the direct (two-dimensional) 2D integral expression which is also derived in this paper, it is shown that, the computational time for the calculation of the Green's function can be substantially reduced using the present solution, instead of the direct 2D integral method.     

Penny-shaped crack in transversely isotropic piezoelectric materials

Using a method of potential functions introduced successively to integrate the field equations of three-dimensional problems for transversely isotropic piezoelectric materials, we obtain the so-called general solution in which the displacement components and electric potential functions are represented by a singular function satisfying some special partial differential equations of 6th order. In order to analyse the mechanical-electric coupling behaviour of penny-shaped crack for above materials, another form of the general solution is obtained under cylindrical coordinate system by introducing three quasi-harmonic functions into the general equations obtained above. It is shown that both the two forms of the general solutions are complete. Furthermore, the mechanical-electric coupling behaviour of penny-shaped crack in transversely isotropic piezoelectric media is analysed under axisymmetric tensile loading case, and the crack-tip stress field and electric displacement field are obtained. The results show that the stress and the electric displacement components near the crack tip have (r ^−1/2) singularity. The project supported by the Natural Science Foundation of Shaanxi Province, China     

Entropy Solutions for Nonlinear Degenerate Problems

We consider a class of elliptic-hyperbolic degenerate equations g(u)-Db(u) +\divgf(u) = fg(u)-\Delta b(u) +\divg\phi (u) =f with Dirichlet homogeneous boundary conditions and a class of elliptic-parabolic-hyperbolic degenerate equations g(u)_t-Db(u) +\divgf(u) = fg(u)_t-\Delta b(u) +\divg\phi (u) =f with homogeneous Dirichlet conditions and initial conditions. Existence of entropy solutions for both problems is proved for nondecreasing continuous functions g and b vanishing at zero and for a continuous vectorial function J satisfying rather general conditions. Comparison and uniqueness of entropy solutions are proved for g and b continuous and nondecreasing and for J continuous.     

The derivation of a thick and thin plate formulation without ad hoc assumptions

For the plate formulation considered in this paper, appropriate three-dimensional elasticity solution representations for isotropic materials are constructed. No a priori assumptions for stress or displacement distributions over the thickness of the plate are made. The strategy used in the derivation is to separate functions of the thickness variable z from functions of the coordinates x and y lying in the midplane of the plate. Real and complex 3-dimensional elasticity solution representations are used to obtain three types of functions of the coordinates x, y and the corresponding differential equations. The separation of the functions of the thickness coordinate can be done by separately considering homogeneous and nonhomogeneous boundary conditions on the upper and lower faces of the plate. One type of the plate solutions derived involves polynomials of the thickness coordinate z. The other two solution forms contain trigonometric and hyperbolic functions of z, respectively. Both bending and stretching (or in-plane) solutions are included in the derivation.     

Singularity formation in a collisionless plasma could occur only at high velocities

The global existence problem is studied for regular solutions of the relativistic Vlasov-Maxwell equations. If it is assumed that the plasma density vanishes a priori for velocities near the speed of light, then regular solutions with arbitrary initial data exist in all of space and time. This assumption is either postulated for a solution or is arranged for all solutions through a modification of the equations themselves.     

Intermittency and Regularity Issues in 3<Emphasis Type="Italic">D</Emphasis> Navier-Stokes Turbulence

Two related open problems in the theory of 3D Navier-Stokes turbulence are discussed in this paper. The first is the phenomenon of intermittency in the dissipation field. Dissipation-range intermittency was first discovered experimentally by Batchelor and Townsend over fifty years ago. It is characterized by spatio-temporal binary behaviour in which long, quiescent periods in the velocity signal are interrupted by short, active ‘events’ during which there are violent fluctuations away from the average. The second and related problem is whether solutions of the 3D Navier-Stokes equations develop finite time singularities during these events. This paper shows that Leray’s weak solutions of the three-dimensional incompressible Navier-Stokes equations can have a binary character in time. The time-axis is split into ‘good’ and ‘bad’ intervals: on the ‘good’ intervals solutions are bounded and regular, whereas singularities are still possible within the ‘bad’ intervals. An estimate for the width of the latter is very small and decreases with increasing Reynolds number. It also decreases relative to the lengths of the good intervals as the Reynolds number increases. Within these ‘bad’ intervals, lower bounds on the local energy dissipation rate and other quantities, such as ||u(·, t)||_∞ and ||∇u(·, t)||_∞, are very large, resulting in strong dynamics at sub-Kolmogorov scales. Intersections of bad intervals for n≧1 are related to the potentially singular set in time. It is also proved that the Navier-Stokes equations are conditionally regular provided, in a given ‘bad’ interval, the energy has a lower bound that is decaying exponentially in time.Final version 17 March 05. Original version November 03.     

Three dimensional analytical solution for finite circular cylinders subjected to indirect tensile test

This paper derives a new three-dimensional (3-D) analytical solution for the indirect tensile tests standardized by ISRM (International Society for Rock Mechanics) for testing rocks, and by ASTM (American Society for Testing and Materials) for testing concretes. The present solution for solid circular cylinders of finite length can be considered as a 3-D counterpart of the classical two dimensional (2-D) solutions by Hertz in 1883 and by Hondros in 1959. The contacts between the two steel diametral loading platens and the curved surfaces of a cylindrical specimen of length H and diameter D are modeled as circular-to-circular Hertz contact and straight-to-circular Hertz contact for ISRM and ASTM standards respectively. The equilibrium equations of the linear elastic circular cylinder of finite length are first uncoupled by using displacement functions, which are then expressed in infinite series of some combinations of Bessel functions, hyperbolic functions, and trigonometric functions. The applied tractions are expanded in Fourier–Bessel series and boundary conditions are used to yield a system of simultaneous equations. For typical rock cylinders of 54 mm diameter subjected to ISRM indirect tensile tests, the contact width is in the order of 2 mm (or a contact angle of 4°) whereas for typical asphalt cylinders of 101.6 mm diameter subjected to ASTM indirect tensile tests the contact width is about 10 mm (or a contact angle of 12°). For such contact conditions, 50 terms in both Fourier and Fourier–Bessel series expansions are found sufficient in yielding converged solutions. The maximum hoop stress is always observed within the central portion on a circular section close to the flat end surfaces. The difference in the maximum hoop stress between the 2-D Hondros solution and the present 3-D solution increases with the aspect ratio H/D as well as Poisson’s ratio ν. When contact friction is neglected, the effect of loading platen stiffness on tensile stress in cylinders is found negligible. For the aspect ratio of H/D = 0.5 recommended by ISRM and ASTM, the error in tensile strength may be up to 15% for both typical rocks and asphalts, whereas for longer cylinders with H/D up to 2 the error ranges from 15% for highly compressible materials, and to 60% for nearly incompressible materials. The difference in compressive radial stress between the 2-D Hertz solution or 2-D Hondros solution and the present 3-D solution also increases with Poisson’s ratio and aspect ratio H/D. In summary, the 2-D solution, in general, underestimates the maximum tensile stress and cannot predict the location of the maximum hoop stress which typically locates close to the end surfaces of the cylinder.     

NEW THEORY FOR EQUATIONS OF NON-FUGHSIAN TYPE REPRESENTATION THEOREM OF TREE SERIES SOLUTION (Ⅰ)

In the analytic theory of differential equations the exact explicit analytic solution has not been obtained for equations of the non-Fuchsian type (Poincare's problem). The new theory proposed in this paper for the first time affords a general method of finding exact analytic expres-sion for irregular integrals.By discarding the assumption of formal solution of classical theory,our method consists in deriving a cor-respondence relation from the equation itself and providing the analytic structure of irregular integrals naturally by the residue theorem. Irregular integrals are made up of three parts: noncontracted part,represented by ordinary recursion series,all-and semi-contracted part by the so-called tree series. Tree series solutions belong to analytic function of the new kind with recursion series as the special case only.     

Application of Uniform Asymptotics to the Second Painlevé Transcendent

In this work we propose a new method for investigating connection problems for the class of nonlinear second‐order differential equations known as the Painlevé equations. Such problems can be characterized by the question as to how the asymptotic behaviours of solutions are related as the independent variable is allowed to pass towards infinity along different directions in the complex plane. Connection problems have been previously tackled by a variety of methods. Frequently these are based on the ideas of isomonodromic deformation and the matching of WKB solutions. However, the implementation of these methods often tends to be heuristic in nature and so the task of rigorising the process is complicated. The method we propose here develops uniform approximations to solutions. This removes the need to match solutions, is rigorous, and can lead to the solution of connection problems with minimal computational effort. Our method relies on finding uniform approximations of differ ential equations of the generic form as the complex‐valued parameter . The details of the treatment rely heavily on the locations of the zeros of the function F in this limit. If they are isolated, then a uniform approximation to solutions can be derived in terms of Airy functions of suitable argument. On the other hand, if two of the zeros of F coalesce as , then an approximation can be derived in terms of parabolic cylinder functions. In this paper we discuss both cases, but illustrate our technique in action by applying the parabolic cylinder case to the “classical” connection problem associated with the second Painlevé transcendent. Future papers will show how the technique can be applied with very little change to the other Painlevé equations, and to the wider problem of the asymptotic behavio ur of the general solution to any of these equations. (Accepted May 15, 1997)     

Potential Method in the Linear Theory of Viscoelastic Materials with Voids

In the present paper the linear theory of viscoelasticity for Kelvin–Voigt materials with voids is considered and some basic results of the classical theory of elasticity are generalized. Indeed, the basic properties of plane harmonic waves are established. The explicit expression of fundamental solution of the system of equations of steady vibrations is constructed by means of elementary functions. The Green’s formulas in the considered theory are obtained. The uniqueness theorems of the internal and external basic boundary value problems (BVPs) are proved. The representation of Galerkin type solution is obtained and the completeness of this solution is established. The formulas of integral representations of Somigliana type of regular vector and regular (classical) solution are obtained. The Sommerfeld-Kupradze type radiation conditions are established. The basic properties of elastopotentials and singular integral operators are given. Finally, the existence theorems for classical solutions of the internal and external basic BVPs of steady vibrations are proved by using of the potential method (boundary integral method) and the theory of singular integral equations.