In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a h...In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.展开更多
3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be effi...3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.展开更多
The oil droplet velocity in an aero engine bearing chamber can determine the initial film state which is the fundament for lubrication design and heat analysis. This paper studied the droplet motion in a respective ae...The oil droplet velocity in an aero engine bearing chamber can determine the initial film state which is the fundament for lubrication design and heat analysis. This paper studied the droplet motion in a respective aero engine bearing chamber and obtained a fast method to calculate the droplet velocity by an analytical method. Comparing the velocity results calculated by the fast method with those from the literatures by the numerical method under different operating conditions, the method proposed in this paper is confirmed to be fast and reliable. The effects of operating conditions on droplet velocity are obtained at the same time. This study contribute to follow-up research work on droplet deposition properties in aero engine bearing chambers展开更多
非量测相机以价格低、体积小、使用灵活等优势被广泛应用于高精度测量工作中,但相机标定结果对测量精度影响较大,针对现有相机标定方法存在着精度不够或标定效率不高等问题,文章提出一种联合加速分割检测特征(features from accelerated...非量测相机以价格低、体积小、使用灵活等优势被广泛应用于高精度测量工作中,但相机标定结果对测量精度影响较大,针对现有相机标定方法存在着精度不够或标定效率不高等问题,文章提出一种联合加速分割检测特征(features from accelerated segment test,FAST)算法和双像光束法平差的相机标定方法。首先利用FAST算法对像片标志点自动提取,利用直接线性变换方法解得相机内、外参数初值;然后基于固定基线长度约束的双像光束法平差模型解算相机相关参数,并通过分类阈值方法提高模型收敛效率。结果表明:该方法相较于普通双像光束法平差的精度有所提升,其标定后的中误差达到0.0064 mm,且实现了相机标定流程的半自动化,提高了作业效率,有望应用于实际场景的相机标定作业。展开更多
目的推动印刷行业朝更高效的方向迈进,提高印刷设备的易操作性和视觉识别性。方法将Kano和FAST(Function Analysis System Technique)模型引入瓦楞纸激光打印设备设计的前期应用需求分析中,通过问卷的方法获取用户的基本要求,并划分为...目的推动印刷行业朝更高效的方向迈进,提高印刷设备的易操作性和视觉识别性。方法将Kano和FAST(Function Analysis System Technique)模型引入瓦楞纸激光打印设备设计的前期应用需求分析中,通过问卷的方法获取用户的基本要求,并划分为几个子类型,进而建立Kano的二维功能属性模型。采用FAST法建立功能树,辅助使用Kano模型,从而更精准地分析用户需求,并更好地根据其需求进行优化设计。结果综合运用设计原理,针对性地挖掘瓦楞纸激光印刷设备在造型识别性、操作易用性、生产安全性上存在的问题,进而输出更优解。结论该设计方法的引入有助于为同类型的印刷设备设计提供参考,并引起更多相关厂家的重视,推动印刷行业向更积极的方向发展。展开更多
Microseismic(MS)event locations are vital aspect of MS monitoring technology used to delineate the damage zone inside the surrounding rock mass.However,complex geological conditions can impose significantly adverse ef...Microseismic(MS)event locations are vital aspect of MS monitoring technology used to delineate the damage zone inside the surrounding rock mass.However,complex geological conditions can impose significantly adverse effects on the final location results.To achieve a high-accuracy location in a complex cavern-containing structure,this study develops an MS location method using the fast marching method(FMM)with a second-order difference approach(FMM2).Based on the established velocity model with three-dimensional(3D)discrete grids,the realization of the MS location can be achieved by searching the minimum residual between the theoretical and actual first arrival times.Moreover,based on the calculation results of FMM2,the propagation paths from the MS sources to MS sensors can be obtained using the linear interpolation approach and the Runge–Kutta method.These methods were validated through a series of numerical experiments.In addition,our proposed method was applied to locate the recorded blasting and MS events that occurred during the excavation period of the underground caverns at the Houziyan hydropower station.The location results of the blasting activities show that our method can effectively reduce the location error compared with the results based on the uniform velocity model.Furthermore,the obtained MS location was verified through the occurrence of shotcrete fractures and spalling,and the monitoring results of the in-situ multipoint extensometer.Our proposed method can offer a more accurate rock fracture location and facilitate the delineation of damage zones inside the surrounding rock mass.展开更多
In 2D fast multipole method for scattering problems,square quadrature rule is used to discretize the Bessel integral identity for diagonal expansion of 2D Helmholtz kernel,and numerical integration error is introduced...In 2D fast multipole method for scattering problems,square quadrature rule is used to discretize the Bessel integral identity for diagonal expansion of 2D Helmholtz kernel,and numerical integration error is introduced. Taking advantage of the relationship between Euler-Maclaurin formula and trapezoidal quadrature rule,and the relationship between trapezoidal and square quadrature rule,sharp computable bound with analytical form on the error of numerical integration of Bessel integral identity by square quadrature rule is derived in this paper. Numerical experiments are presented at the end to demonstrate the accuracy of the sharp computable bound on the numerical integration error.展开更多
This paper describes formulation and implementation of the fast multipole boundary element method （FMBEM） for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the...This paper describes formulation and implementation of the fast multipole boundary element method （FMBEM） for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.展开更多
Fast computation of the landing footprint of a space-to-ground vehicle is a basic requirement for the deployment of parking orbits, as well as for enabling decision makers to develop real-time programs of transfer tra...Fast computation of the landing footprint of a space-to-ground vehicle is a basic requirement for the deployment of parking orbits, as well as for enabling decision makers to develop real-time programs of transfer trajectories. In order to address the usually slow computational time for the determination of the landing footprint of a space-to-ground vehicle under finite thrust, this work proposes a method that uses polynomial equations to describe the boundaries of the landing footprint and uses back propagation(BP) neural networks to quickly determine the landing footprint of the space-to-ground vehicle. First, given orbital parameters and a manoeuvre moment, the solution model of the landing footprint of a space-to-ground vehicle under finite thrust is established. Second, given arbitrary orbital parameters and an arbitrary manoeuvre moment, a fast computational model for the landing footprint of a space-to-ground vehicle based on BP neural networks is provided.Finally, the simulation results demonstrate that under the premise of ensuring accuracy, the proposed method can quickly determine the landing footprint of a space-to-ground vehicle with arbitrary orbital parameters and arbitrary manoeuvre moments. The proposed fast computational method for determining a landing footprint lays a foundation for the parking-orbit configuration and supports the design of real-time transfer trajectories.展开更多
We propose two error control techniques for numerical integrations in fast multiscale collocation methods for solving Fredholm integral equations of the second kind with weakly singular kernels. Both techniques utiliz...We propose two error control techniques for numerical integrations in fast multiscale collocation methods for solving Fredholm integral equations of the second kind with weakly singular kernels. Both techniques utilize quadratures for singular integrals using graded points. One has a polynomial order of accuracy if the integrand has a polynomial order of smoothness except at the singular point and the other has exponential order of accuracy if the integrand has an infinite order of smoothness except at the singular point. We estimate the order of convergence and computational complexity of the corresponding approximate solutions of the equation. We prove that the second technique preserves the order of convergence and computational complexity of the original collocation method. Numerical experiments are presented to illustrate the theoretical estimates.展开更多
The Unsteady Vortex Lattice Method(UVLM) is a medium-fidelity aerodynamic tool that has been widely used in aeroelasticity and flight dynamics simulations. The most timeconsuming step is the evaluation of the induced ...The Unsteady Vortex Lattice Method(UVLM) is a medium-fidelity aerodynamic tool that has been widely used in aeroelasticity and flight dynamics simulations. The most timeconsuming step is the evaluation of the induced velocity. Supposing that the number of bound and wake lattices is N and the computational cost is O (N2), we present an OeNT Dipole Panel Fast Multipole Method(DPFMM) for the rapid evaluation of the induced velocity in UVLM. The multipole expansion coefficients of a quadrilateral dipole panel have been derived in spherical coordinates, whose accuracy is the same as that of the Biot-Savart kernel at the same truncation degree P.Two methods(the loosening method and the shrinking method) are proposed and tested for space partitioning volumetric panels. Compared with FMM for vortex filaments(with three harmonics),DPFMM is approximately two times faster for N2 [103,106]. The simulation time of a multirotor(N~104) is reduced from 100 min(with unaccelerated direct solver) to 2 min(with DPFMM).展开更多
This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations s...This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.展开更多
This paper is dedicated to applying the Fourier amplitude sensitivity test(FAST)method to the problem of mixed extension and inflation of a circular cylindrical tube in the presence of residual stresses.The metafuncti...This paper is dedicated to applying the Fourier amplitude sensitivity test(FAST)method to the problem of mixed extension and inflation of a circular cylindrical tube in the presence of residual stresses.The metafunctions and the Ishigami function are considered in the sensitivity analysis(SA).The effects of the input variables on the output variables are investigated,and the most important parameters of the system under the applied pressure and axial force such as the axial stretch and the azimuthal stretch are determined.展开更多
A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditiona...A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.展开更多
A modified slow-fast analysis method is presented for the periodically excited non-autonomous dynamical system with an order gap between the exciting frequency and the natural frequency.By regarding the exciting term ...A modified slow-fast analysis method is presented for the periodically excited non-autonomous dynamical system with an order gap between the exciting frequency and the natural frequency.By regarding the exciting term as a slow-varying parameter,a generalized autonomous fast subsystem can be defined,the equilibrium branches as well as the bifurcations of which can be employed to account for the mechanism of the bursting oscillations by combining the transformed phase portrait introduced.As an example,a typical periodically excited Hartley model is used to demonstrate the validness of the method,in which the exciting frequency is far less than the natural frequency.The equilibrium branches and their bifurcations of the fast subsystem with the variation of the slow-varying parameter are presented.Bursting oscillations for two typical cases are considered,which reveals that,fold bifurcation may cause the the trajectory to jump between different equilibrium branches,while Hopf bifurcation may cause the trajectory to oscillate around the stable limit cycle.展开更多
基金supported by the Project of the National Key R&D Program(Grant No.2021YFA1000202)National Natural Science Foundation of China(Grant Nos.12120101001,12001326 and 12171283)+2 种基金Natural Science Foundation of Shandong Province(Grant Nos.ZR2021ZD03,ZR2020QA032 and ZR2019ZD42)China Postdoctoral Science Foundation(Grant Nos.BX20190191 and 2020M672038)the Startup Fund from Shandong University(Grant No.11140082063130)。
文摘In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.
基金supported by NSFC(Nos.41274120,41404085,and 41504084)
文摘3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.
基金supported by the Natural Science Foundation of China under Grant No.51275411
文摘The oil droplet velocity in an aero engine bearing chamber can determine the initial film state which is the fundament for lubrication design and heat analysis. This paper studied the droplet motion in a respective aero engine bearing chamber and obtained a fast method to calculate the droplet velocity by an analytical method. Comparing the velocity results calculated by the fast method with those from the literatures by the numerical method under different operating conditions, the method proposed in this paper is confirmed to be fast and reliable. The effects of operating conditions on droplet velocity are obtained at the same time. This study contribute to follow-up research work on droplet deposition properties in aero engine bearing chambers
文摘目的推动印刷行业朝更高效的方向迈进,提高印刷设备的易操作性和视觉识别性。方法将Kano和FAST(Function Analysis System Technique)模型引入瓦楞纸激光打印设备设计的前期应用需求分析中,通过问卷的方法获取用户的基本要求,并划分为几个子类型,进而建立Kano的二维功能属性模型。采用FAST法建立功能树,辅助使用Kano模型,从而更精准地分析用户需求,并更好地根据其需求进行优化设计。结果综合运用设计原理,针对性地挖掘瓦楞纸激光印刷设备在造型识别性、操作易用性、生产安全性上存在的问题,进而输出更优解。结论该设计方法的引入有助于为同类型的印刷设备设计提供参考,并引起更多相关厂家的重视,推动印刷行业向更积极的方向发展。
基金the Key Program of National Natural Science Foundation of China(52039007)for providing financial support.
文摘Microseismic(MS)event locations are vital aspect of MS monitoring technology used to delineate the damage zone inside the surrounding rock mass.However,complex geological conditions can impose significantly adverse effects on the final location results.To achieve a high-accuracy location in a complex cavern-containing structure,this study develops an MS location method using the fast marching method(FMM)with a second-order difference approach(FMM2).Based on the established velocity model with three-dimensional(3D)discrete grids,the realization of the MS location can be achieved by searching the minimum residual between the theoretical and actual first arrival times.Moreover,based on the calculation results of FMM2,the propagation paths from the MS sources to MS sensors can be obtained using the linear interpolation approach and the Runge–Kutta method.These methods were validated through a series of numerical experiments.In addition,our proposed method was applied to locate the recorded blasting and MS events that occurred during the excavation period of the underground caverns at the Houziyan hydropower station.The location results of the blasting activities show that our method can effectively reduce the location error compared with the results based on the uniform velocity model.Furthermore,the obtained MS location was verified through the occurrence of shotcrete fractures and spalling,and the monitoring results of the in-situ multipoint extensometer.Our proposed method can offer a more accurate rock fracture location and facilitate the delineation of damage zones inside the surrounding rock mass.
基金the National Natural Science Foundation of China (No. 11074170)the Independent Research Program of State Key Laboratory of Machinery System and Vibration (SKLMSV) (No. MSV-MS-2008-05)the Visiting Scholar Program of SKLMSV (No. MSV-2009-06)
文摘In 2D fast multipole method for scattering problems,square quadrature rule is used to discretize the Bessel integral identity for diagonal expansion of 2D Helmholtz kernel,and numerical integration error is introduced. Taking advantage of the relationship between Euler-Maclaurin formula and trapezoidal quadrature rule,and the relationship between trapezoidal and square quadrature rule,sharp computable bound with analytical form on the error of numerical integration of Bessel integral identity by square quadrature rule is derived in this paper. Numerical experiments are presented at the end to demonstrate the accuracy of the sharp computable bound on the numerical integration error.
基金Project supported by the National Natural Science Foundation of China(No.11074170)the State Key Laboratory Foundation of Shanghai Jiao Tong University(No.MSVMS201105)
文摘This paper describes formulation and implementation of the fast multipole boundary element method （FMBEM） for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.
基金supported by the National Natural Science Foundation of China (61603398)。
文摘Fast computation of the landing footprint of a space-to-ground vehicle is a basic requirement for the deployment of parking orbits, as well as for enabling decision makers to develop real-time programs of transfer trajectories. In order to address the usually slow computational time for the determination of the landing footprint of a space-to-ground vehicle under finite thrust, this work proposes a method that uses polynomial equations to describe the boundaries of the landing footprint and uses back propagation(BP) neural networks to quickly determine the landing footprint of the space-to-ground vehicle. First, given orbital parameters and a manoeuvre moment, the solution model of the landing footprint of a space-to-ground vehicle under finite thrust is established. Second, given arbitrary orbital parameters and an arbitrary manoeuvre moment, a fast computational model for the landing footprint of a space-to-ground vehicle based on BP neural networks is provided.Finally, the simulation results demonstrate that under the premise of ensuring accuracy, the proposed method can quickly determine the landing footprint of a space-to-ground vehicle with arbitrary orbital parameters and arbitrary manoeuvre moments. The proposed fast computational method for determining a landing footprint lays a foundation for the parking-orbit configuration and supports the design of real-time transfer trajectories.
基金The NNSF (10371137 and 10201034) of Chinathe Foundation (20030558008) of Doctoral Program of National Higher Education, Guangdong Provincial Natural Science Foundation (1011170) of China and the Advanced Research Foundation of Zhongshan UniversityThe US National Science Foundation (9973427 and 0312113)NSF (10371122) of China and the Chinese Academy of Sciences under the program of "Hundred Distinguished Young Chinese Scientists."
文摘We propose two error control techniques for numerical integrations in fast multiscale collocation methods for solving Fredholm integral equations of the second kind with weakly singular kernels. Both techniques utilize quadratures for singular integrals using graded points. One has a polynomial order of accuracy if the integrand has a polynomial order of smoothness except at the singular point and the other has exponential order of accuracy if the integrand has an infinite order of smoothness except at the singular point. We estimate the order of convergence and computational complexity of the corresponding approximate solutions of the equation. We prove that the second technique preserves the order of convergence and computational complexity of the original collocation method. Numerical experiments are presented to illustrate the theoretical estimates.
文摘The Unsteady Vortex Lattice Method(UVLM) is a medium-fidelity aerodynamic tool that has been widely used in aeroelasticity and flight dynamics simulations. The most timeconsuming step is the evaluation of the induced velocity. Supposing that the number of bound and wake lattices is N and the computational cost is O (N2), we present an OeNT Dipole Panel Fast Multipole Method(DPFMM) for the rapid evaluation of the induced velocity in UVLM. The multipole expansion coefficients of a quadrilateral dipole panel have been derived in spherical coordinates, whose accuracy is the same as that of the Biot-Savart kernel at the same truncation degree P.Two methods(the loosening method and the shrinking method) are proposed and tested for space partitioning volumetric panels. Compared with FMM for vortex filaments(with three harmonics),DPFMM is approximately two times faster for N2 [103,106]. The simulation time of a multirotor(N~104) is reduced from 100 min(with unaccelerated direct solver) to 2 min(with DPFMM).
基金supported by the National Natural Science Foundation of China (11172291)the National Science Foundation for Post-doctoral Scientists of China (2012M510162)the Fundamental Research Funds for the Central Universities (KB2090050024)
文摘This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.
文摘This paper is dedicated to applying the Fourier amplitude sensitivity test(FAST)method to the problem of mixed extension and inflation of a circular cylindrical tube in the presence of residual stresses.The metafunctions and the Ishigami function are considered in the sensitivity analysis(SA).The effects of the input variables on the output variables are investigated,and the most important parameters of the system under the applied pressure and axial force such as the axial stretch and the azimuthal stretch are determined.
基金Project(Z110803)supported by the State Key Laboratory of Geomechanics and Geotechnical Engineering,ChinaProject(2008AA062303)supported by the National High Technology Research and Development Program of China
文摘A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.
基金supported by the National Natural Science Foundation of China(Grants11632008 and 11872189)
文摘A modified slow-fast analysis method is presented for the periodically excited non-autonomous dynamical system with an order gap between the exciting frequency and the natural frequency.By regarding the exciting term as a slow-varying parameter,a generalized autonomous fast subsystem can be defined,the equilibrium branches as well as the bifurcations of which can be employed to account for the mechanism of the bursting oscillations by combining the transformed phase portrait introduced.As an example,a typical periodically excited Hartley model is used to demonstrate the validness of the method,in which the exciting frequency is far less than the natural frequency.The equilibrium branches and their bifurcations of the fast subsystem with the variation of the slow-varying parameter are presented.Bursting oscillations for two typical cases are considered,which reveals that,fold bifurcation may cause the the trajectory to jump between different equilibrium branches,while Hopf bifurcation may cause the trajectory to oscillate around the stable limit cycle.
