During transient electromagnetic method (TEM) exploration of a copper mine, we detected the late-channel TEM signal reversal phenomenon (a voltage change from positive to negative) caused by the influence of the i...During transient electromagnetic method (TEM) exploration of a copper mine, we detected the late-channel TEM signal reversal phenomenon (a voltage change from positive to negative) caused by the influence of the induced polarization (IP) effect, which affects the depth and precision of the TEM detection. The conventional inversion method is inefficient because it is difficult to process the data. In this paper, the Cole-Cole model is adopted to analyze the effect of Dc resistivity, chargeability, time constant, and frequency exponent on the TEM response in an homogeneous half space model. Singular Value Decomposition (SVD) is used to invert the measured TEM data, and the Dc resistivity, chargeability, time constant and frequency exponent were extracted from the measured TEM data in the mine area. The extracted parameters are used for interpreting the detection result as a supplement. This reveals why the TEM data acquired in the area has a low resolution. It was found that the DC resistivity and time constant do not significantly change the results, however, the chargeability and frequency exponent have a significant effect. Because of these influences, the SVD method is more accurate than the conventional method in the apparent resistivity profile. The area of the copper mine is confined accurately based on the SVD inverted data. The conclusion has been verified by drill and is identical to the practical geological situation.展开更多
An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge...An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge. In this study, a comparison between the electromagnetic fields of an infinitesimal point charge and a dipole has been presented. First, the time-domain potential function of a point source in a 3D conductive medium is derived. Then the electric and magnetic fields in a 3D homogeneous lossless space are derived via the relation between the potential and field. The field differences between the infinitesimal point charge and the dipole in the step-off time, far-source, and near-source zones are analyzed, and the accuracy of the solutions from these sources is investigated. It is also shown that the field of the infinitesimal point charge in the near-source zone is different from that of the dipole, whereas the far-source zone fields of these two sources are identical. The comparison of real and simulated data shows that the infinitesimal point charge represents the real source better than the divole source.展开更多
We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular...We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.展开更多
We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method...We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized: the real number of the electromagnetic field solution based on the Gaver-Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the "ground excitation-stratum measurement method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drill- hole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research.展开更多
A finite element method with boundary element method (FEM-BEM) is presented for computing electromagnetic induction. The features of an edge element method including the volume and surface edge element method are inve...A finite element method with boundary element method (FEM-BEM) is presented for computing electromagnetic induction. The features of an edge element method including the volume and surface edge element method are investigated in depth. Surface basis functions of edge elements to an arbitrary shape of target are derived according to the geometrical property of basis functions and applied to discretize the surface integral equation for 3-D general targets. The proposed model is presented to compute resonant frequencies and surface current of underground unexplored ordnance (UXO), and then the electromagnetic responses of single target with different frequencies and positions of sensor are simulated and results are validated by experiments.展开更多
Tunneling machines, or excavators, are large and good conductors and affect the reliability of data gathering and interpretation in advanced detection using transient electromagnetic methods. In our experiment, we use...Tunneling machines, or excavators, are large and good conductors and affect the reliability of data gathering and interpretation in advanced detection using transient electromagnetic methods. In our experiment, we used a coincident-loop and central loop type of configuration, where the coil plane l) vertical to and 2) parallel to the working face. A SIROTEM instrument at different locations was used to observe the transient electromagnetic responses of the excavator and to analyze the response amplitudes. The result shows that the tunneling machine affects the advanced detection data and is related to the way the coil is coupled. When the excavator is 6 m from the observatory, the interference of tunneling machine can be ignored.展开更多
The main design factors which effect on operating speed of solenoid actuator for valve operation are mass of plunger, electromagnetic motive force, inductance and return spring, and these factors are not independent b...The main design factors which effect on operating speed of solenoid actuator for valve operation are mass of plunger, electromagnetic motive force, inductance and return spring, and these factors are not independent but related with each other in view point of design and electromagnetic theory. It is impossible to increase the operating speed by only change the value of any one design factor. The change of any one value results in change of any value related it in various design factors. Permanent magnets are as assistant materials which make higher flux density in air gap. Electromagnetic motive force in controlled only by current. This paper presents a speed increasing design method of solenoid actuator using a solenoid, by some governing equations which are composed of electromagnetic theory and empirical knowledge, and proved the propriety by experiments.展开更多
Solenoid valve is one of the executive parts of weft insertion control system. According to the response characteristics of the solenoid valve, an improved design becomes a necessity. Firstly, the numerical model was ...Solenoid valve is one of the executive parts of weft insertion control system. According to the response characteristics of the solenoid valve, an improved design becomes a necessity. Firstly, the numerical model was established after analyzing the solenoid valve during its start-up and shut-down. Comparing the simulation data with the practical data, it is verified that the numerical simulation model has a high feasibility. Secondly, excitation voltage and spring pre-compression were adjusted respectively, and the response rules after adjusting were investigated. The research of the study shows: the response time tends to be inverse proportional to the excitation voltage during start-up, and it becomes a constant value with the increase of the excitation voltage; the response time is proportional to the spring pre-compression when the solenoid valve starts up, it is inverse proportional to spring pre-compression when the solenoid valve shuts down. And the total response time is a constant value with the increase of the spring pre-compression. Therefore, the value of the excitation voltage and the spring pre-compression should be selected when the curve is becoming flatten. The results of the research can provide the reference to the further development of the solenoid valve.展开更多
In the present paper,from the second order partial differential equations for solving the magnetotelluric(MT) fields of general anisotropic medium,we first obtained the second order partial differential equations for ...In the present paper,from the second order partial differential equations for solving the magnetotelluric(MT) fields of general anisotropic medium,we first obtained the second order partial differential equations for some anisotropic media with special conductivity(e.g.diagonal anisotropy,transverse anisotropy,azimuthal anisotropy,etc.) by simplifying the electrical conductivity tensor of anisotropic medium.And then we obtained the analytic solutions to MT fields for the case of transverse and azimuthal anisotropy through converting the conductivity parameter based on that of diagonal anisotropy.We further discussed the influence of the selection of integral limit and step length on precision in solving the analytic solutions for MT fields of isotropic medium.Finally,we presented the MT responses of two transverse and azimuthal anisotropic media as well as some applications of the analytic solutions to MT fields of anisotropic medium.展开更多
Flexible magnetoelectric (ME) materials have been studied for new applications such as memory, energy harvesters, and magnetic field sensors. Herein, with the widely studied and progressive advantages of ME phenomen...Flexible magnetoelectric (ME) materials have been studied for new applications such as memory, energy harvesters, and magnetic field sensors. Herein, with the widely studied and progressive advantages of ME phenomena in the multiferroic field, we demonstrate a new approach for utilizing flexible ME materials as gate dielectric layers in ME organic field-effect transistors (ME-OFET) that can be used for sensing a magnetic field and extracting the ME properties of the gate dielectric itself. The magnetoelectric nanohybrid gate dielectric layer comprises sandwiched stacks of magnetostrictive CoFe2O4 nanoparticles and a highly piezoelectric poly(vinylidene fluoride-co-trifluoroethylene) layer. While varying the magnetic field applied to the ME gate dielectric, the ME effect in the functional gate dielectric modulates the channel conductance of the ME-OFET owing to a change in the effective gate field. The clear separation of the ME responses in the gate dielectric layer of ME-OFET from those of the other parameters was demonstrated using the AC gate biasing method and enabled the extraction of the ME coefficient of ME materials. Additionally, the device shows high stability after cyclic bending of 10,000 cycles at a banding radius of 1.2 cm. The device has significant potential for not only the extraction of the intrinsic characterization of ME materials but also the sensing of a magnetic field in integrated flexible electronic systems.展开更多
Concrete forms of resonant response (ER) for a strong electromagnetic (EM) wave beam (photon flux) propagating in a static magnetic field to a standing gravitational wave (gravitons) are given, and the corresponding p...Concrete forms of resonant response (ER) for a strong electromagnetic (EM) wave beam (photon flux) propagating in a static magnetic field to a standing gravitational wave (gravitons) are given, and the corresponding perturbation solutions and resonant conditions are obtained. It is found that perturbed EM fields (PEMFs) contain three new components with frequencies Io,* w,l and ωPg respectively. In the case of ωe?ωg, the PEMFs are manifested as the EM wave beams with frequency ωe and a standing EM wave with ωg. The former and the background EM wave beam (BE-MWB) have the same propagating direction, while in the case of ωg?ωe, all PEMFs are expressed as the standing EM waves with frequency ωg. The resonant response occurs in two cases of ωe = 1/2 ωg andωe, = ωg only. Then not only the first order perturbed energy fluxes (PEFs) propagating in the same and opposite directions of the BEMWB can be generated, but also radial and tangential PEFs which are perpendicular to the above directions can be produced. This effect might provide a new way for the EM detection of the gravitational waves (GWs). Moreover, the possible schemes of displaying perturbed effects induced by the standing GW withh = 10-33 - 10-35 and λg = 0.1 m at the level of the single photon avalanche and in a typicla laboratory dimension are reviewed.展开更多
Conventional attractive magnetic force models (proportional to the coil current squared and inversely proportional to the gap squared) cannot simulate the nonlinear responses of magnetic bearings in the presence of el...Conventional attractive magnetic force models (proportional to the coil current squared and inversely proportional to the gap squared) cannot simulate the nonlinear responses of magnetic bearings in the presence of electromagnetic losses,flux leakage or saturation of iron.In this paper,based on results from an experimental set-up designed to study magnetic force,a novel parametric model is presented in the form of a nonlinear polynomial with unknown coefficients.The parameters of the proposed model are identified using the weighted residual method.Validations of the model identified were performed by comparing the results in time and frequency domains.The results show a good correlation between experiments and numerical simulations.展开更多
基金supported by the National Technology R&D Program in the 11th Five year Plan of China(No.2007BAQ00168-1-1)the National Natural Science Foundation of China(No. 41103052/D0309)the Shanxi Province Excellent Graduate Innovation Program(No. 20113038)
文摘During transient electromagnetic method (TEM) exploration of a copper mine, we detected the late-channel TEM signal reversal phenomenon (a voltage change from positive to negative) caused by the influence of the induced polarization (IP) effect, which affects the depth and precision of the TEM detection. The conventional inversion method is inefficient because it is difficult to process the data. In this paper, the Cole-Cole model is adopted to analyze the effect of Dc resistivity, chargeability, time constant, and frequency exponent on the TEM response in an homogeneous half space model. Singular Value Decomposition (SVD) is used to invert the measured TEM data, and the Dc resistivity, chargeability, time constant and frequency exponent were extracted from the measured TEM data in the mine area. The extracted parameters are used for interpreting the detection result as a supplement. This reveals why the TEM data acquired in the area has a low resolution. It was found that the DC resistivity and time constant do not significantly change the results, however, the chargeability and frequency exponent have a significant effect. Because of these influences, the SVD method is more accurate than the conventional method in the apparent resistivity profile. The area of the copper mine is confined accurately based on the SVD inverted data. The conclusion has been verified by drill and is identical to the practical geological situation.
基金supported by Chinese National Programs for Fundamental Research and Development(No.2012CB416605)the National Natural Science Foundation of China(No.41174090)Development Project of National Key Scientific Equipment(No.ZDYZ2012-1-05-04)
文摘An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge. In this study, a comparison between the electromagnetic fields of an infinitesimal point charge and a dipole has been presented. First, the time-domain potential function of a point source in a 3D conductive medium is derived. Then the electric and magnetic fields in a 3D homogeneous lossless space are derived via the relation between the potential and field. The field differences between the infinitesimal point charge and the dipole in the step-off time, far-source, and near-source zones are analyzed, and the accuracy of the solutions from these sources is investigated. It is also shown that the field of the infinitesimal point charge in the near-source zone is different from that of the dipole, whereas the far-source zone fields of these two sources are identical. The comparison of real and simulated data shows that the infinitesimal point charge represents the real source better than the divole source.
基金supported by the Young Scientists Fund of the National Natural Science Foundation of China(No.41304082)the China Postdoctoral Science Foundation(No.2016M590731)+2 种基金the Young Scientists Fund of the Natural Science Foundation of Hebei Province(No.D2014403011)the Program for Young Excellent Talents of Higher Education Institutions of Hebei Province(No.BJ2016046)the Geological survey project of China Geological Survey(No.1212011121197)
文摘We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.
基金supported by the Young Scientists Fund of the National Natural Science Foundation of China(No.41304082)the Young Scientists Fund of the Natural Science Foundation of Hebei Province(No.D2014403011)the Geological survey project of China Geological Survey(No.12120114090201)
文摘We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized: the real number of the electromagnetic field solution based on the Gaver-Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the "ground excitation-stratum measurement method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drill- hole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research.
文摘A finite element method with boundary element method (FEM-BEM) is presented for computing electromagnetic induction. The features of an edge element method including the volume and surface edge element method are investigated in depth. Surface basis functions of edge elements to an arbitrary shape of target are derived according to the geometrical property of basis functions and applied to discretize the surface integral equation for 3-D general targets. The proposed model is presented to compute resonant frequencies and surface current of underground unexplored ordnance (UXO), and then the electromagnetic responses of single target with different frequencies and positions of sensor are simulated and results are validated by experiments.
基金support received from the National Basic Research Program of China (No2007CB209400)the National Natural Science Foundation of China (No50774085)the Young Scientists Fund of the School Science Foundation of CUMT (No2008A046)
文摘Tunneling machines, or excavators, are large and good conductors and affect the reliability of data gathering and interpretation in advanced detection using transient electromagnetic methods. In our experiment, we used a coincident-loop and central loop type of configuration, where the coil plane l) vertical to and 2) parallel to the working face. A SIROTEM instrument at different locations was used to observe the transient electromagnetic responses of the excavator and to analyze the response amplitudes. The result shows that the tunneling machine affects the advanced detection data and is related to the way the coil is coupled. When the excavator is 6 m from the observatory, the interference of tunneling machine can be ignored.
文摘The main design factors which effect on operating speed of solenoid actuator for valve operation are mass of plunger, electromagnetic motive force, inductance and return spring, and these factors are not independent but related with each other in view point of design and electromagnetic theory. It is impossible to increase the operating speed by only change the value of any one design factor. The change of any one value results in change of any value related it in various design factors. Permanent magnets are as assistant materials which make higher flux density in air gap. Electromagnetic motive force in controlled only by current. This paper presents a speed increasing design method of solenoid actuator using a solenoid, by some governing equations which are composed of electromagnetic theory and empirical knowledge, and proved the propriety by experiments.
基金supported by the National Natural Science Foundation of China(No.51006090)
文摘Solenoid valve is one of the executive parts of weft insertion control system. According to the response characteristics of the solenoid valve, an improved design becomes a necessity. Firstly, the numerical model was established after analyzing the solenoid valve during its start-up and shut-down. Comparing the simulation data with the practical data, it is verified that the numerical simulation model has a high feasibility. Secondly, excitation voltage and spring pre-compression were adjusted respectively, and the response rules after adjusting were investigated. The research of the study shows: the response time tends to be inverse proportional to the excitation voltage during start-up, and it becomes a constant value with the increase of the excitation voltage; the response time is proportional to the spring pre-compression when the solenoid valve starts up, it is inverse proportional to spring pre-compression when the solenoid valve shuts down. And the total response time is a constant value with the increase of the spring pre-compression. Therefore, the value of the excitation voltage and the spring pre-compression should be selected when the curve is becoming flatten. The results of the research can provide the reference to the further development of the solenoid valve.
基金supported by the National Natural Science Foundation of China(Grant No.40774035)
文摘In the present paper,from the second order partial differential equations for solving the magnetotelluric(MT) fields of general anisotropic medium,we first obtained the second order partial differential equations for some anisotropic media with special conductivity(e.g.diagonal anisotropy,transverse anisotropy,azimuthal anisotropy,etc.) by simplifying the electrical conductivity tensor of anisotropic medium.And then we obtained the analytic solutions to MT fields for the case of transverse and azimuthal anisotropy through converting the conductivity parameter based on that of diagonal anisotropy.We further discussed the influence of the selection of integral limit and step length on precision in solving the analytic solutions for MT fields of isotropic medium.Finally,we presented the MT responses of two transverse and azimuthal anisotropic media as well as some applications of the analytic solutions to MT fields of anisotropic medium.
文摘Flexible magnetoelectric (ME) materials have been studied for new applications such as memory, energy harvesters, and magnetic field sensors. Herein, with the widely studied and progressive advantages of ME phenomena in the multiferroic field, we demonstrate a new approach for utilizing flexible ME materials as gate dielectric layers in ME organic field-effect transistors (ME-OFET) that can be used for sensing a magnetic field and extracting the ME properties of the gate dielectric itself. The magnetoelectric nanohybrid gate dielectric layer comprises sandwiched stacks of magnetostrictive CoFe2O4 nanoparticles and a highly piezoelectric poly(vinylidene fluoride-co-trifluoroethylene) layer. While varying the magnetic field applied to the ME gate dielectric, the ME effect in the functional gate dielectric modulates the channel conductance of the ME-OFET owing to a change in the effective gate field. The clear separation of the ME responses in the gate dielectric layer of ME-OFET from those of the other parameters was demonstrated using the AC gate biasing method and enabled the extraction of the ME coefficient of ME materials. Additionally, the device shows high stability after cyclic bending of 10,000 cycles at a banding radius of 1.2 cm. The device has significant potential for not only the extraction of the intrinsic characterization of ME materials but also the sensing of a magnetic field in integrated flexible electronic systems.
基金We thank Profs. V. N. Rudenko and M. V. Sazhin for their helpful discussions and suggestions, colleagues of Sternberg Astronomical Institute of Moscow University for their hospitality, and Prof. Luo Jun of Huazhong University of Science and Technology fo
文摘Concrete forms of resonant response (ER) for a strong electromagnetic (EM) wave beam (photon flux) propagating in a static magnetic field to a standing gravitational wave (gravitons) are given, and the corresponding perturbation solutions and resonant conditions are obtained. It is found that perturbed EM fields (PEMFs) contain three new components with frequencies Io,* w,l and ωPg respectively. In the case of ωe?ωg, the PEMFs are manifested as the EM wave beams with frequency ωe and a standing EM wave with ωg. The former and the background EM wave beam (BE-MWB) have the same propagating direction, while in the case of ωg?ωe, all PEMFs are expressed as the standing EM waves with frequency ωg. The resonant response occurs in two cases of ωe = 1/2 ωg andωe, = ωg only. Then not only the first order perturbed energy fluxes (PEFs) propagating in the same and opposite directions of the BEMWB can be generated, but also radial and tangential PEFs which are perpendicular to the above directions can be produced. This effect might provide a new way for the EM detection of the gravitational waves (GWs). Moreover, the possible schemes of displaying perturbed effects induced by the standing GW withh = 10-33 - 10-35 and λg = 0.1 m at the level of the single photon avalanche and in a typicla laboratory dimension are reviewed.
文摘Conventional attractive magnetic force models (proportional to the coil current squared and inversely proportional to the gap squared) cannot simulate the nonlinear responses of magnetic bearings in the presence of electromagnetic losses,flux leakage or saturation of iron.In this paper,based on results from an experimental set-up designed to study magnetic force,a novel parametric model is presented in the form of a nonlinear polynomial with unknown coefficients.The parameters of the proposed model are identified using the weighted residual method.Validations of the model identified were performed by comparing the results in time and frequency domains.The results show a good correlation between experiments and numerical simulations.