Most of the carbonate formation are highly heterogeneous with cavities of different sizes, which makes the prediction of cavity-filled reservoir in carbonate rocks difficult. Large cavities in carbonate formations pos...Most of the carbonate formation are highly heterogeneous with cavities of different sizes, which makes the prediction of cavity-filled reservoir in carbonate rocks difficult. Large cavities in carbonate formations pose serious threat to drilling operations. Logging-whiledrilling （LWD） is currently used to accurately identify and evaluate cavities in reservoirs during drilling. In this study, we use the self-adaptive hp-FEM algorithm simulate and calculate the LWD resistivity responses of fracture-cavity reservoir cavities. Compared with the traditional h-FEM method, the self-adaptive hp-FEM algorithm has the characteristics of the self-adaptive mesh refinement and the calculations exponentially converge to highly accurate solutions. Using numerical simulations, we investigated the effect of the cavity size, distance between cavity and borehole, and transmitted frequency on the LWD resistivity response. Based on the results, a method for recognizing cavities is proposed. This research can provide the theoretical basis for the accurate identification and quantitative evaluation of various carbonate reservoirs with cavities encountered in practice.展开更多
The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and ...The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and it is explicit in the time domain. Consequently it is a best mixture of FEM and finite volume method （FVM）. RK-DGFEM can obtain local high-order accuracy by using high-order polynomial basis. Numerical experiments of transverse magnetic （TM） wave propagation in a 2-D resonator are performed. A high-order Lagrange polynomial basis is adopted. Numerical results agree well with analytical solution. And different order Lagrange interpolation polynomial basis impacts on simulation result accuracy are discussed. Computational results indicate that the accuracy is evidently improved when the order of interpolation basis is increased. Finally, L^2 errors of different order polynomial basis in RK-DGFEM are presented. Computational results show that L^2 error declines exponentially as the order of basis increases.展开更多
A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be ...A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be applied to a variety of velocity potential problems. The present paper, however, focused on its application to the problem of water entry of a wedge with varying speed. The continuity of the velocity achieved herein is particularly important for this kind of nonlinear free surface flow problem, because when the time stepping method is used, the free surface is updated through the velocity obtained at each node and the accuracy of the velocity is therefore crucial. Calculation was made for a case when the distance S that the wedge has travelled and time t follow the relationship s=Dtα, where D and α are constants, which is found to lead to a self similar flow field when the effect due to gravity is ignored.展开更多
To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed...To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.展开更多
A general multiple-loop feedback approach for realization of Four-Terminal Floating Nullor RC(FTFN-RC) filter is presented.The proposed filter is constructed by multi-output FTFNs, capacitors and resistors.It can simu...A general multiple-loop feedback approach for realization of Four-Terminal Floating Nullor RC(FTFN-RC) filter is presented.The proposed filter is constructed by multi-output FTFNs, capacitors and resistors.It can simultaneously realize slow-pass, band-pass(if order is even number), and high-pass filter responses.With RC elements grounded and requiring no component matching con-straints, it is fully integrated conveniently.Simulations are performed for the fourth-order Butterworth filter to verify the validity of the circuit.展开更多
Several algorithms were proposed relating to the development of a framework of the perturbation-based stochastic finite element method (PSFEM) for large variation nonlinear dynamic problems. For this purpose, algorith...Several algorithms were proposed relating to the development of a framework of the perturbation-based stochastic finite element method (PSFEM) for large variation nonlinear dynamic problems. For this purpose, algorithms and a framework related to SFEM based on the stochastic virtual work principle were studied. To prove the validity and practicality of the algorithms and framework, numerical examples for nonlinear dynamic problems with large variations were calculated and compared with the Monte-Carlo Simulation method. This comparison shows that the proposed approaches are accurate and effective for the nonlinear dynamic analysis of structures with random parameters.展开更多
Higher-order Time Domain Finite Element Method (TDFEM) based on the nodal inter- polation is proposed for two-dimensional electromagnetic analysis. The detailed algorithms of the method are presented firstly, and then...Higher-order Time Domain Finite Element Method (TDFEM) based on the nodal inter- polation is proposed for two-dimensional electromagnetic analysis. The detailed algorithms of the method are presented firstly, and then the accuracy, CPU time and memory consumption of the higher-order node-based TDFEM are investigated. The high performance of the presented approach is validated by numerical results of the transient responses of Transverse Electric (TE) field and Transverse Magnetic (TM) field in a rectangular waveguide.展开更多
A method for static aeroelastic analysis based on the high-order panel method and modal method is presented. The static aeroelastic characteristics of flexible wings are investigated using this method. Three-dimension...A method for static aeroelastic analysis based on the high-order panel method and modal method is presented. The static aeroelastic characteristics of flexible wings are investigated using this method. Three-dimensional aerodynamic models of flexible wings are constructed based on the geometry of wing configuration, and the modal method is adopted to achieve the fluid-structure coupling. The static aeroelastic characteristics of the AGARD445.6 wing and a low-aspect-ratio wing are investigated in this study. The influences of elastic structural deformation on aerodynamic forces are studied with an emphasis analyzing the aerodynamic coefficients, wing root loads, structural deformation and pressure distribution of different sections, and results are compared with the results from wind-tunnel tests and the elastic results based on experimental aerodynamic forces. It is concluded that aerodynamic forces can be accurately calculated with the high-order panel method. The method presented in this study is feasible, credible and efficient. Comprehensive static aeroelastic characteristics can be provided by the method for early phases of aircraft design.展开更多
In order to obtain the robust high-resolution beamforming, a high order cross sensor processing(CSP) approach is developed. According to the relation ship between the target bearing and the phase difference of each el...In order to obtain the robust high-resolution beamforming, a high order cross sensor processing(CSP) approach is developed. According to the relation ship between the target bearing and the phase difference of each element receiving signal, this method exploits the property that the same diagonal of covariance matrix with the same phase difference and obtains(2M-1)(N-1)virtual elements(N is the original array number) by executing M order CSP. The extended virtual elements can effectively increase the physical aperture of linear array, reduce the main lobe width of beam-forming, and improve the bearing resolution. The CSP method accumulates the data on the same sub-diagonal of the covariance matrix, which can decrease the impact of background noise on beam-forming. The theoretical analysis and experimental results both show that this method has high resolution in bearing estimation, compared with the MUSIC method, which has better robustness under the lower signal-to-noise ratio(SNR).展开更多
In this paper,the discontinuous Galerkin(DG)method combined with localized artificial diffusivity is investigated in the context of numerical simulation of broadband compressible turbulent flows with shocks for under-...In this paper,the discontinuous Galerkin(DG)method combined with localized artificial diffusivity is investigated in the context of numerical simulation of broadband compressible turbulent flows with shocks for under-resolved cases.Firstly,the spectral property of the DG method is analyzed using the approximate dispersion relation(ADR)method and compared with typical finite difference methods,which reveals quantitatively that significantly less grid points can be used with DG for comparable numerical error.Then several typical test cases relevant to problems of compressible turbulence are simulated,including one-dimensional shock/entropy wave interaction,two-dimensional decaying isotropic turbulence,and two-dimensional temporal mixing layers.Numerical results indicate that higher numerical accuracy can be achieved on the same number of degrees of freedom with DG than high order finite difference schemes.Furthermore,shocks are also well captured using the localized artificial diffusivity method.The results in this work can provide useful guidance for further applications of DG to direct and large eddy simulation of compressible turbulent flows.展开更多
The approach of nonconforming finite element method admits users to solve the partial differential equations with lower complexity,but the accuracy is usually low.In this paper,we present a family of highaccuracy nonc...The approach of nonconforming finite element method admits users to solve the partial differential equations with lower complexity,but the accuracy is usually low.In this paper,we present a family of highaccuracy nonconforming finite element methods for fourth order problems in arbitrary dimensions.The finite element methods are given in a unified way with respect to the dimension.This is an effort to reveal the balance between the accuracy and the complexity of finite element methods.展开更多
基金supported by the National Natural Science Foundation of China(No. 41074099)
文摘Most of the carbonate formation are highly heterogeneous with cavities of different sizes, which makes the prediction of cavity-filled reservoir in carbonate rocks difficult. Large cavities in carbonate formations pose serious threat to drilling operations. Logging-whiledrilling （LWD） is currently used to accurately identify and evaluate cavities in reservoirs during drilling. In this study, we use the self-adaptive hp-FEM algorithm simulate and calculate the LWD resistivity responses of fracture-cavity reservoir cavities. Compared with the traditional h-FEM method, the self-adaptive hp-FEM algorithm has the characteristics of the self-adaptive mesh refinement and the calculations exponentially converge to highly accurate solutions. Using numerical simulations, we investigated the effect of the cavity size, distance between cavity and borehole, and transmitted frequency on the LWD resistivity response. Based on the results, a method for recognizing cavities is proposed. This research can provide the theoretical basis for the accurate identification and quantitative evaluation of various carbonate reservoirs with cavities encountered in practice.
文摘The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and it is explicit in the time domain. Consequently it is a best mixture of FEM and finite volume method （FVM）. RK-DGFEM can obtain local high-order accuracy by using high-order polynomial basis. Numerical experiments of transverse magnetic （TM） wave propagation in a 2-D resonator are performed. A high-order Lagrange polynomial basis is adopted. Numerical results agree well with analytical solution. And different order Lagrange interpolation polynomial basis impacts on simulation result accuracy are discussed. Computational results indicate that the accuracy is evidently improved when the order of interpolation basis is increased. Finally, L^2 errors of different order polynomial basis in RK-DGFEM are presented. Computational results show that L^2 error declines exponentially as the order of basis increases.
文摘A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be applied to a variety of velocity potential problems. The present paper, however, focused on its application to the problem of water entry of a wedge with varying speed. The continuity of the velocity achieved herein is particularly important for this kind of nonlinear free surface flow problem, because when the time stepping method is used, the free surface is updated through the velocity obtained at each node and the accuracy of the velocity is therefore crucial. Calculation was made for a case when the distance S that the wedge has travelled and time t follow the relationship s=Dtα, where D and α are constants, which is found to lead to a self similar flow field when the effect due to gravity is ignored.
基金Supported by the National Natural Science Foundation of China under (Grant No.107 72040,50709005 and 50921001)the Major National Science and Technology Projects of China under (Grant No.2008ZX05026-02)the Open Fund of State Key Laboratory of Ocean Engineering
文摘To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.
基金Supported by the Hunan Province Project of Education Department of Financial Aid (No.04C346)
文摘A general multiple-loop feedback approach for realization of Four-Terminal Floating Nullor RC(FTFN-RC) filter is presented.The proposed filter is constructed by multi-output FTFNs, capacitors and resistors.It can simultaneously realize slow-pass, band-pass(if order is even number), and high-pass filter responses.With RC elements grounded and requiring no component matching con-straints, it is fully integrated conveniently.Simulations are performed for the fourth-order Butterworth filter to verify the validity of the circuit.
文摘Several algorithms were proposed relating to the development of a framework of the perturbation-based stochastic finite element method (PSFEM) for large variation nonlinear dynamic problems. For this purpose, algorithms and a framework related to SFEM based on the stochastic virtual work principle were studied. To prove the validity and practicality of the algorithms and framework, numerical examples for nonlinear dynamic problems with large variations were calculated and compared with the Monte-Carlo Simulation method. This comparison shows that the proposed approaches are accurate and effective for the nonlinear dynamic analysis of structures with random parameters.
基金Supported by National Natural Science Foundation of China (No. 60601024)
文摘Higher-order Time Domain Finite Element Method (TDFEM) based on the nodal inter- polation is proposed for two-dimensional electromagnetic analysis. The detailed algorithms of the method are presented firstly, and then the accuracy, CPU time and memory consumption of the higher-order node-based TDFEM are investigated. The high performance of the presented approach is validated by numerical results of the transient responses of Transverse Electric (TE) field and Transverse Magnetic (TM) field in a rectangular waveguide.
基金supported by the National Natural Science Foundation of China (Grant Nos. 60736025, 90716006 and 10902006)the Doctoral Program Foundation of Institutions of Higher Education of China (Grant No. 20091102110015)the Major Programs of China National Space Administration (Grant No. D2120060013)
文摘A method for static aeroelastic analysis based on the high-order panel method and modal method is presented. The static aeroelastic characteristics of flexible wings are investigated using this method. Three-dimensional aerodynamic models of flexible wings are constructed based on the geometry of wing configuration, and the modal method is adopted to achieve the fluid-structure coupling. The static aeroelastic characteristics of the AGARD445.6 wing and a low-aspect-ratio wing are investigated in this study. The influences of elastic structural deformation on aerodynamic forces are studied with an emphasis analyzing the aerodynamic coefficients, wing root loads, structural deformation and pressure distribution of different sections, and results are compared with the results from wind-tunnel tests and the elastic results based on experimental aerodynamic forces. It is concluded that aerodynamic forces can be accurately calculated with the high-order panel method. The method presented in this study is feasible, credible and efficient. Comprehensive static aeroelastic characteristics can be provided by the method for early phases of aircraft design.
基金supported by the National Natural Science Foundation of China(61372180)the Young Talent Frontier Project of Institute of Acoustics of Chinese Academy of Sciences(Y454341261)
文摘In order to obtain the robust high-resolution beamforming, a high order cross sensor processing(CSP) approach is developed. According to the relation ship between the target bearing and the phase difference of each element receiving signal, this method exploits the property that the same diagonal of covariance matrix with the same phase difference and obtains(2M-1)(N-1)virtual elements(N is the original array number) by executing M order CSP. The extended virtual elements can effectively increase the physical aperture of linear array, reduce the main lobe width of beam-forming, and improve the bearing resolution. The CSP method accumulates the data on the same sub-diagonal of the covariance matrix, which can decrease the impact of background noise on beam-forming. The theoretical analysis and experimental results both show that this method has high resolution in bearing estimation, compared with the MUSIC method, which has better robustness under the lower signal-to-noise ratio(SNR).
基金supported by the National Basic Research Program of China(Grant No.2009CB724104)
文摘In this paper,the discontinuous Galerkin(DG)method combined with localized artificial diffusivity is investigated in the context of numerical simulation of broadband compressible turbulent flows with shocks for under-resolved cases.Firstly,the spectral property of the DG method is analyzed using the approximate dispersion relation(ADR)method and compared with typical finite difference methods,which reveals quantitatively that significantly less grid points can be used with DG for comparable numerical error.Then several typical test cases relevant to problems of compressible turbulence are simulated,including one-dimensional shock/entropy wave interaction,two-dimensional decaying isotropic turbulence,and two-dimensional temporal mixing layers.Numerical results indicate that higher numerical accuracy can be achieved on the same number of degrees of freedom with DG than high order finite difference schemes.Furthermore,shocks are also well captured using the localized artificial diffusivity method.The results in this work can provide useful guidance for further applications of DG to direct and large eddy simulation of compressible turbulent flows.
基金supported by National Natural Science Foundation of China (Grant No.11101415)the National Center for Mathematics and Interdisciplinary Sciences,CAS
文摘The approach of nonconforming finite element method admits users to solve the partial differential equations with lower complexity,but the accuracy is usually low.In this paper,we present a family of highaccuracy nonconforming finite element methods for fourth order problems in arbitrary dimensions.The finite element methods are given in a unified way with respect to the dimension.This is an effort to reveal the balance between the accuracy and the complexity of finite element methods.
