Decay of Almost Periodic Solutions¶of Conservation Laws

We consider the asymptotic behavior of solutions of systems of inviscid or viscous conservation laws in one or several space variables, which are almost periodic in the space variables in a generalized sense introduced by Stepanoff and Wiener, which extends the original one of H. Bohr. We prove that if u(x,t) is such a solution whose inclusion intervals at time t, with respect to ?>0, satisfy l _epsiv;(t)/t→0 as t→∞, and such that the scaling sequence u ^T (x,t)=u(T x,T t) is pre-compact as t→∞ in L _loc ¹(? ^{d +1} ₊, then u(x,t) decays to its mean value \(\), which is independent of t, as t→∞. The decay considered here is in L ¹ _loc of the variable ξ≡x/t, which implies, as we show, that \(\) as t→∞, where M _x denotes taking the mean value with respect to x. In many cases we show that, if the initial data are almost periodic in the generalized sense, then so also are the solutions. We also show, in these cases, how to reduce the condition on the growth of the inclusion intervals l _?(t) with t, as t→∞, for fixed ? > 0, to a condition on the growth of l _?(0) with ?, as ?→ 0, which amounts to imposing restrictions only on the initial data. We show with a simple example the existence of almost periodic (non-periodic) functions whose inclusion intervals satisfy any prescribed growth condition as ?→ 0. The applications given here include inviscid and viscous scalar conservation laws in several space variables, some inviscid systems in chromatography and isentropic gas dynamics, as well as many viscous 2 × 2 systems such as those of nonlinear elasticity and Eulerian isentropic gas dynamics, with artificial viscosity, among others. In the case of the inviscid scalar equations and chromatography systems, the class of initial data for which decay results are proved includes, in particular, the L ^∞ generalized limit periodic functions. Our procedures can be easily adapted to provide similar results for semilinear and kinetic relaxations of systems of conservation laws.     

Static statement of the stability problem under creep

In technological processes of rod bending, the critical time is determined [1] by the criterion of unbounded increase A → ∞ in the bent axis amplitude, which is equivalent to the requirement A ? A ₀, where A ₀ is value of the amplitude at the initial time t = 0. In this case, the mathematical models of the process of buckling of rods and plates [2] are constructed in the framework of the theory of small displacements. This contradiction can be removed by the assumption that the critical state is realized for deflections A of the order of several A ₀, i.e., at the time instant corresponding to a sharp increase in displacements. Naturally, this assumption is of local character, because the instant of the transition to the accelerated increase in deflections depends on specific conditions such as, for example, the support conditions, the creep coefficient, the type of the system imperfectness, the value of A ₀, and the eccentricity of the load application.In what follows, we show that, in the case of longitudinal bending (buckling), the time instant directly preceding the beginning of the catastrophic increase in deflections can be determined by the variation in the phase volume of the system.     

Miscible Thermo-Viscous Fingering Instability in Porous Media. Part 1: Linear Stability Analysis

The development of the thermo-viscous fingering instability of miscible displacements in homogeneous porous media is examined. In this first part of the study dealing with stability analysis, the basic equations and the parameters governing the problem in a rectilinear geometry are developed. An exponential dependence of viscosity on temperature and concentration is represented by two parameters, thermal mobility ratio β _T and a solutal mobility ratio β _C , respectively. Other parameters involved are the Lewis number Le and a thermal-lag coefficient λ. The governing equations are linearized and solved to obtain instability characteristics using either a quasi-steady-state approximation (QSSA) or initial value calculations (IVC). Exact analytical solutions are also obtained for very weakly diffusing systems. Using the QSSA approach, it was found that an increase in thermal mobility ratio β _T is seen to enhance the instability for fixed β _C , Le and λ. For fixed β _C and β _T , a decrease in the thermal-lag coefficient and/or an increase in the Lewis number always decrease the instability. Moreover, strong thermal diffusion at large Le as well as enhanced redistribution of heat between the solid and fluid phases at small λ is seen to alleviate the destabilizing effects of positive β _T . Consequently, the instability gets strictly dominated by the solutal front. The linear stability analysis using IVC approach leads to conclusions similar to the QSSA approach except for the case of large Le and unity λ flow where the instability is seen to get even less pronounced than in the case of a reference isothermal flow of the same β _C , but β _T  = 0. At practically, small value of λ, however, the instability ultimately approaches that due to β _C only.     

Centers and Uniform Isochronous Centers of Planar Polynomial Differential Systems

For planar polynomial vector fields of the form

$$\begin{aligned} (-y+X(x,y))\dfrac{\partial }{\partial x}+(x+Y(x,y))\dfrac{\partial }{\partial y}, \end{aligned}$$

where X and Y start at least with terms of second order in the variables x and y, we determine necessary and sufficient conditions under which the origin is a center or a uniform isochronous centers.

     

Spray-Flame Dynamics in a Rich Droplet Array

We numerically study spray-flame dynamics. The initial state of the spray is schematized by alkane droplets located at the nodes of a face-centered 2D-lattice. The droplets are surrounded by a gaseous mixture of alkane and air. The lattice spacing s reduced by the combustion length scale is large enough to consider that the chemical reaction occurs in a heterogeneous medium. The overall spray equivalence ratio is denoted by ?_T, with ?_T = ?_L + ?_G, where ?_G corresponds to the equivalence ratio of the gaseous surrounding mixture at the initial saturated partial pressure, while ?_L is the so-called liquid loading. To model such a heterogenous combustion, the retained chemical scheme is a global irreversible one-step reaction governed by an Arrhenius law, with a modified heat of reaction depending on the local equivalence ratio. ?_T is chosen in the range 0.9 ≤ ?_T ≤ 2. Three geometries (s = 3, s = 6, s = 12) and four liquid loadings, ?_L = 0.3, ?_L = 0.5, ?_L = 0.7, ?_L = 0.85 are studied. In the rich sprays, our model qualitatively retrieves the recent experimental measurements: the rich spray-flames can propagate faster than the single-phase flames with the same overall equivalence ratio. To analyse the conditions for this enhancement, we introduce the concept of “spray Peclet number”, which compares the droplet vaporization time with the combustion propagation time of the single-phase flame spreading in the fresh surrounding mixture.     

Bistable Boundary Reactions in Two Dimensions

In a bounded domain \({\Omega \subset \mathbb R^2}\) with smooth boundary we consider the problem

$\Delta u = 0 \quad {\rm{in }}\, \Omega, \qquad \frac{\partial u}{\partial \nu} = \frac1\varepsilon f(u) \quad {\rm{on }}\,\partial\Omega,$

where ν is the unit normal exterior vector, ε > 0 is a small parameter and f is a bistable nonlinearity such as f(u) = sin(π u) or f(u) = (1 ? u ²)u. We construct solutions that develop multiple transitions from ?1 to 1 and vice-versa along a connected component of the boundary ?Ω. We also construct an explicit solution when Ω is a disk and f(u) = sin(π u).

     

Power-law creep and residual stresses in a carbopol gel

We report on the interplay between creep and residual stresses in a carbopol microgel. When a constant shear stress σ is applied below the yield stress σ _y, the strain is shown to increase as a power law of time, γ(t) = γ ₀ + (t/τ) ^α , with an exponent α = 0.39 ± 0.04 that is strongly reminiscent of Andrade creep in hard solids. For applied shear stresses lower than some typical value σ _c ? 0.2σ _y, the microgel experiences a more complex, anomalous creep behavior, characterized by an initial decrease of the strain, that we attribute to the existence of residual stresses of the order of σ _c that persist after a rest time under a zero shear rate following preshear. The influence of gel concentration on creep and residual stresses are investigated as well as possible aging effects. We discuss our results in light of previous works on colloidal glasses and other soft glassy systems.     

Operator method for solving the plane elasticity problem for a strip with periodic cuts

In the present paper, we use the conformal mapping z/c = ζ?2a sin ζ (a, c?const, ζ = u + iv) of the strip {|v| ≤ v ₀, |u| < ∞} onto the domain D, which is a strip with symmetric periodic cuts. For the domain D, in the orthogonal system of isometric coordinates u, v, we solve the plane elasticity problem. We seek the biharmonic function in the form F = C ψ ₀ + S ψ*₀ + x(C ψ ₁ ? S ψ ₂) + y(C ψ ₂ + S ψ ₁), where C(v) and S(v) are the operator functions described in [1] and ψ ₀(u), …, ψ ₂(u) are the desired functions. The boundary conditions for the function F posed for v = ±v ₀ are equivalent to two operator equations for ψ ₁(u) and ψ ₂(u) and to two ordinary differential equations of first order for ψ ₀(u) and ψ*₀(u) [2]. By finding the functions ψ _j (u) in the form of trigonometric series with indeterminate coefficients and by solving the operator equations, we obtain infinite systems of linear equations for the unknown coefficients. We present an efficient method for solving these systems, which is based on studying stable recursive relations. In the present paper, we give an example of analysis of a specific strip (a = 1/4, v ₀ = 1) loaded on the boundary v = v ₀ by a normal load of intensity p. We find the particular solutions corresponding to the extension of the strip by the longitudinal force X ^∞ and to the transverse and pure bending of the strip due to the transverse force Y ^∞ and the constant moment M ^∞, respectively. We also present the graphs of normal and tangential stresses in the transverse cross-section x = 0 and study the stress concentration effect near the cut bottom.     

The Moments Sliding Vector Field on the Intersection of Two Manifolds

In this work, we consider a special choice of sliding vector field on the intersection of two co-dimension 1 manifolds. The proposed vector field, which belongs to the class of Filippov vector fields, will be called moments vector field and we will call moments trajectory the associated solution trajectory. Our main result is to show that the moments vector field is a well defined, and smoothly varying, Filippov sliding vector field on the intersection \(\Sigma \) of two discontinuity manifolds, under general attractivity conditions of \(\Sigma \). We also examine the behavior of the moments trajectory at first order exit points, and show that it exits smoothly at these points. Numerical experiments illustrate our results and contrast the present choice with other choices of Filippov sliding vector field.     

Synchronization transitions induced by partial time delay in a excitatory–inhibitory coupled neuronal network

Information transmission delays are an inherent factor of neuronal systems as a consequence of the finite propagation speeds and time lapses occurring by both dendritic and synaptic processes. In real neuronal systems, some delay between two neurons is too small and can be ignored, which results in partial time delay. In this paper, we focus on investigating influences of partial time delay on synchronization transitions in a excitatory–inhibitory (E–I) coupled neuronal networks. Here, we suppose time delay between two neurons equals to \(\tau \) with probability \(p_{\mathrm{delay}}\) and investigate effect of partial time delay on synchronization transitions of the neuronal networks by controlling \(\tau \) and \(p_{\mathrm{delay}}\) under three cases. In these three cases, excitatory synapses are always considered to delayed with probability \(p_{\mathrm{delay}}\), while inhibitory synapses are considered to be without delays (case I), delayed with probability \(p_{\mathrm{delay}}\) (case II), and always delayed (case III), respectively. It is revealed that, in the first two cases, partial time delay has little influences on synchronization of the neuronal network for small \(p_{\mathrm{delay}}\), while it could induce synchronization transitions at \(\tau \) around integer multiples of the period of individual neuron T when \(p_{\mathrm{delay}}\) is large enough, while in the case III, partial time delay could induce synchronization transitions at \(\tau \) being around odd integer multiples of T / 2 for small \(p_{\mathrm{delay}}\) and at \(\tau \) being around integer multiples of T for large \(p_{\mathrm{delay}}\). Most interesting observation is that partial time delay could induce frequent synchronization transitions at \(\tau \) being around integer multiples of T / 2 for intermediate \(p_{\mathrm{delay}}\). Moreover, effect of rewiring probability on synchronization transitions induced by partial time delay has been discussed. It is found that synchronization transitions induced by partial time delay are robust to rewiring probability for large \(p_{\mathrm{delay}}\) under the three cases.     

Regularity and Singularities of Optimal Convex Shapes in the Plane

We focus here on the analysis of the regularity or singularity of solutions Ω ₀ to shape optimization problems among convex planar sets, namely:

$J(\Omega_{0})={\rm min} \{J(\Omega), \Omega \quad {\rm convex},\Omega \in \mathcal{S}_{\rm ad}\},$

where \({\mathcal{S}_{\rm ad}}\) is a set of 2-dimensional admissible shapes and \({J:\mathcal{S}_{\rm ad}\rightarrow\mathbb{R}}\) is a shape functional. Our main goal is to obtain qualitative properties of these optimal shapes by using first and second order optimality conditions, including the infinite dimensional Lagrange multiplier due to the convexity constraint. We prove two types of results:

1. i)
   under a suitable convexity property of the functional J, we prove that Ω ₀ is a W ^2,p -set, \({p\in[1, \infty]}\). This result applies, for instance, with p = ∞ when the shape functional can be written as J(Ω) = R(Ω) + P(Ω), where R(Ω) = F(|Ω|, E _f (Ω), λ₁(Ω)) involves the area |Ω|, the Dirichlet energy E _f (Ω) or the first eigenvalue of the Laplace–Dirichlet operator λ₁(Ω), and P(Ω) is the perimeter of Ω;
    

1. ii)
   under a suitable concavity assumption on the functional J, we prove that Ω ₀ is a polygon. This result applies, for instance, when the functional is now written as J(Ω) = R(Ω) ? P(Ω), with the same notations as above.
    

     

Steady Nearly Incompressible Vector Fields in Two-Dimension: Chain Rule and Renormalization

Given bounded vector field \({b : {\mathbb{R}^{d}} \to {\mathbb{R}^{d}}}\), scalar field \({u : {\mathbb{R}^{d}} \to {\mathbb{R}}}\), and a smooth function \({\beta : {\mathbb{R}} \to {\mathbb{R}}}\), we study the characterization of the distribution \({{\rm div}(\beta(u)b)}\) in terms of div b and div(ub). In the case of BV vector fields b (and under some further assumptions), such characterization was obtained by L. Ambrosio, C. De Lellis and J. Malý, up to an error term which is a measure concentrated on the so-called tangential set of b. We answer some questions posed in their paper concerning the properties of this term. In particular, we construct a nearly incompressible BV vector field b and a bounded function u for which this term is nonzero. For steady nearly incompressible vector fields b (and under some further assumptions), in the case when d = 2, we provide complete characterization of div(\({\beta(u)b}\)) in terms of div b and div(ub). Our approach relies on the structure of level sets of Lipschitz functions on \({{\mathbb{R}^{2}}}\) obtained by G. Alberti, S. Bianchini and G. Crippa. Extending our technique, we obtain new sufficient conditions when any bounded weak solution u of \({\partial_t u + b \cdot \nabla u=0}\) is renormalized, that is when it also solves \({\partial_t \beta(u) + b \cdot \nabla \beta(u)=0}\) for any smooth function \({\beta \colon{\mathbb{R}} \to {\mathbb{R}}}\). As a consequence, we obtain new a uniqueness result for this equation.     

Revisiting History Effects in Adverse-Pressure-Gradient Turbulent Boundary Layers

The goal of this study is to present a first step towards establishing criteria aimed at assessing whether a particular adverse-pressure-gradient (APG) turbulent boundary layer (TBL) can be considered well-behaved, i.e., whether it is independent of the inflow conditions and is exempt of numerical or experimental artifacts. To this end, we analyzed several high-quality datasets, including in-house numerical databases of APG TBLs developing over flat-plates and the suction side of a wing section, and five studies available in the literature. Due to the impact of the flow history on the particular state of the boundary layer, we developed three criteria of convergence to well-behaved conditions, to be used depending on the particular case under study. (i) In the first criterion, we develop empirical correlations defining the R e _?? -evolution of the skin-friction coefficient and the shape factor in APG TBLs with constant values of the Clauser pressure-gradient parameter β = 1 and 2 (note that β = δ ^?/τ _w dP _e /dx, where δ ^? is the displacement thickness, τ _w the wall-shear stress and dP _e /dx the streamwise pressure gradient). (ii) In the second one, we propose a predictive method to obtain the skin-friction curve corresponding to an APG TBL subjected to any streamwise evolution of β, based only on data from zero-pressure-gradient TBLs. (iii) The third method relies on the diagnostic-plot concept modified with the shape factor, which scales APG TBLs subjected to a wide range of pressure-gradient conditions. These three criteria allow to ensure the correct flow development of a particular TBL, and thus to separate history and pressure-gradient effects in the analysis.     

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall considering the interaction between squirmer and wall are numerically studied with an immersed boundary-lattice Boltzmann method. The power-law index, Reynolds number, initial orientation angle of squirmer, and initial distance of squirmer from the wall are all taken into account to investigate the swimming characteristics for pusher (β?<?0), neutral squirmer (β?=?0), and puller (β?>?0) (three kinds of swimmer types) near the no-slip boundary. Four new kinds of swimming modes are found. Results show that, for the pushers and pullers, the wall displays an increasing attraction with increasing power-law index n, which differs from the neutral squirmer who always departs from the wall after the first collision with the wall. Both the initial orientation angle and initial distance from the wall only affect the moving situations rather than the moving modes of the squirmers. However, the squirmers depart from the wall as the Reynolds number increases and chaotic orbits appear for some squirmers at Re?=?5. Several typical flow fields are analyzed and the power consumption and torque for different kinds of flows are also studied. It is found that, as the absolute value of β increases, the power consumption generally increases in shear-thinning (n?=?0.4), Newtonian (n?=?1), and shear-thickening (n?=?1.6) fluids. Moreover, the pushers (β?<?0) and the pullers (β?>?0) expend almost the same power if the absolute value of β remains the same. In addition, the power consumption of the squirmers is highly dependent on the power-law index n.     

Does Density Ratio Significantly Affect Turbulent Flame Speed?

In order to experimentally study whether or not the density ratio σ substantially affects flame displacement speed at low and moderate turbulent intensities, two stoichiometric methane/oxygen/nitrogen mixtures characterized by the same laminar flame speed S _L = 0.36 m/s, but substantially different σ were designed using (i) preheating from T _u = 298 to 423 K in order to increase S _L , but to decrease σ, and (ii) dilution with nitrogen in order to further decrease σ and to reduce S _L back to the initial value. As a result, the density ratio was reduced from 7.52 to 4.95. In both reference and preheated/diluted cases, direct images of statistically spherical laminar and turbulent flames that expanded after spark ignition in the center of a large 3D cruciform burner were recorded and processed in order to evaluate the mean flame radius \(\bar {R}_{f}\left (t \right )\) and flame displacement speed \(S_{t}=\sigma ^{-1}{d\bar {R}_{f}} \left / \right . {dt}\) with respect to unburned gas. The use of two counter-rotating fans and perforated plates for near-isotropic turbulence generation allowed us to vary the rms turbulent velocity \(u^{\prime }\) by changing the fan frequency. In this study, \(u^{\prime }\) was varied from 0.14 to 1.39 m/s. For each set of initial conditions (two different mixture compositions, two different temperatures T _u , and six different \(u^{\prime })\), five (respectively, three) statistically equivalent runs were performed in turbulent (respectively, laminar) environment. The obtained experimental data do not show any significant effect of the density ratio on S _t . Moreover, the flame displacement speeds measured at u′/S _L = 0.4 are close to the laminar flame speeds in all investigated cases. These results imply, in particular, a minor effect of the density ratio on flame displacement speed in spark ignition engines and support simulations of the engine combustion using models that (i) do not allow for effects of the density ratio on S _t and (ii) have been validated against experimental data obtained under the room conditions, i.e. at higher σ.     

Planar Vector Fields with a Given Set of Orbits

We determine all the \({\mathcal{C}^1}\) planar vector fields with a given set of orbits of the form y ? y(x) = 0 satisfying convenient assumptions. The case when these orbits are branches of an algebraic curve is also study. We show that if a quadratic vector field admits a unique irreducible invariant algebraic curve \({g(x, y) = \sum_{j=0}^S a_j(x) y^{S-j}= 0}\) with S branches with respect to the variable y, then the degree of the polynomial g is at most 4S.     

Error estimate of the modified point-mass trajectory model of an artillery shell

The article deals with the motion of an axially symmetric spinning artillery shell in the gravity field under the action of the system of aerodynamic forces and moments adopted in ballistics. As the starting point, the system of differential equations of motion of the shell is taken, which is obtained from the original “accurate” system by its linearization in the variables describing the angular motion of the symmetry axis and by additional linearization in the angle between the velocity vector of the center of mass and the vertical plane (l-system). This article examines the system of differential equations of the translational motion and axial rotation of the shell which describes its modified point-mass trajectory model as applied to l-system (m-system). By small parameter methods, an estimate is obtained for the difference of the solution of l-system with given initial data and the solution of m-system with the same initial data for the variables of translational motion and axial rotation. This analytical evaluation is built in such a way that it corresponds with certain numerical estimates for components of the translational motion and axial rotation. It is observed that, under accepted assumptions, m-system and l-system determine the translational motion of the shell with the same order of the error as compared to the original “accurate” nonlinear system of equations of motion of the shell. But m-system does not contain rapidly oscillating variables describing the angular motion of the symmetry axis, and so its numerical integration requires tens of times less computational resources than the numerical integration of l-system. Numerical simulation data are represented.     

A simple criterion for finite time stability with application to impacted buckling of elastic columns

The existing theories of finite-time stability depend on a prescribed bound on initial disturbances and a prescribed threshold for allowable responses. It remains a challenge to identify the critical value of loading parameter for finite time instability observed in experiments without the need of specifying any prescribed threshold for allowable responses. Based on an energy balance analysis of a simple dynamic system, this paper proposes a general criterion for finite time stability which indicates that finite time stability of a linear dynamic system with constant coefficients during a given time interval [0, t _f ] is guaranteed provided the product of its maximum growth rate (determined by the maximum eigen-root p₁ >0) and the duration t _f does not exceed 2, i.e., p₁t _f <2. The proposed criterion (p₁t _f =2) is applied to several problems of impacted buckling of elastic columns: (i) an elastic column impacted by a striking mass, (ii) longitudinal impact of an elastic column on a rigid wall, and (iii) an elastic column compressed at a constant speed (“Hoff problem”), in which the time-varying axial force is replaced approximately by its average value over the time duration. Comparison of critical parameters predicted by the proposed criterion with available experimental and simulation data shows that the proposed criterion is in robust reasonable agreement with the known data, which suggests that the proposed simple criterion (p₁t _f =2) can be used to estimate critical parameters for finite time stability of dynamic systems governed by linear equations with constant coefficients.     

Linearized Stability for Semiflows Generated by a Class of Neutral Equations,with Applications to State-Dependent Delays

We prove a principle of linearized stability for semiflows generated by neutral functional differential equations of the form x′(t) = g(? x _t , x _t ). The state space is a closed subset in a manifold of C ²-functions. Applications include equations with state-dependent delay, as for example x′(t) = a x′(t + d(x(t))) + f (x(t + r(x(t)))) with \({a\in\mathbb{R}, d:\mathbb{R}\to(-h,0), f:\mathbb{R}\to\mathbb{R}, r:\mathbb{R}\to[-h,0]}\).     

Effects of Fuel Lewis Number on Localised Forced Ignition of Globally Stoichiometric Stratified Mixtures: a Numerical Investigation

The influences of fuel Lewis number Le_F on localised forced ignition of globally stoichiometric stratified mixtures have been analysed using three-dimensional compressible Direct Numerical Simulations (DNS) for cases with Le_F ranging from 0.8 to 1.2. The globally stoichiometric stratified mixtures with different values of root-mean-square (rms) equivalence ratio fluctuation (i.e. ?^′= 0.2, 0.4 and 0.6) and the Taylor micro-scale l_? of equivalence ratio ? variation (i.e. l_?/l_f= 2.1, 5.5 and 8.3 with l_f being the Zel’dovich flame thickness of the stoichiometric laminar premixed flame) have been considered for different initial rms values of turbulent velocity u^′. A pseudo-spectral method is used to initialise the equivalence ratio variation following a presumed bi-modal distribution for prescribed values of ?^′ and l_?/l_f for global mean equivalence ratio 〈?〉=1.0. The localised ignition is accounted for by a source term in the energy transport equation that deposits energy for a stipulated time interval. It has been observed that the maximum values of temperature and the fuel reaction rate magnitude increase with decreasing Le_F during the period of external energy deposition. The initial values of Le_F, u^′/S_b(?=1), ?^′ and l_?/l_f have been found to have significant effects on the extent of burning of the stratified mixtures following localised ignition. For a given value of u^′/S_b(?=1), the extent of burning decreases with increasing Le_F. An increase in u^′ leads to a monotonic reduction in the burned gas mass for all values of Le_F in all stratified mixture cases but an opposite trend is observed for the Le_F=0.8 homogeneous mixture. It has been found that an increase in ?^′ has adverse effects on the burned gas mass, whereas the effects of l_?/l_f on the extent of burning are non-monotonic and dependent on ?^′ and Le_F. Detailed physical explanations have been provided for the observed Le_F, u^′/S_b(?=1), ?^′ and l_?/l_f dependences.