Free Convection from a Vertical Plate in a Porous Medium Subjected to a Porous Medium Subjected to a Sudden Change in Surface Heat Flux

An analysis is made of the transient free convection from a vertical flat plate which is embedded in a fluid-saturated porous medium. It is assumed that for time a steady state temperature or velocity has been obtained in the boundary-layer which occurs due to a uniform flux dissipation rate . Then at time the heat flux on the plate is suddenly changed to and maintained at this value for 0$$ " align="middle" border="0"> . An analytical solution has been obtained for the temperature/velocity field for small times in which the transport effects are confined within an inner layer adjacent to the plate. These effects cause a new steady boundary layer. A numerical solution of the full boundary-layer equations is then obtained for the whole transient from to the steady state, firstly by means of a step-by-step method and then by a matching technique. The transition between the two distinct solution methods is always observed to occur very near to the turning point of the plate surface temperature, a time at which the fluid temperature is close to its steady state profile. The solution obtained using the step-by-step method shows excellent agreement with the small time analytical solution. Results are presented to illustrate the occurrence of transients from both small and large increases and decreases in the levels of existing energy inputs.     

An Approximate Solution for Fluid Infiltration near a Circular Opening in Unsaturated Media

A regular perturbation technique is employed to approximate the solution for fluid infiltration from a circular opening into an unsaturated medium. Introducing two empirical constitutive relations and relating the permeability k and water content with pore fluid pressure p, a nonlinear diffusion equation in terms of pore pressure is established. After rearranging the nonlinear diffusion equation, a parameter perturbation on is performed and an approximate solution with an error of is obtained, which correlate to a condition in which = . This approximate solution is verified by a finite difference solution and compared also with a linear solution in which the diffusivity is constant. It is shown that the perturbation solution with terms up to and including first-order can give a reasonably accurate solution for the parameter range for p ₀ selected in this paper. The solution procedure provided in this paper also avoids the numerical problem normally encountered for a small time solution. The solution may also be used to overcome difficulties arising in solution procedure by the similarity transformation (Boltzmann), commonly conducted on diffusion equation, which cannot be applied for a finite wellbore problem.     

On the universality of the velocity profiles of a turbulent flow in an axially rotating pipe

If a fluid enters an axially rotating pipe, it receives a tangential component of velocity from the moving wall, and the flow pattern change according to the rotational speed. A flow relaminarization is set up by an increase in the rotational speed of the pipe. It will be shown that the tangential- and the axial velocity distribution adopt a quite universal shape in the case of fully developed flow for a fixed value of a new defined rotation parameter. By taking into account the universal character of the velocity profiles, a formula is derived for describing the velocity distribution in an axially rotating pipe. The resulting velocity profiles are compared with measurements of Reich [10] and generally good agreement is found.Nomenclature b constant, equation (34) - D pipe diameter - l mixing length - l ₀ mixing length in a non-rotating pipe - N rotation rate,N=Re /Re _D - p pressure - R pipe radius - Re _D flow-rate Reynolds number, - Re rotational Reynolds number, Re =v _w D/ - Re* Reynolds number based on the friction velocity, Re*=v*R/ - (Re*)₀ Reynolds number based on the friction velocity in a non-rotating pipe - Ri Richardson number, equation (10) - r coordinate in radial direction - dimensionless coordinate in radial direction, - v _r ,v ,v _z time mean velocity components - v _r ,v ,v _z velocity fluctations - v _w tangential velocity of the pipe wall - v* friction velocity, - axial mean velocity - v _ZM maximum axial velocity - dimensionless radial distance from pipe wall, - y ⁺ dimensionless radial distance from pipe wall - y ₁ ⁺ constant - Z rotation parameter,Z =v _w/v _* =N Re _D /2Re* - _m eddy viscosity - ( _m )₀ eddy viscosity in a non-rotating pipe - coefficient of friction loss - von Karman constant - ₁ constant, equation (31) - density - dynamic viscosity - kinematic viscosity     

The Experimental Approach to Effective Stress Law of Coal Mass by Effect of Methane

Effective stress law of all kinds of coal samples, including steam coal, fat coal, corking coal, thin coal and anthracite, under pore pressure of gas, is experimentally studied using a newly developed test machine. These samples are taken from Coal Mines in Wuda, Hebi, Yanzhou, Yangquan, Qingshui, and Gujiao in China. The experiment results show that, under pore pressure of gas, the tested coal samples comply with Biots effective stress law, where the Biots coefficient is not a constant, and is bilinear function of volumetric stress () and pore pressure (p), that is, We define four areas according to the numerical feature of , that is, functionless area of pore pressure, normal function area, fracturing function area, and quasi-soil function area. The effective stress law of coal mass introduced by this paper is a constitutive equation in the study of coupled solid and fluid. This has significance in the drainage and outburst of methane in coal seam.     

On the continuum modelling of porous media containing fluid: a molecular viewpoint with particular attention to scale

Mass conservation and linear momentum balance relations for a porous body and any fluid therein, valid at any given length scale in excess of nearest-neighbour molecular separations, are established in terms of local weighted averages of molecular quantities. The mass density field for the porous body at a given scale is used to identify its boundary at this scale, and a porosity field is defined for any pair of distinct length scales. Specific care is paid to the interpretation of the stress tensor associated with each of the body and fluid at macroscopic scales, and of the force per unit volume each exerts on the other. Consequences for the usual microscopic and macroscopic viewpoints are explored.Nomenclature material system; Section 2.1. - porous body (example of a material system); Sections 2.1, 3.1, 4.1 - fluid body (example of a material system); Sections 2.1, 3.1, 4.1 - weighting function; Sections 2.1, 2.3 - ,h weighting function corresponding to spherical averaging regions of radius and boundary mollifying layer of thicknessh; Section 3.2 - Euclidean space; Section 2.1 - V space of all displacements between pairs of points in; Section 2.1 - mass density field corresponding to; (2.3)₁ - ^P , ^f mass density fields for , ; (4.1) - P momentum density field corresponding to; (2.3)₂ - v velocity field corresponding to; (2.4) - S _r (X) interior of sphere of radiusr with centre at pointx; (3.3) - boundary ofany region - region in which ^p > 0 with = ,h; (3.1) - subset of whose points lie at least+h from boundary of ; (3.4) - abbreviated versions of ; Section 3.2, Remark 4 - strict interior of ; (3.7) - analogues of for fluid system ; Section 3.2 - general version of corresponding to any choice of weighting function; (4.6) - interfacial region at scale; (3.8) - ₀ scale of nearest-neighbour separations in ; Section 3.2. Remark 1 - porosity field at scales ( ₁; ₂); (3.9) - pore space at scales ( ₁; ₂); (3.12)     

Flow in Random Porous Media

Flow in a porous medium with a random hydraulic conductivity tensor K(x) is analyzed when the mean conductivity tensor (x) is a non-constant function of position x. The results are a non-local expression for the mean flux vector (x) in terms of the gradient of the mean hydraulic head (x), an integrodifferential equation for (x), and expressions for the two point covariance functions of q(x) and (x). When K(x) is a Gaussian random function, the joint probability distribution of the functions q(x) and (x) is determined.     

Optimal cupolas of uniform strength: Spherical M-shells and axisymmetric T-shells

Summary The first part of this paper is concerned with the optimal design of spherical cupolas obeying the von Mises yield condition. Five different load combinations, which all include selfweight, are investigated. The second part of the paper deals with the optimal quadratic meridional shape of cupolas obeying the Tresca yield condition, considering selfweight plus the weight of a non-carrying uniform cover. It is established that at long spans some non-spherical Tresca cupolas are much more economical than spherical ones.

---

Optimale Kuppeln gleicher Festigkeit: Kugelschalen und axialsymmetrische Schalen

Übersicht Im ersten Teil dieser Arbeit wird der optimale Entwurf sphärischer Kuppeln behandelt, wobei die von Misessche Fließbewegung zugrunde gelegt wird. Fünf verschiedene Lastkombinationen werden untersucht. Der zweite Teil befaßt sich mit der optimalen quadratischen Form des Meridians von Kuppeln, die der Fließbedingung von Tresca folgen.

---

List of Symbols a_k, b_k, c_k, A_k, B_k, C_k coefficients used in series solutions - A, B constants in the nondimensional equation of the meridional curve - normal component of the load per unit area of the middle surface - meridional and circumferential forces per unit width - radial pressure per unit area of the middle surface, - skin weight per unit area of the middle surface, - vertical external load per unit horizontal area, - base radius, - R radius of convergence - s - cupola thickness, - u, w subsidiary functions for quadratic cupolas - vertical component of the load per unit area of middle surface - resultant vertical force on a cupola segment - structural weight of cupola, - combined weight of cupola and skin, - distance from the axis of rotation, - vertical distance from the shell apex, - z auxiliary variable in series solutions - specific weight of structural material of cupola - radius of the middle surface, - uniaxial yield stress - meridional stress, - circumferential stress, - _a, _b, _c, _d, _e subsidiary variables used in evaluating the meridional stress - auxiliary function used in series solutions This paper constitutes the third part of a study of shell optimization which was initiated and planned by the late Prof. W. Prager     

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

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

High Reynolds Number Flow past a Plate Mounted at a Small Angle of Attack

A plane incompressible fluid flow past a plate mounted in a homogeneous stream at a small angle of attack ^* is investigated on the basis of an asymptotic analysis of the Navier-Stokes equations at high Reynolds numbers (Re). In the neighborhood of the leading edge the flow structure is studied in detail. It is found that separation is initiated in a small vicinity of the leading edge at and the length of the slow reverse stream zone is of the order O(1) at . The nonuniqueness of the solution is detected at and the hysteresis phenomenon is explained. It is shown that under certain conditions the solutions obtained also hold for flows past bodies of small thickness.     

Mathematical model of the entry length in the flow of suspensions of solid particles in a gas

A mathematical model consisting of equations of mass and momentum and for the velocity field has been used for computing the entry length of the flow of non-Newtonian fluids in laminar, transition and turbulent regions. Experimental data measured in a vertical flow of a suspension of solid particles in air have been used for verifying the predictions. n flow index for laminar flow - Re Reynolds number defined for the flow of the carrier medium - q exponent for turbulent flow - ratio of core radius with a flat velocity profile to pipe radius - _c ratio of the axial component of local velocity in the core to mean velocity - w mean flow velocity - ratio of axial distance from the pipe entrance to the pipe radius - ratio of the entrance length to the pipe radius - relative mass fraction of particles - ratio of the distance from the pipe wall to the pipe radius - coefficient of pressure loss due to friction     

The Two-Dimensional Attractor of a Differential Equation with State-Dependent Delay

The delay differential equation

Global Attractors: Topology and Finite-Dimensional Dynamics

Many dissipative evolution equations possess a global attractor with finite Hausdorff dimension d. In this paper it is shown that there is an embedding X of into , with N=[2d+2], such that X is the global attractor of some finite-dimensional system on with trivial dynamics on X. This allows the construction of a discrete dynamical system on which reproduces the dynamics of the time T map on and has an attractor within an arbitrarily small neighborhood of X. If the Hausdorff dimension is replaced by the fractal dimension, a similar construction can be shown to hold good even if one restricts to orthogonal projections rather than arbitrary embeddings.     

Strong Solutions for Generalized Newtonian Fluids

We consider the motion of a generalized Newtonian fluid, where the extra stress tensor is induced by a potential with p-structure (p = 2 corresponds to the Newtonian case). We focus on the three dimensional case with periodic boundary conditions and extend the existence result for strong solutions for small times from \tfrac{5}{3}$$ " align="middle" border="0"> (see [16]) to \tfrac{7}{5}.$$ " align="middle" border="0"> Moreover, for we improve the regularity of the velocity field and show that for all 0.$$ " align="middle" border="0"> Within this class of regularity, we prove uniqueness for all \tfrac{7}{5}.$$ " align="middle" border="0"> We generalize these results to the case when p is space and time dependent and to the system governing the flow of electrorheological fluids as long as      

Effect of a Short Roughness Strip on a Turbulent Boundary Layer

The response of a turbulent boundary layer to a short roughness strip is investigated using laser Doppler velocimetry (LDV) and laser induced fluorescence (LIF). Skin friction coefficients are inferred from accurate near-wall measurements. There is an undershoot in , where is the undisturbed smooth wall skin friction coefficient, immediately after the strip. Downstream of the strip, overshoots before relaxing back to unity in an oscillatory manner. The roughness strip has a major effect on the turbulent stresses ; these quantities increase, relative to the undisturbed smooth wall, in the region between the two internal layers originating at the upstream and downstream edges of the strip. The increase in the ratio suggests a decrease in near-wall anisotropy. From the flow visualizations, it is inferred that streamwise vortical structures are weakened immediately downstream of the strip. Consistently, streamwise length scales are also reduced; direct support for this is provided by measured two-point velocity correlations.     

Inhomogeneous ‘longitudinal’ plane waves in a deformed elastic material

By definition, a homogeneous isotropic compressible Hadamard material has the property that an infinitesimal longitudinal homogeneous plane wave may propagate in every direction when the material is maintained in a state of arbitrary finite static homogeneous deformation. Here, as regards the wave, homogeneous means that the direction of propagation of the wave is parallel to the direction of eventual attenuation; and longitudinal means that the wave is linearly polarized in a direction parallel to the direction of propagation. In other words, the displacement is of the form u = ncos k(n · x – ct), where n is a real vector. It is seen that the Hadamard material is the most general one for which a longitudinal inhomogeneous plane wave may also propagate in any direction of a predeformed body. Here, inhomogeneous means that the wave is attenuated, in a direction distinct from the direction of propagation; and longitudinal means that the wave is elliptically polarized in the plane containing these two directions, and that the ellipse of polarization is similar and similarly situated to the ellipse for which the real and imaginary parts of the complex wave vector are conjugate semi-diameters. In other words, the displacement is of the form u = {S exp i(S · x – ct)}, where S is a complex vector (or bivector). Then a Generalized Hadamard material is introduced. It is the most general homogeneous isotropic compressible material which allows the propagation of infinitesimal longitudinal inhomogeneous plane circularly polarized waves for all choices of the isotropic directional bivector. Finally, the most general forms of response functions are found for homogeneously deformed isotropic elastic materials in which longitudinal inhomogeneous plane waves may propagate with a circular polarization in each of the two planes of central circular section of the ⁿ -ellipsoid, where is the left Cauchy-Green strain tensor corresponding to the primary pure homogeneous deformation.     

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.     

Investigation of the Geometrical Dispersion Regime in a Single Fracture Using Positron Emission Projection Imaging

The transport properties of a natural fracture crossing a limestone block of 36 cm × 26 cm × 60 cm is studied using positron emission projection imaging. This non-invasive technique allows to measure the spatial distribution of the activity of a radioactive solution (here irradiated-copper-EDTA solution) within the fracture. The fracture aperture is measured from the spatial distribution of the activity as the fracture is completely filled with the tracer. The experiment consists in injecting the tracer at a constant flow rate in the plane of the fracture filled with an identical non-radioactive solution. Every 10 min, a two-dimensional grey scale image of the concentration field is recorded. The heterogeneity of the tracer distribution increases with time in relation with the spatial heterogeneity of the aperture field, and favours only slightly the region of larger aperture. The correlation length of the aperture distribution is larger than the correlation length of the concentration distribution of the tracer within the sample. Consequently, the concentration distribution cannot be modelled using a classical advection–dispersion equation; the mixing process has not reached a stationary Fickian dispersion regime in the finite size domain of the experiment. Nevertheless, the transversally averaged concentration profiles evaluated along the flow direction x rescale adequately with an advective variable , where is the mean velocity and t the time. This result is explained in the context of the geometrical dispersion regime where the mixing dispersion zone grows proportionally with time. Different approaches are proposed to characterise this anomalous dispersion regime.     

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.     

Wärmeübertragung in einem axial rotierenden,durchströmten Rohr im Bereich des thermischen Einlaufs

Zusammenfassung Der Einfluß der Rotation auf das Temperaturprofil und die Wärmeübergangszahl einer turbulenten Rohrströmung im Bereich des thermischen Einlaufs wird theoretisch untersucht und mit Meßwerten verglichen. Es wird angenommen, daß das Geschwindigkeitsprofil voll ausgebildet ist. Die Rotation hat aufgrund der radial ansteigenden Zentrifugalkräfte einen ausgeprägten Einfluß auf die Unterdrückung der turbulenten Bewegung. Dadurch verschlechtert sich die Wärmeübertragung mit steigender Rotations-Reynoldszahl und die thermische Einlauflänge nimmt beträchtlich zu.

---

Heat transfer in an axially rotating pipe in the thermal entrance region. Part 1: Effect of rotation on turbulent pipe flow

The effects of rotation on the temperature distribution and the heat transfer to a fluid flowing inside a tube are examined by analysis in the thermal entrance region. The theoretical results are compared with experimental findings. The flow is assumed to have a fully developed velocity profile. Rotation was found to have a very marked influence on the suppression of the turbulent motion because of radially growing centrifugal forces. Therefore, a remarkable decrease in heat transfer with increasing rotational Reynolds number can be observed. The thermal entrance length increases remarkably with growing rotational Reynolds number.

---

Formelzeichen a Temperaturleitzahl - C _n , ,C ₁,C ₃ Konstanten - c _p spezifische Wärme bei konstantem Druck - D Rohrdurchmesser - E Funktion nach Gl. (30) - H _n Eigenfunktionen - l hydrodynamischer Mischungsweg - l _q thermischer Mischungsweg - Massenstrom - N=Re /Re Reynoldszahlenverhältnis - Nu Nusseltzahl - Nu Nusseltzahl für die thermisch voll ausgebildete Strömung - Pr Prandtlzahl - Pr _t turbulente Prandtlzahl - Wärmestromdichte - Re _* Schubspannungsreynoldszahl - R _n Eigenfunktionen - Durchfluß-Reynoldszahl - Re v =D/ Rotations-Reynoldszahl - Ri Richardsonzahl - R Rohrradius - r Koordinate in radialer Richtung - dimensionslose Koordinate in radialer Richtung - T Temperatur - T Temperaturschwankung - T _b bulk temperature - mittlere Axialgeschwindigkeit - v Geschwindigkeit - v Geschwindigkeitsschwankung - turbulenter Wärmestrom - dimensionsloser Wandabstand - =1/6 Konstante - Integrationsvariable - Integrationsvariable - , ₁, ₂, dimensionslose Temperaturen - Wärmeleitzahl - _n Eigenwerte - kinematische Viskosität - Dichte - tangentiale Koordinate - , Hilfsfunktionen Indizes m in der Rohrmitte - r radial - w an der Rohrwand - z axial - 0 am Rohreintritt - 0 ohne Rotation - tangential     

Velocity and concentration measurements in a model diesel engine

Laser Doppler anemometry and Rayleigh scattering have been used to quantify the velocity and concentration fields after the start of injection in a model diesel engine motored at 200 rpm in the absence of compression. Fuel injection was simulated by a transient jet of vapour Freon-12 initiated at 40 degrees before top-dead-centre through a nozzle incorporated into the centre of a permanently open intake valve. Swirl was induced by means of 60 degree vanes located in the inlet, port. The piston configurations comprised a flat and a re-entrant piston-bowl.The results indicate that for the two nozzle geometries investigated the mass flux decays faster than momentum with nearly constant decay rates along the centreline. The nozzle with the larger exit diameter and wider jet angle gave rise to slower decay of both mass and momentum with associated lower velocity and concentration fluctuations.List of symbols D ₀ nozzle diameter - r radial coordinate - mean axial velocity - mean axial velocity at the centreline - ₀ mean axial velocity at the nozzle exit - rms of axial velocity fluctuations - mean concentration (mole fraction) - mean concentration at the nozzle exit - rms of concentration fluctuations - x axial coordinate A version of this paper was presented at the ASME Winter Annual Meeting of 1984 and printed in AMD, Vol. 66