STABILITY OF VISCOUS SHOCK WAVES FOR THE ONE-DIMENSIONAL COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH DENSITY-DEPENDENT VISCOSITY

We study the large-time behavior toward viscous shock waves to the Cauchy problem of the one-dimensional compressible isentropic Navier-Stokes equations with density- dependent viscosity. The nonlinear stability of the viscous shock waves is shown for certain class of large initial perturbation with integral zero which can allow the initial density to have large oscillation. Our analysis relies upon the technique developed by Kanel~ and the continuation argument.

NONLINEAR STABILITY OF VISCOUS SHOCK WAVES FOR ONE-DIMENSIONAL NONISENTROPIC COMPRESSIBLE NAVIER–STOKES EQUATIONS WITH A CLASS OF LARGE INITIAL PERTURBATION

We study the nonlinear stability of viscous shock waves for the Cauchy problem of one-dimensional nonisentropic compressible Navier–Stokes equations for a viscous and heat conducting ideal polytropic gas. The viscous shock waves are shown to be time asymptotically stable under large initial perturbation with no restriction on the range of the adiabatic exponent provided that the strengths of the viscous shock waves are assumed to be sufficiently small.The proofs are based on the nonlinear energy estimates and the crucial step is to obtain the positive lower and upper bounds of the density and the temperature which are uniformly in time and space.

POINTWISE ESTIMATES OF SOLUTIONS FOR THE NONLINEAR VISCOUS WAVE EQUATION IN EVEN DIMENSIONS

In this article,we study the pointwise estimates of solutions to the nonlinear viscous wave equation in even dimensions(n≥4).We use the Green’s function method.Our approach is on the basis of the detailed analysis of the Green’s function of the linearized system.We show that the decay rates of the solution for the same problem are different in even dimensions and odd dimensions.It is shown that the solution exhibits a generalized Huygens principle.

THE STABILITY OF THE DELTA WAVE TO PRESSURELESS EULER EQUATIONS WITH VISCOUS AND FLUX PERTURBATIONS

This paper is concerned with the pressureless Euler equations with viscous and flux perturbations.The existence of Riemann solutions to the pressureless Euler equations with viscous and flux perturbations is obtained.We show the stability of the delta wave of the pressureless Euler equations to the perturbations;that is,the limit solution of the pressureless Euler equations with viscous and flux perturbations is the delta wave solution of the pressureless Euler equations as the viscous and flux perturbation simultaneously vanish in the case u_(-)> u_(+).

An ADI Finite Volume Element Method for a Viscous Wave Equation with Variable Coefficients

Based on rectangular partition and bilinear interpolation,we construct an alternating-direction implicit(ADI)finite volume element method,which combined the merits of finite volume element method and alternating direction implicit method to solve a viscous wave equation with variable coefficients.This paper presents a general procedure to construct the alternating-direction implicit finite volume element method and gives computational schemes.Optimal error estimate in L2 norm is obtained for the schemes.Compared with the finite volume element method of the same convergence order,our method is more effective in terms of running time with the increasing of the computing scale.Numerical experiments are presented to show the efficiency of our method and numerical results are provided to support our theoretical analysis.

Deformed two-dimensional rogue waves in the (2+1)-dimensional Korteweg–de Vries equation

Within the(2+1)-dimensional Korteweg–de Vries equation framework,new bilinear B¨acklund transformation and Lax pair are presented based on the binary Bell polynomials and gauge transformation.By introducing an arbitrary functionφ(y),a family of deformed soliton and deformed breather solutions are presented with the improved Hirota’s bilinear method.By choosing the appropriate parameters,their interesting dynamic behaviors are shown in three-dimensional plots.Furthermore,novel rational solutions are generated by taking the limit of the obtained solitons.Additionally,twodimensional(2D)rogue waves(localized in both space and time)on the soliton plane are presented,we refer to them as deformed 2D rogue waves.The obtained deformed 2D rogue waves can be viewed as a 2D analog of the Peregrine soliton on soliton plane,and its evolution process is analyzed in detail.The deformed 2D rogue wave solutions are constructed successfully,which are closely related to the arbitrary functionφ(y).This new idea is also applicable to other nonlinear systems.

TWO-DIMENSIONAL RIEMANN PROBLEMS:FROM SCALAR CONSERVATION LAWS TO COMPRESSIBLE EULER EQUATIONS

In this paper we survey the authors＇ and related work on two-dimensional Riemann problems for hyperbolic conservation laws, mainly those related to the compressible Euler equations in gas dynamics. It contains four sections： 1. Historical review. 2. Scalar conservation laws. 3. Euler equations. 4. Simplified models.

STABILITY OF VISCOUS CONTACT WAVE FOR COMPRESSIBLE NAVIER-STOKES SYSTEM OF GENERAL GAS WITH FREE BOUNDARY

In this paper, we study the large time behavior of solutions to the nonisentropic Navier-Stokes equations of general gas, where polytropic gas is included as a special case, with a free boundary. First we construct a viscous contact wave which approximates to the contact discontinuity, which is a basic wave pattern of compressible Euler equation, in finite time as the heat conductivity tends to zero. Then we prove the viscous contact wave is asymptotic stable if the initial perturbations and the strength of the contact wave are small. This generalizes our previous result [6] which is only for polytropic gas.

DIFFRACTION OF A SOLITARY WAVE BY A THIN WEDGE PART 2. EFFECTS OF VISCOUS BOUNDARY LAYERS

Based on the study on the Mach reflection of a solitary wave in [3] , we continue to investi- gate effects of the boundary layers on the bottom and the vertical side wall. By using matched asymptotic methods, the two-dimensional KdV equation is modified to account for effects of viscosity. Numerical simulation of the problem shows that the effects of side wall are important while the effects of the bottom can be neglected. The results including the side wall's effects agree satisfactorily with those of Melville's experiments. Finally, we establish the simplified concept of the side wall effect and conclude that it repre- sents the physical reason for the discrepancy between the experiments and the previous calculations based on the inviscid fluid flow theory.

Centered simple waves for the two-dimensional pseudo-steady isothermal ?ow around a convex corner

The two-dimensional(2D) pseudo-steady isothermal flow, which is isentropic and irrotational, around a convex corner is studied. The self-similar solutions for the supersonic flow around the convex corner are constructed, where the properties of the centered simple wave are used for the 2D isentropic irrotational pseudo-steady Euler equations. The geometric procedures of the center simple waves are given. It is proven that the supersonic flow turns the convex corner by an incomplete centered expansion wave or an incomplete centered compression wave, depending on the conditions of the downstream state.

粘性波动方程解的逐点估计

本文中研究二维空间带粘性项非线性波动方程解的逐点估计。在经典解整体存在的情况下,解可由格林函数表示,通过对格林函数的详细分析,以及它与非线性项的相互作用,得到解的逐点衰减。

Two-dimensional Wave Equations with Fractal Boundaries

This paper focuses on two cases of two-dimensional wave equations with fractal boundaries. The first case is the equation with classical derivative. The formal solution is obtained. And a definition of the solution is given. Then we prove that under certain conditions, the solution is a kind of fractal function, which is continuous, differentiable nowhere in its domain. Next, for specific given initial position and 3 different initial velocities, the graphs of solutions are sketched. By computing the box dimensions of boundaries of cross-sections for solution surfaces, we evaluate the range of box dimension of the vibrating membrane. The second case is the equation with p-type derivative. The corresponding solution is shown and numerical example is given.

A compact scheme for two-dimensional nonlinear time fractional wave equations

Based on the equivalent integro-differential form of the considered problem, a numerical approach to solving the two-dimensional nonlinear time fractional wave equations(NTFWEs) is considered in this paper. To this end, an alternating direction implicit(ADI) numerical scheme is derived. The scheme is established by combining the secondorder convolution quadrature formula and Crank–Nicolson technique in time and afourth-order difference approach in space. The convergence and unconditional stability of the proposed compact ADI scheme are strictly discussed after a concise solvabilityanalysis. A numerical example is shown to demonstrate the theoretical analysis.

Numerical studies of high-dimensional nonlinear viscous and nonviscous wave equations by using finite difference methods

The numerical solutions of two-dimensional(2D)and three-dimensional(3D)nonlinear viscous and nonviscous wave equations via the unified alternating direction implicit(ADI)finite difference methods(FDMs)are obtained in this paper.By making use of the discrete energy method,it is proven that their numerical solutions converge to exact solutions with an order of two in both time and space with respect to H^(1)-norm.Numerical results confirm that they are relatively accurate and high-resolution,and more successfully simulate the conservation of the energy for nonviscous equations,and the dissipation of the energy for viscous equation.

THE INTERACTIONS BETWEEN WAVE-CURRENTS AND OFFSHORE STRUCTURES WITH CONSIDERATION OF FLUID VISCOSITY

Study of the how held around the large scale offshore structures under the action of waves and viscous currents is of primary importance for the scouring estimation and protection in the vicinity of the structures. But very little has been known in its mechanism when the viscous effects is taken into consideration. As a part of the efforts to tackle the problem, a numerical model is presented for the simulation of the how held around a fixed vertical truncated circular cylinder subjected to waves and viscous currents based on the depth-averaged Reynolds equations and depth-averaged k-epsilon turbulence model. Finite difference method with a suitable iteration defect correct method and an artificial open boundary condition are adopted in the numerical process. Numerical results presented relate to the interactions of a pure incident viscous current with Reynolds number Re = 10(5), a pure incident regular sinusoidal wave, and the coexisting of viscous current and wave with a circular cylinder, respectively. Flow fields associated with the hydrodynamic coefficients of the fixed cylinder, as well as corresponding free surface profiles and wave amplitudes, are discussed. The present method is found to be relatively straightforward, computationally effective and numerically stable for treating the problem of interactions among waves, viscous currents and bodies.

Damping Solitary Wave in a Three-Dimensional Rectangular Geometry Plasma

The solitary waves of a viscous plasma confined in a cuboid under the three types of boundary condition are theoretically investigated in the present paper.By introducing a threedimensional rectangular geometry and employing the reductive perturbation theory,a quasi-Kd V equation is derived in the viscous plasma and a damping solitary wave is obtained.It is found that the damping rate increases as the viscosity coefficient increases,or increases as the length and width of the rectangle decrease,for all kinds of boundary condition.Nevertheless,the magnitude of the damping rate is dominated by the types of boundary condition.We thus observe the existence of a damping solitary wave from the fact that its amplitude disappears rapidly for a → 0and b → 0,or ν→ ＋∞.

基于紧致交错差分格式的二维声波及黏滞声波方程数值模拟

针对常规网格差分难以适用于地震波场数值模拟中复杂介质的问题,首次将紧致交错有限差分格式应用于黏滞声波方程的数值模拟研究并同声波方程的数值模拟进行了模拟精度、频散关系和稳定性分析等方面的比较。理论研究结果表明:当差分精度相同时,紧致交错网格所需节点数要少于常规的中心差分和交错差分格式,计算效率更高;同常规的交错差分与中心差分格式相比,紧致差分的截断误差更小,数值频散也更低,能够适用于粗网格计算;差分精度相同情况下时进行数值模拟,紧致交错格式所需要的时间网格更小,稳定性条件也更为严格;紧致交错差分格式在完全匹配层(perfectly matched layer, PML)条件下,能够对边界反射进行有效吸收。最后,对均匀、水平层状介质以及Marmousi模型进行了黏滞声波方程的数值模拟和波场特征分析,实验结果证明了该方法对于复杂介质的数值模拟的适应性和有效性,并具有较高的模拟精度及计算效率。

基于POD有限元法的粘滞声波方程震源波场重构

叠前逆时偏移是地震勘探中一种流行的地下结构成像方法,其成像条件需要同时刻的震源波场值与检波器波场值.这在实际计算中就需要把正演模拟的所有时刻的震源波场数据全部存储下来,存储需求大.虽然震源波场重构技术可以降低对于波场数据的存储需求,但会引入额外的计算复杂度.为解决这个问题,本文提出了POD有限元法,并将其应用于粘滞震源波场重构.这里的本征正交分解(Proper Orthogonal Decomposition,POD)方法是一种降维方法,能够在降低数据量的同时提供足够的计算精度.数值算例显示,该方法比传统的有限元法更节省存储空间,能够加快重构速度.

粘性依赖于密度的一维等熵可压缩Navier-Stokes方程组粘性激波的非线性稳定性

该文主要研究粘性系数依赖于密度的一维等熵可压缩Navier-Stokes方程组Cauchy问题整体解的大时间渐近行为,主要研究目的是改进文献[7]的结果至γ>1,κ≥0.注意到γ=2,κ=1时,一维等熵可压缩Navier-Stokes方程组对应于Saint-Venant浅水波方程组,该方程组描述了地表浅水运动的规律,在物理学和海洋学中有重要的应用^([1,4,6])。注意到文献^([7])中通过利用Kanel的方法^([19])来推导比容的一致上下界估计,在得出比容的上界时,该方法要求κ<1/2.对该文所研究的问题而言,需要首先利用Kanel’的方法^([19])来推导比容的一致上下界估计.为了扩大κ的取值范围,还需要对比容的上下界作更为精细的能量估计.在得出比容的一致上下界估计之后,可通过精心设计的连续性技巧,将Navier-Stokes方程组的局部解一步步延拓为整体解,并得到对应的大时间渐近行为.

Viscous Shock Wave to an Inflow Problem for Compressible Viscous Gas with Large Density Oscillations

This paper is concerned with the inflow problem for one-dimensional compressible Navier-Stokes equations. For such a problem, Huang, Matsumura, and Shi showed in [4] that there exists viscous shock wave solution to the inflow problem and both the boundary layer solution, the viscous shock wave, and their superposition are time-asymptotically nonlinear stable provided that both the initial perturbation and the boundary velocity are assumed to be sufficiently small. The main purpose of this paper is to show that similar stability results still hold for a class of large initial perturbation which can allow the initial density to have large oscillations. The proofs are given by an elementary energy method and our main idea is to use the smallness of the strength of the viscous shock wave and the boundary velocity to control the possible growth of the solutions induced by the nonlinearity of the compressible Navier-Stokes equations and the inflow boundary condition.The key point in our analysis is to deduce the desired uniform positive lower and upper bounds on the density.