In this study,the wave motion in elastodynamics for unbounded media is modeled using an unsplit-field perfectly matched layer(PML)formulation that is solved by employing an isogeometric analysis(IGA).In the adopted co...In this study,the wave motion in elastodynamics for unbounded media is modeled using an unsplit-field perfectly matched layer(PML)formulation that is solved by employing an isogeometric analysis(IGA).In the adopted combination,the non-uniform rational B-spline(NURBS)functions are employed as basis functions.Moreover,the unbounded and artificial domains,defined in the PML method,are contained in a single patch domain.Based on the proposed scheme,the approximation of the geometry problem is set in a new scheme in which the PML’s absorbing and attenuation properties and the description of traveling waves can be represented.This includes a higher continuity and smoother approximation of the computed domain.As high-order NURBS basis functions are non-interpolatory,a penalty method is present to apply a time-dependent displacement load.The performance of the NURBS-based PML is analyzed through numerical examples for 1D and 2D domains,considering homogeneous and heterogeneous media.Further,we verify the long-time numerical stability of the present method.The developed method can be used to simulate hypothetical stratified domains commonly encountered in soil-structure interaction analyses.展开更多
Human motion modeling is a core technology in computer animation,game development,and humancomputer interaction.In particular,generating natural and coherent in-between motion using only the initial and terminal frame...Human motion modeling is a core technology in computer animation,game development,and humancomputer interaction.In particular,generating natural and coherent in-between motion using only the initial and terminal frames remains a fundamental yet unresolved challenge.Existing methods typically rely on dense keyframe inputs or complex prior structures,making it difficult to balance motion quality and plausibility under conditions such as sparse constraints,long-term dependencies,and diverse motion styles.To address this,we propose a motion generation framework based on a frequency-domain diffusion model,which aims to better model complex motion distributions and enhance generation stability under sparse conditions.Our method maps motion sequences to the frequency domain via the Discrete Cosine Transform(DCT),enabling more effective modeling of low-frequency motion structures while suppressing high-frequency noise.A denoising network based on self-attention is introduced to capture long-range temporal dependencies and improve global structural awareness.Additionally,a multi-objective loss function is employed to jointly optimize motion smoothness,pose diversity,and anatomical consistency,enhancing the realism and physical plausibility of the generated sequences.Comparative experiments on the Human3.6M and LaFAN1 datasets demonstrate that our method outperforms state-of-the-art approaches across multiple performance metrics,showing stronger capabilities in generating intermediate motion frames.This research offers a new perspective and methodology for human motion generation and holds promise for applications in character animation,game development,and virtual interaction.展开更多
The precise tuning of magnetic nanoparticle size and magnetic domains,thereby shaping magnetic properties.However,the dynamic evolution mechanisms of magnetic domain configurations in relation to electromagnetic(EM)at...The precise tuning of magnetic nanoparticle size and magnetic domains,thereby shaping magnetic properties.However,the dynamic evolution mechanisms of magnetic domain configurations in relation to electromagnetic(EM)attenuation behavior remain poorly understood.To address this gap,a thermodynamically controlled periodic coordination strategy is proposed to achieve precise modulation of magnetic nanoparticle spacing.This approach unveils the evolution of magnetic domain configurations,progressing from individual to coupled and ultimately to crosslinked domain configurations.A unique magnetic coupling phenomenon surpasses the Snoek limit in low-frequency range,which is observed through micromagnetic simulation.The crosslinked magnetic configuration achieves effective low-frequency EM wave absorption at 3.68 GHz,encompassing nearly the entire C-band.This exceptional magnetic interaction significantly enhances radar camouflage and thermal insulation properties.Additionally,a robust gradient metamaterial design extends coverage across the full band(2–40 GHz),effectively mitigating the impact of EM pollution on human health and environment.This comprehensive study elucidates the evolution mechanisms of magnetic domain configurations,addresses gaps in dynamic magnetic modulation,and provides novel insights for the development of high-performance,low-frequency EM wave absorption materials.展开更多
The author studies the infinite element method for the boundary value problems of second order elliptic equations on unbounded and multiply connected domains. The author makes a partition of the domain into infinite n...The author studies the infinite element method for the boundary value problems of second order elliptic equations on unbounded and multiply connected domains. The author makes a partition of the domain into infinite number of elements. Without dividing the domain, as usual, into a bounded one and an exterior one, he derives an initial value problem of an ordinary differential equation for the combined stiffness matrix, then obtains the approximate solution with a small amount of computer work. Numerical examples are given.展开更多
By the interpolating inequality and a priori estimates in the weighted space,the existence of the global solutions for the Ginzburg-Landau equation coupled with the BBM equation in an unbounded domain is considered, a...By the interpolating inequality and a priori estimates in the weighted space,the existence of the global solutions for the Ginzburg-Landau equation coupled with the BBM equation in an unbounded domain is considered, and the existence of the maximal attractor is obtained.展开更多
In this paper, we prove the existence of random attractors for a stochastic reaction-diffusion equation with distribution derivatives on unbounded domains. The nonlinearity is dissipative for large values of the state...In this paper, we prove the existence of random attractors for a stochastic reaction-diffusion equation with distribution derivatives on unbounded domains. The nonlinearity is dissipative for large values of the state and the stochastic nature of the equation appears spatially distributed temporal white noise. The stochastic reaction-diffusion equation is recast as a continuous random dynamical system and asymptotic compactness for this demonstrated by using uniform estimates far-field values of solutions. The results are new and appear to be optimal.展开更多
Some specific non-isotropic Jacobi approximations in multiple-dimensions are investigated, which are used for numerical solutions of differential equations on various unbounded domains. The convergence of proposed sch...Some specific non-isotropic Jacobi approximations in multiple-dimensions are investigated, which are used for numerical solutions of differential equations on various unbounded domains. The convergence of proposed schemes are proved. Some efficient algorithms are provided. Numerical results are presented to illustrate the efficiency of this new approach.展开更多
In this paper, we use the concentration-compactness principle together with the Mountain Pass Lemma to get the existence of nontrivial solutions and the existence of infinitely many solutions of the problem need not b...In this paper, we use the concentration-compactness principle together with the Mountain Pass Lemma to get the existence of nontrivial solutions and the existence of infinitely many solutions of the problem need not be compact operators from E to R~1.展开更多
In this paper, a new domain decomposition method based on the natural boundary reduction, which solves wave problems over an unbounded domain, is suggestted. An circular artifcial boundary is introduced. The original ...In this paper, a new domain decomposition method based on the natural boundary reduction, which solves wave problems over an unbounded domain, is suggestted. An circular artifcial boundary is introduced. The original unbounded domain is divided into two subdomains, an internal bounded region and external unbounded region outside the artificial boundary. A Dirichlet-Neumann(D-N) alternating iteration algorithm is constructed. We prove that the algorithm is equavilent to preconditional Richardson iteration method. Numerical studies are performed by finite element method. The numerical results show that the convergence rate of the discrete D-N iteration is independent of the fnite element mesh size.展开更多
In this paper we investigate the existence and stability of periodic solutions(on a half-line R_(+))and almost periodic solutions on the whole line time-axis R to the Boussinesq system on several classes of unbounded ...In this paper we investigate the existence and stability of periodic solutions(on a half-line R_(+))and almost periodic solutions on the whole line time-axis R to the Boussinesq system on several classes of unbounded domains of R^(n) in the framework of interpolation spaces.For the linear Boussinesq system we combine the L^(p)—L^(q)-smoothing estimates and interpolation functors to prove the existence of bounded mild solutions.Then,we prove the existence of periodic solutions by invoking Massera’s principle.We also prove the existence of almost periodic solutions.Then we use the results of the linear Boussinesq system to establish the existence,uniqueness and stability of the small periodic and almost periodic solutions to the Boussinesq system using fixed point arguments and interpolation spaces.Our results cover and extend the previous ones obtained in[13,34,38].展开更多
The following nonlinear hyperbolic equation is discussed in this paper:where A= -? + Iandx∈Rn. The model comes from the transverse deflection equation of an extensible beam. We prove that there exists a unique local ...The following nonlinear hyperbolic equation is discussed in this paper:where A= -? + Iandx∈Rn. The model comes from the transverse deflection equation of an extensible beam. We prove that there exists a unique local solution of the above equation as M depends on x.展开更多
This paper is focused on studying the structure of solutions bounded from below to degenerate elliptic equations with Neumann and Dirichlet boundary conditions in unbounded domains.After establishing the weak maximum ...This paper is focused on studying the structure of solutions bounded from below to degenerate elliptic equations with Neumann and Dirichlet boundary conditions in unbounded domains.After establishing the weak maximum principles,the global boundary Holder estimates and the boundary Harnack inequalities of the equations,we show that all solutions bounded from below are linear combinations of two special solutions(exponential growth at one end and exponential decay at the other)with a bounded solution to the degenerate equations.展开更多
On the basis of introducing the modified Cauchy kernel, we discuss the Hoelder continuity of the Cauchy-type singular integral operator on unbounded domains for regular functions by dividing into the following three c...On the basis of introducing the modified Cauchy kernel, we discuss the Hoelder continuity of the Cauchy-type singular integral operator on unbounded domains for regular functions by dividing into the following three cases: two points are on the boundary of region; one point is on the boundary and another point is in the interior(or exterior) of the region; two points are in the interior (or exterior) of the region.展开更多
In this article,we investigate the longtime behavior for the following nonautonomous nonclassical parabolic equations on unbounded domain ut−∆ut−∆u+λu=f(x,u(x,t−ρ(t)))+g(x,t).Under some suitable conditions on the d...In this article,we investigate the longtime behavior for the following nonautonomous nonclassical parabolic equations on unbounded domain ut−∆ut−∆u+λu=f(x,u(x,t−ρ(t)))+g(x,t).Under some suitable conditions on the delay term f and the non-autonomous forcing term g,we prove the existence of uniform attractors in Banach space CH1(RN)for the multivalued process generated by non-autonomous nonclassical parabolic equations with delays in unbounded domain.展开更多
Since [1] established the Pohozaev identity in bounded domains, this identity has been the principal tool to deal with the non-existence of the equation
In this article, we discuss the approximate method of solving the Riemann-Hilbert boundary value problem for nonlinear uniformly elliptic complex equation of first order (0.1) with the boundary conditions (0.2) in a m...In this article, we discuss the approximate method of solving the Riemann-Hilbert boundary value problem for nonlinear uniformly elliptic complex equation of first order (0.1) with the boundary conditions (0.2) in a multiply connected unbounded domain D, the above boundary value problem will be called Problem A. If the complex Equation (0.1) satisfies the conditions similar to Condition C of (1.1), and the boundary condition (0.2) satisfies the conditions similar to (1.5), then we can obtain approximate solutions of the boundary value problems (0.1) and (0.2). Moreover the error estimates of approximate solutions for the boundary value problem is also given. The boundary value problem possesses many applications in mechanics and physics etc., for instance from (5.114) and (5.115), Chapter VI, [1], we see that Problem A of (0.1) possesses the important application to the shell and elasticity.展开更多
In [1], I. N. Vekua propose the Poincaré problem for some second order elliptic equations, but it can not be solved. In [2], the authors discussed the boundary value problem for nonlinear elliptic equations of se...In [1], I. N. Vekua propose the Poincaré problem for some second order elliptic equations, but it can not be solved. In [2], the authors discussed the boundary value problem for nonlinear elliptic equations of second order in some bounded domains. In this article, the Poincaré boundary value problem for general nonlinear elliptic equations of second order in unbounded multiply connected domains have been completely investigated. We first provide the formulation of the above boundary value problem and corresponding modified well posed-ness. Next we obtain the representation theorem and a priori estimates of solutions for the modified problem. Finally by the above estimates of solutions and the Schauder fixed-point theorem, the solvability results of the above Poincaré problem for the nonlinear elliptic equations of second order can be obtained. The above problem possesses many applications in mechanics and physics and so on.展开更多
We consider fully nonlinear equations of the formF(z,u,Du,D<sup>2</sup>u) = F(x,y,u,u<sub>z</sub>,u<sub>y</sub>,u<sub>z</sub>z,u<sub>z</sub>y,u<sub>...We consider fully nonlinear equations of the formF(z,u,Du,D<sup>2</sup>u) = F(x,y,u,u<sub>z</sub>,u<sub>y</sub>,u<sub>z</sub>z,u<sub>z</sub>y,u<sub>y</sub>y) = 0 (1)in unbounded open subset G = R<sup>2</sup>\Ωof the plane R<sup>2</sup>,where F is a real continuous function on U = G×R×R<sup>2</sup>×R<sup>3</sup> and Ω= Ω<sub>i</sub>,Ω<sub>i</sub> is a simply connected region (i=1,2,"",N) . We assume the function F hascontinuous partial der ivatives F<sub>u<sub>z</sub>z</sub>, F<sub>u<sub>z</sub>y</sub>, F<sub>u<sub>y</sub>y</sub>. on U.For a real function r C( G) a real function u(x,y) is called a solution of (1) satisfyingu = r on G,(2)if there exists a constant P0】2 such that u C<sup>1</sup> ( ) W<sup>2.p</sup><sub>Loc</sub>0 (G) satisfies (1) almost everywhere and (2)in the common sense.The method for treating the above exterior Dirichlet problem in a given unbounded region is as fol-展开更多
A coupling method of finite element and infinite large element is proposed for the numerical solution of an eigenvalue problem in unbounded domains in this paper. With some conditions satisfied, the considered problem...A coupling method of finite element and infinite large element is proposed for the numerical solution of an eigenvalue problem in unbounded domains in this paper. With some conditions satisfied, the considered problem is proved to have discrete spectra. Several numerical experiments are presented. The results demonstrate the feasibility of the proposed method.展开更多
文摘In this study,the wave motion in elastodynamics for unbounded media is modeled using an unsplit-field perfectly matched layer(PML)formulation that is solved by employing an isogeometric analysis(IGA).In the adopted combination,the non-uniform rational B-spline(NURBS)functions are employed as basis functions.Moreover,the unbounded and artificial domains,defined in the PML method,are contained in a single patch domain.Based on the proposed scheme,the approximation of the geometry problem is set in a new scheme in which the PML’s absorbing and attenuation properties and the description of traveling waves can be represented.This includes a higher continuity and smoother approximation of the computed domain.As high-order NURBS basis functions are non-interpolatory,a penalty method is present to apply a time-dependent displacement load.The performance of the NURBS-based PML is analyzed through numerical examples for 1D and 2D domains,considering homogeneous and heterogeneous media.Further,we verify the long-time numerical stability of the present method.The developed method can be used to simulate hypothetical stratified domains commonly encountered in soil-structure interaction analyses.
基金supported by the National Natural Science Foundation of China(Grant No.72161034).
文摘Human motion modeling is a core technology in computer animation,game development,and humancomputer interaction.In particular,generating natural and coherent in-between motion using only the initial and terminal frames remains a fundamental yet unresolved challenge.Existing methods typically rely on dense keyframe inputs or complex prior structures,making it difficult to balance motion quality and plausibility under conditions such as sparse constraints,long-term dependencies,and diverse motion styles.To address this,we propose a motion generation framework based on a frequency-domain diffusion model,which aims to better model complex motion distributions and enhance generation stability under sparse conditions.Our method maps motion sequences to the frequency domain via the Discrete Cosine Transform(DCT),enabling more effective modeling of low-frequency motion structures while suppressing high-frequency noise.A denoising network based on self-attention is introduced to capture long-range temporal dependencies and improve global structural awareness.Additionally,a multi-objective loss function is employed to jointly optimize motion smoothness,pose diversity,and anatomical consistency,enhancing the realism and physical plausibility of the generated sequences.Comparative experiments on the Human3.6M and LaFAN1 datasets demonstrate that our method outperforms state-of-the-art approaches across multiple performance metrics,showing stronger capabilities in generating intermediate motion frames.This research offers a new perspective and methodology for human motion generation and holds promise for applications in character animation,game development,and virtual interaction.
基金supported by the National Natural Science Foundation of China(22265021,52231007,and 12327804)the Aeronautical Science Foundation of China(2020Z056056003)Jiangxi Provincial Natural Science Foundation(20232BAB212004).
文摘The precise tuning of magnetic nanoparticle size and magnetic domains,thereby shaping magnetic properties.However,the dynamic evolution mechanisms of magnetic domain configurations in relation to electromagnetic(EM)attenuation behavior remain poorly understood.To address this gap,a thermodynamically controlled periodic coordination strategy is proposed to achieve precise modulation of magnetic nanoparticle spacing.This approach unveils the evolution of magnetic domain configurations,progressing from individual to coupled and ultimately to crosslinked domain configurations.A unique magnetic coupling phenomenon surpasses the Snoek limit in low-frequency range,which is observed through micromagnetic simulation.The crosslinked magnetic configuration achieves effective low-frequency EM wave absorption at 3.68 GHz,encompassing nearly the entire C-band.This exceptional magnetic interaction significantly enhances radar camouflage and thermal insulation properties.Additionally,a robust gradient metamaterial design extends coverage across the full band(2–40 GHz),effectively mitigating the impact of EM pollution on human health and environment.This comprehensive study elucidates the evolution mechanisms of magnetic domain configurations,addresses gaps in dynamic magnetic modulation,and provides novel insights for the development of high-performance,low-frequency EM wave absorption materials.
基金This work was supported by the China State Major Key Project for Basic Researches Science Fund of the Ministry of Education
文摘The author studies the infinite element method for the boundary value problems of second order elliptic equations on unbounded and multiply connected domains. The author makes a partition of the domain into infinite number of elements. Without dividing the domain, as usual, into a bounded one and an exterior one, he derives an initial value problem of an ordinary differential equation for the combined stiffness matrix, then obtains the approximate solution with a small amount of computer work. Numerical examples are given.
文摘By the interpolating inequality and a priori estimates in the weighted space,the existence of the global solutions for the Ginzburg-Landau equation coupled with the BBM equation in an unbounded domain is considered, and the existence of the maximal attractor is obtained.
文摘In this paper, we prove the existence of random attractors for a stochastic reaction-diffusion equation with distribution derivatives on unbounded domains. The nonlinearity is dissipative for large values of the state and the stochastic nature of the equation appears spatially distributed temporal white noise. The stochastic reaction-diffusion equation is recast as a continuous random dynamical system and asymptotic compactness for this demonstrated by using uniform estimates far-field values of solutions. The results are new and appear to be optimal.
文摘Some specific non-isotropic Jacobi approximations in multiple-dimensions are investigated, which are used for numerical solutions of differential equations on various unbounded domains. The convergence of proposed schemes are proved. Some efficient algorithms are provided. Numerical results are presented to illustrate the efficiency of this new approach.
文摘In this paper, we use the concentration-compactness principle together with the Mountain Pass Lemma to get the existence of nontrivial solutions and the existence of infinitely many solutions of the problem need not be compact operators from E to R~1.
文摘In this paper, a new domain decomposition method based on the natural boundary reduction, which solves wave problems over an unbounded domain, is suggestted. An circular artifcial boundary is introduced. The original unbounded domain is divided into two subdomains, an internal bounded region and external unbounded region outside the artificial boundary. A Dirichlet-Neumann(D-N) alternating iteration algorithm is constructed. We prove that the algorithm is equavilent to preconditional Richardson iteration method. Numerical studies are performed by finite element method. The numerical results show that the convergence rate of the discrete D-N iteration is independent of the fnite element mesh size.
基金financially supported by the Vietnam National Foundation for Science and Technology Development under grant number 101.02-2021.04financially supported by Vietnam Ministry of Education and Training under Project B2022-BKA-06.
文摘In this paper we investigate the existence and stability of periodic solutions(on a half-line R_(+))and almost periodic solutions on the whole line time-axis R to the Boussinesq system on several classes of unbounded domains of R^(n) in the framework of interpolation spaces.For the linear Boussinesq system we combine the L^(p)—L^(q)-smoothing estimates and interpolation functors to prove the existence of bounded mild solutions.Then,we prove the existence of periodic solutions by invoking Massera’s principle.We also prove the existence of almost periodic solutions.Then we use the results of the linear Boussinesq system to establish the existence,uniqueness and stability of the small periodic and almost periodic solutions to the Boussinesq system using fixed point arguments and interpolation spaces.Our results cover and extend the previous ones obtained in[13,34,38].
文摘The following nonlinear hyperbolic equation is discussed in this paper:where A= -? + Iandx∈Rn. The model comes from the transverse deflection equation of an extensible beam. We prove that there exists a unique local solution of the above equation as M depends on x.
文摘This paper is focused on studying the structure of solutions bounded from below to degenerate elliptic equations with Neumann and Dirichlet boundary conditions in unbounded domains.After establishing the weak maximum principles,the global boundary Holder estimates and the boundary Harnack inequalities of the equations,we show that all solutions bounded from below are linear combinations of two special solutions(exponential growth at one end and exponential decay at the other)with a bounded solution to the degenerate equations.
基金the National Natural Science Foundation of China(Grant Nos.1140116211301136)the ScienceFoundation of Hebei Province(Grant No.A2014208158)
文摘On the basis of introducing the modified Cauchy kernel, we discuss the Hoelder continuity of the Cauchy-type singular integral operator on unbounded domains for regular functions by dividing into the following three cases: two points are on the boundary of region; one point is on the boundary and another point is in the interior(or exterior) of the region; two points are in the interior (or exterior) of the region.
基金Supported by the 2018 research funding of higher education of Gansu province project[2018B-075]
文摘In this article,we investigate the longtime behavior for the following nonautonomous nonclassical parabolic equations on unbounded domain ut−∆ut−∆u+λu=f(x,u(x,t−ρ(t)))+g(x,t).Under some suitable conditions on the delay term f and the non-autonomous forcing term g,we prove the existence of uniform attractors in Banach space CH1(RN)for the multivalued process generated by non-autonomous nonclassical parabolic equations with delays in unbounded domain.
基金This work is supported in port by the Foundation of Zhongshan University Advanced Research Center.
文摘Since [1] established the Pohozaev identity in bounded domains, this identity has been the principal tool to deal with the non-existence of the equation
文摘In this article, we discuss the approximate method of solving the Riemann-Hilbert boundary value problem for nonlinear uniformly elliptic complex equation of first order (0.1) with the boundary conditions (0.2) in a multiply connected unbounded domain D, the above boundary value problem will be called Problem A. If the complex Equation (0.1) satisfies the conditions similar to Condition C of (1.1), and the boundary condition (0.2) satisfies the conditions similar to (1.5), then we can obtain approximate solutions of the boundary value problems (0.1) and (0.2). Moreover the error estimates of approximate solutions for the boundary value problem is also given. The boundary value problem possesses many applications in mechanics and physics etc., for instance from (5.114) and (5.115), Chapter VI, [1], we see that Problem A of (0.1) possesses the important application to the shell and elasticity.
文摘In [1], I. N. Vekua propose the Poincaré problem for some second order elliptic equations, but it can not be solved. In [2], the authors discussed the boundary value problem for nonlinear elliptic equations of second order in some bounded domains. In this article, the Poincaré boundary value problem for general nonlinear elliptic equations of second order in unbounded multiply connected domains have been completely investigated. We first provide the formulation of the above boundary value problem and corresponding modified well posed-ness. Next we obtain the representation theorem and a priori estimates of solutions for the modified problem. Finally by the above estimates of solutions and the Schauder fixed-point theorem, the solvability results of the above Poincaré problem for the nonlinear elliptic equations of second order can be obtained. The above problem possesses many applications in mechanics and physics and so on.
基金The project is supported by the Natural Science Foundation of Fujian,P.R.China.
文摘We consider fully nonlinear equations of the formF(z,u,Du,D<sup>2</sup>u) = F(x,y,u,u<sub>z</sub>,u<sub>y</sub>,u<sub>z</sub>z,u<sub>z</sub>y,u<sub>y</sub>y) = 0 (1)in unbounded open subset G = R<sup>2</sup>\Ωof the plane R<sup>2</sup>,where F is a real continuous function on U = G×R×R<sup>2</sup>×R<sup>3</sup> and Ω= Ω<sub>i</sub>,Ω<sub>i</sub> is a simply connected region (i=1,2,"",N) . We assume the function F hascontinuous partial der ivatives F<sub>u<sub>z</sub>z</sub>, F<sub>u<sub>z</sub>y</sub>, F<sub>u<sub>y</sub>y</sub>. on U.For a real function r C( G) a real function u(x,y) is called a solution of (1) satisfyingu = r on G,(2)if there exists a constant P0】2 such that u C<sup>1</sup> ( ) W<sup>2.p</sup><sub>Loc</sub>0 (G) satisfies (1) almost everywhere and (2)in the common sense.The method for treating the above exterior Dirichlet problem in a given unbounded region is as fol-
基金Supported partial by the Natinal Science Foundation of China under Grant No. 10401020 and Grant No. 10471073.
文摘A coupling method of finite element and infinite large element is proposed for the numerical solution of an eigenvalue problem in unbounded domains in this paper. With some conditions satisfied, the considered problem is proved to have discrete spectra. Several numerical experiments are presented. The results demonstrate the feasibility of the proposed method.