Reflection and refraction of attenuated waves at boundary of elastic solid and porous solid saturated with two immiscible viscous fluids

The propagation of elastic waves is studied in a porous solid saturated with two immiscible viscous fluids.The propagation of three longitudinal waves is represented through three scalar potential functions.The lone transverse wave is presented by a vector potential function.The displacements of particles in different phases of the aggregate are defined in terms of these potential functions.It is shown that there exist three longitudinal waves and one transverse wave.The phenomena of reflection and refraction due to longitudinal and transverse waves at a plane interface between an elastic solid half-space and a porous solid half-space saturated with two immiscible viscous fluids are investigated.For the presence of viscosity in pore-fluids,the waves refracted to the porous medium attenuate in the direction normal to the interface.The ratios of the amplitudes of the reflected and refracted waves to that of the incident wave are calculated as a nonsingular system of linear algebraic equations.These amplitude ratios are used to further calculate the shares of different scattered waves in the energy of the incident wave.The modulus of the amplitude and the energy ratios with the angle of incidence are computed for a particular numerical model.The conservation of the energy across the interface is verified.The effects of variations in non-wet saturation of pores and frequencies on the energy partition are depicted graphically and discussed.     

Boundary layer flow of a Jeffrey fluid with convective boundary conditions

The boundary layer stretched flow of a Jeffrey fluid subject to the convective boundary conditions was investigated. The governing dimensionless problems were computed by using the homotopy analysis approach. Convergence of the derived solutions was checked and the influence of embedded parameters was analyzed by plotting graphs. It was noticed that the velocity increases with an increase in the Deborah number. Furthermore, it was found that the temperature is also an increasing function of the Biot number. We further found that for fixed values of other parameters, the local Nusselt number increases by increasing the suction parameter and Deborah number. Numerical values of the skin friction coefficient and local Nusselt numbers were computed and examined. Copyright © 2011 John Wiley & Sons, Ltd.     

An extended theory to predict the onset of viscous instabilities for miscible displacements in porous media

Many enhanced oil recovery schemes involve the displacement of oil by a miscible fluid. Whether a displacement is stable or unstable has a profound effect on how efficiently a solvent displaces oil within a reservoir. That is, if viscous fingers are present, the displacement efficiency and, hence, the economic return of the recovery scheme is seriously impaired bacause of macroscopic bypassing of the oil. As a consequence, it is of interest to be able to predict the boundary which separates stable displacements from those which are unstable.This paper presents a dimensionless scaling group for predicting the onset of hydrodynamic instability of a miscible displacement in porous media. An existing linear perturbation analysis was extended in order to obtain the scaling group. The new scaling group differs from those obtained in previous studies because it takes into account a variable unperturbed concentration profile, both transverse dimensions of the porous medium, and both the longitudinal and the transverse dispersion coefficient.It has been shown that stability criteria derived in the literature are special cases of the general condition given here. Therefore, the stability criterion obtained in this study should be used for a displacement conducted under arbitrary conditions. The stability criterion is verified by comparing it with miscible displacement experiments carried out in a Hele-Shaw cell. Moreover, a comparison of the theory with some porous medium experiments from the literature also supports the validity of the theory.Nomenclature c solvent concentration - C _g fractional glycerine volume - D molecular diffusion coefficient, cm²/s - D _L longitudinal dispersion coefficient, cm²/s - D _T transverse dispersion coefficient, cm²/s - g gravitational acceleration, cm/s² - h distance between the plates, cm - I _sr dimensionless scaling group - k permeability, cm² - L _x width of the porous medium, cm - L _y height of the porous medium, cm - t time, s - u velocity in thex direction, cm/s - v velocity in they direction, cm/s - V displacement velocity, cm/s - w velocity in thez direction, cm/s - z length of the graded viscosity bank, cm - eigenvalue in thex direction - eigenvalue in they direction - wave number - viscosity, poise - density, g/cc - time constant, s^-1 - porosity     

Plane waves in a semi-infinite fluid saturated porous medium

The field equations governing the propagation of waves in an incompressible liquid-saturated porous medium are investigated and a general solution is presented. It has been revealed that coupled longitudinal and transverse waves propagate in the porous medium. The propagation of transverse waves in the fluid phase is completely due to the interaction between the solid and fluid phases. The dispersion relationship and attenuation features are discussed. Unlike other investigations, all explicit forms of the arguments are derived. The reflection of the plane harmonic waves at the plane, traction-free boundary, which shows the influence of the dissipation on the velocity, and the attenuation coefficients of the reflected waves is studied. It is of interest that pore pressure is produced in the process of reflection, even in the case of the incidence of transverse waves.     

Propagation of acceleration waves in incompressible saturated porous solids

Within the framework of the incompressible porous media model, the propagation properties of acceleration waves in liquid-filled porous solids is discussed. The incompressibility of the two constituents in the model forces the amplitudes of the longitudinal waves in the skeleton and in the liquid to satisfy a certain relation. The two propagation speeds are presented by examination for the existence of acceleration waves and only longitudinal and transverse waves are realizable in the incompressible two-phase porous materials.     

Spin-up in a saturated porous medium

It is shown that, on the Brinkman model, spin-up is confined to boundary layers whose thickness is of order k ^1/2, and the spin-up is established in a time of order k/, where k, , and denote permeability, density, porosity and dynamic viscosity, respectively.     

Boundary conditions for stokes flows near a porous membrane

A theoretical development is carried out to model the boundary conditions for Stokes flows near a porous membrane, which, in general, allows non-zero slip as well as normal flow at the surface. Two types of models are treated: an infinitesimally thin plate with a periodic array of circular apertures and a series of parallel slits. For Stokes flows, the mean normal flux and slip velocity are proportional to the pressure difference across the membrane and the average shear stress at the membrane, respectively. The appropriate proportionality constants which depend on the membrane geometry are calculated as functions of the porosity. An interesting feature of the results is that the slip at the membrane has, in general, a direction different from that of the applied shear for these models.     

Hydrodynamics of viscous fluid through porous slit with linear absorption

The exact solutions for the viscous fluid through a porous slit with linear ab-sorption are obtained. The Stokes equation with non-homogeneous boundary conditions is solved to get the expressions for the velocity components, pressure distribution, wall shear stress, fractional absorption, and leakage flux. The volume flow rate and mean flow rate are found to be useful in obtaining a convenient form of the longitudinal velocity component and pressure difference. The points of the maximum velocity components for a fixed axial distance are identified. The value of the linear absorption parameter is ran-domly chosen, and the rest available data of the rat kidney to the tabulate pressure drop and fractional absorption are incorporated. The effects of the linear absorption, uniform absorption, and flow rate parameters on the flow properties are discussed by graphs. It is found that forward flow occurs only if the volume flux per unit width is greater than the absorption velocity throughout the length of the slit, otherwise back flow may occur. The leakage flux increases with the increase in the linear absorption parameter. Streamlines are drawn to help the analysis of the flow behaviors during the absorption of the fluid flow through the renal tubule and purification of blood through an artificial kidney.     

流体饱和多孔介质黏弹性动力人工边界

基于Biot流体饱和多孔介质本构方程,分别考察具有辐射阻尼性质的外行柱面波和球 面波在圆柱面和球面人工边界上引起的法向、切向应力的表达式. 在应力表达形式上,固相 介质和孔隙流体的法向和切向应力都是由两项组成,它们分别与质点的位移和速度成正比, 因此,可在人工边界的法向和切向设置连续分布的并联弹簧------黏滞阻尼器,用来模拟人工边 界以外的无限域介质对来自有限计算域的外行波动的能量吸收作用,从而形成了流体饱和多 孔介质的黏弹性动力人工边界. 流体饱和多孔介质的黏弹性动力人工边界可方便地与大型通 用软件结合,用于分析饱和土中复杂的结构-地基动力相互作用问题. 算例表明流体饱和多 孔介质黏弹性动力人工边界具有较好的精度和稳定性.     

Rayleigh waves with impedance boundary conditions in anisotropic solids

The paper is concerned with the propagation of Rayleigh waves in an elastic half-space with impedance boundary conditions. The half-space is assumed to be orthotropic and monoclinic with the symmetry plane _x3=0$x_3=0$. The main aim of the paper is to derive explicit secular equations of the wave. For the orthotropic case, the secular equation is obtained by employing the traditional approach. It is an irrational equation. From this equation, a new version of the secular equation for isotropic materials is derived. For the monoclinic case, the method of polarization vector is used for deriving the secular equation and it is an algebraic equation of eighth-order. When the impedance parameters vanish, this equation coincides with the secular equation of Rayleigh waves with traction-free boundary conditions.     

Building generalized open boundary conditions for fluid dynamics problems

This paper deals with the design of an efficient open boundary condition (OBC) for fluid dynamics problems. Such problematics arise, for instance, when one solves a local model on a fine grid that is nested in a coarser one of greater extent. Usually, the local solution U^loc is computed from the coarse solution U^ext, thanks to an OBC formulated as , where B_h and B_H are discretizations of the same differential operator (B_h being defined on the fine grid and B_H on the coarse grid). In this paper, we show that such an OBC cannot lead to the exact solution, and we propose a generalized formulation , where g is a correction term. When B_h and B_H are discretizations of a transparent operator, g can be computed analytically, at least for simple equations. Otherwise, we propose to approximate g by a Richardson extrapolation procedure. Numerical test cases on a 1D Laplace equation and on a 1D shallow water system illustrate the improved efficiency of such a generalized OBC compared with usual ones. Copyright © 2012 John Wiley & Sons, Ltd.     

Horizontal boundary-layer natural convection in a porous medium saturated with a gas

Boundary-layer analysis is performed for free convection flow over a hot horizontal surface embedded in a porous medium saturated with a gas of variable properties. The variable gas properties are accounted for via the assumption that thermal conductivity and dynamic viscosity are proportional to temperature. A similarity solution is shown to exist for the case of constant surface temperature. Numerical results for the stream function, horizontal velocity, and temperature profiles within the boundary layer as well as for the mass of entrained gas, surface slip velocity, and heat transfer rate at different values of the wall-temperature parameter are presented. Asymptotic solutions for large heating are also available to support the numerical work.     

Boundary conditions control for a shallow‐water model

A variational data assimilation technique was used to estimate optimal discretization of interpolation operators and derivatives in the nodes adjacent to the rigid boundary. Assimilation of artificially generated observational data in the shallow‐water model in a square box and assimilation of real observations in the model of the Black sea are discussed. It is shown in both experiments that controlling the discretization of operators near a rigid boundary can bring the model solution closer to observations as in the assimilation window and beyond the window. This type of control also allows to improve climatic variability of the model. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling high speed flow of a compressible fluid in a saturated porous medium

A modification to the Forchheimer-Brinkman equation, for the modelling of high speed flow of a compressible fluid in a dense saturated porous medium, is proposed. The modified equation is applied to a flow in which choking can occur.     

A model for solute transport following an abrupt pressure impact in saturated porous media

A mathematical model is developed of an abrupt pressure impact applied to a compressible fluid with solute, flowing through saturated porous media. Nondimensional forms of the macroscopic balance equations of the solute mass and of the fluid mass and momentum lead to dominant forms of these equations. Following the onset of the pressure change, we focus on a sequence of the first two time intervals at which we obtain reduced forms of the balance equations. At the very first time period, pressure is proven to be distributed uniformly within the affected domain, while solute remains unaffected. During the second time period, the momentum balance equation for the fluid conforms to a wave form, while the solute mass balance equation conforms to an equation of advective transport. Fluid's nonlinear wave equation together with its mass balance equation, are separately solved for pressure and velocity. These are then used for the solution of solute's advective transport equation. The 1-D case, conforms to a pressure wave equation, for the solution of fluid's pressure and velocity. A 1-D analytical solution of the transport problem, associates these pressure and velocity with an exponential power which governs solute's motion along its path line.     

Boundary condition formulation for viscous fluid flow problems

In the solution of the Navier-Stokes equations by difference methods in infinite regions, the question arises as to the nature of the approximate boundary conditions at those portions of the computational region boundary where these conditions are not determined directly by the formulation of the basic problem. In certain cases of practical importance, these boundary conditions may be obtained by coupling the N-S equations with equations which are similar to the boundary-layer equations.In the present paper, we propose boundary conditions for the case of viscous incompressible fluid flow. Their application is illustrated for the problem of flow past the leading edge of a semi-infinite flat plate.The author wishes to thank I. Yu. Brailovskaya and L. A. Chudov for helpful suggestions in the course of this investigation.     

Advective fluxes in multiphase porous media under nonisothermal conditions

We consider the case in which more than one fluid phase occupies the void space of a porous medium. The advective flux law is formulated for a fluid phase, under nonisothermal conditions and with the presence of solutes in the fluid phases. The derivation of the flux laws is based on an approximated version of the averaged balance equation for linear momentum. Taking into account momentum transfer through the interface between the fluid phases, leads to coupling between the flow in adjacent phases. Fluxes are also shown to depend on the surface tension at the interface between the adjacent fluid phases. Since the latter depends on temperature and solute concentration in the two phases, the advective flux is shown to depend on both temperature and solute concentration gradients in the two phases. A preliminary order of magnitude analysis gives conditions under which the coupling phenomena are not negligible. The approach is applied to the unsaturated zone, as a typical example of a multiphase porous medium.The main conclusion is that the well known Darcy law for single phase flow, may have to be modified for a multi fluid phase system, especially when temperature and solute concentration are not uniform.     

Natural convection flow of a viscous incompressible fluid in a rectangular porous cavity heated from below with cold sidewalls

We consider the flow, which is induced by differential heating on the boundaries of a porous cavity heated from below. In particular we allow the sidewalls to have the same cold temperature as the upper surface, and thus the problem is a variant of the Darcy-Bénard convection problem, but one where there is flow at all non-zero Grashof numbers. Attention is focused on how the flow and heat transfer is affected by variations in the cavity aspect ratio, the Grashof number and the Darcy number. The flow becomes weaker as the Darcy number decreases from the pure fluid limit towards the Darcy-flow limit. In addition the number of cells which form in the cavity varies primarily with the aspect ratio and is always even due to the symmetry imposed by the cold sidewalls.     

Comparison of outflow boundary conditions for subsonic aeroacoustic simulations

Aeroacoustics simulations require much more precise boundary conditions than classical aerodynamics. Two classes of non‐reflecting boundary conditions for aeroacoustics are compared in the present work: the characteristic analysis‐based methods and the Tam and Dong approach. In the characteristic methods, waves are identified and manipulated at the boundaries, whereas the Tam and Dong approach use modified linearized Euler equations in a buffer zone near outlets to mimic a non‐reflecting boundary. The principles of both approaches are recalled, and recent characteristic methods incorporating the treatment of transverse terms are discussed. Three characteristic techniques—the original Navier–Stokes characteristic boundary conditions (NSCBC) of Poinsot and Lele and two versions of the modified method of Yoo and Im—are compared with the Tam and Dong method for four typical aeroacoustics problems: vortex convection on a uniform flow, vortex convection on a shear flow, acoustic propagation from a monopole, and acoustic propagation from a dipole. Results demonstrate that the Tam and Dong method generally provides the best results and is a serious alternative solution to characteristic methods even though its implementation might require more care than the usual NSCBC approaches. Copyright © 2011 John Wiley & Sons, Ltd.