A first Integral of Navier­Stokes equations for two-dimensional flows

As wellknown, Bernoulli's equation is obtained as the first integral of Euler's equations in the absence of vorticity. Even in case of non-vanishing vorticity, a first integral from Euler's equations is obtained by using the so called Clebsch transformation [1] for inviscid flows. In contrast to this, a generalisation of this procedure towards viscous flows has not been established so far. In the present paper a first integral of Navier-Stokes equations is constructed in the case of two-dimensional flow by making use of an alternative representation of the fields in terms of complex coordinates and introducing a potential representation for the pressure. The associated boundary conditions are also considered. The first integral is a suitable tool for the development of new analytical methods and numerical codes in fluid dynamics. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the use of potential fields in fluid mechanics

In classical fluid mechanics, potential fields have been employed to enable the integration of the equations of motion. As is well known, Bernoulli's equation is obtained as a first integral of Euler's equations in the absence of vorticity and viscosity if the velocity vector is perceived as the gradient of a scalar potential. The so-called Clebsch transformation [1] involving three scalar potentials allows for a further extension to flows with non-vanishing vorticity; the resulting equations turn out to be self-adjoint, allowing for a variational formulation. All attempts in classic literature, however, are restricted to inviscid flows and the finding of a potential representation enabling the integration of the Navier-Stokes equations remains desirable. Progress on this topic was reported by [3, 4] who constructed a first integral of the two-dimensional incompressible Navier-Stokes equations by making use of an auxiliary potential field and a representation of the fields in terms of complex coordinates. The new formulation proved to be useful in numerical applications and moreover, replacing the scalar potential by a tensor potential, the theory can be successfully generalised to encompass three-dimensional Navier-Stokes flow. Related to the first integral a finite element method was presented in [2] based on a formulation involving the velocities and the first order derivatives of the introduced potential. This way the dynamic boundary condition could be incorporated elegantly and the system of equations fitted into the first order system least-squares methodology. However, a promising alternative approach results if one considers the streamfunction and a slightly modified potential field as independent variables. This new approach involves Laplacian operators rather than mixed derivatives and allows for a convenient embodiment of the Neumann conditions on the streamfunction that is in contrast to the original stream function / potential formulation [4]. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Equivalence transformations of the Clebsch equations

We consider the Clebsch equations describing the motion of an ideal incompressible fluid in the absence of mass forces. These equations contain an arbitrary element, namely a function of three variables. We find the equivalence transformations of the variables which act on the arbitrary element and preserve the structure of the equations. The equivalence transformation is pointed out that takes the arbitrary element to zero. Some examples are given.     

Clifford分析中于特征流形上奇异积分方程的正则化

借助于多元复分析的思想，此文利用Clifford分析中于特征流形上奇异积分的两种Poincare＇-Bretrand置换公式，研究特征流形上奇异积分方程的Fredhlom理论，找到了它的正则化算子．     

Continuous and Discrete Clebsch Variational Principles

The Clebsch method provides a unifying approach for deriving variational principles for continuous and discrete dynamical systems where elements of a vector space are used to control dynamics on the cotangent bundle of a Lie group via a velocity map. This paper proves a reduction theorem which states that the canonical variables on the Lie group can be eliminated, if and only if the velocity map is a Lie algebra action, thereby producing the Euler–Poincaré (EP) equation for the vector space variables. In this case, the map from the canonical variables on the Lie group to the vector space is the standard momentum map defined using the diamond operator. We apply the Clebsch method in examples of the rotating rigid body and the incompressible Euler equations. Along the way, we explain how singular solutions of the EP equation for the diffeomorphism group (EPDiff) arise as momentum maps in the Clebsch approach. In the case of finite-dimensional Lie groups, the Clebsch variational principle is discretized to produce a variational integrator for the dynamical system. We obtain a discrete map from which the variables on the cotangent bundle of a Lie group may be eliminated to produce a discrete EP equation for elements of the vector space. We give an integrator for the rotating rigid body as an example. We also briefly discuss how to discretize infinite-dimensional Clebsch systems, so as to produce conservative numerical methods for fluid dynamics.      

Elastoviscoplastic fluid flow in non-circular tubes: Transversal field and interplay of elasticity and plasticity

An analytical study of elastoviscoplastic fluid flow in tubes of non-circular cross section is presented. The constitutive structure of the fluid is described by a linear frame invariant combination of the Phan-Thien−Tanner model of viscoelastic fluids and the Bingham model of plastic fluids. Non-circular tube cross sections are modeled by the shape factor method a one-to-one mapping of the circular base contour into a wide spectrum family of arbitrary tube contours. Field variables are expanded into asymptotic series in terms of the elasticity measure, the Weissenberg number We, coupled with an asymptotic expansion in terms of the geometrical mapping parameter ε leading to a set of hierarchical momentum balance equations which are solved successively up to and including the third order in We when the secondary field appears for the first time. The computational algorithm developed is applied to the study of the non-rectilinear flow in tubes with triangular and square cross sections. We find that the presence of the yield stress dampens the intensity of the purely viscoelastic vortices, the higher the yield stress the lower the intensity of the vortices in the cross-section, and the further away the vortices are from the center of the cross section as compared to the purely viscoelastic vortices. The results also evidence that viscoelasticity increases the axial flow for given viscoplastic conditions and pressure drop, and consequently increases the rate of flow, a phenomenon that may find applications in optimizing material transportation.     

Clebsch transformation of hydrodynamic equations

A geometrïc meaning of the Clebsch transformation of the hydrodynamic equations for an ideal fluid is established, as well as its relation to the reduction of differential 1-forms to the Darboux normal form. The corresponding kinematic conservation laws are described. The well-known conservation law for the number of quasiparticles is considered as a particular example.     

The method of potential functions in problems of the theory of elasticity for bodies with a defect

The method of potential functions using a Fourier transformation in the class of slowly increasing distributions, corresponding to the classical method of complex potentials, is proposed for solving well-known problems of the theory of elasticity for bodies with a defect. It is shown that when a Fourier transformation with respect to all the spatial variables is used, the solution of the dynamic problem of the theory of elasticity can also be represented in terms of a jump in the stresses and displacements at the defect. The correctness of the transformed problem is considered (in terms of an analogue of the Lopatinskii condition). The solution of the system of Helmholtz equations, to which the system of Lamé equations is reduced in the case of the two-dimensional dynamic problem, is expressed in terms of the jump in the stresses and displacements at the defect as a result of solving the corresponding singular integral equations.     

Exact solutions of some mixed problems of uncoupled thermoelasticity for a truncated hollow circular cone with a groove along the generatrix

An elastic body of finite dimensions in the form of a truncated hollow circular cone with a groove along the generatrix is considered. The uncoupled problem of thermoelasticity is formulated for this body for different types of boundary conditions on all the surfaces. These are the conditions for specifying the displacements or sliding clamping on surfaces with fixed angular coordinates and the conditions for specifying the stresses on surfaces with a fixed radial coordinate (shear stresses are assumed to be zero). It is assumed that the temperature is a specified function of all the spherical coordinates. Some auxiliary functions, related to the displacements, are introduced first, and equations for these functions are then derived using Lamé's equations. A finite integral Fourier transformation with respect to one of the angular variables is then employed. After this, by solving certain Sturm-Liouville problems, a new integral transformation is constructed and is applied to the equations with respect to the other angular variable. As a result a one-dimensional system of differential equations is obtained, to solve which an integral Mellin transformation is employed in a special way. Finally, exact solutions of some problems of thermoelasticity are constructed in series for this body.     

Optimal fields for integral equations

The theory of optimal fields is developed for optimal control problems in which the state variables are determined by integral equations. The Hilbert integral is considered and the Hamilton-Jacobi equations are derived. The results obtained contain two maximum principles as special cases; one reflecting the special character of field theory and one corresponding to the results previously obtained by use of variations.     

A Lagrangian formulation for steady three-dimensional Newton-Busemann flow

A new Lagrangian formulation for steady three dimensional inviscid flow over rigid bodies is developed. First, the continuity equation is eliminated by the use of two stream functions. This is followed by a transformation to new independent variables, two of which are these stream functions and the third one is a Lagrangian time distinct from the Eulerian time. This Lagrangian formulation with the use of Lagrangian time requires only three independent variables and allows the free boundary problem of flow with shock wave to be rendered a fixed boundary one thereby making it easier to solve. In the Newtonian limit the governing equations reduce to the geodesic equations of the body surface. For flow past two-dimensional bodies, bodies of revolution, and conical bodies and wings, the problems are solved to within quadrature. All known Newtonian flow solutions are found to be special cases of the present theory.     

An asymptotic method of solving transient dynamic contact problems of the theory of elasticity for a strip

An asymptotic method is proposed for solving transient dynamic contact problems of the theory of elasticity for a thin strip. The solution of problems by means of the integral Laplace transformation (with respect to time) and the Fourier transformation (with respect to the longitudinal coordinate) reduces to an integral equation in the form of a convolution of the first kind in the unknown Laplace transform of contact stresses under the punch. The zeroth term of the asymptotic form of the solution of the integral equation for large values of the Laplace parameter is constructed in the form of the superposition of solutions of the corresponding Wiener-Hopf integral equations minus the solution of the corresponding integral equation on the entire axis. In solving the Wiener-Hopf integral equations, the symbols of the kernel of the integral equation in the complex plane is presented in special form — in the form of uniform expansion in terms of exponential functions. The latter enables integral equations of the second kind to be obtained for determining the Laplace-Fourier transform of the required contact stresses, which, in turn, is effectively solved by the method of successive approximations. After Laplace inversion of the zeroth term of the asymptotic form of the solution of the integral equations, the asymptotic solution of the transient dynamic contact problem is determined. By way of example, the asymptotic solution of the problem of the penetration of a plane punch into an elastic strip lying without friction on a rigid base is given. Formulae are derived for the active elastic resistance force on the punch of a medium preventing the penetration of the punch, and the law of penetration of the punch into the elastic strip is obtained, taking into account the elastic stress wave reflected from the strip face opposite the punch and passing underneath it.     

Penrose transformation and complex integral geometry

A survey is given of results on the representation of solutions of systems of massless equations in terms of solutions of the Cauchy-Riemann equations on the space of Penrose twisters. A detailed introduction to the theory of twistors is given. The basic integral transform — the Penrose transformation — is investigated by means of complex integral geometry. Other approaches to the theory of twistors are discussed and compared. In addition, the Yang-Mills equation is considered from the point of view of nonlinear Cauchy-Riemann equations.     

Transient response of dissimilar piezoelectric layers with multiple interacting interface cracks

In this study, we examine the dynamic behavior of two bonded dissimilar piezoelectric layers containing multiple interfacial cracks subjected to electro-mechanical impact loading. The problem was formulated through Fourier transformation into singular integral equations in which the unknown variables are the jumps of displacement and electric potential across the crack surface in the Laplace transform domain. The resulting integral equations together with the corresponding single-valued conditions are solved numerically for the densities of electro-elastic dislocations on a crack surface. The dynamic field intensity factors and dynamic energy release rate (DERR) history are obtained for both permeable and impermeable crack. The stress field is also obtained for the interface crack under impact loads. The results show that the field intensity factors at the crack tips and dynamic energy release rate depend on the interfacial crack geometry, electromechanical coupling and the electric boundary conditions on the crack surface.     

The Far-Field Equations in Linear Elasticity for Disconnected Rigid Bodies and Cavities

In this paper the far-field equations in linear elasticity for scattering from disjoint rigid bodies and cavities are considered. The direct scattering problem is formulated in differential and integral form. The boundary integral equations are constructed using a combination of single- and double-layer potentials. Using a Fredholm type theory it is proved that these boundary integral equations are uniquely solvable. Assuming that the incident field is produced by a superposition of plane incident waves in all directions of propagation and polarization it is established that the scattered field is also expressed as the superposition of the corresponding scattered fields. A pair of integral equations of the first kind which hold independently of the boundary conditions are constructed for the far-field region. The properties of the Herglotz functions are used to derive solvability conditions for the far-field equations. It is also proved that the far-field operators, in terms of which we can express the far-field equations, are injective and have dense range. An analytical example for spheres illuminates the theoretical results.     

A Note on a Singular Integral Equation

A direct method is presented for solving a singular integralequation which is a generalization of one occurring in wave-guidetheory and in the theory of dislocations. The present methodavoids elaborate application of the complex variable methodsnormally associated with such singular integral equations.     

Difference Potentials Method and its Applications

The purpose of this article is to acquaint the reader with the general concepts and capabilities of the Difference Potentials Method (DPM). DPM is used for the numerical solution of boundary-value and some other problems in difference and differential formulations. Difference potentials and DPM play the same role in the theory of solutions of linear systems of difference equations on multi-dimensional non-regular meshes as the classical Cauchy integral and the method of singular integral equations do in the theory of analytical functions (solutions Cauchy-Riemann system). The application of DPM to the solution of problems in difference formulation forms the first aspect of the method. The second aspect of the DPM implementation is the discretization and numerical solution of the Calderon-Seeley boundary pseudo-differential equations. The latter are equivalent to elliptical differential equations with variable coefficients in the domain; they are written making no use of fundamental solutions and integrals. Because of this fact ordinary methods for discretization of integral equations cannot be applied in this case. Calderon-Seeley equations have probably not been used for computations before the theory of DPM appeared. This second aspect for the implementation of DPM is that it does not require difference approximation on the boundary conditions of the original problem. The latter circumstance is just the main advantage of the second aspect in comparison with the first one. To begin with, we put forward and justify the main constructions and applications of DPM for problems connected with the Laplace equation. Further, we also outline the general theory and applications: both those already realized and anticipated.     

Remarks on the Linearization of Differential Equations

The Hamiltonian formalism and the theory of canonical transformationsare used in this paper, first of all, to show that, given anordinary non-linear differential equation it is always possiblein principle to find a variable transformation reducing it toa linear equation, or a system of linear equations. The proofgiven is not to be construed as a general practical method forfinding this transformation; it merely shows that such a transformationmust always exist. It is suggested that this may also hold true for partial differentialequations. The conjecture is made plausible, in two cases, bythe use of canonical transformation procedures for linearizingsimple non-linear partial differential equations—one beinga slight generalization of Burger's equation and the other anequation in three independent variables reminiscent of the Eulerequations for fluid flow.     

Topological analysis of an integrable system related to the rolling of a ball on a sphere

A new integrable system describing the rolling of a rigid body with a spherical cavity on a spherical base is considered. Previously the authors found the separation of variables for this system on the zero level set of a linear (in angular velocity) first integral, whereas in the general case it is not possible to separate the variables. In this paper we show that the foliation into invariant tori in this problem is equivalent to the corresponding foliation in the Clebsch integrable system in rigid body dynamics (for which no real separation of variables has been found either). In particular, a fixed point of focus type is possible for this system, which can serve as a topological obstacle to the real separation of variables.     

A General Theory of Resonance Between Weakly Coupled Oscillatory Systems

Resonant oscillations are investigated between several fullynon-linear oscillatory systems which are subject to a weak coupling.The equations of motion are obtained from a Hamiltonian whichis first expressed in terms of angle action variables, and thenaveraged. It is shown that this results in a reduced set ofequations, the reduction depending on the number of resonancesbetween the various systems. For example a resonance in a systemwith two degrees of freedom results in equations which are mathematicallyequivalent to those for one degree of freedom. The theory isillustrated by application to the forced oscillations of a simplependulum and to the resonant interaction between the two modesof oscillation of a double pendulum.