Dynamics and Control of Initialized Fractional-Order Systems

Due to the importance of historical effects in fractional-order systems,this paper presents a general fractional-order system and control theorythat includes the time-varying initialization response. Previous studieshave not properly accounted for these historical effects. Theinitialization response, along with the forced response, forfractional-order systems is determined. The scalar fractional-orderimpulse response is determined, and is a generalization of theexponential function. Stability properties of fractional-order systemsare presented in the complex w-plane, which is a transformation of thes-plane. Time responses are discussed with respect to pole positions inthe complex w-plane and frequency response behavior is included. Afractional-order vector space representation, which is a generalizationof the state space concept, is presented including the initializationresponse. Control methods for vector representations of initializedfractional-order systems are shown. Finally, the fractional-orderdifferintegral is generalized to continuous order-distributions whichhave the possibility of including all fractional orders in a transferfunction.     

Computation of Fractional Order Derivative and Integral via Power Series Expansion and Signal Modelling

The three techniques of s-to-z transform, power series expansion (PSE) and signal modelling are combined to develop a new procedure for efficiently computing the fractional order derivatives and integrals of discrete-time signals. A mapping function between the s-plane and the z-plane is first chosen, and then a PSE of this mapping function raised to fractional order is performed to get the desired infinite impulse response of the ideal digital fractional operator. Finally, the desired impulse response is modelled as the impulse response of a linear invariant system whose rational transfer function is determined using deterministic signal modelling techniques. Three non-iterative techniques, namely Padé, Prony and Shanks’ methods have been considered in this paper. Using Al-Alaoui’s rule as s-to-z transform, computation examples show that both Prony and Shanks’ method can achieve more accurate fractional differentiation and integration than Padé method which is equivalent to continued fraction expansion technique.     

Exact calculation of the flow of a staggered tandem cascade with moving second blade row

Summary  The plane flow around a tandem cascade of flat plates is calculated by means of conformal mapping. The blades of the two rows are perpendicular to each other. The first row is stationary, the second row moves with constant velocity. The conformal mapping will be constructed by a “mapping flow”. The blades of one row are stream lines and those of the other row are potential lines of the flow. By conformal mapping, the physical flow around the tandem cascade of the physical ζ-plane is converted into a flow between infinitely long straight walls in the z-plane, each wall corresponding to one of the blades. The conditions far upstream and far downstream of the cascade are represented by source-vortices. In the z-plane, the boundary conditions may be easily fulfilled by reflection and repetition of the source-vortices, and the flow may be calculated by well-known methods. The physical flow searched for is then obtained by inverse mapping. Received 24 July 2000; accepted for publication 6 December 2000     

On the Computation of the Third Order Terms of the Series Defining the Center Manifold for a Scalar Delay Differential Equation

When computing the third order terms of the series of powers of the function whose graph is the center manifold, at an equilibrium point of a scalar delay differential equation with a single constant delay r > 0, some problems occur at the term w_2,1z²[`(z)].{w_{2,1}z^2\overline{z}.} More precisely, in order to determine the values at 0, respectively −r of the function w _2,1(.), an algebraic system of equations must be solved. We show that the two equations are dependent, hence the system has an infinity of solutions. Then we show how we can overcome this lack of uniqueness and provide a formula for w _2,1(0).     

Three-dimensional wake structure measurement using a modified PIV technique

Three-dimensional vortical structures have been measured in a circular-cylinder wake using particle imaging velocimetry (PIV) for the Reynolds number range of 2×10³ to 1×10⁴. The PIV was modified, compared with the conventional one, in terms of its light sheet arrangement to capture reliably streamwise vortices. While in agreement with previous reports, the presently measured spanwise structures complement the data in the literature in the streamwise evolution of the near-wake spanwise vortex in size, strength, streamwise and lateral convection velocities, shedding new light upon vigorous interactions between oppositely signed spanwise structures. The longitudinal vortices display mushroom patterns in the (x, z)-plane in the immediate proximity to the cylinder. Their most likely inclination in the (x, y)-plane is inferred from the measurements in different (x, z)-planes. The longitudinal vortices in the (y, z)-plane show alternate change in sign, though not discernible at x/d > 15. They decay in the maximum vorticity and circulation rapidly from x/d = 5 to 10 and slowly for x/d > 10, and are further compared with the spanwise vortices in size, strength and rate of decay.     

Statistical mechanics of temporary polymer networks I. The equilibrium theory

In part I of this work (the present article) the equilibrium state of temporary polymer networks is treated in the framework of thermodynamics and statistical mechanics. The network is described as an open system. Thereby we use a modified spring-bead model in which the beads represent junctions that decay and reform thus adding a viscous component to the assumed elastic behaviour of the permanent network. The relevant statistical equation — analogous to Liouville's equation — is solved. The grand-canonical probability density function and two of three equations of state are derived. Explicit formulae are given for several relevant probabilities. For instance the probabilityw (z)dz that a network chain connecting two junctions has a contour length betweenz andz +dz is given by the Wien type formulaw(z) =A z ³ exp {–B z} whereA andB do not depend onz.     

A numerical study of various rheological polydispersity measures

Model calculations were performed in order to investigate the sensitivity of various rheological polydispersity parameters for variations in the moments of the molar mass distribution (MMD) of linear polymers. Molar mass distributions were generated with the Gaussian and the Generalised exponential distribution functions, using a fixed weight average molar mass M _w and variable M _w /M _n and M _z /M _n . Assuming linear entangled polymeric chains, the linear viscoelastic properties were predicted by calculating the stress relaxation modulus of the consecutive monodisperse fractions with the BSW relaxation time spectrum and blending these curves with the double reptation blending rule. BSW relaxation parameters appropriate for polypropylene were used.  It was found that both the zero-shear viscosity and the so-called cross-over frequency, at which and are equal, depend mostly on M _w but also significantly on both M _w /M _n and M _z /M _w . By contrast, the steady-state compliance depends mainly on M _z /M _w , its functional dependence on moments of the MMD being best described by the Ferry equation.  None of the polydispersity parameters PI (from the modulus cross-over), MODSEP (the modulus separation) or PDR (from the shape of the flow curve), as introduced in literature depends solely on the polydispersity M _w /M _n . PI is the most sensitive indicator for this purpose. Finally, the parameters ER ( at a fixed low value of , MODSEP en DRI (from the shape of the flow curve) are shown to be good indicators for the weight (M _z /M _w ) of the high molar mass tail of the molar mass distribution. Received: 5 May 1998 Accepted: 30 July 1998     

Solution for the thermoelastic problem in vertically inhomogeneous media

The fundamental transient-thermoelastic problem with body forces and a heat source in vertically inhomogeneous media is investigated by a method presented in this paper. The basic equations in Fourier transforms and Laplace transform are obtained in the form of two sets of first order linear ordinary differential equations inz, Eq. (7). Furthermore, forN-layered media, the general solution in the transformed spaces of thej-th layer is given for fully connected interface between layers, Eq. (11). Finally, under general condition, a closed-form solution for the quasi-static transient displacements, stresses, temperature in the body can be obtained by the convolution theorems for the two integral transforms. In the final solution, the Green's functions can be expressed in terms of Hankel transforms of order zero and unity as well as inverse Laplace transform, and come out rather neatly. Comprehensive Institute of Geotechnical Investigation and Surveying, Ministry of Urban and Rural Construction and Environmental Protection     

A survey of exact solutions of inviscid field equations in the theory of shear flow instabilty

This paper presents a survey of undisturbed flows that take one or another of the field equations of inviscid shear flow instability theory (e.g. theRayleigh equation,Taylor-Goldstein-Haurwitz equation or theKuo equation) to a differential equation satisfied by aknown transcendental function forarbitrary complex values of the parameters. Some mean velocity profiles having this feature are well known. Thus, piecewise linear mean velocity profiles take theRayleigh equations to a constant-coefficient differential equation and the exponential mean velocity profile takes theRayleigh equation to theGauss hypergeometric equation. Less well known is the fact that a variety of mean velocity profiles take theRayleigh equation to a differential equation due toKarl Heun. These profiles include: (i) the sinusoidal profile; (ii) the hyperbolic tangent profile (an example pointed out byMiles (1963)); (iii) the profile in the form of the square of a hyperbolic secant (theBickley jet); and, (iv) a skewed velocity profile in which each component has the form of a quadratic function of the variable exp(–z/l) (in whichz is the cross-stream coordinate andl is a length scale). In all of these cases, one or another author has previously identified aregular neutral mode solution of theRayleigh equation and has expressed that solution in the form of elementary functions. Such regular neutral modes apparently represent cases in which the solution ofHeun's equation (which is normally an infinite series) truncates to a single term. The survey concludes by noting that the parabolic mean velocity profile takes theRayleigh equation to thedifferential equation of the spheroidal wave function.Dedicated to Mårten T. Landahl on the occasion of his sixty-fifth birthday by a former apprentice as a token of his respect and gratitude.     

On the generalized Cauchy function and new Conjecture on its exterior singularities

This article studies on Cauchy’s function f (z) and its integral, (2πi)J[ f (z)] ≡■C f (t)dt/(t z) taken along a closed simple contour C, in regard to their comprehensive properties over the entire z = x + iy plane consisted of the simply connected open domain D + bounded by C and the open domain D outside C. (1) With f (z) assumed to be C n (n < ∞-times continuously differentiable) z ∈ D + and in a neighborhood of C, f (z) and its derivatives f (n) (z) are proved uniformly continuous in the closed domain D + = [D + + C]. (2) Cauchy’s integral formulas and their derivatives z ∈ D + (or z ∈ D ) are proved to converge uniformly in D + (or in D = [D +C]), respectively, thereby rendering the integral formulas valid over the entire z-plane. (3) The same claims (as for f (z) and J[ f (z)]) are shown extended to hold for the complement function F(z), defined to be C n z ∈ D and about C. (4) The uniform convergence theorems for f (z) and F(z) shown for arbitrary contour C are adapted to find special domains in the upper or lower half z-planes and those inside and outside the unit circle |z| = 1 such that the four general- ized Hilbert-type integral transforms are proved. (5) Further, the singularity distribution of f (z) in D is elucidated by considering the direct problem exemplified with several typ- ical singularities prescribed in D . (6) A comparative study is made between generalized integral formulas and Plemelj’s formulas on their differing basic properties. (7) Physical sig- nificances of these formulas are illustrated with applicationsto nonlinear airfoil theory. (8) Finally, an unsolved inverse problem to determine all the singularities of Cauchy function f (z) in domain D , based on the continuous numerical value of f (z) z ∈ D + = [D + + C], is presented for resolution as a conjecture.     

Thermodynamische Konsistenz zwischen Dichteanomalie und Regelation von Wasser

This paper shows the derivation of thermodynamic relations between the anomal behaviour of the density of fluid water and the anomal slope of the melting line in thep, T-plane. From these relations anomalities of the constant pressure curves in theT, s-plane between the saturation and the melting lines are revealed. Phase boundaries in theT, s-plane are calculated. A melting table from the triplepoint up to 1000 bar is included.     

On the r-value of textured sheet metals

The r-value of a sheet metal is a measure of plastic anisotropy frequently used for prediction of performance in deep-drawing. It has also figured prominently in the literature for validation of theories where the predicted angular dependence of r is compared with the measured dependence. As plastic anisotropy in sheet metals is caused mainly by the preferred orientations of grains within the polycrystalline metal, it is natural to ask how r would depend on the orientation distribution function (ODF) w which defines the crystallographic texture of the polycrystal. In this paper a general formula relating r to w is derived for textured sheet metals whose plastic flow behavior is governed by a plastic potential f(σ, w), the anisotropic part of which depends linearly on the texture coefficients; here σ denotes the deviator of the Cauchy stress. Specific forms of this formula for orthorhombic sheets of cubic and of hexagonal metals are explicitly given.     

A new analytical solution to axisymmetric Biot＇s consolidation of a finite soil layer

A new analytical method is presented to study the axisymmetric Biot＇s consolidation of a finite soil layer. Starting from the governing equations of axisymmetric Blot＇s consolidation, and based on the property of Laplace transform, the relation of basic variables for a point of a finite soil layer is established between the ground surface （z= 0） and the depth z in the Laplace and Hankel transform domains. Combined with the boundary conditions of the finite soil layer, the analytical solution of any point in the transform domain can be obtained. The actual solution in the physical domain can be obtained by inverse Laplace and Hankel transforms. A numerical analysis for the axisymmetric consolidation of a finite soil layer is carried out.     

Vibrations of a polygonal plate weakened by a circular hole with two rectilinear cuts

We consider free bending vibrations of a finite isotropic plate bounded on the outside by a polygon L ₂ (e.g., a square, an ellipse, a circle, etc.) and on the inside by the contour L ₁ formed by the circle of radius r with two rectilinear cuts located symmetrically on the axis Ox. The plate is rigidly clamped along the entire outer contour L ₂, and the inner contour L ₁ is free. We reduce the solution of the plate vibration problem to the integration of a fourth-order differential equation [2, 4, 5]. The method used in this paper is well known in the literature [2, 4, 5, 11] for simple simply connected domains, but the case under study (a doubly connected domain with cuts) has not yet been considered, because no mapping functions z = λ(ξ) have been known for complicated domains (doubly connected domain with cuts). The author [6–8] is the first in the scientific world to find such mapping functions. The obtained theoretic solution is illustrated by numerical examples.     

Closed-form solution for a mode-III interfacial edge crack between two bonded dissimilar elastic strips

A closed-form solution is obtained for the problem of a mode-III interfacial edge crack between two bonded semi-infinite dissimilar elastic strips. A general out-of-plane displacement potential for the crack interacting with a screw dislocation or a line force is constructed using conformal mapping technique and existing dislocation solutions. Based on this displacement potential, the stress intensity factor (SIF, K_III) and the energy release rate (ERR, G_III) for the interfacial edge crack are obtained explicitly. It is shown that, in the limiting special cases, the obtained results coincide with the results available in the literature. The present solution can be used as the Green’s function to analyze interfacial edge cracks subjected to arbitrary anti-plane loadings. As an example, a formula is derived correcting the beam theory used in evaluation of SIF (K_III) and ERR (G_III) of bimaterials in the double cantilever beam (DCB) test configuration.     

Evaluation of molecular weight distribution from dynamic moduli

A method to evaluate molecular weight distribution (MWD) from dynamic moduli is presented here. It relies on the least-square fitting of the dynamic data to a model whose parameters depend on the MWD. In particular, the analytical solution for the relaxation modulus previously obtained from the double reptation model, with the Tuminello step relaxation function and the Generalized Exponential Function (GEX) describing the MWD (Nobile and Cocchini 2000), has been used. A Finite Element Approximation (FEA) has been applied to calculate dynamic moduli from the relaxation modulus as a function of MWD. The sensitiveness of the GEX-double reptation dynamic moduli on the model parameters has also been investigated and the results show that large changes of the M_w/M_n ratio weakly affect the dynamic moduli, while small changes of the M_z/M_w ratio significantly deform the dynamic moduli curves. The use of rheological data to obtain MWD, by the model used in this paper, will, therefore, be able to give rather well defined M_z/M_w ratios, while more uncertainty will be presented in the M_w/M_n results. The so-called GEX-rheological model for the dynamic moduli was applied to fit the experimental data of different polymers in order to obtain the best-fit parameters of the MWD of these polymers, without the need for the inversion of the double reptation integral equation. The stability of the results has been confirmed through the evaluation of the 90% confidence intervals for the first molecular weight averages. Finally, concerning the M_w and M_z values, the predictions obtained from the dynamic moduli measurements differ by less than 10% from those obtained from GPC measurements while, as expected, more uncertainty is present in the M_n predictions. Received: 6 February 2000 Accepted: 22 August 2000     

Solution of the heat transfer of laminar forced-convection in non-circular pipes

Summary The heat transfer problems of forced-convection in non-circular pipes have many engineering applications. In this a paper a formal solution is given when the mapping function z=w()= a, which maps conformally the cross-section of the channel onto the unit circle in the -plane is known. The expression for the average velocity, average temperature, mixed mean temperature, heat transfer rate and the Nusselt number have been expressed in terms of the constants a _n .     

Stresses in a half space due to Newtonian gravitation

An efficient general solution is obtained for the problem of the elastic half space z > 0 with a traction-free surface experiencing gravitational attraction to an arbitrarily shaped body located in z < 0. Many components of the stress field can be written down immediately if the potential of the attracting body is known. Results are given for the case of attraction to a uniform sphere.     

Measurement of velocity profiles in a rectangular microchannel with aspect ratio <Emphasis Type="Italic">α</Emphasis> = 0.35

In this work, we measured 14 horizontal velocity profiles along the vertical direction of a rectangular microchannel with aspect ratio α = h/w = 0.35 (h is the height of the channel and w is the width of the channel) using microPIV at Re = 1.8 and 3.6. The experimental velocity profiles are compared with the full 3D theoretical solution, and also with a Poiseuille parabolic profile. It is shown that the experimental velocity profiles in the horizontal and vertical planes are in agreement with the theoretical profiles, except for the planes close to the wall. The discrepancies between the experimental data and 3D theoretical results in the center vertical plane are less than 3.6%. But the deviations between experimental data and Poiseuille’s results approaches 5%. It indicates that 2D Poiseuille profile is no longer a perfect theoretical approximation since α = 0.35. The experiments also reveal that, very near the hydrophilic wall (z = 0.5–1 μm), the measured velocities are significantly larger than the theoretical velocity based on the no-slip assumption. A proper discussion on some physical effects influencing the near wall velocity measurement is given.     

Thin plate bending problem of dissimilar strips with two bond lines

Summary A solution to the thin plate bending problem of partially bonded dissimilar strips with two bond lines is presented. The two strips are symmetrically bonded with respect to the interface which is on theX-axis. The complex stress functions approach together with the rational mapping function technique are used in the analysis. A concentrated bending moment applied at each strip is considered. Distributions of bending and torsional moments, as well as the stress intensity of debonding (SID) at the debonding tips are obtained, and the debonding extension is investigated.