Measurement of the relaxation time in Darcy flow

The inertia of a liquid flowing through a porous medium is normally ignored, but if the acceleration is great, it may be important. The relaxation time, defined so that it alone accounts for the inertia, has been determined experimentally with a simple oscillator. A U-Tube is provided with a porous plug and filled with a liquid. During pendulation of the liquid, the frequency and the damping define the relaxation time. The measured value of the relaxation time is about 10 times the theoretical estimate derived from Navier-Stokes equation.Symbols E modulus of elasticity - E _D dissipated energy - E _k kinetic energy - g acceleration of gravity - G pressure gradient - h height - K ₀ permeability - L length of porous plug - n porosity - P dissipated power - pressure - R half the tube length - R _c radius of the tube bend - r radial coordinate - r _o radius of the tube - s coordinate along a streamline in the tube - t time - v flux per unit area - it relaxation time - , auxiliary variables - , v dynamic and kinematic viscosity - , velocity potential for inviscid flow and gravity potential - dissipation function - displacement of the liquid - , _o frequency of damped and undamped oscillations     

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

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

Transient oscillatory and constantly accelerated non-Newtonian flow in a porous medium

This paper looks at the magnetohydrodynamic (MHD) analysis for transient flow of an Oldroyd-B fluid in a porous medium. The presented analysis takes into account the modified Darcy's law. The flow is induced due to constantly accelerated and oscillating plate. Expressions for the corresponding velocity field and the adequate tangential stress are determined by means of the Fourier sine transform. The influence of various parameters of interest on the velocity and tangential stress has been shown and discussed. A comparison for different kinds of fluids is also provided.     

Unsteady flow of ground-water in an aquifer system with viscoelastic properties

The transient flow behaviour of groundwater in an aquifer-aquitard system with viscoelastic properties is studied. On the basis of previous work (Hantush, Neuman, Brutsaert, Corapcioglu), new partial-differential integral equations are assembled. The well-known equations (Hantush, Brutsaert) are special cases of these equations, which describe the flow of slightly compressible groundwater in aquifer-aquitard layers with viscoelastic properties.Analytical solutions of the coupled partial differential-integral equations are obtained by Laplace transform. The viscoelastic properties enhance the heterogeneities of an elastic aquifer system with delay and feed. Numerical inversion of the Laplace transformed drawdown is in good agreement with the analytical solutions obtained.     

The inertial effect on the natural convection flow within a fluid-saturated porous medium

The general momentum equation for fluid flow within a porous medium is supposedly valid for any fluid-porous medium configuration. One of the main concerns of using the general equations refers to the inclusion of both inertia terms, namely, the convective inertia term and the Forchheimer term. In this study, we go beyond the important discussion about the correctness of including both terms in the general momentum equations by focusing upon the effect of the convective inertia term on the heat transfer results. The fluid-porous medium system considered here is a cavity bounded by solid surfaces with vertical walls maintained at constant but different temperatures. The natural convection problem is solved numerically, and the results are compared with a general theory developed by using the method of scale analysis. It is demonstrated that the convective inertia term effect upon the heat transfer results is minor for 0.01 ≤ Pr ≤ 1, 10 ≤ Ra_D ≤ 10⁴, 10⁻⁸ ≤ Da ≤ 10⁻², and porosities 0.4 and 0.8. It is also shown that, contrary to the general belief, the convective inertial effect upon the heat transfer within the cavity is minimized when the Prandtl number is reduced.     

Stokes' first problem for the fourth order fluid in a porous half

In this study, the flow of a fourth order fluid in a porous half space is modeled. By using the modified Darcy's law, the flow over a suddenly moving flat plate is studied numerically. The influence of various parameters of interest on the velocity profile is revealed.     

Solution for the oscillation of a newtonian fluid in a co-axial viscometer

When one cup of a co-axial viscometer oscillates,the measured moment on theother(stationary)cup shows a phase lag,partly due to the inertial effect of the fluidwithin the gap between the cups.In this paper such an effect is illustrated by a newexact solution of the Navier-Stokes equation,which is derived herein by a scheme ofreducing it to a two-point boundary value problem for ODEs,Our numerical resultsindicate that,as the Womersley number a or the dimensionless gap widthδincreases,the fluid velocity profile within the gap gradually deviates from the linear one andtransits to that of the boundary layer type,with the result that the moment decreasesin the magnitude and lags behind in the phase.With the advantage of high accuracyand excellent stability,the scheme proposed herein can readily be extended to solveother linear periodic problems.     

Modelling of flow pattern in a shallow aquifer affected by reservoirs

A compartmental model is developed to estimate flow parameters of a shallow aquifer affected by water loads in surface reservoirs and to evaluate its nonsteady flow distribution. The method incorporates temporal piezometric head measurements and sampling of water for dissolved chemicals and isotopes analyses. Each compartment is governed by a set of equations describing the conservation of linear momentum and mass balance expressions for water, isotopes, and dissolved chemicals. The number of compartmental balance expressions always must be greater than that of the unknown flow parameters associated with each compartment. An optimization method is described to yield spatial distribution of aquifer storativity, transmissivity, porosity, leakage, and compliance coefficients and fluxes leaking into an aquifer's lower boundary. Future predictions of an aquifer's piezometric head distribution in a compartmental system is formulated on the basis of the estimated flow parameters and the leakage components.Compartmental modelling which incorporates concentrations of environmental tracers, may yield efficiency in computing resources and accuracy enhancement for predicting an aquifer's flow regime.     

Oblique stagnation-point flow of an incompressible visco-elastic fluid towards a stretching surface

An analysis is made of the steady two-dimensional stagnation-point flow of an incompressible viscoelastic fluid over a flat deformable surface when the surface is stretched in its own plane with a velocity cx, where x is the distance from the stagnation-point and c is a positive constant. It is shown that for a viscoelastic fluid of short memory (obeying Walters’ B model), a boundary layer is formed when the stretching velocity of the surface is less than ax, where ax+2by is the inviscid free-stream velocity and y is the distance normal to the plate, a and b being constants and the velocity at a point increases with increase in the elasticity of the fluid. On the other hand an inverted boundary layer is formed when the surface stretching velocity exceeds ax and the velocity decreases with increase in the elasticity of the fluid. A novel result of the analysis is that the flow near the stretching surface is that corresponding to an inviscid stagnation-point flow when a=c. Temperature distribution in the boundary layer is found in three cases, namely: (i) the sheet with constant surface temperature (CST); (ii) the sheet with variable surface temperature (VST) and (iii) the sheet with prescribed quadratic power law surface heat flux (PHF) for various values of non-dimensional parameters. It is found that in all the three cases when a/c>1, temperature at a point decreases with increase in the elasticity of the fluid and when a/c<1, temperature at a point increases with increase in the elasticity of the fluid. Further temperature at a point decreases with increase in the radiation parameter and wall temperature parameter.     

Effective Radius and Underpressure of an Air Extraction Well in a Heterogeneous Porous Medium

We consider the compressible steady state air flow in a porous medium caused by an extraction well and governed by Darcy's law. For a homogeneous soil matrix we have derived formulas (in 2-D and 3-D) to determine the effective radius of a single well depending on the well position and the depth of the domain. For inhomogeneous case (in 2-D) the influence of the heterogeneity, well position and the depth of the water table on the effective radius and on the pressure at the well is studied.     

Heat transfer in oblique stagnation-point flow of an incompressible viscous fluid towards a stretching surface

Steady two-dimensional oblique stagnation-point flow of an incompressible viscous fluid over a flat deformable sheet is investigated when the sheet is stretched in its own plane with a velocity proportional to the distance from the stagnation-point. It is shown that the flow has a boundary layer structure for values of a/c (> 1), where ax+2by and cx are the x-component of the free stream velocity and the stretching velocity of the plate respectively, x being the distance from the stagnation-point. On the other hand when a/c < 1, the flow has an inverted boundary layer structure. It is also observed that the velocity at a point increases with increase in the free stream shear. For a fixed value of a/c, the streamlines becomes more and more oblique towards the left of the stagnation-point with increase in b/c where b > 0. On the other hand the streamlines become increasingly oblique to the right of the stagnation-point with increase in |b/c| when b < 0. For a fixed value of the Prandtl number Pr, temperature at a point decreases with increase in a/c. Further for a given value of a/c, the surface heat flux increases with increase in Pr.     

Transient mixed convection flow of a second grade viscoelastic fluid past an inclined backward facing step

The transient mixed convection of a second-grade viscoelastic fluid past an inclined backward facing step was studied numerically. The combined effects of the Reynolds number, the elastic effect, the inclined angle of the flow channel on the reattachment length, and the phenomena of heat transfer are examined during the development of the flow field. The Gauss-Seidal method with successive over relaxation was implemented to solve the stream-vorticity and energy equations. The results indicate that the reattachment length increases to the maximum as the inclined angle increases up to 150° or 180°. At these cases, the point of reattachment is close to the point of the local maximum value of Nu_x or is overshooting it. It is observed that the reattachment length increases as the Reynolds number increases or the elastic coefficient decreases. In the meantime, the contact point of the isotherm on the upper plate moves upward and is close to upstream flow as the inclined angle is around 150°.     

The numerical solution of the unsteady natural convection flow in a square cavity at high Rayleigh number using SADI method

The unsteady natural convection flow in a square cavity at high Rayleigh number Ra=10 ⁷ and 2×10 ⁷ has been computed using cubic spline integration. The required solutions to the two dimensional Navier-Stokes and energy equations have been obtained using two alternate numerical formulations on non-uniform grids. The main features of the transient flow have been briefly discussed. The results obtained by using the present method are in good agreement with the theoretical predictions ^[1,2].The steady state results have been compared with accurate solutions presented recently for Ra=10 ⁷.     

Transient free‐surface flow of a viscoelastic fluid in a narrow channel

The interplay between inertia and elasticity is examined for transient free‐surface flow inside a narrow channel. The lubrication theory is extended for the flow of viscoelastic fluids of the Oldroyd‐B type (consisting of a Newtonian solvent and a polymeric solute). While the general formulation accounts for non‐linearities stemming from inertia effects in the momentum conservation equation, and the upper‐convected terms in the constitutive equation, only the front movement contributes to non‐linear coupling for a flow inside a straight channel. In this case, it is possible to implement a spectral representation in the depthwise direction for the velocity and stress. The evolution of the flow field is obtained locally, but the front movement is captured only in the mean sense. The influence of inertia, elasticity and viscosity ratio is examined for pressure‐induced flow. The front appears to progress monotonically with time. However, the velocity and stress exhibit typically a strong overshoot upon inception, accompanied by a plug‐flow behaviour in the channel core. The flow intensity eventually diminishes with time, tending asymptotically to Poiseuille conditions. For highly elastic liquids the front movement becomes oscillatory, experiencing strong deceleration periodically. A multiple‐scale solution is obtained for fluids with no inertia and small elasticity. Comparison with the exact (numerical) solution indicates a wide range of validity for the analytical result. Copyright © 2004 John Wiley & Sons, Ltd.     

Interaction of differently shaped bodies in a potential flow of perfect compressible fluid: Axisymmetric internal problem

An approach is proposed to investigate an axisymmetric system consisting of an infinite thin elastic cylindrical shell that contains a potential flow of perfect compressible fluid and a periodically vibrating spherical inclusion. The approach emerged as part of a project devoted to developing methods to bring plugged oil wells back into production by the Vibration Theory Department of the S. P. Timoshenko Institute of Mechanics. This mathematical approach allows transforming the general solutions to equations of mathematical physics from one coordinate system to another to obtain an exact analytical solution (in the form of Fourier series) to interaction problems for systems of rigid and elastic bodies __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 9, pp. 16–31, September 2006.     

Analytical quantification of coefficients in the Ergun equation for fluid friction in a packed bed

A new analytical derivation for momentum transport during laminar flow through granular porous media is discussed and some of its implied results described. In the very low Reynolds number regime fully developed laminar flow is assumed and in the higher laminar Reynolds number regime the Forchheimer (non-Darcy) effect is modelled through reference to form drag induced by the solid constituents of the porous medium. The results are compared to the Ergun equation, which is empirically based on experimental measurements, and the correspondence is shown to be remarkably close.     

Propogation of waves in a linear elastic fluid

The non-classical symmetry method is used to determine particular forms of the arbitrary velocity and forcing terms in a linear wave equation used to model the propogation of waves in a linear elastic fluid. The behaviour of solutions derived using the non-classical symmetry method are discussed. Solutions satisfy a given initial profile and wave velocity. For some solutions the arbitrary forcing terms and wave velocity can be written in terms of the initial wave profile. Relationships between the arbitrary forcing, arbitrary velocity and the solution are derived.     

The in-plane stiffnesses of a honeycomb core including the thickness effect

Summary The subject of the consideration is the contribution of a regular honeycomb core to the effective in-plane stiffnesses of a sandwich structure. Due to the coupling of the core displacements with those of the sandwich face sheets, the stiffness contribution of the core is not proportional to its total thickness, as could be expected. The corresponding thickness effect is investigated by means of an appropriate closed-form approach. In doing so, the total elastic core strain energy is calculated based on an adequately chosen displacement representation. Further on, the resultant effective stiffnesses are derived as a function of the total core thickness. A comparative computation of the effective stiffnesses by finite element analysis gives good agreement. Received 5 August 1997; accepted for publication 23 October 1997     

Unsteady flow in the Ekman layer of an elastico-viscous liquid

The transient flow in the Ekman layer of an elastico-viscous liquid near a flat plate is discussed. Initially the fluid and the plate were rotating together and the plate then suddently starts moving with a uniform velocity in its own plane relative to the rotating frame of reference. It is shown that the ultimate steady state is reached through decay of inertial oscillations whose frequency decreases with increase in the elastic parameter.