The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the...The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the existing numerical schemes for it show some deficiencies. Based on the original hyperbolic mild-slope equation, a nonlinear dispersion relation is introduced in present paper to effectively take the nonlinear effect of waves into account and a new numerical scheme is proposed. The weakly nonlinear dispersion relation and the improved numerical scheme are applied to the simulation of wave transformation over an elliptic shoal. Numerical tests show that the improvement of the numerical scheme makes efficient the solution to the hyperbolic mild-slope equation, A comparison of numerical results with experimental data indicates that the results obtained by use of the new scheme are satisfactory.展开更多
The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However, its computational method in Cartesian coordinates often renders the model inaccur...The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However, its computational method in Cartesian coordinates often renders the model inaccurate in areas with irregular shorelines, such as estuaries and harbors. Based on the hyperbolic mild-slope equation in Cartesian coordinates, the numerical model in orthogonal curvilinear coordinates is developed. The transformed model is discretized by the finite difference method and solved by the ADI method with space-staggered grids. The numerical predictions in curvilinear co- ordinates show good agreemenl with the data obtained in three typical physical expedments, which demonstrates that the present model can be used to simulate wave propagation, for normal incidence and oblique incidence, in domains with complicated topography and boundary conditions.展开更多
Based on the hyperbolic mild-slope equations derived by Copeland (1985), a numerical model is established in unstag- gered grids. A composite 4 th-order Adam-Bashforth-Moulton (ABM) scheme is used to solve the model i...Based on the hyperbolic mild-slope equations derived by Copeland (1985), a numerical model is established in unstag- gered grids. A composite 4 th-order Adam-Bashforth-Moulton (ABM) scheme is used to solve the model in the time domain. Terms involving the first order spatial derivatives are differenced to O ( Δx )4accuracy utilizing a five-point formula. The nonlinear dispersion relationship proposed by Kirby and Dalrymple (1986) is used to consider the nonlinear effect. A numerical test is performed upon wave propagating over a typical shoal. The agreement between the numerical and the experimental results validates the present model. Biodistribution and applications are also summarized.展开更多
We are concerned with the large-time behavior of 3D quasilinear hyperbolic equations with nonlinear damping.The main novelty of this paper is two-fold.First,we prove the optimal decay rates of the second and third ord...We are concerned with the large-time behavior of 3D quasilinear hyperbolic equations with nonlinear damping.The main novelty of this paper is two-fold.First,we prove the optimal decay rates of the second and third order spatial derivatives of the solution,which are the same as those of the heat equation,and in particular,are faster than ones of previous related works.Second,for well-chosen initial data,we also show that the lower optimal L^(2) convergence rate of the k(∈[0,3])-order spatial derivatives of the solution is(1+t)^(-(2+2k)/4).Therefore,our decay rates are optimal in this sense.The proofs are based on the Fourier splitting method,low-frequency and high-frequency decomposition,and delicate energy estimates.展开更多
In this paper,we present a novel spatial reconstruction scheme,called AENO,that results from a special averaging of the ENO polynomial and its closest neighbour,while retaining the stencil direction decided by the ENO...In this paper,we present a novel spatial reconstruction scheme,called AENO,that results from a special averaging of the ENO polynomial and its closest neighbour,while retaining the stencil direction decided by the ENO choice.A variant of the scheme,called m-AENO,results from averaging the modified ENO(m-ENO)polynomial and its closest neighbour.The concept is thoroughly assessed for the one-dimensional linear advection equation and for a one-dimensional non-linear hyperbolic system,in conjunction with the fully discrete,high-order ADER approach implemented up to fifth order of accuracy in both space and time.The results,as compared to the conventional ENO,m-ENO and WENO schemes,are very encouraging.Surprisingly,our results show that the L_(1)-errors of the novel AENO approach are the smallest for most cases considered.Crucially,for a chosen error size,AENO turns out to be the most efficient method of all five methods tested.展开更多
High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of th...High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.展开更多
In this paper, we provide an explicit expression for the full Dirichlet-to-Neumann map corresponding to a radial potential for a hyperbolic differential equation in 3-dimensional. We show that the Dirichlet-Neumann op...In this paper, we provide an explicit expression for the full Dirichlet-to-Neumann map corresponding to a radial potential for a hyperbolic differential equation in 3-dimensional. We show that the Dirichlet-Neumann operators corresponding to a potential radial have the same properties for hyperbolic differential equations as for elliptic differential equations. We numerically implement the coefficients of the explicit formulas. Moreover, a Lipschitz type stability is established near the edge of the domain by an estimation constant. That is necessary for the reconstruction of the potential from Dirichlet-to-Neumann map in the inverse problem for a hyperbolic differential equation.展开更多
We consider a strongly non-linear degenerate parabolic-hyperbolic problem with p(x)-Laplacian diffusion flux function. We propose an entropy formulation and prove the existence of an entropy solution.
In this paper,the oscillation of solutions of hyperbolic partial functional differential equations is studied,and oscillatory criteria of solutions with three kinds of boundary conditions are obtained.
An efficient numerical model for wave refraction, diffraction and reflection is presented in this paper. In the model the modified time-dependent mild-slope equation is transformed into an evolution equation and an im...An efficient numerical model for wave refraction, diffraction and reflection is presented in this paper. In the model the modified time-dependent mild-slope equation is transformed into an evolution equation and an improved ADI method involving a relaxation factor is adopted to solve it. The method has the advantage of improving the numerical stability and convergence rate by properly determining the relaxation factor. The range of the relaxation factor making the differential scheme unconditionally stable is determined by stability analysis. Several verifications are performed to examine the accuracy of the present model. The numerical results coincide with the analytic solutions or experimental data very well and the computer time is reduced.展开更多
A new form of hyperbolic mild slope equations is derived with the inclusion of the amphtude dispersion of nonlinear waves. The effects of including the amplitude dispersion effect on the wave propagation are discussed...A new form of hyperbolic mild slope equations is derived with the inclusion of the amphtude dispersion of nonlinear waves. The effects of including the amplitude dispersion effect on the wave propagation are discussed. Wave breaking mechanism is incorporated into the present model to apply the new equations to surf zone. The equations are solved nu- merically for regular wave propagation over a shoal and in surf zone, and a comparison is made against measurements. It is found that the inclusion of the amplitude dispersion can also improve model' s performance on prediction of wave heights around breaking point for the wave motions in surf zone.展开更多
In this paper, oscillatory properties for solutions of the systems of certain quasilinear impulsive delay hyperbolic equations with nonlinear diffusion coefficient are investigated. A sufficient criterion for oscillat...In this paper, oscillatory properties for solutions of the systems of certain quasilinear impulsive delay hyperbolic equations with nonlinear diffusion coefficient are investigated. A sufficient criterion for oscillations of such systems is obtained.展开更多
Sufficient conditions are obtained for the oscillation of solutions of the systems of quasilinear hyperbolic differential equation with deviating arguments under nonlinear boundary condition.
New hyperbolic mild slope equations for random waves are developed with the inclusion of amplitude dispersion. The frequency perturbation around the peak frequency of random waves is adopted to extend the equations fo...New hyperbolic mild slope equations for random waves are developed with the inclusion of amplitude dispersion. The frequency perturbation around the peak frequency of random waves is adopted to extend the equations for regular waves to random waves. The nonlinear effect of amplitude dispersion is incorporated approximately into the model by only considering the nonlinear effect on the carrier waves of random waves, which is done by introducing a representative wave amplitude for the carrier waves. The computation time is gready saved by the introduction of the representative wave amplitude. The extension of the present model to breaking waves is also considered in order to apply the new equations to surf zone. The model is validated for random waves propagate over a shoal and in surf zone against measurements.展开更多
A proper orthogonal decomposition(POD) method was successfully used in the reduced-order modeling of complex systems.In this paper,we extend the applications of POD method,namely,apply POD method to a classical fini...A proper orthogonal decomposition(POD) method was successfully used in the reduced-order modeling of complex systems.In this paper,we extend the applications of POD method,namely,apply POD method to a classical finite element(FE) formulation for second-order hyperbolic equations with real practical applied background,establish a reduced FE formulation with lower dimensions and high enough accuracy,and provide the error estimates between the reduced FE solutions and the classical FE solutions and the implementation of algorithm for solving reduced FE formulation so as to provide scientific theoretic basis for service applications.Some numerical examples illustrate the fact that the results of numerical computation are consistent with theoretical conclusions.Moreover,it is shown that the reduced FE formulation based on POD method is feasible and efficient for solving FE formulation for second-order hyperbolic equations.展开更多
Nonlinear effect is of importance to waves propagating from deep water to shallow water. The non-linearity of waves is widely discussed due to its high precision in application. But there are still some problems in de...Nonlinear effect is of importance to waves propagating from deep water to shallow water. The non-linearity of waves is widely discussed due to its high precision in application. But there are still some problems in dealing with the nonlinear waves in practice. In this paper, a modified form of mild-slope equation with weakly nonlinear effect is derived by use of the nonlinear dispersion relation and the steady mild-slope equation containing energy dissipation. The modified form of mild-slope equation is convenient to solve nonlinear effect of waves. The model is tested against the laboratory measurement for the case of a submerged elliptical shoal on a slope beach given by Berkhoff et al. The present numerical results are also compared with those obtained through linear wave theory. Better agreement is obtained as the modified mild-slope equation is employed. And the modified mild-slope equation can reasonably simulate the weakly nonlinear effect of wave propagation from deep water to coast.展开更多
In this paper, oscillatory properties of all solutions for neutral type impulsive hyperbolic equations with several delays under the Robin boundary condition are investigated and several new sufficient conditions for ...In this paper, oscillatory properties of all solutions for neutral type impulsive hyperbolic equations with several delays under the Robin boundary condition are investigated and several new sufficient conditions for oscillation are presented.展开更多
Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with...Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with curved boundary or irregular coastline. The proposed wave-induced current model includes a nearshore current module established through orthogonal curvilinear transformation form of shallow water equations and a wave module based on the curvilinear parabolic approximation wave equation. The wave module actually serves as the driving force to provide the current module with required radiation stresses. The Crank-Nicolson finite difference scheme and the alternating directions implicit method are used to solve the wave and current module, respectively. The established surf zone currents model is validated by two numerical experiments about longshore currents and rip currents in basins with rip channel and breakwater. The numerical results are compared with the measured data and published numerical results.展开更多
In this article, a proper orthogonal decomposition (POD) method is used to study a classical splitting positive definite mixed finite element (SPDMFE) formulation for second- order hyperbolic equations. A POD redu...In this article, a proper orthogonal decomposition (POD) method is used to study a classical splitting positive definite mixed finite element (SPDMFE) formulation for second- order hyperbolic equations. A POD reduced-order SPDMFE extrapolating algorithm with lower dimensions and sufficiently high accuracy is established for second-order hyperbolic equations. The error estimates between the classical SPDMFE solutions and the reduced-order SPDMFE solutions obtained from the POD reduced-order SPDMFE extrapolating algorithm are provided. The implementation for solving the POD reduced-order SPDMFE extrapolating algorithm is given. Some numerical experiments are presented illustrating that the results of numerical computation are consistent with theoretical conclusions, thus validating that the POD reduced-order SPDMFE extrapolating algorithm is feasible and efficient for solving second-order hyperbolic equations.展开更多
基金This subject was financially supported by the National Natural Science Foundation of China(Grant No. 59839330 and No.59976047) by the Visiting Scholal Foundation of State Key Hydraulic Lab.of High Speed Flows of Dalian University of Technology.
文摘The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the existing numerical schemes for it show some deficiencies. Based on the original hyperbolic mild-slope equation, a nonlinear dispersion relation is introduced in present paper to effectively take the nonlinear effect of waves into account and a new numerical scheme is proposed. The weakly nonlinear dispersion relation and the improved numerical scheme are applied to the simulation of wave transformation over an elliptic shoal. Numerical tests show that the improvement of the numerical scheme makes efficient the solution to the hyperbolic mild-slope equation, A comparison of numerical results with experimental data indicates that the results obtained by use of the new scheme are satisfactory.
基金supported by the National Basic Research Program of China ( Grant No.2006CB403302)the National Natural Science Foundation of China (Grant Nos .50839001 and 50709004)the Scientific Research Foundation of the Higher Education Institutions of Liaoning Province (Grant No.2006T018)
文摘The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However, its computational method in Cartesian coordinates often renders the model inaccurate in areas with irregular shorelines, such as estuaries and harbors. Based on the hyperbolic mild-slope equation in Cartesian coordinates, the numerical model in orthogonal curvilinear coordinates is developed. The transformed model is discretized by the finite difference method and solved by the ADI method with space-staggered grids. The numerical predictions in curvilinear co- ordinates show good agreemenl with the data obtained in three typical physical expedments, which demonstrates that the present model can be used to simulate wave propagation, for normal incidence and oblique incidence, in domains with complicated topography and boundary conditions.
基金This work is financially supported by the National Natural Science Foundation of China (50409015).
文摘Based on the hyperbolic mild-slope equations derived by Copeland (1985), a numerical model is established in unstag- gered grids. A composite 4 th-order Adam-Bashforth-Moulton (ABM) scheme is used to solve the model in the time domain. Terms involving the first order spatial derivatives are differenced to O ( Δx )4accuracy utilizing a five-point formula. The nonlinear dispersion relationship proposed by Kirby and Dalrymple (1986) is used to consider the nonlinear effect. A numerical test is performed upon wave propagating over a typical shoal. The agreement between the numerical and the experimental results validates the present model. Biodistribution and applications are also summarized.
基金partially supported by the National Nature Science Foundation of China(12271114)the Guangxi Natural Science Foundation(2023JJD110009,2019JJG110003,2019AC20214)+2 种基金the Innovation Project of Guangxi Graduate Education(JGY2023061)the Key Laboratory of Mathematical Model and Application(Guangxi Normal University)the Education Department of Guangxi Zhuang Autonomous Region。
文摘We are concerned with the large-time behavior of 3D quasilinear hyperbolic equations with nonlinear damping.The main novelty of this paper is two-fold.First,we prove the optimal decay rates of the second and third order spatial derivatives of the solution,which are the same as those of the heat equation,and in particular,are faster than ones of previous related works.Second,for well-chosen initial data,we also show that the lower optimal L^(2) convergence rate of the k(∈[0,3])-order spatial derivatives of the solution is(1+t)^(-(2+2k)/4).Therefore,our decay rates are optimal in this sense.The proofs are based on the Fourier splitting method,low-frequency and high-frequency decomposition,and delicate energy estimates.
文摘In this paper,we present a novel spatial reconstruction scheme,called AENO,that results from a special averaging of the ENO polynomial and its closest neighbour,while retaining the stencil direction decided by the ENO choice.A variant of the scheme,called m-AENO,results from averaging the modified ENO(m-ENO)polynomial and its closest neighbour.The concept is thoroughly assessed for the one-dimensional linear advection equation and for a one-dimensional non-linear hyperbolic system,in conjunction with the fully discrete,high-order ADER approach implemented up to fifth order of accuracy in both space and time.The results,as compared to the conventional ENO,m-ENO and WENO schemes,are very encouraging.Surprisingly,our results show that the L_(1)-errors of the novel AENO approach are the smallest for most cases considered.Crucially,for a chosen error size,AENO turns out to be the most efficient method of all five methods tested.
文摘High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.
文摘In this paper, we provide an explicit expression for the full Dirichlet-to-Neumann map corresponding to a radial potential for a hyperbolic differential equation in 3-dimensional. We show that the Dirichlet-Neumann operators corresponding to a potential radial have the same properties for hyperbolic differential equations as for elliptic differential equations. We numerically implement the coefficients of the explicit formulas. Moreover, a Lipschitz type stability is established near the edge of the domain by an estimation constant. That is necessary for the reconstruction of the potential from Dirichlet-to-Neumann map in the inverse problem for a hyperbolic differential equation.
文摘We consider a strongly non-linear degenerate parabolic-hyperbolic problem with p(x)-Laplacian diffusion flux function. We propose an entropy formulation and prove the existence of an entropy solution.
文摘In this paper,the oscillation of solutions of hyperbolic partial functional differential equations is studied,and oscillatory criteria of solutions with three kinds of boundary conditions are obtained.
文摘An efficient numerical model for wave refraction, diffraction and reflection is presented in this paper. In the model the modified time-dependent mild-slope equation is transformed into an evolution equation and an improved ADI method involving a relaxation factor is adopted to solve it. The method has the advantage of improving the numerical stability and convergence rate by properly determining the relaxation factor. The range of the relaxation factor making the differential scheme unconditionally stable is determined by stability analysis. Several verifications are performed to examine the accuracy of the present model. The numerical results coincide with the analytic solutions or experimental data very well and the computer time is reduced.
基金the National Natural Science Foundation of China (Grant Nos .50479053 and10672034)the Programfor Changjiang Scholars and Innovative Research Teamin Universitythe foundation for doctoral degree education of the Education Ministry of China
文摘A new form of hyperbolic mild slope equations is derived with the inclusion of the amphtude dispersion of nonlinear waves. The effects of including the amplitude dispersion effect on the wave propagation are discussed. Wave breaking mechanism is incorporated into the present model to apply the new equations to surf zone. The equations are solved nu- merically for regular wave propagation over a shoal and in surf zone, and a comparison is made against measurements. It is found that the inclusion of the amplitude dispersion can also improve model' s performance on prediction of wave heights around breaking point for the wave motions in surf zone.
基金Supported by Hunan Provincial NSF(05jj400008)of China under Grant.
文摘In this paper, oscillatory properties for solutions of the systems of certain quasilinear impulsive delay hyperbolic equations with nonlinear diffusion coefficient are investigated. A sufficient criterion for oscillations of such systems is obtained.
基金This work is supported in part by NNSF of China(10571126)and in part by Program for New Century Excellent Talents in University.
文摘Sufficient conditions are obtained for the oscillation of solutions of the systems of quasilinear hyperbolic differential equation with deviating arguments under nonlinear boundary condition.
基金supported by the National Natural Science Foundation of China(Grant Nos.50479053and10672034)the Program for Changjiang Scholars and Innovative Research Teamin University,and thefoundationfordoctoral degree education of the Education Ministry of China
文摘New hyperbolic mild slope equations for random waves are developed with the inclusion of amplitude dispersion. The frequency perturbation around the peak frequency of random waves is adopted to extend the equations for regular waves to random waves. The nonlinear effect of amplitude dispersion is incorporated approximately into the model by only considering the nonlinear effect on the carrier waves of random waves, which is done by introducing a representative wave amplitude for the carrier waves. The computation time is gready saved by the introduction of the representative wave amplitude. The extension of the present model to breaking waves is also considered in order to apply the new equations to surf zone. The model is validated for random waves propagate over a shoal and in surf zone against measurements.
基金supported by the National Science Foundation of China (11061009,40821092)the National Basic Research Program (2010CB428403,2009CB421407,2010CB951001)Natural Science Foundation of Hebei Province (A2010001663)
文摘A proper orthogonal decomposition(POD) method was successfully used in the reduced-order modeling of complex systems.In this paper,we extend the applications of POD method,namely,apply POD method to a classical finite element(FE) formulation for second-order hyperbolic equations with real practical applied background,establish a reduced FE formulation with lower dimensions and high enough accuracy,and provide the error estimates between the reduced FE solutions and the classical FE solutions and the implementation of algorithm for solving reduced FE formulation so as to provide scientific theoretic basis for service applications.Some numerical examples illustrate the fact that the results of numerical computation are consistent with theoretical conclusions.Moreover,it is shown that the reduced FE formulation based on POD method is feasible and efficient for solving FE formulation for second-order hyperbolic equations.
文摘Nonlinear effect is of importance to waves propagating from deep water to shallow water. The non-linearity of waves is widely discussed due to its high precision in application. But there are still some problems in dealing with the nonlinear waves in practice. In this paper, a modified form of mild-slope equation with weakly nonlinear effect is derived by use of the nonlinear dispersion relation and the steady mild-slope equation containing energy dissipation. The modified form of mild-slope equation is convenient to solve nonlinear effect of waves. The model is tested against the laboratory measurement for the case of a submerged elliptical shoal on a slope beach given by Berkhoff et al. The present numerical results are also compared with those obtained through linear wave theory. Better agreement is obtained as the modified mild-slope equation is employed. And the modified mild-slope equation can reasonably simulate the weakly nonlinear effect of wave propagation from deep water to coast.
文摘In this paper, oscillatory properties of all solutions for neutral type impulsive hyperbolic equations with several delays under the Robin boundary condition are investigated and several new sufficient conditions for oscillation are presented.
基金supported by the National Natural Science Foundation of China (Grant Nos. 50839001 and 50979036)
文摘Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with curved boundary or irregular coastline. The proposed wave-induced current model includes a nearshore current module established through orthogonal curvilinear transformation form of shallow water equations and a wave module based on the curvilinear parabolic approximation wave equation. The wave module actually serves as the driving force to provide the current module with required radiation stresses. The Crank-Nicolson finite difference scheme and the alternating directions implicit method are used to solve the wave and current module, respectively. The established surf zone currents model is validated by two numerical experiments about longshore currents and rip currents in basins with rip channel and breakwater. The numerical results are compared with the measured data and published numerical results.
基金supported by the National Science Foundation of China(11271127,11361035)Science Research of Guizhou Education Department(QJHKYZ[2013]207)Natural Science Foundation of Inner Mongolia(2012MS0106)
文摘In this article, a proper orthogonal decomposition (POD) method is used to study a classical splitting positive definite mixed finite element (SPDMFE) formulation for second- order hyperbolic equations. A POD reduced-order SPDMFE extrapolating algorithm with lower dimensions and sufficiently high accuracy is established for second-order hyperbolic equations. The error estimates between the classical SPDMFE solutions and the reduced-order SPDMFE solutions obtained from the POD reduced-order SPDMFE extrapolating algorithm are provided. The implementation for solving the POD reduced-order SPDMFE extrapolating algorithm is given. Some numerical experiments are presented illustrating that the results of numerical computation are consistent with theoretical conclusions, thus validating that the POD reduced-order SPDMFE extrapolating algorithm is feasible and efficient for solving second-order hyperbolic equations.