Scheme of Investigation of the Spectrum of a Family of Perturbed Operators and Its Application to Spectral Problems in Thick Junctions

We develop a scheme for the investigation of the asymptotic behavior of eigenvalues and eigenvectors of a family of self-adjoint compact operators {A: > 0} that act in different spaces and lose their compactness in the limit case 0. We prove the Hausdorff convergence of the spectrum of the operator A to the spectrum of the limit operator A₀, obtain asymptotic estimates for this convergence both to points of the discrete spectrum and to points of the essential spectrum of the operator A₀, and prove asymptotic estimates for eigenvectors of A. This scheme is applied to the investigation of the asymptotic behavior of eigenvalues and eigenfunctions of the Neumann problem in a thick singularly degenerate junction that consists of two domains connected by an -periodic system of thin rods of fixed length.     

Transcendentally small transversality in the rapidly forced pendulum

The rapidly forced pendulum equation with forcing sin((t/), where =<₀^p,p = 5, for ₀, sufficiently small, is considered. We prove that stable and unstable manifolds split and that the splitting distanced(t) in the ( ,t) plane satisfiesd(t) = sin(t/) sech(/2) +O( ₀ exp(–/2)) (2.3a) and the angle of transversal intersection,, in thet = 0 section satisfies 2 tan/2 = 2S _s = (/2) sech(/2) +O(( ₀ /) exp(–/2)) (2.3b) It follows that the Melnikov term correctly predicts the exponentially small splitting and angle of transversality. Our method improves a previous result of Holmes, Marsden, and Scheuerle. Our proof is elementary and self-contained, includes a stable manifold theorem, and emphasizes the phase space geometry.     

Propagation of a spherical wave in nonlinearly compressible and elastoplastic materials

The problem of spherical wave propagation in soil under the action of an intense uniformly decreasing load ₀(t) applied to the boundary of a cavity with radius r₀ is considered. Soil with a high stress level is modeled either by ideally nonlinearly compressible or elastoplastic material, taking account of linear irreversible unloading for the material. In contrast to [1–7], in order to describe material movement use is made of strain theory [8] with determining functions = (), _i=_i(_i), where , _i, , _i are the first and second invariants of strain and stress tensors. During material loading these functions are presented in the form of polynomials ()=(_i+₂¦¦), _ii)=(_i-_2i)_i, in which constant coefficients _i, _i=1, 2) are determined by experiment, taking account of the triaxial stressed state of soil. Solution of the problem is constructed by an analytically reversible method, with prescribed shape for the shock-wave (SW) surface in the form of a second-degree polynomial relating to time t and a numerical method of characteristics for a prescribed arbitrarily decreasing load _i(t). On the basis of the analytical equations obtained, calculations are carried out for material parameters (including loading profile) in a computer and stresses and mass velocity of plastic and elastoplastic materials are compared.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 95–100, July–August, 1986.The authors express their sincere thanks to Kh. A. Rakhmatulin for discussing the results of this work.     

On laminar flow through a uniformly porous pipe

Numerous investigations ([1] and [4–9]) have been made of laminar flow in a uniformly porous circular pipe with constant suction or injection applied at the wall. The object of this paper is to give a complete analysis of the numerical and theoretical solutions of this problem. It is shown that two solutions exist for all values of injection as well as the dual solutions for suction which had been noted by previous investigators. Analytical solutions are derived for large suction and injection; for large suction a viscous layer occurs at the wall while for large injection one solution has a viscous layer at the centre of the channel and the other has no viscous layer anywhere. Approximate analytic solutions are also given for small values of suction and injection.

Nomenclature

General r distance measured radially - z distance measured along axis of pipe - u velocity component in direction of z increasing - v velocity component in direction of r increasing - p pressure - density - coefficient of kinematic viscosity - a radius of pipe - V velocity of suction at the wall - r ²/a ² - R wall or suction Reynolds number, Va/ - f() similarity function defined in (6) - u ₀() eigensolution - U(0) a velocity at z=0 - K an arbitrary constant - B _K Bernoulli numbers Particular Section 5 perturbation parameter, –2/R - ² a constant, –K - x / - g(x) f()/ Section 6 perturbation parameter, –R/2 - ² a constant, –K - g() f() - g _c ()=g() near centre of pipe - * point where g()=0 Section 7 2/R - ² K - t (1–)/ - w(t, ) [1–f(t)]/ - ₀, ₁ constants - g() f()– ₀ - ₀/ - ₀ a constant - * point where f()=0     

Global Behavior of a Nonlinear Quasiperiodic Mathieu Equation

In this paper, we investigate the interaction of subharmonicresonances in the nonlinear quasiperiodic Mathieu equation,x + [ + (cos ₁ t + cos ₂ t)] x + x³ = 0.We assume that 1 and that the coefficient of the nonlinearterm, , is positive but not necessarily small.We utilize Lie transform perturbation theory with elliptic functions –rather than the usual trigonometric functions – to study subharmonic resonances associated with orbits in 2m:1 resonance with a respective driver. In particular, we derive analytic expressions that place conditions on (, , ₁, ₂) at which subharmonic resonance bands in a Poincaré section of action space begin to overlap. These results are used in combination with Chirikov's overlap criterion to obtain an overview of the O() global behavior of equation (1) as a function of and ₂ with ₁, , and fixed.     

A Turbulence Model Sensitivity Study for CH4/H2 Bluff-Body Stabilized Flames

The objective of this work is to assess the performances of different turbulence models in predicting turbulent diffusion flames in conjunction with the flamelet model.The k– model, the Explicit Algebraic Stress Model (EASM) and the k– model withvaried anisotropy parameter C (LEA k– model)are first applied to the inert turbulent flow over a backward-facing step, demonstrating the quality of the turbulence models. Following this, theyare used to simulate the CH₄/H₂ bluff-body flame studied by the University of Sydney/Sandia.The numerical results are compared to experimental values of the mixture fraction, velocity field, temperature and constituent mass fractions.The comparisons show that the overall result depends on the turbulence model used, and indicate that theEASM and the LEA k– models perform better than the k– model and mimic most of the significant flow features.     

Homogenization of the Robin problem in a thick multilevel junction

We consider a mixed boundary-value problem for the Poisson equation in a plane two-level junction that is the union of a domain ₀ and a large number 2N of thin rods with variable thickness of order = (N ^–1). The thin rods are divided into two levels, depending on their length. In addition, the thin rods from each level are -periodically alternated. We investigate the asymptotic behavior of the solution as 0 under the Robin conditions on the boundaries of the thin rods. By using some special extension operators, a convergence theorem is proved.Published in Neliniini Kolyvannya, Vol. 7, No. 3, pp. 336–355, July–September, 2004.     

Application of the Chapman-Enskog method to the case of a binary two-temperature gas mixture

An interesting property of the flows of a binary mixture of neutral gases for which the molecular mass ratio =m/M1 is that within the limits of the applicability of continuum mechanics the components of the mixture may have different temperatures. The process of establishing the Maxwellian equilibrium state in such a mixture divides into several stages, which are characterized by relaxation times _i which differ in order of magnitude. First the state of the light component reaches equilibrium, then the heavy component, after which equilibrium between the components is established [1]. In the simplest case the relaxation times differ from one another by a factor of ^*.Here the mixture component temperature difference relaxation time _T /, where is the relaxation time for the light component. If 1, 1, so that _T ~1, then for the characteristic hydrodynamic time scale t~1 the relative temperature difference will be of order unity. In the absence of strong external force fields the component velocity difference is negligibly small, since its relaxation time _vt1.In the case of a fully ionized plasma the Chapman-Enskog method is quite easily extended to the case of the two-temperature mixture [3], since the Landau collision integral is used, which decomposes directly with respect to . In the Boltzmann cross collision integral, the quantity appears in the formulas relating the velocities before and after collision, which hinders the decomposition of this integral with respect to , which is necessary for calculating the relaxation terms in the equations for temperatures differing from zero in the Euler approximation [4] (the transport coefficients are calculated considerably more simply, since for their determination it is sufficient to account for only the first (Lorentzian [5]) terms of the decomposition of the cross collision integrals with respect to ). This led to the use in [4] for obtaining the equations of the considered continuum mixture of a specially constructed model kinetic equation (of the Bhatnagar-Krook type) which has an undetermined degree of accuracy.In the following we use the Boltzmann equations to obtain the equations of motion of a two-temperature binary gas mixture in an approximation analogous to that of Navier-Stokes (for convenience we shall term this approximation the Navier-Stokes approximation) to determine the transport coefficients and the relaxation terms of the equations for the temperatures. The equations in the Burnett approximation, and so on, may be obtained similarly, although this derivation is not useful in practice.     

The influence of blending polystyrenes of narrow molecular weight distribution on melt creep flow and creep recovery in elongation

Creep and creep recovery experiments in elongation were performed with melts of anionically polymerized polystyrenes (PS) and with their blends at a temperature of 150 °C. For stresses ₀ < 10 000 N/m² the samples with narrow molecular weight distribution show linear viscoelastic behavior up to the maximum Hencky strain = 3.5, achievable in a newly developed elongational rheometer for polymer melts. The compliances,D (t), of the blends are linear-viscoelastic only up to a strain limit _L . For strains beyond _L the compliance of each blend depends on the stress ₀. For a series of binary blends, prepared from the same components of narrow MWD, the linear-viscoelastic limit _L seems to be independent of the mixing ratio and stress. _L seems to be a function only of the molecular weights of the original components, the blends investigated were made from.Paper presented at the Annual Conference of the German Society of Rheology at Ulm, March 7–10, 1983.     

Analytical evaluation of stress-strain test data

An analytical model for deducing the actual stress-strain properties from laboratory test results is discussed. As an illustration, an elastic bilinear material is used for unconfined cylindrical compression test conditions, as simulated with a finite element analysis. The results obtained are applicable for assisting in evaluating measured strength and stiffness properties of some clay soils, concrete test cylinders, concrete cores, and rock cores.The quantitative results of this study can be used for interpreting measured stress-strain data for unconfined compression test conditions. The error in measured results is shown to be influenced by Poisson's ratio, length-to-diameter ratio of the specimen, end condition, and ratio of inelastic modulus to initial elastic modulus. Curves for adjusting the measured results to the theoretical results are presented.Nomenclature D specimen diameter - E _i initial elastic stiffness modulus - E _y elastic stiffness modulus beyond the yield stress, plastic or inelastic modulus - L specimen length - axial strain - _av average strain - _g gage length strain - _y yield strain - Poisson's ratio - compressive stress - _av average stress - _t theoretical compressive stress - _y yield stress - _ym measured stress at the yield strain     

Dielectric properties of heterogeneous mixtures with a polar constituent

Summary After defining the boundaries for the dielectric constant of a heterogeneous mixture, the behaviour of such a mixture is studied as a function of the frequency, when one of its components is polar. Deviations from a semicircle are to be expected for the function _m =f( _m ) even when the dielectric properties of the polar constituent can be described with a semicircular Cole-Cole-arc. The relaxation time of the mixture is shorter than that of the polar constituent.     

Thermodynamically complete equations of state for α-, ε-, and γ-iron

The numerical model of phase transition in iron in stress waves described in [1] contains equations of state with a limited range of applicability. They do not consider thermal excitation of conduction electrons and the presence of and — -triple point on the phase equilibrium curve, the effect of which should appear in shock loading of porous or preheated specimens. The present study will offer thermodynamically complete equations of state for the -, -, -phases of iron, free of these shortcomings.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 112–114, May–June, 1986.     

The extension of Beltrami's ellipsoid to anisotropic bodies

Summary The spectral decomposition of the compliance, stiffness, and failure tensors for transversely isotropic materials was studied and their characteristic values were calculated using the components of these fourth-rank tensors in a Cartesian frame defining the principal material directions. The spectrally decomposed compliance and stiffness or failure tensors for a transversely isotropic body (fiber-reinforced composite), and the eigenvalues derived from them define in a simple and efficient way the respective elastic eigenstates of the loading of the material. It has been shown that, for the general orthotropic or transversely isotropic body, these eigenstates consist of two double components, ₁ and ₂ which are shears (₂ being a simple shear and ₁, a superposition of simple and pure shears), and that they are associated with distortional components of energy. The remaining two eigenstates, with stress components ₃, and ₄, are the orthogonal supplements to the shear subspace of ₁ and ₂ and consist of an equilateral stress in the plane of isotropy, on which is superimposed a prescribed tension or compression along the symmetry axis of the material. The relationship between these superimposed loading modes is governed by another eigenquantity, the eigenangle .The spectral type of decomposition of the elastic stiffness or compliance tensors in elementary fourth-rank tensors thus serves as a means for the energy-orthogonal decomposition of the energy function. The advantage of this type of decomposition is that the elementary idempotent tensors to which the fourth-rank tensors are decomposed have the interesting property of defining energy-orthogonal stress states. That is, the stress-idempotent tensors are mutually orthogonal and at the same time collinear with their respective strain tensors, and therefore correspond to energy-orthogonal stress states, which are therefore independent of each other. Since the failure tensor is the limiting case for the respective _x, which are eigenstates of the compliance tensor S, this tensor also possesses the same remarkable property.An interesting geometric interpretation arises for the energy-orthogonal stress states if we consider the projections of _x in the principal₃D stress space. Then, the characteristic state ₂ vanishes, whereas stress states ₁, ₃ and ₄ are represented by three mutually orthogonal vectors, oriented as follows: The ₃ and ₄ lie on the principal diagonal plane (₃₁₂) with subtending angles equaling (–/2) and (-), respectively. On the positive principal ₃-axis, is the eigenangle of the orthotropic material, whereas the ₁-vector is normal to the (₃₁₂)-plane and lies on the deviatoric -plane. Vector ₂ is equal to zero.It was additionally conclusively proved that the four eigenvalues of the compliance, stiffness, and failure tensors for a transversely isotropic body, together with value of the eigenangle , constitute the five necessary and simplest parameters with which invariantly to describe either the elastic or the failure behavior of the body. The expressions for the _x-vector thus established represent an ellipsoid centered at the origin of the Cartesian frame, whose principal axes are the directions of the ₁-, ₃- and ₄-vectors. This ellipsoid is a generalization of the Beltrami ellipsoid for isotropic materials.Furthermore, in combination with extensive experimental evidence, this theory indicates that the eigenangle alone monoparametrically characterizes the degree of anisotropy for each transversely isotropic material. Thus, while the angle for isotropic materials is always equal to _i = 125.26° and constitutes a minimum, the angle || progressively increases within the interval 90–180° as the anisotropy of the material is increased. The anisotropy of the various materials, exemplified by their ratiosE _L/2G_L of the longitudinal elastic modulus to the double of the longitudinal shear modulus, increases rapidly tending asymptotically to very high values as the angle approaches its limits of 90 or 180°.     

A solution of the E-ε model for nearly homogeneous turbulence with a mean shear

An analytical solution of the E- model for the downstream evolution of a stationary and nearly homogeneous turbulent shear flow is presented. In case that the turbulent time scale has adjusted itself to the time scale imposed by the shear, an asymptotic solution can be derived from the full solution, which shows that both E and increase downstream exponentially. By comparing this asymptotic solution with experimental data a value for the unknown constant c _l , in the -equation, is derived. Moreover, we find an expression for the downstream development of the variance of a scalar, which is also compared with experimental data. The analytical solution shows that a homogeneous shear flow with a uniform velocity gradient can only be obtained if the shear is sufficiently small. In the experiments this condition is not always satisfied. A discussion is given on how a nearly homogeneous shear flow can be obtained over a limited downstream interval by changing the initial conditions in E and , and a comparison is made with experimental data. Finally it is shown that better transverse homogeneity can be obtained by taking an exponential velocity profile instead of a linear profile.     

Near-wall turbulence models for 3D boundary layers

Calculations of the three-dimensional boundary layer in an S shaped duct are performed with various – models. Three different near-wall models are used for the – model, of which one is using a new set of near-wall damping functions deduced from direct numerical simulations of turbulent channel flow available in the literature. The results show that it is possible to obtain damping functions giving better agreement, especially for and , with direct simulation data and experiments than with damping functions deduced from trial and error.     

Turbulent flow investigation inside and outside plain-orifice atomizers with rounded orifice inlets

The performances of three linear eddy viscosity models (LEVM) and one algebraic Reynolds stress model (ARSM) on the simulation of the internal and external flows in the plain-orifice atomizers with rounded orifice inlets are evaluated. The validity of the computational model is first assessed through the testing of a backward facing step flow, a sudden expansion pipe flow and a liquid column collapsing problem. Then the atomizer internal and external flows are analyzed by comparing the computed discharge coefficients with available experimental data and by comparing the turbulence intensity profiles at the orifice exit. The results are also illustrated by the fluid/air interface plot. It is found that the turbulence models investigated exhibit zonal behaviors, i.e., none of the models investigated performs well throughout the entire flow field. It is worthwhile to note that the standard k-model is not necessarily the worst among the models investigated. In average, the ARSM model gives better results as compared to the standard k-model and the low Reynolds number models. The turbulence strength has a significant influence on the global characteristics of the flow field. The models with better predictions of the turbulence kinetic energy, such as Gatski–Speziales ARSM model and Nagano–Hishidas low Reynolds number model, can yield better predictions of the global characteristics of the flow field, e.g., the reattachment lengths for the backward-facing step flow and the sudden expansion pipe flow, and the discharge coefficient for the atomizer flow.     

Fluctuating flow in a rotating fluid

Summary Fluctuating flow of a viscous fluid rotating over a disk whose angular velocity oscillates about a nonzero mean is investigated. Initially the disk and the fluid rotate in the same sense with different angular velocities ₁ and ₂ ( ₂> ₁) and at a particular instant of time, the angular velocity of the disk becomes ₁[1+ sin( )]. The problem is solved as an initial boundary value problem and it is found that for small values of the results of analytical and numerical methods are in excellent agreement. The effect of frequency parameter on surface skin frictions has been analysed for various values of angular velocity ratio s and amplitude parameter .

---

Fluktuierende Strömung in einer rotierenden Flüssigkeit

Übersicht Untersucht wird die fluktuierende Strömung einer viskosen Flüssigkeit, die über einer Scheibe, deren Winkelgeschwindigkeit um einen von Null verschiedenen Mittelwert schwankt, rotiert. Anfangs drehen sich die Scheibe und die Flüssigkeit gleichsinnig, aber mit verschiedenen Winkelgeschwindigkeiten ₁ und ₂ ( ₂> ₁). Zu einem Anfangszeitpunkt geht die Winkelgeschwindigkeit der Scheibe über in ₁[1+ sin ( )]. Die Aufgabe wird als Anfangs-/Randwertproblem gelöst. Für kleine Werte stimmen die analytischen und numerischen Ergebnisse hervorragend überein. Für verschiedene Werte des Winkelgeschwindigkeitsverhältnisses und des Amplitudenparameters wurde der Einfluß des Frequenzparameters auf die Reibspannungen an der Scheibe untersucht.

     

Flow model for velocity distribution in fixed porous beds under isothermal conditions

In a brief survey of the previous work the limitations of the modified Darcy equation and of the vectorial form of the Ergun equation are discussed. To include the effect of wall friction on the flows the viscous resistance term is added to the vectorial form of the Ergun equation. Using the generalized Ergun equation a one-dimensional formulation is presented for flow of fluids through packed beds taking into account the variation of porosity along the radial direction. It is found that there is a reasonable agreement between the numerical and the experimental results and it is observed that the variation of porosity with radial position has greater influence on channeling of velocity near the walls. For the assumption of constant porosity the velocity profiles exhibit similar nature as the plug flow profiles with a thin boundary layer near the wall.

---

Modell der Geschwindigkeitsverteilung in einem isotherm durchströmten Festbett

Zusammenfassung In der vorliegenden Arbeit werden eingangs die Anwendbarkeitsgrenzen der modifizierten Darcy-Gleichung und der in vektorieller Form geschriebenen Ergun-Gleichung diskutiert. Um Einflüsse der Wandreibung auf eine Strömung mit in der Ergun-Gleichung berücksichtigen zu können, wird ein Reibungsterm hinzugefügt. Die so generalisierte Gleichung kann benutzt werden, um die eindimensional gerichtete Strömung durch eine Kugel schüttung zu berechnen. Eine radiale Veränderung der Schüttungsporosität ist dabei mit in die Betrachtung eingeschlossen. Das nichtlineare Grenzwertproblem wird numerisch gelöst und mit experimentellen Daten aus der Literatur verglichen. Die mit Meßwerten zufriedenstellend übereinstimmenden Rechenergebnisse zeigen, daß die radiale Porositätsverteilung in einem Festbett einen erheblichen Einfluß auf die Durchströmungsgeschwindigkeit in Wandnähe ausübt; die Berechnungen geben die Strömungsrandgängigkeit wieder. Wird die Bettporosität als unveränderlich angenommen, erhält man pfropfenströmungsähnliche Geschwindigkeitsprofile mit einer dünnen Wandgrenzschicht, in welcher die Geschwindigkeit auf den Wert null abfällt.

---

Nomenclature A Tridiagonal matrix defined in Eq. (20) - a Bed radius - d_p Particle diameter - f₁ 150 (1–)²/(³d _p ² ) Darcy resistance term - f₂ 1,75(1–)/(³d_p) Parameter of resistance due to inertial effects - ¯f₁ 150(1–)²/³ - ¯f₂ 1,75(1–)/₃ - G Column vector defined in Eq. (20) - k Permeability, /f₁ - L Length of the bed - P Pressure - r Radial co-ordinate - R_p Reynolds number based on particle diameter, v₀d_p/ - , v_z Superficial velocity vector, axial component - v_1z Average superficial velocity defined in Eq. (20) - V Absolute magnitude of velocity - ¯v The average velocity - v₀ The velocity at the centre of the tube - X Column vector defined in Eq. (20) - r^* Dimensionless radial co-ordinate, r/a - p^* Dimensionless pressure, p/v ₀ ² - v _z ^* Dimensionless axial component of velocity, v_z/v₀ - ¯v^* Dimensionless average velocity defined in Eq. (20) - z^* Dimensionless axial co-ordinate, z/L Greek letters Ratio of tube radius to particle diameter, a/d_p - Porosity or void fraction - ₀ Porosity at the axis of the container - Dynamical viscosity - Kinematic viscosity - p Density - Distance from the wall of the container, defined in Eq. (16)     

Numerical solution of quasilinear singularly perturbed ordinary differential equation without turning points

In this paper we consider a quasilinear second order ordinary diferential equation with a small parameter ε, Firstly an approximate problem is constructed. Then an iterative procedure is developed. Finally we give an algorithm whose accuracy is good for arbitrary ε>0.