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     

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.     

On the Ellipticity of the Middle Surface of a Shell and its Application to the Asymptotic Analysis of ‘Membrane’ Shells

We consider a surface S = (), where ² is a bounded, connected, open set with a smooth boundary and : ³ is a smooth map; let () denote the components of the two-dimensional linearized strain tensor of S and let 0 with length 0 > 0. We assume the the norm ,|| ()||0, in the space V0() = { H¹() × H¹() × L²(); = 0 on 0 } is equivalent to the usual product norm on this space. We then establish that this assumption implies that the surface S is uniformly elliptic and that we necessarily have 0 = .     

On shock layer method for the hypersonic flow problems

Chernyi’s series method[1] is not proper for the case that(γ-l)/(γ+l)<<2/(γ+1)×M²sin²β (γ=c_p/c_v-adiabatic index number, M-Much number, β-shock incidence). In this paper, we only suppose that in the neighbour of the shock, there exists a shock layer in which the density of the gas is very big, but we do not remove the case that (γ-1)/(γ+1)<<2/(γ+1)M²sin²β.     

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.     

On forced vibration of anisotropic cylinders

Summary The dynamic response of a circular cylinder with thick walls of transverse curvilinear isotropy subjected to a uniformly distributed pressure varying periodically with time is analyzed by means of the Laplace transformation, and the exact solution is obtained in closed form. The previously obtained solutions for forced vibrations with isotropy, and free vibrations with transverse curvilinear isotropy are included as special cases of the general results reported here.Nomenclature t time - r, , z cylindrical coordinates - _ii components of normal strain - _ii components of normal stress - u radial displacement - c _ij elastic constant - mass density - c ² c ₁₁/ - ² c ₂₂/c ₁₁ - a, b inner, outer radius of the cylinder - , A, B constants - forced angular frequency - function defined by (9) - p, real, complex variables - constant defined by (14) - real number - , Lamé elastic constants - J (x) Bessel function of first kind - Y (x) Bessel function of second kind - I (x) modified Bessel function of first kind - K (x) modified Bessel function of second kind     

Construction of exact numerical solutions of the stationary traveling wave type for viscous thin films

An effective numerical procedure, based on the Galerkin method, for finding solutions of the stationary traveling wave type in the complete formulation is proposed for the case of viscous liquid films. Examples of a viscous film flowing freely down a vertical surface have been calculated. The calculations have been made for various values of the dimensionless surface tension , including =0. The method makes it possible to predict a number of bifurcations that occur as decreases. The existence of numerous families of stationary traveling waves when 1 was demonstrated in [6]. The present study shows that as 1 all but one of these families of wave solutions disappear. The shape of the periodic and solitary waves and the pressure distribution in the film are found for various . When =0 and the wave number is fairly small, the periodic solution has a singularity, as predicted in [14]: at the crest of the wave a corner point appears; the angle between the tangents at this point =140–150. The method proposed can be used to calculate other wavy film flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 94–100, May–June, 1990.     

Transport in ordered and disordered porous media IV: Computer generated porous media for three-dimensional systems

In the method of volume averaging, the difference between ordered and disordered porous media appears at two distinct points in the analysis, i.e. in the process of spatial smoothing and in the closure problem. In theclosure problem, the use of spatially periodic boundary conditions isconsistent with ordered porous media and the fields under consideration when the length-scale constraint,r ₀L is satisfied. For disordered porous media, spatially periodic boundary conditions are an approximation in need of further study.In theprocess of spatial smoothing, average quantities must be removed from area and volume integrals in order to extractlocal transport equations fromnonlocal equations. This leads to a series of geometrical integrals that need to be evaluated. In Part II we indicated that these integrals were constants for ordered porous media provided that the weighting function used in the averaging process contained thecellular average. We also indicated that these integrals were constrained by certain order of magnitude estimates for disordered porous media. In this paper we verify these characteristics of the geometrical integrals, and we examine their values for pseudo-periodic and uniformly random systems through the use of computer generated porous media.

Nomenclature

Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - a _i i=1, 2, 3 gaussian probability distribution used to locate the position of particles - I unit tensor - L general characteristic length for volume averaged quantities, m - L characteristic length for , m - L characteristic length for , m - characteristic length for the -phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1, 2, 3 lattice vectors, m - m convolution product weighting function - m _v special convolution product weighting function associated with the traditional volume average - n _i i=1, 2, 3 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - r position vector, m - r _m support of the weighting functionm, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume,, m³ - x positional vector locating the centroid of an averaging volume, m - x ₀ reference position vector associated with the centroid of an averaging volume, m - y position vector locating points relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - /L, small parameter in the method of spatial homogenization - standard deviation ofa _i - _r standard deviation ofr - _r intrinsic phase average of      

Shear softening and thixotropic properties of wheat flour doughs in dynamic testing at high shear strain

Shear softening and thixotropic properties of wheat flour doughs are demonstrated in dynamic testing with a constant stress rheometer. This behaviour appears beyond the strictly linear domain (strain amplitude ₀ 0.2%),G,G and |^*| decreasing with ₀, the strain response to a sine stress wave yet retaining a sinusoidal shape. It is also shown thatG recovers progressively in function of rest time. In this domain, as well as in the strictly linear domain, the Cox-Merz rule did not apply but() and | ^*())| may be superimposed by using a shift factor, its value decreasing in the former domain when ₀ increases. Beyond a strain amplitude of about 10–20%, the strain response is progressively distorted and the shear softening effects become irreversible following rest.     

The dynamic performance of the Weissenberg rheogoniometer III. Large amplitude oscillatory shearing — harmonic analysis

The harmonic content of the nonlinear dynamic behaviour of 1% polyacrylamide in 50% glycerol/water was studied using a standard Model R 18 Weissenberg Rheogoniometer. The Fourier analysis of the Oscillation Input and Torsion Head motions was performed using a Digital Transfer Function Analyser.In the absence of fluid inertia effects and when the amplitude of the (fundamental) Oscillation Input motion _I is much greater than the amplitudes of the Fourier components of the Torsion Head motion _Tn empirical nonlinear dynamic rheological propertiesG _n (, ⁰),G _n (, ⁰) and/or _n (, ⁰), _n (, ⁰) may be evaluated without a-priori-knowledge of a rheological constitutive equation. A detailed derivation of the basic equations involved is presented.Cone and plate data for the third harmonic storage modulus (dynamic rigidity)G ₃ (, ⁰), loss modulusG ₃ (, ⁰) and loss angle ₃ (, ⁰) are presented for the frequency range 3.14 × 10^–2 1.25 × 10² rad/s at two strain amplitudes, _CP ⁰ = 2.27 and 4.03. Composite cone and plate and parallel plates data for both the third and fifth harmonic dynamic viscosities ₃ (, ⁰), _S (, ⁰) and dynamic rigiditiesG ₃ (, ⁰),G ₅ (, ⁰) are presented for strain amplitudes in the ranges 1.10 _CP ⁰ 4.03 and 1.80 _PP ⁰ 36 for a single frequency, = 3.14 × 10^–1 rad/s. Good agreement was obtained between the results from both geometries and the absence of significant fluid inertia effects was confirmed by the superposition of the data for different gap widths.     

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     

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

Mass transport and reaction in catalyst pellets

When heterogeneous chemical reactions take place in porous catalysts, mass transport can occur by bulk diffusion, Knudsen diffusion, and convective transport. Previous studies of these phenomena have been largely based on Maxwell's dusty gas model with the convective transport or Darcy flow added to the diffusive transport. This is done in order to satisfy one of the limiting conditions encountered in the study of flow in porous media. A more fundamental approach consists of the use of the method of volume averaging and the general form of the species momentum equation. For an N-component system, this leads to N independent flux relations to be used in conjunction with the volume-averaged species continuity equations.Roman Letters _A (t) surface area of a species body, m² - a _v interfacial area per unit volume, m^-1 - A _e area of entrances and exits for the -phase contained within the averaging volume, m² - A _K area of the - surface contained within the averaging volume, m² - b _A species A body force, N/kg - b mass average body force, N/kg - B inverse tortuosity tensor for bulk diffusion - c total molar concentration, moles/m³ - c _A species A molar concentration, moles/m³ - _A surface concentration of species A, moles/m² - C_A² intrinsic phase average molar concentration, moles/m³ - c _A – C_A², spatial deviation concentration, moles/m³ - c _A mean molecular speed for species A, m/s - binary diffusion coefficient, m²/s - D _A ^{K, eff} Knudsen diffusion coefficient for species A, m²/s - f vector that maps P _A into P _A , m - g gravitational vector, m/s² - G second order tensor that maps N _A into N _A for free molecule flow conditions - H inverse tortuosity tensor for Knudsen diffusion - I unit tensor - j _A c _A u _A ^* , molar diffusive flux, moles/m²s - K Darcy's Law permeability tensor, m² - L macroscopic length scale, m - L _D diffusive length, m - l characteristic length for the -phase, m - l _A mean free path for species A, m - M _A molecular weight of species A, kg/mole - n outwardly directed unit normal vector - n _K unit normal vector directed from the -phase toward the -phase - n outwardly directed unit normal vector at the entrances and exits of the -phase contained within the averaging volume - N _A c _A v _A molar flux of species A, moles/m²s - N _A intrinsic phase average of the species A molar flux, moles/m²s - \~N _A spatial deviation of the molar flux of species A, moles/m²s - p total pressure, N/m² - P p + , total pressure over and above the hydrostatic pressure, N/m² - P _A partial pressure of species A, N/m² - p _A intrinsic phase average partial pressure, N/m² - P_A – p _A, spatial deviation partial pressure, N/m² - P _A p_A + _AA partial pressure of species A over and above the hydrostatic pressure of species A, N/m² - p _ab diffusive force exerted by species B on species A, N/m³ - universal gas constant, N m/moles K - R _A molar rate of production of species A owing to homogeneous chemical reaction, moles/m³s - molar rate of production of species A owing to heterogeneous chemical reaction, moles/m²s - r _A mass rate of production of species A owing to homogeneous chemical reaction, kg/m³s - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - t _A species stress vector, N/m² - T _A species stress tensor, N/m² - T total stress tensor, N/m² - T temperature, K - T spatial average temperature, K - u _A v _A – v, mass diffusion velocity, m/su _A ^* v_A – v^*, molar diffusion velocity, m/s - u _o velocity of the rigid, solid phase relative to some inertial frame, m/s - v _A species velocity, m/s - v mass average velocity, m/s - v ^* molar average velocity, m/s - v _A ^* species velocity of those molecules of species A generated by chemical reaction, m/s - _A (t) volume of a species A body, m³ - averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ - v phase average, mass average velocity, m/s - w arbitrary velocity vector, m/s - x _A c _A /c mole fraction of species A - X _A intrinsic phase average mole fraction - X _A – X _A , spatial deviation mole fraction Greek Letters V/V volume fraction of the -phase - _A sum of all terms in the species A momentum equation that are small compared to the diffusive force, N/m³ - viscosity of the -phase, Ns/m² - _A mass density of species A, kg/m³ - total mass density, kg/m³ - _a species viscous stress tensor, N/m² - total viscous stress tensor, N/m² - tortuosity factor - total body force potential function, Nm/kg - _a species body force potential function, Nm/kg - 3.1416 - _{a a} / mass fraction of species A     

Determination of permeability tensors for two-phase flow in homogeneous porous media: Theory

In this paper we continue previous studies of the closure problem for two-phase flow in homogeneous porous media, and we show how the closure problem can be transformed to a pair of Stokes-like boundary-value problems in terms of pressures that have units of length and velocities that have units of length squared. These are essentially geometrical boundary value problems that are used to calculate the four permeability tensors that appear in the volume averaged Stokes' equations. To determine the geometry associated with the closure problem, one needs to solve the physical problem; however, the closure problem can be solved using the same algorithm used to solve the physical problem, thus the entire procedure can be accomplished with a single numerical code.Nomenclature a a vector that maps V onto , m^-1. - A a tensor that maps V onto . - A area of the - interface contained within the macroscopic region, m². - A area of the -phase entrances and exits contained within the macroscopic region, m². - A area of the - interface contained within the averaging volume, m². - A area of the -phase entrances and exits contained within the averaging volume, m². - Bo Bond number (= (=(–)g²/). - Ca capillary number (= v/). - g gravitational acceleration, m/s². - H mean curvature, m^-1. - I unit tensor. - permeability tensor for the -phase, m². - viscous drag tensor that maps V onto V. - ^* dominant permeability tensor that maps onto v , m². - ^* coupling permeability tensor that maps onto v , m². - characteristic length scale for the -phase, m. - l characteristic length scale representing both and , m. - L characteristic length scale for volume averaged quantities, m. - n unit normal vector directed from the -phase toward the -phase. - n unit normal vector representing both n and n . - n unit normal vector representing both n and n . - P pressure in the -phase, N/m². - p superficial average pressure in the -phase, N/m². - p intrinsic average pressure in the -phase, N/m². - p –p , spatial deviation pressure for the -phase, N/m². - r ₀ radius of the averaging volume, m. - r position vector, m. - t time, s. - v fluid velocity in the -phase, m/s. - v superficial average velocity in the -phase, m/s. - v intrinsic average velocity in the -phase, m/s. - v –v , spatial deviation velocity in the -phase, m/s. - V volume of the -phase contained within the averaging volmue, m³. - averaging volume, m³. Greek Symbols V /, volume fraction of the -phase. - viscosity of the -phase, Ns/m². - density of the -phase, kg/m³. - surface tension, N/m. - (v +v ^T ), viscous stress tensor for the -phase, N/m².     

New calculations on the Faraday effect in wave guides

Summary The propagation of electromagnetic waves is investigated theoretically for a round wave guide, containing a gyroelectrie-gyromagnetic medium with gyroaxis parallel to the guide in the form of a cylindrical shell of thickness, adjacent to the wall of the guide. An equation is set up, permitting to compute the change in the propagation constant due to the presence of the shell, including terms proportional to ². Assuming only the presence of gyromagnetism, the change ₁ of first order in for TE-waves is determined and is found to be the same fpr right- and left-circular polarization. The second order difference ₂ ⁺ — ₂ ^- for the two senses of polarization, however, appears to have a non-vanishing value which, just like ₁ can be expressed in terms of the radius of the guide, the frequency, the dielectric constant and the elements of the gyromagnetic permeability tensor which characterize the medium of the shell.     

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.     

On physical criteria for the selection of wood for soundboards of musical instruments

In order to develop criteria for the physical evaluation of wood for soundboards of musical instruments, measurements were made of dynamic Young's modulusE, static Young's modulusE, internal frictionQ ^–1 in longitudinal direction, and specific gravity for numerous species of broad-leaved wood. From the results obtained, including those of our previous paper on coniferous wood [1], it was found that the suitability of wood for soundboards could be evaluated by the quantity ofQ ^–1/(E/), and that there were very high correlations betweenQ ^–1/(E/) andE/, and betweenE andE, regardless of wood species. Consequently, it becomes possible to select practically any wood suitable for soundboards by using the value ofE/, which can be measured easily, and it was derived that the relation betweenE/ andQ ^–1 of wood could be expressed by an exponential equation regardless of wood species.     

Instability of rarefaction shocks in systems of conservation laws

It is well-known that rarefaction shocks are unstable solutions of nonlinear hyperbolic conservation laws. Indeed, for scalar equations rarefaction shocks are unstable in the class of smooth solutions, but for systems one can only say in general that rarefaction shocks are unstable in the larger class of weak solutions. (Here unstable refers to a lack of continuous dependence upon perturbations of the initial data.) Since stability in the class of weak solutions is not well understood, ([T, TE]), entropy considerations have played a leading role in ruling out shocks that violate the laws of physics. However, for non-strictly hyperbolic systems the analogy with the equations of gas dynamics breaks down, and general entropy or admissibility criteria for the variety of shocks which appear, (see, e.g., [IMPT]), are not known. In this paper we address the question of when the instability of a shock can be demonstrated within the class of smooth solutions alone. We show by elementary constructions that this occurs whenever there exists an alternative solution to the Riemann problem with the same shock data which consists entirely of rarefaction waves and contact discontinuities with at least one non-zero rarefaction wave. We show that for 2×2 strictly hyperbolic, genuinely nonlinear systems the condition is both necessary and sufficient. We show too that for the full 3×3 (Euler) equations of gas dynamics with polytropic equations of state, rarefaction shocks of moderate strength are unstable in the class of smooth solutions if and only if the adiabatic gas constant satisfies 1 < < 5/3 (see Theorem 8). More precisely, there is a constant y _*, 0 < y _* < 1, depending only on , such that if y _* p _lp _rp _l for 1-shocks, and if y _* p _rP _lp _r for 3-shocks (where p _r and p _l denote the pressures on both sides of the rarefaction shock), then the shock is unstable if and only if 1 < < 5/3. Thus for such shocks, the theory of the Riemann problem for polytropic gases in the range 1 < < 5/3 can be rigorously developed with a knowledge of the smooth solutions alone by using stability under smoothing as an admissibility criterion, rather than by using the classical entropy inequalities.     

Transport in ordered and disordered porous media III: Closure and comparison between theory and experiment

In this paper we examine the closure problem associated with the volume averaged form of the Stokes equations presented in Part II. For both ordered and disordered porous media, we make use of a spatially periodic model of a porous medium. Under these circumstances the closure problem, in terms of theclosure variables, is independent of the weighting functions used in the spatial smoothing process. Comparison between theory and experiment suggests that the geometrical characteristics of the unit cell dominate the calculated value of the Darcy's law permeability tensor, whereas the periodic conditions required for thelocal form of the closure problem play only a minor role.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A interfacial area of the- interface associated with the local closure problem, m² - A _p surface area of a particle, m² - b vector used to represent the pressure deviation, m^–1 - B ⁰ B+I, a second order tensor that maps v _m ontov - B second-order tensor used to represent the velocity deviation - d _p 6V _p/A_p, effective particle diameter, m - d a vector related to the pressure, m - D a second-order tensor related to the velocity, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor calculated on the basis of a spatially periodic model, m² - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p characteristic length for the volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - characteristic length (pore scale) for the-phase - _i i=1, 2, 3 lattice vectors, m - weighting function - m(-y) , convolution product weighting function - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - n unit normal vector pointing from the-phase toward the -phase - p pressure in the-phase, N/m² - p _m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p – p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function - r position vector, m - r position vector locating points in the-phase, m. - V averaging volume, m³ - B volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v _m superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - v traditional superficial volume averaged velocity, m/s - v – v _m , spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the -phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * , weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

Effect of the Nature of the Interaction between a Gas and a Surface on the Relations for the Knudsen Layer in the Case of Strong Evaporation

The effect of the temperature accommodation coefficient _T on the relations at the Knudsen layer edge is investigated for strong evaporation using the moment method. An explicit expression for the dimensionless density as a function of the temperature and the Mach number M is obtained for 0 < _T < 1. For _T = 0 the entire solution is obtained in explicit form. It is shown that for = 0 and a condensation coefficient << 1 the temperature outside the Knudsen layer changes sharply as M varies from 0 to a certain value much less than unity after which the temperature ceases to depend on . For the model of specular reflection of the molecules from the surface the density and the temperature outside the Knudsen layer are found in explicit form as functions of the Mach number.