An exact solution for non-darcy free convection from a horizontal line source of heat

In this paper we consider the free convection from a horizontal line source of heat which is embedded in an unbounded porous medium saturated with a fluid at rest under gravity. The convective fluid and the porous medium are in local thermal equilibrium.

---

Eine exakte Lösung der nicht-darcy'schen freien Konvektion von einer horizontalen, linienförmigen Wärmequelle

Zusammenfassung In dem Aufsatz wird die freie Konvektion von einer horizontalen, linienförmigen Wärmequelle untersucht, die in ein unbegrenztes poröses Medium eingebettet ist. Die Poren des porösen Mediums sind mit einem Fluid gefüllt, das unter Schwerkrafteinfluß ruht. Das strömende Fluid und das poröse Medium sind örtlich im thermischen Gleichgewicht.

---

Nomenclature c _p specific heat of convective fluid - F _o parameter,=/(vl>g - g acceleration due to gravity - k thermal conductivity of the saturated porous medium - l typical length scale of body - Q heat flux per unit length of a line source - Ra Rayleigh number, =gQl/²c_p - Ra _x local Rayleigh number, =xg Qx/ a²c_p - T temperature - T temperature of ambient fluid - u, x andy components of velocity - x, y coordinates vertically upwards and normal to axis of plume - X, Y non-dimensional coordinates vertically upwards and normal to axis of plume Greek symbols equivalent themal diffusivity - coefficient of thermal expansion - similarity variable - non-dimensional temperature - x permeability of porous medium - viscosity of convective fluid - v kinematic viscosity of convective fluid - density of convective fluid - stream function - non-dimensional stream function - the Forchheimer's coefficient     

Wave‐field structure in active bubble systems in shock tubes with “discontinuities” in cross section

The wave structure in active bubble media in shock tubes with sudden changes of profiles in the form of discontinuities in cross section and a onephase liquid waveguide is analyzed numerically. In axisymmetric formulation, the paper studies wave amplification due to reflection from a wall and focusing at the buttend of a rigid rod aligned coaxially with the channel. In this configuration, the amplification effect results from twodimensional cumulation of the shock wave after it leaves the annular channel and reaches the buttend of the rod. A Mach configuration forms in the focus spot. The geometrical characteristics of the shock tube allow one to control (to some extent) the amplification coefficient and the coordinates of the focus spot. In particular, it is shown that the wave can be focused near the second discontinuity of cross section — a rigid wall (in the region of passage through the interface to the waveguide) — and intensified upon reflection. If the waveguide radius is equal to the height of the Mach stem, the reflected wave has a maximum amplitude.     

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.     

Asymptotic Modeling of Nonlinear Wave Processes in Shock‐Loaded Elastoplastic Materials

Nonlinear wave processes in shockloaded elastoplastic materials are modeled. A comparison of the results obtained with experimental data and numerical solutions of exact systems of dynamic equations shows that the model equations proposed qualitatively describe the stressdistribution evolution in both the elasticflow and plasticflow regions and can be used to solve one and twodimensional problems of pulsed deformation and fracture of elastoplastic media.     

Symmetry-Breaking Bifurcations of Charged Drops

It has been observed experimentally that an electrically charged spherical drop of a conducting fluid becomes nonspherical (in fact, a spheroid) when a dimensionless number X inversely proportional to the surface tension coefficient is larger than some critical value (i.e., when <_c). In this paper we prove that bifurcation branches of nonspherical shapes originate from each of a sequence of surface-tension coefficients ), where ₂=_c. We further prove that the spherical drop is stable for any >₂, that is, the solution to the system of fluid equations coupled with the equation for the electrostatic potential created by the charged drop converges to the spherical solution as t provided the initial drop is nearly spherical. We finally show that the part of the bifurcation branch at =₂ which gives rise to oblate spheroids is linearly stable, whereas the part of the branch corresponding to prolate spheroids is linearly unstable.     

Asymptotic behavior of solutions of nonlinear elliptic equations

We study and obtain formulas for the asymptotic behavior as ¦x¦ of C ² solutions of the semilinear equation u=f(x, u), x (*) where is the complement of some ball in ⁿ and f is continuous and nonlinear in u. If, for large x, f is nearly radially symmetric in x, we give conditions under which each positive solution of (*) is asymptotic, as ¦x¦, to some radially symmetric function. Our results can also be useful when f is only bounded above or below by a function which is radially symmetric in x or when the solution oscillates in sign. Examples when f has power-like growth or exponential growth in the variables x and u usefully illustrate our results.     

Eine analytische Näherungslösung für die eingefrorene laminare Grenzschichtströmung eines binären Gemisches

Zusammenfassung Für die eingefrorene laminare Grenzschichtströmung eines teilweise dissoziierten binären Gemisches entlang einer stark gekühlten ebenen Platte wird eine analytische Näherungslösung angegeben. Danach läßt sich die Wandkonzentration als universelle Funktion der Damköhler-Zahl der Oberflächenreaktion angeben. Für das analytisch darstellbare Konzentrationsprofil stellt die Damköhler-Zahl den Formparameter dar. Die Wärmestromdichte an der Wand bestehend aus einem Wärmeleitungs- und einem Diffusionsanteil wird angegeben und diskutiert. Das Verhältnis beider Anteile läßt sich bei gegebenen Randbedingungen als Funktion der Damköhler-Zahl ausdrücken.

---

An analytical approximation for the frozen laminar boundary layer flow of a binary mixture

An analytical approximation is derived for the frozen laminar boundary layer flow of a partially dissociated binary mixture along a strongly cooled flat plate. The concentration at the wall is shown to be a universal function of the Damkohler-number for the wall reaction. The Damkohlernumber also serves as a parameter of shape for the concentration profile which is presented in analytical form. The heat transfer at the wall depending on a conduction and a diffusion flux is derived and discussed. The ratio of these fluxes is expressed as a function of the Damkohler-number if the boundary conditions are known.

---

Formelzeichen A Atom - A₂ Molekül - C Konstante in Gl. (20) - c₁=1/(2C) Konstante in Gl. (35) - c_p spezifische Wärme bei konstantem Druck - D binärer Diffusionskoeffizient - Ec=u ² /(2h_f) Eckert-Zahl - h spezifische Enthalpie - h_t=h+u²/2 totale spezifische Enthalpie - h _A ⁰ spezifische Dissoziationsenthalpie - K_w Reaktionsgeschwindigkeitskonstante der heterogenen Wandreaktion - 1= /( ) Champman-Rubesin-Parameter - Le=Pr/Sc Lewis-Zahl - M Molmasse - p statischer Druck - Pr= c_pf/ Prandtl-Zahl - q_w Wärmestromdichte an der Wand - q_cw, q_dw Wärmeleitungsbzw. Diffusionsanteil der Wärmestromdichte an der Wand - universelle Gaskonstante - R=/(2M_a) individuelle Gaskonstante der molekularen Komponente - Re_x= u x/ Reynolds-Zahl - Sc=/( D) Schmidt-Zahl - T absolute Temperatur - T_d=h _A ⁰ /R charakteristische Dissoziationstemperatur - u, v x- und y-Komponenten der Geschwindigkeit - U=u/u normierte x-Komponente der Geschwindigkeit - x, y Koordinaten parallel und senkrecht zur Platte Griechische Symbole - =A/ Dissoziationsgrad - Grenzschichtdicke - ₂ Impulsverlustdicke - Damköhler-Zahl der Oberflächenreaktion - =T/T normierte Temperatur - =y/ normierter Wandabstand - Wärmeleitfähigkeit - dynamische Viskosität - , ^* Ähnlichkeitskoordinaten - Dichte - Schubspannung Indizes A auf ein Atom bezogen - M auf ein Molekül bezogen - f auf den eingefrorenen Zustand bezogen - w auf die Wand bezogen - auf den Außenrand der Grenzschicht bezogen     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

The Thermodynamic Significance of the Local Volume Averaged Temperature

Since the temperature is not an additive function, the traditional thermodynamic point of view suggests that the volume integral of the temperature has no precise physical meaning. This observation conflicts with the customary analysis of non-isothermal catalytic reactors, heat pipes, driers, geothermal processes, etc., in which the volume averaged temperature plays a crucial role. In this paper we identify the thermodynamic significance of the volume averaged temperature in terms of a simple two-phase heat transfer process. Given the internal energy as a function of the point temperature and the density

Representation of the rheological properties of polymer melts in terms of complex fluidity

In dynamic rheological experiments melt behavior is usually expressed in terms of complex viscosity ^* () or complex modulusG ^* (). In contrast, we attempted to use the complex fluidity ^* () = 1/µ ^* () to represent this behavior. The main interest is to simplify the complex-plane diagram and to simplify the determination of fundamental parameters such as the Newtonian viscosity or the parameter of relaxation-time distribution when a Cole-Cole type distribution can be applied. ^* () complex shear viscosity - () real part of the complex viscosity - () imaginary part of the complex viscosity - G ^* () complex shear modulus - G() storage modulus in shear - G() loss modulus in shear - J ^* () complex shear compliance - J() storage compliance in shear - J() loss compliance in shear - shear strain - rate of strain - angular frequency (rad/s) - shear stress - loss angle - ^* () complex shear fluidity - () real part of the complex fluidity - () imaginary part of the complex fluidity - ₀ zero-viscosity - ₀ average relaxation time - h parameter of relaxation-time distribution     

Numerical analysis of axisymmetric buckling of conical shells

Nonlinear boundaryvalue problems of axisymmetric buckling of conical shells under a uniformly distributed normal pressure are solved by the shooting method. The problems are formulated for a system of six firstorder ordinary differential equations with independent rotation and displacement fields. Simply supported and clamped cases are considered. Branching solutions of the boundaryvalue problems are studied for different pressures and geometrical parameters of the shells. The nonmonotonic and discontinuous curves of equilibrium states obtained show that collapse, i.e., snapthrough instability is possible. For a simply supported shell, multivalued solutions are obtained for both external and internal pressure. For a clamped thinwalled shell, theoretical results are compared with experimental data.     

Qualitative behavior of dissipative wave equations on bounded domains

The qualitative behavior of solutions of the mixed problem u_tt = u-a(x)u_t in IR x , u=0 on IR x , is studied in the case when a>0 and IRⁿ is bounded. Roughly speaking, if aa_min>0, then solutions decay at least as fast as exp t( –1/2a_min), with the possible exception of a finite dimensional set of smooth solutions whose existence is associated with a phenomenon of overdamping. If a_max is sufficiently small, depending on , then no overdamping occurs.Partially supported by NSF grant NSF GP 34260.This work was partially supported by the National Science Foundation under Grant No. GP 34260     

Analytical Solution for the Effect of Viscous Dissipation on Mixed Convection in Saturated Porous Media

A new analytical solution is introduced for the effect of viscous dissipation on mixed convection flow and heat transfer about an isothermal vertical wall embedded in Darcy and non-Darcy porous media with uniform free stream velocity. The effect of viscous dissipation on mixed convection in both regimes has been analyzed for both the aiding and opposing flows using Gebhart number, Ge _x =gx/c _p. The governing parameters are Re, Ra, Pe and Ge _x . The case of Re=0 corresponds to Darcy mixed convection region and Re/Pe is identified as the mixed convection governing parameter, Ra=0 leading to pure forced convection. A good agreement was found between the numerical and analytical solutions. It was found from the Nusselt number results that viscous dissipation lowers the heat transfer rate in both Darcy and Forchheimer flow regimes for aiding as well as opposing flows.     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

Averaging of a Layered Elastic Medium with Low Dynamic Dissipation at the Interlayer Boundary

Averaged relations for a layered medium with dynamic dissipation at the interlayer boundary are constructed for dynamic problems of longitudinal shear and twodimensional theory of elasticity. For low dissipation, the averaged boundaryvalue problem is shown to be singularly perturbed. The degeneration of the boundaryvalue problem is studied qualitatively.     

Exact Solutions of the Hydrodynamic Equations Derived from Partially Invariant Solutions

The paper proposes a heuristic approach to constructing exact solutions of the hydrodynamic equations based on the specificity of these equations. A number of systems of hydrodynamic equations possess the following structure: they contain a reduced system of n equations and an additional equation for an extra function w. In this case, the reduced system, in which w = 0, admits a Lie group G. Taking a certain partially invariant solution of the reduced system with respect to this group as a seed:rdquo; solution, we can find a solution of the entire system, in which the functional dependence of the invariant part of the seed solution on the invariants of the group G has the previous form. Implementation of the algorithm proposed is exemplified by constructing new exact solutions of the equations of rotationally symmetric motion of an ideal incompressible liquid and the equations of concentrational convection in a plane boundary layer and thermal convection in a rotating layer of a viscous liquid.     

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     

A general model for the effective viscosity of pseudoplastic and dilatant fluids

Summary A new and very general expression is proposed for correlation of data for the effective viscosity of pseudoplastic and dilatant fluids as a function of the shear stress. Most of the models which have been proposed previously are shown to be special cases of this expression. A straightforward procedure is outlined for evaluation of the arbitrary constants.

---

Zusammenfassung Eine neue und sehr allgemeine Formel wird für die Korrelation der Werte der effektiven Viskosität von strukturviskosen und dilatanten Flüssigkeiten in Abhängigkeit von der Schubspannung vorgeschlagen. Die meisten schon früher vorgeschlagenen Methoden werden hier als Spezialfälle dieser Gleichung gezeigt. Ein einfaches Verfahren für die Auswertung der willkürlichen Konstanten wird beschrieben.

---

Nomenclature b arbitrary constant inSisko model (eq. [5]) - n arbitrary exponent in eq. [1] - x independent variable - y(x) dependent variable - y ₀(x) limiting behavior of dependent variable asx 0 - y(x) limiting behavior of dependent variable asx - z original dependent variable - arbitrary constant inSisko model (eq. [5]) andBird-Sisko model (eq. [6]) - arbitrary exponent in eqs. [2] and [8] - effective viscosity = shear stress/rate of shear - _A effective viscosity at = _A - _B empirical constant in eqs. [2] and [8] - ₀ limiting value of effective viscosity as 0 - ₀() limiting behavior of effective viscosity as 0 - limiting value of effective viscosity as - () limiting behavior of effective viscosity as - rate of shear - arbitrary constant inBird-Sisko model (eq.[6]) - shear stress - _A arbitrary constant in eqs. [2] and [8] - ₀ shear stress at inBingham model - _1/2 shear stress at = ( ₀ + )/2 With 8 figures     

New integral estimates for deformations in terms of their nonlinear strains

If u is a bi-Lipschitzian deformation of a bounded Lipschitz domain in ⁿ (n2), we show that the L ^P norm (p1, pn) of a certain nonlinear strain function e(u) associated with u dominates the distance in L ^q (q= np/(n–p) if p if p>n) from u to a suitably chosen rigid motion of ⁿ . This work extends that of F. John, who proved corresponding estimates for p}>1 under the hypothesis that u has uniformly small strain. We also obtain a bound for the oscillation of Du in L ². These estimates are apparently the first to apply with no a priori pointwise hypotheses upon the strain of u. In ³ the integral e(u) ² d³ is dominated by typical hyperelastic energy functionals proposed in the literature for modeling the behavior of rubber; thus the case n=3, p=2 gives the first general bound for the deformations of such materials in terms of the associated nonlinear elastic work.     

Similar solutions and the asymptotics of filtration equations

In this paper we consider the asymptotic behavior of solutions of the quasilinear equation of filtration as t. We prove that similar solutions of the equation u _t = (u )_xx asymptotically represent solutions of the Cauchy problem for the full equation u _t = [(u)]_xx if (u) is close to u for small u.