A POSTERIORI ERROR ESTIMATE OF THE DSD METHOD FOR FIRST-ORDER HYPERBOLIC EQUATIONS

A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.

Numerical Analysis of Upwind Difference Schemes for Two-Dimensional First-Order Hyperbolic Equations with Variable Coefficients

In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.

THE OPTIMAL LARGE TIME BEHAVIOR OF3D QUASILINEAR HYPERBOLIC EQUATIONS WITH NONLINEAR DAMPING

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.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

AENO:a Novel Reconstruction Method in Conjunction with ADER Schemes for Hyperbolic Equations

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.

Efficient Sparse-Grid Implementation of a Fifth-Order Multi-resolution WENO Scheme for Hyperbolic Equations

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.

Dirichlet-to-Neumann Map for a Hyperbolic Equation

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.

Linear Quadratic Optimal Control for Systems Governed by First-Order Hyperbolic Partial Differential Equations

This paper focuses on linear-quadratic(LQ)optimal control for a class of systems governed by first-order hyperbolic partial differential equations(PDEs).Different from most of the previous works,an approach of discretization-then-continuousization is proposed in this paper to cope with the infinite-dimensional nature of PDE systems.The contributions of this paper consist of the following aspects:(1)The differential Riccati equations and the solvability condition of the LQ optimal control problems are obtained via the discretization-then-continuousization method.(2)A numerical calculation way of the differential Riccati equations and a practical design way of the optimal controller are proposed.Meanwhile,the relationship between the optimal costate and the optimal state is established by solving a set of forward and backward partial difference equations(FBPDEs).(3)The correctness of the method used in this paper is verified by a complementary continuous method and the comparative analysis with the existing operator results is presented.It is shown that the proposed results not only contain the classic results of the standard LQ control problem of systems governed by ordinary differential equations as a special case,but also support the existing operator results and give a more convenient form of computation.

Oscillatory Criteria for a Class of Boundary Value Problem of Nonlinear Hyperbolic Equations

目的研究一类具有连续偏差变元的双曲偏泛函微分方程边值问题解的振动性．方法利用平均化方法，将多维边值问题解的振动性问题转化为常微分方程及其不等式的一维振动问题进行讨论．结果与结论推广了已有的一类具有离散偏差变元的双曲方程边值问题解的振动性的结果，得到了一类具有连续偏差变元的双曲偏泛函微分方程在两类不同边界条件下解的振动准则．

Comparative Studies between Picard’s and Taylor’s Methods of Numerical Solutions of First Ordinary Order Differential Equations Arising from Real-Life Problems

To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.

Existence of Entropy Solution for Degenerate Parabolic-Hyperbolic Problem Involving p(x)-Laplacian with Neumann Boundary Condition

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.

ON THE OSCILLATION OF SOLUTIONS OF HYPERBOLIC PARTIAL FUNCTIONAL DIFFERENTIAL EQUATIONS

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.

Hyperbolic Mild Slope Equations with Inclusion of Amplitude Dispersion Effect:Regular Waves

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.

Oscillation Theorem of Systems of Quasilinear Impulsive Delay Hyperbolic Equations

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.

OSCILLATION OF SOLUTIONS OF THE SYSTEMS OF QUASILINEAR HYPERBOLIC EQUATION UNDER NONLINEAR BOUNDARY CONDITION

Sufficient conditions are obtained for the oscillation of solutions of the systems of quasilinear hyperbolic differential equation with deviating arguments under nonlinear boundary condition.

Hyperbolic Mild Slope Equations with Inclusion of Amplitude Dispersion Effect:Random Waves

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 REDUCED FE FORMULATION BASED ON POD METHOD FOR HYPERBOLIC EQUATIONS

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.

New Numerical Scheme for Simulation of Hyperbolic Mild-Slope Equation

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.

OSCILLATION CRITERIA OF NEUTRAL TYPE IMPULSIVE HYPERBOLIC EQUATIONS

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.

A POD REDUCED-ORDER SPDMFE EXTRAPOLATING ALGORITHM FOR HYPERBOLIC EQUATIONS

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.