Dynamic response of an axially loaded rigid sphere embedded in a saturated poroelastic medium

In this paper, an analytical solution for the response of a rigid sphere embedded in a full space poroelastic medium subjected to a dynamic lateral load is derived. The solution is obtained using Biots theory for acoustic waves. In this solution, the displacements of the solid skeleton and the pore pressure are expressed in terms of three scalar potentials. These potentials correspond to the wave velocities of the slow and fast compressional waves and to the shear wave. The governing equation for the dynamic motion is expressed in the frequency domain using Fourier transformation. Different boundary and load conditions were investigated. Curves showing variation in the fluid pressure and solid displacements with the loads frequency were plotted in non-dimensional forms.     

Shear wave dispersion and attenuation in periodic systems of alternating solid and viscous fluid layers

The attenuation and dispersion of elastic waves in fluid-saturated rocks due to the viscosity of the pore fluid is investigated using an idealized exactly solvable example of a system of alternating solid and viscous fluid layers. Waves in periodic layered systems at low frequencies are studied using an asymptotic analysis of Rytov’s exact dispersion equations. Since the wavelength of shear waves in fluids (viscous skin depth) is much smaller than the wavelength of shear or compressional waves in solids, the presence of viscous fluid layers necessitates the inclusion of higher terms in the long-wavelength asymptotic expansion. This expansion allows for the derivation of explicit analytical expressions for the attenuation and dispersion of shear waves, with the directions of propagation and of particle motion being in the bedding plane. The attenuation (dispersion) is controlled by the parameter which represents the ratio of Biot’s characteristic frequency to the viscoelastic characteristic frequency. If Biot’s characteristic frequency is small compared with the viscoelastic characteristic frequency, the solution is identical to that derived from an anisotropic version of the Frenkel–Biot theory of poroelasticity. In the opposite case when Biot’s characteristic frequency is greater than the viscoelastic characteristic frequency, the attenuation/dispersion is dominated by the classical viscoelastic absorption due to the shear stiffening effect of the viscous fluid layers. The product of these two characteristic frequencies is equal to the squared resonant frequency of the layered system, times a dimensionless proportionality constant of the order 1. This explains why the visco-elastic and poroelastic mechanisms are usually treated separately in the context of macroscopic (effective medium) theories, as these theories imply that frequency is small compared to the resonant (scattering) frequency of individual pores.     

Peristaltic motion of a third-order fluid in a planar channel

In order to determine the characteristics of the peristaltic transport of shear thinning non-Newtonian materials, the motion of a third-order fluid in a planar channel having walls that are transversely displaced by an infinite, harmonic traveling wave of large wavelength and negligibly small Reynolds number was analyzed using a perturbation expansion in terms of a variant of the Deborah number. Within the range of validity of this analysis, we found the pumping rate of a shear-thinning fluid is less than that for a Newtonian fluid having a shear viscosity the same as the lower-limiting viscosity of the nonNewtonian material. Also, the space of variables for which trapping of a bolus of fluid occurs is reduced for the shear-thinning fluid investigated here.     

Long wavelength approximation to peristaltic motion of a power law fluid

Peristaltic motion of a power law fluid in a two-dimensional channel is studied. Assuming that the wavelength of the peristaltic wave is large in comparison to the mean half-width of the channel, a solution for the stream function is obtained as an asymptotic expansion in terms of slope parameter. Expressions for axial pressure gradient and shear stress are derived. The effect of flow behaviour indexn on the streamline pattern and shear stress is studied and the phenomenon of trapping is discussed.     

Body waves in fractured porous media saturated by two immiscible newtonian fluids

A study of body waves in fractured porous media saturated by two fluids is presented. We show the existence of four compressional and one rotational waves. The first and third compressional waves are analogous to the fast and slow compressional waves in Biot's theory. The second compressional wave arises because of fractures, whereas the fourth compressional wave is associated with the pressure difference between the fluid phases in the porous blocks. The effects of fractures on the phase velocity and attenuation coefficient of body waves are numerically investigated for a fractured sandstone saturated by air and water phases. All compressional waves except the first compressional wave are diffusive-type waves, i.e., highly attenuated and do not exist at low frequencies.Now at Izmir Institute of Technology, Faculty of Engineering, Gaziosmanpasa Bulvari, No.16, Cankaya, Izmir, Turkey.     

Compressional waves in fluid-saturated porous solid containing a penny-shaped crack

Diffraction of normal compression waves by a penny-shaped crack in a fluid-saturated porous medium is investigated. Two wave types are considered, namely, compressional wave of the first kind, and the second kind. The former, also known as fast wave, propagates primarily through the solid, whereas the latter or slow wave, propagates mainly in the fluid. Each wave propagates in the medium along with induced wave of the same type in the companion constituent of the material. Application of Biot’s theory in conjunction with integral transform technique reduces the problem to a mixed boundary-value problem whose solution is in turn governed by a Fredholm integral equation of the second kind. Near-field and far-field solutions are obtained in terms of the dynamic stress-intensity factor and the scattering cross section, respectively. They are of particular importance to the linear elastic fracture mechanics (LEFM) and in the scattering theory of elastic waves. The mode I stress-intensity factors are computed numerically for a set of selected material property values, and shown graphically for various mass density and viscosity-to-permeability ratios. The obtained results reveal significant impact of the presence of pore fluid upon the stress-intensity factors, both magnitudes and frequencies at their peak values. The influence of the fluid is also observed from the calculated scattering cross sections of the scattered far-field. Accuracy of the present solution procedure is verified by comparing the numerical results with existing results in the limiting case of dry elastic materials.     

孔隙地层中声波诱导的电磁场测井响应模拟

声电效应测井过程中,当充流体井中发射的声波传播至井外时,在孔隙地层中引起含带电离 子的孔隙流体相对于固体骨架流动,导致电磁场. 在忽略转换电磁场对声场影响的条件下,先采用Biot-Rosenbaum模型计算孔隙地层中的 声场,再计算声波诱导的电场和磁场. 计算结果表明,有伴随纵波波群和横波波群的电磁场、有伴随斯通利波的电磁场,还有微弱的 临界折射电磁波. 计算还表明: 接收器偏离井轴越远,电场和磁场的强度越大; 在井壁处,磁场强度的波形与轴向电场强度的波形相似.     

Reflection and transmission of elastic waves by a moving slab

Scattering of elastic waves by a moving slab is considered. Two cases corresponding to strong contact and good lubrication are treated. It is shown that the motion introduces new effective compressional and shear wave velocities in the moving slab. The amplitude of the reflection and transmission coefficients are given for various angles of incidence, frequencies and velocities of motion.     

Numerical simulation of a fluid-saturated soil under a dynamically loaded pile

The behaviour of the soil under a dynamically loaded pile toe is studied. The soil is modelled as a fluid saturated porous continuum. The constitutive behaviour of the solid skeleton is described by the elasto-plastic model of Drücker-Prager. The wave propagation is simulated with a dynamical finite-element program.A two-phase model of soil gives useful information about effective stress and pore pressure in the soil. In saturated soil the main wave under the pile toe propagates more downards than in dry soil, due to the higher compressional stiffness in saturated soil. The plastic zone under the pile toe propagates with the velocity of the fast compressional wave. The pore fluid influences the plasticity strongly and can be expected to affect pile driving too.The distribution of effective stress and pore pressure under the pile toe depends on the permeability of the soil and cannot be calculated uniquely from a single-phase calculation. Therefore, a nonlinear soil cannot be modelled correctly as a conventional single-phase material.     

Dynamic shear and compressional behavior of polydimethylsiloxanes: Ultrasonic and low frequency characterization

The dynamic rheological properties of a series of entangled polydimethylsiloxanes are investigated. Rheometrical shear measurements are first reported for three different molecular weight samples, together with ultrasonic tests, thus allowing to obtain dynamic shear master curves over ten decades. Winter‘s model for slightly polydisperse polymers is used and works well in this case. This provides an example of the applicability of this model. The propagation of longitudinal waves in the MHz range for the same PDMSs is studied next. The results are analyzed and combined with static measurements of the compressibility. They indicate that different relaxation mechanisms have to be considered in shear and compression, and that the ratio between volume and shear viscosity is frequency dependent, with a preponderance for shear effects at higher frequencies. Shear and compressional relaxational mechanisms are also well separated. Moreover, it is verified that the shear and compressional moduli are independent of the molecular weight in the transition region, above a certain frequency. This frequency corresponds to a wavelength comparable to the distance between entanglements, in the case of shear waves. Received: 9 December 1997 Accepted: 30 March 1998     

Interaction of shear flows of an ideal incompressible fluid in a channel

The problem of the decay of an arbitrary discontinuity for the equations describing plane-parallel shear flows of an ideal fluid in a narrow channel is considered. The class of particular solutions corresponding to fluid flows with piecewise constant vorticity is studied. In this class, the existence of self-similar solutions describing all possible unsteady wave configurations resulting from the nonlinear interaction of the specified shear flows is established. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 6, pp. 34–47, November–December, 2006.     

Improvements to a smooth particle hydrodynamics simulator for investigating submarine landslide generated waves

Submarine landslides can exhibit complex rheologies, including a finite yield stress and shear thinning, yet are often simulated numerically using a Newtonian fluid rheology and simplistic boundary conditions. Here we present improvements made to a Smoothed Particle Hydrodynamics simulator to allow the accurate simulation of submarine landslide generated waves. The improvements include the addition of Bingham and Herschel-Bulkley rheologies, which better simulate the behavior of submarine mudflows. The interaction between the base of the slide and the slope is represented more accurately through the use of a viscous stress boundary condition. This condition treats the interface between the seafloor and the slide as a fluid boundary layer with a user-defined viscosity and length scale. Modifications to the pressure and density calculations are described that improve their stability for landslide generated wave scenarios. An option for pressure decomposition is introduced to prevent particle locking under high pressure. This facilitates the application of this simulator to landslide scenarios beneath significant water depths. Additional modifications to the reaveraging and renormalization routines improve the stability of the free surface and fluid density. We present the mathematical formulations of these improvements alongside commentary on their performance and applicability to landslide generated wave modeling. The modifications are verified against analytical fluid flow solutions and a wave generation experiment.     

Squirt-Flow in Fluid-Saturated Porous Media: Propagation of Rayleigh Waves

Propagation of attenuated waves is studied in a squirt-flow model of porous solid permeated by two different pore regimes saturated with same viscous fluid. Presence of soft compliant microcracks embedded in the grains of stiff porous rock defines the double-porosity formation. Microcracks and pores respond differently to the compressional effect of a propagating wave, which induces the squirt-flow from microcracks to pores. Elastodynamics of constituent particles in porous aggregate is represented through a single-porosity formulation, which involves the frequency-dependent complex moduli. This formulation is deduced as a special case of double-porosity formation allowing the wave-induced flow of pore-fluid. This squirt-flow model of porous solid supports the attenuated propagation of two compressional waves and one shear wave. Superposition of these body waves, subject to stress-free surface, defines the propagation of Rayleigh wave. This wave is governed by a complex irrational dispersion equation, which is solved numerically after rationalising into an algebraic equation. For existence of Rayleigh wave, a complex solution of the dispersion equation should represent a leaky wave, which decays for propagation along any direction in the semi-infinite medium. A numerical example is solved to analyse the effects of squirt-flow on phase velocity, attenuation and polarisation of the Rayleigh waves, for different combinations of parameters. Numerical results suggest the existence of an additional (second) Rayleigh wave in the squirt-flow model of dissipative porous solids.     

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法. 鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究. 先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解. 对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析. 分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成. 当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导. 当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导. 随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.      

Exact solutions to one-dimensional transient response of incompressible fluid-saturated single-layer porous media

Based on the Biot theory of porous media,the exact solutions to onedimensional transient response of incompressible saturated single-layer porous media under four types of boundary conditions are developed.In the procedure,a relation between the solid displacement u and the relative displacement w is derived,and the well-posed initial conditions and boundary conditions are proposed.The derivation of the solution for one type of boundary condition is then illustrated in detail.The exact solutions for the other three types of boundary conditions are given directly.The propagation of the compressional wave is investigated through numerical examples.It is verified that only one type of compressional wave exists in the incompressible saturated porous media.     

Simulation of peristaltic flow of chyme in small intestine for couple stress fluid

In this article couple stress fluid have been considered for the peristaltic flow of chyme in intestine. Problem under consideration have been formulated assuming that two non-periodic sinusoidal wave of different wavelength propagate with same speed c along the outer wall of the tube. Governing equations have been simplified under the assumptions of long wavelength and low Reynolds number approximation (such assumption are consistent that Re (Reynold number) is very small and long wavelength approximation also exists in the small intestine). Exact solutions have been evaluated for velocity and pressure rise. Physical behaviour of different parameter of couple stress fluid have been presented graphically for velocity, pressure rise, pressure gradient and frictional forces. The stream lines are also made against different parameters.     

Model of a strong discontinuity for the equations of spatial long waves propagating in a free-boundary shear flow

The long-wave equations describing three-dimensional shear wave motion of a free-surface ideal fluid are rearranged to a special form and used to describe discontinuous solutions. Relations at the discontinuity front are derived, and stability conditions for the discontinuity are formulated. The problem of determining the flow parameters behind the discontinuity front from known parameters before the front and specified velocity of motion of the front are investigated. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 4, pp. 206–213, July–August, 2008.     

Nonorthogonal stagnation-point flow over a lubricated surface

An attempt is made to study a steady two-dimensional flow of a viscous incompressible fluid incident at some angle onto a plate lubricated with a thin layer of a power-law fluid. Similar and nonsimilar solutions of the governing partial differential equations are obtained numerically by imposing the continuity of velocity and shear stress at the interface layer between the fluid and the lubricant. The Keller box method is applied to obtain the solutions. The limiting cases for full and no-slip conditions are compared.     

Wave forces on a submerged horizontal plate–Part II: Solitary and cnoidal waves

This paper is the companion to Part I under the same title, and is mainly concerned with wave loads due to nonlinear waves of solitary and cnoidal type propagating over a submerged, horizontal and thin plate. Following the development of the nonlinear model (via the Level I Green–Naghdi theory) for the flow of an incompressible and inviscid fluid given in Part I, the wave-induced loads on the submerged, fixed (and rigid) plate are calculated, and results are compared with the available laboratory data, and with linear solutions of the problem. Dependence of the loads on wave conditions (wave height and wavelength) and plate characteristics (submergence depth and plate width) are studied for both the solitary and cnoidal wave cases.