On the perturbational global attractivity of nonautomous delay differential equations

Consider the perturbed nonautonomous linear delay differential equation x(t) = - a(t)x(t-τ) + F(t, x₁, t ⩾ 0 where x₁(s)=x(t+s) for −δ≤s≤0. Suppose that a(t) ∈ C([0,∞), (0,∞)), τ≥0,F:[0, ∞) x C[−δ,0] → R is a continuous functions and F(t,0) ≡ 0. Here C[−δ,0] is the space of continuous functions Φ: [−δ,0] → R with ∥Φ∥<H for the norm | Φ |, where |·| is any norm in R and 0<H≤+∞. Most of the known papers [1–5,7] have been concerned with the local or global asymptotic behavior of the zero solution of Eq. (*) when a(t) is independent of t i. e., a(t) is autonomous. The aim in this paper is to derive the sufficient conditions for the global attractivity of the zero solution of of Eq. (*) When a(t) is nonautomous. Our results, which extend and improve the known results, are even “sharp”. At the same time, the method used in this paper can be applicable to the perturbed nonlinear equation. Project supported by the Natural Science Foundation of Hunan     

Measure-theoretic foundations of statistical mechanics

The basic formulas of classical equilibrium statistical mechanics are derived from well-known theorems in measure theory and ergodic theory. The method used is a generalization of the methods of Khinchin and Grad and deals with several, in fact a “complete set”, of “invariants” or “integrals of the motion”. Most of the results are simple corollaries of Birkhoff's ergodic theorem, and since time-averages are used, the whole approach is characterized by an absence of statistical “ensembles” and probability notions. In the course of the development a “generalized temperature” is introduced, and a generalization of the second law of thermodynamics is derived. Formulas for the “microcanonical”, “canonical”, and “grand canonical” distributions appear as special cases of the general theory.     

One approximate solution of the Nekrasov problem

An approximate solution ω = A[ω, μ] of the nonlinear integral Nekrasov equation is obtained by successive replacement of the kernel of the integral operator by a close one. The solution is sought not directly at the bifurcation point μ₁ = 3 of the linearized equation ω = μL[ω] but at the point μ = 1 at which operator A[ω, μ], remaining nonlinear in ω, is linear in μ. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 50–56, November–December, 2007.     

ROSSBY WAVES WITH THE CHANGE OF β

In this paper, the change of the Rossby parameter β with latitude is considered and the parameter γ≡-dβ/dy=2sinφ/a² is introduced and the β-plane approximation is extended into f=f₀+β₀y-γ₀y²/2 which includes the parameter γ. Such approximation closes further to practice especially in the high latitude regions.We give emphasis to the research of the effect of the parameter γ on the Rossby waves. It is seen that the effect of the parameter γ is remarkable in the high latitude regions. It can produce the Rossby waves caused by the pure parameter γ. And the phase speed formula of Rossby waves with the change of ft is generally given, which is degenerated into the well-known Rossby formula when γ₀=0. The researches also point out that when the change of β is regarded, even if the basic current u is a linear function of y the unstable modes can also take place. However,the parameter γ usually plays a stable part in the Rossby waves and it does affect the longitudinal scale and the structure of constant phase lines(trough-ridge lines)of Rossby waves and slow down the growing or decaying of Rossby waves.     

The COVID-19 pandemic: model-based evaluation of non-pharmaceutical interventions and prognoses

For the electrically actuated microbeam subjected to a combination of DC and AC voltage loadings, we studied the nonlinear dynamical responses and internal resonance analytically. The flexible boundary condition is considered. The modal interaction due to three-to-one internal resonance between the first and second modes is highlighted. The method of multiple scales is employed to get the modulation equation, which describes the amplitude and phases of the involving modes. Then, the equilibrium solutions of the modulation equation are calculated and their stability is examined. The frequency and force response curves reflecting the primary resonance (the first mode) are presented. We find that the first and second modes in the proposed model are coupled, and hence the energy transfer can occur between the involving modes. Moreover, it is found that the response of the system may encounter Hopf bifurcation for some specific parameters. As for simulation, the Galerkin scheme is applied to verify the analytical results in terms of time history, phase-plane portrait, Fourier spectrum, and Poincare section. The simulation results are in good agreement with the analytical ones.     

Mini-Tensile Experiments of Clock-Rolled Zirconium Plate

We present our efforts to measure the tensile strength of clock-rolled pure zirconium in the through-thickness (TT) direction of the plate. Although the plate is too thin to produce standard ASTM tensile samples in the TT orientation, such measurements are relevant to benchmarking our constitutive models of hardening and texture evolution. We have designed a fixture and sample to perform tensile tests on our 9 mm thick plate. The sample is a double-ligament mini-tensile sample: 8 × 8 × 1 mm overall; each ligament has a gage section of 1 × 1 × 3 mm. In contrast, our standard “macro” tensile sample is a flat dogbone with a gage section of 3 × 1.5 × 25 mm. We validate our design by comparing the results of mechanical tests performed on samples of both geometries. Although the hardening response is nearly identical, the flow stress of the miniature samples is offset by +25 MPa at the onset of plastic yield. We present our efforts to resolve the origin of this offset.     

Special case of a traveling wave in one model of a multicomponent medium

The Kuropatenko model for a multicomponent medium whose components are polytropic gases is considered. It is assumed that, as x → ±∞, the multicomponent medium is in a homogeneous state with constant gas-dynamic parameters — velocity, pressure, and temperature. For the traveling wave flows, conditions similar to the Hugoniot conditions are obtained and used to uniquely determine the flow parameters for x → −∞ from the flow parameters x → +∞ and traveling wave velocity. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 39–47, July–August, 2009.     

Applications of the signs of Melnikov's function

The existence and stability of peridic solutions for the two-dimensional system x′=f(x)+g(x, a), O<ε<1,a∈R whoseunperturbed system is Hamlitonan be decided by using the signs of Melnikov's function. Theresults can be applied to the comstuction of phase portralts in the bifiacation set of codimension two bifurcations of flows with double aero eignvahues.     

Existence of the minimal positive solution of some nonlinear elliptic systems when the nonlinearity is the sum of a sublinear and a superlinear term

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｉｓｐａｐｅｒ,ｗｅｃｏｎｓｉｄｅｒｔｈｅｅｌｌｉｐｔｉｃｓｙｓｔｅｍ(1λ)　-Δｕ=ｆ(λ,ｘ,ｕ)-ｖ　　(ｉｎΩ),-Δｖ=δｕ-γｖ(ｉｎΩ),ｕ=ｖ=0(ｏｎΩ),ｗｈｅｒｅΩｉｓａｓｍｏｏｔｈｂｏｕｎｄｅｄｄｏｍａｉｎｉｎＲＮ(Ｎ≥2)ａｎｄλｉｓａｒｅａｌｐａｒａｍｅｔｅｒ.Ｔｈｅｓｏｌｕｔｉｏｎｓ(ｕ,ｖ)ｏｆｔｈｉｓｓｙｓｔｅｍｒｅｐｒｅｓｅｎｔｓｔｅａｄｙｓｔａｔｅｓｏｌｕｔｉｏｎｓｏｆｒｅａｃｔｉｏｎｄｉｆｆｕｓｉｏｎｓｙｓｔｅｍｓｄｅｒｉｖｅｄｆｒｏｍｓｅｖｅｒａｌａｐ…     

Continuity of a Deformation as a Function of its Cauchy-Green Tensor

If the Riemann-Christoffel tensor associated with a field C of class 𝒞² of positive-definite symmetric matrices of order 3 vanishes in a simply connected open subset Ω⊂ℝ³, then this field is the Cauchy-Green tensor field associated with a deformation Θ of class 𝒞³ of the set Ω, and Θ is uniquely determined up to isometries of ℝ³. Let denote the equivalence class formed by all such deformations, and let denote the mapping defined in this fashion. We establish here that the mapping ℱ is continuous, for certain natural metrizable topologies. (Accepted October 14, 2002) Published online March 12, 2003 Communicated by S. S. Antman     

Viscosity of concentrated noncolloidal bidisperse suspensions

At the same solid volume fraction (Φ) the relative viscosity (η _r ) of a concentrated noncolloidal bidisperse suspension of hard spherical particles is lower than that of a monodisperse suspension. In this paper a semi-analytical viscosity model of noncolloidal bidisperse suspensions is derived using an integration method. In this model the random loose packing density obtained by computer simulation is taken as the limit of solid volume fraction Φ _m which depends upon both the diameter ratio (λ) of large to small particles and the volume fraction of large particles (ξ=Φ _l /Φ). This model shows that at high solid volume fraction, Φ > 0.40, both λ and ξ significantly influence η _r . For example, at Φ=0.5, it predicts that for monodisperse suspensions η _r =70, while for bidisperse suspensions (λ=2 and ξ=0.7) η _r =40. Comparison shows that, at high solid volume fraction (0.4–0.5), the relative viscosity predicted by this model is in good agreement with that predicted by the work of Shapiro and Probstein (1992) and of Patlazhan (1993), but is higher than that predicted by the work of others. Received: 27 February 2001 Accepted: 25 April 2001     

Turbulence Modeling Effects on the Prediction of Equilibrium States of Buoyant Shear Flows

The effects of turbulence modeling on the prediction of equilibrium states of turbulent buoyant shear flows were investigated. The velocity field models used include a two-equation closure, a Reynolds-stress closure assuming two different pressure-strain models and three different dissipation rate tensor models. As for the thermal field closure models, two different pressure-scrambling models and nine different temperature variance dissipation rate ɛ_τ) equations were considered. The emphasis of this paper is focused on the effects of the ɛ_τ-equation, of the dissipation rate models, of the pressure-strain models and of the pressure-scrambling models on the prediction of the approach to equilibrium turbulence. Equilibrium turbulence is defined by the time rate of change of the scaled Reynolds stress anisotropic tensor and heat flux vector becoming zero. These conditions lead to the equilibrium state parameters, given by /ɛ, _τ/ɛ_τ, , Sk/ɛ and G/ɛ, becoming constant. Here, and _τ are the production of turbulent kinetic energy k and temperature variance , respectively, ɛ and ɛ_τ are their respective dissipation rates, R is the mixed time scale ratio, G is the buoyant production of k and S is the mean shear gradient. Calculations show that the ɛ_τ-equation has a significant effect on the prediction of the approach to equilibrium turbulence. For a particular ɛ_τ-equation, all velocity closure models considered give an equilibrium state if anisotropic dissipation is accounted for in one form or another in the dissipation rate tensor or in the ɛ-equation. It is further found that the models considered for the pressure-strain tensor and the pressure-scrambling vector have little or no effect on the prediction of the approach to equilibrium turbulence. Received 21 April 2000 and accepted 21 February 2001     

Breakdown of the Reynolds Analogy in a Stagnation Region Under Inflow Disturbances

A systematic analysis is performed for the Reynolds analogy breakdown at stagnation-region flow and heat transfer in the presence of inflow disturbances. The Reynolds analogy breakdown between momentum and energy transfers in a stagnation region is scrutinized by varying the Reynolds number (5000≤Re≤20000), the amplitude (0.00075≤A≤0.003) and the length scale (λ/δ=10.6). A spanwise sinusoidal variation is superimposed on the velocity component normal to the wall. Self-similarity solutions are obtained with trigonometric series expansions. The Reynolds analogy criterion demonstrates that the rate of change of skin friction is different from that of wall heat transfer. Different evolutions of the rates of skin friction and wall heat transfer are due to the difference between 〈s'v'〉 and 〈v'T'〉. An in-depth analysis on 〈s'v'〉 and 〈 v'T'〉 is performed by analysis using disturbance correlations based on the fluctuating velocity transport equations in vorticity form. It is found that the pressure fluctuations, the wall blocking and the Lamb vectors are responsible for the breakdown of the Reynolds analogy. A direct comparison is made between momentum and energy balances associated with the three responsible mechanisms. A common finding is that their profiles are changed significantly at a location where the evolution of the streamwise vortex is strong. Received 12 May 2000 and accepted 6 March 2001     

DNS of Mixing and Reaction of Two Species in a Turbulent Channel Flow: A Validation of the Conditional Moment Closure

We consider the chemical reaction in a turbulent flow for the case that the time scale of turbulence and the time scale of the reaction are comparable. This process is complicated by the fact that the reaction takes place intermittently at those locations where the species are adequately mixed. This is known as spatial segregation. Several turbulence models have been proposed to take the effect of spatial segregation into account. Examples are the probability density function (PDF) and the conditional moment closure (CMC) models. The main advantage of these models is that they are able to parameterize the effects of turbulent mixing on the chemical reaction rate. As a price several new unknown terms appear in these models for which closure hypothesis must be supplied. Examples are the conditional dissipation 〈 χ ∣ φ 〉, the conditional diffusion 〈 κ ∇² φ ∣ u, φ 〉 and the conditional velocity 〈 u ∣ φ 〉. In the present study we investigate these unknown terms that appear in the PDF and CMC model by means of a direct numerical simulation (DNS) of a fully developed turbulent flow in a channel geometry. We present the results of two simulations in which a scalar is released from a continuous line source. In the first we consider turbulent mixing without chemical reaction and in the second we add a binary reaction. The results of our simulations agree very well with experimental data for the quantities on which information is available. Several closure hypotheses that have been proposed in the literature, are considered and validated with help of our simulation results. This revised version was published online in July 2006 with corrections to the Cover Date.     

Thermal convection in a porous layer: Effects of anisotropy and surface boundary conditions

The theory describing the onset of convection in a homogeneous porous layer bounded above and below by isothermal surfaces is extended to consider an upper boundary which is partly permeable. The general boundary condition p + λ ∂p/∂n = constant is applied at the top surface and the flow is investigated for various λ in the range 0 ⩽ λ < ∞. Estimates of the magnitude and horizontal distribution of the vertical mass and heat fluxes at the surface, the horizontally-averaged heat flux (Nusselt number) and the fraction of the fluid which recirculates within the layer are found for slightly supercritical conditions. Comparisons are made with the two limiting cases λ → ∞, where the surface is completely impermeable, and λ = 0, where the surface is at constant pressure. Also studied are the effects of anisotropy in permeability, ξ = K _H /K _V , and anisotropy is thermal conductivity, η = k _H /k _V , both parameters being ratios of horizontal to vertical quantities. Quantitative results are given for a wide variety of the parameters λ, ξ and η. In the limit ξ/η → 0 there is no recirculation, all fluid being converted out of the top surface, while in the limit ξ/η → ∞ there is full recirculation.     

Slow Dynamics for the Cahn‐Hilliard Equation in Higher Space Dimensions: The Motion of Bubbles

It is known that the Van der Waals‐Cahn‐Hilliard (W‐C‐H) dynamics can be approximated by a Quasi‐static Stefan problem with surface tension. It turns out that the Stefan problem has a manifold of equilibria equal in dimension to that of the domain Ω: any sphere of fixed radius with interface contained in the domain is an equilibrium (indistinguishable from the point of view of the perimeter functional). We resolve this degeneracy by showing that at the W‐C‐H level this manifold is replaced by a quasi‐invariant stable manifold, on which the typical solution preserves its “bubble” like shape until it reaches the boundary. Moreover, we show that the “bubble” moves superslowly. We also obtain an equation that determines those special spheres that correspond to equilibria at the W‐C‐H level. Our work establishes the phenomenon of superslow motion in higher space dimensions in the class of single interface solutions. (Accepted February 12, 1996)     

Improvement of cold-wire response for measurement of temperature dissipation

This work aimed at improving fine-scale measurements using cold-wire anemometry. The dissipation ɛ _θ of the temperature variance was measured on the axis of a heated turbulent round jet. The measurements were performed with a constant current anemometer (CCA) operating fine Pt–10%Rh wires at very low overheat. The CCA developed for this purpose allowed the use of the current injection method in order to estimate the time constant of the wire. In the first part of the paper, it is shown that the time constants obtained for two wire diameters −d=1.2 and d=0.58 μm – compare well with those measured at the same time using two other methods (laser excitation and pulsed wire). Moreover, for these two wires, the estimated time constants were in good agreement with those obtained from a semi-empirical relation. In the second part of the paper, a compensation procedure – post-processing filtering – was developed in order to improved the frequency response of the cold-wire probes. The measurements carried out on the axis of the jet (Re _D =16 500, Re _λ  ≃ 167) showed that the frequency response of the 1.2 μm wire was significantly improved. In fact, the spectral characteristics of the compensated signal obtained with the 1.2 μm wire compared fairly well with those from the 0.58 μm wire. Moreover, the results indicated that the compensation procedure must be applied when the cut-off frequency of the cold-wire f _c is lower than two times the Kolmogorov frequency f _K. In the case where f _c ≃ 0.6f _K, the compensation procedure can reduce the error in the estimate of ɛ _θ by more than 20%. When f _c ≃ 2f _K, the effect of the compensation is reduced to about 5%. Received: 3 November 2000/Accepted: 23 March 2001     

Existence of Singular Self-Similar Solutions of the Three-Dimensional Euler Equations in a Bounded Domain

A self-similar solution of the three-dimensional (3d) incompressible Euler equations is defined byu(x,t) =U(y)/t*-t) ^{α, y = x/(t* ~ t)β,α,β>} 0, whereU(y) satisfiesζU + βy. ΔU + U. VU + VP = 0,divU = 0. For α = β = 1/2, which is the limiting case of Leray’s self-similar Navier—Stokes equations, we prove the existence of(U,P) ε H³(Ω,R³ X R) in a smooth bounded domain Ω, with the inflow boundary data of non-zero vorticity. This implies the possibility that solutions of the Euler equations blow up at a timet = t*, t* < +∞.     

The propagation mechanism of high speed turbulent deflagrations

The propagation regimes of combustion waves in a 30 cm by 30 cm square cross–sectioned tube with an obstacle array of staggered vertical cylindrical rods (with BR=0.41 and BR=0.19) are investigated. Mixtures of hydrogen, ethylene, propane, and methane with air at ambient conditions over a range of equivalence ratios are used. In contrast to the previous results obtained in circular cross–sectioned tubes, it is found that only the quasi–detonation regime and the slow turbulent deflagration regimes are observed for ethylene–air and for propane–air. The transition from the quasi–detonation regime to the slow turbulent deflagration regime occurs at (where D is the tube “diameter” and is the detonation cell size). When , the quasi–detonation velocities that are observed are similar to those in unobstructed smooth tubes. For hydrogen–air mixtures, it is found that there is a gradual transition from the quasi–detonation regime to the high speed turbulent deflagration regime. The high speed turbulent deflagration regime is also observed for methane–air mixtures near stoichiometric composition. This regime was previously interpreted as the “choking” regime in circular tubes with orifice plate obstacles. Presently, it is proposed that the propagation mechanism of these high speed turbulent deflagrations is similar to that of Chapman–Jouguet detonations and quasi-detonations. As well, it is observed that there exists unstable flame propagation at the lean limit where . The local velocity fluctuates significantly about an averaged velocity for hydrogen–air, ethylene–air, and propane–air mixtures. Unstable flame propagation is also observed for the entire range of high speed turbulent deflagrations in methane–air mixtures. It is proposed that these fluctuations are due to quenching of the combustion front due to turbulent mixing. Quenched pockets of unburned reactants are swept downstream, and the subsequent explosion serves to overdrive the combustion front. The present study indicates that the dependence on the propagation mechanisms on obstacle geometry can be exploited to elucidate the different complex mechanisms of supersonic combustion waves. Received 5 November 2001 / Accepted 12 June 2002 / Published online 4 November 2002 Correspondence to: J. Chao (e-mail: jenny.chao@mail.mcgill.ca) An abridged version of this paper was presented at the 18th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Seattle, USA, from July 29 to August 3, 2001.     

The method of composite expansions applied to boundary layer problems in symmetric bending of the spherical shells

In this paper,the method of composite expansions which was proposed by W. Z. Chien (1948)[5]is extended to investigate two-parameter boundary layer problems.For the problems of symmetric deformations of the spherical shells under the action of uniformly distribution load q, its nonlinear equilibrium equations can be written as follows: where ε and δ are undetermined parameters.If δ=1 and ε is a small parameter, the above-mentioned problem is called first boundary layer problem; if ε is a small parameter, and δ is a small parameter, too, the above-mentioned problem is called second boundary layer problem.For the above-mentioned problems, however, we assume that the constants ε, δ and p satisfy the following equation: εp=1-ε In defining this condition by using the extended method of composite expansions, we find the asymptotic solution of the above-mentioned problems with the clamped boundary conditions.