Ultrasonic velocities of a set of saturated sandstone samples were measured at simulated in-situ pressures in the laboratory.The samples were obtained from the W formation of the WXS Depression and covered low to near...Ultrasonic velocities of a set of saturated sandstone samples were measured at simulated in-situ pressures in the laboratory.The samples were obtained from the W formation of the WXS Depression and covered low to nearly high porosity and permeability ranges.The brine and four different density oils were used as pore fluids,which provided a good chance to investigate fluid viscosity-induced velocity dispersion.The analysis of experimental observations of velocity dispersion indicates that（1）the Biot model can explain most of the small discrepancy（about 2–3%）between ultrasonic measurements and zero frequency Gassmann predictions for high porosity and permeability samples saturated by all the fluids used in this experiment and is also valid for medium porosity and permeability samples saturated with low viscosity fluids（less than approximately 3 mP·S）and（2）the squirt flow mechanism dominates the low to medium porosity and permeability samples when fluid viscosity increases and produces large velocity dispersions as high as about 8%. The microfracture aspect ratios were also estimated for the reservoir sandstones and applied to calculate the characteristic frequency of the squirt flow model,above which the Gassmann’ s assumptions are violated and the measured high frequency velocities cannot be directly used for Gassmann’s fluid replacement at the exploration seismic frequency band for W formation sandstones.展开更多
Mesoscopic fluid flow is the major cause of wave attenuation and velocity dispersion at seismic frequencies in porous rocks. The Johnson model provides solutions for the frequency-dependent quality factor and phase ve...Mesoscopic fluid flow is the major cause of wave attenuation and velocity dispersion at seismic frequencies in porous rocks. The Johnson model provides solutions for the frequency-dependent quality factor and phase velocity in partially saturated porous media with pore patches of arbitrary shapes. We use the Johnson model to derive approximations for the quality factor Q at the high and low frequency limit, and obtain the approximate equation for Qmin based on geophysical and geometric parameters. A more accurate equation for Qmin is obtained after correcting for the linear errors between the exact and approximate Q values. The complexity of the pore patch shape affects the maximum attenuation of Qmin and the transition frequency ftr; furthermore, the effect on ftr is stronger than that on Qmin. Numerical solutions to Biot’s equation are computationally intensive; thus, we build an equivalent viscoelastic model on the basis of the Zener model, which well approximates the wave attenuation and dispersion in porous rocks in the seismic band.展开更多
There are complex heterogeneous entities in the underground medium,and the heterogeneous scale has a substantial impact on wave propagation.In this study,we used a set of 11 samples of glass beads as high-velocity het...There are complex heterogeneous entities in the underground medium,and the heterogeneous scale has a substantial impact on wave propagation.In this study,we used a set of 11 samples of glass beads as high-velocity heterogeneous bodies to evaluate the impact of such heterogeneous bodies on the propagation of P-wave.We vary the heterogeneous scale by changing the diameter of the glass beads from 0.18 to 11 mm while keeping the same volume proportion(10%)of the beads for the set of 11 samples.The pulse transmission method was used to record measurements at the ultrasonic frequencies of 0.34,0.61,and 0.84 MHz in the homogeneous matrix.The relationship between P-wave fi eld features and heterogeneity scale,P-wave velocity,and the multiple of the wave number and heterogeneous scale(ka)was observed in the laboratory,which has sparked widespread interest and research.Heterogeneous scale affects P-wave propagation,and its wave field changes are complex.The waveform,amplitude,and velocity of the recorded P-waves correlate with the heterogeneous scale.For the forward scattering while large-scale heterogeneities,noticeable direct and diff racted waves are observed in the laboratory,which indicates that the infl uence of direct and diff racted waves cannot be ignored for large-scale heterogeneities.The relationship between velocity and ka shows frequency dependence;the reason is that the magnitude of change in velocity caused by wave number is diff erent from that caused by heterogeneous scale.According to the change in the recorded waveform,amplitude variation,or the relationship between the velocity measured at diff erent frequencies and the heterogeneous scale,the identifi ed turning points of the ray approximation are all around ka=10.When ka is less than 1,the velocity changes slowly and gradually approaches the eff ective medium velocity.The ray velocity measured for heterogeneous media with large velocity perturbations in the laboratory is signifi cantly smaller than the velocity predicted by the perturbation theory.展开更多
Scratch test and friction test were performed to evaluate the internal and external interface behaviors of TiAISiN coating, respectively. The critical compressive and shearing stress of coating failure during scratch ...Scratch test and friction test were performed to evaluate the internal and external interface behaviors of TiAISiN coating, respectively. The critical compressive and shearing stress of coating failure during scratch test were calculated and the values are 30.84 MPa and 4.98 MPa respectively. The average friction coefficients of TiAISiN coat- ing against 2Crl2Ni4Mo3VNbN steel are 0.70 （sliding speed 50 m/rain）, 0.63 （sliding speed 100 m/min）, and 0.81 （sliding speed 150 m/min）. The elements diffusion was analyzed by EDS. A1 and Si element of coating material dif- fuse to the steel disc, except Ti element. The oxidation decreases with the increase of sliding speed, but the adhesion increases with the increase of sliding speed. More A1 element diffuses to the steel disc at the high sliding speed, but the diffusion of Si element keeps almost constant at dlfferent sliding speeds.展开更多
This paper proposes a weak Galerkin finite element method to solve incompressible quasi-Newtonian Stokes equations. We use piecewise polynomials of degrees k + 1(k 0) and k for the velocity and pressure in the interio...This paper proposes a weak Galerkin finite element method to solve incompressible quasi-Newtonian Stokes equations. We use piecewise polynomials of degrees k + 1(k 0) and k for the velocity and pressure in the interior of elements, respectively, and piecewise polynomials of degrees k and k + 1 for the boundary parts of the velocity and pressure, respectively. Wellposedness of the discrete scheme is established. The method yields a globally divergence-free velocity approximation. Optimal priori error estimates are derived for the velocity gradient and pressure approximations. Numerical results are provided to confirm the theoretical results.展开更多
Kinetic Monte Carlo methods provide a powerful computational tool for the simulation of microscopic processes such as the diffusion of interacting particles on a surface, at a detailed atomistic level. However such al...Kinetic Monte Carlo methods provide a powerful computational tool for the simulation of microscopic processes such as the diffusion of interacting particles on a surface, at a detailed atomistic level. However such algorithms are typically computationatly expensive and are restricted to fairly small spatiotemporal scales. One approach towards overcoming this problem was the development of coarse-grained Monte Carlo algorithms. In recent literature, these methods were shown to be capable of efficiently describing much larger length scales while still incorporating information on microscopic interactions and fluctuations. In this paper, a coarse-grained Langevin system of stochastic differential equations as approximations of diffusion of interacting particles is derived, based on these earlier coarse-grained models. The authors demonstrate the asymptotic equivalence of transient and long time behavior of the Langevin approximation and the underlying microscopic process, using asymptotics methods such as large deviations for interacting particles systems, and furthermore, present corresponding numerical simulations, comparing statistical quantities like mean paths, auto correlations and power spectra of the microscopic and the approximating Langevin processes. Finally, it is shown that the Langevin approximations presented here are much more computationally efficient than conventional Kinetic Monte Carlo methods, since in addition to the reduction in the number of spatial degrees of freedom in coarse-grained Monte Carlo methods, the Langevin system of stochastic differential equations allows for multiple particle moves in a single timestep.展开更多
The study of soil and rigid body system interactions is very important for the exploration of the Moon and Mars worldwide.The discrete element method(DEM)is a relatively accurate simulation method to study dry sand so...The study of soil and rigid body system interactions is very important for the exploration of the Moon and Mars worldwide.The discrete element method(DEM)is a relatively accurate simulation method to study dry sand soil mechanical properties.However,it is not suitable for bodies that are in mutual contact,connected due to constraints or have complex inertia properties due to their geometry.An efficient combination of the two-dimensional discrete element and multibody dynamic modeling method is proposed to solve the problem,in which the contacts and frictions among the granular spheres and the multibody system,including the smooth and rough rigid bodies,are taken into account.In this work,the soil field is modeled by a two-dimensional DEM,and the dynamics of the constrained rigid body system are modeled by the Cartesian method.A detection algorithm is developed to address the interactions between spherical discrete elements and roughly shaped rigid bodies.The advantage of this coupled method is that it enables the simultaneous capture of both responses.Finally,the program is verified by simulation experiments of the three-ball collision and the collision among the rectangular bars and the three balls.Based on this,the movement of the toothed wheel in the granular matter is analyzed,and the results show that the wheel with six teeth and 30°inclination has the fastest forward speed.In extraterrestrial objects,the wheel grip worsens,but the forward speed first increases and then decreases with decreasing gravity acceleration and loads on wheels,which proves that the coupled two-dimensional DEM and multibody dynamic program is effective in solving engineering problems.展开更多
基金sponsored by the National Natural Science Foundation of China(Grant Nos.40830423and40904029)CNOOC Zhanjiang Research Project(Contract No.Z2008SLZJ-FN0158)the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education Ministry
文摘Ultrasonic velocities of a set of saturated sandstone samples were measured at simulated in-situ pressures in the laboratory.The samples were obtained from the W formation of the WXS Depression and covered low to nearly high porosity and permeability ranges.The brine and four different density oils were used as pore fluids,which provided a good chance to investigate fluid viscosity-induced velocity dispersion.The analysis of experimental observations of velocity dispersion indicates that（1）the Biot model can explain most of the small discrepancy（about 2–3%）between ultrasonic measurements and zero frequency Gassmann predictions for high porosity and permeability samples saturated by all the fluids used in this experiment and is also valid for medium porosity and permeability samples saturated with low viscosity fluids（less than approximately 3 mP·S）and（2）the squirt flow mechanism dominates the low to medium porosity and permeability samples when fluid viscosity increases and produces large velocity dispersions as high as about 8%. The microfracture aspect ratios were also estimated for the reservoir sandstones and applied to calculate the characteristic frequency of the squirt flow model,above which the Gassmann’ s assumptions are violated and the measured high frequency velocities cannot be directly used for Gassmann’s fluid replacement at the exploration seismic frequency band for W formation sandstones.
基金sponsored by the National Science and Technology Major Project(Grant No.2011ZX05025-001-07)
文摘Mesoscopic fluid flow is the major cause of wave attenuation and velocity dispersion at seismic frequencies in porous rocks. The Johnson model provides solutions for the frequency-dependent quality factor and phase velocity in partially saturated porous media with pore patches of arbitrary shapes. We use the Johnson model to derive approximations for the quality factor Q at the high and low frequency limit, and obtain the approximate equation for Qmin based on geophysical and geometric parameters. A more accurate equation for Qmin is obtained after correcting for the linear errors between the exact and approximate Q values. The complexity of the pore patch shape affects the maximum attenuation of Qmin and the transition frequency ftr; furthermore, the effect on ftr is stronger than that on Qmin. Numerical solutions to Biot’s equation are computationally intensive; thus, we build an equivalent viscoelastic model on the basis of the Zener model, which well approximates the wave attenuation and dispersion in porous rocks in the seismic band.
基金supported by the National Science and Technology Major Project of China(No.2017ZX05005-004).
文摘There are complex heterogeneous entities in the underground medium,and the heterogeneous scale has a substantial impact on wave propagation.In this study,we used a set of 11 samples of glass beads as high-velocity heterogeneous bodies to evaluate the impact of such heterogeneous bodies on the propagation of P-wave.We vary the heterogeneous scale by changing the diameter of the glass beads from 0.18 to 11 mm while keeping the same volume proportion(10%)of the beads for the set of 11 samples.The pulse transmission method was used to record measurements at the ultrasonic frequencies of 0.34,0.61,and 0.84 MHz in the homogeneous matrix.The relationship between P-wave fi eld features and heterogeneity scale,P-wave velocity,and the multiple of the wave number and heterogeneous scale(ka)was observed in the laboratory,which has sparked widespread interest and research.Heterogeneous scale affects P-wave propagation,and its wave field changes are complex.The waveform,amplitude,and velocity of the recorded P-waves correlate with the heterogeneous scale.For the forward scattering while large-scale heterogeneities,noticeable direct and diff racted waves are observed in the laboratory,which indicates that the infl uence of direct and diff racted waves cannot be ignored for large-scale heterogeneities.The relationship between velocity and ka shows frequency dependence;the reason is that the magnitude of change in velocity caused by wave number is diff erent from that caused by heterogeneous scale.According to the change in the recorded waveform,amplitude variation,or the relationship between the velocity measured at diff erent frequencies and the heterogeneous scale,the identifi ed turning points of the ray approximation are all around ka=10.When ka is less than 1,the velocity changes slowly and gradually approaches the eff ective medium velocity.The ray velocity measured for heterogeneous media with large velocity perturbations in the laboratory is signifi cantly smaller than the velocity predicted by the perturbation theory.
基金Supported by China Postdoctoral Science Foundation(No. 20110490380 and No. 20110490383)Dongfang Turbine Co, Ltd (No. 2011GZ011)State Key Laboratory of Tribology, Tsinghua University (No. SKLT10A01)
文摘Scratch test and friction test were performed to evaluate the internal and external interface behaviors of TiAISiN coating, respectively. The critical compressive and shearing stress of coating failure during scratch test were calculated and the values are 30.84 MPa and 4.98 MPa respectively. The average friction coefficients of TiAISiN coat- ing against 2Crl2Ni4Mo3VNbN steel are 0.70 （sliding speed 50 m/rain）, 0.63 （sliding speed 100 m/min）, and 0.81 （sliding speed 150 m/min）. The elements diffusion was analyzed by EDS. A1 and Si element of coating material dif- fuse to the steel disc, except Ti element. The oxidation decreases with the increase of sliding speed, but the adhesion increases with the increase of sliding speed. More A1 element diffuses to the steel disc at the high sliding speed, but the diffusion of Si element keeps almost constant at dlfferent sliding speeds.
基金supported by Major Research Plan of National Natural Science Foundation of China (Grant No. 91430105)
文摘This paper proposes a weak Galerkin finite element method to solve incompressible quasi-Newtonian Stokes equations. We use piecewise polynomials of degrees k + 1(k 0) and k for the velocity and pressure in the interior of elements, respectively, and piecewise polynomials of degrees k and k + 1 for the boundary parts of the velocity and pressure, respectively. Wellposedness of the discrete scheme is established. The method yields a globally divergence-free velocity approximation. Optimal priori error estimates are derived for the velocity gradient and pressure approximations. Numerical results are provided to confirm the theoretical results.
基金supported by the National Science Foundation (Nos.DMS-0413864,DMS-0715125,DMS-0715125)the CDI-Type II Award (No.NSF-CMMI-0835673)+1 种基金the European Commision Marie-Curie Grant (No.FP6-517911)the Swedish Research Council
文摘Kinetic Monte Carlo methods provide a powerful computational tool for the simulation of microscopic processes such as the diffusion of interacting particles on a surface, at a detailed atomistic level. However such algorithms are typically computationatly expensive and are restricted to fairly small spatiotemporal scales. One approach towards overcoming this problem was the development of coarse-grained Monte Carlo algorithms. In recent literature, these methods were shown to be capable of efficiently describing much larger length scales while still incorporating information on microscopic interactions and fluctuations. In this paper, a coarse-grained Langevin system of stochastic differential equations as approximations of diffusion of interacting particles is derived, based on these earlier coarse-grained models. The authors demonstrate the asymptotic equivalence of transient and long time behavior of the Langevin approximation and the underlying microscopic process, using asymptotics methods such as large deviations for interacting particles systems, and furthermore, present corresponding numerical simulations, comparing statistical quantities like mean paths, auto correlations and power spectra of the microscopic and the approximating Langevin processes. Finally, it is shown that the Langevin approximations presented here are much more computationally efficient than conventional Kinetic Monte Carlo methods, since in addition to the reduction in the number of spatial degrees of freedom in coarse-grained Monte Carlo methods, the Langevin system of stochastic differential equations allows for multiple particle moves in a single timestep.
基金supported by the Key Program of National Natural Science Foundation of China(Grant No.11932001)the General Program of National Natural Science Foundation of China(Grant Nos.11772186 and 11772188)for which the authors are grateful.This research was also supported by the Key Laboratory of Hydrodynamics(Ministry of Education).
文摘The study of soil and rigid body system interactions is very important for the exploration of the Moon and Mars worldwide.The discrete element method(DEM)is a relatively accurate simulation method to study dry sand soil mechanical properties.However,it is not suitable for bodies that are in mutual contact,connected due to constraints or have complex inertia properties due to their geometry.An efficient combination of the two-dimensional discrete element and multibody dynamic modeling method is proposed to solve the problem,in which the contacts and frictions among the granular spheres and the multibody system,including the smooth and rough rigid bodies,are taken into account.In this work,the soil field is modeled by a two-dimensional DEM,and the dynamics of the constrained rigid body system are modeled by the Cartesian method.A detection algorithm is developed to address the interactions between spherical discrete elements and roughly shaped rigid bodies.The advantage of this coupled method is that it enables the simultaneous capture of both responses.Finally,the program is verified by simulation experiments of the three-ball collision and the collision among the rectangular bars and the three balls.Based on this,the movement of the toothed wheel in the granular matter is analyzed,and the results show that the wheel with six teeth and 30°inclination has the fastest forward speed.In extraterrestrial objects,the wheel grip worsens,but the forward speed first increases and then decreases with decreasing gravity acceleration and loads on wheels,which proves that the coupled two-dimensional DEM and multibody dynamic program is effective in solving engineering problems.
