Wall pressure fluctuations of a turbulent separated and reattaching flow affected by an unsteady wake

The spatio-temporal characteristics of the wall-pressure fluctuations in separated and reattaching flows over a backward-facing step were investigated through pressure-velocity joint measurements carried out using multiple-arrayed microphones and split-film probes. A spoke-wheel-type wake generator was installed upstream of the backward-facing step. The flow structure at the effective forcing frequency (St _f=0.2) was found to be well organized in terms of wall pressure spectrum, cross-correlation, wavenumber-frequency spectrum, and wavelet auto-correlation. Introduction of the unsteady wake (St _f=0.2) reduced the reattachment length by 10%. In addition, the unsteady wake enhanced the turbulence intensity near the separation edge and, as a consequence, enhanced the quadrupole sound sources; however, the turbulence intensity near the reattachment region was weakened and the overall flow noise was attenuated. The greater organization of the flow structure induced by the unsteady wake led to a weakening of the dipole sound sources, which are the dominant sound sources in this system. The dipole sound sources generated by wall pressure fluctuations were calculated using Curles integral formula.Abbreviations AR Aspect ratio - SBF Spatial box filtering Roman symbols C _p Wall pressure fluctuation coefficient, p/0.5U ² - H Step height of backward-facing step (mm) - H _s Shape factor (H _s = ^*/) - R _s Distance from acoustic source point to observation point (m) - Re _H Reynolds number, U H/ - St The reduced frequency, fH/U - St _f Normalized forcing frequency by unsteady wake, f _p H/U - T Vortex shedding period (s) - U Free-stream velocity (m/s) - a Speed of sound (m/s) - f Frequency (Hz) - f _p Wake passing frequency (Hz) - k Turbulent kinetic energy (m²/s²) - k _x Streamwise wave number (1/m) - k _z Spanwise wave number (1/m) - l _j Cosine of angle - p Instantaneous wall pressure (Pa) - p _rms Root-mean-square of wall pressure (Pa) - p _SBF Spatial box filtered wall pressure (Pa) - p _d Dipole sound source (Pa) - p _w Conditionally-averaged wall pressure (Pa) - q Dynamic pressure, 0.5U ² (Pa) - r Distance from origin to observation point (mm) - u _c Convection velocity (m/s) - u_max Root-mean-square of streamwise velocity (m/s) - x _R Time-mean reattachment length (mm) Greek symbols _p Forward-flow time fraction - Auto-correlation of pressure at x ₀ - Two-dimensional cross-correlation of pressure with streamwise separation interval , spanwise separation interval , and time delay , at (x ₀, z ₀) - Boundary layer thickness (mm, _99%) - ^* Displacement thickness (mm, ) - _ij Kroneckers delta function - Phase angle (°) - Wavelength (mm) - Momentum thickness (mm, ) - Angle between vertical axis and observation point (°) - Density (kg/m³) - Time delay (s) - Streamwise separation interval (m) - Spanwise separation interval (m) - _p (f; x ₀) Autospectrum of pressure measured at x ₀ (Pa² s) - _pp (, ; x ₀) Streamwise cross spectrum of pressure at x ₀ (Pa² s) - _pp (, , ; x ₀, z ₀) Streamwise and spanwise cross spectrum of pressure at (x ₀, z ₀) (Pa² s) - _pp (kx, ; x ₀) Streamwise wavenumber-frequency spectrum of pressure at x ₀ (Pa² s) - _pp (kx, kz, ; x ₀, z ₀) Two-dimensional wavenumber-frequency spectrum of pressure at (x ₀, z ₀) (Pa² s)     

The value of the critical exponent for reaction-diffusion equations in cones

Let D R ^N be a cone with vertex at the origin i.e., D = (0, )x where S ^N–1 and x D if and only if x = (r, ) with r=¦x¦, . We consider the initial boundary value problem: u _t = u+u ^p in D×(0, T), u=0 on Dx(0, T) with u(x, 0)=u ₀(x) 0. Let ₁ denote the smallest Dirichlet eigenvalue for the Laplace-Beltrami operator on and let ₊ denote the positive root of (+N–2) = ₁. Let p ^* = 1 + 2/(N + ₊). If 1 < p < p ^*, no positive global solution exists. If p>p ^*, positive global solutions do exist. Extensions are given to the same problem for u _t=+¦x¦ u ^p .This research was supported in part by the Air Force Office of Scientific Research under Grant # AFOSR 88-0031 and in part by NSF Grant DMS-8 822 788. The United States Government is authorized to reproduce and distribute reprints for governmental purposes not withstanding any copyright notation therein.     

Nonlocal dispersion in media with continuously evolving scales of heterogeneity

General nonlocal diffusive and dispersive transport theories are derived from molecular hydrodynamics and associated theories of statistical mechanical correlation functions, using the memory function formalism and the projection operator method. Expansion approximations of a spatially and temporally nonlocal convective-dispersive equation are introduced to derive linearized inverse solutions for transport coefficients. The development is focused on deriving relations between the frequency-and wave-vector-dependent dispersion tensor and measurable quantities. The resulting theory is applicable to porous media of fractal character.Nomenclature C _v (t) particle velocity correlation function - C _v ,(t) particle fluctuation velocity correlation function - C _j (x,t) current correlation function - D(x,t) dispersion tensor - D(x,t) fluctuation dispersion tensor - f ₀(x,p) equilibrium phase probability distribution function - f(x, p;t) nonequilibrium phase probability distribution function - G(x,t) conditional probability per unit volume of finding a particle at (x,t) given it was located elsewhere initially - (k,t) Fourier transform ofG(x,t) - G(x,t) fluctuation conditional probability per unit volume of finding a particle at (x,t) given it was located elsewhere initially - k wave vector - K(t) memory function - L Liouville operator - m mass - p(t) particle momentum coordinate - P = (0)( , (0)) projection operator - Q =I-P projection operator - s real Laplace space variable - S(k, ) time-Fourier transform of(k,t) - t time - v(t) particle velocity vector - v(t) particle fluctuation velocity vector - V phase space velocity - time-Fourier variable - ^(itn)(k) frequency moment of(k,t) - x(t) particle displacement coordinate - x(t) particle displacement fluctuation coordinate - friction coefficient - (t) normalized correlation function General Functions () Dirac delta function - () Gamma function Averages ₀ Equilibrium phase-space average - Nonequilibrium phase-space average - (,) L ² inner product with respect tof ₀     

Heat and mass transfer in laminar flow between parallel porous plates

Summary The problem of heat transfer in a two-dimensional porous channel has been discussed by Terrill [6] for small suction at the walls. In [6] the heat transfer problem of a discontinuous change in wall temperature was solved. In the present paper the solution of Terrill for small suction at the walls is revised and the whole problem is extended to the cases of large suction and large injection at the walls. It is found that, for all values of the Reynolds number R, the limiting Nusselt number Nu increases with increasing R.Nomenclature stream function - 2h channel width - x, y distances measured parallel and perpendicular to the channel walls respectively - U velocity of fluid at x=0 - V constant velocity of fluid at the wall - =y/h nondimensional distance perpendicular to the channel walls - f() function defined in equation (1) - coefficient of kinematic viscosity - R=Vh/ suction Reynolds number - density of fluid - C _p specific heat at constant pressure - K thermal conductivity - T temperature - x=x ₀ position where temperature of walls changes - T ₀, T ₁ temperature of walls for x<x ₀, x>x ₀ respectively - = (T – T ₁)/T ₀ – T ₁) nondimensional temperature - =x/h nondimensional distance along channel - R ^* = Uh/v channel Reynolds number - Pr = C _p/K Prandtl number - _n eigenvalues - B _n() eigenfunctions - B _n ⁽⁰⁾ , () eigenfunctions for R=0 - B ₀ ⁽ⁱ⁾ , B ₀ ⁽ⁱⁱ⁾ ... change in eigenfunctions when R0 and small - K _n constants given by equation (13) - h heat transfer coefficient - Nu Nusselt number - _m mean temperature - C _n constants given by equation (18) - perturbation parameter - B _0i () perturbation approximations to B ₀() - Q = B ₀/ ₀ derivative of eigenfunction with respect to eigenvalue - z nondimensional distance perpendicular to the channel walls - F(z) function defined by (54)     

Theory of single well tests in a leaky aquifer. Analytical solution for non steady flow

The evaluation of a pump test or a slug test in a single well that completely penetrates a leaky aquifer does not yield a unique relation between the hydraulic properties of the aquifer, independent of the testing conditions. If the flow is transient, the drawdown is characterized by a single similarity parameter that does not distinguish between the storativity and the leakage factor. If the flow is quasi stationary, the drawdown is characterized by a single similarity parameter that does not distinguish between the transmissivity and the leakage factor. The general non steady solution, which is derived in closed form, is characterized bythree similarity parameters.Nomenclature a e 0.8905 = auxiliary parameter - b thickness of the aquifer - b _c thickness of the semipervious stratum - B() auxiliary function - f(s),g(s) auxiliary functions in the complex plane - F(t),G(t) auxiliary functions of time - h undisturbed level of the phreatic surface - K conductivity of the aquifer - K _c conductivity of the semipervious stratum - m ₀ leakage factor - m dimensionless leakage factor - N(s) auxiliary function in the complex plane - Q _w (t) discharge flux - Q steady discharge flux - Q ₀ constant discharge flux during limited time - Q(t) dimensionless discharge flux - r ₀ radius of the well - r radial coordinate - r dimensionless radial coordinate - s complex variable - s ₀ pole - S storativity of the aquifer - S _n n'th part of an integration contour - t time - t dimensionless time - T transmissivity of the aquifer - ,,,,, dimensionless parameters - Euler's number - dummy variable - ₁(), ₂() auxiliary functions - (r, t) drawdown - ₀(t) drawdown in the well - (r, t) dimensionless drawdown - ₀(t) dimensionless drawdown in the well     

The secondary flow induced around a sphere in an oscillating stream of elastico-viscous liquid

The present paper is devoted to the theoretical study of the secondary flow induced around a sphere in an oscillating stream of an elastico-viscous liquid. The boundary layer equations are derived following Wang's method and solved by the method of successive approximations. The effect of elasticity of the liquid is to produce a reverse flow in the region close to the surface of the sphere and to shift the entire flow pattern towards the main flow. The resistance on the surface of the sphere and the steady secondary inflow increase with the elasticity of the liquid.Nomenclature a radius of the sphere - b ^ik contravariant components of a tensor - e contravariant components of the rate of strain tensor - F() see (47) - G total nondimensional resistance on the surface of the sphere - g _ik covariant components of the metric tensor - f, g, h secondary flow components introduced in (34) - k ₀ measure of relaxation time minus retardation time (elastico-viscous parameter) - K =k ₀ ²/V ₀ ² , nondimensional parameter characterizing the elasticity of the liquid - n measure of the ratio of the boundary layer thickness and the oscillation amplitude - N, T defined in (44) - p arbitrary isotropic pressure - p _ik covariant components of the stress tensor - p ^ik contravariant components of the stress tensor associated with the change of shape of the material - R =V ₀ a/v, the Reynolds number - S =a/V ₀, the Strouhall number - r, , spherical polar coordinates - u, v, w r, , component of velocity - t time - V(, t) potential velocity distribution around the sphere - V ₀ characteristic velocity - u, v, t, y, P nondimensional quantities defined in (15) - reciprocal of s - density - defined in (32) - defined in (42) - ₀ limiting viscosity for very small changes in deformation velocity - complex conjugate of - oscillation frequency - = ₀/, the kinematic coefficient of viscosity - , defined in (52) - (, y) stream function defined in (45) - =(NT/2n)^1/2 y - /t convective time derivative (1) _ik      

On nonlinear problems of mixed type: A qualitative theory using infinite-dimensional center manifolds

We consider the equation a(y)u_xx+div_y(b(y)_yu)+c(y)u=g(y, u) in the cylinder (–l,l)×, being elliptic where b(y)>0 and hyperbolic where b(y)<0. We construct self-adjoint realizations in L₂() of the operatorAu= (1/a) div_y(b_yu)+(c/a) in the case ofb changing sign. This leads to the abstract problem u_xx+Au=g(u), whereA has a spectrum extending to + as well as to –. For l= it is shown that all sufficiently small solutions lie on an infinite-dimensional center manifold and behave like those of a hyperbolic problem. Anx-independent cross-sectional integral E=E(u, u_x) is derived showing that all solutions on the center manifold remain bounded forx ±. For finitel, all small solutionsu are close to a solution on the center manifold such that u(x)-(x) Ce ^-(1-|x|) for allx, whereC and are independent ofu. Hence, the solutions are dominated by hyperbolic properties, except close to the terminal ends {±1}×, where boundary layers of elliptic type appear.     

Nonlinear rheological behavior of a concentrated spherical silica suspension

Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G^* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< ^–1; = linear viscoelastic relaxation time) and shear-thickening for larger (>^–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and _– (t, ). In the shear-thinning regime, the and dependence of ₊(t,) and _–(t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson     

Theoretical investigation of the effect of an internal heat source on turbulent liquid metal heat transfer in a channel with uniform heat flux

Summary The effect of an internal heat source on the heat transfer characteristics for turbulent liquid metal flow between parallel plates is studied analytically. The analysis is carried out for the conditions of uniform internal heat generation, uniform wall heat flux, and fully established temperature and velocity profiles. Consideration is given both to the uniform or slug flow approximation and the power law approximation for the turbulent velocity profile. Allowance is made for turbulent eddying within the liquid metal through the use of an idealized eddy diffusivity function. It is found that the Nusselt number is unaffected by the heat source strength when the velocity profile is assumed to be uniform over the channel cross section. In the case of a 1/7-power velocity expression, the Nusselt numbers are lower than those in the absence of internal heat generation, and decrease with diminishing eddy conduction. Nusselt numbers, in the absence of an internal heat source, are compared with existing calculations, and indications are that the present results are adequate for preliminary design purposes.Nomenclature A hydrodynamic parameter - a half height of channel - a ₁ a constant, 1+0.01 Pr Re ^0.9 - a ₂ a constant, 0.01 Pr Re ^0.9 - C _p specific heat at constant pressure - D _h hydraulic diameter of channel, 4a - h heat transfer coefficient, q _w/(t _w–t _b) - I ₁ integral defined by (17) - I ₂ integral defined by (18) - k diffusivity parameter, (1+0.01 Pr Re ^0.9)^1/2 - m exponent in power velocity expression - Nu Nusselt number, hD _h/ - Nu ₀ Nusselt number in absence of internal heat generation - Pr Prandtl number, / - Q heat generation rate per volume - q _w wall heat flux - Re Reynolds number for channel, 2/ - s ratio of heat generation rate to wall heat flux, Qa/q _w - T dimensionless temperature, (t _w–t)/(t _w–t _b) - t fluid temperature, t _w wall temperature, t _b fluid bulk temperature - u fluid velocity in x direction, , fluid mean velocity - x longitudinal coordinate measured from channel entrance - x ⁺ dimensionless longitudinal coordinate, 2(x/a)/Pr Re - y transverse coordinate measured from channel centerline - z transverse coordinate measured from channel wall, a–y - molecular diffusivity of heat, /C _p - dummy variable of integration - dummy variable of integration - _H eddy diffusivity of heat - _M eddy diffusivity of momentum - dummy variable of integration - fluid thermal conductivity - _T dimensionless diffusivity, Pr ( _H/) - fluid kinematic viscosity - dummy variable of integration - fluid density - dummy variable of integration - ratio of eddy diffusivity for heat transfer to that for momentum transfer, _H/ _M - average value of - dimensionless velocity distribution, u/     

Translatory accelerating motion of a circular disk in a viscous fluid

The exact solution of the equation of motion of a circular disk accelerated along its axis of symmetry due to an arbitrarily applied force in an otherwise still, incompressible, viscous fluid of infinite extent is obtained. The fluid resistance considered in this paper is the Stokes-flow drag which consists of the added mass effect, steady state drag, and the effect of the history of the motion. The solutions for the velocity and displacement of the circular disk are presented in explicit forms for the cases of constant and impulsive forcing functions. The importance of the effect of the history of the motion is discussed.Nomenclature a radius of the circular disk - b one half of the thickness of the circular disk - C dimensionless form of C ₁ - C ₁ magnitude of the constant force - D fluid drag force - f(t) externally applied force - F() dimensionaless form of applied force - F ₀ initial value of F - g gravitational acceleration - H() Heaviside step function - k magnitude of impulsive force - K dimensionless form of k - M a dimensionless parameter equals to (1+37#x03C0;_s/4_f) - S displacement of disk - t time - t ₁ time of application of impulsive force - u velocity of the disk - V dimensionless velocity - V ₀ initial velocity of V - V _t terminal velocity - parameter in (13) - parameter in (13) - (t) Dirac delta function - ratio of b/a - () function given in (5) - dynamical viscosity of the fluid - kinematic viscosity of the fluid - _f fluid density - _s mass density of the circular disk - dimensionless time - _i dimensionless form of t _i - dummy variable - dummy variable     

Two-phase brine mixtures in the geothermal context and the polymer flood model

Two-phase mixtures of hot brine and steam are important in geothermal reservoirs under exploitation. In a simple model, the flows are described by a parabolic equation for the pressure with a derivative coupling to a pair of wave equations for saturation and salt concentration. We show that the wave speed matrix for the hyperbolic part of the coupled system is formally identical to the corresponding matrix in the polymer flood model for oil recovery. For the class ofstrongly diffusive hot brine models, the identification is more than formal, so that the wave phenomena predicted for the polymer flood model will also be observed in geothermal reservoirs.Roman Symbols A,B coefficient matrices (5) - c(x,t) salt concentration (primary dependent variable) - C(p, s, c, q _t) wave speed matrix (6) - f source term (5) - g acceleration due to gravity (constant) - h _b(p, c) brine specific enthalpy - h _v(p) vapour specific enthalpy - j conservation flux (1) - k absolute permeability (constant) - k _b(s), k_v(s) relative permeabilities of the brine and vapour phases - K conductivity - p(x,t) pressure (primary dependent variable) - q volume flux (Darcy velocity) (3) - s(x,t) brine saturation (primary dependent variable) - t time (primary independent variable) - T=T _sat(p) saturation temperature - u _b(p, c) brine specific internal energy - u _m T rock matrix specific internal energy - u _v(p) vapour specific internal energy - U(x, t) shock velocity - x space (primary independent variable) Greek Symbols porosity (constant) - _b(p, c) brine dynamic viscosity - _v(p) vapour dynamic viscosity - (p, s, c) conservation density (1) - _b(p, c) brine density - _v(p) vapour density Suffixes b brine - m rock matrix - t total - v vapour - S salt - M mass - E energy     

Limiting analytical solutions for complex saturated and unsaturated transport problems

Non-linear diffusion and velocity-dependent dispersion problems are under consideration. The necessary and sufficient conditions allowing the comparison of solutions to the two dimensional convection-dispersion equations with different coefficients are obtained. These conditions provide a framework within which solutions to the complex non-linear problems mentioned above can be estimated by solutions to the problems possessing analytical solvability.Nomenclature c(x, y, t) concentration of solute in solution,ML ^–3 - C(h)=d/dh moisture capacity function - D,D _ij hydrodynamic dispersion coefficient, a second order tensor,L ² T ^–1 - D _L longitudinal hydrodynamic dispersion coefficient,L ² T ^–1 - D _m molecular diffusion coefficient,L ² T ^–1 - D _T transverse hydrodynamic coefficient,L ² T ^–1 - G flow domain for the unsaturated flow problem - G _z , G _w flow domain and complex potential domain, respectively, for the hydrodynamic dispersion problem - h piezometric head,L - I _n given mass flux normal to the boundary,MLT ^–1 - k hydraulic conductivity,LT ^–1 - K(h) unsaturated hydraulic conductivity,LT ^–1 - L continuously differentiable function with respect to all arguments - m porosity - n(x,t) outer normal vector to the boundary - t time,T - V(x, y, t) seepage velocity vector withV=V,LT ^–1 - x Cartesian coordinate system - x horizontal coordinate,L - y vertical coordinate (elevation),L - (x),(x,t) given functions in initial and boundary conditions (3), (4) - ₁(,) angle between vectors ₁c andV - boundary of the flow domain - _L , _T longitudinal and transverse dispersivities, respectively,L - water mass density,ML ^–3 - v _i components of a unit vector in the direction of the outward normal to the boundary - =–kh velocity potential - =/m - stream function defined such thatw=+i is the complex potential - =/m     

Viscous properties of calcium carbonate filled polymer melts

Summary The viscous properties of calcium carbonate filled polyethylene and polystyrene melts were examined. The relative vircosity _r defined in the previous paper gave an asymtptotic value( _r)_l in the range of the shear stress below 10⁵ dyne/cm².( _r)_l of the calcium carbonate filled system was higher than that of the glass beads or glass balloons filled system at the same volume fraction of the filler. Maron-Pierce equation with ₀ = 0.44 was able to approximate the( _r)_l — relationship. However, it was deduced here that the high value of( _r)_l of calcium carbonyl filled system was due to the apparent increase of and this increase was attributed to the fixed polymer layer formed on the powder particle. By assuming the particle as a sphere with a diameter of 2 µm, the thickness of the fixed polymer layer was estimated as about 0.17 µm. The yield stress estimated from the Casson's plots increased exponentially with.

---

Zusammenfassung Es wurden die viskosen Eigenschaften von Polyäthylen-und Polystyrol-Schmelzen untersucht, die mit Kalziumkarbonat-Teilchen gefüllt waren. Für die relative Viskosität _r, wie sie in einer vorangegangenen Veröffentlichung definiert worden war, ergab sich bei Schubspannungen unterhalb 10⁵ dyn/cm² ein asymptotischer Wert( _r)_l. Dieser war bei den mit Kalziumkarbonat gefüllten Schmelzen höher als bei Schmelzen, die bis zur gleichen Volumenkonzentration mit Glaskugeln oder Glasballons gefüllt waren. Die ( _r) _l —-Abhängigkeit ließ sich durch eine Gleichung nachMaron und Pierce mit ₀ = 0,44 beschreiben. Es wurde jedoch geschlossen, daß der hohe( _r)_l-Wert der mit Kalziumkarbonat gefüllten Schmelzen auf eine scheinbare Zunahme von zurückzuführen ist, verursacht durch eine feste Polymerschicht auf der Teilchenoberfläche. Unter Annahme kugelförmiger Teilchen mit einem Durchmesser von 2 µm ließ sich die zugeordnete Schichtdicke zu 0,17 µm abschätzen. Die mittels der Casson-Beziehung geschätzte Fließspannung ergab eine exponentielle-Abhängigkeit.

---

---

With 7 figures and 1 table     

On some nearly viscometric flows of viscoelastic liquids

Superposition of oscillatory shear imposed from the boundary and through pressure gradient oscillations and simple shear is investigated. The integral fluid with fading memory shows flow enhancement effects due to the nonlinear structure. Closed-form expressions for the change in the mass transport rate are given at the lowest significant order in the perturbation algorithm. The elasticity of the liquid plays as important a role in determining the enhancement as does the shear dependent viscosity. Coupling of shear thinning and elasticity may produce sharp increases in the flow rate. The interaction of oscillatory shear components may generate a steady flow, either longitudinal or orthogonal, resulting in increases in flow rates akin to resonance, and due to frequency cancellation, even in the absence of a mean gradient. An algorithm to determine the constitutive functions of the integral fluid of order three is outlined.Nomenclature A _n Rivlin-Ericksen tensor of order . - A _k Non-oscillatory component of the first order linear viscoelastic oscillatory velocity field induced by the kth wave in the pressure gradient - d Half the gap between the plates - e _x, e _z Unit vectors in the longitudinal and orthogonal directions, respectively - G(s) Relaxation modulus - G History of the deformation - Stress response functional - I() Enhancement defined as the ratio of the frequency dependent part of the discharge to the frequencyindependent part of it at the third order - I ^*() Enhancement defined as the ratio of the increase in discharge due to oscillations to the total discharge without the oscillations - k Power index in the relaxation modulus G(s) - k _i ^–1 Relaxation times in the Maxwell representation of the quadratic shear relaxation modulus (s ₁, s ₂) - m _i ^–1, n _i ^–1 Relaxation times in the Maxwell representations of the constitutive functions ₁(s ₁,s ₂,s ₃) and ₄ (s ₁, s ₂,s ₃), respectively - P Constant longitudinal pressure gradient - p Pressure field - _mx ,⁽³⁾ _nz ,⁽³⁾ Mean volume transport rates at the third order in the longitudinal and orthogonal directions, respectively - ₀,⁽³⁾, ₁,⁽³⁾ Frequency independent and dependent volume transport rates, respectively, at the third order - s = t- Difference between present and past times t and      

On properties of solutions of functional differential equations continuously differentiable on (0, +∞)

We establish the saddle-point property of the system of functional differential equations (t) = Ax(t) + Bx((t)) + C ((t)) + f (x(t), x((t))), (0) = 0.Translated from Neliniini Kolyvannya, Vol. 7, No. 3, pp. 302–310, July–September, 2004.     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

Note on the cancellation of contaminating noise in the measurement of turbulent wall pressure fluctuations

A method is presented for the cancellation of wide band contaminating noise occurring within internal flow configurations such as rectangular channels and pipes. Facility generated noise within these flow systems contaminates the turbulent wall pressure signature at low frequencies thus preventing the possible extraction of useful information. The proposed methodology utilizes the signals from two flush mounted wall pressure transducers. A first estimate for the one-point spectral density is obtained using a least mean square algorithm. A secondary correction to this estimate is obtained by taking advantage of the planar homogeneity of the turbulence. The application of the technique is demonstrated in a fully developed turbulent channel flow for which a more than 40 dB cancellation is obtained at low frequencies. In this low frequency range, the power spectral density is shown to have an approximate quadratic dependence, substantiating past theoretical predictions reported in the literature.List of symbols d transducer diameter - d ⁺ non-dimensional transducer diameter, - f frequency (Hz) - h channel half-height - l spanwise separation of transducers - p _N (t) signal due to contaminating noise - p _q (t) system output signal - p _R (t) turbulent wall pressure signal at reference transducer - p ₁(t) turbulent wall pressure signal at primary transducer - R _h channel Reynolds number, - s _R (t) pressure signal from reference transducer - s ₁ (t) pressure signal from primary transducer - t time - U channel centerline velocity, (maximum velocity) - u ^* shear velocity, - W(f) optimized filter function used for cancellation (Fourier transform) - correction factor, (noise to signal ratio) - ² coherence function - _z lateral turbulence macro-scale - _N correlation length scale of contaminating noise - v kinematic viscosity - fluid density - _wall wall shear stress - _pN (f) auto spectral density of p _N (t) - _pR (f) auto spectral density of p _R (t) - ₁ p auto spectral density of p ₁(t) - _sR (f) auto spectral density of s _R (t) - ₁ s(f) auto spectral density of s ₁(t) - _1ssr cross spectral density of s ₁ (t) with s _R (t)     

Transient flow towards a well in an aquifer including the effect of fluid inertia

Transient propagation of weak pressure perturbations in a homogeneous, isotropic, fluid saturated aquifer has been studied. A damped wave equation for the pressure in the aquifer is derived using the macroscopic, volume averaged, mass conservation and momentum equations. The equation is applied to the case of a well in a closed aquifer and analytical solutions are obtained to two different flow cases. It is shown that the radius of influence propagates with a finite velocity. The results show that the effect of fluid inertia could be of importance where transient flow in porous media is studied.List of symbols b Thickness of the aquifer, m - c ₀ Wave velocity, m/s - k Permeability of the porous medium, m² - n Porosity of the porous medium - p( ,t) Pressure, N/m² - Q Volume flux, m³/s - r Radial coordinate, m - r _w Radius of the well, m - s Transform variable - S Storativity of the aquifer - S _d(r, t) Drawdown, m - t Time, s - T Transmissivity of the aquifer, m²/s - ( ,t) Velocity of the fluid, m/s - Coordinate vector, m - z Vertical coordinate, m - Coefficient of compressibility, m²/N - Coefficient of fluid compressibility, m²/N - Relaxation time, s - (r, t) Hydraulic potential, m - Dynamic viscosity of the fluid, Ns/m² - Dimensionless radius - Density of the fluid, Ns²/m⁴ - (, ) Dimensionless drawdown - Dimensionless time - , x Dummy variables - ₀, ₁ Auxilary functions     

Uniform Convergence for Approximate Traveling Waves in Linear Reaction–Diffusion–Hyperbolic Systems

In this paper we study linear reaction–hyperbolic systems of the form , (i = 1, 2, ..., n) for x > 0, t > 0 coupled to a diffusion equation for p ₀ = p ₀(x, y, θ, t) with “near-equilibrium” initial and boundary data. This problem arises in a model of transport of neurofilaments in axons. The matrix (k _ij ) is assumed to have a unique null vector with positive components summed to 1 and the v _j are arbitrary velocities such that . We prove that as the solution converges to a traveling wave with velocity v and a spreading front, and that the convergence rate in the uniform norm is , for any small positive α.     

Thermal stresses in rectangular plates

Summary A method of determining the thermal stresses in a flat rectangular isotropic plate of constant thickness with arbitrary temperature distribution in the plane of the plate and with no variation in temperature through the thickness is presented. The thermal stress have been obtained in terms of Fourier series and integrals that satisfy the differential equation and the boundary conditions. Several examples have been presented to show the application of the method.Nomenclature x, y rectangular coordinates - _x, _y direct stresses - _xy shear stress - ø Airy's stress function - E Young's modulus of elasticity - coefficient of thermal expansion - T temperature - ² Laplace operator: - ⁴ biharmonic operator - 2a length of the plate - 2b width of the plate - a/b aspect ratio - a _mr, b_ms, c_nr, d_ns Fourier coefficients defined in equation (6) - _m=m/a m=1, 2, 3, ... _n=n/2a n=1, 3, 5, ... - _r=r/b r=1, 2, 3, ... _s=s/2b s=1, 3, 5, ... - A _m, B_m, C_n, D_n, E_r, F_r, G_s, H_s Fourier coefficients - K _rand L _s Fourier coefficients defined in equation (20) - direct stress at infinity - T ₁(x, y) temperature distribution symmetrical in x and y - T ₂(x, y) temperature distribution symmetrical in x and antisymmetrical in y - T ₃(x, y) temperature distribution antisymmetrical in x and symmetrical in y - T ₄(x, y) temperature distribution antisymmetrical in x and y