A case of strong nonlinearity: Intermittency in highly turbulent flows

It has long been suspected that flows of incompressible fluids at large or infinite Reynolds number (namely at small or zero viscosity) may present finite time singularities. We review briefly the theoretical situation on this point. We discuss the effect of a small viscosity on the self-similar solution to the Euler equations for inviscid fluids. Then we show that single-point records of velocity fluctuations in the Modane wind tunnel display correlations between large velocities and large accelerations in full agreement with scaling laws derived from Leray's equations (1934) for self-similar singular solutions to the fluid equations. Conversely, those experimental velocity–acceleration correlations are contradictory to the Kolmogorov scaling laws.     

Inverse laminar boundary-layer problems with assigned shear. The mechul function revisited

An inverse method is presented which accurately determines the pressure distribution for assigned wall shear in a two-dimensional, laminar, incompressible boundary layer. The method reformulates the mechul function scheme of Cebeci and Keller to produce a stable solution in the marching direction and to increase accuracy in the normal direction. In the reformulation a modified pressure gradient parameter variation in the normal direction is used in conjunction with three-point backward differences for streamwise derivatives and fourth-order accurate splines for normal derivatives. The resulting spline-finite difference equations are solved by Newton-Raphson iteration together with partial pivoting. Numerical solutions are presented for self-similar and non self-similar flows and compared with published results.     

Self-similar collapse of a circular cavity of a power-law liquid

Using the lubrication approximation we investigate the self-similar axisymmetric flow of a power-law liquid towards a central circular cavity. It is shown that this problem has a self-similar solution of the second kind. The self-similarity exponent is found by solving a non-linear eigenvalue problem arising from the requirement that the integral curve that represents the solution must join the appropriate singular points in the phase plane of the governing equation. The eigenvalues for different values of the rheological index are computed. Numerical integration of the equations allows us to determine the shape of the solution in terms of the physical variables. We make a detailed analysis of the influence of the rheology on the properties of the solutions.     

A new numerical technique for the analysis of parametrically excited nonlinear systems

A new computational scheme using Chebyshev polynomials is proposed for the numerical solution of parametrically excited nonlinear systems. The state vector and the periodic coefficients are expanded in Chebyshev polynomials and an integral equation suitable for a Picard-type iteration is formulated. A Chebyshev collocation is applied to the integral with the nonlinearities reducing the problem to the solution of a set of linear algebraic equations in each iteration. The method is equally applicable for nonlinear systems which are represented in state-space form or by a set of second-order differential equations. The proposed technique is found to duplicate the periodic, multi-periodic and chaotic solutions of a parametrically excited system obtained previously using the conventional numerical integration schemes with comparable CPU times. The technique does not require the inversion of the mass matrix in the case of multi degree-of-freedom systems. The present method is also shown to offer significant computational conveniences over the conventional numerical integration routines when used in a scheme for the direct determination of periodic solutions. Of course, the technique is also applicable to non-parametrically excited nonlinear systems as well.     

Global existence of solutions to the periodic initial value problems for two-dimensional Newton-Boussinesq equations

A class of periodic initial value problems for two-dimensional Newton- Boussinesq equations are investigated in this paper. The Newton-Boussinesq equations are turned into the equivalent integral equations. With iteration methods, the local existence of the solutions is obtained. Using the method of a priori estimates, the global existence of the solution is proved.     

Some problems in the theory of plane jets of conducting fluid

Using the boundary-layer equations as a basis, the author considers the propagation of plane jets of conducting fluid in a transverse magnetic field (noninductive approximation).The propagation of plane jets of conducting fluid is considered in several studies [1–12]. In the first few studies jet flow in a nonuniform magnetic field is considered; here the field strength distribution along the jet axis was chosen in order to obtain self-similar solutions. The solution to such a problem given a constant conductivity of the medium is given in [1–3] for a free jet and in [4] for a semibounded jet; reference [5] contains a solution to the problem of a free jet allowing for the dependence of conductivity on temperature. References [6–8] attempt an exact solution to the problem of jet propagation in any magnetic field. An approximate solution to problems of this type can be obtained by using the integral method. References [9–10] contain the solution obtained by this method for a free jet propagating in a uniform magnetic field.The last study [10] also gives a comparison of the exact solution obtained in [3] with the solution obtained by the integral method using as an example the propagation of a jet in a nonuniform magnetic field. It is shown that for scale values of the jet velocity and thickness the integral method yields almost-exact values. In this study [10], the propagation of a free jet is considered allowing for conduction anisotropy. The solution to the problem of a free jet within the asymptotic boundary layer is obtained in [1] by applying the expansion method to the small magnetic-interaction parameter. With this method, the problem of a turbulent jet is considered in terms of the Prandtl scheme. The Boussinesq formula for the turbulent-viscosity coefficient is used in [12].This study considers the dynamic and thermal problems involved with a laminar free and semibounded jet within the asymptotic boundary layer, propagating in a magnetic field with any distribution. A system of ordinary differential equations and the integral condition are obtained from the initial partial differential equations. The solution of the derived equations is illustrated by the example of jet propagation in a uniform magnetic field. A similar solution is obtained for a turbulent free jet with the turbulent-exchange coefficient defined by the Prandtl scheme.     

Solution to the problem of a swirling jet from an annular orifice

The problem of the propagation of a laminar immersed fan jet with swirling was considered in [1–3]. In [1], the jet source scheme was used to find a self-similar solution for a weakly swirling jet. An attempt to solve by an integral method the analogous problem for a jet emanating from a slit of finite size was made in [2]. In [3], the equations of motion for a jet with arbitrary swirling were reduced under a number of assumptions to the equations that describe the flow of a flat immersed jet. This paper gives the numerical solution to the problem of the propagation of a radial jet emanating with arbitrary swirling from a slit of finite size and an analytic solution for the main section of the jet.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 49–54, March–April, 1991.     

Concentration of electric field near electrodes on piezoelectric layer

This paper aims at the study of electric field concentration near the edge of the electrodes for a piezoelectric layer. For simplicity, attention is focussed on a piezoelectric layer that is bonded to a rigid, conductive substrate and then covered by a pair of conductive electrodes on the surface. The analysis is performed within the scope of linear piezoelectricity. Fourier transform technique is applied to reduce the boundary value problem to a pair of dual integral equations. An analytical solution to the integral equations is obtained by using an iteration method. Numerical results for distributions of stresses and electric fields are displayed and discussed.     

An analytical theory of heated duct flows in supersonic combustors

One-dimensional analytical theory is developed for supersonic duct flow with variation of cross section, wall friction, heat addition, and relations between the inlet and outlet flow parameters are obtained. By introducing a selfsimilar parameter, effects of heat releasing, wall friction, and change in cross section area on the flow can be normalized and a self-similar solution of the flow equations can be found. Based on the result of self-similar solution, the sufficient and necessary condition for the occurrence of thermal choking is derived. A relation of the maximum heat addition leading to thermal choking of the duct flow is derived as functions of area ratio, wall friction, and mass addition, which is an extension of the classic Rayleigh flow theory, where the effects of wall friction and mass addition are not considered. The present work is expected to provide fundamentals for developing an integral analytical theory for ramjets and scramjets.     

定态火焰在可燃预混气中产生的压力波

火焰在可燃预混气中传播时，在火焰面前方产生一道压力波。忽略点火及火焰的初期加速，仅考虑火焰达到稳定传播速度的情况。用Ｏｐｅｎｈｅｉｍ自相似解分析流场，得到相应的控制方程及定解条体；用自适应步长的四阶Ｒｕｎｇｅ－Ｋｕｔｔａ法对方程积分，讨论了流场压力波结构及弱激波近似声波解；认为火焰为间断面，能量释放在火焰面后瞬时完成。利用火焰面两侧的能量关系，得到了火焰位置、燃速及对应Ｃ－Ｊ条件的火焰位置、Ｃ－Ｊ燃速。     

Self-similar problems of elastic contact for non-convex punches

Self-similar problems of contact for non-convex punches are considered. The non-convexity of the punch shapes introduces differences from the traditional self-similar contact problems when punch profiles are convex and their shapes are described by homogeneous functions. First, three-dimensional Hertz type contact problems are considered for non-convex punches whose shapes are described by parametric-homogeneous functions. Examples of such functions are numerous including both fractal Weierstrass type functions and smooth log-periodic sine functions. It is shown that the region of contact in the problems is discrete and the solutions obey a non-classical self-similar law. Then the solution to a particular case of the contact problem for an isotropic linear elastic half-space when the surface roughness is described by a log-periodic function, is studied numerically, i.e. the contact problem for rough punches is studied as a Hertz type contact problem without employing additional assumptions of the multi-asperity approach. To obtain the solution, the method of non-linear boundary integral equations is developed. The problem is solved only on the fundamental domain for the parameter of self-similarity because solutions for other values of the parameter can be obtained by renormalization of this solution. It is shown that the problem has some features of chaotic systems, namely the global character of the solution is independent of fine distinctions between parametric-homogeneous functions describing roughness, while the stress field of the problem is sensitive to small perturbations of the punch shape.     

Motion of a strong shock front in a two-dimensional gas layer with moving internal boundary

We examine the problem of planar one-dimensional motion of a strong shock wave with moving internal boundary in which the initial position of the front, its intensity, the mass of the gas involved in the motion, and the energy contained in this gas are known. The problem is not self-similar and its exact solution, which involves working with partial differential equations, presents serious difficulties. In the following we determine the law of shock-front motion in this problem via the method of [1], which makes it possible to find a system of ordinary differential equations for the problem. The method is based on an initial specification of the power-law coupling between the dimensionless Lagrangian and Eulerian variables and replacement of the energy equation by this coupling and the energy integral. The solution is sought in the first approximation.     

Boundary-layer flow of Reiner–Philippoff fluids past a stretching wedge

This paper presents a numerical analysis of the steady boundary-layer flow of a Reiner–Philippoff fluid induced by a 90° stretching wedge in a variable free stream. The governing partial differential equations are converted into a set of two ordinary differential equations by the use of a similarity transformation. The flow is therefore governed by a stretching velocity parameter λ and two non-Newtonian fluid parameters γ and μ₀. The variation of the skin friction, as well as other flow characteristics, as a function of the governing parameters is presented graphically and tabulated. A stability analysis has also been performed for this self-similar flow based on linear disturbances to the steady similarity solutions. The results presented in this paper reveal that there are no multiple (dual) solutions for the present problem and the unique solution is stable.     

Hydraulic fracture crack propagation driven by a non-Newtonian fluid

The problem of hydraulic fracture crack propagation in a porous medium is considered. The fracture is driven by an incompressible viscous fluid with a power-law rheology of the pseudoplastic type. The fluid seepage is described by an equation generalizing the Darcy law in the hydraulic approximation. It is shown that the system of governing equations has a power-law self-similar solution, whereas, in the limiting cases of low and high fluid saturation in the porous medium, there are some families of power-law or exponential self-similar solutions. The complete self-similar solution is constructed. The effect of the nonlinear rheology of the fracturing fluid on the behavior of the solution is studied. The problem is solved analytically for an arbitrary boundary condition at the crack inlet when the viscous stresses in the non-Newtonian fluid are close to a constant.     

On thermoelastostatics of composites with nonlocal properties of constituents II. Estimation of effective material and field parameters

One considers a linear thermoelastic composite medium, which consists of a homogeneous matrix containing a statistically homogeneous random set of ellipsoidal uncoated or coated heterogeneities. It is assumed that the stress–strain constitutive relations of constituents are described by the nonlocal integral operators, whereas the equilibrium and compatibility equations remain unaltered as in classical local elasticity. The general integral equations connecting the stress and strain fields in the point being considered and the surrounding points are obtained. The method is based on a centering procedure of subtraction from both sides of a known initial integral equation their statistical averages obtained without any auxiliary assumptions such as, e.g., effective field hypothesis implicitly exploited in the known centering methods. In a simplified case of using of the effective field hypothesis for analyzing composites with one sort of heterogeneities, one proves that the effective moduli explicitly depend on both the strain and stress concentrator factor for one heterogeneity inside the infinite matrix and does not directly depend on the elastic properties (local or nonlocal) of heterogeneities. In such a case, the Levin’s (1967) formula in micromechanics of composites with locally elastic constituents is generalized to their nonlocal counterpart. A solution of a volume integral equation for one heterogeneity subjected to inhomogeneous remote loading inside an infinite matrix is proposed by the iteration method. The operator representation of this solution is incorporated into the new general integral equation of micromechanics without exploiting of basic hypotheses of classical micromechanics such as both the effective field hypothesis and “ellipsoidal symmetry” assumption. Quantitative estimations of results obtained by the abandonment of the effective field hypothesis are presented.     

Variational iteration method for two-dimensional steady slip flow in micro-channels

In this paper, the two-dimensional steady slip flow in microchannels is investigated. Research on micro flow, especially on micro slip flow, is very important for designing and optimizing the micro electromechanical system (MEMS). The Navier-Stokes equations for two-dimensional steady slip flow in microchannels are reduced to a nonlinear third-order differential equation by using similarity solution. The variational iteration method (VIM) is used to solve this nonlinear equation analytically. Comparison of the result obtained by the present method with numerical solution reveals that the accuracy and fast convergence of the new method.     

Orthogonal polynomials and determinant formulas of function-valued Padé-type approximation using for solution of integral equations

To solve Fredholm integral equations of the second kind, a generalized linear functional is introduced and a new function-valued Padé-type approximation is defined. By means of the power series expansion of the solution, this method can construct an approximate solution to solve the given integral equation. On the basis of the orthogonal polynomials, two useful determinant expressions of the numerator polynomial and the denominator polynomial for Padé-type approximation are explicitly given. Project supported by the National Natural Science Foundation of China (No.10271074)     

Development of laminar jet flows with zero excess impulse

The problem of the development of a laminar jet of viscous incompressible fluid with zero excess impulse (the wake of a hydrodynamic motor) was investigated for the first time by Birkhoff and Zarantonello [1], who found a self-similar solution to the dynamical problem for the case of a two-dimensional laminar wake. The problem of the development of turbulent wakes of hydrodynamic motors in the near and far flow regions was solved by Ginevskii [2] on the basis of an integral method. In the present paper, the method of asymptotic expansions is used on the basis of the boundary layer equations to solve nonself-similar problems of the development of laminar jet flows of a viscous incompressible fluid with zero excess impulse. The obtained solution takes into account the influence of the details of the source (finite size of the body and its geometry) and the value of the Prandtl number on the velocity and temperature distribution. In the case of a laminar axi-symmetric wake, a self-similar solution is obtained to the thermal problem, the solution being valid in a wide range of Prandtl numbers.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 27–33, May–June, 1984.     

The influence of magnetic field upon the collapse of a cylindrical shock wave

In the present paper, the problem of propagation of collapsing cylindrical shock wave in an ideal gas permeated by a transverse magnetic field with infinite electrical conductivity is investigated. Here it is assumed that the medium ahead of the shock front is uniform and at rest. Also, its counter pressure concerning the motion of the wave front is neglected. This problem admits a self similar solution of second kind. The similarity exponent has been computed by solving a nonlinear eigenvalue problem and integrating numerically the self-similar equations for various values of adiabatic heat exponent and Cowling number. Numerical computations have been performed to determine the flow field behind the shock wave. The influence of magnetic field strength and adiabatic heat exponent on the flow parameters for various cases is presented.     

Calculation of the boundary layer on the electrically conducting wall of a plane channel

There are presently available quite a large number of works devoted to the study of the motion of an electrically conducting fluid in boundary layers formed on electrodes or on the nonconducting walls of various MHD devices. However, the methods of solving the boundary layer equations in these studies are based on various simplifying assumptions which allow the problem to be reduced to the solution of a system of ordinary differential equations. Thus, in [1] there is imposed on the flow the special magnetic fieldH1/x, which enables the problem to be reduced to the self-similar form, while in the studies of other authors [2, 3] either the solution is sought in the form of expansions in x, or it is assumed that the problem is locally self-similar [4]. In the present paper we construct the solution of the MHD boundary layer equations which is obtained by one of the numerical methods which has long been used for solving the boundary layer equations for a nonconducting fluid.