The Lie Theory of Two-Variable Hypergeometric Functions

We develop a Lie-algebraic method that associates with each of the 34 distinct second-order hypergeometric functions in two variables a canonical system of partial differential equations. The special functions arise by partial separation of variables in these simple systems. Some consequences are a demonstration that all such functions appear as solutions of the 4-variable wave equation and a classification of the possible imbeddings. In each case the functions are characterized by first- and second-order operators in the enveloping algebra of the conformal symmetry algebra for the wave equation. In some cases the 3-variable wave and heat equations and the 2-variable Helmholtz equation also arise. This intimate relationship between Horn functions and some fundamental equations of mathematical physics shows that these functions are more interesting than was previously recognized and permits use of the powerful tools of Lie theory and separation of variables to obtain properties of the functions.     

Construction of splitting schemes based on transition operator approximation

The stability analysis of approximate solutions to unsteady problems for partial differential equations is usually based on the use of the canonical form of operator-difference schemes. Another possibility widely used in the analysis of methods for solving Cauchy problems for systems of ordinary differential equations is associated with the estimation of the norm of the transition operator from the current time level to a new one. The stability of operator-difference schemes for a first-order model operator-differential equation is discussed. Primary attention is given to the construction of additive schemes (splitting schemes) based on approximations of the transition operator. Specifically, classical factorized schemes, componentwise splitting schemes, and regularized operator-difference schemes are related to the use of a certain multiplicative transition operator. Additive averaged operator-difference schemes are based on an additive representation of the transition operator. The construction of second-order splitting schemes in time is discussed. Inhomogeneous additive operator-difference schemes are constructed in which various types of transition operators are used for individual splitting operators.     

A model-coupling framework for nearshore waves,currents, sediment transport,and seabed morphology

This paper presents a framework for synchronously coupling wave, current, sediment transport, and seabed morphology for the accurate simulation of multi-physics coastal ocean processes. The governing equations, which represent models that are commonly adopted in practical simulations, are discretized using finite-difference methods. The resulting system is validated against analytical solutions. In order to test the performance of the proposed framework and the numerical methods, dam-break flow over a mobile-bed and evolution of a wave-driven sand dune are simulated. The interactions among waves, currents, and seabed morphology are illustrated.     

Detailed simulation of the pulsating detonation wave in the shock-attached frame

The paper is devoted to the numerical investigation of the stability of propagation of pulsating gas detonation waves. For various values of the mixture activation energy, detailed propagation patterns of the stable, weakly unstable, irregular, and strongly unstable detonation are obtained. The mathematical model is based on the Euler system of equations and the one-stage model of chemical reaction kinetics. The distinctive feature of the paper is the use of a specially developed computational algorithm of the second approximation order for simulating detonation wave in the shock-attached frame. In distinction from shock capturing schemes, the statement used in the paper is free of computational artifacts caused by the numerical smearing of the leading wave front. The key point of the computational algorithm is the solution of the equation for the evolution of the leading wave velocity using the second-order grid-characteristic method. The regimes of the pulsating detonation wave propagation thus obtained qualitatively match the computational data obtained in other studies and their numerical quality is superior when compared with known analytical solutions due to the use of a highly accurate computational algorithm.     

Vector domain decomposition schemes for parabolic equations

A new class of domain decomposition schemes for finding approximate solutions of timedependent problems for partial differential equations is proposed and studied. A boundary value problem for a second-order parabolic equation is used as a model problem. The general approach to the construction of domain decomposition schemes is based on partition of unity. Specifically, a vector problem is set up for solving problems in individual subdomains. Stability conditions for vector regionally additive schemes of first- and second-order accuracy are obtained.     

On a certain class of a nonlinear evolution partial differential equations

An existence result and a priori bound for the solution of a second-order nonlinear parabolic equation are established. Also a generalized tanh-function method is used for constructing exact travelling wave solutions for the nonlinear diffusion equation of Fisher type originated from the considered partial differential equation. And new multiple soliton solutions are obtained.     

On the convergence of solutions for a class of nonconformal FEM schemes for quasilinear elliptic equations

We consider versions of the nonconformal finite element method for the approximation to a second-order quasilinear elliptic equation in divergence form. For the construction of grid schemes, we use an approach used earlier for the nonstationary convection-diffusion equation and based on the Galerkin-Petrov limit approximation to the mixed statement of the original problem. The accuracy of solutions of nonconformal schemes with triangular linear finite elements is estimated in the absence of interior penalty terms, which are usually used in methods close to DG-methods for the stabilization of the scheme solution.     

有限深两层流中内孤立波的高阶解*

本文研究有限水深两层流中孤立波的三阶近似理论,并考虑了自由表面对孤立波的影响,运用坐标变形方法得到了三阶内孤立波的发展方程,求得波速的解析表达式。对方程进行了数值计算,得到了几种参数下三阶解曲线,指出自由表面对波型和波速的影响是二阶的。计算表明三阶解对一阶、二阶解有明显的改进,使其更加接近试验结果。     

带吸收边界条件的声波方程显式差分格式的稳定性分析

１．引言 在进行无界或半无界区域上各种波动方程的数值求解时,需引进入工边界以限制计算范围,在这些边界上应加相应的人工边界条件．这种边界条件应保证所求得的有界区域上的解很好地逼近原来无界区域上的解．对波动方程来说,就是在边界上人工反射应尽可能地小,使之对区域内部解的影响在允许的误差范围以内.因而它们被称为无反射边界条件或吸收边界条件．这种边界条件还应保证所形成的有界区域上的微分方程定解问题是适定的．这也是各种数值方法稳定的必要条件。 近二十多年来,已发展了声波方程的各种类型的吸收边界条件,其中以Ｃｌ…     

Characteristic-based schemes for dispersive waves I. The method of characteristics for smooth solutions

In order to embark on the development of numerical schemes for stiff problems, we have studied a model of relaxing heat flow. To isolate those errors unavoidably associated with discretization, a method of characteristics is developed, containing three free parameters depending on the stiffness ratio. It is shown that such “decoupled” schemes do not take into account the interaction between the wave families and hence result in incorrect wave speeds. We also demonstrate that schemes can differ by up to two orders of magnitude in their rms errors even while maintaining second-order accuracy. We show that no method of characteristics solution can be better than second-order accurate. Next, we develop “coupled” schemes which account for the interactions, and here we obtain two additional free parameters. We demonstrate how coupling of the two wave families can be introduced in simple ways and how the results are greatly enhanced by this coupling. Finally, numerical results for several decoupled and coupled schemes are presented, and we observe that dispersion relationships can be a very useful qualitative tool for analysis of numerical algorithms for dispersive waves. © 1993 John Wiley & Sons, Inc.     

Linear second order in time energy stable schemes for hydrodynamic models of binary mixtures based on a spatially pseudospectral approximation

We develop two linear, second order energy stable schemes for solving the governing system of partial differential equations of a hydrodynamic phase field model of binary fluid mixtures. We first apply the Fourier pseudo-spectral approximation to the partial differential equations in space to obtain a semi-discrete, time-dependent, ordinary differential and algebraic equation (DAE) system, which preserves the energy dissipation law at the semi-discrete level. Then, we discretize the DAE system by the Crank-Nicolson (CN) and the second-order backward differentiation/extrapolation (BDF/EP) method in time, respectively, to obtain two fully discrete systems. We show that the CN method preserves the energy dissipation law while the BDF/EP method does not preserve it exactly but respects the energy dissipation property of the hydrodynamic model. The two new fully discrete schemes are linear, unconditional stable, second order accurate in time and high order in space, and uniquely solvable as linear systems. Numerical examples are presented to show the convergence property as well as the efficiency and accuracy of the new schemes in simulating mixing dynamics of binary polymeric solutions.     

二阶Camassa-Holm方程行波解

对二阶Camassa-Holm方程行波解的情况进行了讨论.利用解的唯一性,得到了如下结论:二阶CH方程的行波解唯一存在,但不具有u(x,t)=kem(x-ct)形式.还为二阶CH方程行波解的研究提供了一种新途径和方法.     

A novel convenient finite difference method for shallow water waves derived by fifth-order Kortweg and De-Vries-type equation

In this research, two new finite difference schemes are derived and presented for estimating a solution to the fifth-order Kortweg and De-Vries equation. The global conservation law on any time–space regions which yields a three-level linear-implicit algorithm with a diagonal system is exactly preserved by the intended finite difference method. Theoretically verified and numerically proved, the created schemes are unconditionally stable and have the second-order accuracy both in time and space. The obtained results guarantee that the novel idea offers a new aspect to analyze the wave behavior.     

Constructing exact solutions to discrete systems with the trial function method

Based on the homogenous balance method and the trial function method, several trial function methods composed of exponential functions are proposed and applied to nonlinear discrete systems. With the.help of symbolic computation system, the new exact solitary wave solutions to discrete nonlinear mKdV lattice equation, discrete nonlinear （2 ＋ 1） dimensional Toda lattice equation, Ablowitz-Ladik-lattice system are constructed.The method is of significance to seek exact solitary wave solutions to other nonlinear discrete systems.     

A CELL-CENTERED ALE METHOD WITH HLLC-2D RIEMANN SOLVER IN 2D CYLINDRICAL GEOMETRY

This paper presents a second-order direct arbitrary Lagrangian Eulerian (ALE) method for compressible flow in two-dimensional cylindrical geometry.This algorithm has half-face fluxes and a nodal velocity solver,which can ensure the compatibility between edge fluxes and the nodal flow intrinsically.In two-dimensional cylindrical geometry,the control vol-ume scheme and the area-weighted scheme are used respectively,which are distinguished by the discretizations for the source term in the momentum equation.The two-dimensional second-order extensions of these schemes are constructed by employing the monotone up-wind scheme of conservation law (MUSCL) on unstructured meshes.Numerical results are provided to assess the robustness and accuracy of these new schemes.     

Analytical solution of a spatially variable coefficient advection–diffusion equation in up to three dimensions

Analytical solutions are provided for the two- and three-dimensional advection–diffusion equation with spatially variable velocity and diffusion coefficients. We assume that the velocity component is proportional to the distance and that the diffusion coefficient is proportional to the square of the corresponding velocity component. There is a simple transformation which reduces the spatially variable equation to a constant coefficient problem for which there are available a large number of known analytical solutions for general initial and boundary conditions. These solutions are also solutions to the spatially variable advection–diffusion equation. The special form of the spatial coefficients has practical relevance and for divergent free flow represent corner or straining flow. Unlike many other analytical solutions, we use the transformation to obtain solutions of the spatially variable coefficient advection–diffusion equation in two and three dimensions. The analytical solutions, which are simple to evaluate, can be used to validate numerical models for solving the advection–diffusion equation with spatially variable coefficients. For numerical schemes which cannot handle flow stagnation points, we provide analytical solution to the spatially variable coefficient advection–diffusion equation for two-dimensional corner flow which contains an impermeable flow boundary. The impermeable flow boundary coincides with a streamline along which the fluid velocity is finite but the concentration vanishes. This example is useful for validating numerical schemes designed to predict transport around a curved boundary.     

A new high-order energy-preserving scheme for the modified Korteweg-de Vries equation

In this paper, a new high-order energy-preserving scheme is proposed for the modified Korteweg-de Vries equation. The proposed scheme is constructed by using the Hamiltonian boundary value methods in time, and Fourier pseudospectral method in space. Exploiting this method, we get second-order and fourth-order energy-preserving integrators. The proposed schemes not only have high accuracy, but also exactly conserve the total mass and energy in the discrete level. Finally, single solitary wave and the interaction of two solitary waves are presented to illustrate the effectiveness of the proposed schemes.     

Flux-splitting schemes for parabolic equations with mixed derivatives

Difference schemes of required quality are often difficult to construct as applied to boundary value problems for parabolic equations with mixed derivatives. Specifically, difficulties arise in the design of monotone difference schemes and unconditionally stable locally one-dimensional splitting schemes. In parabolic problems, certain opportunities are offered by restating the problem in question so that the quantities to be determined are fluxes (directional derivatives). The original problem is then rewritten as a boundary value one for a system of equations in flux variables. Weighted schemes for parabolic equations in flux coordinates are examined. Unconditionally stable locally one-dimensional flux schemes that are first- and second-order accurate in time are constructed for a parabolic equation without mixed derivatives. A feature of systems in flux variables for equations with mixed derivatives is that the terms with time derivatives are coupled with each other.     

Splitting positive definite mixed element method for viscoelasticity wave equation

A splitting positive definite mixed finite element method is proposed for second-order viscoelasticity wave equation. The proposed procedure can be split into three independent symmetric positive definite integro-differential sub-system and does not need to solve a coupled system of equations. Error estimates are derived for both semidiscrete and fully discrete schemes. The existence and uniqueness for semidiscrete scheme are proved. Finally, a numerical example is provided to illustrate the efficiency of the method.     

Explicit Peakon solutions to a family of wave-breaking equations

The singular traveling wave solutions of a general 4-parameter family equation which unifies the Camass-Holm equation, the Degasperis-Procesi equation and the Novikov equation are investigated in this paper. At first, we obtain the explicit peakon solutions for one of its specific case that $a=(p+2)c$, $b=(p+1)c$ and $c=1$, which is referred to a generalized Camassa-Holm-Novikov (CHN) equation, by reducing it to a second-order ordinary differential equation (ODE) and solving its associated first-order integrable ODE. By observing the characteristics of peakon solutions to the CHN equation, we construct the peakon solutions for the general 4-parameter breaking wave equation. It reveals that singularities of the peakon solutions come up only when the solutions attain singular points of the equation, which might be a universal principal for all singular traveling wave solutions for wave breaking equations.