ON A CLASS OF METHOD FOR SOLVING PROBLEMS WITH RANDOM BOUNDARY NOTCHES AND/OR CRACKS(Ⅱ) COMPUTATIONS FOR BOUNDARY NOTCHES

This work is a continuation of the discussion of [1], On a class of method for solving problems with random boundary notches and/or cracks, (I) by C. Ouyang (Appl. Math. & Mech., vol. 1, No. 2, 1980). Here computations for boundary notches are made by using the theory and formulas presented in [1]. In the computation modification is also made for the boundary conditions in parametric plane in [1]. Numerical results for examples show that within ranges of parameter considered in the paper, for example L, the present method in quite workable in practical computations.     

Structure of filtration pseudoturbulence

The problem of large-scale (pseudoturbulent) motion of a fluid in a nonuniform porous medium was formulated in [1]. Since in practice the local porosity (x) is unknown, it may be considered a continuous random point function. The difference of the local values of the porosity (x) from the mean value for the medium as a whole leads to the occurrence of random pseudoturbulent motions of the filtering fluid, which are superposed on the mean filtration flow. The characteristics of the large-scale filtration motion in a medium with this sort of random porosity were considered in detail in [1], where the formal solution is presented for the resulting equations for two-point correlations, based on the use of considerations of spatial invariance. Also presented is a qualitative discussion of the effect of pseudoturbulence of the filtering medium on the transport processes in the medium.We note that the considered problem of pseudoturbulence of a filtering fluid in a nonuniform porous medium does not have anything in common with the statistical problem of the motion of small fluid elements in a broken porous space. The latter problem is interesting in connection with the analysis of the convective diffusion processes in a porous body (both uniform and nonuniform) and, beginning with [3], has been considered in several studies, including [2].In the present study we have used a method for solving the problem which is significantly different in comparison with that of [1], based on the representation of the variation of the local porosity from point to point as a random process with independent increments. This method has the advantage that it permits expressing the required correlation functions in the form of quadratures with arbitrary values of the characteristic parameters. In the following, for simplicity we consider the axisymmetric problem under the assumption that the two-point correlation of the deviations of the local values of the porosity from the mean are representable in the form of an isotropie Gaussian function of the distance between the points. The explicit expressions for the correlations are also written in some approximation and the physical consequences resulting from these assumptions are discussed.     

Again discussing about singular perturbation of general boundary value problem for higher order elliptic equation containing two-parameter

In this paper using the method of The Two-Variable Expansion Procedure [11] we again discuss the construction of asymptotic expression of solution of general boundary value problem for higher order ellitptic equation containing two-parameter whose boundary condition is more general than [1]. We give asymptotic expression of solution as well as the estimation corresponding to the remainder term.     

Electrical parameters of a channel with finite electrodes with allowance for the electrode potential drop

Spatial problems involving the electric field in an MHD channel were formulated in [1] with allowance for the electrode potential drop. It was assumed that the electrode layer had a small thickness, so that relationships on the boundary of the layer could be applied to the surface of the electrode. It was assumed that the electrode potential drop ° could be represented as a function of the current density jn at the electrode in the form of a known function ° =f (j_n) determined experimentally or deduced from the appropriate electrode-layer theory. An approximate method was then put forward for solving such problems by reducing them to the determination of the electric field from a known distribution of the magnetic field and the gas-dynamic parameters. It was shown that when =°/ E is small (E is the characteristic induced or applied potential difference), the solution can be sought in the form of series in powers of . In the zero-order approximation, the electric field is determined without taking into account the electrode processes. The first approximation gives a correction of the order of . The quantity °, which is present in the boundary conditions on the electrode in the first-order approximation, is determined from the current density calculated in the zero-order approximation.One of the problems discussed in [1] was concerned with the electric current in a channel with one pair of symmetric electrodes. Its solution was found in the first approximation in the form of the integral Keldysh-Sedov formula. In this paper we report an analysis of the solution for ° taken in the form of a step function.     

Calculation of heat transfer accompanying hypersonic three-dimensional flow of air over slender blunt-nosed cones

The effective length method [1, 2] has been used to make systematic calculations of the heat transfer for laminar and turbulent boundary layers on slender blunt-nosed cones at small angles of attack ( + 5° in a separationless hypersonic air stream dissociating in equilibrium (half-angles of the cones 0 20°, angles of attack 0 15°, Mach numbers 5 M 25). The parameters of the gas at the outer edge of the boundary layer were taken equal to the inviscid parameters on the surface of the cones. Analysis of the results leads to simple approximate dependences for the heat transfer coefficients.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 173–177, September–October, 1981.     

Heat transfer in laminar flow in a uniformly porous channel with an applied transverse magnetic field

Summary The steady laminar flow of an incompressible, viscous, and electrically conducting fluid between two parallel porous plates with equal permeability has been discussed by Terrill and Shrestha [6]. In this paper, using the solution of [6] for the velocity field, the heat transfer problems of (i) uniform wall temperature and (ii) uniform heat flux at wall are solved.For small suction Reynolds numbers we find that the Nusselt number, with increasing Reynolds number, increases for case (i) and decreases for (ii).Nomenclature stream function - 2h channel width - x, y distances measured parallel, perpendicular to the channel walls - U velocity of fluid in the x direction at x=0 - V constant velocity of suction at the wall - nondimensional distance, y/h - nondimensional distance, x/h - f() function defined in (1) - density - coefficient of kinematic viscosity - R suction Reynolds number, V h/ - Re channel Reynolds number, 4U h/ - B ₀ magnetic induction - electrical conductivity - M Hartmann number, B ₀ h(/)^1/2 - K constant defined in (3) - A constant defined in (5) - 4R/Re - q local heat flux per unit area at the wall - k thermal conductivity - T temperature of the fluid - X –1/ ln(1–) - C _p specific heat at constant pressure - j current density - Pr Prandtl number, C _p/k - P mass transfer Péclet number, R Pr - Pe mass transfer Péclet number, P/ - T ₀ temperature at x=0 - T _H() temperature in the fully developed region - T _h(X, ) temperature in the entrance region - Y _n () eigenfunctions, uniform wall temperature - _n eigenvalues - e() function defined by (24) - B _n 2/3 _n ² - A _n constants defined by (28) - a _2m constants defined by (30) - F _n () eigenfunctions, uniform wall heat flux - a _n , b _n , c _n , d _n , e _n constants defined by (45) and (48) - S a parameter, U ²/q - h ₁ heat transfer coefficient - T _m mean temperature - Nu Nusselt number - Nu _T Nusselt number, uniform wall temperature - Nu _q Nusselt number, uniform wall heat flux     

Hypersonic flow past blunt-edged delta wings

A method is proposed for calculating hypersonic ideal-gas flow past blunt-edged delta wings with aspect ratios = 100–200. Systematic wing flow calculations are carried out on the intervals 6 M 20, 0 20, 60 80; the results are analyzed in terms of hypersonic similarity parameters.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 175–179, September–October, 1990.     

Effect of Porous Coatings on Stability of Hypersonic Boundary Layers

The influence of ultrasoundabsorbing coatings on stability of hypersonic boundary layers is considered. Two types of coatings were used in experiments: feltmetal with a random porous microstructure and a sheet perforated by blind cylindrical microchannels. The experiments were performed in a wind tunnel at a Mach number M = 5.95 on sharp cones with a 7° apex halfangle. Evolution of natural disturbances and artificially induced wave packets in the boundary layer was studied with the help of hotwire anemometry. Spatial characteristics of artificial disturbances were obtained. It is demonstrated that such coatings exert a stabilizing effect on secondmode disturbances.     

A non-dualistic unified field theory of gravitation,electromagnetism and spin

The wisdom of classicalunified field theories in the conceptual framework of Weyl, Eddington, Einstein and Schrödinger has often been doubted and in particular there does not appear to be any empirical reason why the Einstein-Maxwell (E-M) theory needs to be geometrized. The crux of the matter is, however not whether the E-M theory is aesthetically satisfactory but whether it answers all the modern questions within the classical context. In particular, the E-M theory does not provide a classical platform from which the Dirac equation can be derived in the way Schrödinger's equation is derived from classical mechanics via the energy equation and the Correspondence Principle. The present paper presents a non-dualistic unified field theory (UFT) in the said conceptual framework as propounded by M. A. Tonnelat. By allowing the metric formds ²=g dx ^v x ^v and the non-degenerate two-formF=(1/2> _l) dx ^vdx ^vto enter symmetrically into the theory we obtain a UFT which contains Einstein's General Relativity and the Born-Infeld electrodynamics as special cases. Above all, it is shown that the Dirac equation describing the electron in an external gravito-electromagnetic field can be derived from the non-dualistic Einstein equation by a simple factorization if the Correspondence Principle is assumed.     

Heat Transfer on the Surface on a Spherically Blunted Cone Exposed to a Supersonic Spatial Flow with Injection of a Coolant Gas

The heattransfer processes in a supersonic spatial flow around a spherically blunted cone with allowance for heat overflow along the longitudinal and circumferential coordinates and injection of a coolant gas are studied numerically. The prospects of using highly heatconducting materials and injection of a coolant gas for reduction of the maximum temperatures at the body surface are demonstrated. The solutions of the direct and inverse problems in one, two, and threedimensional formulations for different shell materials are compared. The error of the thinwall method in determining the heat flux on the heatloaded boundary of the body is estimated.     

On the sensitive dynamical system and the transition from the apparently deterministic process to the completely random process

The detailed analysis of the dynamical process of coin tossing is made. Through calculations, it is illustrated how and why the result is extremely sensitive to the initial conditions. It is also shown that, as the initial height of the mass center of the coin increases, the final configuration, i.e. head or tail, becomes more and more sensitive to the initial parameters (the initial velocity angular velocity, and the initial orientation), the coefficient of the air drag, and the energy absorption factor of the surface on which the coin bounces. If we keep the head upward initially but allow a small range for the change of some other initial parameters, the frequency that the final configuration is head, would be 1 if the initial height h of the mass center is sufficiently small, and would be clo to 1/2 if h is sufficiently large. An interesting question is how this frequency changes continuously from 1 to 1/2 as h increases. Detailed calculations show that such a transition is very similar to the transition from laminar to turbulent flows. A basic difference between the transition stage and the completely random stage is indicated: In the completely random stage, the deterministic process of the individual case is extremely sensitive to the initial conditions and the dynamical parameters, out the statistical properties of the ensemble are insensitive to the small changes of the initial conditions and the dynamical parameters. On the contrary, in the transition stage, both the deterministic process of the individual case and the statistical properties of the ensemble are sensitive to the initial conditions and the dynamical parameters. The mechanism for this feature of the transition stage is the existence of the long-train structure in the parameter space. The illuminations of this analysis on some other random phenomena are discussed.     

Heat transport in a parabolic tube of circular cross-section

The local temperature has been determined for a viscous liquid flowing through a paraboloidal tube. Wall temperature and inlet temperature have been considered constant. The liquid flow was considered as creeping flow and its velocity distribution was determined by solving the biharmonic differential equation of the stream function. The local temperature was evaluated numerically from the analytical results.

---

Wärmetransport im Paraboloidrohr

Zusammenfassung Es wird die lokale Temperatur in einer viskosen Strömung durch ein Paraboloidrohr bestimmt. Dabei wird konstante Wand- und Einlauftemperatur angenommen. Die kriechende Strömungsgeschwindigkeit wurde aus der Lösung der biharmonischen Differentialgleichung der Stromfunktion bestimmt. Die lokale Temperatur wurde aus den analytischen Ergebnissen für einige Paraboloidrohre numerisch bestimmt.

---

Nomenclature ₁ F ₁ confluent hypergeometric function - diffusivity - T(, , ) temperature - T _w temperature at the paraboloidal wall - T _i temperature at the inlet - u(, ) flow velocity of viscous liquid in -direction - volumetric flow - eigenvalues of confluent hypergeometric function - streamfunction - _o wall of paraboloidal tube - _o inlet of paraboloidal tube - , , paraboloidal coordinates     

General expressions for unsteady forces in a lattice of profiles

A large number of studies have been devoted to the unsteady flow of a viscid incompressible fluid past a lattice of thin profiles and the determination of the resulting aerodynamic forces and moments. For example, in the particular case of the motion of a lattice with stagger with zero phase shift of the oscillations between neighboring profiles, Haskind [1] determined the unsteady lift force and moment. Popescu [2] suggested expressions for the force and moment in the case when =0 and =0, using the method of conformal mapping. Samoilovich [3] obtained equations for the unsteady lift force and moment by the method of the acceleration potential for phase shift =0 and = of the oscillations between neighboring profiles. Musatov [4] used an electronic digital computer to calculate the overall unsteady aerodynamic characteristics of a grid by the vortex method, taking into account the amplitude of the oscillations and the initial circulation for =m (m1). Gorelov [5] determined the coefficients of the over-all unsteady aerodynamic force and moment of a profile in a lattice with the stagger and any value of =m. He used a method based on the unsteady flow past an isolated profile with subsequent account for the interference of the profiles in the lattice.In the following we find general expressions for the unsteady lift force and moment acting on a lattice moving in an incompressible fluid with the constant velocity U. These formulas generalize the known formulas for the isolated profile [6]. The profiles of a staggered grid (Section 1) are considered to be thin and slightly curved, and perform oscillations with a phase shift of the oscillations between neighboring profiles. The method of separation of singularities is used to obtain the solution in closed form. The coefficients of the expansion of the complex velocity in a series in the derivatives of a function are calculated. An integral equation relative to the unknown tangential velocity component in the wake is derived (Section 2), and its analytic solution is given (Section 3). For =0 the solution coincides with the solution obtained earlier in [7]. Expressions are obtained for the forces and moments (Section 4) in the form of four terms. The first two terms determine the force and moment for motion with constant circulation, and the last two determine these characteristics for motion with variable circulation. The suction force acting at the leading edges of the profiles is found in a general form. Particular cases of closely and widely spaced lattices are considered. Computational results are presented.     

Stability of elastic bodies under omnilateral compression

Conclusions We have analyzed here the stability of the equilibrium of a simply connected isotropic compressible body with the elastic potential of arbitrary form and under uniform omnilateral deformation. A survey has been given here of earlier results obtained by other authors. The basic celations have been stated in a general form covering the theory of finite subcritical strains and two variants of the theory of small subcritical strains. For the latter theory new relations have been rigorously derived from which perturbations of tracking surface loads can be calculated, on the basis of corresponding expressions in the theory of finite subcritical strains. It has been proven that the sufficient conditions for the applicability of the static method of analysis are satisfied when the same boundary conditions are given over the entire body surface, as well as in several cases of different boundary conditions given at different segments of the boundary surface. It has been shown in a general form, for the theory of finite subcritical strains and for two variants of the theory of small subcritical strains, that the equilibrium of an elastic body under omnilateral deformation is stable, if a tracking load, is given over the entire boundary surface. As an example of problems with different boundary conditions at different segments of the boundary surface, we have considered the conventional problem concerning the stability of a bar on hinge supports and under uniform omnilateral deformation. It has been rigorously proven that in this case the equilibrium is stable when tracking loads are given at the lateral surfaces and is unstable when dead loads are given at the lateral surfaces. These conclusions apply to the theory of finite subcritical strains as well as to the theory of small subcritical strains, and they represent the complete version pertaining to compressible bodies.Institute of Mechanics, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Prikladnaya Mekhanika, Vol. 12, No. 6, pp. 3–27, June, 1976.     

An unsteady-state thermal boundary layer on a flat isothermal plate

The Kármán-Polhausen integral method is used to investigate the problem of an unsteady-state thermal boundary layer on an isothermal plate with a stepwise change in the conditions of flow around the plate; analytical expressions are obtained for the thickness of the thermal boundary layer. A dependence is found for the rate of movement of the boundary between the steady-state and unsteady-state regions of the solution on the Prandtl number. A similar problem was solved in [1, 2] for a dynamic layer, Goodman [3] discusses the more partial problem of an unsteady-state thermal boundary layer under steady-state flow conditions. Rozenshtok [4] considers the problem in an adequate statement but, unfortunately, he permitted errors of principle to enter into the writing of the system of characteristic equations; this led to absolutely invalid results. In an evaluation of the advantages and shortcomings of the integral method under consideration, given in [4], it must only be added that the method is applicable to problems in which the initial conditions differ from zero since, in this case, approximation of the velocity and temperature profiles by polynomials is not admissible.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 64–69, July–August, 1970.     

Electrode layers in flows of an inviscid plasma with good conductivity

The structure of the electromagnetic electrode layers that are produced in flows across a magnetic field by a completely ionized and inviscid plasma with good conductivity and a high magnetic Reynolds number is examined in a linear approximation. Flow past a corrugated wall and flow in a plane channel of slowly varying cross section with segmented electrodes are taken as specific examples. The possibility is demonstrated of the formation of nondissipative electrode layers with thicknesses on the order of the Debye distance or electron Larmor radius and of dissipative layers with thicknesses on the order of the skin thickness, as calculated from the diffusion rate in a magnetic field [2].In plasma flow in a transverse magnetic field, near the walls, along with the gasdynamie boundary layers, which owe their formation to viscosity, thermal conductivity, etc. (because of the presence of electromagnetic fields, their structures may vary considerably from that of ordinary gasdynamic layers), proper electromagnetic boundary layers may also be produced. An example of such layers is the Debye layer in which the quasi-neutrality of the plasma is upset. No less important, in a number of cases, is the quasi-neutral electromagnetic boundary layer, in which there is an abrupt change in the frozen-in parameter k=B/p (B is the magnetic field and p is the density of the medium). This layer plays a special role when we must explicitly allow for the Hall effect and the related formation of a longitudinal electric field (in the direction of the veloeiryv of the medium). We will call this the magnetic layer. The magnetic boundary layer can be dissipative as well as noudissipative (see below). The dissipative magnetic layer has been examined in a number of papers: for an incompressible medium with a given motion law in [1], for a compressible medium with good conductivity in [2], and with poor conductivity in [3]. In the present paper, particular attention will be devoted to nondissipative magnetic boundary layers.     

Flow over a body with development of a steady separation zone as re → ∞

This paper discusses formulation of the total problem of flow of an incompressible liquid over a body, with formation of a closed stationary separation zone as Re . The scheme used is based on the method of matched asymptotic expansions [1]. Following [1], it is postulated that the separated zone is developed (i.e., it is not infinitely fragmented and does not vanish as Re ), and the flow inside it has a definite degree of regularity with respect to Re. With these hypotheses we can use the Prandtl-Batchelor theorem [2], which states that, in the limit as Re , a region of circulating flow becomes vortex flow of an inviscid liquid with constant vorticity . Therefore, a basis for constructing matched asymptotic expansions is the vortex-potential problem (the problem of determining a stream function , satisfying the equation = 0 in the region of translational motion and the equation = in a certain region, unknowna priori, of circulating motion). In the general case the solution of the vortex-potential problem depends on two parameters: the total pressure p_o and the vorticity in the separated zone. These parameters appear in the condition for matching the solutions of the first and second boundary-layer approximations (at the boundary of the separated zone for the end Re values) with the corresponding solutions for the inviscid flow. It is shown in the present paper that the conditions for matching the cyclic boundary layer with the external translational flow are the same additional relations which allow us to close the total problem. Thus, in using the method of matched asymptotic expansions to solve the problem of flow over a body with closed stationary separation zones one must simultaneously consider no less than two approximations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 28–37, March–April, 1978.The authors thank G. Yu. Stepanov for discussion of the paper and valuable comments.     

Some problems of nonisothermal steady flow of a viscous fluid

The paper presents solutions to the problems of plane Couette flow, axial flow in an annulus between two infinite cylinders, and flow between two rotating cylinders. Taking into account energy dissipation and the temperature dependence of viscosity, as given by Reynolds's relation =₀ exp (–T) (₀, =const). Two types of boundary conditions are considered: a) the two surfaces are held at constant (but in general not equal) temperatures; b) one surface is held at a constant temperature, the other surface is insulated.Nonisothermal steady flow in simple conduits with dissipation of energy and temperature-dependent viscosity has been studied by several authors [1–11]. In most of these papers [1–6] viscosity was assumed to be a hyperbolic function of temperature, viz. =_m 1/1+²(T–T_m.Under this assumption the energy equation is linear in temperature and can he easily integrated. Couette flow with an exponential viscosity-temperature relation. =₀ e ^–T (₀, =const), (0.1) was studied in [7, 8]. Couette flow with a general (T) relation was studied in (9).Forced flow in a plane conduit and in a circular tube with a general (T) relation was studied in [10]. In particular, it has been shown in [10] that in the case of sufficiently strong dependence of viscosity on temperature there can exist a critical value of the pressure gradient, such that a steady flow is possible only for pressure gradients below this critical value.In a previous work [11] the authors studied Polseuille flow in a circular tube with an exponential (T) relation. This thermohydrodynamic problem was reduced to the problem of a thermal explosion in a cylindrical domain, which led to the existence of a critical regime. The critical conditions for the hydrodynamic thermal explosion and the temperature and velocity profiles were calculated.In this paper we treat the problems of Couette flow, pressureless axial flow in an annulus, and flow between two rotating cylinders taking into account dissipation and the variation of viscosity with temperature according to Reynolds's law (0.1). The treatment of the Couette flow problem differs from that given in [8] in that the constants of integration are found by elementary methods, whereas in [8] this step involved considerable difficulties. The solution to the two other problems is then based on the Couette problem.     

An application of non-singularity BEM to plate bending, buckling and natural vibration problems

In this paper the fundamental solution of the singular governing equation of plate static bending is taken as the Green's function, which can satisfy the governing equation precisely in the plate region. Based on the principle of superposition, let the function values on the plate boundary, induced by a set of the Green's function sources (including the known sources in the plate region and the unknown sources in the fictitious region), satisfy the prescribed conditions on specially chosen boundary matching points, and the corresponding semi-analytical and semi-numerical solution can be obtained, which is free from the restraint of boundary forms and boundary conditions. The more matching points there are on the boundary, the better the accuracy of results is. Finally, in static bending problems a set of linear algebraic equations has to be computed; in buckling problems the minimum value of buckling eigenvalue equation has to be found; in natural vibration problems the eigenvalues of the frequency equation have to be calculated. Numerical examples are given and the results are compared with those by the analytical method and other methods. It can be seen that they are very close to each other.     

Effect of a longitudinal magnetic field on the characteristics of a stabilized channel arc

The results of calculations of the temperature profiles and volt-ampere characteristics of a long cylindrical argon arc in a longitudinal uniform magnetic field are presented. The calculation was made for the following parameters: pressure p =0.1–10.0 atm; temperatures T = 1000-20,000°K; magnetic field induction B =0-10 T; diameter of cylindrical channel d = 1.0 cm. It is shown that for strongly radiating arcs (p1.0 atm) the temperature profiles become more inflated with an increase in the magnetic field, while for weakly radiating arcs (p 0.1 atm) the appearance of loops in the volt-ampere characteristics is typical for certain conditions (14,000T20,000°K, B1.0 T), indicating the impossibility of arcing under these conditions.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 147–153, March–April, 1975.