浅水长波近似方程组的拟小波解

以浅水长波近似方程组为例,提出了拟小波方法求解(1 1)维非线性偏微分方程组数值解,该方程用拟小波离散格式离散空间导数,得到关于时间的常微分方程组,用四阶Runge-K utta方法离散时间导数,并将其拟小波解与解析解进行比较和验证.     

一种改进的截断展开法求非线性发展方程的精确解

给出了一种改进的截断展开法,利用此方法借助于计算机符号计算求得了Burgers方程和浅水长波近似方程组的精确解,其中包括孤子解,并讨论其具体应用.改进后的方法与以前的方法相比能得到方程的更多形式的精确解.所给出的改进的截断展开法也可以用来研究其它非线性发展方程的孤子解,是求非线性发展方程精确解的一种有效的直接方法.     

解高维广义对称正则长波方程的Fourier谱方法

１引言对称正则长波方程（ＳＲＬＷＥ）是正则化长波方程（ＲＬＷＥ）的一种对称叙述［１］用于描述弱非线性作用下空间变换的离子声波传播．［１］得到了方程组（１．１）的双曲正割平方孤立波解、四个不变量和数值结果、明显地,从（１．１）中消去ρ,得到一类正则长波方程（ＲＬＷＥ）代替（１．２）中第三项、第四项对ｔ的导数为对ｘ的导数,得到Ｂｏｕｓｓｉｎｅｓｑ方程．［２］对一类广义对称正则长波方程组提出了谱方法,证明了古典光滑解的存在性和唯一性,建立了近似解的收敛性和误差估计。［３］研究了高维对称正则长波方程整体…     

赤道β-平面近似的无量纲线性浅水方程组的边值问题

应用分层理论所提供的方法,研究了赤道β-平面近似的无量纲线性浅水方程组的边值问题的适定性以及相应的识别方法.对于在超平面{x=f(y,t)}和{y=g(x,t)}的边值问题,分别给出了它们局部解的解空间构造.     

浅水方程组的数值方法

构造了浅水方程组的二阶精度的TVD格式。格式由简单的TVD Runge-Kutta型时间离散和有坡度限制的空间对称离散格式组成。数值耗散项用局部棱柱化河道流的特征变量构造。格式的主要优点是能够计算天然河道中浅水方程组的弱解并且构造简单。格式能够求出天然河道或非平底部渠道中的精确静水解。给出了渠道溃坝问题数值解与解析解的比较,验证格式精度高。实际天然河道型梯级水库溃坝的数值实验表明格式稳定,适应性强。     

长水波近似方程组的动力学行为和辛几何算法

长水波近似方程组作为一类重要的非线性方程有着许多广泛的应用前景,特别是在浅水非线性色散波的研究中具有重要意义.给出了长水波近似方程组的动力学行为,并基于Hamilton空间体系的多辛理论研究了长水波近似方程组的数值解法,讨论了利用Preissmann方法构造离散多辛格式的途径,并构造了一种典型的半隐式的多辛格式,格式满足多辛守恒律.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

非线性方程组新的再生核方法

给出了二阶非线性方程组边值问题新的数值方法.该方法基于再生核和最小二乘法.首先建立再生核直积空间,接着证明了近似解的一致收敛性,避免了耗时的施密特正交化过程.数值算例展现出该算法简单、有效.     

长水波近似方程组的新精确解

依据齐次平衡法的思想 ,首先提出了求非线性发展方程精确解的新思路 ,这种方法通过改变待定函数的次序 ,优势是使求解的复杂计算得到简化 .应用本文的思路 ,可得到某些非线性偏微分方程的新解 .其次我们给出了长水波近似方程组的一些新精确解 ,其中包括椭圆周期解 ,我们推广了有关长波近似方程的已有结果 .     

一种解非线性方程组的不精确牛顿——兰索斯方法

本文提出一种不完全线搜索技术的不精确牛顿—克雷洛夫(Newton-Krylov)子空间方法解对称非线性方程组,其中克雷洛夫子空间方法采用的是兰索斯(Lanczos)类分解技术.迭代方向是通过使用兰索斯方法近似求解非线性方程组的牛顿方程获得的.在合理的假设条件下,分析了算法的全局收敛性和局部超线性收敛速率.最后,数值结果显示了该算法的有效性.     

一类一维的辐射流体力学方程组激波的存在性

该文主要讨论一维空间中一类辐射流体力学方程组的激波. 由Rankine-Hugoniot条件及熵条件得此问题可表述为关于辐射流体力学方程组带自由边界的初边值问题. 首先通过变量代换, 将其自由边界转换为固定边界, 然后研究关于此非线性方程组的一个初边值问题解的存在唯一性. 为此先构造了此问题的一个近似解, 然后分别通过Picard迭代与Newton迭代对此非线性问题构造近似解序列. 通过一系列估计与紧性理论得到此近似解序列的收敛性, 其极限即为原辐射热力学方程组的一个激波.     

Homogenization Method in the Problem of Long Wave Propagation from a Localized Source in a Basin over an Uneven Bottom

In the framework of the linearized shallow water equations, the homogenization method for wave type equations with rapidly oscillating coefficients that generally cannot be represented as periodic functions of the fast variables is applied to the Cauchy problem for the wave equation describing the evolution of the free surface elevation for long waves propagating in a basin over an uneven bottom. Under certain conditions on the function describing the basin depth, we prove that the solution of the homogenized equation asymptotically approximates the solution of the original equation. Model homogenized wave equations are constructed for several examples of one-dimensional sections of the real ocean bottom profile, and their numerical and asymptotic solutions are compared with numerical solutions of the original equations.     

An Efficient Method for Computing Hyperbolic Systems with Geometrical Source Terms Having Concentrations

We propose a simple numerical method for calculating both unsteady and steady state solution of hyperbolic system with geometrical source terms having concentrations. Physical problems under consideration include the shallow water equations with topography,and the quasi one-dimensional nozzle flows. We use the interface value, rather than the cell-averages, for the source terms, which results in a well-balanced scheme that can capture the steady state solution with a remarkable accuracy. This method approximates the source terms via the numerical fluxes produced by an (approximate) Riemann solver for the homogeneous hyperbolic systems with slight additional computation complexity using Newton‘s iterations and numerical integrations. This method solves well the subor super-critical flows, and with a transonic fix, also handles well the transonic flows over the concentration. Numerical examples provide strong evidence on the effectiveness of this new method for both unsteady and steady state calculations.     

Analysis of the Nonlinear Shallow Water Equations Over Nonplanar Topography

The role played by the beach bottom profile on coastal inundation phenomena is analyzed here by means of approximate analytical solutions of the nonlinear shallow water equations (NSWEs) over uneven bottoms. These are obtained by only using the assumptions of small waves at the seaward boundary and small topographic forcing. Our work, built on the Carrier and Greenspan [ 1 ] hodographic transformation and on the solution of the boundary value problem (BVP) for the NSWEs proposed by Antuono and Brocchini [ 2 ], focuses on the propagation of nonlinear non-breaking waves over quasi-planar beaches. Since the terms associated with the perturbed bottom only appear in the second-order perturbed solutions, the breaking conditions for the planar-beach bathymetry also predict well the breaking occurring on the nonplanar beaches analyzed here. The most important results, concerning the shoreline position and the near-shoreline velocity, are given for both pulse-like and periodic input waves propagating over two types of nonplanar bathymetries. The solution proposed here is a fundamental benchmark for any numerical and theoretical analyzes concerned with estimates of wave run-up on beaches of complex shape.     

Numerical solutions of two-dimensional nonlinear integral analysis

In this paper, the approximate solutions for two different type of two-dimensional nonlinear integral equations: two-dimensional nonlinear Volterra-Fredholm integral equations and the nonlinear mixed Volterra-Fredholm integral equations are obtained using the Laguerre wavelet method. To do this, these two-dimensional nonlinear integral equations are transformed into a system of nonlinear algebraic equations in matrix form. By solving these systems, unknown coefficients are obtained. Also, some theorems are proved for convergence analysis.Some numerical examples are presented and results are compared with the analytical solution to demonstrate the validity and applicability of the proposed method.     

An efficient splitting technique for two‐layer shallow‐water model

We consider the numerical approximation of the weak solutions of the two‐layer shallow‐water equations. The model under consideration is made of two usual one‐layer shallow‐water model coupled by nonconservative products. Because of the nonconservative products of the system, which couple both one‐layer shallow‐water subsystems, the usual numerical methods have to consider the full model. Of course, uncoupled numerical techniques, just involving finite volume schemes for the basic shallow‐water equations, are very attractive since they are very easy to implement and they are costless. Recently, a stable layer splitting technique was introduced [Bouchut and Morales de Luna, M2AN Math Model Numer Anal 42 (2008), 683–698]. In the same spirit, we exhibit new splitting technique, which is proved to be well balanced and non‐negative preserving. The main benefit issuing from the here derived uncoupled method is the ability to correctly approximate the solution of very severe benchmarks. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1396–1423, 2015     

Interaction of Shallow Water Waves

In this paper, we consider the Riemann problem and interaction of elementary waves for a nonlinear hyperbolic system of conservation laws that arises in shallow water theory. This class of equations includes as a special case the equations of classical shallow water equations. We study the bore and dilatation waves and their properties, and show the existence and uniqueness of the solution to the Riemann problem. Towards the end, we discuss numerical results for different initial data along with all possible interactions of elementary waves. It is noticed that in contrast to the p -system, the Riemann problem is solvable for arbitrary initial data, and its solution does not contain vacuum state.     

The C 1 wavelet method for numerical solutions of the Neumann problem

In this article we use the C ¹ wavelet bases on Powell-Sabin triangulations to approximate the solution of the Neumann problem for partial differential equations. The C ¹ wavelet bases are stable and have explicit expressions on a three-direction mesh. Consequently, we can approximate the solution of the Neumann problem accurately and stably. The convergence and error estimates of the numerical solutions are given. The computational results of a numerical example show that our wavelet method is well suitable to the Neumann boundary problem.