The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LW...The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.展开更多
Aimed at the demand of contingency return at any time during the near-moon phase in the manned lunar landing missions,a fast calculation method for three-impulse contingency return trajectories is proposed.Firstly,a t...Aimed at the demand of contingency return at any time during the near-moon phase in the manned lunar landing missions,a fast calculation method for three-impulse contingency return trajectories is proposed.Firstly,a three-impulse contingency return trajectory scheme is presented by combining the Lambert transfer and maneuver at the special point.Secondly,a calculation model of three-impulse contingency return trajectories is established.Then,fast calculation methods are proposed by adopting the high-order Taylor expansion of differential algebra in the twobody trajectory dynamics model and perturbed trajectory dynamics model.Finally,the performance of the proposed methods is verified by numerical simulation.The results indicate that the fast calculation method of two-body trajectory has higher calculation efficiency compared to the semi-analytical calculation method under a certain accuracy condition.Due to its high efficiency,the characteristics of the three-impulse contingency return trajectories under different contingency scenarios are further analyzed expeditiously.These findings can be used for the design of contingency return trajectories in future manned lunar landing missions.展开更多
The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a sp...The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.展开更多
Poverty incidence is the key index that needs to be measured in the poverty exit examination and evaluation of 832 poverty-stricken counties and 128 000 poverty-stricken villages. In this paper, based on the statement...Poverty incidence is the key index that needs to be measured in the poverty exit examination and evaluation of 832 poverty-stricken counties and 128 000 poverty-stricken villages. In this paper, based on the statement of general concept and conventional calculation method of the poverty incidence, the calculation method of poverty incidence in the exit evaluation of poverty-stricken counties （also including poverty-stricken township and poverty-stricken villages） was investigated through the view of the third-party evaluation. In addition to considering the previous ＂number of planned poverty remaining population＂, the method also needed to give consideration to the exiting mistaken population, evaluation missing population. Based on the case in Yuanyang County, Honghe Hani and Yi Autonomous Prefecture, Yunnan Province, the poverty incidences of 10 exit planning villages by the end of 2017 in Yuanyang County were investigated and estimated, and suggestions were proposed to promote the precise poverty alleviation and poverty relief of the county.展开更多
In order to study the calculation methods of bending behavior of Chinese reinforced concrete beams from 1912 to 1949, tests on the mechanical performance of 66 rebars from different modem Chinese concrete buildings, t...In order to study the calculation methods of bending behavior of Chinese reinforced concrete beams from 1912 to 1949, tests on the mechanical performance of 66 rebars from different modem Chinese concrete buildings, the concrete compressive strength of 12 modem Chinese concrete buildings, and the concrete cover thickness of 9 modem Chinese concrete buildings are carried out; and the actual material properties and structural conformations of modem Chinese concrete buildings are obtained. Then, the comparison on calculation methods of bending behavior including the original Chinese calculation method, the present Chinese calculation method, the present American calculation method and the present European calculation method is studied. The results show that the original Chinese calculation method of bending behavior is based on the allowable stress calculation method, and the design safety factors are 3.55 to 4. 0. In term of the calculation area of longitudinal rebars of reinforced concrete beams, without considering earthquake action, the original Chinese structural calculation method is safer than the present Chinese structural calculation method, the present European structural calculation method, and the present American structural calculation method. The results can provide support for the structural safety assessments of modem Chinese reinforced concrete buildings.展开更多
In re-entry, the drilling riser hanging to the holding vessel takes on a free hanging state, waiting to be moved from the initial random position to the wellhead. For the re-entry, dynamics calculation is often done t...In re-entry, the drilling riser hanging to the holding vessel takes on a free hanging state, waiting to be moved from the initial random position to the wellhead. For the re-entry, dynamics calculation is often done to predict the riser motion or evaluate the structural safety. A dynamics calculation method based on Flexible Segment Model （FSM） is proposed for free hanging marine risers. In FSM, a riser is discretized into a series of flexible segments. For each flexible segment, its deflection feature and external forces are analyzed independently. For the whole riser, the nonlinear governing equations are listed according to the moment equilibrium at nodes. For the solution of the nonlinear equations, a linearization iteration scheme is provided in the paper. Owing to its flexibility, each segment can match a long part of the riser body, which enables that good results can be obtained even with a small number of segments. Moreover, the linearization iteration scheme can avoid widely used Newton-Rapson iteration scheme in which the calculation stability is influenced by the initial points. The FSM-based dynamics calculation is timesaving and stable, so suitable for the shape prediction or real-time control of free hanging marine risers.展开更多
The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper,...The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.展开更多
The method to predict roll deformation precisely and efficiently is vital for the strip shape control of a six-high rolling mill. Traditional calculation methods of roll deformation, such as the finite element method ...The method to predict roll deformation precisely and efficiently is vital for the strip shape control of a six-high rolling mill. Traditional calculation methods of roll deformation, such as the finite element method and the influence function method, have been widely used due to their accuracies. However, the required calculation time is too long to be applied to the realtime control. Therefore, a rapid calculation method for predicting roll deformation of a six-high rolling mill was proposed, which employed the finite difference method to calculate the roll deflection and used a polynomial to describe the nonlinear relationship between roll flattening and roll contact pressure. Furthermore, a new correction strategy was proposed in the iteration, where the roll center flattening and the roll flattening deviation were put forward and corrected simultaneously in the iteration process according to the static equilibrium of roll. Finally, by the comparison with traditional methods, the proposed method was proved to be more efficient and it was suitable for the online calculation of the strip shape control.展开更多
Debris flows are the one type of natural disaster that is most closely associated with hu- man activities. Debris flows are characterized as being widely distributed and frequently activated. Rainfall is an important ...Debris flows are the one type of natural disaster that is most closely associated with hu- man activities. Debris flows are characterized as being widely distributed and frequently activated. Rainfall is an important component of debris flows and is the most active factor when debris flows oc- cur. Rainfall also determines the temporal and spatial distribution characteristics of the hazards. A reasonable rainfall threshold target is essential to ensuring the accuracy of debris flow pre-warning. Such a threshold is important for the study of the mechanisms of debris flow formation, predicting the characteristics of future activities and the design of prevention and engineering control measures. Most mountainous areas have little data regarding rainfall and hazards, especially in debris flow forming re- gions. Therefore, both the traditional demonstration method and frequency calculated method cannot satisfy the debris flow pre-warning requirements. This study presents the characteristics of pre-warning regions, included the rainfall, hydrologic and topographic conditions. An analogous area with abundant data and the same conditions as the pre-warning region was selected, and the rainfall threshold was calculated by proxy. This method resolved the problem of debris flow pre-warning in ar- eas lacking data and provided a new approach for debris flow pre-warning in mountainous areas.展开更多
All-vertical-piled wharf is a kind of high-piled wharf, but it is extremely different from the traditional ones in some aspects, such as the structural property, bearing characteristics, failure mechanism, and static ...All-vertical-piled wharf is a kind of high-piled wharf, but it is extremely different from the traditional ones in some aspects, such as the structural property, bearing characteristics, failure mechanism, and static or dynamic calculation methods. In this paper, the finite element method （FEM） and theoretical analysis method are combined to analyze the structural property, bearing behavior and failure mode of the all-vertical-piled wharf in offshore deep water, and to establish simplified calculation methods determining the horizontal static ultimate bearing capacity and the dynamic response for the all-vertical-piled wharf. Firstly, the bearing capability and failure mechanism for all-vertical-piled wharf are studied by use of FEM, and the failure criterion is put forward for all-vertical-piled wharf based on the ＇plastic hinge＇. According to the failure criterion and P-Y curve method, the simplified calculation method of the horizontal static ultimate bearing capacity for all-vertical-piled wharf is proposed, and it is verified that the simplified method is reasonable by comparison with the FEM. Secondly, the displacement dynamic magnification factor for the all-vertical-piled wharf under wave cyclic loads and ship impact loads is calculated by the FEM and the theory formula based on the single degree of freedom （SDOF） system. The results obtained by the two methods are in good agreement with each other, and the simplified calculation method of the displacement dynamic magnification factor for all-vertical-piled wharf under dynamic loads is proposed. Then the simplified calculation method determining the dynamic response for the all-vertical-piled wharf is proposed in combination with P-Y curve method. That is, the dynamic response of the structure can be obtained through the static calculation results of P-Y curve method multiplied by the displacement dynamic magnification factor. The feasibility of the simplified dynamic response method is verified by comparison with the FEM under different conditions.展开更多
The Levenberg-Marquardt method, the best algorithm to obtain the least-square solution of nonlinear equations, is applied to calculate the stable phase equilibria. It can get the best combination between robustness an...The Levenberg-Marquardt method, the best algorithm to obtain the least-square solution of nonlinear equations, is applied to calculate the stable phase equilibria. It can get the best combination between robustness and speed of the calculations. Its application to ternary AI-Si-Mg system is executed in detail. The calculated phase equilibria agree well with the experimental results. Furthermore, the Levenberg-Marquardt method is not sensitive to the initial values.展开更多
In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can be...In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling （LWD）, in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.展开更多
On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the...On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software ＂Auto mapping system of submarine topography and geomorphology MBChart＂. This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China＇ s next Five-year Plan.展开更多
This paper presents a rapid and simple risk calculation method for large and complex engineering systems, the simulated maximum entropy method （SMEM）, which is based on integration of the advantages of the Monte Car...This paper presents a rapid and simple risk calculation method for large and complex engineering systems, the simulated maximum entropy method （SMEM）, which is based on integration of the advantages of the Monte Carlo and maximum entropy methods, thus avoiding the shortcoming of the slow convergence rate of the Monte Carlo method in risk calculation. Application of SMEM in the calculation of reservoir flood discharge risk shows that this method can make full use of the known information under the same conditions and obtain the corresponding probability distribution and the risk value. It not only greatly improves the speed, compared with the Monte Carlo method, but also provides a new approach for the risk calculation in large and complex engineering systems.展开更多
For a long time, because of the lack of investment capital and enough attentions, the overall constructions of rural power grid were far behind than the urban power grid in Chongqing Jiangbei Power Company. The low vo...For a long time, because of the lack of investment capital and enough attentions, the overall constructions of rural power grid were far behind than the urban power grid in Chongqing Jiangbei Power Company. The low voltage problems were highlighted in the rural power grid due to the characteristics of rural power grid. Using the distribution network flow calculation method, we evaluated the low voltage problems of the rural power grid which belongs to Chongqing Jiangbei Power Company. In addition, we collected the data of distribution transformers in electricity consumption peak period. Some practical management strategies were proposed by the analysis and evaluation of potential and appeared low voltage problems.展开更多
The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise w...The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.展开更多
A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercriti...A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of s...Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.展开更多
基金supported by the Major Research Project on Scientific Instrument Development of the National Natural Science Foundation of China(42327901)National Natural Science Foundation of China(42030806,42074120,41904104,423B2405).
文摘The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.
基金co-supported by the National Natural Science Foundation of China(No.12072365)the Technology Innovation Team of Manned Space Engineering,China。
文摘Aimed at the demand of contingency return at any time during the near-moon phase in the manned lunar landing missions,a fast calculation method for three-impulse contingency return trajectories is proposed.Firstly,a three-impulse contingency return trajectory scheme is presented by combining the Lambert transfer and maneuver at the special point.Secondly,a calculation model of three-impulse contingency return trajectories is established.Then,fast calculation methods are proposed by adopting the high-order Taylor expansion of differential algebra in the twobody trajectory dynamics model and perturbed trajectory dynamics model.Finally,the performance of the proposed methods is verified by numerical simulation.The results indicate that the fast calculation method of two-body trajectory has higher calculation efficiency compared to the semi-analytical calculation method under a certain accuracy condition.Due to its high efficiency,the characteristics of the three-impulse contingency return trajectories under different contingency scenarios are further analyzed expeditiously.These findings can be used for the design of contingency return trajectories in future manned lunar landing missions.
文摘The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.
文摘Poverty incidence is the key index that needs to be measured in the poverty exit examination and evaluation of 832 poverty-stricken counties and 128 000 poverty-stricken villages. In this paper, based on the statement of general concept and conventional calculation method of the poverty incidence, the calculation method of poverty incidence in the exit evaluation of poverty-stricken counties （also including poverty-stricken township and poverty-stricken villages） was investigated through the view of the third-party evaluation. In addition to considering the previous ＂number of planned poverty remaining population＂, the method also needed to give consideration to the exiting mistaken population, evaluation missing population. Based on the case in Yuanyang County, Honghe Hani and Yi Autonomous Prefecture, Yunnan Province, the poverty incidences of 10 exit planning villages by the end of 2017 in Yuanyang County were investigated and estimated, and suggestions were proposed to promote the precise poverty alleviation and poverty relief of the county.
基金The National Natural Science Foundation of China(No.51138002)the Foundation for the Author of National Excellent Doctoral Dissertation of PR China(No.201452)the Open Fund of Shanghai Key Laboratory of Engineering Structure Safety(No.2015-KF06)
文摘In order to study the calculation methods of bending behavior of Chinese reinforced concrete beams from 1912 to 1949, tests on the mechanical performance of 66 rebars from different modem Chinese concrete buildings, the concrete compressive strength of 12 modem Chinese concrete buildings, and the concrete cover thickness of 9 modem Chinese concrete buildings are carried out; and the actual material properties and structural conformations of modem Chinese concrete buildings are obtained. Then, the comparison on calculation methods of bending behavior including the original Chinese calculation method, the present Chinese calculation method, the present American calculation method and the present European calculation method is studied. The results show that the original Chinese calculation method of bending behavior is based on the allowable stress calculation method, and the design safety factors are 3.55 to 4. 0. In term of the calculation area of longitudinal rebars of reinforced concrete beams, without considering earthquake action, the original Chinese structural calculation method is safer than the present Chinese structural calculation method, the present European structural calculation method, and the present American structural calculation method. The results can provide support for the structural safety assessments of modem Chinese reinforced concrete buildings.
基金supported by the National Natural Science Foundation of China (Grant No. 51009092)the Doctoral Foundation of Education Ministry of China (Grant No. 20090073120013)the Scientific Research Foundation of State Education Ministry for the Returned Overseas Chinese Scholars
文摘In re-entry, the drilling riser hanging to the holding vessel takes on a free hanging state, waiting to be moved from the initial random position to the wellhead. For the re-entry, dynamics calculation is often done to predict the riser motion or evaluate the structural safety. A dynamics calculation method based on Flexible Segment Model （FSM） is proposed for free hanging marine risers. In FSM, a riser is discretized into a series of flexible segments. For each flexible segment, its deflection feature and external forces are analyzed independently. For the whole riser, the nonlinear governing equations are listed according to the moment equilibrium at nodes. For the solution of the nonlinear equations, a linearization iteration scheme is provided in the paper. Owing to its flexibility, each segment can match a long part of the riser body, which enables that good results can be obtained even with a small number of segments. Moreover, the linearization iteration scheme can avoid widely used Newton-Rapson iteration scheme in which the calculation stability is influenced by the initial points. The FSM-based dynamics calculation is timesaving and stable, so suitable for the shape prediction or real-time control of free hanging marine risers.
基金Sponsored by the National Natural Science Foundation of China (10672080)
文摘The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.
基金This work was financially supported by the National Natural Science Foundation of China （51674028）, and Fundamental Research Funds for the Central Universities （FRF-IC- 16-001）.
文摘The method to predict roll deformation precisely and efficiently is vital for the strip shape control of a six-high rolling mill. Traditional calculation methods of roll deformation, such as the finite element method and the influence function method, have been widely used due to their accuracies. However, the required calculation time is too long to be applied to the realtime control. Therefore, a rapid calculation method for predicting roll deformation of a six-high rolling mill was proposed, which employed the finite difference method to calculate the roll deflection and used a polynomial to describe the nonlinear relationship between roll flattening and roll contact pressure. Furthermore, a new correction strategy was proposed in the iteration, where the roll center flattening and the roll flattening deviation were put forward and corrected simultaneously in the iteration process according to the static equilibrium of roll. Finally, by the comparison with traditional methods, the proposed method was proved to be more efficient and it was suitable for the online calculation of the strip shape control.
基金supported by the National Natural Science Foundation of China(Nos.40830742 and 40901007)
文摘Debris flows are the one type of natural disaster that is most closely associated with hu- man activities. Debris flows are characterized as being widely distributed and frequently activated. Rainfall is an important component of debris flows and is the most active factor when debris flows oc- cur. Rainfall also determines the temporal and spatial distribution characteristics of the hazards. A reasonable rainfall threshold target is essential to ensuring the accuracy of debris flow pre-warning. Such a threshold is important for the study of the mechanisms of debris flow formation, predicting the characteristics of future activities and the design of prevention and engineering control measures. Most mountainous areas have little data regarding rainfall and hazards, especially in debris flow forming re- gions. Therefore, both the traditional demonstration method and frequency calculated method cannot satisfy the debris flow pre-warning requirements. This study presents the characteristics of pre-warning regions, included the rainfall, hydrologic and topographic conditions. An analogous area with abundant data and the same conditions as the pre-warning region was selected, and the rainfall threshold was calculated by proxy. This method resolved the problem of debris flow pre-warning in ar- eas lacking data and provided a new approach for debris flow pre-warning in mountainous areas.
基金financially supported by the Education Commission Fund of Chongqing(Grant No.KJ1500518)the Science Commission Fund of Chongqing(Grant No.cstc2016jcyjA0123)the Open Fund of Key Laboratory of Regulation Technology for Inland Waterway in Transportation Industry of Chongqing Jiaotong University(Grant No.NHHD-201506)
文摘All-vertical-piled wharf is a kind of high-piled wharf, but it is extremely different from the traditional ones in some aspects, such as the structural property, bearing characteristics, failure mechanism, and static or dynamic calculation methods. In this paper, the finite element method （FEM） and theoretical analysis method are combined to analyze the structural property, bearing behavior and failure mode of the all-vertical-piled wharf in offshore deep water, and to establish simplified calculation methods determining the horizontal static ultimate bearing capacity and the dynamic response for the all-vertical-piled wharf. Firstly, the bearing capability and failure mechanism for all-vertical-piled wharf are studied by use of FEM, and the failure criterion is put forward for all-vertical-piled wharf based on the ＇plastic hinge＇. According to the failure criterion and P-Y curve method, the simplified calculation method of the horizontal static ultimate bearing capacity for all-vertical-piled wharf is proposed, and it is verified that the simplified method is reasonable by comparison with the FEM. Secondly, the displacement dynamic magnification factor for the all-vertical-piled wharf under wave cyclic loads and ship impact loads is calculated by the FEM and the theory formula based on the single degree of freedom （SDOF） system. The results obtained by the two methods are in good agreement with each other, and the simplified calculation method of the displacement dynamic magnification factor for all-vertical-piled wharf under dynamic loads is proposed. Then the simplified calculation method determining the dynamic response for the all-vertical-piled wharf is proposed in combination with P-Y curve method. That is, the dynamic response of the structure can be obtained through the static calculation results of P-Y curve method multiplied by the displacement dynamic magnification factor. The feasibility of the simplified dynamic response method is verified by comparison with the FEM under different conditions.
基金This research is supported by the State Key Fundamental Research Project(G2000067202-1).
文摘The Levenberg-Marquardt method, the best algorithm to obtain the least-square solution of nonlinear equations, is applied to calculate the stable phase equilibria. It can get the best combination between robustness and speed of the calculations. Its application to ternary AI-Si-Mg system is executed in detail. The calculated phase equilibria agree well with the experimental results. Furthermore, the Levenberg-Marquardt method is not sensitive to the initial values.
基金supported by the National Natural Science Foundation of China(No.41204094)Science Foundation of China University of Petroleum,Beijing(No.2462015YQ0506)
文摘In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling （LWD）, in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.
基金This study was supported by Projects under contract Nos DY105 China's 0cean-03-01-01 and DY105-03-01-07the National Natural Science Foundation of China under contract No.40506017the Youth Foundation of Marine High-tech Project of China under contract No.2002AA616010.
文摘On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software ＂Auto mapping system of submarine topography and geomorphology MBChart＂. This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China＇ s next Five-year Plan.
基金supported by the National Water Pollution Control and Management Technology Major Projects(Grant No. 2009ZX07423-001)the National Natural Science Foundation of China (Grants No.51179069and 40971300)the Fundamental Research Funds for the Central Universities (Grants No.10QX43,09MG16,and 10QG23)
文摘This paper presents a rapid and simple risk calculation method for large and complex engineering systems, the simulated maximum entropy method （SMEM）, which is based on integration of the advantages of the Monte Carlo and maximum entropy methods, thus avoiding the shortcoming of the slow convergence rate of the Monte Carlo method in risk calculation. Application of SMEM in the calculation of reservoir flood discharge risk shows that this method can make full use of the known information under the same conditions and obtain the corresponding probability distribution and the risk value. It not only greatly improves the speed, compared with the Monte Carlo method, but also provides a new approach for the risk calculation in large and complex engineering systems.
文摘For a long time, because of the lack of investment capital and enough attentions, the overall constructions of rural power grid were far behind than the urban power grid in Chongqing Jiangbei Power Company. The low voltage problems were highlighted in the rural power grid due to the characteristics of rural power grid. Using the distribution network flow calculation method, we evaluated the low voltage problems of the rural power grid which belongs to Chongqing Jiangbei Power Company. In addition, we collected the data of distribution transformers in electricity consumption peak period. Some practical management strategies were proposed by the analysis and evaluation of potential and appeared low voltage problems.
文摘The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.
基金the National Science Council ot Taiwan,China for funding this research(Project no.:NSC 94-2218-E-035-011)
文摘A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
基金the National Science Council of Taiwan for funding this research (NSC 96-2221-E-019-061).
文摘Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.
