A technique for evaluation of three-dimensional behavior in turbulent boundary layers using computer augmented hydrogen bubble-wire flow visualization

A system is described which allows the recreation of the three-dimensional motion and deformation of a single hydrogen bubble time-line in time and space. By digitally interfacing dualview video sequences of a bubble time-line with a computer-aided display system, the Lagrangian motion of the bubble-line can be displayed in any viewing perspective desired. The u and v velocity history of the bubble-line can be rapidly established and displayed for any spanwise location on the recreated pattern. The application of the system to the study of turbulent boundary layer structure in the near-wall region is demonstrated.List of Symbols Reynolds number based on momentum thickness u /v - t⁺ nondimensional time - u shear velocity - u local streamwise velocity, x-direction - u ⁺ nondimensional streamwise velocity - v local normal velocity, -direction - x ⁺ nondimensional coordinate in streamwise direction - ⁺ nondimensional coordinate normal to wall - ₊ ^wire wire nondimensional location of hydrogen bubble-wire normal to wall - z ⁺ nondimensional spanwise coordinate - momentum thickness - v kinematic viscosity - _W wall shear stress     

Non-persistence of Homoclinic Connections for Perturbed Integrable Reversible Systems

The dynamics of an analytic reversible vector field (X,) is studied in with one real parameter close to 0; X=0 is a fixed point. The differential D_x (0,0) generates an oscillatory dynamics with a frequency of order 1—due to two simple, opposite eigenvalues lying on the imaginary axis—and it also generates a slow dynamics which changes from a hyperbolic type—eigenvalues are —to an elliptic type—eigenvalues are —as passes trough 0. The existence of reversible homoclinic connections to periodic orbits is known for such vector fields. In this paper we study a particular subclass of such vector fields, obtained by small reversible perturbations of the normal form. We give an explicit condition on the perturbation, generically satisfied, which prevents the existence of a homoclinic connections to 0 for the perturbed system. The normal form system of any order admits a reversible homoclinic connection to 0, which then does not survive under perturbation of higher order. It will be seen that normal form essentially decouples the hyperbolic and elliptic part of the linearization to any chosen algebraic order. However, this decoupling does not persist arbitrary reversible perturbation, which finally causes the appearance of small amplitude oscillations.     

Flow of second-order fluids over an enclosed rotating disc

Summary This paper is devoted to a study of the flow of a second-order fluid (flowing with a small mass rate of symmetrical radial outflow m, taken negative for a net radial inflow) over a finite rotating disc enclosed within a coaxial cylinderical casing. The effects of the second-order terms are observed to depend upon two dimensionless parameters ₁ and ₂. Maximum values ₁ and ₂ of the dimensionless radial distances at which there is no recirculation, for the cases of net radial outflow (m>0) and net radial inflow (m<0) respectively, decrease with an increase in the second-order effects [represented by T(=₁+₂)]. The velocities at ₁ and ₂ as well as at some other fixed radii have been calculated for different T and the associated phenomena of no-recirculation/recirculation discussed. The change in flow phenomena due to a reversal of the direction of net radial flow has also been studied. The moment on the rotating disc increases with T.Nomenclature , , z coordinates in a cylindrical polar system - z ₀ distance between rotor and stator (gap length) - =/z ₀, dimensionless radial distance - =z/z ₀, dimensionless axial distance - _s = _s/z₀, dimensionless disc radius - V =(u, v, w), velocity vector - dimensionless velocity components - uniform angular velocity of the rotor - , p fluid density and pressure - P =p/(₂ z ₀₂ ² , dimensionless pressure - ₁, ₂, ₃ kinematic coefficients of Newtonian viscosity, elastico-viscosity and cross-viscosity respectively - ₁, _{2 2}/z ₀ ² , resp. ₃/z ₀ ² , dimensionless parameters representing the ratio of second-order and inertial effects - m = , mass rate of symmetrical radial outflow - l a number associated with induced circulatory flow - Rm =m/(z ₀₁), Reynolds number of radial outflow - R _l =l/(z ₀₁), Reynolds number of induced circulatory flow - Rz =z ₀ ² /₁, Reynolds number based on the gap - ₁, ₂ maximum radii at which there is no recirculation for the cases Rm>0 and Rm<0 respectively - ₁(T), ₂(T) ₁ and ₂ for different T - U _1(T) ⁽⁺⁾ = dimensionless radial velocity, Rm>0 - V _1(T) ⁽⁺⁾ = , dimensionless transverse velocity, Rm>0 - U _2(T) ^(–) = , dimensionless radial velocity, Rm=–Rn<0, m=–n - V _2(T) ^(–) = , dimensionless transverse velocity, Rm<0 - C _m moment coefficient     

Der Spannungs- und Verschiebungszustand drehsymmetrischer Membrane mit beliebig gekrümmter Meridiankurve

Zusammenfassung Die Einführung von Zylinderkoordinaten (x, r, ) in die Gleichgewichtsbedingungen der Schnittkräfte bzw. in die Beziehungen zwischen Verzerrung und Verschiebungen am differentialen Schalenabschnitt ermöglicht die Berechnung des Spannungs- und Verschiebungszustandes von drehsymmetrischen Membranen mit beliebig gekrümmter Meridiankurve auf die Integration einer einfachen, linearen partiellen Differentialgleichung zweiter Ordnung für eine charakteristische FunktionF bzw. zurückzuführen. Eine geschlossene Lösung und damit eine Darstellung der Schnittkräfte und Verschiebungen durch explizite Formeln ist bei harmonischer Belastung cosn für zwei Funktionsgruppen=x ² und=x ^–3 möglich. Im Sonderfall der drehsymmetrischen und der antimetrischen Belastung mitn=0 undn=1 gelten die Gleichungen der Schnitt- und Verschiebungsgrößen für eine beliebige Meridianfunktion=(). Die Betrachtungen der Randbedingungen offener Schalen bei harmonischer Belastung geben über die infinitesimalen Deformationen einer drehsymmetrischen Membran mit überall negativer Krümmung Aufschluß.     

Multiaxial elongation of polyisobutylene with various and changing strain rate ratios

The multiaxial elongational rheometer equipped with rotary clamps is modified such that in addition to simple, equibiaxial and multiaxial elongations also tests with new modes of elongation can be performed. As an example, polyisobutylene is elongated with a ratio of the principal strain rates of and magnitudes of the maximum strain rate , 0.04 and 0.08 s^–1. As a test result, the first elongational viscosityµ ₁ (t) is obtained which follows closely the linear viscoelastic shear viscosity . In contrast, the second elongational viscosityµ ₂ (t) remains below . By means of a further modification of the rheometer, the test modes can be varied during the deformation period. This allows one to investigate the influence of a well-defined rheological pre-history on the following rheological behaviour. As an example a variation ofm = 0.5 2 was performed. The measured normal-stress differences superpose from the single steps of deformation similar to the linear viscoelastic prediction.Dedicated to Prof. F. R. Schwarzl on the occasion of his 60th birthday     

Homologies of Intersections of Pseudomanifolds with Local Homotopic Conditions on Links of Strata

In terms of local homotopic properties of the links of strata of an n-dimensional PL-pseudomanifold X, we obtain a sufficient condition for the natural homomorphisms of the jth intersection homology groups with perversity multiindices and to be isomorphisms for all j i, where i < n – 3.     

Strömungstechnische Gesichtspunkte der Flußbelastung durch Abwasser aus Kläranlagen

Zusammenfassung Für ein im Durchlauf betriebenes System bestehend aus einem Fluß (Vorfluter) und den angeschlossenen Kläranlagen wird eine Methode zur Bestimmung der Vorfluterbelastung durch die eingeleiteten Klärwässer angegeben. Die Methode erfaßt mit Rücksicht auf die Anwendung des Verursacherprinzips im Gewässerschutz die Belastung durch jede Kläranlage für sich, und zwar in Abhängigkeit von der Wasserführung, den Emissionsraten der betreffenden Kläranlage und dem Selbstreinigungsvermögen von den organischen Stoffen aus der betreffenden Kläranlage. Die abhängigen Veränderlichen sind mit der Fließgeschwindigkeit gewichtete Mittelwerte von Schmutzstoffdichten über den Vorfluterquerschnitt. Im Falle konstanter Vorflutertemperatur und zeitunabhängiger Struktur der Klärwässer ergeben sich beispielsweise für die abhängigen Veränderlichen einfache analytische Darstellungen, welche sich als spezielle Formen des -Theorems erweisen. Es wird gezeigt, bei einem unendlich langen Vorfluter mit konstantem Volumenstrom stromabwärts der Klärwassereinleitungen stimmen die erwähnten gewichteten Mittelwerte mit den entsprechenden ungewichteten stromabwärts der Klärwassereinleitungen überein. Die entwickelte Methode kann leicht erweitert werden, um den Sauerstoffschwund im Vorfluter durch jede Kläranlage für sich zu bestimmen.

---

Fluid mechanical aspects of river pollution by effluents from waste treatment plants

The pollution of a river by effluent inflows from waste treatment plants is modeled under steady-state conditions. With respect to modern policies of environmental protection the method describes the river pollution by each plant separately, depending on the flow conditions, the emission rates of the plant and the microbiological decomposition of the biodegradable matter from the plant. Each dependent variable is a weighted cross-sectional mean of a density of organic matter. If the water temperature is constant and the composition of each effluent is independent of time the method gives simple analytic expressions for the dependent variables, which prove to be special versions of the -theorem. It is shown for an infinitely long river of constant volume rate of flow downstream of the effluent inflows: the weighted means mentioned agree with the corresponding nonweighted downstream of the effluent inflows. The present paper can easily the extended to determine the oxygen deficit in the river due to each plant.

---

Bezeichnungen a Anzahl der Kläranlagen - D(t_b) Kennzahl, Einführung in 4.3 - ^eA Emissionsrate der abbaubaren or ganischen Verschmutzung aus der -ten Kläranlage - ^eU Emissionsrate der nichtabbaubaren organischen Verschmutzung aus der -ten Kläranlage - Vorfluterquerschnitt, Einführung in Gl. (4) - F Flächeninhalt von - dF Betrag eines Flächenelements, Einführung in Gl. (6) - J_A Diffusionsstromdichten, Einführung in Gl. (2) bzw. Gl. (3) - L Anzahl der Stromstrecken - M Gesamtmasse der abbaubaren or- ganischen Verschmutzung in den N Teilchen, Einführung in Gl. (17) - N Anzahl der verschmutzten Flußwasserteilchen, welche die -te Nahfeldvermischungszone während des Zeitintervalles t_a t_b für immer verlassen - P(x, t, x, t_c) Teilchendichte, Einführung in Gl. (11) und Gl.(12) - Q Selbstreinigungsvermögen, Einführung in Gl.(26) - t Zeitpunkt, Einführung in Gl.(11) - t, t_b Intervallgrenzen, Einführung in 4.1 - t_c Zeitpunkt, Einführung in Gl.(11) - t Zeitdifferenz, Einführung im Anschluß an Gl.(10) - t^* charakteristische Zeit, Einführung in 4.3 - Strömungsgeschwindigkeit Komponente von ¯b in Richtung der zu Tal weisenden Oberflächennormalen eines Vorfluterquerschnitts, Einführung in Gl. (5) und Gl. (6) - Volumenstrom, Einführung in Gl. (7) - x Ortsvektor - x Ortsvektor eines bestimmten markierten Teilchens zur Zeit t_c, Einführung in Gl.(11) - x längs der Stromachse gemessene Längenkoordinate - x x-Koordinate des Vorfluterquerschnitts durch x - x,x₊₁ x-Koordinaten der Vorfluterquerschnitte, welche die -te Stromstrecke stromaufwärts bzw. stromabwärts begrenzen. Einführung in 4.2. - transformierte Variable, Einführung in Gl.(65) - Zeitvariable - (t_b) Kennzahl, Einführung in 4.3. - Masse der abbaubaren organischen Verschmutzung in dem markierten Teilchen, Einführung in Gl.(14) - , Integrationsvariablen, Einführung in Gl.(38) bzw. Gl.(28) - _A durch die -te Kläranlage bedingte Dichte der abbaubaren organischen Verschmutzung - _U durch die -te Kläranlage bedingte Dichte der nichtabbaubaren organischen Verschmutzung - Mittelwerte von bzw· , Einführung in Gl.(31) bzw. Gl.(8) - _m -Wert zu einem Maximum, Einführung in Gl.(31) - Verhältnis zweier Mittelwerte, Einführung in Gl.(64) - stochastischer Mittelwert einer Zufallsgröße Y - Y Schwankung einer Zufallsgröße Y um den stochastisehen Mittelwert - Mittlung über den Vorfluterquerschnitt Der saubere Vorfluter sei definiert durch Standardwerte für Mindestanforderungen an die Flußwasserqualität. Vorschläge für solche Standardwerte werden in jüngster Zeit unter Berücksichtigung des Umweltschutzes ausführlich diskutiert ([1]; [2], S.- K 13 -).     

The Extinction Behavior of the Solutions for a Class of Reaction-Diffusion Equations

1　ＩｎｔｒｏｄｕｃｔｉｏｎａｎｄＬｅｍｍａｓＴｈｅｒｅａｒｅｍａｎｙｒｅｓｕｌｔｓａｂｏｕｔｅｘｉｓｔｅｎｃｅ (ｇｌｏｂａｌｏｒｌｏｃａｌ)ａｎｄａｓｙｍｐｔｏｔｉｃｂｅｈａｖｉｏｒｏｆｓｏｌｕｔｉｏｎｓｆｏｒｒｅａｃｔｉｏｎ＿ｄｉｆｆｕｓｉｏｎｅｑｕａｔｉｏｎｓ[1- 9].Ｂｙｔｈｅａｉｄｓｏｆｒｅｓｕｌｔｓ[2 ,3]ｏｆｅｑｕａｔｉｏｎ ｕ/ ｔ=Δｕ-λ｜ｕ｜γ- 1ｕｗｉｔｈｉｎｉｔｉａｌ＿ｂｏｕｎｄａｒｙｖａｌｕｅｓ,ｐａｐｅｒ [4 ]ｓｔｕｄｉｅｄｔｈｅｐｒｏｂｌｅｍｏｆ ｕ/ ｔ=Δｕ-λ｜ｅβｔｕ｜γ- …     

Singular solutions of some quasilinear elliptic equations

We study isolated singularities of the quasilinear equation in an open set of ^N , where 1 < p N, p -1 q < N(p — 1)/ (N -p). We prove that, for any positive solution, if a singularity at the origin is not removable then either or u(x)/(x) any positive constant as x 0 where is the fundamental solution of the p-harmonic equation: . Global positive solutions are also classified.     

Invertible Contractions and Asymptotically Stable ODE’S that are not \mathcal{C}^{1}-Linearizable

We present an example of a contraction diffeomorphism in infinite dimensions that is not -linearizable, and we construct a regular ordinary differential equation in a Hilbert space whose time-one map is that diffeomorphism. With this we have an example of an asymptotically stable ODE that is not -conjugate to its linear part.     

Traveling Waves of a Singularly Perturbed System of Integral–Differential Equations Arising from Neuronal Networks

A nonlinear nonlocal model arising from synaptically coupled neuronal networks with two integral terms is considered. The existence and stability of several traveling wave solutions are established by using ideas in differential equations and functional analysis. Steady-state solutions of some inhomogeneous integral–differential equations are also investigated. We consider several types of kernel functions: (I) positive functions, such as and , where ρ>0 is a constant; (II) nonnegative kernels with compact supports, for examples, (i) 1$$" align="middle" border="0"> , and (ii) {\pi\over 2}$$" align="middle" border="0"> ; (III) Mexican hat type kernel functions, such as and , where A>B>0 and a>b>0 are constants.Dedicated to Professor Yulin Zhou and Professor Boling Guo on the Occassions of their birthdays.     

Bifurcation Phenomena from Standing Pulse Solutions in Some Reaction-Diffusion Systems

Bifurcation phenomena from standing pulse solutions of the problem is considered. (>0) is a sufficiently small parameter and is a positive one. It is shown that there exist two types of destabilization of standing pulse solutions when decreases. One is the appearance of travelling pulse solutions via the static bifurcation and the other is that of in-phase breathers via the Hopf bifurcation. Furthermore which type of destabilization occurs first with decreasing is discussed for the piecewise linear nonlinearities f and g.     

Dynamic properties of some polymer solutions subjected to a steady shear superimposed on an oscillatory shear flow

Summary Results of measurements of dynamic moduli in parallel superposition of a steady and oscillatory shear flow were compared with different theories.Most of the rheological models which are based on the assumption that the material properties (spectrum of relaxation times) are unaffected by the presence of the steady rates of shear give results incompatible with our experimental results described earlier. From these models the best fit is given by the WJFLMB-model.A much better fit is given by theories based on the assumption that the relaxation spectrum is cut-off at a critical value of the relaxation time depending on the prescribed steady rate of shear. Under certain conditions these theories explain qualitatively the linear relations between the frequency ₀, at whichG (, ) = 0 and the steady rate of shear , as we found experimentally for several polymer solutions.

---

Zusammenfassung Die Ergebnisse von Messungen der dynamischen Moduln bei Superposition einer stationären und einer oszillierenden Scherströmungsk omponente werden mit verschiedenen Theorien verglichen.Die meisten rheologischen Modelle sind auf die Annahme gegründet, daß die Stoffeigenschaften (das Relaxationsspektrum) durch das Vorliegen einer stationären Scherströmung nicht verändert werden. Solche Theorien liefern Ergebnisse, die mit unseren früher beschriebenen Ergebnissen unvereinbar sind. Die beste Übereinstimmung liefert noch das WJFLMB-Modell.Eine wesentlich genauere Anpassung ermöglichen diejenigen Theorien, die auf der Annahme gründen, daß das Relaxationsspektrum bei einem kritischen Wert der Relaxationszeit abgeschnitten wird und daß dieser Wert von der überlagerten stationären Scherströmung abhängt. Unter gewissen Bedingungen erklären diese Theorien zumindest qualitativ die lineare Beziehung zwischen der Frequenz ₀, bei welcherG (, ) = 0 wird, und der zugeordneten Schergeschwindigkeit , wie wir sie bei verschiedenen Experimenten gefunden haben.

---

---

With 4 figures     

Nonlinear Taylor stability of viscoelastic fluids

Using Stuart's energy method, the torque on the inner cylinder, for a second order fluid, in the supercritical regime is calculated. It is found that when the second normal stress difference is negative, the flow is more stable than for a Newtonian fluid and the torque is reduced. If the second normal stress difference is positive, then the flow is more stable and there is no torque reduction. Experimental data related to the present work are discussed.Nomenclature a amplitude of the fundamentals - A _ij ⁽¹⁾ , A _ij ⁽²⁾ first and second Rivlin-Ericksen tensors - d r ₂–r ₁ - D d/dx - E - F - g _ij metric tensor - G torque on the inner cylinder in the supercritical regime - h height of the cylinders - k ₀ /d ² - k ₁ /d ² - I ₁ - I ₂ - I ₃ - I ₄ - r ₁, r ₂ radii of inner and outer cylinders respectively - r ₀ 1/2(r ₁+r ₂) - R Reynolds number ₁ r ₁ d/ ₀ - R _c critical Reynolds number - T Taylor number r _{1 1} ² d ^{3 2}/ ₀ ² *) - T _c critical Taylor number - u ₁, v ₁, w ₁ Fundamentals of the disturbance - u _i , v _i , w _i , (i>1) harmonics - mean velocity (not laminar velocity) - u –u ₁/ar _{1 1} - v v ₁/Rar _{1 1} - x (r–r ₀)/d - , material constants - 0 viscosity - wave number d - density - ₁ angular velocity of inner cylinder - tilde denotes complex conjugate     

Dimensional analysis of pore scale and field scale immiscible displacement

A basic re-examination of the traditional dimensional analysis of microscopic and macroscopic multiphase flow equations in porous media is presented. We introduce a macroscopic capillary number which differs from the usual microscopic capillary number Ca in that it depends on length scale, type of porous medium and saturation history. The macroscopic capillary number is defined as the ratio between the macroscopic viscous pressure drop and the macroscopic capillary pressure. can be related to the microscopic capillary number Ca and the LeverettJ-function. Previous dimensional analyses contain a tacit assumption which amounts to setting = 1. This fact has impeded quantitative upscaling in the past. Our definition for , however, allows for the first time a consistent comparison between macroscopic flow experiments on different length scales. Illustrative sample calculations are presented which show that the breakpoint in capillary desaturation curves for different porous media appears to occur at 1. The length scale related difference between the macroscopic capillary number for core floods and reservoir floods provides a possible explanation for the systematic difference between residual oil saturations measured in field floods as compared to laboratory experiment.     

Einflu\ der TemperaturabhÄngigkeit der Stoffwerte auf den WÄrmeübergang an der laminar überströmten ebenen Platte

Zusammenfassung Der Einflu\ der TemperaturabhÄngigkeit der Stoffwerte auf den WÄrmeübergang ist zur Zeit noch nicht befriedigend geklÄrt. Da Messungen des WÄrmeübergangs stets mit mehr oder weniger gro\en Fehlern behaftet sind, und da der Einflu\ der TemperaturabhÄngigkeit für die einzelnen Stoffwerte nicht getrennt untersucht werden kann, ist es experimentell praktisch nicht möglich, Korrekturfaktoren zu bestimmen. Es bleibt daher nur der Versuch, rechnerisch dieses Problem anzugehen, wobei es wichtig ist, Verfahren zu benutzen, die nicht zu lange Rechenzeiten erfordern, um eine gro\e Anzahl von Parameterstudien durchzuführen.über eine solche Rechnung für eine lÄngsangeströmte, ebene Platte bei laminarer Strömung wird berichtet. Dabei konnte auch der Dissipationsterm berücksichtigt werden.Bereits aus der Form der zu lösenden Differentialgleichungen ergibt sich, da\ die VerhÄltnisse der Stoffwerte einzeln berücksichtigt werden müssen, da sich die Stoffwerte in den Differentialgleichungen nicht zu einer Prandtl-Zahl zusammenfassen lassen. Die Ergebnisse der Parameterstudien werden diskutiert und mit den bisher verwendeten Korrekturfaktoren verglichen.

---

Influence of temperature dependent physical properties on the heat transfer in the laminar boundary layer of a parallel affluxed flat plate

The influence of temperature dependent properties on the heat transfer is not jet satisfying claryfied. Measurements of heat transfer contain more ore less great faults, and it is practically not possible to investigate the influence of the temperature dependence for particular physical properties separately. Therefore, the only possibility is to examine this problem numerically. Hereby it is important to use methods which don't require too much computing time in order to study a great number of parameters.In the following article such a numerical method for the calculation of a laminar boundary layer will be introduced. By this method it is also possible to consider the dissipation.Already the form of the differential equations, which are to be solved, shows, that the influence of a variable property fluid has to be considered separately, because it is not possible to combine all physical properties with the Pr-number. The results of the parameter studies will be discussed and will be compared with the well known correction terms.

---

Bezeichnungen C dimensionslose Stoffwertkorrektur - c _p spezifische WÄrmekapazitÄt - Eckertsche Kenngrö\e - dimensionslose Stromfunktion - erste Ableitung der dimensionslosen Stromfunktion - n Exponent des ZÄhigkeitstermes (Gl. (43)) - Nusseltsche Kennzahl - Pr=c_p/ Prandtlsche Kennzahl - WÄrmestromdichte - Reynoldsche Kennzahl - r Recovery-Faktor (Gl. (50)) - T Temperatur - t _E Eigentemperatur - u,u Strömungsgeschwindigkeit inx-Richtung - Strömungsgeschwindigkey-Richtung - WÄrmeübergangskoeffizient - dimensionslosey-Koordinat - dimensionslose Temperatur - Dichte - dynamische ViskositÄt - v kinematische ViskositÄt - WÄrmeleitfÄhigkeit Indizes au\erhalb der Grenzschicht - w an der Wand Vorgetragen auf der Sitzung des GVC-Fachausschusses WÄrme- und Stoffübertragung am 30. 3. 1981 in Freiburg     

Geometric Analysis of One Singular Cauchy Problem

We consider the Cauchy problem , x(0) = 0, where a ₀₀₀ = 0, a ₀₀₁ = 0, and a ₀₀₂ = 0, and prove the existence of continuously differentiable solutions x(0,] with required asymptotic properties.     

Stability of rarefaction waves in viscous media

We study the time-asymptotic behavior of weak rarefaction waves of systems of conservation laws describing one-dimensional viscous media, with strictly hyperbolic flux functions. Our main result is to show that solutions of perturbed rarefaction data converge to an approximate, Burgers rarefaction wave, for initial perturbations w ₀ with small mass and localized as w ₀(x)= The proof proceeds by iteration of a pointwise ansatz for the error, using integral representations of its various components, based on Green's functions. We estimate the Green's functions by careful use of the Hopf-Cole transformation, combined with a refined parametrix method. As a consequence of our method, we also obtain rates of decay and detailed pointwise estimates for the error.This pointwise method has been used successfully in studying stability of shock and constant-state solutions. New features in the rarefaction case are time-varying coefficients in the linearized equations and error waves of unbounded mass (log (t)). These diffusion waves have amplitude (t ^-1/2logt) in linear degenerate transversal fields and (t ^-1/2logt) in genuinely nonlinear transversal fields, a distinction which is critical in the stability proof.     

Hydro-mechanical aspects of the sand production problem

This paper examines the hydro-mechanical aspect of the sand production problem and sets the basic frame of the corresponding mathematical modelling. Accordingly, piping and surface erosion effects are studied on the basis of mass balance and particle transport considerations as well as Darcy's law. The results show that surface erosion is accompanied by high changes of porosity and permeability close to the free surface. Quantities which can be measured in experiment, like the amount of produced solids or fluid discharge, can be used in an inverse way to determine the constitutive parameters of the problem.Notation dV Volume element - dV _s Volume of solids pt - dV _v Volume of voids - dV _ff Volume of fluid phase - dV _fs Volume of fluidized-particles - Volume of mixture - dM _s Mass of solids - dM _ff Mass of fluid phase - d _{M fs} Mass of fluidized-particles - Mass of mixture - s Density of solids - f Density of fluid - ff Density of fluid phase - fs Density of fluidized-particles - Density of mixture - _i ^ff Velocity of fluid - _i ^fs Velocity of fluidized-particles - _i ^s Velocity of solids - Velocity of mixture - q ^ff Volume-discharge of fluid - q ^fs Volume-discharge of fluidized-particles - Volume-discharge of mixture - m ^ff Mass-discharge of fluid - m ^fs Mass-discharge of fluidized-particles - Mass-discharge of mixture - _er Rate of mass-eroded - _dep Rate of mass-deposited - Mass generation term - dS _i Surface element - Pore-surface element - D _IJ Tensor of mechanical dispersion - x _i Location - t Time - Porosity - c Transport concentration - c _cr Critical value ofc - p Fluid-pressure - k Permeability coefficient - _k Kinematic viscosity - Spatial frequency of erosion starter points     

The Most Unstable Profiles of a Plane-Parallel Flow in a Channel

It is well known (after Rayleigh) that a plan-parallel flow in a channel can be unstable only if the basic velocity profile U(z) possesses inflection points. The profile determines (via the Rayleigh equation) the maximal increment c _i of small perturbations and 'eigenvalues' c (see Equation (1)). The increment and the imaginary parts c _i of the eigenvalues c provide a quantitative characterization of the basic stability properties of the flow. Here we find some best possible bounds for these values. The bounds are determined by the following parameters: wave number ; enstrophy of the basic flow; width of the channel L. A similar approach can be applied to models of atmosphere, ocean, plasma etc.