Efficient local structure‐preserving schemes for the RLW‐type equation

The multisymplectic schemes have been used in numerical simulations for the RLW‐type equation successfully. They well preserve the local geometric property, but not other local conservation laws. In this article, we propose three novel efficient local structure‐preserving schemes for the RLW‐type equation, which preserve the local energy exactly on any time‐space region and can produce richer information of the original problem. The schemes will be mass‐ and energy‐preserving as the equation is imposed on appropriate boundary conditions. Numerical experiments are presented to verify the efficiency and invariant‐preserving property of the schemes. Comparisons with the existing nonconservative schemes are made to show the behavior of the energy affects the behavior of the solution.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1678–1691, 2017     

Fourth‐order compact scheme for the one‐dimensional sine‐Gordon equation

Finite difference scheme to the generalized one‐dimensional sine‐Gordon equation is considered in this paper. After approximating the second order derivative in the space variable by the compact finite difference, we transform the sine‐Gordon equation into an initial‐value problem of a second‐order ordinary differential equation. Then Padé approximant is used to approximate the time derivatives. The resulting fully discrete nonlinear finite‐difference equation is solved by a predictor‐corrector scheme. Both Dirichlet and Neumann boundary conditions are considered in our proposed algorithm. Stability analysis and error estimate are given for homogeneous Dirichlet boundary value problems using energy method. Numerical results are given to verify the condition for stability and convergence and to examine the accuracy and efficiency of the proposed algorithm. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Local energy‐ and momentum‐preserving schemes for Klein‐Gordon‐Schrödinger equations and convergence analysis

In this article, we obtain local energy and momentum conservation laws for the Klein‐Gordon‐Schrödinger equations, which are independent of the boundary condition and more essential than the global conservation laws. Based on the rule that the numerical methods should preserve the intrinsic properties as much as possible, we propose local energy‐ and momentum‐preserving schemes for the equations. The merit of the proposed schemes is that the local energy/momentum conservation law is conserved exactly in any time‐space region. With suitable boundary conditions, the schemes will be charge‐ and energy‐/momentum‐preserving. Nonlinear analysis shows LEP schemes are unconditionally stable and the numerical solutions converge to the exact solutions with order . The theoretical properties are verified by numerical experiments. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1329–1351, 2017     

An enhanced‐physics‐based scheme for the NS‐α turbulence model

We study a new enhanced‐physics‐based numerical scheme for the NS‐alpha turbulence model that conserves both energy and helicity. Although most turbulence models (in the continuous case) conserve only energy, NS‐alpha is one of only a very few that also conserve helicity. This is one reason why it is becoming accepted as the most physically accurate turbulence model. However, no numerical scheme for NS‐alpha, until now, conserved both energy and helicity, and thus the advantage gained in physical accuracy by modeling with NS‐alpha could be lost in a computation. This report presents a finite element numerical scheme, and gives a rigorous analysis of its conservation properties, stability, solution existence, and convergence. A key feature of the analysis is the identification of the discrete energy and energy dissipation norms, and proofs that these norms are equivalent (provided a careful choice of filtering radius) in the discrete space to the usual energy and energy dissipation norms. Numerical experiments are given to demonstrate the effectiveness of the scheme over usual (helicity‐ignoring) schemes. A generalization of this scheme to a family of high‐order NS‐alpha‐deconvolution models, which combine the attractive physical properties of NS‐alpha with the high accuracy gained by combining α‐filtering with van Cittert approximate deconvolution. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Spatial behaviour of solutions of the dual‐phase‐lag heat equation

In this paper we study the spatial behaviour of solutions of some problems for the dual‐phase‐lag heat equation on a semi‐infinite cylinder. The theory of dual‐phase‐lag heat conduction leads to a hyperbolic partial differential equation with a third derivative with respect to time. First, we investigate the spatial evolution of solutions of an initial boundary‐value problem with zero boundary conditions on the lateral surface of the cylinder. Under a boundedness restriction on the initial data, an energy estimate is obtained. An upper bound for the amplitude term in this estimate in terms of the initial and boundary data is also established. For the case of zero initial conditions, a more explicit estimate is obtained which shows that solutions decay exponentially along certain spatial‐time lines. A class of non‐standard problems is also considered for which the temperature and its first two time derivatives at a fixed time T are assumed proportional to their initial values. Copyright © 2004 John Wiley & Sons, Ltd.     

New energy identities and super convergence analysis of the energy conserved splitting FDTD methods for 3D Maxwell's equations

The energy‐conserved splitting finite‐difference time‐domain (EC‐S‐FDTD) method has recently been proposed to solve the Maxwell equations with second order accuracy while numerically keep the L² energy conservation laws of the equations. In this paper, the EC‐S‐FDTD scheme for the 3D Maxwell equations is proved to be energy‐conserved and unconditionally stable in the discrete H¹ norm. The EC‐S‐FDTD scheme is of second‐order accuracy both in time step and spatial steps, which suggests the super‐convergence of this scheme in the discrete H¹ norm. And the divergence of the electric field of the EC‐S‐FDTD scheme in the discrete L² norm is second‐order accurate. Numerical experiments confirm our theoretical analysis. Copyright © 2012 John Wiley & Sons, Ltd.     

On the decay and blow‐up of solution for coupled nonlinear wave equations with nonlinear damping and source terms

In this work, we consider a nonlinear coupled wave equations with initial‐boundary value conditions and nonlinear damping and source terms. Under suitable assumptions on the damping terms and source terms and initial data in the stable set, we obtain that the decay estimates of the energy function is exponential or polynomial by using Nakao's method. By using the energy method, we obtain the blow‐up result of solution with some positive or nonpositive initial energy. Copyright © 2016 John Wiley & Sons, Ltd.     

Dynamical behavior of fractional‐order energy‐saving and emission‐reduction system and its discretization

This paper attempts to explore dynamical behavior and mathematical properties of the three‐dimensional fractional‐order energy‐saving and emission‐reduction system. Theoretically, the conditions of local stability of fractional‐order system's equilibrium points are obtained. Numerical investigations on the dynamics of this system are carried out, and the existence of the asymptotically stable attractor is found. Combined with the fractional‐order subsystem, we discuss the relationship between energy‐saving and emission‐reduction and economic growth, and carbon emissions and economic growth. Furthermore, we discretize the fractional‐order system and give necessary and sufficient conditions of its stabilization. It is shown that the stability of the discretization system is impacted by the system's fractional parameter. Numerical simulations show the richer dynamical behavior of the fractional‐order system and verify the theoretical results. Recommendations for Resource Managers

• The impact of carbon emissions on economic growth is one of the main reasons for energy‐saving and emission‐reduction.
• Control measures on people's low‐carbon life through government intervention are required to protect the natural environment.
• New energy‐saving and emission‐reduction technologies should be implemented to achieve sustainable social and economic development.

     

Radiation conditions for the linear water‐wave problem in periodic channels

We study the well‐posedness of the linearized water‐wave problem in a periodic channel with fixed or freely floating compact bodies. Floquet–Bloch–Gelfand‐transform techniques lead to a generalized spectral problem with quadratic dependence on a complex parameter, and the asymptotics of the solutions at infinity can be described using Floquet waves. These are constructed from Jordan chains, which are related with the eigenvalues of the quadratic pencil and which are calculated in the paper in some typical cases. Posing proper radiation conditions requires a careful study of spaces of incoming and outgoing waves, especially in the threshold situation. This is done with the help of a certain skew‐Hermitian form q , which is closely related to the Umov–Poynting vector of energy transportation. Our radiation conditions make the problem operator into a Fredholm operator of index zero and provides natural (energy) classification of outgoing/incoming waves. They also lead to a novel, most natural properties and interpretation of the scattering matrix, which becomes unitary and symmetric even at threshold.     

A second‐order energy stable backward differentiation formula method for the epitaxial thin film equation with slope selection

In this article, we study a new second‐order energy stable Backward Differentiation Formula (BDF) finite difference scheme for the epitaxial thin film equation with slope selection (SS). One major challenge for higher‐order‐in‐time temporal discretizations is how to ensure an unconditional energy stability without compromising numerical efficiency or accuracy. We propose a framework for designing a second‐order numerical scheme with unconditional energy stability using the BDF method with constant coefficient stabilizing terms. Based on the unconditional energy stability property that we establish, we derive an stability for the numerical solution and provide an optimal convergence analysis. To deal with the highly nonlinear four‐Laplacian term at each time step, we apply efficient preconditioned steepest descent and preconditioned nonlinear conjugate gradient algorithms to solve the corresponding nonlinear system. Various numerical simulations are presented to demonstrate the stability and efficiency of the proposed schemes and solvers. Comparisons with other second‐order schemes are presented.     

Recurrent solutions of the linearly coupled complex cubic‐quintic Ginzburg‐Landau equations

In this paper, we will establish the bounded solutions, periodic solutions, quasiperiodic solutions, almost periodic solutions, and almost automorphic solutions for linearly coupled complex cubic‐quintic Ginzburg‐Landau equations, under suitable conditions. The main difficulty is the nonlinear terms in the equations that are not Lipschitz‐continuity, traditional methods cannot deal with the difficulty in our problem. We overcome this difficulty by the Galerkin approach, energy estimate method, and refined inequality technique.     

Global non‐existence of solutions of a class of wave equations with non‐linear damping and source terms

In this paper we consider the non‐linear wave equation a,b>0, associated with initial and Dirichlet boundary conditions. We prove, under suitable conditions on α,β,m,p and for negative initial energy, a global non‐existence theorem. This improves a result by Yang (Math. Meth. Appl. Sci. 2002; 25 :825–833), who requires that the initial energy be sufficiently negative and relates the global non‐existence of solutions to the size of Ω. Copyright © 2004 John Wiley & Sons, Ltd.     

A fourth‐order compact finite difference scheme for solving an N‐carrier system with Neumann boundary conditions

In this study, we develop a fourth‐order compact finite difference scheme for solving a model of energy exchanges in a generalized N‐carrier system with heat sources and Neumann boundary conditions, which extends the concept of the well‐known parabolic two‐step model for microheat transfer. By using the matrix analysis, the compact finite difference numerical scheme is shown to be unconditionally stable. The accuracy of the solution obtained by the scheme is tested by a numerical example. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Riesz basis property of root vectors of non‐self‐adjoint operators generated by aircraft wing model in subsonic airflow

This paper is the third in a series of several works devoted to the asymptotic and spectral analysis of a model of an aircraft wing in a subsonic air flow. This model has been developed in the Flight Systems Research Center of UCLA and is presented in the works by Balakrishnan. The model is governed by a system of two coupled integro‐differential equations and a two‐parameter family of boundary conditions modeling the action of the self‐straining actuators. The differential parts of the above equations form a coupled linear hyperbolic system; the integral parts are of the convolution type. The system of equations of motion is equivalent to a single operator evolution–convolution equation in the energy space. The Laplace transform of the solution of this equation can be represented in terms of the so‐called generalized resolvent operator, which is an operator‐valued function of the spectral parameter. This generalized resolvent operator is a finite‐meromorphic function on the complex plane having the branch cut along the negative real semi‐axis. Its poles are precisely the aeroelastic modes and the residues at these poles are the projectors on the generalized eigenspaces. In the first two papers (see [33, 34]) and in the present one, our main object of interest is the dynamics generator of the differential parts of the system. This generator is a non‐self‐adjoint operator in the energy space with a purely discrete spectrum. In the first paper, we have shown that the spectrum consists of two branches, and have derived their precise spectral asymptotics with respect to the eigenvalue number. In the second paper, we have derived the asymptotical approximations for the mode shapes. Based on the asymptotical results of the first two papers, in the present paper, we (a) prove that the set of the generalized eigenvectors of the aforementioned differential operator is complete in the energy space; (b) construct the set of vectors which is biorthogonal to the set of the generalized eigenvectors in the case when there might be not only eigenvectors but associate vectors as well; and (c) prove that the set of the generalized eigenvectors forms a Riesz basis in the energy space. To prove the main result of the paper, we made use of the Nagy–Foias functional model for non‐self‐adjoint operators. The results of all three papers will be important for the reconstruction of the solution of the original initial‐boundary‐value problem from its Laplace transform in the forthcoming papers. Copyright © 2000 John Wiley & Sons, Ltd.     

Energy decay estimates and infinite blow‐up phenomena for a strongly damped semilinear wave equation with logarithmic nonlinear source

This paper deals with the energy decay estimates and infinite blow‐up phenomena for a strongly damped semilinear wave equation with logarithmic nonlinear source term under null Dirichlet boundary condition. By constructing a new family of potential wells, together with logarithmic Sobolev inequality and perturbation energy technique, we establish sufficient conditions to guarantee the solution exists globally or occurs infinite blow‐up and derive the polynomial or exponential energy decay estimates under some appropriate conditions.     

The Cauchy problem of a fluid‐particle interaction model with external forces

In this paper, we consider the Cauchy problem of a fluid‐particle interaction model with external forces. We first construct the asymptotic profile of the system. The global existence and uniqueness theorem for the solution near the profile is given. Finally, optimal decay rate of the solution to the background profile is obtained by combining the decay rate analysis of a linearized equation with energy estimates for the nonlinear terms. The main method used in this paper is the energy method combining with the macro‐micro decomposition.     

Homogenization of the Signorini boundary‐value problem in a thick junction and boundary integral equations for the homogenized problem

We consider a mixed boundary‐value problem for the Poisson equation in a thick junction Ω_ε which is the union of a domain Ω₀ and a large number of ε—periodically situated thin cylinders. The non‐uniform Signorini conditions are given on the lateral surfaces of the cylinders. The asymptotic analysis of this problem is done as ε→0, i.e. when the number of the thin cylinders infinitely increases and their thickness tends to zero. We prove a convergence theorem and show that the non‐uniform Signorini boundary conditions are transformed in the limiting variational inequalities in the region that is filled up by the thin cylinders as ε→0. The convergence of the energy integrals is proved as well. The existence and uniqueness of the solution to this non‐standard limit problem is established. This solution can be constructed by using a penalty formulation and successive iteration. For some subclass, these problems can be reduced to an obstacle problem in Ω₀ and an appropriate postprocessing. The equations in Ω₀ finally are also treated with boundary integral equations. Copyright © 2010 John Wiley & Sons, Ltd.     

On resonant interactions of gravity‐capillary waves without energy exchange

We consider resonant triad interactions of gravity‐capillary waves and investigate in detail special resonant triads that exchange no energy during their interactions so that the wave amplitudes remain constant in time. After writing the resonance conditions in terms of two parameters (or two angles of wave propagation), we first identify a region in the two‐dimensional parameter space, where resonant triads can be always found, and then describe the variations of resonant wavenumbers and wave frequencies over the resonance region. Using the amplitude equations recovered from a Hamiltonian formulation for water waves, it is shown that any resonant triad inside the resonance region can interact without energy exchange if the initial wave amplitudes and relative phase satisfy the two conditions for fixed point solutions of the amplitude equations. Furthermore, it is shown that the symmetric resonant triad exchanging no energy forms a transversely modulated traveling wave field, which can be considered a two‐dimensional generalization of Wilton ripples.     

Persistence of undercompressive phase boundaries for isothermal Euler equations including configurational forces and surface tension

The persistence of subsonic phase boundaries in a multidimensional Van der Waals fluid is analyzed. The phase boundary is considered as a sharp free boundary that connects liquid and vapor bulk phase dynamics given by the isothermal Euler equations. The evolution of the boundary is driven by effects of configurational forces as well as surface tension. To analyze this problem, the equations and trace conditions are linearized such that one obtains a general hyperbolic initial boundary value problem with higher‐order boundary conditions. A global existence theorem for the linearized system with constant coefficients is shown. The proof relies on the normal mode analysis and a linear form in suitable spaces that is defined using an associated adjoint problem. Especially, the associated adjoint problem satisfies the uniform backward in time Kreiss–Lopatinski? condition. A new energy‐like estimate that also includes surface energy terms leads finally to the uniqueness and regularity for the found solutions of the problem in weighted spaces. Copyright © 2016 John Wiley & Sons, Ltd.     

Error analysis for solving the Korteweg‐de Vries equation by a Legendre pseudo‐spectral method

A Legendre pseudo‐spectral method is proposed for the Korteweg‐de Vries equation with nonperiodic boundary conditions. Appropriate base functions are chosen to get an efficient algorithm. Error analysis is given for both semi‐discrete and fully discrete schemes. The numerical results confirm to the theoretical analysis. © (2000) John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 513–534, (2000)