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 non- singular 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,

Heat and Mass Transfer in MHD Visco-Elastic Fluid Flow through a Porous Medium over a Stretching Sheet with Chemical Reaction

This paper presents the study of convective heat and mass transfer characteristics of an incompressible MHD visco-elastic fluid flow immersed in a porous medium over a stretching sheet with chemical reaction and thermal stratification effects. The resultant governing boundary layer equations are highly non-linear and coupled form of partial differential equations, and they have been solved by using fourth order Runge-Kutta integration scheme with Newton Raphson shooting method. Numerical computations are carried out for the non-dimensional physical parameters. Here a numerical has been carried out to study the effect of different physical parameters such as visco-elasticity, permeability of the porous medium, magnetic field, Grashof number, Schmidt number, heat source parameter and chemical reaction parameter on the flow, heat and mass transfer characteristics.

Visco-elastic fluid flow past an infinite vertical porous plate in the presence of first-order chemical reaction

An analysis has been developed to study the unsteady free convection flow of an incompressible visco-elastic fluid on a continuously moving vertical porous plate in the presence of a first-order chemical reaction. The governing equations are solved numerically using an implicit finite difference technique. The obtained numerical solutions are compared with the analytical solutions. The velocity profiles are presented. A parametric analysis is performed to illustrate the influences of the visco-elastic parameter, the dimensionless chemical reaction parameter, and the plate moving velocity on the steady state velocity profiles, the time dependent friction coefficient, the Nusselt number, and the Sherwood number.

Mixed Convection MHD Flow of Viscoelastic Fluid in a Porous Medium past a Hot Vertical Plate

The boundary layer flow of a steady incompressible and visco-elastic fluid with short memory (obeying Walters’ B fluid model) passing over a hot vertical porous plate has been investigated in the presence of transverse magnetic field. The momentum and energy equations are reduced to couple non-linear partial differential equations along with the boundary conditions by using a suitable similarity transformation. These partial differential equations are transformed to a system of coupled non-linear ordinary differential equations by employing a perturbation technique. The system is solved by developing a suitable numerical procedure such as implicit finite difference scheme along with Newton’s linearization method. The computational results for the flow quantities have presented graphically for the effects of thermal radiation, viscous dissipation, heat generation/absorption, visco-elasticity, Hartmann number and the permeability parameter. Results demonstrated that Prandtl number has more pronouncing effect on the temperature distribution rather than the viscosity parameter as well as the thermal radiation parameter. Further the velocity gradient changes significantly due to the presence of temperature dependent variable viscosity.

P-SV wave elastic impedance and fluid identification factor in weakly anisotropic media

The P-SV wave reflection coefficients in VTI and HTI media were obtained by approximation of the Jilek（2002a and b）equation in orthotropic anisotropic media.An approximate equation for P-SV wave elastic impedance can be derived from the combination of the new coefficients with S-wave elastic impedance（Duffaut et al.,2000）.On this basis, the fluid identification factor in weakly anisotropic media was constructed and used to identify the Castagna and Smith（1994）lithologic combination and achieved good results. Finally,we specifically analyzed the anisotropic parameter impacts P-SV wave elastic impedance and fluid factor trends.

Fluid sensitivity study of elastic parameters in lowmedium porosity and permeability reservoir rocks

In this article,based on the acoustic measurements of core samples obtained from the low to medium porosity and permeability reservoirs in the WXS Depression,the densities and P and S wave velocities of these core samples were obtained.Then based on these data,a series of elastic parameters were computed.From the basic theory and previous pore fluid research results,we derived a new fluid identification factor（F）.Using the relative variations,Ag/w and Ao/w,of the elastic parameters between gas and water saturated samples and between oil and water saturated samples,λρ,σHSFIF,Kρ,λρ-2μρ,and F as quantitative indicators,we evaluate the sensitivity of the different fluid identification factors to identify reservoir fluids and validate the effects by crossplots.These confirm that the new fluid identification factor（F） is more sensitive for distinguishing oil and water than the traditional method and is more favorable for fliud identification in low to medium porosity and permeability reservoirs.

Mechanical Analysis of Viscous-Elastic Fluid Acting on Residual Oil in the Micro Pore

In order to analyze the normal deviatoric stress that viscous-elastic fluid acting on the residual oil under the situation of different flooding conditions and different permeabilities, Viscous-elastic fluid flow equation is established in the micro pore by choosing the continuity equation, motion equation and the upper-convected Maxwell constitutive equation, the flow field is computed by using numerical analysis, the forces that driving fluid acting on the residual oil in micro pore are got, and the influence of flooding conditions, pore width and viscous-elasticity of driving fluid on force is compared and analyzed. The results show that: the more viscous-elasticity of driving fluid increases, the greater the normal deviatoric stress acting on the residual oil increases;using constant pressure gradient flooding, the lager the pore width is, the greater normal deviatoric stress acting on the residual oil will be.

Heat and Mass Transfer in Visco-Elastic Fluid through Rotating Porous Channel with Hall Effect

This paper examined the hydromagnetic boundary layer flow of viscoelastic fluid with heat and mass transfer in a vertical channel with rotation and Hall current. A constant suction and injection is applied to the plates. A strong magnetic field is applied in the direction normal to the plates. The entire system rotates with uniform angular velocity (Ω), about the axis perpendicular to the plates. The governing equations are solved by perturbation technique to obtain an analytical result for velocity, temperature, concentration distributions and results are presented graphically for various values of viscoelastic parameter (K2), Prandtl number (Pr), Schmidt number (Sc), radiation parameter (R), heat generation parameter (Qh) and Hall parameter (m).

Numerical Calculation of the Horizontal Stress Difference of Viscous-Elastic Fluid Acting on the Residual Oil in Micro Pore

In order to analyze the microscopic theory of viscous-elastic fluid flooding residual oil, the flow equation of polymer solution in the micro pore can be derived by selecting upper-convected Maxwell constitutive equation, continuity equation and motion equation. Then, the flow velocity field and stress field can be calculated under the boundary condition, and with the theory of stress tensor, the horizontal stress difference of polymer solution acting on the residual oil can be calculated. The results show that the greater the elasticity of viscous-elastic fluid is, the wider the flow channel is, the greater the horizontal stress difference is. The force acting on residual oil by viscous-elastic fluid can be increased by increasing the concentration of the polymer solution.

RESEARCH ON COMPLICATED DAM FOUNDATION DURING EXCAVATION WITH THREE DIMENSION VISCO ELASTIC PLASTIC CONDENSED FEM

On the basis of the two geological factors rock structure and ground stress environment,a visco elastic plastic model is established to analyze the rock stability of dam foundation and dam abutment during excavation by a kind of FE condensed method.Rock mechanics is applied in analysing the dynamic process of displacements,stresses,yielding destruction of construction base level,soft interface of dam foundation.Results of the FE analysis indicate that theories and methods in this paper are reasonable and reliable.

Hydroelastic Response of A Circular Plate in Waves on A Two-Layer Fluid of Finite Depth

The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interfacial wave modes, of a given wave frequency, on a two-layer fluid of finite and constant depth. In linear potential-flow theory, with the aid of angular eigenfunction expansions, the diffraction potentials can be expressed by the Bessel functions. A system of simultaneous equations is derived by matching the velocity and the pressure between the open-water and the platecovered regions, while incorporating the edge conditions of the plate. Then the complex nested series are simplified by utilizing the orthogonality of the vertical eigenfunctions in the open-water region. Numerical computations are presented to investigate the effects of different physical quantities, such as the thickness of the plate, Young’s modulus, the ratios of the densities and of the layer depths, on the dispersion relations of the flexural-gravity waves for the two-layer fluid. Rapid convergence of the method is observed, but is slower at higher wave frequency. At high frequency, it is found that there is some energy transferred from the interfacial mode to the surface mode.

ANALYSIS OF COUPLED-MODE FLUTTER OF PIPES CONVEYING FLUID ON THE ELASTIC FOUNDATION

The governing equation of solid-liquid couple vibration of pipe conveying fluid on the elastic foundation was derived. The critical velocity and complex frequency of pipe conveying fluid on Winkler elastic foundation and two-parameter foundation were calculated by po,ver series method. Compared,with pipe without considering elastic foundation, the numerical results show that elastic foundation can increase the critical flow velocity of static instability and dynamic instability of pipe. And the increase of foundation parameters may increase the critical flow velocity of static instability and dynamic instability of pipe, thereby delays the occurrence of divergence and flutter instability of pipe. For higher mass ratio beta, in the combination of certain foundation parameters, pipe behaves the phenomenon of restabilization and redivergence after the occurrence of static instability, and then coupled-mode flutter takes place.

An explicit finite element－finite difference method for analyzing the effect of visco－elastic local topography on the earthquake motion

An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite elementdiscrete model is formed by using the artificial boundary and finite element method, and the dynamic equationsof local nodes in the discrete model are obtained according to the theory of the special finite element method similar to the finite difference method, and then the explicit step-by-step integration formulas are presented by usingthe explicit difference method for solving the visco-elastic dynamic equation and Generalized Multi-transmittingBoundary. The method has the advantages of saving computing time and computer memory space, and it is suitable for any case of topography and has high computing accuracy and good computing stability.

Fluids discrimination by ray-path elastic impedance inversion: A successful case from Sulige tight gas field

Existing seismic prediction methods struggle to effectively discriminate between fluids in tight gas reservoirs,such as those in the Sulige gas field in the Ordos Basin,where porosity and permeability are extremely low and the relationship between gas and water is complicated.In this paper,we have proposed a comprehensive seismic fluid identification method that combines ray-path elastic impedance(REI)inversion with fluid substitution for tight reservoirs.This approach is grounded in geophysical theory,forward modeling,and real data applications.We used geophysics experiments in tight gas reservoirs to determine that Brie's model is better suited to calculate the elastic parameters of mixed fluids than the conventional Wood’s model.This yielded a more reasonable and accurate fluid substitution model for tight gas reservoirs.We developed a forward model and carried out inversion of REI.which reduced the non-uniqueness problem that has plagued elastic impedance inversion in the angle domain.Our well logging forward model in the ray-path domain with different fluid saturations based on a fluid substitution model proved that REI identifies fluids more accurately when the ray parameters are large.The distribution of gas saturation can be distinguished from the crossplot of REI(p=0.10)and porosity.The inverted ray-path elastic impedance profile was further used to predict the porosity and gas saturation profile.Our new method achieved good results in the application of 2D seismic data in the western Sulige gas field.

On Axially Symmetric Vibrations of Fluid Filled Poroelastic Spherical Shells

Employing Biot’s theory of wave propagation in liquid saturated porous media, waves propagating in a hollow poroelastic closed spherical shell filled with fluid are studied. The frequency equation of axially symmetric vibrations for a pervious and an impervious surface is obtained. Free vibrations of a closed spherical shell are studied as a particular case when the fluid is vanished. Frequency as a function of ratio of thickness to inner radius is computed in absence of dissipation for two types of poroelastic materials each for a pervious and an impervious surface. Results of previous works are obtained as a particular case of the present study.

Vibration of Visco-Elastic Parallelogram Plate with Parabolic Thickness Variation

The main objective of the present investigation is to study the vibration of visco-elastic parallelogram plate whose thickness varies parabolically. It is assumed that the plate is clamped on all the four edges and that the thickness varies parabolically in one direction i.e. along length of the plate. Rayleigh-Ritz technique has been used to determine the frequency equation. A two terms deflection function has been used as a solution. For visco-elastic, the basic elastic and viscous elements are combined. We have taken Kelvin model for visco-elasticity that is the combination of the elastic and viscous elements in parallel. Here the elastic element means the spring and the viscous element means the dashpot. The assumption of small deflection and linear visco-elastic properties of “Kelvin” type are taken. We have calculated time period and deflection at various points for different values of skew angles, aspect ratio and taper constant, for the first two modes of vibration. Results are supported by tables. Alloy “Duralumin” is considered for all the material constants used in numerical

ONE-DIMENSIONAL CONSOLIDATION OF LAYERED AND VISCO-ELASTIC SOILS UNDER ARBITRARY LOADING

Based on the layered visco-elastic soil model, according to the Terzaghi's one dimensional consolidation theory, by the method of Laplace transform and matrix transfer technique, the problems about the consolidation of layered and saturated visco-elastic soils under arbitrary loading were solved. Through deductions, the general solution, in the terms of layer thickness, the modulus and the coefficients of permeability and Laplacian transform's parameters was obtained. The strain and deformation of the layered and saturated visco-elastic soils under arbitrary loading can be calculated by Laplace inversion. According to the results of several numerical examples, the consolidation of visco-elastic soils logs behind that of elastic soils. The development of effective stress and the displacement is vibrant process under cyclic loading. Finally, an engineering case is studied and the results prove that the methods are very effective.

Effect of Non-Homogeneity on Thermally Induced Vibration of Orthotropic Visco-Elastic Rectangular Plate of Linearly Varying Thickness

The analysis presented here is to study the effect of non-homogeneity on thermally induced vibration of orthotropic visco-elastic rectangular plate of linearly varying thickness. Thermal vibrational behavior of non-homogeneous rectangular plates of variable thickness having clamped boundary conditions on all the four edges is studied. For non–homogeneity of the plate material, density is assumed to vary linearly in one direction. Using the method of separation of variables, the governing differential equation is solved. An approximate but quite convenient frequency equation is derived by using Rayleigh-Ritz technique with a two-term deflection function. Time period and deflection at different points for the first two modes of vibration are calculated for various values of temperature gradients, non- homogeneity constant, taper constant and aspect ratio. Comparison studies have been carried out with non-homogeneous visco-elastic rectangular plate to establish the accuracy and versatility.

Forward modeling of ocean-bottom cable data and wave-mode separation in fluid–solid elastic media with irregular seabed

In marine seismic exploration, ocean-bottom cable techniques accurately record the multicomponent seismic wavefield; however, the seismic wave propagation in fluid–solid media cannot be simulated by a single wave equation. In addition, when the seabed interface is irregular, traditional finite-difference schemes cannot simulate the seismic wave propagation across the irregular seabed interface. Therefore, an acoustic–elastic forward modeling and vector-based P-and S-wave separation method is proposed. In this method, we divide the fluid–solid elastic media with irregular interface into orthogonal grids and map the irregular interface in the Cartesian coordinates system into a horizontal interface in the curvilinear coordinates system of the computational domain using coordinates transformation. The acoustic and elastic wave equations in the curvilinear coordinates system are applied to the fluid and solid medium, respectively. At the irregular interface, the two equations are combined into an acoustic–elastic equation in the curvilinear coordinates system. We next introduce a full staggered-grid scheme to improve the stability of the numerical simulation. Thus, separate P-and S-wave equations in the curvilinear coordinates system are derived to realize the P-and S-wave separation method.

FORMATION OF SINGULARITY FOR COMPRESSIBLE VISCOELASTICITY

The formation of singularity and breakdown of classical solutions to the three- dimensional compressible viscoelasticity and inviscid elasticity are considered. For the compressible inviscid elastic fluids, the finite-time formation of singularity in classical solu- tions is proved for certain initial data. For the compressible viscoelastic fluids, a criterion in term of the temporal integral of the velocity gradient is obtained for the breakdown of smooth solutions.