Generalized diagonalization of matrices over quaternion field

Ｉｎｒｅｃｅｎｔｙｅａｒｓ,ａｐｐｌｉｃａｔｉｏｎｓｏｆｑｕａｔｅｒｎｉｏｎｍａｔｒｉｃｅｓａｒｅｂｅｃｏｍｉｎｇｍｏｒｅａｎｄｍｏｒｅｉｍｐｏｒｔａｎｔａｎｄｅｘｔｅｎｓｉｖｅｉｎｒｉｇｉｄｍｅｃｈａｎｉｃｓ,ｑｕａｎｔｕｍｍｅｃｈａｎｉｃｓ,ｃｏｎｔｒｏｌｔｈｅｏｒｙａｎｄｈｅｌｉｃａｌｔｅｃｈｎｏｌｏｇｙ［１～３］．Ｗｉｔｈｔｈｅｒａｐｉｄｄｅｖｅｌｏｐｍｅｎｔｏｆｔｈｅａｂｏｖｅｄｉｓｃｉｐｌｉｎｅｓ,ｉｔｉｓｇｅｔｔｉｎｇｍｏｒｅａｎｄｍｏｒｅｎｅｃｅｓｓａｒｙｆｏｒｕｓｔｏｆｕｒｔｈ…     

The double velocity correlation function of homogeneous turbulence with constant mean velocity gradient

In this article, as the velocity gradient is taken as a constant value, we obtain the solutions of the equation of fluctuation velocity after Fourier transform. When the mean velocity gradient is small, they represent the picture of eddies, of which the homogeneous turbulence (both isotropic and nonisotropic) of the final period is composed. By using the eddies of these types at different times, we may compose the steady turbulent field with the constant velocity gradient and this field may represent the turbulent field in the central part of the channel flow or pipe flow approximately. Then we may obtain the double velocity correlation function of this turbulent field, which involves both longitudinal correlation coefficient and the transversal correlation coefficient . We compare these theoretical coefficients with the experimental data of these coefficients at initial period and final period of isotropic homogeneous turbulence. And then we obtain the relationship between the turbulent double velocity correlation coefficient and the mean velocity gradient. Finally, we get the expressions of the Reynolds stress and the eddy viscosity coefficient.     

The flow regime in the turbulent near wake of a hemisphere

The work presented is a wind tunnel study of the near wake region behind a hemisphere immersed in three different turbulent boundary layers. In particular, the effect of different boundary layer profiles on the generation and distribution of near wake vorticity and on the mean recirculation region is examined. Visualization of the flow around a hemisphere has been undertaken, using models in a water channel, in order to obtain qualitative information concerning the wake structure.List of symbols C _p pressure coefficient, - D diameter of hemisphere - n vortex shedding frequency - p pressure on model surface - p ₀ static pressure - Re Reynolds number, - St Strouhal number, - U, V, W local mean velocity components - mean freestream velocity inX direction - U ^* shear velocity, - u, v, w velocity fluctuations inX, Y andZ directions - X Cartesian coordinate in longitudinal direction - Y Cartesian coordinate in lateral direction - Z Cartesian coordinate in direction perpendicular to the wall - it* boundary layer displacement thickness, - diameter of model surface roughness - elevation angleI - O boundary layer momentum thickness, - _w wall shearing stress - dynamic viscosity of fluid - density of fluid - streamfunction - _x longitudinal component of vorticity, - _y lateral component of vorticity, - _z vertical component of vorticity, This paper was presented at the Ninth symposium on turbulence, University of Missouri-Rolla, October 1–3, 1984     

The existence of limit cycles for system\dot x=Q(x,y), Y=P(x)

In [1], by a transformation on the Lienard equation system

On the operator of symmetric differentiation on a compact Riemannian manifold with boundary

We state a particular case of one of the theorems which we shall prove. Let Ω be a bounded open set in ℝ ⁿ with smooth boundary and let σ=(σ _ij )be a symmetric second-order tensor with components σ _ij εH ^k(Ω) for some (positive or negative) integer k; H ^k are Sobolev spaces on Ω. Then we have for some u _i εH ^k +1(Ω),i=1,...,n, if and only if (if k<0, the integral is in fact a duality) for any symmetric tensor (ω with components and such that ). Some applications in the theory of elasticity are also given.     

Forced convective flow over a flat plate in non-Newtonian power law fluids

A new forced convection parameter

Transient response of a finite crack in an orthotropic strip under the in-plane shear impact

The transient response of a central crack in an orthotropic strip under the in-plane shear impact loading is studied by using the dual integral equation method proposed by Copson and Sih. The general formula for the shear stress intensity factor near the crack tip is derived. Numerical results of with in various cases are obtained by solving the second kind Fredholm integral equation and by performing the inverse Laplace transform.     

A class of two-level explicit difference schemes for solving three dimensional heat conduction equation

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｉｓｐａｐｅｒｄｅａｌｓｗｉｔｈｔｈｅｉｎｉｔｉａｌ＿ｂｏｕｎｄａｒｙｖａｌｕｅｐｒｏｂｌｅｍｏｆｔｈｒｅｅ＿ｄｉｍｅｎｓｉｏｎａｌｈｅａｔｃｏｎｄｕｃｔｉｏｎｅｑｕａｔｉｏｎｉｎｔｈｅｒｅｇｉｏｎＤ :0≤ｘ,ｙ ,ｚ≤Ｌ ,0 ≤ｔ≤Ｔ ｕ ｔ= 2 ｕ ｘ2 2 ｕ ｙ2 2 ｕ ｚ2 ,ｕ｜ｘ=0 =ｆ1(ｙ,ｚ,ｔ) ,　ｕ｜ｘ=Ｌ =ｆ2 (ｙ ,ｚ,ｔ) ,ｕ｜ｙ=0 =ｇ1(ｚ,ｘ,ｔ) ,　ｕ｜ｙ=Ｌ =ｇ2 (ｚ,ｘ,ｔ) ,ｕ｜ｚ=0 =ｈ1(ｘ ,ｙ ,ｔ) ,　ｕ｜ｚ=Ｌ =ｈ2 (ｘ ,ｙ ,ｔ) ,ｕ｜ｔ=0 =φ(ｘ ,ｙ,ｚ) .(1 )(2 )…     

A technique for evaluation of three-dimensional behavior in turbulent boundary layers using computer augmented hydrogen bubble-wire flow visualization

A system is described which allows the recreation of the three-dimensional motion and deformation of a single hydrogen bubble time-line in time and space. By digitally interfacing dualview video sequences of a bubble time-line with a computer-aided display system, the Lagrangian motion of the bubble-line can be displayed in any viewing perspective desired. The u and v velocity history of the bubble-line can be rapidly established and displayed for any spanwise location on the recreated pattern. The application of the system to the study of turbulent boundary layer structure in the near-wall region is demonstrated.List of Symbols Reynolds number based on momentum thickness u /v - t⁺ nondimensional time - u shear velocity - u local streamwise velocity, x-direction - u ⁺ nondimensional streamwise velocity - v local normal velocity, -direction - x ⁺ nondimensional coordinate in streamwise direction - ⁺ nondimensional coordinate normal to wall - ₊ ^wire wire nondimensional location of hydrogen bubble-wire normal to wall - z ⁺ nondimensional spanwise coordinate - momentum thickness - v kinematic viscosity - _W wall shear stress     

The extended Jordan's lemma and the relation between laplace transform and fourier transform

I.IntroductionItiswell4n0wnthatthecontourintergrationofcomp1exvariableftinctioncanmakealotintegrationveryconvenient.Jordan'slemmahasaveryimportantstatusintheonec0mplexvariableintegration,anditisveryusefulforavarityofintegration.Withthetheoryoffunctionsofo…     

Two-phase flow in heterogeneous porous media I: The influence of large spatial and temporal gradients

In order to capture the complexities of two-phase flow in heterogeneous porous media, we have used the method of large-scale averaging and spatially periodic models of the local heterogeneities. The analysis leads to the large-scale form of the momentum equations for the two immiscible fluids, a theoretical representation for the large-scale permeability tensor, and a dynamic, large-scale capillary pressure. The prediction of the permeability tensor and the dynamic capillary pressure requires the solution of a large-scale closure problem. In our initial study (Quintard and Whitaker, 1988), the solution to the closure problem was restricted to the quasi-steady condition and small spatial gradients. In this work, we have relaxed the constraint of small spatial gradients and developed a dynamic solution to the closure problem that takes into account some, but not all, of the transient effects that occur at the closure level. The analysis leads to continuity and momentum equations for the-phase that are given by

Asymptotic method for singular perturbation problem of ordinary difference equations

This paper is taken up for the following difference equation problem(P,)(L,y)_k≡εy(k 1) a(k,ε)y(k) b(k,ε)y(k-1)=f(k,ε)(1≤k≤N-1),B_1y≡-y(0) c_1y(1)=a,B_2y≡-c_2y(N-1) y(N)=βwhereεis a small parameter,c_1,c_2,a,βconstants and a(k,ε),b(k,ε),f(k,ε)(1≤k≤N)functions of k andε.Firstly,the case with constant coefficients isconsidered.Secondly,a general method based on extended transformation is given tohandle(P.)where the coefficients may be variable and uniform asymptotic expansionsare obtained Finally,a numerical example is provided to illustrate the proposed method.     

Numerical study of a nonlinear integrodifferential equation

I.IntroductionItiswell-knobal.nthatKorteweg-deVriesequationisacanonicalmodeltodescribethebalanceofthenonlineareffectandthedispersiveeffectofaphysicalsystem.Thisequationpossessestheso-called'soliton"solution,whichhasbeenfoundnumericallybyZabuskyandKruskall'].Ho-c'Jlever,sometimesthebalanceofnonlinearityanddispersionofasystemmayleadtoa,integroditTerentialequationinsteadofadifferentialequation.Forinstance,inthestudyofvortexbreakdownofanunboundedrotatingfluidLeibovich12]derivedfollowingnonline…     

On nonlinear hyperbolic equation in unbounded domain

The following nonlinear hyperbolic equation is discussed in this paper: where The model comes from the transverse deflection equation of an extensible beam. We prove that there exists a unique local solution of the above equation as M depends on x.     

The numerical stabilities of multiderivative block method

In [1], a class of multiderivative block methods (MDBM) was studied for the numerical solutions of stiff ordinary differential equations. This paper is aimed at solving the problem proposed in [1] that what conditions should be fulfilled for MDBMs in order to guarantee the A-stabilities. The explicit expressions of the polynomials and in the stability functions are given. Furthermore, we prove . With the aid of symbolic computations and the expressions of diagonal Pade' approximations, we obtained the biggest block size k of the A-stable MDBM for any given l (the order of the highest derivatives used in MDBM, l≥1)     

Plane stress dynamic fields near a propagating crack-tip in a power-law material

Based on the plastic-dynamic equations, the asymptotic behaviour of the near-tip fields for a plane stress tensile crack propagating in a power-law material has been studied in this paper. It is shown that the stress and strain singularities are, respectively, of the order and , whereA is a constant which is related to the size of plastic region,r is the distance to the crack tip,n is the power-law exponent. Projects sponsored by the National Science Foundation.     

The unsteady aerodynamics of a first stage stator vane row

The fundamental unsteady aerodynamics on a vane row of an axial flow research compressor stage are experimentally investigated, demonstrating the effects of airfoil camber and steady loading. In particular, the rotor wake generated unsteady surface pressure distributions on the first stage vane row are quantified over a range of operating conditions. These cambered airfoil unsteady data are correlated with predictions from a flat plate cascade inviscid flow model. At the design point, the unsteady pressure difference coefficient data exhibit good correlation with the nonseparated predictions, with the aerodynamic phase lag data exhibiting fair trendwise correlation. The quantitative phase lag differences are associated with the camber of the airfoil. An aft suction surface flow separation region is indicated by the steady state surface static pressure data as the aerodynamic loading is increased. This separation affects the increased incidence angle unsteady pressure data.List of symbols b airfoil semi-chord - C airfoil chord - C _p dynamic pressure coefficient, - _p static pressure coefficient, - i incidence angle - k reduced frequency, - N number of rotor revolutions - p dynamic pressure difference - static pressure difference, - S stator vane circumferential spacing - U _t rotor blade tip speed - u longitudinal perturbation velocity - V absolute velocity - V _axial absolute axial velocity - v transverse perturbation velocity - x _sep location of separation point - inlet angle - inlet air density - blade passing angular frequency     

Behaviour of the Reynolds stress on rough walls

Quadrant analysis is used to study the contributions, associated to the four quadrants of the (u,v) plane, to the production of the turbulent shear stress on rough walls. The measurements are described for a fully developed turbulent flow between two rough plates with varying the parameterλ _z(span/height ratio of roughness elements). The application of this technique indicates that the (Q ₂) events (ejections) and (Q ₄) events (sweeps) cause an intense production of the Reynolds stress —ρ as compared with the (Q ₁) and (Q ₃) quadrants. The (Q ₂) contribution to the Reynolds stress depends on the geometry factorλ _z. Variation of the parameterλ _zaffects the distributions of and forH⩽3, whereH is a particular threshold. Comparison with boundary layer flow shows that the region 0.2⩽y/h⩽0.7 is characterized by a Reynolds stress production, independent of the flow nature. The third moment of the longitudinal velocity fluctuations is found to be sensitive to the surface roughness.     

Optimal cupolas of uniform strength: Spherical M-shells and axisymmetric T-shells

Summary The first part of this paper is concerned with the optimal design of spherical cupolas obeying the von Mises yield condition. Five different load combinations, which all include selfweight, are investigated. The second part of the paper deals with the optimal quadratic meridional shape of cupolas obeying the Tresca yield condition, considering selfweight plus the weight of a non-carrying uniform cover. It is established that at long spans some non-spherical Tresca cupolas are much more economical than spherical ones.

---

Optimale Kuppeln gleicher Festigkeit: Kugelschalen und axialsymmetrische Schalen

Übersicht Im ersten Teil dieser Arbeit wird der optimale Entwurf sphärischer Kuppeln behandelt, wobei die von Misessche Fließbewegung zugrunde gelegt wird. Fünf verschiedene Lastkombinationen werden untersucht. Der zweite Teil befaßt sich mit der optimalen quadratischen Form des Meridians von Kuppeln, die der Fließbedingung von Tresca folgen.

---

List of Symbols a_k, b_k, c_k, A_k, B_k, C_k coefficients used in series solutions - A, B constants in the nondimensional equation of the meridional curve - normal component of the load per unit area of the middle surface - meridional and circumferential forces per unit width - radial pressure per unit area of the middle surface, - skin weight per unit area of the middle surface, - vertical external load per unit horizontal area, - base radius, - R radius of convergence - s - cupola thickness, - u, w subsidiary functions for quadratic cupolas - vertical component of the load per unit area of middle surface - resultant vertical force on a cupola segment - structural weight of cupola, - combined weight of cupola and skin, - distance from the axis of rotation, - vertical distance from the shell apex, - z auxiliary variable in series solutions - specific weight of structural material of cupola - radius of the middle surface, - uniaxial yield stress - meridional stress, - circumferential stress, - _a, _b, _c, _d, _e subsidiary variables used in evaluating the meridional stress - auxiliary function used in series solutions This paper constitutes the third part of a study of shell optimization which was initiated and planned by the late Prof. W. Prager     

Existence of Weak Solutions to the Equations of Non-Stationary Motion of Non-Newtonian Fluids with Shear Rate Dependent Viscosity

In the present paper we prove the existence of weak solutions to the equations of non-stationary motion of an incompressible fluid with shear rate dependent viscosity in a cylinder Q = Ω × (0,T), where denotes an open set. For the power-low model with we are able to construct a weak solution with ∇ · u = 0.