An unsymmetrical van der pol equation with stable subharmonics

Possible stable subharmonic solutions of the equation

ÿ − k(1 + 2cy − y²)ÿ + Y = bkμ cos μt, c > 0

, klarge, are discussed by the techniques used by J.E. Littlewood for van der Pol's equation in Acta Math. 97 (1957), that is the case of the above equation with c = 0 and

, k large. Their variation as c increases is also considered briefly.     

Non-linear oscillations with multiple forcing terms

The problem of determining the transient response of a non-linear oscillator of the form ü + u = εƒ(u,u) + E(t) is studied by the method of multiple time scales, using the symbolic computation system MACSYMA. when the excitation E(t) consists of a finite number of harmonic forcing terms. Here ε is a small parameter and ƒ(u,u) is a non-linear function of its arguments. In particular, the Van der Pol and Duffing oscillators are studied in detail. It is found that when the forcing frequencies are not close to each other or close to the primary resonance of the system, then the response of the system is analogous to the behavior when only one forcing term is present. However, when the forcing frequencies are close to each other or close to the primary resonance, then the behavior is quite different, exhibiting certain oscillations not observed in the case of one forcing term.     

Group properties of utt=[F(u)ux]x

The group properties and the associated Lie algebra are developed for the subject quasilinear wave equation, for arbitrary f[fεC²(R), f > 0, f ≠ 0]. From the resulting information sets of explicit invariant solutions are constructed for wave propagation in gases and for the transonic equation.     

Inertial instability of the Stewartson -layer

Inertial stability of a vertical shear layer (Stewartson E^1/4-layer) on the sidewall of a cylindrical tank with respect to stationary axisymmetric perturbations is inverstigated by means of a linear theory. The stability is determined by two non-dimensional parameters, the Rossby number Ro = U/2ΩL and Ekman number E = v/ΩH², where U and L = (E/4)1/4H are the characteristic velocity and width of the shear layer, respectively, Ω the angular velocity of the basic rotation, v the kinematic viscosity and H the depth of the tank.

For a given Ekman number, the flow is more unstable for larger values of the Rossby number. For E = 10⁻⁴, which is a typical value of the Ekman number realized in rotating tank experiments, the critical Rossby number Ro_c for instability and the critical axial wavenumber m_c non-dimensionalized by L⁻¹ are found to be 1.3670 and 8.97, respectively. The value of Ro_c increases and that of m_c decreases with increasing E.     

An application of Nagumo''s lemma to some singularly perturbed systems

The existence and asymptotic behavior as ε → 0⁺ of periodic, almost periodic, and bounded solutions of the differential system x = f(t, x, y, ε), Ωy′ = g(t, x, y, ε), are considered where x, f; are n-vectors, y, g are m-vectors and Ω = diag{ε^h₁}…, ε^h_m for integral h_i, h₁ h₂ …, h_m. The principal tools are a lemma of Nagumo which allows the construction of appropriate upper and lower solutions and the asymptotic theory of singularly perturbed linear differential systems.     

Flow in a cylinder driven by rotating split-disk endwalls

Flow of an incompressible viscous fluid contained in a cylindrical vessel (radius R, height H) is considered. Each of the cylinder endwalls is split into two parts which rotate steadily about the central axis with different rotation rates: the inner disk (r < r₁) rotating at Ω₁, and the outer annulus (r₁ < r < R) rotating at Ω₂. Numerical solutions to the axisymmetric Navier-Stokes equations are secured for small system Ekman numbers E ( v/(ΩH²)). In the linear regime, when the Rossby number Ro

Full-size image (1K)

, the numerical results are shown to be compatible with the theoretical prediction as well as the available experimental measurements. Emphasis is placed on the results in the nonlinear regime in which Ro is finite. Details of the structures of azimuthai and meridional flows are presented by the numerical results. For a fixed Ekman number, the gross features of the flow remain qualitatively unchanged as Ro increases. The meridional flows are characterized by two circulation cells. The shear layer is a region of intense axial flow toward the endwall and of vanishing radial velocity. The thicknesses of the shear layer near r = r₁ and the Ekman layer on the endwall scale with E and E , respectively. The numerical results are consistent with these scalings.     

On a Korteweg–de Vries-like equation with higher degree non-linearity

We deal with a non-linear partial differential equation which has been widely investigated owing to its applications in quantum field theory, as well as plasma and solid-state physics. It is the matter of a third order KdV-like equation with higher degree non-linearity in the coefficient of the transport term; it can be derived from a Lagrangian or an Hamiltonian density. In the current literature specific attention has been devoted to the search for traveling-wave solutions, depending upon a positive parameter v, which assesses the speed of the solitary wave. The velocity v is always assumed to be constant, as its dependence on the wave-amplitude is neglected in the mathematical model. In this context, Coffey [On series expansions giving closed-form solutions of Korteweg–de Vries-like equations, SIAM J. Appl. Math. 50 (6) (1990) 1580–1592] exploits an algebraic recursive technique to obtain these solutions in closed form for particular values of v. The aim of this paper is to extend these results by showing that closed-form solutions are achievable for every value of v: to this purpose we supply a proper mathematical framework for these issues by taking into account a suitable special function, namely an elliptic function in the sense of Weierstraß. Furthermore we obtain two classes of the so-called kink solutions, see [M.W. Coffey, On series expansions giving closed-form solutions of Korteweg–de Vries-like equations, SIAM J. Appl. Math. 50 (6) (1990) 1580–1592; B. Dey, Domain wall solutions of KdV-like equations with higher order non-linearity, J. Phys. A 19 (1) (1986) L9–L12], as well as an exponential development of the general solution, for which we prove the convergence. Eventually we show how to implement the resulting functions by means of a symbolic manipulation program.     

First order piecewise linear systems with random parametric excitation

The Fokker-Planck equation is used to develop a general method for finding the spectral density and other properties of first order systems governed by stochastic differential equations of the form

dx/dt + f(x)[1 + m(t)] = n(t)

, where f(x) is piecewise linear and m(t) and n(t) represent stationary Gaussian white noise. The method is similar to one used by the authors to deal with the case m(t) = 0, but is complicated by the possible existence of irregular (singular) points of the Fokker-Planck equation. Graphical results for some special cases are presented.     

Experimental study of turbulence in isothermal and stably stratified homogeneous shear flows: Mean flow measurements

The evolution of freestream turbulence under the combined action of linear shear and stable linear temperature profile is investigated. The experiment is carried out in a small, open circuit, low-speed test cell that uses air as working fluid. The temperature gradient formed at the entrance to the test section by means of an array of 24 horizontal, differentially heated elements is varied to get a maximum Brunt-Vaisala frequency N_o[=({g/T_m}{∂T/∂y})^1/2] of 3.1⁻¹. Linear velocity profiles are produced using screens of variable mesh size. The Reynolds number Re_M based on centre-line velocity and mesh size is varied from 80 to 175. Isothermal studies are carried out in four different experiments with varying velocity gradients. The effect of inlet turbulence level on growth of turbulence is studied in these flows by keeping the shear parameter Sh (=(x/u)(∂u/∂y)) constant. The range of shear parameters considered is 2.5–7.0. Shear and stratification combined produce a maximum gradient Richardson number Ri_g (= N_o²/(∂u/∂y)²) of 0.0145. Results have been presented in terms of evolution of variance of velocity fluctuations, Reynolds shear stress and temperature fluctuations. Measurements show the following: In isothermal flows the growth rate of turbulence quantities depends on both shear parameter and inlet turbulence level. There are distinct stages in the evolution of the flow and that can be identified by the power-law exponent of growth of turbulence. Shear is seen to promote the growth of turbulence and accelerate it towards a fully developed equilibrium state. Stratification initially suppresses the growth of turbulence and, hence, enhances the degree of underdevelopment. Under these conditions shear becomes active and subsequently enhances the growth rate of turbulence quantities.     

Asymptotic soliton train solutions of Kaup–Boussinesq equations

Asymptotic soliton trains arising from a ‘large and smooth’ enough initial pulse are investigated by the use of the quasiclassical quantization method for the case of Kaup–Boussinesq shallow water equations. The parameter varying along the soliton train is determined by the Bohr–Sommerfeld quantization rule which generalizes the usual rule to the case of ‘two potentials’ h₀(x) and u₀(x) representing initial distributions of height and velocity, respectively. The influence of the initial velocity u₀(x) on the asymptotic stage of the evolution is determined. Excellent agreement of numerical solutions of the Kaup–Boussinesq equations with predictions of the asymptotic theory is found.     

Exact analysis of Burgers equation on semiline with flux condition at the origin

The initial boundary value problem for the Burgers equation in the domain x 0, t > 0 with flux boundary condition at x = 0 has been solved exactly. The behaviour of the solution as t tends to infinity is studied and the “asymptotic profile at infinity” is obtained. In addition, the uniqueness of the solution of the initial boundary value problem is proved and its inviscid limit as → 0 is obtained.     

Weakly nonlocal envelope solitary waves: numerical calculations for the Klein-Gordon (φ) equation

“Weakly nonlocal” solitary waves differ from ordinary solitary waves by possessing small amplitude, oscillatory “wings” that extend indefinitely from the large amplitude “core”. Such generalized solitary waves have been discovered in capillarygravity water waves, particle physics models, and geophysical Rossby waves. In this work, we present explicit calculations of weakly nonlocal envelope solitary waves. Each is a sine wave modulated by a slowly-varying “envelope” that itself propagates at the group velocity. Our example is the cubically nonlinear Klein-Gordon equation, which is a model in particle physics (φ⁴ theory) and in electrical engineering (with a different sign). Both cases have weakly nonlocal“breather” solitons. Via the Lorentz invariance, each breather generates a one-parameter family of nonlocal envelope solitary waves. The φ⁴ breather was described and calculated in earlier work. This generates envelope solitons which have “wings” that are (mostly) proportional to the second harmonic of the sinusoidal factor. In this article, we calculate breathers and envelope solitary waves for the second, electrical engineering case. Since these, unlike the φ⁴ waves, contain only odd harmonics, the envelope solitary waves are nonlocal only via the third harmonic.     

Coherent fine scale eddies in turbulence transition of spatially-developing mixing layer

To investigate the relationship between characteristics of the coherent fine scale eddy and a laminar–turbulent transition, a direct numerical simulation (DNS) of a spatially-developing turbulent mixing layer with Re_ω,0 = 700 was conducted. On the onset of the transition, strong coherent fine scale eddies appears in the mixing layer. The most expected value of maximum azimuthal velocity of the eddy is 2.0 times Kolmogorov velocity (u_k), and decreases to 1.2u_k, which is an asymptotic value in the fully-developed state, through the transition. The energy dissipation rate around the eddy is twice as high compared with that in the fully-developed state. However, the most expected diameter and eigenvalues ratio of strain rate acting on the coherent fine scale eddy are maintained to be 8 times Kolmogorov length (η) and :β:γ = −5:1:4 in the transition process. In addition to Kelvin–Helmholtz rollers, rib structures do not disappear in the transition process and are composed of lots of coherent fine scale eddies in the fully-developed state instead of a single eddy observed in early stage of the transition or in laminar flow.     

Impingement heat transfer at a circular cylinder due to a submerged slot jet of water

An experimental investigation was carried out on the heat transfer due to a submerged slot jet of water impinging on a circular cylinder in crossflow. The cylinder diameter and the slot width are of the same order of magnitude, specifically D_s = 2.0 and 3.0 mm and D_c = 2.5 and 3.0 mm. The experimental apparatus allowed variation of the slot width, the cylinder diameter, and the distance from nozxle exit to heater. Conditions of impingement from the bottom (ascending flow) were taken into consideration as well as impingement from above (descending flow). The Nusselt number was determined as a function of Reynolds and Prandtl numbers in the range 1.5 × 10³ < Re < 2.0 × 10⁴, 2.7 < Pr < 7.0, and 1.5 ≤ z/D_s ≤ 10. The experimental data were correlated with a simple equation that fits 90% of the data with a precision of 20%.     

A survey of experimental heat transfer data for nucleate pool boiling of liquid metals and a new correlation

The experimental data for heat transfer during nucleate pool boiling of saturated liquid metals on plain surfaces are surveyed and a new correlation is presented. The correlation is h = Cq^0.7p_r^m, where C and m are, respectively, 13.7 and 0.22 p_r < 0.001 and 6.9 and 0.12 for p_r > 0.001 (h is in W/m² K and q in W/m²). This correlation has been verified with data for K, Na, Cs, Li, and Hg from 17 sources over the reduced pressure (p_r) range of 4.3 × 10⁻⁶ to 1.8 × 10⁻². The correlation of Subbotin et al. was found unsatisfactory, but a modified correlation was developed that also gives good agreement with most of the data.     

The horizontal intrusion layer of melt in a saturated porous medium

This article presents a theory of how the melt region advances as an intrusion layer along the top boundary of a solid phase-change material that is heated from the side. The phase-change material fills the pores of a solid matrix. We show that the thickness of the horizontal melt layer increases as x^3/5, where x is the horizontal distance measured by from the leading edge of the layer. The total length of the intrusion layer increases as t^3/4, and as T_max^5/4. Finite-difference simulations of convection melting in the Darcy-Rayleigh number range of 200–800 agree with the theoretical results. We also show that in a rectangular porous medium heated from the side, the size of the entire melt region is dominated by the melting contributed by the horizontal intrusion layer, if the time is great enough so that the group (Ste Fo)^3/4 is greater than 1.     

A general correlation for the local loss coefficient in Newtonian axisymmetric sudden expansions

Results from numerical simulations and guidance from an approximated corrected-theory, developed by Oliveira and Pinho (1997), (Oliveira, P.J. and Pinho, F.T. 1997. Pressure drop coefficient of laminar Newtonian flow in axisymmetric sudden expansions. Int. J. Heat and Fluid flow 18, 518–529) have been used to arrive at a correlation expressing the irreversible loss coefficient for laminar Newtonian flow in axisymmetric sudden expansions. The correlation is valid for the ranges 1.5 < D₂/D₁ < 4 and 0.5 < Re < 200 with errors of less than 5%, except for 25 < Re < 100 where the error could be as much as 7%. The recirculation bubble length is also presented for the same range of conditions and the pressure recovery coefficient was calculated for Reynolds numbers above 15.     

Equations differentielles non lineaires: Differentiabilite des solutions par rapport aux conditions initiales et problemes de la stabilite

In this paper, we examine the problem of the asymptotic stability of the solutions of a differential equation (E), Y = X(Y, t), from the point of view of the differentiability of the solutions with respect to the initial conditions. The method allows us to deal with cases in which the mapping X is not differentiable everywhere, and in which the variational equation of (E) is not defined in the usual meaning. This method can be carried on for the research of periodic solutions. We give two examples.

Zusammenfassung

In diesem Bericht wird das Problem der asymptotischen Stabilität der Lösungen einer Differential Gleichung (E), Y = X(Y, t), vom Standpunkt der Differenzierbarkeit der Lösungen in Bezug auf die Anfangsbedingungen untersucht. Das Verfahren erlaubt die Fälle zu behandeln, worin die Abbildung X nicht überall differenzierbar ist, und worin auch die Variations Gleichung nicht im üblichen Sinne definiert ist. Das Verfahren kann erweitert werden zur Untersuchung periodischer Lösungen. Zwei Beispiele werden vorgeführt.

Abstract

B aдннoй paбaтe иccлeдyeтcя вoпoc oб accимeтpиХecкoй ycтoйХибocти peшeний диффepeHциaлънoгo ypaвнeния (E), Y = X(Y, t), c тoХки зpeния диффepeнциpyeмocти peшeний пo oтнoшeнию к нaХaлъным ycлoвиям. Иcпoлъзyeмый мeтoд дaeт вoзмoжнocтъ иccлeдoвaтъ cлyХaи, кoгдa пpeoбpaзoвaниe X нe вcюдy диффepeнциaлнoe. Bapиaциoннoe ypaвнeниe для oпpeдeляeтcя нecкoлъкo инaХe Хeм oбыХнo. Дaнный мeтoд мoжeт иcпoлъзoвaтъcя для иccлeдoвa*ncy;ия пepиoдиХecкиx peшeний. Дaютcя двa пpимepa.