On the unsteady rotational flow of a fractional second grade fluid through a circular cylinder

Here the velocity field and the associated tangential stress corresponding to the rotational flow of a generalized second grade fluid within an infinite circular cylinder are determined by means of the Laplace and finite Hankel transforms. At time t=0 the fluid is at rest and the motion is produced by the rotation of the cylinder around its axis. The solutions that have been obtained are presented under series form in terms of the generalized G-functions. The similar solutions for ordinary second grade and Newtonian fluids are obtained from general solution for β→1, respectively, β→1 and α ₁→0. Finally, the influences of the pertinent parameters on the fluid motion, as well as a comparison between models, is underlined by graphical illustrations.     

A unified approach to closed-form solutions of moving heat-source problems

A closed-form model for the computation of temperature distribution in an infinitely extended isotropic body with a time-dependent moving-heat sources is discussed. The temperature solutions are presented for the sources of the forms: (i) ₀₁(t)=₀ exp(−λt), (ii) ₀₂(t) =₀(t/t ^*)exp(−λt), and ₀₃(t)=₀[1+a cos(ωt)], where λ and ω are real parameters and t ^* characterizes the limiting time. The reduced (or dimensionless) temperature solutions are presented in terms of the generalized representation of an incomplete gamma function Γ(α,x;b) and its decomposition C _Γ and S _Γ. The solutions are presented for moving, -point, -line, and -plane heat sources. It is also demonstrated that the present analysis covers the classical temperature solutions of a constant strength source under quasi-steady state situations. Received on 13 June 1997     

On the unsteady rotational flow of a generalized Maxwell fluid through a circular cylinder

In this paper, the velocity field and the associated tangential stress corresponding to the rotational flow of a generalized Maxwell fluid within an infinite circular cylinder are determined by means of the Laplace and finite Hankel transforms. Initially, the fluid is at rest, and the motion is produced by the rotation of the cylinder about its axis with a unsteady angular velocity. The solutions that have been obtained are presented under series form in terms of the generalized G _a,b,c (, t)-functions. The similar solutions for the ordinary Maxwell and Newtonian fluids, performing the same motion, are obtained as special cases, when β → 1, respectively β → 1 and λ → 0, from general solutions. Finally, the solutions that have been obtained are compared by graphical illustrations, and the influence of the pertinent parameters on the fluid motion is also underlined by graphical illustrations.     

Exact analytical solutions for axial flow of a fractional second grade fluid between two coaxial cylinders

The velocity field and the adequate shear stress corresponding to the longitudinal flow of a fractional second grade fluid, between two infinite coaxial circular cylinders, are determined by applying the Laplace and finite Hankel transforms. Initially the fluid is at rest, and at time t = 0⁺, the inner cylinder suddenly begins to translate along the common axis with constant acceleration. The solutions that have been obtained are presented in terms of generalized G functions. Moreover, these solutions satisfy both the governing differential equations and all imposed initial and boundary conditions. The corresponding solutions for ordinary second grade and Newtonian fluids are obtained as limiting cases of the general solutions. Finally, some characteristics of the motion, as well as the influences of the material and fractional parameters on the fluid motion and a comparison between models, are underlined by graphical illustrations.     

Some exact solutions of the oscillatory motion of a generalized second grade fluid in an annular region of two cylinders

The velocity field and the associated shear stress corresponding to the longitudinal oscillatory flow of a generalized second grade fluid, between two infinite coaxial circular cylinders, are determined by means of the Laplace and Hankel transforms. Initially, the fluid and cylinders are at rest and at t = 0＋ both cylinders suddenly begin to oscillate along their common axis with simple harmonic motions having angular frequencies Ω1 and Ω2. The solutions that have been obtained are presented under integral and series forms in terms of the generalized G and R functions and satisfy the governing differential equation and all imposed initial and boundary conditions. The respective solutions for the motion between the cylinders, when one of them is at rest, can be obtained from our general solutions. Furthermore, the corresponding solutions for the similar flow of ordinary second grade fluid and Newtonian fluid are also obtained as limiting cases of our general solutions. At the end, the effect of different parameters on the flow of ordinary second grade and generalized second grade fluid are investigated graphically by plotting velocity profiles.     

RETRACTED ARTICLE: Flow of fractional Maxwell fluid between coaxial cylinders

This paper deals with the study of unsteady flow of a Maxwell fluid with fractional derivative model, between two infinite coaxial circular cylinders, using Laplace and finite Hankel transforms. The motion of the fluid is produced by the inner cylinder that, at time t = 0⁺, is subject to a time-dependent longitudinal shear stress. Velocity field and the adequate shear stress are presented under series form in terms of the generalized G and R functions. The solutions that have been obtained satisfy all imposed initial and boundary conditions. The corresponding solutions for ordinary Maxwell and Newtonian fluids are obtained as limiting cases of general solutions. Finally, the influence of the pertinent parameters on the fluid motion as well as a comparison between the three models is underlined by graphical illustrations.     

On the unsteady rotational flow of fractional Oldroyd-B fluid in cylindrical domains

This paper concerned with the unsteady rotational flow of fractional Oldroyd-B fluid, between two infinite coaxial circular cylinders. To solve the problem we used the finite Hankel and Laplace transforms. The motion is produced by the inner cylinder that, at time t=0⁺, is subject to a time-dependent rotational shear. The solutions that have been obtained, presented under series form in terms of the generalized G functions, satisfy all imposed initial and boundary conditions. The corresponding solutions for ordinary Oldroyd-B, fractional and ordinary Maxwell, fractional and ordinary second grade, and Newtonian fluids, performing the same motion, are obtained as limiting cases of general solutions.     

Wave Packets, Signaling and Resonances in a Homogeneous Waveguide

We solve the initial-boundary-value linear stability problem for small localised disturbances in a homogeneous elastic waveguide formally by applying a combined Laplace – Fourier transform. An asymptotic evaluation of the solution, expressed as an inverse Laplace – Fourier integral, is carried out by means of the mathematical formalism of absolute and convective instabilities. Wave packets, triggered by perturbations localised in space and finite in time, as well as responses to sources localised in space, with the time dependence satisfying e^−iωt + O(e^−ɛt ), for t → ∞, where Im ω₀ = 0 and ω > 0 , that is, the signaling problem, are treated. For this purpose, we analyse the dispersion relation of the problem analytically, and by solving numerically the eigenvalue stability problem. It is shown that due to double roots in a wavenumber k of the dispersion relation function D(k, ω), for real frequencies ω, that satisfy a collision criterion, wave packets with an algebraic temporal decay and signaling with an algebraic temporal growth, that is, temporal resonances, are present in a neutrally stable homogeneous waveguide. Moreover, for any admissible combination of the physical parameters, a homogeneous waveguide possesses a countable set of temporally resonant frequencies. Consequences of these results for modelling in seismology are discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

Scale-adaptive simulation of flow past wavy cylinders at a subcritical Reynolds number

The Xu & Yan scale-adaptive simulation (XYSAS) model is employed to simulate the flows past wavy cylinders at Reynolds number 8 × 10 3.This approach yields results in good agreement with experimental measurements.The mean flow field and near wake vortex structure are replicated and compared with that of a corresponding circular cylinder.The effects of wavelength ratios λ/D m from 3 to 7,together with the amplitude ratios a /D m of 0.091 and 0.25,are fully investigated.Owing to the wavy configuration,a maximum reduction of Strouhal number and root-meansquare (r.m.s) fluctuating lift coefficients are up to 50% and 92%,respectively,which means the vortex induced vibration (VIV) could be effectively alleviated at certain larger values of λ/D m and a /D m.Also,the drag coefficients can be reduced by 30%.It is found that the flow field presents contrary patterns with the increase of λ/D m.The free shear layer becomes much more stable and rolls up into mature vortex only further downstream when λ/D m falls in the range of 5-7.The amplitude ratio a /D m greatly changes the separation line,and subsequently influences the wake structures.     

Vorticities in a LES Model for 3D Periodic Turbulent Flows

This paper is devoted to the study of a LES model to simulate turbulent 3D periodic flow. We focus our attention on the vorticity equation derived from this LES model for small values of the numerical grid size δ. We obtain entropy inequalities for the sequence of corresponding vorticities and corresponding pressures independent of δ, provided the initial velocity u₀ is in L_x² while the initial vorticity ω₀ = ∇ × u₀ is in L_x¹. When δ tends to zero, we show convergence, in a distributional sense, of the corresponding equations for the vorticities to the classical 3D equation for the vorticity.     

Analysis of regular and chaotic dynamics of the Euler-Bernoulli beams using finite difference and finite element methods

Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions(symmetric and non-symmetric)are studied in this work.Reliability of the obtained results is verified by the finite difference method(FDM)and the finite element method(FEM)with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes(regular and non-regular).The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.     

Relaxation patterns of long,linear, flexible,monodisperse polymers: BSW spectrum revisited

Theoretical predictions for the dynamic moduli of long, linear, flexible, monodisperse polymers are summarized and compared with experimental observations. Surprisingly, the predicted 1/2 power scaling of the long-time modes of the relaxation spectrum is not found in the experiments. Instead, scaling with a power of about 1/4 extends all the way up to the longest relaxation times near τ/τ _max = 1. This is expressed in the empirical relaxation time spectrum of Baumgaertel-Schausberger-Winter, denoted as “BSW spectrum,” and justifies a closer look at the properties of the BSW spectrum. Working with the BSW spectrum, however, is made difficult by the fact that hypergeometric functions occur naturally in BSW-based rheological material functions. BSW provides no explicit solutions for the dynamic moduli, G ^′(ω), G ^″(ω), or the relaxation modulus G(t). To overcome this problem, close approximations of simple analytical form are shown for these moduli. With these approximations, analysis of linear viscoelastic data allows the direct determination of BSW parameters.     

Estimating the viscoelastic moduli of complex fluids using the generalized Stokes–Einstein equation

We obtain the linear viscoelastic shear moduli of complex fluids from the time-dependent mean square displacement, <Δr ²(t)>, of thermally-driven colloidal spheres suspended in the fluid using a generalized Stokes–Einstein (GSE) equation. Different representations of the GSE equation can be used to obtain the viscoelastic spectrum, G˜(s), in the Laplace frequency domain, the complex shear modulus, G ^*(ω), in the Fourier frequency domain, and the stress relaxation modulus, G _r (t), in the time domain. Because trapezoid integration (s domain) or the Fast Fourier Transform (ω domain) of <Δr ²(t)> known only over a finite temporal interval can lead to errors which result in unphysical behavior of the moduli near the frequency extremes, we estimate the transforms algebraically by describing <Δr ²(t)> as a local power law. If the logarithmic slope of <Δr ²(t)> can be accurately determined, these estimates generally perform well at the frequency extremes. Received: 8 September 2000/Accepted: 9 March 2000     

On resolution to Wu＇s conjecture on Cauchy function＇s exterior singularities

This is a series of studies on Wu’s conjecture and on its resolution to be presented herein. Both are devoted to expound all the comprehensive properties of Cauchy’s function f(z) (z = x + iy) and its integral J[f(z) ] ≡(2πi) -1 C f(t)(t z) -1dt taken along the unit circle as contour C,inside which(the open domain D+) f(z) is regular but has singularities distributed in open domain Doutside C. Resolution is given to the inverse problem that the singularities of f(z) can be determined in analytical form in terms of the values f(t) of f(z) numerically prescribed on C(|t| = 1) ,as so enunciated by Wu’s conjecture. The case of a single singularity is solved using complex algebra and analysis to acquire the solution structure for a standard reference. Multiple singularities are resolved by reducing them to a single one by elimination in principle,for which purpose a general asymptotic method is developed here for resolution to the conjecture by induction,and essential singularities are treated with employing the generalized Hilbert transforms. These new methods are applicable to relevant problems in mathematics,engineering and technology in analogy with resolving the inverse problem presented here.     

The crumbling of a viscous prism with an inclined free surface

Summary We study the two-dimensional instantaneous Stokes flow driven by gravity in a viscous triangular prism supported by a horizontal rigid substrate and a vertical wall. The oblique side of the prism, inclined at an angle α with respect to the substrate, is a fluid-air interface, where the stresses are zero and surface tension is neglected. We develop the stream function ψ in polar coordinates (r,θ) centered at the vertex of α and split it into an inhomogeneous part, which accounts for gravity effects, and a homogeneous part, which is expressed as a series expansion. The inhomogeneous part and the first term of the expansion may be envisioned, respectively, as self-similar solutions of the first kind and of the second kind for r→0, each one holding in complementary α domains with a crossover at α _c =21.47^∘, which we study in some detail. The coefficients of the series are calculated by truncating the expansion and using the method of direct collocation with a suitable set of points at the wall. The solution strictly holds for t=0, because later the free surface ceases to be a plane; nevertheless, it provides a good approximation for the early time evolution of the fluid profile, as shown by the comparison with numerical simulations. For 0<α<45^∘, the vertex angle remains constant and the edge remains strictly at rest; the transition at α _c manifests itself through a change in the rate of growth of the curvature. The time at which the edge starts to move (waiting time) cannot be calculated since the instantaneous solution ceases to be valid. For α>45^∘, the instantaneous local solution indicates that the surface near the vertex is launched against the substrate so that the edge starts to move immediately with a power law dependence on time t. However, due to the high value of the exponent, the vertex may seem to be at rest for some finite time even in this case. Received 29 August 1997; accepted for publication 21 January 1998     

Dynamics of a Non-Autonomous ODE System Occurring in Coagulation Theory

We consider a constant coefficient coagulation equation with Becker–D?ring type interactions and power law input of monomers J ₁(t) = α t ^ω, with α > 0 and . For this infinite dimensional system we prove solutions converge to similarity profiles as t and j converge to infinity in a similarity way, namely with either or constants, where is a function of t only. This work generalizes to the non-autonomous case a recent result of da Costa et al. (2004). Markov Processes Relat. Fields 12, 367–398. and provides a rigorous derivation of formal results obtained by Wattis J. Phys. A: Math. Gen. 37, 7823–7841. The main part of the approach is the analysis of a bidimensional non-autonomous system obtained through an appropriate change of variables; this is achieved by the use of differential inequalities and qualitative theory methods. The results about rate of convergence of solutions of the bidimensional system thus obtained are fed into an integral formula representation for the solutions of the infinite dimensional system which is then estimated by an adaptation of methods used by da Costa et al. (2004). Markov Processes Relat. Fields 12, 367–398.      

A Collocation-Based Algorithm for Computing the Pull-In Range of a Class of Phase-Locked Loops

The pull-in range (ω_p) of a phase-locked loop (PLL) is defined as the maximum value of loop detuning ω_0s for which pull-in occurs from anywhere on the PLL's phase plane. That is, pull-in is guaranteed from anywhere on the phase plane if ω_0s < ω_p. Simple approximation is available for computing ω_p for the high gain PLL where saddle-node bifurcation occurs at ω_0s = ω_p. Unlike the high gain case, a simple approximation for ω_p is not available for the low gain case where bifurcation from a separatrix cycle occurs at ω_0s = ω_p. The vector field model for a class of second-order PLLs is shown to have rotational properties, which imply the existence of a separatrix cycle. The external stability of this separatrix cycle is an indicator of the type of bifurcation (saddle-node or separatrix cycle) which terminates the limit cycle associated with the PLL's stable false lock state and the PLL pulls-in (i.e. achieve phase lock). A formula is given for determining the separatrix cycle's stability, which indicates that these paratrix cycle is externally stable for small values of closed loop gain. A collocation-based algorithm is presented for computing the PLL's separatrix cycle and the value of pull-in range frequency ω_0s = ω_p at which a stable separatrix cycle exists.     

Intermediate Asymptotics Beyond Homogeneity and Self-Similarity: Long Time Behavior for u t = Δϕ(u)

We investigate the long time asymptotics in L¹₊(R) for solutions of general nonlinear diffusion equations u_t = Δϕ(u). We describe, for the first time, the intermediate asymptotics for a very large class of non-homogeneous nonlinearities ϕ for which long time asymptotics cannot be characterized by self-similar solutions. Scaling the solutions by their own second moment (temperature in the kinetic theory language) we obtain a universal asymptotic profile characterized by fixed points of certain maps in probability measures spaces endowed with the Euclidean Wasserstein distance d₂. In the particular case of ϕ(u) ~ u^m at first order when u ~ 0, we also obtain an optimal rate of convergence in L¹ towards the asymptotic profile identified, in this case, as the Barenblatt self-similar solution corresponding to the exponent m. This second result holds for a larger class of nonlinearities compared to results in the existing literature and is achieved by a variation of the entropy dissipation method in which the nonlinear filtration equation is considered as a perturbation of the porous medium equation.     

Oscillatory Couette flow in the presence of an inclined magnetic field

Oscillatory MHD Couette flow of electrically conducting fluid between two parallel plates in a rotating system in the presence of an inclined magnetic field is considered when the upper plate is held at rest and the lower plate oscillates non-torsionally . An exact solution of the governing equations has been obtained by using Laplace transform technique. Asymptotic behavior of the solution is analyzed for M ² ≪1, K ² ≪1 and ω ≪1 and for large M ², K ² and ω. Numerical results of velocities are depicted graphically and the frictional shearing stresses are presented in tables. It is found that a thin boundary layer is formed near the lower plate, for large values of rotation parameter K ², Hartman number M ² and frequency parameter ω. The thickness of this boundary layer increases with increase in inclination of the magnetic field with the axis of rotation.     

An investigation of electromagnetic wave propagation in plasma by shock tube

This paper presents the electromagnetic wave propagation characteristics in plasma and the attenuation coefficients of the microwave in terms of the parameters he, v, w, L, wb. The φ800 mm high temperature shock tube has been used to produce a uniform plasma. In order to get the attenuation of the electromagnetic wave through the plasma behind a shock wave, the microwave transmission has been used to measure the relative change of the wave power. The working frequency is f = (2-35)GHz (ω=2πf, wave length A =15cm-8mm). The electron density in the plasma is ne = (3&#215;10^10-1&#215;10^14) cm^-3. The collision frequency v = (1&#215;10^8-6&#215;10^10) Hz. The thickness of the plasma layer L = (2-80)cm. The electron circular frequency ωb=eBo/me, magnetic flux density B0 = (0-0.84)T. The experimental results show that when the plasma layer is thick (such as L/λ≥10), the correlation between the attenuation coefficients of the electromagnetic waves and the parameters ne,v,ω, L determined from the measurements are in good agreement with the theoretical predictions of electromagnetic wave propagations in the uniform infinite plasma. When the plasma layer is thin (such as when both L and A are of the same order), the theoretical results are only in a qualitative agreement with the experimental observations in the present parameter range, but the formula of the electromagnetic wave propagation theory in an uniform infinite plasma can not be used for quantitative computations of the correlation between the attenuation coefficients and the parameters ne,v,ω, L. In fact, if ω&lt;ωp, v^2&lt;&lt;ω^2, the power attenuations K of the electromagnetic waves obtained from the measurements in the thin-layer plasma are much smaller than those of the theoretical predictions. On the other hand, if ω&gt;ωp, v^2&lt;&lt;ω^2 (just v≈f), the measurements are much larger than the theoretical results. Also, we have measured the electromagnetic wave power attenuation value under the magnetic field and without a magnetic field. The result indicates that the value measured under the magnetic field shows a distinct improvement.