Estimation of Chaotic and Regular (Stick–Slip and Slip–Slip) Oscillations Exhibited by Coupled Oscillators with Dry Friction

In this paper, we present a novel approach to quantify regular or chaotic dynamics of either smooth or non-smooth dynamical systems. The introduced method is applied to trace regular and chaotic stick–slip and slip–slip dynamics. Stick–slip and slip–slip periodic and chaotic trajectories are analyzed (for the investigated parameters, a stick–slip dynamics dominates). Advantages of the proposed numerical technique are given.     

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.     

Vlasov–Maxwell–Boltzmann Diffusive Limit

Inspired by the work (Bastea et al. in J Stat Phys 1011087–1136, 2000) for binary fluids, we study the diffusive expansion for solutions around Maxwellian equilibrium and in a periodic box to the Vlasov–Maxwell–Boltzmann system, the most fundamental model for an ensemble of charged particles. Such an expansion yields a set of dissipative new macroscopic PDEs, the incompressible Vlasov–Navier–Stokes–Fourier system and its higher order corrections for describing a charged fluid, where the self-consistent electromagnetic field is present. The uniform estimate on the remainders is established via a unified nonlinear energy method and it guarantees the global in time validity of such an expansion up to any order.     

In situ force-balance tensiometry

 Although a fundamental physical parameter, surface tension is difficult to measure. Common tensiometry inaccuracy comes from failure to control air–liquid–solid contact conditions, or account for liquid meniscus geometry and buoyancy corrections. This paper describes an in situ tensiometry technique, based on withdrawal of a thin-walled tube from the liquid interface, that enforces a known air–liquid–solid contact condition. This technique can be pursued at any level of experimental hygiene. Experimental results for filtered tap water, an alcohol–water solution, and a surfactant–water solution show that results repeatable to three significant digits are obtained with modest effort for a variety of geometrical parameters. Received: 7 October 1997/Accepted: 23 April 1998     

From the Boltzmann Equation to an Incompressible Navier–Stokes–Fourier System

We establish a Navier–Stokes–Fourier limit for solutions of the Boltzmann equation considered over any periodic spatial domain of dimension two or more. We do this for a broad class of collision kernels that relaxes the Grad small deflection cutoff condition for hard potentials and includes for the first time the case of soft potentials. Appropriately scaled families of DiPerna–Lions renormalized solutions are shown to have fluctuations that are compact. Every limit point is governed by a weak solution of a Navier–Stokes–Fourier system for all time.     

Development and validation of a thermo-mechanical finite element model of the steel quenching process including solid–solid phase changes

A computational thermo-metallographic and thermoelastoplastic model for the analysis of the quenching process is developed and validated. The diffusive transfor-mations are modeled according to the Johnson–Mehl–Avrami–Kolmogorov model and the Scheil’s additivity rule. Two different models are investigated for the non-diffusive transformation—the Koistinen–Marburger model and the Yu model. A large displacement formulation is assumed for the deformation analysis, modeling the plastic behavior of the material according to the Prandtl–Reuss model. Two different bilinear hardening models—the isotropic and the kinematic hardening model—are used and compared. The model allows to evaluate the transient stress and strain distributions during the quenching process, the final phases and hardness distributions, and to predict the residual stress and the final deformation of the processed part. A good agreement between computational results and reference data is found     

Solving initial–boundary-value creep problems

The Galerkin–Bubnov method with global approximations is used to find approximate solutions to initial–boundary-value creep problems. It is shown that this approach allows obtaining solutions available in the literature. The features of how the solutions of initial–boundary-value problems for oneand three-dimensional models are found are analyzed. The approximate solutions found by the Galerkin–Bubnov method with global approximations is shown to be invariant to the form of the equations of the initial–boundary-value problem. It is established that solutions of initial–boundary-value creep problems can be classified according to the form of operators in the mathematical problem formulation     

Temperature dependence of the interfacial tension of PS/PMMA, PS/PE, and PMMA/PE blends

The effect of temperature on the interfacial tension for PS/PMMA, PS/PE, and PMMA/PE was measured using the imbedded fiber retraction method. Interfacial tensions for PS/PMMA, PS/PE, and PMMA/PE were measured over temperature ranges of 160–250 °C, 140–220 °C, and 140–220 °C, respectively. The interfacial tension was found to follow a dependence of 3.6–0.013 T dyn/cm, 7.6–0.051 T dyn/cm and 11.8–0.017 T dyn/cm for PS/PMMA, PS/PE, and PMMA/PE, respectively. Comparison of the data with the mean field theory of Helfand and Sapse were made; however, a simple linear fit to the data described the temperature dependence in the experimental window as well as the predictions of the mean field theory. Received: 6 July 1999 Accepted: 23 March 2000     

Characteristics of liquid sheets formed by splash plate nozzles

An experimental study was conducted to identify the effect of viscosity on the characteristics of liquid sheets formed by a splash plate nozzle. Various mixtures of corn syrup and water are used to obtain viscosities in the range 1–170 mPa.s. Four different splash plates with nozzle diameters of 0.5, 0.75, 1, and 2 mm, with a constant plate angle of 55° were tested. Liquid sheets formed under various operating conditions were directly visualized. The sheet atomization process for the range of parameters studied here is governed by two different mechanisms: Rayleigh–Plateau (R–P) and Rayleigh–Taylor (R–T) instabilities. R–P occurs at the rim and R–T occurs on the thin sheet. The rim instability can be laminar or turbulent, depending on the jet Reynolds number. The R–T instability of the sheet is observed at the outer edges of the radially spreading sheet, where the sheet is the thinnest. It can also occur inside the sheet, due to formation of holes and ruptures.     

On a Generalized Nonlinear K–ε Model and the Use of Extended Thermodynamics in Turbulence

A resent extension of the nonlinear K–ε model is critically discussed from a basic theoretical standpoint. While it was said in the paper that this model was formulated to incorporate relaxation effects, it will be shown that the model is incapable of describing one of the most basic such turbulent flows as is obvious but is described for clarity. It will be shown in detail that this generalized nonlinear K–ε model yields erroneous results for the Reynolds stress tensor when the mean strains are set to zero in a turbulent flow – the return-to-isotropy problem which is one of the most elementary relaxational turbulent flows. It is clear that K–ε type models cannot describe relaxation effects. While their general formalism can describe relaxation effects, the nonlinear K–ε model – which the paper is centered on – cannot. The deviatoric part of the Reynolds stress tensor is predicted to be zero when it actually only gradually relaxes to zero. Since this model was formulated by using the extended thermodynamics, it too will be critically assessed. It will be argued that there is an unsubstantial physical basis for the use of extended thermodynamics in turbulence. The role of Material Frame-Indifference and the implications for future research in turbulence modeling are also discussed. Received 19 February 1998 and accepted 23 October 1998     

Prediction of separation flows around a 6：1 prolate spheroid using RANS/LES hybrid approaches

This paper presents hybrid Reynolds-averaged Navier–Stokes (RANS) and large-eddy-simulation (LES) methods for the separated flows at high angles of attack around a 6:1 prolate spheroid. The RANS/LES hybrid methods studied in this work include the detached eddy simulation (DES) based on Spalart–Allmaras (S–A), Menter’s k–ω shear-stress-transport (SST) and k–ω with weakly nonlinear eddy viscosity formulation (Wilcox–Durbin+, WD+) models and the zonal-RANS/LES methods based on the SST and WD+ models. The switch from RANS near the wall to LES in the core flow region is smooth through the implementation of a flow-dependent blending function for the zonal hybrid method. All the hybrid methods are designed to have a RANS mode for the attached flows and have a LES behavior for the separated flows. The main objective of this paper is to apply the hybrid methods for the high Reynolds number separated flows around prolate spheroid at high-incidences. A fourth-order central scheme with fourth-order artificial viscosity is applied for spatial differencing. The fully implicit lower–upper symmetric-Gauss–Seidel with pseudo time sub-iteration is taken as the temporal differentiation. Comparisons with available measurements are carried out for pressure distribution, skin friction, and profiles of velocity, etc. Reasonable agreement with the experiments, accounting for the effect on grids and fundamental turbulence models, is obtained for the separation flows. The project supported by the National Natural Science Foundation of China (10502030 and 90505005).     

Global Solvability in W 2 2,1 -Weighted Spaces of the Two-Dimensional Navier–Stokes Problem in Domains with Strip-Like Outlets to Infinity

The time-dependent Navier–Stokes system is studied in a two-dimensional domain with strip-like outlets to infinity in weighted Sobolev function spaces. It is proved that under natural compatibility conditions there exists a unique solution with prescribed fluxes over cross-sections of outlets to infinity which tends in each outlet to the corresponding time-dependent Poiseuille flow. The obtained results are proved for arbitrary large norms of the data (in particular, for arbitrary fluxes) and globally in time. The authors are supported by EC FP6 MC–ToK programme SPADE2, MTKD–CT–2004–014508.     

Measurement and chemical kinetic model predictions of detonation cell size in methanol–oxygen mixtures

In this study, detonation cell sizes of methanol–oxygen mixtures are experimentally measured at different initial pressures and compositions. Good agreement is found between the experiment data and predictions based on the chemical length scales obtained from a detailed chemical kinetic model. To assess the detonation sensitivity in methanol–oxygen mixtures, the results are compared with those of hydrogen–oxygen and methane–oxygen mixtures. Based on the cell size comparison, it is shown that methanol–oxygen is more detonation sensitive than methane–oxygen but less sensitive than hydrogen–oxygen.     

On Lyapunov stability of impulsive Takagi–Sugeno fuzzy systems

We analyze the Lyapunov stability of impulsive Takagi–Sugeno fuzzy systems. Using the direct Lyapunov method, we establish sufficient conditions for the stability of these systems. We show that these conditions can be expressed in terms of a system of linear matrix inequalities. As an example, we consider an impulsive fuzzy control in a two-species “predator–prey” model. Translated from Neliniini Kolyvannya, Vol. 11, No. 4, pp. 481–494, October–December, 2008.     

The Low Mach Number Limit for the Full Navier–Stokes–Fourier System

We study the low Mach number asymptotic limit for solutions to the full Navier–Stokes–Fourier system, supplemented with ill-prepared data and considered on an arbitrary time interval. Convergencetowards the incompressible Navier–Stokes equations is shown.     

Elastic properties of Al<Subscript>2</Subscript>O<Subscript>3</Subscript>–NiAl: a modified version of Hashin–Shtrikman bounds

The modified nonlinear relations for the estimation of elastic constants of Al₂O₃–NiAl composite material are developed. The concept of microstructure and interconnectivity of phases at the interface is used. Hashin–Shtrikman relations are described in their actual form and modified version of Hashin–Shtrikman relations for bulk and shear moduli are discussed. These relations for elastic and mechanical properties are applied mainly for Al₂O₃–NiAl composite material. Theoretical predictions using modified relations are compared with Hashin–Shtrikman bounds and experimental results of elastic properties for Al₂O₃–NiAl matrix-inclusion-based composite. It is found that the predicted values of elastic and mechanical properties using modified relations are quite close to the experimental results.     

Conformational evolution under steady shear flow: a comparison between cyclic and linear block copolymers

Cyclic block copolymer is a special type of block copolymer having no free ends. Comparison of the dynamic behavior between cyclic and linear block copolymers enables an understanding of the role of chain ends in dynamics of the latter. In relation to this point, analysis was made on the conformational dynamics for a cyclic bead-spring type diblock copolymer chain, AoB, under the steady shear flow. Further comparison was made on the conformational behavior of the AoB chain and that of two symmetric linear triblock copolymer chains, A–B–A and B–A–B. For these chains, the mobility was set to be higher for the A segments than the B segments. Thus, for the AoB chain under the steady shear flow, the segments of the A block exhibit less orientational anisotropy than those of the B block. This orientational contrast is enhanced for the A–B–A chain partly because the constraint for the motion of the segments is less near the chain ends than near the chain center. Nevertheless, for the B–A–B chain, the segmental orientation over the A block becomes more anisotropic than that over the B block. Detailed analysis shows that this result is attributable to a high orientational correlation for the segments of two end B blocks, in particularly for those near the block junctions. The correlated B segments exert a tensile force on the A block thereby significantly enhancing the orientational anisotropy of the A segments.     

Aspect-ratio dependence of transient Taylor vortices close to threshold

We perform a detailed numerical study of transient Taylor vortices arising from the instability of cylindrical Couette flow with the exterior cylinder at rest for radius ratio η = 0.5 and variable aspect ratio Γ. The result of Abshagen et al. (J Fluid Mech 476:335–343, 2003) that onset transients apparently evolve on a much smaller time–scale than decay transients is recovered. It is shown to be an artefact of time scale estimations based on the Stuart–Landau amplitude equation which assumes frozen space dependence while full space–time dependence embedded in the Ginzburg–Landau formalism needs to be taken into account to understand transients already at moderate aspect ratio. Sub-critical pattern induction is shown to explain the apparently anomalous behaviour of the system at onset while decay follows the Stuart–Landau prediction more closely. The dependence of time scales on boundary effects is studied for a wide range of aspect ratios, including non-integer ones, showing general agreement with the Ginzburg–Landau picture able to account for solutions modulated by Ekman pumping at the disks bounding the cylinders.      

Non-Ergodicity of the Nosé–Hoover Thermostatted Harmonic Oscillator

The Nosé–Hoover thermostat is a deterministic dynamical system designed for computing phase space integrals for the canonical Gibbs distribution. Newton’s equations are modified by coupling an additional reservoir variable to the physical variables. The correct sampling of the phase space according to the Gibbs measure is dependent on the Nosé–Hoover dynamics being ergodic. Hoover presented numerical experiments to show that the Nosé–Hoover dynamics are non-ergodic when applied to the harmonic oscillator. In this article, we prove that the Nosé–Hoover thermostat does not give an ergodynamical system for the one- dimensional harmonic oscillator when the “mass” of the reservoir is large. Our proof of non-ergodicity uses KAM theory to demonstrate the existence of invariant tori for the Nosé–Hoover dynamical system that separate phase space into invariant regions. We present numerical experiments motivated by our analysis that seem to show that the dynamical system is not ergodic even for a moderate thermostat mass.