In this paper,we establish global classical solutions of semilinear wave equations with small compact supported initial data posed on the product space R^(3)×T.The semilinear nonlinearity is assumed to be of the ...In this paper,we establish global classical solutions of semilinear wave equations with small compact supported initial data posed on the product space R^(3)×T.The semilinear nonlinearity is assumed to be of the cubic form.The main ingredient here is the establishment of the L^(2)-L^(∞)decay estimates and the energy estimates for the linear problem,which are adapted to the wave equation on the product space.The proof is based on the Fourier mode decomposition of the solution with respect to the periodic direction,the scaling technique,and the combination of the decay estimates and the energy estimates.展开更多
Research on seismic anisotropy and attenuation plays a significant role in exploration geophysics. To enhance the imaging quality for complicated structures, we develop several effective improvements for anisotropic a...Research on seismic anisotropy and attenuation plays a significant role in exploration geophysics. To enhance the imaging quality for complicated structures, we develop several effective improvements for anisotropic attenuation effects in reverse-time migration (Q-RTM) on surface and vertical seismic profiling (VSP) acquisition geometries. First, to suppress pseudo-shear wave artifact and numerical instability of the commonly used anisotropic pseudo-acoustic wave equations, an optimized pure P-wave dispersion relation is derived and the corresponding pure-mode wave equation is solved by combining the finite-difference and Possion methods. Second, a simplified anisotropic pure-mode visco-acoustic wave equation (PVAWE) based on standard linear solid model is established. Third, a time-dispersion correlation strategy is applied to improve the modeling accuracy. Fourth, we extend a target-oriented scheme to anisotropic attenuated modeling and imaging. Instead of the conventional wavefield modeling and RTM, the proposed approach can extract available wavefield information near the target regions and produce high imaging resolution for target structures. Last, both anisotropic surface and VSP Q-RTMs are executed by combining optimized PVAWE, time-dispersion correlation and target-oriented algorithm. Modeling examples demonstrate the advantages of our schemes. Moreover, our modified Q-compensated imaging workflow can be regarded as a supplement to the classical anisotropic RTM.展开更多
Two(3+1)-dimensional shallow water wave equations are studied by using residual symmetry and the consistent Riccati expansion(CRE) method. Through localization of residual symmetries, symmetry reduction solutions of t...Two(3+1)-dimensional shallow water wave equations are studied by using residual symmetry and the consistent Riccati expansion(CRE) method. Through localization of residual symmetries, symmetry reduction solutions of the two equations are obtained. The CRE method is applied to the two equations to obtain new B?cklund transformations from which a type of interesting interaction solution between solitons and periodic waves is generated.展开更多
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this...In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations: the first-order (velocity-stress) and the second- order (displacement-stress) separate elastic wave equa- tions, with the first-order (velocity-stress) and the second- order (displacement-stress) full (or mixed) elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.展开更多
The solutions to a linear wave equation can satisfy the principle of superposition, i.e., the linear superposition of two or more known solutions is still a solution of the linear wave equation. We show in this articl...The solutions to a linear wave equation can satisfy the principle of superposition, i.e., the linear superposition of two or more known solutions is still a solution of the linear wave equation. We show in this article that many nonlinear wave equations possess exact traveling wave solutions involving hyperbolic, triangle, and exponential functions, and the suitable linear combinations of these known solutions can also constitute linear superposition solutions to some nonlinear wave equations with special structural characteristics. The linear superposition solutions to the generalized KdV equation K(2,2,1), the Oliver water wave equation, and the k(n, n) equation are given. The structure characteristic of the nonlinear wave equations having linear superposition solutions is analyzed, and the reason why the solutions with the forms of hyperbolic, triangle, and exponential functions can form the linear superposition solutions is also discussed.展开更多
A complete discrimination system for the fourth order polynomial is given. As an application, we have reduced a (1+1)-dimensional dispersive long wave equation with general coefficients to an elementary integral fo...A complete discrimination system for the fourth order polynomial is given. As an application, we have reduced a (1+1)-dimensional dispersive long wave equation with general coefficients to an elementary integral form and obtained its all possible exact travelling wave solutions including rational function type solutions, solitary wave solutions, triangle function type periodic solutions and Jacobian elliptic functions double periodic solutions. This method can be also applied to many other similar problems.展开更多
Cylinder pressure fluctuation during combustion process of internal combustion engine is closely related to combustion noise and knock.The current studies are based on cylinder pressure test to obtain information on c...Cylinder pressure fluctuation during combustion process of internal combustion engine is closely related to combustion noise and knock.The current studies are based on cylinder pressure test to obtain information on combustion noise and knock,but there is little for simulation of combustion pressure fluctuation.Based on effects of combustion process in the combustion chamber on cylinder pressure by using wave equation,the mechanism of pressure fluctuation during combustion is researched with three-dimensional acoustic wave equation and flow field model of KIVA program.The cylinder pressure fluctuation curve,temperature field and acoustic field are obtained from the coupled calculation of the wave equation and KIVA program.The frequency spectrum analysis is taken with the cylinder pressure oscillation of cylinder pressure measured and calculated.The calculation result is consistent with the experimental result.This indicates that the cylinder pressure fluctuation can be correctly calculated with the wave equation.Analysis of calculation results of temperature field and acoustic field shows that sound field changes faster than flame propagates,and distribution of sound field is more complicated.Combustion pressure oscillation in the diesel engine is under highly unstable condition.This indicates that the combination of cylinder pressure fluctuation model and combustion model is an effective method to study the pressure oscillations and a new method to study the combustion noise and knock.展开更多
A complete classical symmetry classification and a nonclassical symmetry classification of a class of nonlinear wave equations are given with three arbitrary parameter functions. The obtained results show that such no...A complete classical symmetry classification and a nonclassical symmetry classification of a class of nonlinear wave equations are given with three arbitrary parameter functions. The obtained results show that such nonlinear wave equations admit richer classical and nonclassical symmetries, leading to the conservation laws and the reduction of the wave equations. Some exact solutions of the considered wave equations for particular cases are derived.展开更多
Based on splitting multi-symplectic structures, a new multi-symplectic scheme is proposed and applied to a nonlinear wave equation. The explicit multi-symplectic scheme of the nonlinear wave equation is obtained, and ...Based on splitting multi-symplectic structures, a new multi-symplectic scheme is proposed and applied to a nonlinear wave equation. The explicit multi-symplectic scheme of the nonlinear wave equation is obtained, and the corresponding multi-symplectic conservation property is proved. The backward error analysis shows that the explicit multi-symplectic scheme has good accuracy. The sine-Gordon equation and the Klein-Gordon equation are simulated by an explicit multi-symplectic scheme. The numerical results show that the new explicit multi-symplectic scheme can well simulate the solitary wave behaviors of the nonlinear wave equation and approximately preserve the relative energy error of the equation.展开更多
In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models...In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order (M= 1) and the fourth order (M = 2) recursion formulas.展开更多
The mixed covolume method for the regularized long wave equation is devel- oped and studied. By introducing a transfer operator γh, which maps the trial function space into the test function space, and combining the ...The mixed covolume method for the regularized long wave equation is devel- oped and studied. By introducing a transfer operator γh, which maps the trial function space into the test function space, and combining the mixed finite element with the finite volume method, the nonlinear and linear Euler fully discrete mixed covolume schemes are constructed, and the existence and uniqueness of the solutions are proved. The optimal error estimates for these schemes are obtained. Finally, a numerical example is provided to examine the efficiency of the proposed schemes.展开更多
We investigate (2+1)-dimensional generalized modified dispersive water wave (GMDWW) equation by utilizing the bifurcation theory of dynamical systems. We give the phase portraits and bifurcation analysis of the plane ...We investigate (2+1)-dimensional generalized modified dispersive water wave (GMDWW) equation by utilizing the bifurcation theory of dynamical systems. We give the phase portraits and bifurcation analysis of the plane system corresponding to the GMDWW equation. By using the special orbits in the phase portraits, we analyze the existence of the traveling wave solutions. When some parameter takes special values, we obtain abundant exact kink wave solutions, singular wave solutions, periodic wave solutions, periodic singular wave solutions, and solitary wave solutions for the GMDWW equation.展开更多
The spectral distribution exp( ), where {} are the eigenvalues of the negative Laplacian -△=- in the (x^1,x^2)-plane, is studied for a variety of domains, where -∞< t <∞ and i=(1/2)(-1) . The dependence of (...The spectral distribution exp( ), where {} are the eigenvalues of the negative Laplacian -△=- in the (x^1,x^2)-plane, is studied for a variety of domains, where -∞< t <∞ and i=(1/2)(-1) . The dependence of (t)on the connectivity of a domain and the boundary conditions are analyzed. Particular attention is given to a general bounded domain Ω in R^2 with a smooth boundary Ω, where a finite number of piecewise smooth Dirichlet, Neumann and Robin boundary conditions on the piecewise smooth parts Γj(j = 1,……,n) of Ω are considered such that Some geometrical properties of Ω(e.g., the area of Ω, the total lengths of the boundary, the curvature of its boundary, etc.) are determined, from the asymptotic expansions of (t) for |t| → 0.展开更多
This paper is concerned with the initial-boundary value problem for damped wave equations with a nonlinear convection term in the multi-dimensional half space R n + : u tt u + u t + divf (u) = 0, t 〉 0, x = (x...This paper is concerned with the initial-boundary value problem for damped wave equations with a nonlinear convection term in the multi-dimensional half space R n + : u tt u + u t + divf (u) = 0, t 〉 0, x = (x 1 , x ′ ) ∈ R n + := R + × R n 1 , u(0, x) = u 0 (x) → u + , as x 1 → + ∞ , u t (0, x) = u 1 (x), u(t, 0, x ′ ) = u b , x ′ = (x 2 , x 3 , ··· , x n ) ∈ R n 1 . (I) For the non-degenerate case f ′ 1 (u + ) 〈 0, it was shown in [10] that the above initialboundary value problem (I) admits a unique global solution u(t, x) which converges to the corresponding planar stationary wave φ(x 1 ) uniformly in x 1 ∈ R + as time tends to infinity provided that the initial perturbation and/or the strength of the stationary wave are sufficiently small. And in [10] Ueda, Nakamura, and Kawashima proved the algebraic decay estimates of the tangential derivatives of the solution u(t, x) for t → + ∞ by using the space-time weighted energy method initiated by Kawashima and Matsumura [5] and improved by Nishihkawa [7]. Moreover, by using the same weighted energy method, an additional algebraic convergence rate in the normal direction was obtained by assuming that the initial perturbation decays algebraically. We note, however, that the analysis in [10] relies heavily on the assumption that f ′ (u) 〈 0. The main purpose of this paper isdevoted to discussing the case of f ′ 1 (u b ) ≥ 0 and we show that similar results still hold for such a case. Our analysis is based on some delicate energy estimates.展开更多
The numerical solution for a type of quasilinear wave equation is studied.The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved.The...The numerical solution for a type of quasilinear wave equation is studied.The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved.The error of the difference solution is estimated.The theoretical results are controlled on a numerical example.展开更多
The spheroidal wave functions are found to have extensive applications in many branches of physics and mathematics. We use the perturbation method in supersymmetric quantum mechanics to obtain the analytic ground eige...The spheroidal wave functions are found to have extensive applications in many branches of physics and mathematics. We use the perturbation method in supersymmetric quantum mechanics to obtain the analytic ground eigenvalue and the ground eigenfunction of the angular spheroidal wave equation at low frequency in a series form. Using this approach, the numerical determinations of the ground eigenvalue and the ground eigenfunction for small complex frequencies are also obtained.展开更多
The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, w...The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, which uses an orthogonal function system with a weight function as the basis function. Compared with the conventional moving least-squares (MLS) approximation, the algebraic equation system in the IMLS approximation is not ill-conditioned, and can be solved directly without deriving the inverse matrix. Because there are fewer coefficients in the IMLS than in the MLS approximation, fewer nodes are selected in the IEFG method than in the element-free Galerkin method. Thus, the IEFG method has a higher computing speed. In the IEFG method, the Galerkin weak form is employed to obtain a dis- cretized system equation, and the penalty method is applied to impose the essential boundary condition. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the wave equations and the boundary-initial conditions depend on time, the scal- ing parameter, number of nodes and the time step length are considered for the convergence study.展开更多
A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly t...A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly to the second-order elastic wave equation.In view of this,based on the first-order CPML absorbing boundary condition,we propose a new CPML(NCPML)boundary which can be directly applied to the second-order wave equation.We first systematically extend the first-order CPML technique into second-order wave equations,neglecting the space-varying characteristics of the partial damping coefficient in the complex-frequency domain,avoiding the generation of convolution in the time domain.We then transform the technique back to the time domain through the inverse Fourier transform.Numerical simulation indicates that the space-varying characteristics of the attenuation factor have little influence on the absorption effect and increase the memory at the same time.A number of numerical examples show that the NCPML proposed in this study is effective in simulating elastic wave propagation,and this algorithm is more efficient and requires less memory allocation than the conventional PML absorbing boundary.展开更多
In this paper, we make use of the auxiliary equation and the expanded mapping methods to find the new exact periodic solutions for (2+1)-dimensional dispersive long wave equations in mathematical physics, which are...In this paper, we make use of the auxiliary equation and the expanded mapping methods to find the new exact periodic solutions for (2+1)-dimensional dispersive long wave equations in mathematical physics, which are expressed by Jacobi elliptic functions, and obtain some new solitary wave solutions (m → 1). This method can also be used to explore new periodic wave solutions for other nonlinear evolution equations.展开更多
文摘In this paper,we establish global classical solutions of semilinear wave equations with small compact supported initial data posed on the product space R^(3)×T.The semilinear nonlinearity is assumed to be of the cubic form.The main ingredient here is the establishment of the L^(2)-L^(∞)decay estimates and the energy estimates for the linear problem,which are adapted to the wave equation on the product space.The proof is based on the Fourier mode decomposition of the solution with respect to the periodic direction,the scaling technique,and the combination of the decay estimates and the energy estimates.
基金supported by the National Key R&D Program of China(2021YFA0716902)National Natural Science Foundation of China(42004119,42174156)+1 种基金the Fundamental Research Funds for the Central Universities,CHD(300102261306)the National Engineering Research Center of Offshore Oil and Gas Exploration,No.6 Courtyard,Taiyanggong South Street,Chaoyang District,Beijing,100028.
文摘Research on seismic anisotropy and attenuation plays a significant role in exploration geophysics. To enhance the imaging quality for complicated structures, we develop several effective improvements for anisotropic attenuation effects in reverse-time migration (Q-RTM) on surface and vertical seismic profiling (VSP) acquisition geometries. First, to suppress pseudo-shear wave artifact and numerical instability of the commonly used anisotropic pseudo-acoustic wave equations, an optimized pure P-wave dispersion relation is derived and the corresponding pure-mode wave equation is solved by combining the finite-difference and Possion methods. Second, a simplified anisotropic pure-mode visco-acoustic wave equation (PVAWE) based on standard linear solid model is established. Third, a time-dispersion correlation strategy is applied to improve the modeling accuracy. Fourth, we extend a target-oriented scheme to anisotropic attenuated modeling and imaging. Instead of the conventional wavefield modeling and RTM, the proposed approach can extract available wavefield information near the target regions and produce high imaging resolution for target structures. Last, both anisotropic surface and VSP Q-RTMs are executed by combining optimized PVAWE, time-dispersion correlation and target-oriented algorithm. Modeling examples demonstrate the advantages of our schemes. Moreover, our modified Q-compensated imaging workflow can be regarded as a supplement to the classical anisotropic RTM.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11975156 and 12175148)。
文摘Two(3+1)-dimensional shallow water wave equations are studied by using residual symmetry and the consistent Riccati expansion(CRE) method. Through localization of residual symmetries, symmetry reduction solutions of the two equations are obtained. The CRE method is applied to the two equations to obtain new B?cklund transformations from which a type of interesting interaction solution between solitons and periodic waves is generated.
基金partially supported by China National Major Science and Technology Project (Subproject No:2011ZX05024-001-03)
文摘In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations: the first-order (velocity-stress) and the second- order (displacement-stress) separate elastic wave equa- tions, with the first-order (velocity-stress) and the second- order (displacement-stress) full (or mixed) elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.
文摘The solutions to a linear wave equation can satisfy the principle of superposition, i.e., the linear superposition of two or more known solutions is still a solution of the linear wave equation. We show in this article that many nonlinear wave equations possess exact traveling wave solutions involving hyperbolic, triangle, and exponential functions, and the suitable linear combinations of these known solutions can also constitute linear superposition solutions to some nonlinear wave equations with special structural characteristics. The linear superposition solutions to the generalized KdV equation K(2,2,1), the Oliver water wave equation, and the k(n, n) equation are given. The structure characteristic of the nonlinear wave equations having linear superposition solutions is analyzed, and the reason why the solutions with the forms of hyperbolic, triangle, and exponential functions can form the linear superposition solutions is also discussed.
文摘A complete discrimination system for the fourth order polynomial is given. As an application, we have reduced a (1+1)-dimensional dispersive long wave equation with general coefficients to an elementary integral form and obtained its all possible exact travelling wave solutions including rational function type solutions, solitary wave solutions, triangle function type periodic solutions and Jacobian elliptic functions double periodic solutions. This method can be also applied to many other similar problems.
基金supported by National Natural Science Foundation of China (Grant No. 50706032 and No. 50476065)PhD Programs Foundation of Ministry of Education of China (Grant No. 20070056011)
文摘Cylinder pressure fluctuation during combustion process of internal combustion engine is closely related to combustion noise and knock.The current studies are based on cylinder pressure test to obtain information on combustion noise and knock,but there is little for simulation of combustion pressure fluctuation.Based on effects of combustion process in the combustion chamber on cylinder pressure by using wave equation,the mechanism of pressure fluctuation during combustion is researched with three-dimensional acoustic wave equation and flow field model of KIVA program.The cylinder pressure fluctuation curve,temperature field and acoustic field are obtained from the coupled calculation of the wave equation and KIVA program.The frequency spectrum analysis is taken with the cylinder pressure oscillation of cylinder pressure measured and calculated.The calculation result is consistent with the experimental result.This indicates that the cylinder pressure fluctuation can be correctly calculated with the wave equation.Analysis of calculation results of temperature field and acoustic field shows that sound field changes faster than flame propagates,and distribution of sound field is more complicated.Combustion pressure oscillation in the diesel engine is under highly unstable condition.This indicates that the combination of cylinder pressure fluctuation model and combustion model is an effective method to study the pressure oscillations and a new method to study the combustion noise and knock.
基金Project supported by the National Natural Science Foundation of China(Nos.11071159 and11301259)the Shanghai Key Projects(No.12510501700)+1 种基金the Scientific Research of College of Inner Mongolia(No.NJZZ14053)the Natural Science Foundation of Inner Mongolia(Nos.2013MS0105and 2014MS0114)
文摘A complete classical symmetry classification and a nonclassical symmetry classification of a class of nonlinear wave equations are given with three arbitrary parameter functions. The obtained results show that such nonlinear wave equations admit richer classical and nonclassical symmetries, leading to the conservation laws and the reduction of the wave equations. Some exact solutions of the considered wave equations for particular cases are derived.
基金Project supported by the National Natural Science Foundation of China(Nos.11161017,11071251,and 10871099)the National Basic Research Program of China(973 Program)(No.2007CB209603)+1 种基金the Natural Science Foundation of Hainan Province(No.110002)the Scientific Research Foun-dation of Hainan University(No.kyqd1053)
文摘Based on splitting multi-symplectic structures, a new multi-symplectic scheme is proposed and applied to a nonlinear wave equation. The explicit multi-symplectic scheme of the nonlinear wave equation is obtained, and the corresponding multi-symplectic conservation property is proved. The backward error analysis shows that the explicit multi-symplectic scheme has good accuracy. The sine-Gordon equation and the Klein-Gordon equation are simulated by an explicit multi-symplectic scheme. The numerical results show that the new explicit multi-symplectic scheme can well simulate the solitary wave behaviors of the nonlinear wave equation and approximately preserve the relative energy error of the equation.
基金National Basic Research Program of China Under Grant No. 2007CB714200National Natural Science Foundation of China Under Grant No. 90715038
文摘In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order (M= 1) and the fourth order (M = 2) recursion formulas.
基金supported by the National Natural Science Fundation of China (No. 11061021)the Science Research of Inner Mongolia Advanced Education (Nos. NJ10006, NJ10016, and NJZZ12011)the National Science Foundation of Inner Mongolia (Nos. 2011BS0102 and 2012MS0106)
文摘The mixed covolume method for the regularized long wave equation is devel- oped and studied. By introducing a transfer operator γh, which maps the trial function space into the test function space, and combining the mixed finite element with the finite volume method, the nonlinear and linear Euler fully discrete mixed covolume schemes are constructed, and the existence and uniqueness of the solutions are proved. The optimal error estimates for these schemes are obtained. Finally, a numerical example is provided to examine the efficiency of the proposed schemes.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11361069 and 11775146).
文摘We investigate (2+1)-dimensional generalized modified dispersive water wave (GMDWW) equation by utilizing the bifurcation theory of dynamical systems. We give the phase portraits and bifurcation analysis of the plane system corresponding to the GMDWW equation. By using the special orbits in the phase portraits, we analyze the existence of the traveling wave solutions. When some parameter takes special values, we obtain abundant exact kink wave solutions, singular wave solutions, periodic wave solutions, periodic singular wave solutions, and solitary wave solutions for the GMDWW equation.
文摘The spectral distribution exp( ), where {} are the eigenvalues of the negative Laplacian -△=- in the (x^1,x^2)-plane, is studied for a variety of domains, where -∞< t <∞ and i=(1/2)(-1) . The dependence of (t)on the connectivity of a domain and the boundary conditions are analyzed. Particular attention is given to a general bounded domain Ω in R^2 with a smooth boundary Ω, where a finite number of piecewise smooth Dirichlet, Neumann and Robin boundary conditions on the piecewise smooth parts Γj(j = 1,……,n) of Ω are considered such that Some geometrical properties of Ω(e.g., the area of Ω, the total lengths of the boundary, the curvature of its boundary, etc.) are determined, from the asymptotic expansions of (t) for |t| → 0.
基金The research of Fan Lili was supported by two grants from the National Natural Science Foundation of China (10871151 10925103)+1 种基金the research of Liu Hongxia was supported by National Natural Science Foundation of China (10871082)the research of Yin Hui was supported by National Natural Sciences Foundation of China (10901064)
文摘This paper is concerned with the initial-boundary value problem for damped wave equations with a nonlinear convection term in the multi-dimensional half space R n + : u tt u + u t + divf (u) = 0, t 〉 0, x = (x 1 , x ′ ) ∈ R n + := R + × R n 1 , u(0, x) = u 0 (x) → u + , as x 1 → + ∞ , u t (0, x) = u 1 (x), u(t, 0, x ′ ) = u b , x ′ = (x 2 , x 3 , ··· , x n ) ∈ R n 1 . (I) For the non-degenerate case f ′ 1 (u + ) 〈 0, it was shown in [10] that the above initialboundary value problem (I) admits a unique global solution u(t, x) which converges to the corresponding planar stationary wave φ(x 1 ) uniformly in x 1 ∈ R + as time tends to infinity provided that the initial perturbation and/or the strength of the stationary wave are sufficiently small. And in [10] Ueda, Nakamura, and Kawashima proved the algebraic decay estimates of the tangential derivatives of the solution u(t, x) for t → + ∞ by using the space-time weighted energy method initiated by Kawashima and Matsumura [5] and improved by Nishihkawa [7]. Moreover, by using the same weighted energy method, an additional algebraic convergence rate in the normal direction was obtained by assuming that the initial perturbation decays algebraically. We note, however, that the analysis in [10] relies heavily on the assumption that f ′ (u) 〈 0. The main purpose of this paper isdevoted to discussing the case of f ′ 1 (u b ) ≥ 0 and we show that similar results still hold for such a case. Our analysis is based on some delicate energy estimates.
文摘The numerical solution for a type of quasilinear wave equation is studied.The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved.The error of the difference solution is estimated.The theoretical results are controlled on a numerical example.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10875018 and 10773002)
文摘The spheroidal wave functions are found to have extensive applications in many branches of physics and mathematics. We use the perturbation method in supersymmetric quantum mechanics to obtain the analytic ground eigenvalue and the ground eigenfunction of the angular spheroidal wave equation at low frequency in a series form. Using this approach, the numerical determinations of the ground eigenvalue and the ground eigenfunction for small complex frequencies are also obtained.
基金supported by the National Natural Science Foundation of China (11171208)Shanghai Leading Academic Discipline Project (S30106)
文摘The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, which uses an orthogonal function system with a weight function as the basis function. Compared with the conventional moving least-squares (MLS) approximation, the algebraic equation system in the IMLS approximation is not ill-conditioned, and can be solved directly without deriving the inverse matrix. Because there are fewer coefficients in the IMLS than in the MLS approximation, fewer nodes are selected in the IEFG method than in the element-free Galerkin method. Thus, the IEFG method has a higher computing speed. In the IEFG method, the Galerkin weak form is employed to obtain a dis- cretized system equation, and the penalty method is applied to impose the essential boundary condition. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the wave equations and the boundary-initial conditions depend on time, the scal- ing parameter, number of nodes and the time step length are considered for the convergence study.
基金supported by the National Science and Technology Major Special Sub-project of China(No.2016ZX05024-001-008)the National Natural Science Foundation Joint Fund Prcject of China(No.U1562215).
文摘A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly to the second-order elastic wave equation.In view of this,based on the first-order CPML absorbing boundary condition,we propose a new CPML(NCPML)boundary which can be directly applied to the second-order wave equation.We first systematically extend the first-order CPML technique into second-order wave equations,neglecting the space-varying characteristics of the partial damping coefficient in the complex-frequency domain,avoiding the generation of convolution in the time domain.We then transform the technique back to the time domain through the inverse Fourier transform.Numerical simulation indicates that the space-varying characteristics of the attenuation factor have little influence on the absorption effect and increase the memory at the same time.A number of numerical examples show that the NCPML proposed in this study is effective in simulating elastic wave propagation,and this algorithm is more efficient and requires less memory allocation than the conventional PML absorbing boundary.
基金Project supported by the Anhui Key Laboratory of Information Materials and Devices (Anhui University),China
文摘In this paper, we make use of the auxiliary equation and the expanded mapping methods to find the new exact periodic solutions for (2+1)-dimensional dispersive long wave equations in mathematical physics, which are expressed by Jacobi elliptic functions, and obtain some new solitary wave solutions (m → 1). This method can also be used to explore new periodic wave solutions for other nonlinear evolution equations.