In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hy...In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.展开更多
A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order...A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.展开更多
Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference ...Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference techniques. The existence and uniqueness of the weak solution are proved without any assumptions on choice of the spacetime meshes. Basic error estimates in L-infinity (L-2) norm, that is maximum-norm in time, L-2-norm in space are obtained. The numerical results are given in the last part and the analysis between theoretic and experimental results are obtained.展开更多
Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved thro...Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.展开更多
This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their...This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.展开更多
Finite element method (FEM) is an efficient numerical tool for the solution of partial differential equations (PDEs). It is one of the most general methods when compared to other numerical techniques. PDEs posed in a ...Finite element method (FEM) is an efficient numerical tool for the solution of partial differential equations (PDEs). It is one of the most general methods when compared to other numerical techniques. PDEs posed in a variational form over a given space, say a Hilbert space, are better numerically handled with the FEM. The FEM algorithm is used in various applications which includes fluid flow, heat transfer, acoustics, structural mechanics and dynamics, electric and magnetic field, etc. Thus, in this paper, the Finite Element Orthogonal Collocation Approach (FEOCA) is established for the approximate solution of Time Fractional Telegraph Equation (TFTE) with Mamadu-Njoseh polynomials as grid points corresponding to new basis functions constructed in the finite element space. The FEOCA is an elegant mixture of the Finite Element Method (FEM) and the Orthogonal Collocation Method (OCM). Two numerical examples are experimented on to verify the accuracy and rate of convergence of the method as compared with the theoretical results, and other methods in literature.展开更多
Deformation behavior of slab at the straightening stage during continuous casting was simulated by the explicit dynamic finite element method,and the stress distribution along the width direction of the slab and its c...Deformation behavior of slab at the straightening stage during continuous casting was simulated by the explicit dynamic finite element method,and the stress distribution along the width direction of the slab and its change regularity at slab center during continuous casting were obtained.The influence of distribution and change of stress on the propagation of longitudinal cracks on slab surface was discussed.The results show that the tensional stress appears on slab surface at the inner arc side and the compressive stress appears on slab surface at the outer arc side at stages 6-8 in straightening zone during continuous casting.Longitudinal cracks generally appear on slab top surface and do not appear on slab bottom surface,which are also observed in industry.展开更多
A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretizatio...A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming (P1)2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.展开更多
On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these meth...On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these methods, the transverse displacement is approximated by conforming (bi)linear macroelements or (bi)quadratic elements, and the rotation by conforming (bi)linear elements. The shear stress can be locally computed from transverse displacement and rotation. Uniform in plate thickness, optimal error bounds are obtained for the transverse displacement, rotation, and shear stress in their natural norms. Numerical results are presented to illustrate the theoretical results.展开更多
The two-dimensional problem of generalized thermoelastic diffusion material with thermal and diffusion relaxation times is investigated in the context of the Lord-Shulman theory. As an application of the problem, a pa...The two-dimensional problem of generalized thermoelastic diffusion material with thermal and diffusion relaxation times is investigated in the context of the Lord-Shulman theory. As an application of the problem, a particular type of thermal source is considered and the problem is solved numerically by using a finite element method. The components of displacement, stress, temperature distribution chemical potential and mass concentration are obtained. The resulting quantities are depicted graphically for a special model. An appreciable effect of relaxation times is observed on various resulting quantities.展开更多
The wireline formation tester (WFT) is an important tool for formation evaluation, such as calculating the formation pressure and permeability, identifying the fluid type, and determining the interface between oil a...The wireline formation tester (WFT) is an important tool for formation evaluation, such as calculating the formation pressure and permeability, identifying the fluid type, and determining the interface between oil and water. However, in a low porosity and low permeability formation, the supercharge pressure effect exists, since the mudcake has a poor sealing ability. The mudcake cannot isolate the hydrostatic pressure of the formation around the borehole and the mud seeps into the formations, leading to inaccurate formation pressure measurement. At the same time, the tool can be easily stuck in the low porosity/low permeability formation due to the long waiting and testing time. We present a method for determining the minimum testing time for the wireline formation tester. The pressure distribution of the mudcake and the formation were respectively calculated with the finite element method (FEM). The radius of the influence of mud pressure was also computed, and the minimum testing time in low porosity/low permeability formations was determined within a range of values for different formation permeabilities. The determination of the minimum testing time ensures an accurate formation pressure measurement and minimizes possible accidents due to long waiting and testing time.展开更多
By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved hav...By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved having third-order pseudo- symplectic scheme respectively for general Hamiltonian systems, and they both keep energy conservative. The finite element methods are proved to be symplectic as well as energy conservative for linear Hamiltonian systems. The numerical results are in agree-ment with theory.展开更多
In this article, we consider a two-dimensional symmetric space-fractional diffusion equation in which the space fractional derivatives are defined in Riesz potential sense. The well-posed feature is guaranteed by ener...In this article, we consider a two-dimensional symmetric space-fractional diffusion equation in which the space fractional derivatives are defined in Riesz potential sense. The well-posed feature is guaranteed by energy inequality. To solve the diffusion equation, a fully discrete form is established by employing Crank-Nicolson technique in time and Galerkin finite element method in space. The stability and convergence are proved and the stiffness matrix is given analytically. Three numerical examples are given to confirm our theoretical analysis in which we find that even with the same initial condition, the classical and fractional diffusion equations perform differently but tend to be uniform diffusion at last.展开更多
A space-time finite element method,discontinuous in time but continuous in space,is studied to solve the nonlinear forward-backward heat equation.A linearized technique is introduced in order to obtain the error estim...A space-time finite element method,discontinuous in time but continuous in space,is studied to solve the nonlinear forward-backward heat equation.A linearized technique is introduced in order to obtain the error estimates of the approximate solutions.And the numerical simulations are given.展开更多
Interdigitated finger capacitance of a continuous-wave terahertz photomixer is calculated using the finite element method.For the frequently used electrode width(0.2 μm) and gap width(1.8 μm),the finger capacita...Interdigitated finger capacitance of a continuous-wave terahertz photomixer is calculated using the finite element method.For the frequently used electrode width(0.2 μm) and gap width(1.8 μm),the finger capacitance increases quasi-quadratically with the number of electrodes increasing.The quasi-quadratic dependence can be explained by a sequence of lumped capacitors connected in parallel.For a photomixer composed of 10 electrodes and 9 photoconductive gaps,the finger capacitance increases as the gap width increases at a small electrode width,and follows the reverse trend at a large electrode width.For a constant electrode width,the finger capacitance first decreases and then slightly increases as the gap broadens until the smallest finger capacitance is formed.We also investigate the finger capacitances at different electrode and gap configurations with the 8 μm × 8 μm photomixer commonly used in previous studies.These calculations lead to a better understanding of the finger capacitance affected by the finger parameters,and should lead to terahertz photomixer optimization.展开更多
In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complic...In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.展开更多
We’ll consider the model of two-phase compressible miscible displacement in porous media which includes molecular diffusion and dispersion in one dimensional space. Time-discretization procedure is established and an...We’ll consider the model of two-phase compressible miscible displacement in porous media which includes molecular diffusion and dispersion in one dimensional space. Time-discretization procedure is established and analysed. The optimal error estimate in L2 norm is proved by introducing a new interpolation operator R.展开更多
In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wa...In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wave number domain and with a 2.5D approach.Three-dimensional numerical models formulated in the time/space domain are less frequently used,mainly due to their high computational cost.Notwithstanding,these models present very attractive characteristics,such as the possibility of considering nonlinear behaviors or the modelling of excess pore pressure and non-homogeneous and non-periodic geometries in the longitudinal direction of the track.In this study,two 3D numerical approaches formulated in the time/space domain are compared and experimentally validated.The first one consists of a finite element approach and the second one of a finite difference approach.The experimental validation in an actual case situated in Carregado(Portugal)shows an acceptable fitting between the numerical results and the actual measurements for both models.However,there are some differences among them.This study therefore includes some recommendations for their use in practical soil dynamics and geotechnical engineering.展开更多
文摘In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.
基金supported by the National Natural Science Foundation of China (No. 10601022)NSF ofInner Mongolia Autonomous Region of China (No. 200607010106)513 and Science Fund of InnerMongolia University for Distinguished Young Scholars (No. ND0702)
文摘A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.
文摘Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference techniques. The existence and uniqueness of the weak solution are proved without any assumptions on choice of the spacetime meshes. Basic error estimates in L-infinity (L-2) norm, that is maximum-norm in time, L-2-norm in space are obtained. The numerical results are given in the last part and the analysis between theoretic and experimental results are obtained.
基金Project supported by the National Basic Research Program of China (973 program) (No.G1999032804)
文摘Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.
文摘This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.
文摘Finite element method (FEM) is an efficient numerical tool for the solution of partial differential equations (PDEs). It is one of the most general methods when compared to other numerical techniques. PDEs posed in a variational form over a given space, say a Hilbert space, are better numerically handled with the FEM. The FEM algorithm is used in various applications which includes fluid flow, heat transfer, acoustics, structural mechanics and dynamics, electric and magnetic field, etc. Thus, in this paper, the Finite Element Orthogonal Collocation Approach (FEOCA) is established for the approximate solution of Time Fractional Telegraph Equation (TFTE) with Mamadu-Njoseh polynomials as grid points corresponding to new basis functions constructed in the finite element space. The FEOCA is an elegant mixture of the Finite Element Method (FEM) and the Orthogonal Collocation Method (OCM). Two numerical examples are experimented on to verify the accuracy and rate of convergence of the method as compared with the theoretical results, and other methods in literature.
基金Project(50634030) supported by the National Natural Science Foundation of ChinaProject(20090042120005) supported by the Doctorate Foundation of the Ministry of Education of ChinaProject(2006CB605208-1) supported by the State Basic Research Program of China
文摘Deformation behavior of slab at the straightening stage during continuous casting was simulated by the explicit dynamic finite element method,and the stress distribution along the width direction of the slab and its change regularity at slab center during continuous casting were obtained.The influence of distribution and change of stress on the propagation of longitudinal cracks on slab surface was discussed.The results show that the tensional stress appears on slab surface at the inner arc side and the compressive stress appears on slab surface at the outer arc side at stages 6-8 in straightening zone during continuous casting.Longitudinal cracks generally appear on slab top surface and do not appear on slab bottom surface,which are also observed in industry.
基金supported by the National Natural Science Foundation of China(No.10771150)the National Basic Research Program of China(No.2005CB321701)+1 种基金the Program for New Century Excellent Talents in University(No.NCET-07-0584)the Natural Science Foundation of Sichuan Province(No.07ZB087)
文摘A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming (P1)2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.
基金Project supported by China Postdoctoral Science Foundation (20100481488), Key Fund Project of Advanced Research of the Weapon Equipment (9140A33040512JB3401).
基金supported by NSFC(11571266,91430106,11171168,11071132)NSFC-RGC(China-Hong Kong)(11661161017)
文摘On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these methods, the transverse displacement is approximated by conforming (bi)linear macroelements or (bi)quadratic elements, and the rotation by conforming (bi)linear elements. The shear stress can be locally computed from transverse displacement and rotation. Uniform in plate thickness, optimal error bounds are obtained for the transverse displacement, rotation, and shear stress in their natural norms. Numerical results are presented to illustrate the theoretical results.
文摘The two-dimensional problem of generalized thermoelastic diffusion material with thermal and diffusion relaxation times is investigated in the context of the Lord-Shulman theory. As an application of the problem, a particular type of thermal source is considered and the problem is solved numerically by using a finite element method. The components of displacement, stress, temperature distribution chemical potential and mass concentration are obtained. The resulting quantities are depicted graphically for a special model. An appreciable effect of relaxation times is observed on various resulting quantities.
文摘The wireline formation tester (WFT) is an important tool for formation evaluation, such as calculating the formation pressure and permeability, identifying the fluid type, and determining the interface between oil and water. However, in a low porosity and low permeability formation, the supercharge pressure effect exists, since the mudcake has a poor sealing ability. The mudcake cannot isolate the hydrostatic pressure of the formation around the borehole and the mud seeps into the formations, leading to inaccurate formation pressure measurement. At the same time, the tool can be easily stuck in the low porosity/low permeability formation due to the long waiting and testing time. We present a method for determining the minimum testing time for the wireline formation tester. The pressure distribution of the mudcake and the formation were respectively calculated with the finite element method (FEM). The radius of the influence of mud pressure was also computed, and the minimum testing time in low porosity/low permeability formations was determined within a range of values for different formation permeabilities. The determination of the minimum testing time ensures an accurate formation pressure measurement and minimizes possible accidents due to long waiting and testing time.
基金Project supported by the National Natural Science Foundation of China (No.10471038)
文摘By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved having third-order pseudo- symplectic scheme respectively for general Hamiltonian systems, and they both keep energy conservative. The finite element methods are proved to be symplectic as well as energy conservative for linear Hamiltonian systems. The numerical results are in agree-ment with theory.
文摘In this article, we consider a two-dimensional symmetric space-fractional diffusion equation in which the space fractional derivatives are defined in Riesz potential sense. The well-posed feature is guaranteed by energy inequality. To solve the diffusion equation, a fully discrete form is established by employing Crank-Nicolson technique in time and Galerkin finite element method in space. The stability and convergence are proved and the stiffness matrix is given analytically. Three numerical examples are given to confirm our theoretical analysis in which we find that even with the same initial condition, the classical and fractional diffusion equations perform differently but tend to be uniform diffusion at last.
文摘A space-time finite element method,discontinuous in time but continuous in space,is studied to solve the nonlinear forward-backward heat equation.A linearized technique is introduced in order to obtain the error estimates of the approximate solutions.And the numerical simulations are given.
基金Project supported by the National High Technology Research and Development Program of China (Grant No. 2011AAxxx2008A)Hundred Talent Program of the Chinese Academy of Sciences (Grant No. J08-029)the Main Direction Program of Knowledge Innovation of the Chinese Academy of Sciences (Grant No. YYYJ-1123-4)
文摘Interdigitated finger capacitance of a continuous-wave terahertz photomixer is calculated using the finite element method.For the frequently used electrode width(0.2 μm) and gap width(1.8 μm),the finger capacitance increases quasi-quadratically with the number of electrodes increasing.The quasi-quadratic dependence can be explained by a sequence of lumped capacitors connected in parallel.For a photomixer composed of 10 electrodes and 9 photoconductive gaps,the finger capacitance increases as the gap width increases at a small electrode width,and follows the reverse trend at a large electrode width.For a constant electrode width,the finger capacitance first decreases and then slightly increases as the gap broadens until the smallest finger capacitance is formed.We also investigate the finger capacitances at different electrode and gap configurations with the 8 μm × 8 μm photomixer commonly used in previous studies.These calculations lead to a better understanding of the finger capacitance affected by the finger parameters,and should lead to terahertz photomixer optimization.
基金National Natural Science Foundation of China Under Grant No. 50178065
文摘In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.
基金This work was supported by National Science Foundation and China State Major Key Project for Basic Research
文摘We’ll consider the model of two-phase compressible miscible displacement in porous media which includes molecular diffusion and dispersion in one dimensional space. Time-discretization procedure is established and analysed. The optimal error estimate in L2 norm is proved by introducing a new interpolation operator R.
文摘In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wave number domain and with a 2.5D approach.Three-dimensional numerical models formulated in the time/space domain are less frequently used,mainly due to their high computational cost.Notwithstanding,these models present very attractive characteristics,such as the possibility of considering nonlinear behaviors or the modelling of excess pore pressure and non-homogeneous and non-periodic geometries in the longitudinal direction of the track.In this study,two 3D numerical approaches formulated in the time/space domain are compared and experimentally validated.The first one consists of a finite element approach and the second one of a finite difference approach.The experimental validation in an actual case situated in Carregado(Portugal)shows an acceptable fitting between the numerical results and the actual measurements for both models.However,there are some differences among them.This study therefore includes some recommendations for their use in practical soil dynamics and geotechnical engineering.