To investigate the diffusion reaction between Ti/Al solid diffusion couple, Ti/Al alternate foils formed by hot pressing were annealed at 525, 550, 575 and 600 °C for time ranging from 1 to 40 h. The experimental...To investigate the diffusion reaction between Ti/Al solid diffusion couple, Ti/Al alternate foils formed by hot pressing were annealed at 525, 550, 575 and 600 °C for time ranging from 1 to 40 h. The experimental results show that TiAl3 was the only observed phase at Ti/Al interface. The interface thermodynamics favored the preferential formation of TiAl3 in Ti/Al couple. The growth of TiAl3 layer occurred mainly towards Al foil side and exhibited a parabolic law. Using the interdiffusion coefficients calculated based on the contribution of grain boundary diffusion, the growth of TiAl3 was simulated numerically with the finite difference method, and the simulated results were in good agreement with the experimental ones.展开更多
A data gathering system is designed for the interferometric fiber optic gyroscope (IFOG) of land strapdown inertial system. IFOG is tested and the testing curve is given. The test data of IFOG are analyzed with Allan ...A data gathering system is designed for the interferometric fiber optic gyroscope (IFOG) of land strapdown inertial system. IFOG is tested and the testing curve is given. The test data of IFOG are analyzed with Allan variance method and each error coefficient is identified. Furthermore, a random drift error model for IFOG is built by the method of time series analysis. The conclusion provides supports for improving IFOG design and compensating for errors of IFOG in practice.展开更多
We apply the newly proposed double absorbing boundary condition(DABC)(Hagstrom et al., 2014) to solve the boundary reflection problem in seismic finite-difference(FD) modeling. In the DABC scheme, the local high...We apply the newly proposed double absorbing boundary condition(DABC)(Hagstrom et al., 2014) to solve the boundary reflection problem in seismic finite-difference(FD) modeling. In the DABC scheme, the local high-order absorbing boundary condition is used on two parallel artificial boundaries, and thus double absorption is achieved. Using the general 2D acoustic wave propagation equations as an example, we use the DABC in seismic FD modeling, and discuss the derivation and implementation steps in detail. Compared with the perfectly matched layer(PML), the complexity decreases, and the stability and fl exibility improve. A homogeneous model and the SEG salt model are selected for numerical experiments. The results show that absorption using the DABC is considerably improved relative to the Clayton–Engquist boundary condition and nearly the same as that in the PML.展开更多
We present a finite difference (FD) method for the simulation of seismic wave fields in fractured medium with an irregular (non-fiat) free surface which is beneficial for interpreting exploration data acquired in ...We present a finite difference (FD) method for the simulation of seismic wave fields in fractured medium with an irregular (non-fiat) free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations (in a curvilinear coordinate system) in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry (HTI). Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures.展开更多
In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid colu...In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.展开更多
For thousands of years, the North China Plain has been one of the most productive agricultural regions of the country, and the future of this region will be determined in large part by how global climatic change impac...For thousands of years, the North China Plain has been one of the most productive agricultural regions of the country, and the future of this region will be determined in large part by how global climatic change impact regional conditions and by actions taken to mitigate or adapt to impacts of climate change. It is of great importance to estimate the effects of carbon sequestration measures taken to mitigate or adapt to impacts of climate change with proper economic outcome models. This paper aims to measure the effects of the most commonly used carbon sequestration measures, fertilization and irrigation, on agricultural production in the North China Plain.展开更多
Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network mode...Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network model and measured data,a 3-D fracture network dynamic modeling method based on error analysis was proposed.Firstly,errors of four fracture volume density estimation methods(proposed by ODA,KULATILAKE,MAULDON,and SONG)and that of four fracture size estimation methods(proposed by EINSTEIN,SONG and TONON)were respectively compared,and the optimal methods were determined.Additionally,error index representing the deviation between fracture network model and measured data was established with integrated use of fractal dimension and relative absolute error(RAE).On this basis,the downhill simplex method was used to build the dynamic modeling method,which takes the minimum of error index as objective function and dynamically adjusts the fracture density and size parameters to correct the error index.Finally,the 3-D fracture network model could be obtained which meets the requirements.The proposed method was applied for 3-D fractures simulation in Miao Wei hydropower project in China for feasibility verification and the error index reduced from 2.618 to 0.337.展开更多
Nonuniform flow distribution along the radial direction usually exists in a Z-flow type radial flow adsorber,which will decrease the utilization of adsorbent and the switching time and may result in operating safety p...Nonuniform flow distribution along the radial direction usually exists in a Z-flow type radial flow adsorber,which will decrease the utilization of adsorbent and the switching time and may result in operating safety problems in cryogenic air separation.In order to improve the uniformity of the flow distribution along the radial direction in the adsorber,a differential equation is derived through pressure drop analysis in the Z-flow type radial adsorber with a cone in the middle of the central pipe.The differential equation determines the ideal cross-sectional radii of the cone along the axis.The result shows that the cross-sectional radius of the cone should gradually decrease from 0.3 m to zero along the axis to ensure that the process air is distributed uniformly in the Z-flow type radial flow adsorber and the shape of the cone is a little convex.The flow distribution without the cone in the central pipe is compared under different bed porosities.It is demonstrated that the proposed differential equation can provide theoretical support for designing Z-flow type radial flow adsorbers.展开更多
This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial...This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.展开更多
The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient e...The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient electromagnetic field. We used multiple transmitting loops instead of the traditional single transmitting loop to load the transmitting loop into Cartesian grids. We improved the method for calculating the z-component of the magnetic field based on the characteristics of full space. Then, we established the full- space 3D geoelectrical model using geological data for coalmines. In addition, the transient electromagnetic responses of water-filled goaves of variable shape at different locations were simulated by using the finite-difference time-domain (FDTD) method. Moreover, we evaluated the apparent resistivity results. The numerical modeling results suggested that the resistivity differences between the coal seam and its roof and floor greatly affect the distribution of apparent resistivity, resulting in nearly circular contours with the roadway head at the center. The actual distribution of apparent resistivity for different geoelectrical models of water in goaves was consistent with the models. However, when the goal water was located in one side, a false low-resistivity anomaly would appear on the other side owing to the full-space effect but the response was much weaker. Finally, the modeling results were subsequently confirmed by drilling, suggesting that the proposed method was effective.展开更多
Ultrasonic leaky Lamb waves are sensitive to defects and debonding in multilayer media. In this study, we use the finite-difference method to simulate the response of flexural waves in the presence of defects owing to...Ultrasonic leaky Lamb waves are sensitive to defects and debonding in multilayer media. In this study, we use the finite-difference method to simulate the response of flexural waves in the presence of defects owing to casing corrosion and rough fluctuations at the cement-formation interface. The ultrasonic obliquely incidence could effectively stimulate the flexural waves. The defects owing to casing corrosion change the amplitude of the early- arrival flexural wave, which gradually decrease with increasing defect thickness on the exterior walls and is the lowest when the defect length and wavelength were comparable. The scattering at the defects decreases the energy of flexural waves in the casing that leaks directly to fluids. For rough cement-formation interface, the early-arrival flexural waves do not change, whereas the late-arrival flexural waves have reduced amplitude owing to the scattering at rough interface.展开更多
In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transform...In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transformed coordinate system is derived. The wave field is modeled using the finite-difference method in the transformed coordinate system. The model calculation shows that this method is able to model the seismic wave field with fluctuating surface topography and achieve good results. Finally, the energy curves of the direct and reflected waves are analyzed to show that surface topography has a great influence on the seismic wave's dynamic properties.展开更多
D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration a...D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration are used to solve seismic wave equation, which means difference scheme for space domain and analytic integration for time domain. Both the principle and algorithm of this method are introduced in the paper. Based on the theory, the numerical examples prove that this hybrid method can lead to higher accuracy than the traditional finite difference method and the solution is very close to the exact one. Also the seismic modeling examples show the good performance of this method even in the case of complex surface conditions and complicated structures.展开更多
The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of eros...The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of erosion and deposition processes during debris flow movement, no entrainment assumption is unreasonable. The numerical model considering these processes is used for simulating debris flow. Raster grid networks of a digital elevation model in GIS provide a uniform grid system to describe complex topography. As the raster grid can be used as the finite difference mesh, the numerical model is solved numerically using the Leap-frog finite difference method. Finally, the simulation results can be displayed by GIS easily and used to debris flow evaluation. To illustrate this approach, the proposed methodology is applied to the Yohutagawa debris flow that occurred on 2oth October 2010, in Amami- Oshima area, Japan. The simulation results that reproduced the movement, erosion and deposition are in good agreement with the field investigation. The effectiveness of the dam in this real-ease is also verified by this approach. Comparison with the results were simulated by other models, shows that the present coupled model is more rational and effective.展开更多
The mathematical modeling for evaluation of the sphericity error is proposed with minimum radial separation center. To obtain the minimum sphericity error from the form data, a geometric approximation technique was de...The mathematical modeling for evaluation of the sphericity error is proposed with minimum radial separation center. To obtain the minimum sphericity error from the form data, a geometric approximation technique was devised. The technique regarded the least square sphere center as the initial center of the concentric spheres containing all measurement points, and then the center was moved gradually to reduce the radial separation till the minimum radial separation center was got where the constructed concentric spheres conformed to the minimum zone condition. The method was modeled firstly, then the geometric approximation process was analyzed, and finally,the software for data processing was programmed. As evaluation example, five steel balls were measured and the measurement data were processed with the developed program. The average iteration times of the approximation technique is 4.2, and on average the obtained sphericity error is 0. 529μm smaller than the least square solution,with accuracy increased by 7. 696%.展开更多
Numerical analyses of earthquake effects on the deformation, stability, and load transfer of a slope covered by deposits are traditionally based on the assumption that the slope is a continuum. It would be problem...Numerical analyses of earthquake effects on the deformation, stability, and load transfer of a slope covered by deposits are traditionally based on the assumption that the slope is a continuum. It would be problematic, however, to extend these approaches to the simulation of the slide, collapse and disintegration of the deposits under seismic loading. Contrary to this, a discrete element method (DEM) provides a means to consider large displacement and rotation of the non-continuum. To take the advantages of both methods of continuum and non- continuum analyses, seismic responses of a slope covered by deposits are studied by coupling a twodimensional (a-D) finite difference method and a 2-D DEM, with the bedrock being modelled by the finite difference grids and the deposits being represented by disks. A smooth transition across the boundaries of the continuous/discontinuous domains is obtained by imposing the compatibility condition and equilibrium condition along their interfaces. In the course of computation, the same time-step value is chosen for both continuous and discontinuous domains. The free-field boundaries are adopted for lateral grids of bedrock domain to eliminate the radiation damping effect. When the static equilibrium under gravity load is obtained, dynamic calculation begins under excitation of the seismic wave input from the continuum model bottom. In this way, responses to the earthquake of a slope covered by deposits are analyzed dynamically. Combined with field monitoring data, deformation and stability of the slope are discussed. The effects of the relevant parameters of spectrum characteristic, duration, andpeak acceleration of seismic waves are further investigated and explained from the simulations.展开更多
The overcurrent (OC) protection limit is set usually accorging to a OC protection setting table on digital integrated protection equipment in mine explode isolation high voltage (HV) vacuum switch. For digital integra...The overcurrent (OC) protection limit is set usually accorging to a OC protection setting table on digital integrated protection equipment in mine explode isolation high voltage (HV) vacuum switch. For digital integrated protection equipment, OC protection setting table must be converted to be a microcomputer algorithm. This paper first intro-duced a method of the fitting OC protection setting table to be OC relay inverse time characteristics equations using MATLAB least square fitting. On the basis of analyzing these fitting equations, a notion, “integral limit rate” was put forward initially and a OC in-verse time digital algorithm was developed. MATLAB simulation results and a digital signal processor (DSP) based digital integrated protection equipment running test indicate that this algorithm has less calculation amount, less taking up memory, high control accuracy, implements the no-grade setting of OC delay values, suits for all kinds of low-middle mi-crocomputer system implementation.展开更多
We further develop the lattice Boltzmann (LB) model [Physica A 382 (2007) 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook (BGK) collision term in the LB equation, which...We further develop the lattice Boltzmann (LB) model [Physica A 382 (2007) 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook (BGK) collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference (FLFD) scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including (i) The thermal Couette flow; (ii) One- and two-dlmenslonal FLiemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.展开更多
When modeling wave propagation in infinite space, it is necessary to have stable absorbing boundaries to effectively eliminate spurious reflections from the truncation boundaries. The SH wave equations for Perfectly M...When modeling wave propagation in infinite space, it is necessary to have stable absorbing boundaries to effectively eliminate spurious reflections from the truncation boundaries. The SH wave equations for Perfectly Matched Layers (PML) are deduced and their Crank-Nicolson scheme are presented in this paper. We use the second-, sixth-, and tenth-order finite difference and pseudo-spectral algorithms to compute the spatial derivatives. Two numerical models, a homogeneous isotropic medium and a multi-layer model with a cave, are designed to investigate how the absorbing boundary width and the algorithms determine PML effects. Numerical results show that, for PML, the low-order finite difference algorithms have fairly good absorbing effects when the absorbing boundary is thin, whereas, high-order algorithms always have good absorption when the boundary is thick. Finally, we discuss the reflection coefficient and point out its shortcomings, which is why we use the SNR to quantitatively scale the PML effects,展开更多
基金Project (50771041) supported by the National Natural Science Foundation of ChinaProject (05-0350) supported by the New Century Excellent Talents in University, China
文摘To investigate the diffusion reaction between Ti/Al solid diffusion couple, Ti/Al alternate foils formed by hot pressing were annealed at 525, 550, 575 and 600 °C for time ranging from 1 to 40 h. The experimental results show that TiAl3 was the only observed phase at Ti/Al interface. The interface thermodynamics favored the preferential formation of TiAl3 in Ti/Al couple. The growth of TiAl3 layer occurred mainly towards Al foil side and exhibited a parabolic law. Using the interdiffusion coefficients calculated based on the contribution of grain boundary diffusion, the growth of TiAl3 was simulated numerically with the finite difference method, and the simulated results were in good agreement with the experimental ones.
文摘A data gathering system is designed for the interferometric fiber optic gyroscope (IFOG) of land strapdown inertial system. IFOG is tested and the testing curve is given. The test data of IFOG are analyzed with Allan variance method and each error coefficient is identified. Furthermore, a random drift error model for IFOG is built by the method of time series analysis. The conclusion provides supports for improving IFOG design and compensating for errors of IFOG in practice.
基金supported by the National Nature Science Foundation of China(Grant No.U1262208)the Important National Science & Technology Specific Projects(Grant No.2011ZX05019-008)
文摘We apply the newly proposed double absorbing boundary condition(DABC)(Hagstrom et al., 2014) to solve the boundary reflection problem in seismic finite-difference(FD) modeling. In the DABC scheme, the local high-order absorbing boundary condition is used on two parallel artificial boundaries, and thus double absorption is achieved. Using the general 2D acoustic wave propagation equations as an example, we use the DABC in seismic FD modeling, and discuss the derivation and implementation steps in detail. Compared with the perfectly matched layer(PML), the complexity decreases, and the stability and fl exibility improve. A homogeneous model and the SEG salt model are selected for numerical experiments. The results show that absorption using the DABC is considerably improved relative to the Clayton–Engquist boundary condition and nearly the same as that in the PML.
基金sponsored by the Knowledge Innovation Program of the Chinese Academy of Sciences No.KZCX2-YW-132)the Important National Science and Technology Specific Projects(No.2008ZX05008-006)the National Natural Science Foundation of China Nos.41074033,40721003,40830315,and 40874041)
文摘We present a finite difference (FD) method for the simulation of seismic wave fields in fractured medium with an irregular (non-fiat) free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations (in a curvilinear coordinate system) in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry (HTI). Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures.
基金supported by NSFC(No.41174118)one of the major state S&T special projects(No.2008ZX05020-004)+1 种基金a Postdoctoral Fellowship of China(No.2013M530106)China Scholarship Council(No.2010644006)
文摘In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.
基金Supported by the National Scientific Foundation of China(7087311840801231)+5 种基金the Chinese Academy of Sciences(KZCX2-YW-305-2KZCX2-YW-326-1)the Ministry of Science and Technology of China(2006DFB9192012008BAC43B012008BAK47B022008BAK50B06)~~
文摘For thousands of years, the North China Plain has been one of the most productive agricultural regions of the country, and the future of this region will be determined in large part by how global climatic change impact regional conditions and by actions taken to mitigate or adapt to impacts of climate change. It is of great importance to estimate the effects of carbon sequestration measures taken to mitigate or adapt to impacts of climate change with proper economic outcome models. This paper aims to measure the effects of the most commonly used carbon sequestration measures, fertilization and irrigation, on agricultural production in the North China Plain.
基金Project(51321065)supported by the Innovative Research Groups of the National Natural Science Foundation of ChinaProject(2013CB035904)supported by the National Basic Research Program of China(973 Program)Project(51439005)supported by the National Natural Science Foundation of China
文摘Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network model and measured data,a 3-D fracture network dynamic modeling method based on error analysis was proposed.Firstly,errors of four fracture volume density estimation methods(proposed by ODA,KULATILAKE,MAULDON,and SONG)and that of four fracture size estimation methods(proposed by EINSTEIN,SONG and TONON)were respectively compared,and the optimal methods were determined.Additionally,error index representing the deviation between fracture network model and measured data was established with integrated use of fractal dimension and relative absolute error(RAE).On this basis,the downhill simplex method was used to build the dynamic modeling method,which takes the minimum of error index as objective function and dynamically adjusts the fracture density and size parameters to correct the error index.Finally,the 3-D fracture network model could be obtained which meets the requirements.The proposed method was applied for 3-D fractures simulation in Miao Wei hydropower project in China for feasibility verification and the error index reduced from 2.618 to 0.337.
基金Supported by the Major State Basic Research Development Program of China(2011CB706501)the National Natural Science Foundation of China(51176164)
文摘Nonuniform flow distribution along the radial direction usually exists in a Z-flow type radial flow adsorber,which will decrease the utilization of adsorbent and the switching time and may result in operating safety problems in cryogenic air separation.In order to improve the uniformity of the flow distribution along the radial direction in the adsorber,a differential equation is derived through pressure drop analysis in the Z-flow type radial adsorber with a cone in the middle of the central pipe.The differential equation determines the ideal cross-sectional radii of the cone along the axis.The result shows that the cross-sectional radius of the cone should gradually decrease from 0.3 m to zero along the axis to ensure that the process air is distributed uniformly in the Z-flow type radial flow adsorber and the shape of the cone is a little convex.The flow distribution without the cone in the central pipe is compared under different bed porosities.It is demonstrated that the proposed differential equation can provide theoretical support for designing Z-flow type radial flow adsorbers.
文摘This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.
基金supported by the National Key Scientific Instrument and Equipment Development Project(No.2011YQ03013307)the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education InstitutionsKey Laboratory of Coal Resources Exploration and Comprehensive Utilization,Ministry of Land and Resources
文摘The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient electromagnetic field. We used multiple transmitting loops instead of the traditional single transmitting loop to load the transmitting loop into Cartesian grids. We improved the method for calculating the z-component of the magnetic field based on the characteristics of full space. Then, we established the full- space 3D geoelectrical model using geological data for coalmines. In addition, the transient electromagnetic responses of water-filled goaves of variable shape at different locations were simulated by using the finite-difference time-domain (FDTD) method. Moreover, we evaluated the apparent resistivity results. The numerical modeling results suggested that the resistivity differences between the coal seam and its roof and floor greatly affect the distribution of apparent resistivity, resulting in nearly circular contours with the roadway head at the center. The actual distribution of apparent resistivity for different geoelectrical models of water in goaves was consistent with the models. However, when the goal water was located in one side, a false low-resistivity anomaly would appear on the other side owing to the full-space effect but the response was much weaker. Finally, the modeling results were subsequently confirmed by drilling, suggesting that the proposed method was effective.
基金supported by the Research and Development of Key Instruments and Technologies for Deep Resources Prospecting(No.ZDYZ2012-1-07)the National Natural Science Foundation of China(Nos.41204099,11134011,and 11274341)
文摘Ultrasonic leaky Lamb waves are sensitive to defects and debonding in multilayer media. In this study, we use the finite-difference method to simulate the response of flexural waves in the presence of defects owing to casing corrosion and rough fluctuations at the cement-formation interface. The ultrasonic obliquely incidence could effectively stimulate the flexural waves. The defects owing to casing corrosion change the amplitude of the early- arrival flexural wave, which gradually decrease with increasing defect thickness on the exterior walls and is the lowest when the defect length and wavelength were comparable. The scattering at the defects decreases the energy of flexural waves in the casing that leaks directly to fluids. For rough cement-formation interface, the early-arrival flexural waves do not change, whereas the late-arrival flexural waves have reduced amplitude owing to the scattering at rough interface.
基金This research is sponsored by the Scientific Research Project of the China Geological Survey "Basic Theory, Special Collection and Special Process Method Research on Metal Mineral Seismic Exploration" (Project Number: 2000201 0002146).
文摘In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transformed coordinate system is derived. The wave field is modeled using the finite-difference method in the transformed coordinate system. The model calculation shows that this method is able to model the seismic wave field with fluctuating surface topography and achieve good results. Finally, the energy curves of the direct and reflected waves are analyzed to show that surface topography has a great influence on the seismic wave's dynamic properties.
基金This project is sponsored by the Specialized Prophasic Basic Research of the"973"Programme,contract No:2001cca02300
文摘D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration are used to solve seismic wave equation, which means difference scheme for space domain and analytic integration for time domain. Both the principle and algorithm of this method are introduced in the paper. Based on the theory, the numerical examples prove that this hybrid method can lead to higher accuracy than the traditional finite difference method and the solution is very close to the exact one. Also the seismic modeling examples show the good performance of this method even in the case of complex surface conditions and complicated structures.
基金finanicial support from the Global Environment Research Fund of Japan(S-8)from Grants-in-Aid for Scientific Research(Scientific Research(B),22310113,G.Chen)from Japan Society for the Promotion of Science
文摘The objective of this study is to incorporate a numerical model with GIS to simulate the movement, erosion and deposition of debris flow across the three dimensional complex terrain. In light of the importance of erosion and deposition processes during debris flow movement, no entrainment assumption is unreasonable. The numerical model considering these processes is used for simulating debris flow. Raster grid networks of a digital elevation model in GIS provide a uniform grid system to describe complex topography. As the raster grid can be used as the finite difference mesh, the numerical model is solved numerically using the Leap-frog finite difference method. Finally, the simulation results can be displayed by GIS easily and used to debris flow evaluation. To illustrate this approach, the proposed methodology is applied to the Yohutagawa debris flow that occurred on 2oth October 2010, in Amami- Oshima area, Japan. The simulation results that reproduced the movement, erosion and deposition are in good agreement with the field investigation. The effectiveness of the dam in this real-ease is also verified by this approach. Comparison with the results were simulated by other models, shows that the present coupled model is more rational and effective.
基金Supported by National Natural Science Foundation of China(No.50175081) andTianjin Municipal Science and Technology Commission (No.0431835116).
文摘The mathematical modeling for evaluation of the sphericity error is proposed with minimum radial separation center. To obtain the minimum sphericity error from the form data, a geometric approximation technique was devised. The technique regarded the least square sphere center as the initial center of the concentric spheres containing all measurement points, and then the center was moved gradually to reduce the radial separation till the minimum radial separation center was got where the constructed concentric spheres conformed to the minimum zone condition. The method was modeled firstly, then the geometric approximation process was analyzed, and finally,the software for data processing was programmed. As evaluation example, five steel balls were measured and the measurement data were processed with the developed program. The average iteration times of the approximation technique is 4.2, and on average the obtained sphericity error is 0. 529μm smaller than the least square solution,with accuracy increased by 7. 696%.
基金the National Basic Research Program of China (Grant No. 2008CB425802)
文摘Numerical analyses of earthquake effects on the deformation, stability, and load transfer of a slope covered by deposits are traditionally based on the assumption that the slope is a continuum. It would be problematic, however, to extend these approaches to the simulation of the slide, collapse and disintegration of the deposits under seismic loading. Contrary to this, a discrete element method (DEM) provides a means to consider large displacement and rotation of the non-continuum. To take the advantages of both methods of continuum and non- continuum analyses, seismic responses of a slope covered by deposits are studied by coupling a twodimensional (a-D) finite difference method and a 2-D DEM, with the bedrock being modelled by the finite difference grids and the deposits being represented by disks. A smooth transition across the boundaries of the continuous/discontinuous domains is obtained by imposing the compatibility condition and equilibrium condition along their interfaces. In the course of computation, the same time-step value is chosen for both continuous and discontinuous domains. The free-field boundaries are adopted for lateral grids of bedrock domain to eliminate the radiation damping effect. When the static equilibrium under gravity load is obtained, dynamic calculation begins under excitation of the seismic wave input from the continuum model bottom. In this way, responses to the earthquake of a slope covered by deposits are analyzed dynamically. Combined with field monitoring data, deformation and stability of the slope are discussed. The effects of the relevant parameters of spectrum characteristic, duration, andpeak acceleration of seismic waves are further investigated and explained from the simulations.
文摘The overcurrent (OC) protection limit is set usually accorging to a OC protection setting table on digital integrated protection equipment in mine explode isolation high voltage (HV) vacuum switch. For digital integrated protection equipment, OC protection setting table must be converted to be a microcomputer algorithm. This paper first intro-duced a method of the fitting OC protection setting table to be OC relay inverse time characteristics equations using MATLAB least square fitting. On the basis of analyzing these fitting equations, a notion, “integral limit rate” was put forward initially and a OC in-verse time digital algorithm was developed. MATLAB simulation results and a digital signal processor (DSP) based digital integrated protection equipment running test indicate that this algorithm has less calculation amount, less taking up memory, high control accuracy, implements the no-grade setting of OC delay values, suits for all kinds of low-middle mi-crocomputer system implementation.
基金Supported by the Science Foundations of LCP and CAEP under Grant Nos. 2009A0102005 and 2009B0101012National Natural Science Foundation of China under Grant Nos. 11075021, 11074300, and 11074303+3 种基金National Basic Research Program (973 Program) under Grant No. 2007CB815105Fundamental Research Funds for the Central University under Grant No. 2010YS03Technology Support Program of LangFang under Grant Nos. 2010011029/30/31Science Foundation of NCIAE under Grant No. 2008-ky-13
文摘We further develop the lattice Boltzmann (LB) model [Physica A 382 (2007) 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook (BGK) collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference (FLFD) scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including (i) The thermal Couette flow; (ii) One- and two-dlmenslonal FLiemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.
基金supported jointly by the 973 Program (Grant No.2007CB209505)the National Natural Science Fund (Grant No.40704019,40674061)+1 种基金the School Basic Research Fund of Tsinghua University (JC2007030)PetroChina Innovation Fund (Grant No.060511-1-1)
文摘When modeling wave propagation in infinite space, it is necessary to have stable absorbing boundaries to effectively eliminate spurious reflections from the truncation boundaries. The SH wave equations for Perfectly Matched Layers (PML) are deduced and their Crank-Nicolson scheme are presented in this paper. We use the second-, sixth-, and tenth-order finite difference and pseudo-spectral algorithms to compute the spatial derivatives. Two numerical models, a homogeneous isotropic medium and a multi-layer model with a cave, are designed to investigate how the absorbing boundary width and the algorithms determine PML effects. Numerical results show that, for PML, the low-order finite difference algorithms have fairly good absorbing effects when the absorbing boundary is thin, whereas, high-order algorithms always have good absorption when the boundary is thick. Finally, we discuss the reflection coefficient and point out its shortcomings, which is why we use the SNR to quantitatively scale the PML effects,