One of the most important parameters for oceanic internal waves (IWs) is their amplitude. We have developed a method to retrieve the IW amplitude from nautical X-Band radar images based on the KdV equation for continu...One of the most important parameters for oceanic internal waves (IWs) is their amplitude. We have developed a method to retrieve the IW amplitude from nautical X-Band radar images based on the KdV equation for continuous stratified finite depth system. We have also tested the method of measuring the amplitude of IWs from X-Band radar backscatter image sequences acquired on June 2009 in the northeastern South China Sea. The method was applied in several radar images. Experiments show that the retrieval amplitudes are consistent with the in-situ observational amplitudes of IWs by using the towed thermistor chain and conductivity-temperature-depth (CTD) profile. The uncertainty of the method is also discussed.展开更多
In order to forecast the distribution of crest amplitudes and the occurrence of freak waves in a short crested coastal sea,a novel transformed linear simulation method is initially proposed in this paper.A Hermite tra...In order to forecast the distribution of crest amplitudes and the occurrence of freak waves in a short crested coastal sea,a novel transformed linear simulation method is initially proposed in this paper.A Hermite transformation model expressed as a monotonic cubic polynomial serves as the foundation for the novel simulation technique.The wave crest amplitude exceedance probabilities of two sea states-one with a directional wave spectrum based on the measured wave elevation data at the Yura coast and the other with a typical directional JONSWAP wave spectrum-have been predicted using the novel simulation method that has been proposed.The likelihood that a particular critical wave crest amplitude will be exceeded is directly correlated with the probability that freak waves will occur.It is shown that the novel simulation approach suggested can provide predictions that are more precise than those obtained from the Rayleigh crest amplitude distribution model,the Jahns and Wheeler crest amplitude distribution model,or the conventional linear simulation method.This study also demonstrated that the nonlinear simulation method is less effective than the novel simulation method in terms of efficiency.展开更多
A new theoretical solution is presented here for the dynamic characteristics of a buoyant jet due to opposing small amplitude waves. The conservation equations of mass, tangential moment^n and vertical momentum are so...A new theoretical solution is presented here for the dynamic characteristics of a buoyant jet due to opposing small amplitude waves. The conservation equations of mass, tangential moment^n and vertical momentum are solved by the integral method which encompasses the Gaussian profiles of velocity and density. The action of waves is incorporated into the equations of motion as an external force and a new exact solution is obtained to predict the trajectory, velocity distribution and boundary thickness of the buoyant jet over an arbitrary lateral cross section. It is found that the velocity along the centerline is inversely proportional to the ratio of the momentum of the wave to the buoyant jet. The averaged bound- ary width varies with the fluctuation of the boundary width, the distance from the orifice and the velocity correction function. Owing to the motion of waves, the fluctuation of the boundary width is proportional to the wave steepness.展开更多
The wave propagation in the one-dimensional complex Ginzbur-Landau equation (CGLE) is studied by considering a wave source at the system boundary. A special propagation region, which is an island-shaped zone surroun...The wave propagation in the one-dimensional complex Ginzbur-Landau equation (CGLE) is studied by considering a wave source at the system boundary. A special propagation region, which is an island-shaped zone surrounded by the defect turbulence in the system parameter space, is observed in our numerical experiment. The wave signal spreads in the whole space with a novel amplitude wave pattern in the area. The relevant factors of the pattern formation, such as the wave speed, the maximum propagating distance and the oscillatory frequency, are studied in detail. The stability and the generality of the region are testified by adopting various initial conditions. This finding of the amplitude pattern extends the wave propagation region in the parameter space and presents a new signal transmission mode, and is therefore expected to be of much importance.展开更多
We theoretically analyze the organizing filament of small amplitude scroll waves in general excitable media by perturbation method and explicitly give the expressions of coefficients in Keener theory. In particular fo...We theoretically analyze the organizing filament of small amplitude scroll waves in general excitable media by perturbation method and explicitly give the expressions of coefficients in Keener theory. In particular for the excitable media with equal diffusion, we obtain a close system for the motion of the filament. With an example of the Oregonator model, our results are in good agreement with those simulated by Winfree.展开更多
A new method of detecting the vertical bearing capacity for single-pile with high strain is discussed in this paper. A heavy hammer or a small type of rocket is used to strike the pile top and the detectors are used ...A new method of detecting the vertical bearing capacity for single-pile with high strain is discussed in this paper. A heavy hammer or a small type of rocket is used to strike the pile top and the detectors are used to record vibra- tion graphs. An expression of higher degree of strain (deformation force) is introduced. It is testified theoretically that the displacement, velocity and acceleration cannot be obtained by simple integral acceleration and differential velocity when long displacement and high strain exist, namely when the pile phase generates a whole slip relative to the soil body. That is to say that there are non-linear relations between them. It is educed accordingly that the force P and displacement S are calculated from the amplitude of wave train and (dynamic) P-S curve is drew so as to determine the yield points. Further, a method of determining the vertical bearing capacity for single-pile is dis- cussed. A static load test is utilized to check the result of dynamic test and determine the correlative constants of dynamic-static P(Q)-S curve.展开更多
In this paper we study the problem of generation of surface waves produced due to a) rolling of the plate and b) presence of a line source in front of a fixed vertical plate. The amplitudes of radiated waves at larg...In this paper we study the problem of generation of surface waves produced due to a) rolling of the plate and b) presence of a line source in front of a fixed vertical plate. The amplitudes of radiated waves at large distance from the plate, in both cases, are obtained by a suitable application of Green's integral theorem. These are then studied graphically for various values of the ice cover parameter.展开更多
This paper gives an overall discussion about water level change on slopes under wave action, including wave runup, wave rundown and wave up-down amplitude, and a suggested formula for their calculation.
The waves driven by the wind do not move on the water as ordinarily done by sailboats. Indeed, the movement of the waves driven by the wind is more complex than the sailboats’ translation movement that we know. The m...The waves driven by the wind do not move on the water as ordinarily done by sailboats. Indeed, the movement of the waves driven by the wind is more complex than the sailboats’ translation movement that we know. The movement of the wave in our particular case results from the chain-job done by wind’s stress and gravity forces: material is collected upstream (erosion phenomenon) and then deposited on the wave’s summit by the wind. This material deposited on the summit of the wave by the wind is then removed and dispatched on the downstream side of the wave by gravity forces. As always happens in any chain-job: if the wind works faster than gravity forces, great accumulation of material will occur at the summit of the wave that will lead to an increase in its (the wave in this case) height. If conversely the wind works more slowly, a deficit in material delivery will occur and gravity force goes directly to remove material on the wave’s summit and lead to a decrease in its height. In terms of Mechanics, we know that the main obstacle that can seriously disturb the work of the wind is the unavailability of water or so its viscosity. Given the complexity of the process to be studied, it seemed necessary for us to make a use of modulational instability theories such as the standard NLSE in order to better understand the contribution of wind and water viscosity to modulations of driven waves’ amplitudes (or phases): modulations which sometimes can accidentally trigger unpredictable rogue waves.展开更多
Based on the Burgers equation and Manley-Rowe equation, the derivation about nonlinear interaction of the acoustic waves has been done in this paper. After nonlinear interaction among the low-frequency weak waves and ...Based on the Burgers equation and Manley-Rowe equation, the derivation about nonlinear interaction of the acoustic waves has been done in this paper. After nonlinear interaction among the low-frequency weak waves and the pump wave, the analytical solutions of acoustic waves' amplitude in the field are deduced. The relationship between normalized energy of high-frequency and the change of acoustic energy before and after the nonlinear interaction of the acoustic waves is analyzed. The experimental results about the changes of the acoustic energy are presented. The study shows that new frequencies are generated and the energies of the low-frequency are modulated in a long term by the pump waves, which leads the energies of the low-frequency acoustic waves to change in the pulse trend in the process of the nonlinear interaction of the acoustic waves. The increase and decrease of the energies of the low-frequency are observed under certain typical conditions, which lays a foundation for practical engineering applications.展开更多
Interracial internal waves in a three-layer density-stratified fluid are investigated using a singular perturbation method, and third-order asymptotic solutions of the velocity potentials and third-order Stokes wave s...Interracial internal waves in a three-layer density-stratified fluid are investigated using a singular perturbation method, and third-order asymptotic solutions of the velocity potentials and third-order Stokes wave solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory. As expected, the third-order solutions describe the third-order nonlinear modification and the third-order nonlinear interactions between the interracial waves. The wave velocity depends on not only the wave number and the depth of each layer but also on the wave amplitude.展开更多
Fractures in oil and gas reservoirs have been the topic of many studies and have attracted reservoir research all over the world. Because of the complexities of the fractures, it is difficult to use fractured reservoi...Fractures in oil and gas reservoirs have been the topic of many studies and have attracted reservoir research all over the world. Because of the complexities of the fractures, it is difficult to use fractured reservoir core samples to investigate true underground conditions. Due to the diversity of the fracture parameters, the simulation and evaluation of fractured rock in the laboratory setting is also difficult. Previous researchers have typically used a single material, such as resin, to simulate fractures. There has been a great deal of simplifying of the materials and conditions, which has led to disappointing results in application. In the present study, sandstone core samples were selected and sectioned to simulate fractures, and the changes of the compressional and shear waves were measured with the gradual increasing of the fracture width. The effects of the simulated fracture width on the acoustic wave velocity and amplitude were analyzed. Two variables were defined: H represents the amplitude attenuation ratio of the compressional and shear wave, and x represents the transit time difference value of the shear wave and compressional wave divided by the transit time of the compressional wave. The effect of fracture width on these two physical quantities was then analyzed. Finally, the methods of quantitative evaluation for fracture width with H and x were obtained. The experimental results showed that the rock fractures linearly reduced the velocity of the shear and compressional waves. The effect of twin fractures on thecompressional velocity was almost equal to that of a single fracture which had the same fracture width as the sum of the twin fractures. At the same time, the existence of fractures led to acoustic wave amplitude attenuations, and the compressional wave attenuation was two times greater than that of the shear wave. In this paper, a method was proposed to calculate the fracture width with x and H, then this was applied to the array acoustic imaging logging data.The application examples showed that the calculated fracture width could be compared with fractures on the electric imaging logs. These rules were applied in the well logs to effectively evaluate the fractures, under the case of no image logs, which had significance to prospecting and development of oil and gas in fractured reservoirs.展开更多
This work analyzes the quality of crustal tilt and strain observations during 2014, which were acquired from 269 sets of ground tiltmeters and 212 sets of strainmeters. In terms of data quality, the water tube tiltmet...This work analyzes the quality of crustal tilt and strain observations during 2014, which were acquired from 269 sets of ground tiltmeters and 212 sets of strainmeters. In terms of data quality, the water tube tiltmeters presented the highest rate of excellent quality,approximately 91%, and the pendulum tiltmeters and ground strainmeters yielded rates of81% and 78%, respectively. This means that a total of 380 sets of instruments produced high-quality observational data suitable for scientific investigations and analyses.展开更多
A simple approach to calculate the amplitude component of a wave front propagating in space from a hologram is proposed. It is able to calculate the amplitude distribution on a plane at any distance rapidly using a st...A simple approach to calculate the amplitude component of a wave front propagating in space from a hologram is proposed. It is able to calculate the amplitude distribution on a plane at any distance rapidly using a standard GPU. This is useful for analyzing and reconstructing the 3D image encoded on a hologram.展开更多
An unsteady load calculation method for the support configuration of a monopile-supported offshore wind turbine is developed based on the Fluent software platform.Firstly,the water wave is generated by imposing the in...An unsteady load calculation method for the support configuration of a monopile-supported offshore wind turbine is developed based on the Fluent software platform.Firstly,the water wave is generated by imposing the inlet boundary conditions according to the exact potential flow solution.Then the wave evolution is simulated by solving the unsteady incompressible Navier-Stokes(N-S)equations coupled with the volume of fluid method.For the small amplitude wave with reasonable wave parameters,the numerical wave result agrees well with that of the given wave model.Finally,a monopile support configuration is introduced and a CFD-based load calculation method is established to accurately calculate the unsteady load under the combined action of wave and wind.The computed unsteady wave load on a small-size monopile support located in the small amplitude wave flow coincides with that of the Morison formula.The load calculations are also performed on a large-size monopile support and a monopile-supported offshore wind turbine under the combined action of small amplitude wave and wind.展开更多
In this paper,several new constant-amplitude and variable-amplitude wave solutions(namely,traveling wave solutions) of a generalized nonlinear Schrdinger equation are investigated by using the extended homogeneous b...In this paper,several new constant-amplitude and variable-amplitude wave solutions(namely,traveling wave solutions) of a generalized nonlinear Schrdinger equation are investigated by using the extended homogeneous balance method,where the balance method is applied to solve the Riccati equation and the reduced nonlinear ordinary differential equation,respectively.In addition,stability analysis of those solutions are also conducted by regular phase plane technique.展开更多
基金Supported by the Knowledge Innovation Program of Chinese Academy of Sciences (No. KZCX-YW-12-04)the National Natural Science Foundation of China (No. 41030855)+1 种基金the National High Technology Research and Development Program of China (863 Program) (No. 2008AA09A403)the Marine Public Welfare Project of China (No. 201105032)
文摘One of the most important parameters for oceanic internal waves (IWs) is their amplitude. We have developed a method to retrieve the IW amplitude from nautical X-Band radar images based on the KdV equation for continuous stratified finite depth system. We have also tested the method of measuring the amplitude of IWs from X-Band radar backscatter image sequences acquired on June 2009 in the northeastern South China Sea. The method was applied in several radar images. Experiments show that the retrieval amplitudes are consistent with the in-situ observational amplitudes of IWs by using the towed thermistor chain and conductivity-temperature-depth (CTD) profile. The uncertainty of the method is also discussed.
基金financially supported by the Chinese State Key Laboratory of Ocean Engineering(Grant No.GKZD010068/084).
文摘In order to forecast the distribution of crest amplitudes and the occurrence of freak waves in a short crested coastal sea,a novel transformed linear simulation method is initially proposed in this paper.A Hermite transformation model expressed as a monotonic cubic polynomial serves as the foundation for the novel simulation technique.The wave crest amplitude exceedance probabilities of two sea states-one with a directional wave spectrum based on the measured wave elevation data at the Yura coast and the other with a typical directional JONSWAP wave spectrum-have been predicted using the novel simulation method that has been proposed.The likelihood that a particular critical wave crest amplitude will be exceeded is directly correlated with the probability that freak waves will occur.It is shown that the novel simulation approach suggested can provide predictions that are more precise than those obtained from the Rayleigh crest amplitude distribution model,the Jahns and Wheeler crest amplitude distribution model,or the conventional linear simulation method.This study also demonstrated that the nonlinear simulation method is less effective than the novel simulation method in terms of efficiency.
文摘A new theoretical solution is presented here for the dynamic characteristics of a buoyant jet due to opposing small amplitude waves. The conservation equations of mass, tangential moment^n and vertical momentum are solved by the integral method which encompasses the Gaussian profiles of velocity and density. The action of waves is incorporated into the equations of motion as an external force and a new exact solution is obtained to predict the trajectory, velocity distribution and boundary thickness of the buoyant jet over an arbitrary lateral cross section. It is found that the velocity along the centerline is inversely proportional to the ratio of the momentum of the wave to the buoyant jet. The averaged bound- ary width varies with the fluctuation of the boundary width, the distance from the orifice and the velocity correction function. Owing to the motion of waves, the fluctuation of the boundary width is proportional to the wave steepness.
文摘The wave propagation in the one-dimensional complex Ginzbur-Landau equation (CGLE) is studied by considering a wave source at the system boundary. A special propagation region, which is an island-shaped zone surrounded by the defect turbulence in the system parameter space, is observed in our numerical experiment. The wave signal spreads in the whole space with a novel amplitude wave pattern in the area. The relevant factors of the pattern formation, such as the wave speed, the maximum propagating distance and the oscillatory frequency, are studied in detail. The stability and the generality of the region are testified by adopting various initial conditions. This finding of the amplitude pattern extends the wave propagation region in the parameter space and presents a new signal transmission mode, and is therefore expected to be of much importance.
文摘We theoretically analyze the organizing filament of small amplitude scroll waves in general excitable media by perturbation method and explicitly give the expressions of coefficients in Keener theory. In particular for the excitable media with equal diffusion, we obtain a close system for the motion of the filament. With an example of the Oregonator model, our results are in good agreement with those simulated by Winfree.
文摘A new method of detecting the vertical bearing capacity for single-pile with high strain is discussed in this paper. A heavy hammer or a small type of rocket is used to strike the pile top and the detectors are used to record vibra- tion graphs. An expression of higher degree of strain (deformation force) is introduced. It is testified theoretically that the displacement, velocity and acceleration cannot be obtained by simple integral acceleration and differential velocity when long displacement and high strain exist, namely when the pile phase generates a whole slip relative to the soil body. That is to say that there are non-linear relations between them. It is educed accordingly that the force P and displacement S are calculated from the amplitude of wave train and (dynamic) P-S curve is drew so as to determine the yield points. Further, a method of determining the vertical bearing capacity for single-pile is dis- cussed. A static load test is utilized to check the result of dynamic test and determine the correlative constants of dynamic-static P(Q)-S curve.
基金Supported by the Department of Science and Technology,Government of India for their financial support of this work through the SERC Fast Track Scheme for Young Scientist(No.SR/FTP/MS-037/2011)
文摘In this paper we study the problem of generation of surface waves produced due to a) rolling of the plate and b) presence of a line source in front of a fixed vertical plate. The amplitudes of radiated waves at large distance from the plate, in both cases, are obtained by a suitable application of Green's integral theorem. These are then studied graphically for various values of the ice cover parameter.
文摘This paper gives an overall discussion about water level change on slopes under wave action, including wave runup, wave rundown and wave up-down amplitude, and a suggested formula for their calculation.
文摘The waves driven by the wind do not move on the water as ordinarily done by sailboats. Indeed, the movement of the waves driven by the wind is more complex than the sailboats’ translation movement that we know. The movement of the wave in our particular case results from the chain-job done by wind’s stress and gravity forces: material is collected upstream (erosion phenomenon) and then deposited on the wave’s summit by the wind. This material deposited on the summit of the wave by the wind is then removed and dispatched on the downstream side of the wave by gravity forces. As always happens in any chain-job: if the wind works faster than gravity forces, great accumulation of material will occur at the summit of the wave that will lead to an increase in its (the wave in this case) height. If conversely the wind works more slowly, a deficit in material delivery will occur and gravity force goes directly to remove material on the wave’s summit and lead to a decrease in its height. In terms of Mechanics, we know that the main obstacle that can seriously disturb the work of the wind is the unavailability of water or so its viscosity. Given the complexity of the process to be studied, it seemed necessary for us to make a use of modulational instability theories such as the standard NLSE in order to better understand the contribution of wind and water viscosity to modulations of driven waves’ amplitudes (or phases): modulations which sometimes can accidentally trigger unpredictable rogue waves.
基金financially supported by the Scientific Research Fund of Heilongjiang Provincial Education Department(Grant No.12541132)the Natural Science Youth Foundation of Heilongjiang Province of China(Grant No.QC2015082)
文摘Based on the Burgers equation and Manley-Rowe equation, the derivation about nonlinear interaction of the acoustic waves has been done in this paper. After nonlinear interaction among the low-frequency weak waves and the pump wave, the analytical solutions of acoustic waves' amplitude in the field are deduced. The relationship between normalized energy of high-frequency and the change of acoustic energy before and after the nonlinear interaction of the acoustic waves is analyzed. The experimental results about the changes of the acoustic energy are presented. The study shows that new frequencies are generated and the energies of the low-frequency are modulated in a long term by the pump waves, which leads the energies of the low-frequency acoustic waves to change in the pulse trend in the process of the nonlinear interaction of the acoustic waves. The increase and decrease of the energies of the low-frequency are observed under certain typical conditions, which lays a foundation for practical engineering applications.
基金supported by the Natural Science Foundation of Inner Mongolia,China(Grant No 200711020116)Open Fund of the Key Laboratory of Ocean Circulation and Waves,Chinese Academy of Sciences(Grant No KLOCAW0805)+1 种基金the Key Program of the Scientific Research Plan of Inner Mongolia University of Technology,China(Grant No ZD200608)National Science Fund for Distinguished Young Scholars of China(Grant No 40425015)
文摘Interracial internal waves in a three-layer density-stratified fluid are investigated using a singular perturbation method, and third-order asymptotic solutions of the velocity potentials and third-order Stokes wave solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory. As expected, the third-order solutions describe the third-order nonlinear modification and the third-order nonlinear interactions between the interracial waves. The wave velocity depends on not only the wave number and the depth of each layer but also on the wave amplitude.
基金supported in part by the National Natural Science Foundation of China (Grant No. 41174096)the Graduate Innovation Fund of Jilin University (Project No. 2016103)
文摘Fractures in oil and gas reservoirs have been the topic of many studies and have attracted reservoir research all over the world. Because of the complexities of the fractures, it is difficult to use fractured reservoir core samples to investigate true underground conditions. Due to the diversity of the fracture parameters, the simulation and evaluation of fractured rock in the laboratory setting is also difficult. Previous researchers have typically used a single material, such as resin, to simulate fractures. There has been a great deal of simplifying of the materials and conditions, which has led to disappointing results in application. In the present study, sandstone core samples were selected and sectioned to simulate fractures, and the changes of the compressional and shear waves were measured with the gradual increasing of the fracture width. The effects of the simulated fracture width on the acoustic wave velocity and amplitude were analyzed. Two variables were defined: H represents the amplitude attenuation ratio of the compressional and shear wave, and x represents the transit time difference value of the shear wave and compressional wave divided by the transit time of the compressional wave. The effect of fracture width on these two physical quantities was then analyzed. Finally, the methods of quantitative evaluation for fracture width with H and x were obtained. The experimental results showed that the rock fractures linearly reduced the velocity of the shear and compressional waves. The effect of twin fractures on thecompressional velocity was almost equal to that of a single fracture which had the same fracture width as the sum of the twin fractures. At the same time, the existence of fractures led to acoustic wave amplitude attenuations, and the compressional wave attenuation was two times greater than that of the shear wave. In this paper, a method was proposed to calculate the fracture width with x and H, then this was applied to the array acoustic imaging logging data.The application examples showed that the calculated fracture width could be compared with fractures on the electric imaging logs. These rules were applied in the well logs to effectively evaluate the fractures, under the case of no image logs, which had significance to prospecting and development of oil and gas in fractured reservoirs.
基金supported by Special Foundation of Earthquake Science(201408006)Director Foundation of Institute of Seismology,China Earthquake Administration(201516214)
文摘This work analyzes the quality of crustal tilt and strain observations during 2014, which were acquired from 269 sets of ground tiltmeters and 212 sets of strainmeters. In terms of data quality, the water tube tiltmeters presented the highest rate of excellent quality,approximately 91%, and the pendulum tiltmeters and ground strainmeters yielded rates of81% and 78%, respectively. This means that a total of 380 sets of instruments produced high-quality observational data suitable for scientific investigations and analyses.
基金the UK Engineering and Physical Sciences Research Council (EPSRC) for the support through the Platform Grant in Liquid Crystal Photonics
文摘A simple approach to calculate the amplitude component of a wave front propagating in space from a hologram is proposed. It is able to calculate the amplitude distribution on a plane at any distance rapidly using a standard GPU. This is useful for analyzing and reconstructing the 3D image encoded on a hologram.
基金supported partly by the National Basic Research Program of China("973"Program)(No.2014CB046200)the National Natural Science Foundation of China(No.11372135)the NUAA Fundamental Research Funds(No.NS2013005)
文摘An unsteady load calculation method for the support configuration of a monopile-supported offshore wind turbine is developed based on the Fluent software platform.Firstly,the water wave is generated by imposing the inlet boundary conditions according to the exact potential flow solution.Then the wave evolution is simulated by solving the unsteady incompressible Navier-Stokes(N-S)equations coupled with the volume of fluid method.For the small amplitude wave with reasonable wave parameters,the numerical wave result agrees well with that of the given wave model.Finally,a monopile support configuration is introduced and a CFD-based load calculation method is established to accurately calculate the unsteady load under the combined action of wave and wind.The computed unsteady wave load on a small-size monopile support located in the small amplitude wave flow coincides with that of the Morison formula.The load calculations are also performed on a large-size monopile support and a monopile-supported offshore wind turbine under the combined action of small amplitude wave and wind.
基金supported by the National NSF of China(11571088)NSF of Zhejiang Province(LY13A010020)Program(HNUEYT2013)
文摘In this paper,several new constant-amplitude and variable-amplitude wave solutions(namely,traveling wave solutions) of a generalized nonlinear Schrdinger equation are investigated by using the extended homogeneous balance method,where the balance method is applied to solve the Riccati equation and the reduced nonlinear ordinary differential equation,respectively.In addition,stability analysis of those solutions are also conducted by regular phase plane technique.