During the production,the fluid in the vicinity of the directional well enters the wellbore with different rates,leading to non-uniform flux distribution along the directional well.However,in all existing studies,it i...During the production,the fluid in the vicinity of the directional well enters the wellbore with different rates,leading to non-uniform flux distribution along the directional well.However,in all existing studies,it is oversimplified to a uniform flux distribution,which can result in inaccurate results for field applications.Therefore,this paper proposes a semi-analytical model of a directional well based on the assumption of non-uniform flux distribution.Specifically,the direction well is discretized into a carefully chosen series of linear sources,such that the complex well trajectory can be captured and the nonuniform flux distribution along the wellbore can be considered to model the three-dimensional flow behavior.By using the finite difference method,we can obtain the numerical solutions of the transient flow within the wellbore.With the aid of Green's function method,we can obtain the analytical solutions of the transient flow from the matrix to the wellbore.The complete flow behavior of a directional well is perfectly represented by coupling the above two types of transient flow.Subsequently,on the basis of the proposed model,we conduct a comprehensive analysis of the pressure transient behavior of a directional well.The computation results show that the flux variation along the direction well has a significant effect on pressure responses.In addition,the directional well in an infinite reservoir may exhibit the following flow regimes:wellbore afterflow,transition flow,inclined radial flow,elliptical flow,horizontal linear flow,and horizontal radial flow.The horizontal linear flow can be observed only if the formation thickness is much smaller than the well length.Furthermore,a dip region that appears on the pressure derivative curve indicates the three-dimensional flow behavior near the wellbore.展开更多
Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditi...Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.展开更多
In this study, a morphodynamic numerical model is established with the Regional Ocean Modeling System(ROMS)to investigate the transient behavior of sand waves under realistic sea conditions. The simulation of sand wav...In this study, a morphodynamic numerical model is established with the Regional Ocean Modeling System(ROMS)to investigate the transient behavior of sand waves under realistic sea conditions. The simulation of sand wave evolution comprises two steps: 1) a regional-scale model is configured first to simulate the ocean hydrodynamics, i.e., tides and tidal currents, and 2) the transient behavior of sand waves is simulated in a small computational domain under the time-variant currents extracted from the large model. The evolution of sand waves on the continental shelf in the Beibu Gulf is specifically investigated. The numerical results of the two-year evolution of sand waves under normal sea conditions compare well with the field survey data. The transient behavior of sand waves in individual months shows that the sand waves are more stable in April and October than that in other months, which can be selected as the windows for seabed operations. The effects of sediment properties, including settling velocity, critical shear stress and surface erosion rate, on sand wave evolution are also analyzed. Then, the typhoon-induced currents are further superimposed on the tidal currents as the extreme weather conditions. Sand waves with the average wavelength generally have more active behavior than smaller or larger sand waves. The characteristics of the evolution of sand waves in an individual typhoon process are quite different for different hydrodynamic combinations. For the storm conditions, i.e., the real combination and maximum combination cases, the sand waves experience a significant migration together with a damping in height due to the dominant suspended sediment transport. For the mild conditions, i.e., the pure tidal current and minimum combination cases, the sand waves migrate less, but the heights continue growing due to the dominant bedload transport.展开更多
The dynamical equations for a inertial reciprocating machine excited by two rotating eccentric weights were built by the matrix methodology for establishing dynamical equations of discrete systems. A mathematical mode...The dynamical equations for a inertial reciprocating machine excited by two rotating eccentric weights were built by the matrix methodology for establishing dynamical equations of discrete systems. A mathematical model of electromechanical coupling system for the machine was formed by combining the dynamical equations with the state equations of the two motors. The computer simulation to the model was performed for several values of the damping coefficient or the motor power, respectively. The substance of transient behavior of the machine is unveiled by analyzing the results of the computer simulation, and new methods are presented for diminishing the transient amplitude of the vibrating machine and improving the transient behavior. The reliable mathematical model is provided for intelligent control of the transient behavior and engineering design of the equipment.展开更多
A mathematical model of electromechanical coupling system for a planar inertial vibrating machine is built by setting up dynamical equations of discrete systems with a matrix methodology proposed. The substance of the...A mathematical model of electromechanical coupling system for a planar inertial vibrating machine is built by setting up dynamical equations of discrete systems with a matrix methodology proposed. The substance of the transient behavior of the machine is unveiled by analyzing the results of the computer simulation to the model, and new methods are presented for diminishing the transient amplitude of the vibrating machine and improving the transient behavior. The reliable mathematical model is provided for intelligent control of the transient behavior of the equipment.展开更多
The loss of hydrocarbon production caused by the dynamic behavior of the inner boundary and propped fractures under long-term production conditions has been widely reported in recent studies.However,the quantitative r...The loss of hydrocarbon production caused by the dynamic behavior of the inner boundary and propped fractures under long-term production conditions has been widely reported in recent studies.However,the quantitative relationships for the variations of the inner boundary and propped fractures have not been determined and incorporated in the semi-analytical models for the pressure and rate transient analysis.This work focuses on describing the variations of the inner boundary and propped fractures and capturing the typical characteristics from the pressure transient curves.A generalized semi-analytical model was developed to characterize the dynamic behavior of the inner boundary and propped fractures under long-term production conditions.The pressure-dependent length shrinkage coefficients,which quantify the length changes of the inner zone and propped fractures,are modified and incorporated into this multi-zone semi-analytical model.With simultaneous numerical iterations and numerical inversions in Laplace and real-time space,the transient solutions to pressure and rate behavior are determined in just a few seconds.The dynamic behavior of the inner boundary and propped fractures on transient pressure curves is divided into five periods:fracture bilinear flow(FR1),dynamic PFs flow(FR2),inner-area linear flow(FR3),dynamic inner boundary flow(FR4),and outer-area dominated linear flow(FR5).The early hump during FR2 period and a positive upward shift during FR4period are captured on the log-log pressure transient curves,reflecting the dynamic behavior of the inner boundary and propped fractures during the long-term production period.The transient pressure behavior will exhibit greater positive upward trend and the flow rate will be lower with the shrinkage of the inner boundary.The pressure derivative curve will be upward earlier as the inner boundary shrinks more rapidly.The lower permeability caused by the closure of un-propped fractures in the inner zone results in greater upward in pressure derivative curves.If the permeability loss for the dynamic behavior of the inner boundary caused by the closure of un-propped fractures is neglected,the flow rate will be overestimated in the later production period.展开更多
The mathematical system is formulated by four partial differential equations combined with initial- boundary value conditions to describe transient behavior of three-dimensional semiconductor device with heat conducti...The mathematical system is formulated by four partial differential equations combined with initial- boundary value conditions to describe transient behavior of three-dimensional semiconductor device with heat conduction. The first equation of an elliptic type is defined with respect to the electric potential, the successive two equations of convection dominated diffusion type are given to define the electron concentration and the hole concentration, and the fourth equation of heat conductor is for the temperature. The electric potential appears in the equations of electron concentration, hole concentration and the temperature in the formation of the intensity. A mass conservative numerical approximation of the electric potential is presented by using the mixed finite volume element, and the accuracy of computation of the electric intensity is improved one order. The method of characteristic fractional step difference is applied to discretize the other three equations, where the hyperbolic terms are approximated by a difference quotient in the characteristics and the diffusion terms are discretized by the method of fractional step difference. The computation of three-dimensional problem works efficiently by dividing it into three one-dimensional subproblems and every subproblem is solved by the method of speedup in parallel. Using a pair of different grids （coarse partition and refined partition）, piecewise threefold quadratic interpolation, variation theory, multiplicative commutation rule of differential operators, mathematical induction and priori estimates theory and special technique of differential equations, we derive an optimal second order estimate in L2-norm. This numerical method is valuable in the simulation of semiconductor device theoretically and actually, and gives a powerful tool to solve the international problem presented by J. Douglas, Jr.展开更多
The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally gradedint...The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally gradedinterlayer with a crack perpendicular to the boundaries. The impact load is applied on the face of the crack. Fourier transform and Laplace transform methods are used to formulate the present problem in terms of a singular integral equation in Laplace transform domain. Considering variations of parameters such as the nonhomogeneity constant, the thickness ratio and the crack length, the dynamic stress intensity factors （DSIFs） in time domain are studied and some meaningful conclusions are obtained.展开更多
This paper presents a numerical analysis of the dynamic transient behaviors of undersea cables.In this numerical study,the governing equations based on Euler-Bernoulli beam theory are adopted,and they can satisfy many...This paper presents a numerical analysis of the dynamic transient behaviors of undersea cables.In this numerical study,the governing equations based on Euler-Bernoulli beam theory are adopted,and they can satisfy many applications no matter what the magnitude of the cable tension is.The nonlinear coupled equations are solved by a popular central finite difference method,and the numerical results of transient behaviors are presented when several kinds of surrounding conditions,such as different towing speeds of surface vessel,different currents and waves with various frequencies and amplitudes,are exerted.Then a detailed comparison of the results,including the upper end tension and cable shape in time-domain,is made under the above external excitations,and finally the possible reasons for these are further explained.展开更多
This work focuses on transient thermal behavior of radial fins of rectangular,triangular and hyperbolic profiles with temperature-dependent properties.A hybrid numerical algorithm which combines differential transform...This work focuses on transient thermal behavior of radial fins of rectangular,triangular and hyperbolic profiles with temperature-dependent properties.A hybrid numerical algorithm which combines differential transformation(DTM) and finite difference(FDM) methods is utilized to theoretically study the present problem.DTM and FDM are applied to the time and space domains of the problem,respectively.The accuracy of this method solution is checked against the numerical solution.Then,the effects of some applicable parameters were studied comparatively.Since a broad range of governing parameters are investigated,the results could be useful in a number of industrial and engineering applications.展开更多
This paper proposes a fractional-order simplest chaotic system using a bi-stable locally-active memristor.The characteristics of the memristor and transient transition behaviors of the proposed system are analyzed,and...This paper proposes a fractional-order simplest chaotic system using a bi-stable locally-active memristor.The characteristics of the memristor and transient transition behaviors of the proposed system are analyzed,and this circuit is implemented digitally using ARM-based MCU.Firstly,the mathematical model of the memristor is designed,which is nonvolatile,locally-active and bi-stable.Secondly,the asymptotical stability of the fractional-order memristive chaotic system is investigated and some sufficient conditions of the stability are obtained.Thirdly,complex dynamics of the novel system are analyzed using phase diagram,Lyapunov exponential spectrum,bifurcation diagram,basin of attractor,and coexisting bifurcation,coexisting attractors are observed.All of these results indicate that this simple system contains the abundant dynamic characteristics.Moreover,transient transition behaviors of the system are analyzed,and it is found that the behaviors of transient chaotic and transient period transition alternately occur.Finally,the hardware implementation of the fractional-order bi-stable locally-active memristive chaotic system using ARM-based STM32F750 is carried out to verify the numerical simulation results.展开更多
By establishing concept an transient solutions of general nonlinear systems converging to its equilibrium set, long-time behavior of solutions for cellular neural network systems is studied. A stability condition in g...By establishing concept an transient solutions of general nonlinear systems converging to its equilibrium set, long-time behavior of solutions for cellular neural network systems is studied. A stability condition in generalized sense is obtained. This result reported has an important guide to concrete neural network designs.展开更多
The transient crevice corrosion behavior of 304 stainless steel in NaCl solution has been investigated by a multiphysics coupling model.The model considers local electrochemical reactions,transport of diffe rent speci...The transient crevice corrosion behavior of 304 stainless steel in NaCl solution has been investigated by a multiphysics coupling model.The model considers local electrochemical reactions,transport of diffe rent species,and homogeneous reactions.The moving mesh method is used to obtain the geometrical change of the crevice wall with time due to corrosion.The level set method is employed to quantitatively describe the influence of the precipitation process on electrochemical reactions.The transient crevice corrosion morphology,potential and current distributions,and pH and chloride ion concentration distributions are obtained by simulation.The effect of crevice geometry factors on the corrosion process is also discussed.The simulation results are in good agreement with the experiments,showing that the model has high reliability.展开更多
We study the transient behavior of B-ISDN bandwidth allocation strategies in tills paper. First, we present the united state transition equations and obtain the steady-state probability distribution vector by the matr...We study the transient behavior of B-ISDN bandwidth allocation strategies in tills paper. First, we present the united state transition equations and obtain the steady-state probability distribution vector by the matrixgeometric solution method. Second, we treat the problem in a uniform way that establishing the coupled differential equations describing the dynamic flow process of B-ISDN channel with the steady state probability as its initial condition and solve it numerically based on the fourth-order Runge-kutta procedure which allows a good balance between computing time and accuracy. Finally, experiment results and discussion are given to illustrate the effectes of the method presented in this paper and show that the idea and analysis way can be extended to all dynamic traffic controlling and management in B-ISDN systems.展开更多
Compared with vertical and horizontal wells, the solution and computation of transient pressure responses of slanted wells are more complex. Vertical and horizontal wells are both simplified cases of slanted wells at ...Compared with vertical and horizontal wells, the solution and computation of transient pressure responses of slanted wells are more complex. Vertical and horizontal wells are both simplified cases of slanted wells at particular inclination, so the model for slanted wells is more general and more complex than other models for vertical and horizontal wells. Many authors have studied unsteady-state flow of fluids in slanted wells and various solutions have been proposed. However, until now, few of the published results pertain to the computational efficiency. Whether in the time domain or in the Laplace domain, the computation of integration of complex functions is necessary in obtaining pressure responses of slanted wells, while the computation of the integration is complex and time-consuming. To obtain a perfect type curve the computation time is unacceptable even with an aid of high-speed computers. The purpose of this paper is to present an efficient algorithm to compute transient pressure distributions caused by slanted wells in reservoirs. Based on rigorous derivation, the transient pressure solution for slanted wells of any inclination angle is presented. Assuming an infinite-conductivity wellbore, the location of the equivalent-pressure point is determined. More importantly, according to the characteristics of the integrand in a transient pressure solution for slanted wells, the whole integral interval is partitioned into several small integral intervals, and then the method of variable substitution and the variable step-size piecewise numerical integration are employed. The amount of computation is significantly reduced and the computational efficiency is greatly improved. The algorithm proposed in this paper thoroughly solved the difficulty in the efficient and high-speed computation of transient pressure distribution of slanted wells with any inclination angle.展开更多
The self-excited attractors and hidden attractors in a memcapacitive system which has three elements are studied in this paper.The critical parameter of stable and unstable states is calculated by identifying the eige...The self-excited attractors and hidden attractors in a memcapacitive system which has three elements are studied in this paper.The critical parameter of stable and unstable states is calculated by identifying the eigenvalues of Jacobian matrix.Besides,complex dynamical behaviors are investigated in the system,such as coexisting attractors,hidden attractors,coexisting bifurcation modes,intermittent chaos,and multistability.From the theoretical analyses and numerical simulations,it is found that there are four different kinds of transient transition behaviors in the memcapacitive system.Finally,field programmable gate array(FPGA)is used to implement the proposed chaotic system.展开更多
基金the financial support provided by the National Natural Science Foundation of China(No.52104043)。
文摘During the production,the fluid in the vicinity of the directional well enters the wellbore with different rates,leading to non-uniform flux distribution along the directional well.However,in all existing studies,it is oversimplified to a uniform flux distribution,which can result in inaccurate results for field applications.Therefore,this paper proposes a semi-analytical model of a directional well based on the assumption of non-uniform flux distribution.Specifically,the direction well is discretized into a carefully chosen series of linear sources,such that the complex well trajectory can be captured and the nonuniform flux distribution along the wellbore can be considered to model the three-dimensional flow behavior.By using the finite difference method,we can obtain the numerical solutions of the transient flow within the wellbore.With the aid of Green's function method,we can obtain the analytical solutions of the transient flow from the matrix to the wellbore.The complete flow behavior of a directional well is perfectly represented by coupling the above two types of transient flow.Subsequently,on the basis of the proposed model,we conduct a comprehensive analysis of the pressure transient behavior of a directional well.The computation results show that the flux variation along the direction well has a significant effect on pressure responses.In addition,the directional well in an infinite reservoir may exhibit the following flow regimes:wellbore afterflow,transition flow,inclined radial flow,elliptical flow,horizontal linear flow,and horizontal radial flow.The horizontal linear flow can be observed only if the formation thickness is much smaller than the well length.Furthermore,a dip region that appears on the pressure derivative curve indicates the three-dimensional flow behavior near the wellbore.
基金supported by National Natural Science Foundation of China(11101244,11271231)National Tackling Key Problems Program(20050200069)Doctorate Foundation of the Ministry of Education of China(20030422047)
文摘Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.
基金financially supported by the National Natural Science Foundation of China (Grant Nos. 51579232 and 51890913)the Open Funding of State Key Laboratory of Hydraulic Engineering Simulation and Safety (Grant No. HESS-1712)。
文摘In this study, a morphodynamic numerical model is established with the Regional Ocean Modeling System(ROMS)to investigate the transient behavior of sand waves under realistic sea conditions. The simulation of sand wave evolution comprises two steps: 1) a regional-scale model is configured first to simulate the ocean hydrodynamics, i.e., tides and tidal currents, and 2) the transient behavior of sand waves is simulated in a small computational domain under the time-variant currents extracted from the large model. The evolution of sand waves on the continental shelf in the Beibu Gulf is specifically investigated. The numerical results of the two-year evolution of sand waves under normal sea conditions compare well with the field survey data. The transient behavior of sand waves in individual months shows that the sand waves are more stable in April and October than that in other months, which can be selected as the windows for seabed operations. The effects of sediment properties, including settling velocity, critical shear stress and surface erosion rate, on sand wave evolution are also analyzed. Then, the typhoon-induced currents are further superimposed on the tidal currents as the extreme weather conditions. Sand waves with the average wavelength generally have more active behavior than smaller or larger sand waves. The characteristics of the evolution of sand waves in an individual typhoon process are quite different for different hydrodynamic combinations. For the storm conditions, i.e., the real combination and maximum combination cases, the sand waves experience a significant migration together with a damping in height due to the dominant suspended sediment transport. For the mild conditions, i.e., the pure tidal current and minimum combination cases, the sand waves migrate less, but the heights continue growing due to the dominant bedload transport.
文摘The dynamical equations for a inertial reciprocating machine excited by two rotating eccentric weights were built by the matrix methodology for establishing dynamical equations of discrete systems. A mathematical model of electromechanical coupling system for the machine was formed by combining the dynamical equations with the state equations of the two motors. The computer simulation to the model was performed for several values of the damping coefficient or the motor power, respectively. The substance of transient behavior of the machine is unveiled by analyzing the results of the computer simulation, and new methods are presented for diminishing the transient amplitude of the vibrating machine and improving the transient behavior. The reliable mathematical model is provided for intelligent control of the transient behavior and engineering design of the equipment.
文摘A mathematical model of electromechanical coupling system for a planar inertial vibrating machine is built by setting up dynamical equations of discrete systems with a matrix methodology proposed. The substance of the transient behavior of the machine is unveiled by analyzing the results of the computer simulation to the model, and new methods are presented for diminishing the transient amplitude of the vibrating machine and improving the transient behavior. The reliable mathematical model is provided for intelligent control of the transient behavior of the equipment.
基金financial funding of National Natural Science Foundation of China (No.52004307)China National Petroleum Corporation (No.ZLZX2020-02-04)the Science Foundation of China University of Petroleum,Beijing (No.2462018YJRC015)。
文摘The loss of hydrocarbon production caused by the dynamic behavior of the inner boundary and propped fractures under long-term production conditions has been widely reported in recent studies.However,the quantitative relationships for the variations of the inner boundary and propped fractures have not been determined and incorporated in the semi-analytical models for the pressure and rate transient analysis.This work focuses on describing the variations of the inner boundary and propped fractures and capturing the typical characteristics from the pressure transient curves.A generalized semi-analytical model was developed to characterize the dynamic behavior of the inner boundary and propped fractures under long-term production conditions.The pressure-dependent length shrinkage coefficients,which quantify the length changes of the inner zone and propped fractures,are modified and incorporated into this multi-zone semi-analytical model.With simultaneous numerical iterations and numerical inversions in Laplace and real-time space,the transient solutions to pressure and rate behavior are determined in just a few seconds.The dynamic behavior of the inner boundary and propped fractures on transient pressure curves is divided into five periods:fracture bilinear flow(FR1),dynamic PFs flow(FR2),inner-area linear flow(FR3),dynamic inner boundary flow(FR4),and outer-area dominated linear flow(FR5).The early hump during FR2 period and a positive upward shift during FR4period are captured on the log-log pressure transient curves,reflecting the dynamic behavior of the inner boundary and propped fractures during the long-term production period.The transient pressure behavior will exhibit greater positive upward trend and the flow rate will be lower with the shrinkage of the inner boundary.The pressure derivative curve will be upward earlier as the inner boundary shrinks more rapidly.The lower permeability caused by the closure of un-propped fractures in the inner zone results in greater upward in pressure derivative curves.If the permeability loss for the dynamic behavior of the inner boundary caused by the closure of un-propped fractures is neglected,the flow rate will be overestimated in the later production period.
基金supported by the National Natural Science Foundation of China(Grant Nos.11101124 and 11271231)the National Tackling Key Problems Program for Science and Technology(Grant No.20050200069)the Doctorate Foundation of the Ministry of Education of China(Grant No.20030422047)
文摘The mathematical system is formulated by four partial differential equations combined with initial- boundary value conditions to describe transient behavior of three-dimensional semiconductor device with heat conduction. The first equation of an elliptic type is defined with respect to the electric potential, the successive two equations of convection dominated diffusion type are given to define the electron concentration and the hole concentration, and the fourth equation of heat conductor is for the temperature. The electric potential appears in the equations of electron concentration, hole concentration and the temperature in the formation of the intensity. A mass conservative numerical approximation of the electric potential is presented by using the mixed finite volume element, and the accuracy of computation of the electric intensity is improved one order. The method of characteristic fractional step difference is applied to discretize the other three equations, where the hyperbolic terms are approximated by a difference quotient in the characteristics and the diffusion terms are discretized by the method of fractional step difference. The computation of three-dimensional problem works efficiently by dividing it into three one-dimensional subproblems and every subproblem is solved by the method of speedup in parallel. Using a pair of different grids （coarse partition and refined partition）, piecewise threefold quadratic interpolation, variation theory, multiplicative commutation rule of differential operators, mathematical induction and priori estimates theory and special technique of differential equations, we derive an optimal second order estimate in L2-norm. This numerical method is valuable in the simulation of semiconductor device theoretically and actually, and gives a powerful tool to solve the international problem presented by J. Douglas, Jr.
基金the National Science Foundation for Excellent Young Investigators(10325208)the National Natural Science Foundation of China(10432030)the China Postdoctoral Science Foundation(2004036018)
文摘The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally gradedinterlayer with a crack perpendicular to the boundaries. The impact load is applied on the face of the crack. Fourier transform and Laplace transform methods are used to formulate the present problem in terms of a singular integral equation in Laplace transform domain. Considering variations of parameters such as the nonhomogeneity constant, the thickness ratio and the crack length, the dynamic stress intensity factors （DSIFs） in time domain are studied and some meaningful conclusions are obtained.
基金the National High Technology Research and Development Program (863) of China(No. 2008AA092301-1)the National Natural Science Foundation of China(No. 50909061)the Ph.D. Programs Foundation of Ministry of Education of China(No. 20070248103)
文摘This paper presents a numerical analysis of the dynamic transient behaviors of undersea cables.In this numerical study,the governing equations based on Euler-Bernoulli beam theory are adopted,and they can satisfy many applications no matter what the magnitude of the cable tension is.The nonlinear coupled equations are solved by a popular central finite difference method,and the numerical results of transient behaviors are presented when several kinds of surrounding conditions,such as different towing speeds of surface vessel,different currents and waves with various frequencies and amplitudes,are exerted.Then a detailed comparison of the results,including the upper end tension and cable shape in time-domain,is made under the above external excitations,and finally the possible reasons for these are further explained.
文摘This work focuses on transient thermal behavior of radial fins of rectangular,triangular and hyperbolic profiles with temperature-dependent properties.A hybrid numerical algorithm which combines differential transformation(DTM) and finite difference(FDM) methods is utilized to theoretically study the present problem.DTM and FDM are applied to the time and space domains of the problem,respectively.The accuracy of this method solution is checked against the numerical solution.Then,the effects of some applicable parameters were studied comparatively.Since a broad range of governing parameters are investigated,the results could be useful in a number of industrial and engineering applications.
文摘This paper proposes a fractional-order simplest chaotic system using a bi-stable locally-active memristor.The characteristics of the memristor and transient transition behaviors of the proposed system are analyzed,and this circuit is implemented digitally using ARM-based MCU.Firstly,the mathematical model of the memristor is designed,which is nonvolatile,locally-active and bi-stable.Secondly,the asymptotical stability of the fractional-order memristive chaotic system is investigated and some sufficient conditions of the stability are obtained.Thirdly,complex dynamics of the novel system are analyzed using phase diagram,Lyapunov exponential spectrum,bifurcation diagram,basin of attractor,and coexisting bifurcation,coexisting attractors are observed.All of these results indicate that this simple system contains the abundant dynamic characteristics.Moreover,transient transition behaviors of the system are analyzed,and it is found that the behaviors of transient chaotic and transient period transition alternately occur.Finally,the hardware implementation of the fractional-order bi-stable locally-active memristive chaotic system using ARM-based STM32F750 is carried out to verify the numerical simulation results.
文摘By establishing concept an transient solutions of general nonlinear systems converging to its equilibrium set, long-time behavior of solutions for cellular neural network systems is studied. A stability condition in generalized sense is obtained. This result reported has an important guide to concrete neural network designs.
基金financially supported by the National Key Research and Development Program of China(No.2017YFB0702100)the Key Research Program of Frontier Sciences,CAS(No.QYZDY-SSW-JSC012)。
文摘The transient crevice corrosion behavior of 304 stainless steel in NaCl solution has been investigated by a multiphysics coupling model.The model considers local electrochemical reactions,transport of diffe rent species,and homogeneous reactions.The moving mesh method is used to obtain the geometrical change of the crevice wall with time due to corrosion.The level set method is employed to quantitatively describe the influence of the precipitation process on electrochemical reactions.The transient crevice corrosion morphology,potential and current distributions,and pH and chloride ion concentration distributions are obtained by simulation.The effect of crevice geometry factors on the corrosion process is also discussed.The simulation results are in good agreement with the experiments,showing that the model has high reliability.
文摘We study the transient behavior of B-ISDN bandwidth allocation strategies in tills paper. First, we present the united state transition equations and obtain the steady-state probability distribution vector by the matrixgeometric solution method. Second, we treat the problem in a uniform way that establishing the coupled differential equations describing the dynamic flow process of B-ISDN channel with the steady state probability as its initial condition and solve it numerically based on the fourth-order Runge-kutta procedure which allows a good balance between computing time and accuracy. Finally, experiment results and discussion are given to illustrate the effectes of the method presented in this paper and show that the idea and analysis way can be extended to all dynamic traffic controlling and management in B-ISDN systems.
基金financial support from the special fund of China’s central government for the development of local colleges and universities―the project of national first-level discipline in Oil and Gas Engineering, the National Science Fund for Distinguished Young Scholars of China (Grant No. 51125019)the National Program on Key fundamental Research Project (973 Program, Grant No. 2011CB201005)
文摘Compared with vertical and horizontal wells, the solution and computation of transient pressure responses of slanted wells are more complex. Vertical and horizontal wells are both simplified cases of slanted wells at particular inclination, so the model for slanted wells is more general and more complex than other models for vertical and horizontal wells. Many authors have studied unsteady-state flow of fluids in slanted wells and various solutions have been proposed. However, until now, few of the published results pertain to the computational efficiency. Whether in the time domain or in the Laplace domain, the computation of integration of complex functions is necessary in obtaining pressure responses of slanted wells, while the computation of the integration is complex and time-consuming. To obtain a perfect type curve the computation time is unacceptable even with an aid of high-speed computers. The purpose of this paper is to present an efficient algorithm to compute transient pressure distributions caused by slanted wells in reservoirs. Based on rigorous derivation, the transient pressure solution for slanted wells of any inclination angle is presented. Assuming an infinite-conductivity wellbore, the location of the equivalent-pressure point is determined. More importantly, according to the characteristics of the integrand in a transient pressure solution for slanted wells, the whole integral interval is partitioned into several small integral intervals, and then the method of variable substitution and the variable step-size piecewise numerical integration are employed. The amount of computation is significantly reduced and the computational efficiency is greatly improved. The algorithm proposed in this paper thoroughly solved the difficulty in the efficient and high-speed computation of transient pressure distribution of slanted wells with any inclination angle.
基金Project supported by the National Natural Science Foundation of China(Grant No.51377124)the Science Fund for New Star of Youth Science and Technology of Shaanxi Province,China(Grant No.2016KJXX-40).
文摘The self-excited attractors and hidden attractors in a memcapacitive system which has three elements are studied in this paper.The critical parameter of stable and unstable states is calculated by identifying the eigenvalues of Jacobian matrix.Besides,complex dynamical behaviors are investigated in the system,such as coexisting attractors,hidden attractors,coexisting bifurcation modes,intermittent chaos,and multistability.From the theoretical analyses and numerical simulations,it is found that there are four different kinds of transient transition behaviors in the memcapacitive system.Finally,field programmable gate array(FPGA)is used to implement the proposed chaotic system.
