Displacement type unified equation for bending, buckling and vibration of conical shells and its applications

By using Donnell's simplication and starting from the displacement type equations of conical shells, and introducing a displacement functionU(s,,) (In the limit case, it will be reduced to cylindrical shell displacement function introduced by V. S. Vlasov) and a generalized loadq,(s,,),the equations of conical shells are changed into an eighth—order solvable partial differential equation about the displacement functionU(s,,). As a special case, the general bending problem of conical shells on Winkler foundation has been studied. Detailed numerical results and boundary coefficients for edge unit loads are obtained.The project supported by the National Natural Science Foundation of China.     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

Analysis of turbulent boundary layer interaction with an external supersonic flow behind a step

An integral method of analyzing turbulent flow behind plane and axisymmetric steps is proposed, which will permit calculation of the pressure distribution, the displacement thickness, the momentum-loss thickness, and the friction in the zone of boundary layer interaction with an external ideal flow. The characteristics of an incompressible turbulent equilibrium boundary layer are used to analyze the flow behind the step, and the parameters of the compressible boundary layer flow are connected with the parameters of the incompressible boundary layer flow by using the Cowles-Crocco transformation.A large number of theoretical and experimental papers devoted to this topic can be mentioned. Let us consider just two [1, 2], which are similar to the method proposed herein, wherein the parameter distribution of the flow of a plane nearby turbulent wake is analyzed. The flow behind the body in these papers is separated into a zone of isobaric flow and a zone of boundary layer interaction with an external ideal flow. The jet boundary layer in the interaction zone is analyzed by the method of integral relations.The flow behind plane and axisymmetric steps is analyzed on the basis of a scheme of boundary layer interaction with an external ideal supersonic stream. The results of the analysis by the method proposed are compared with known experimental data.Notation x, y longitudinal and transverse coordinates - X, Y transformed longitudinal and transverse coordinates - , ^*, ^** boundary layer thickness, displacement thickness, momentum-loss thickness of a boundary layer - , ^*, ^** layer thickness, displacement thickness, momentum-loss thickness of an incompressible boundary layer - u, velocity and density of a compressible boundary layer - U, velocity and density of the incompressible boundary layer - , stream function of the compressible and incompressible boundary layers - , dynamic coefficient of viscosity of the compressible and incompressible boundary layers - r₁ radius of the base part of an axisymmetric body - r radius - R transformed radius - M Mach number - friction stress - p pressure - a speed of sound - s enthalpy - v Prandtl-Mayer angle - P Prandtl number - P_t turbulent Prandtl number - r₂ radius of the base sting - b step depth - =0 for plane flow - =1 for axisymmetric flow Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 33–40, May–June, 1971.In conclusion, the authors are grateful to M. Ya. Yudelovich and E. N. Bondarev for useful comments and discussions.     

Evaluation of the dynamic-mechanical properties of solids from a cooperative model for stress relaxation

For many solid materials the stress relaxation process obeys the universal relationF = – (d/d lnt)_max = (0.1 ± 0.01) ( ₀ — _i ), regardless of the structure of the material. Here denotes the stress,t the time, ₀ the initial stress of the experiment and _i the internal stress. A cooperative model accounting for the similarity in relaxation behaviour between different materials was developed earlier. Since this model has a spectral character, the concepts of linear viscoelasticity are used here to evaluate the corresponding prediction of the dynamic mechanical properties, i.e. the frequency dependence of the storageE () and lossE () moduli. Useful numerical approximations ofE () andE () are also evaluated. It is noted that the universal relation in stress relaxation had a counterpart in the frequency dependence ofE (). The theoretical prediction of the loss factor for high-density polyethylene is compared with experimental results. The agreement is good.     

The secondary flow induced around a sphere in an oscillating stream of elastico-viscous liquid

The present paper is devoted to the theoretical study of the secondary flow induced around a sphere in an oscillating stream of an elastico-viscous liquid. The boundary layer equations are derived following Wang's method and solved by the method of successive approximations. The effect of elasticity of the liquid is to produce a reverse flow in the region close to the surface of the sphere and to shift the entire flow pattern towards the main flow. The resistance on the surface of the sphere and the steady secondary inflow increase with the elasticity of the liquid.Nomenclature a radius of the sphere - b ^ik contravariant components of a tensor - e contravariant components of the rate of strain tensor - F() see (47) - G total nondimensional resistance on the surface of the sphere - g _ik covariant components of the metric tensor - f, g, h secondary flow components introduced in (34) - k ₀ measure of relaxation time minus retardation time (elastico-viscous parameter) - K =k ₀ ²/V ₀ ² , nondimensional parameter characterizing the elasticity of the liquid - n measure of the ratio of the boundary layer thickness and the oscillation amplitude - N, T defined in (44) - p arbitrary isotropic pressure - p _ik covariant components of the stress tensor - p ^ik contravariant components of the stress tensor associated with the change of shape of the material - R =V ₀ a/v, the Reynolds number - S =a/V ₀, the Strouhall number - r, , spherical polar coordinates - u, v, w r, , component of velocity - t time - V(, t) potential velocity distribution around the sphere - V ₀ characteristic velocity - u, v, t, y, P nondimensional quantities defined in (15) - reciprocal of s - density - defined in (32) - defined in (42) - ₀ limiting viscosity for very small changes in deformation velocity - complex conjugate of - oscillation frequency - = ₀/, the kinematic coefficient of viscosity - , defined in (52) - (, y) stream function defined in (45) - =(NT/2n)^1/2 y - /t convective time derivative (1) _ik      

The value of the critical exponent for reaction-diffusion equations in cones

Let D R ^N be a cone with vertex at the origin i.e., D = (0, )x where S ^N–1 and x D if and only if x = (r, ) with r=¦x¦, . We consider the initial boundary value problem: u _t = u+u ^p in D×(0, T), u=0 on Dx(0, T) with u(x, 0)=u ₀(x) 0. Let ₁ denote the smallest Dirichlet eigenvalue for the Laplace-Beltrami operator on and let ₊ denote the positive root of (+N–2) = ₁. Let p ^* = 1 + 2/(N + ₊). If 1 < p < p ^*, no positive global solution exists. If p>p ^*, positive global solutions do exist. Extensions are given to the same problem for u _t=+¦x¦ u ^p .This research was supported in part by the Air Force Office of Scientific Research under Grant # AFOSR 88-0031 and in part by NSF Grant DMS-8 822 788. The United States Government is authorized to reproduce and distribute reprints for governmental purposes not withstanding any copyright notation therein.     

On the compactness of quasi-conforming element spaces and the convergence of quasi-conforming element method

In this paper, the compactness of quasi-conforming element spaces and the—convergence of quasi-conforming element method are discussed. The well-known Rellich compactness theorem is generalized to the sequences of quasi-conforming element spaces with certain properties, and the generalized Poincare inequality. The generalized Friedrichs inequality and the generalized inequality of Poincare-Friedrichs are proved true for them. The error estimates are also given. It is shown that the quasi-conforming element method is convergent if the quasi-conforming element spaces have the approximability and the strong continuity, and satisfy the rank condition of element and pass the test IPT. As practical examples, 6-parameter, 9-paramenter, 12-paramenter, 15-parameter, 18-parameter and 21-paramenter quasi-conforming elements are shown to be convergent, and their L²²()-errors are O(h), O(h), O(h ² ), O(h ² ), O(h ), and O(h ⁴ ) respectively.     

Two-dimensional transonic vortex flows of an ideal gas

Equations are obtained for two-dimensional transonic adiabatic (nonisoenergetic and nonisoentropic) vortex flows of an ideal gas, using the natural coordinates (=const is the family of streamlines, and =const is the family of lines orthogonal to them). It is not required that the transonic gas flow be close to a uniform sonic flow (the derivation is given without estimates). Solutions are found for equations describing vortex flows inside a Laval nozzle and near the sonic boundary of a free stream.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 105–109, September–October, 1973.     

Hydrodynamic models as applied to the investigation of magnetic plasma stability

In the present paper magnetohydrodynamic models are employed to investigate the stability of an inhomogeneous magnetic plasma with respect to perturbations in which the electric field may be regarded as a potential field (rot E 0). A hydrodynamic model, actually an extension of the well-known Chew-Goldberg er-Low model [1], is used to investigate motions transverse to a strong magnetic field in a collisionless plasma. The total viscous stress tensor is given; this includes, together with magnetic viscosity, the so-called inertial viscosity.Ordinary two-fluid hydrodynamics is used in the case of strong collisions=. It is shown that the collisional viscosity leads to flute-type instability in the case when, collisions being neglected, the flute mode is stabilized by a finite Larmor radius. A treatment is also given of the case when epithermal high-frequency oscillations (not leading immediately to anomalous diffusion) cause instability in the low-frequency (drift) oscillations in a manner similar to the collisional electron viscosity, leading to anomalous diffusion.Notation f particle distribution function - E electric field component - H₀ magnetic field - density - V particle velocity - e charge - m, M electron and ion mass - _i, _e ion and electron cyclotron frequencies - viscous stress tensor - P pressure - r_i Larmor radius - P pressure tensor - t time - frequency - T temperature - collision frequency - collision time - j current density - _i, _e ion and electron drift frequencies - k_x, k_y, k_z wave-vector components - n₀ particle density - g acceleration due to gravity. The authors are grateful to A. A. Galeev for valuable discussion.     

Analysis of suddenly started laminar flow in the entrance region of a circular tube

Suddenly started laminar flow in the entrance region of a circular tube, with constant inlet velocity, is investigated analytically by using integral momentum approach. A closed form solution to the integral momentum equation is obtained by the method of characteristics to determine boundary layer thickness, entrance length, velocity profile, and pressure gradient.Nomenclature M(, , ) a function - N(, , ) a function - p pressure - p* p/1/2U ², dimensionless pressure - Q(, , ) a function - R radius of the tube - r radial distance - Re 2RU/, Reynolds number - t time - U inlet velocity, constant for all time, uniform over the cross section - u velocity in the boundary layer - u* u/U, dimensionless velocity - u ₁ velocity in the inviscid core - x axial distance - y distance perpendicular to the axis of the tube - y* y/R, dimensionless distance perpendicular to the axis - boundary layer thickness - * displacement thickness - /R, dimensionless boundary layer thickness - momentum thickness - absolute viscosity of the fluid - /, kinematic viscosity of the fluid - x/(R Re), dimensionless axial distance - density of the fluid - tU/(R Re), dimensionless time - _w wall shear stress     

Heteroclinic orbits in singular systems: A unifying approach

We consider singularly perturbed systems , such that=f(, _o, 0). _o ^m , has a heteroclinic orbitu(t). We construct a bifurcation functionG(, ) such that the singular system has a heteroclinic orbit if and only ifG(, )=0 has a solution=(). We also apply this result to recover some theorems that have been proved using different approaches.     

A note on the tan δ(T) peak as a glass transition indicator in biosolids

The plasticization of many biosolids can take place over a fairly broad temperature range. The resulting loss of stiffness is primarily expressed by a drastic drop of G(T) whose magnitude is usually higher than G(T) by one or two orders of magnitude. Both G(T) and G(T) have characteristic properties that can vary widely among biomaterials. Consequently, the tan (T) peak need not be a mark of the transition center and it can be observed at temperatures where different materials have undergone a very different degree of plasticization as judged by the magnitude of G(T). This is demonstrated by computer simulations using typical functions that describe G(T) and G(T) at the glass transition region and with published data on the dynamic mechanical behavior of a variety of biosolids.     

Pulsatile blood flow through arterioles

An analytical study was made to examine the effect of vascular deformability on the pulsatile blood flow in arterioles through the use of a suitable mathematical model. The blood in arterioles is assumed to consist of two layers — both Newtonian but with differing coefficients of viscosity. The flow characteristics of blood as well as the resistance to flow have been determined using the numerical computations of the resulting expressions. The applicability of the model is illustrated using numerical results based on the existing experimental data. r, z coordinate system - u, axial/longitudinal velocity component of blood - p pressure exerted by blood - _b density of blood - µ viscosity of blood - t time - , displacement components of the vessel wall - T _t0,T ₀ known initial stresses - density of the wall material - h thickness of the vessel wall - T _t,T stress components of the vessel - K _l,K _r components of the spring coefficient - C _l,C _r components of the friction coefficient - M _a additional mass of the mechanical model - r ₁ outer radius of the vessel - thickness of the plasma layer - r ₁ inner radius of the vessel - circular frequency of the forced oscillation - k wave number - E ₀,E _t, , _t material parameters for the arterial segment - µ _p viscosity of the plasma layer - Q total flux - Q _p flux across the plasma zone - Q _h flux across the core region - Q mean flow rate - resistance to flow - P pressure difference - l length of the segment of the vessel     

Measurements of thermally stratified pipe flow using image-processing techniques

The cross-correlation technique and Laser Induced Fluorescence (LIF) have been adopted to measure the time-dependent and two-dimensional velocity and temperature fields of a stably thermal-stratified pipe flow. One thousand instantaneous and simultaneous velocity and temperature maps were obtained at overall Richardson numberRi = 0 and 2.5, from which two-dimensional vorticity, Reynolds stress and turbulent heat flux vector were evaluated. The quasi-periodic inclined vortices (which connected to the crest) were revealed from successive instantaneous maps and temporal variation of vorticity and temperature. It has been recognized that these vortices are associated with the crest and valley in the roll-up motion.List of symbols A Fraction of the available light collected - C Concentration of fluorescence - D Pipe diameter - I Fluorescence intensity - L Sampling length along the incident beam - I ₀ Intensity of an excitation beam - I _c (T) Calibration curve between temperature and fluorescence intensity - I _ref Reference intensity of fluorescence radiation - Re _b Reynolds number based on bulk velocity,U _b D/v - Ri Overall Richardson number based on velocity difference,gDT/U ² - t Time - t Time interval between the reference and corresponding matrix - T Temperature - T ₁,T ₂ Temperature of lower and upper layer - T ^* Normalized temperature, (T–T ₁)/T - T _c (I) Inverse function of temperature as a function ofI _c - T _ref Reference temperature - T Temperature difference between upper and lower flow,T ₂ –T ₁ - U ₁ Velocity of lower stream - U ₂ Velocity of upper stream - U _b Bulk velocity - U _c Streamwise mean velocity atY/D=0 - U Streamwise velocity difference between upper and lower flow,U ₁–U ₂ - u, v, T Fluctuating component ofU, V, T - U, V Velocity component of X, Y direction - X Streamwise distance from the splitter plate - Y Transverse distance from the centerline of the pipe - Z Spanwise distance from the centerline of the pipe - Quantum yield - Absorptivity - vorticity calculated from a circulation - Kinematic viscosity - circulation     

On nonlinear problems of mixed type: A qualitative theory using infinite-dimensional center manifolds

We consider the equation a(y)u_xx+div_y(b(y)_yu)+c(y)u=g(y, u) in the cylinder (–l,l)×, being elliptic where b(y)>0 and hyperbolic where b(y)<0. We construct self-adjoint realizations in L₂() of the operatorAu= (1/a) div_y(b_yu)+(c/a) in the case ofb changing sign. This leads to the abstract problem u_xx+Au=g(u), whereA has a spectrum extending to + as well as to –. For l= it is shown that all sufficiently small solutions lie on an infinite-dimensional center manifold and behave like those of a hyperbolic problem. Anx-independent cross-sectional integral E=E(u, u_x) is derived showing that all solutions on the center manifold remain bounded forx ±. For finitel, all small solutionsu are close to a solution on the center manifold such that u(x)-(x) Ce ^-(1-|x|) for allx, whereC and are independent ofu. Hence, the solutions are dominated by hyperbolic properties, except close to the terminal ends {±1}×, where boundary layers of elliptic type appear.     

A Direct Method for the Identification of the Permeability Field Based on Flux Assimilation by a Discrete Analog of Darcy's Law

Inverse models to determine the permeability are generally based on existing forward models for the pressure. The permeabilities are adapted in such a way that the calculated pressures match the specified pressures in a number of points. To assimilate a priori knowledge about the flux, we introduce the flux assimilation method, which is based on the vector potential–pressure formulation of Darcy's law. Thanks to an unconventional discretization technique – the edge-based face element method – not only the specified pressures, but also specified information about the flux density can easily be assimilated. A relatively simple, but insightful analytical example illustrates the potential of this method.     

Asymptotic theory of the turbulent boundary layer as a problem of singular perturbations

We consider an asymptotic theory of the turbulent boundary layer [1,2]. In this paper we make an attempt to further develop the mathematical aspects of this theory. We demonstrate the features of this theory by applying it to a problem which is close to the so-called equilibrium turbulent boundary layer with a pressure gradient and blowing.Notation x, y coordinates, parallel and perpendicular to the wall - u velocity component in the x direction - p, ',v pressure, density, and kinematic viscosity coefficient - l' scale of turbulence - tangential stress - u speed at the outer edge of the boundary layer - thickness of the boundary layer - * displacement thickness - ** momentum loss thickness - Cf coefficient of friction - R Reynolds number Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 86–95, September–October, 1971.The author thanks S. S. Kutateladze, A. I. Leont'ev, and G. V. Aronovich for their interest in this effort.     

On the Discretization of Distributed-Parameter Systems with Quadratic and Cubic Nonlinearities

Approximate methods for analyzing the vibrations of an Euler--Bernoulli beam resting on a nonlinear elastic foundation are discussed. The cases of primary resonance ( _n ) and subharmonic resonance of order one-half ( 2 _n ), where is the excitation frequency and _n is the natural frequency of the nth mode of the beam, are investigated. Approximate solutions based on discretization via the Galerkin method are contrasted with direct application of the method of multiple scales to the governing partial-differential equation and boundary conditions. The amplitude and phase modulation equations show that single-mode discretization leads to erroneous qualitative as well as quantitative predictions. Regions of softening (hardening) behavior of the system, the spatial dependence of the response drift, and frequency-response curves are numerically evaluated and compared using both approaches.     

On the effects of uniform high suction on the rotationally symmetric flow of a conducting liquid near a stationary disc in the presence of a transverse magnetic field

In the present paper an attempt has been made to find out effects of uniform high suction in the presence of a transverse magnetic field, on the motion near a stationary plate when the fluid at a large distance above it rotates with a constant angular velocity. Series solutions for velocity components, displacement thickness and momentum thickness are obtained in the descending powers of the suction parameter a. The solutions obtained are valid for small values of the non-dimensional magnetic parameter m (= 4 _e ² H ₀ ² /) and large values of a (a2).Nomenclature a suction parameter - E electric field - E _r , E , E _z radial, azimuthal and axial components of electric field - F, G, H reduced radial, azimuthal and axial velocity components - H magnetic field - H _r , H , H _z radial, azimuthal and axial components of magnetic field - H ₀ uniform magnetic field - H* displacement thickness and momentum thickness ratio, */ - h induced magnetic field - h _r , h , h _z radial, azimuthal and axial components of induced magnetic field - J current density - m nondimensional magnetic parameter - p pressure - P reduced pressure - R Reynolds number - U ₀ representative velocity - V velocity - V _r , V , V _z radial, azimuthal and axial velocity components - w ₀ uniform suction through the disc. - density - electrical conductivity - kinematic viscosity - _e magnetic permeability - a parameter, (/)^1/2 z - a parameter, a - * displacement thickness - momentum thickness - angular velocity     

Sensitivity analysis of the non-linear dynamic response of beams using space-time perturbation modes

The focus of the present work is directed towards the development of an effective reduced basis technique for calculating the sensitivity of the non-linear dynamic structural response of mechanical systems with respect to variations in the design variables.The proposed methodology is formulated within the context of a mixed space-time finite element method, which naturally allows the treatment of initial and boundary value problems. The time dependency of the solutions is implied in the assumed space-time modal shapes, and hence the partial differential equations of motion are directly reduced to a set of non-linear simultaneous equations of a purely algebraic nature.The independent field variables are approximated in terms of perturbations modes or path derivatives with respect to a load control parameter. These modes, extracting information about the kinematic and dynamic behavior of the structural system through the higher order derivatives of the strain and kinetic energies, are appropriate bases for non-linear dynamic problems. The sensitivity derivatives of the field variables are then approximated using a combination of perturbation modes and of their sensitivity derivatives.The resulting computational procedure offers high potential for the effective and numerically efficient sensitivity analysis of dynamic systems exhibiting periodic-in-time response. The proposed methodology is illustrated addressing non-linear beam problems subjected to harmonic loading and the results obtained are compared with those of a full finite-element model.Nomenclature (O, I _i), (is1, 2, 3) Inertial frame of origin O - (P, s _i), (is1, 2, 3) Local frame in the undeformed configuration - (Q, s _i ^*), (is1, 2, 3) Local frame in the deformed configuration - t Time - l Abscissa along the beam reference line - L Beam length - (·)s(·)/t Partial derivative with respect to time - (·)s(·)/l Partial derivative with respect to space - u Position vector of the beam reference line - r Rotation parameters - ds(u, r) Generalized displacement vector - R(r) Rotation tensor associated with r - (r) Tensor defined in equations (4) and (5) - s· Finite rotation vector - as(a _s, a _v) Conformal rotation vector - Angular velocity - ws(u, ) Generalized velocity vector - k Curvature - e Generalized strains - ps(h, l) Generalized momenta - fs(s, m) Generalized sectional stress resultants - f _es(S _e, m _e) Applied external loads - M Inertia tensor