An exact Riemann solver for detonation products

Numerical methods based upon the Riemann Problem are considered for solving the general initial-value problem for the Euler equations applied to real gases. Most of such methods use an approximate solution of the Riemann problem when real gases are involved. These approximate Riemann solvers do not yield always a good resolution of the flow field, especially for contact surfaces and expansion waves. Moreover, approximate Riemann solvers cannot produce exact solutions for the boundary points. In order to overcome these shortcomings, an exact solution of the Riemann problem is developed, valid for real gases. The method is applied to detonation products obeying a 5th order virial equation of state, in the shock-tube test case. Comparisons between our solver, as implemented in Random Choice Method, and finite difference methods, which do not employ a Riemann solver, are given.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Delta Shock Wave Solution of the Riemann Problem for the Non-homogeneous Modified Chaplygin Gasdynamics

The motivation of the present study is to derive the solution of the Riemann problem for modified Chaplygin gas equations in the presence of constant external force. The analysis leads to the fact that in some special circumstances delta shock appears in the solution of the Riemann problem. Also, the Rankine–Hugoniot relations for delta shock wave which are utilized to determine the strength, position and propagation speed of the delta shocks have been derived. Delta shock wave solution to the Riemann problem for the modified Chaplygin gas equation is obtained. It is found that the external force term, appearing in the governing equations, influences the Riemann solution for the modified Chaplygin gas equation.

     

Explicit solution to the exact Riemann problem and application in nonlinear shallow‐water equations

The Riemann solver is the fundamental building block in the Godunov‐type formulation of many nonlinear fluid‐flow problems involving discontinuities. While existing solvers are obtained either iteratively or through approximations of the Riemann problem, this paper reports an explicit analytical solution to the exact Riemann problem. The present approach uses the homotopy analysis method to solve the nonlinear algebraic equations resulting from the Riemann problem. A deformation equation defines a continuous variation from an initial approximation to the exact solution through an embedding parameter. A Taylor series expansion of the exact solution about the embedding parameter provides a series solution in recursive form with the initial approximation as the zeroth‐order term. For the nonlinear shallow‐water equations, a sensitivity analysis shows fast convergence of the series solution and the first three terms provide highly accurate results. The proposed Riemann solver is implemented in an existing finite‐volume model with a Godunov‐type scheme. The model correctly describes the formation of shocks and rarefaction fans for both one and two‐dimensional dam‐break problems, thereby verifying the proposed Riemann solver for general implementation. Copyright © 2007 John Wiley & Sons, Ltd.     

A Riemann solver and upwind methods for a two‐phase flow model in non‐conservative form

We present a theoretical solution for the Riemann problem for the five‐equation two‐phase non‐conservative model of Saurel and Abgrall. This solution is then utilized in the construction of upwind non‐conservative methods to solve the general initial‐boundary value problem for the two‐phase flow model in non‐conservative form. The basic upwind scheme constructed is the non‐conservative analogue of the Godunov first‐order upwind method. Second‐order methods in space and time are then constructed via the MUSCL and ADER approaches. The methods are systematically assessed via a series of test problems with theoretical solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

The Riemann problem for non-ideal isentropic compressible two phase flows

We consider the Riemann problem for a five-equation, two-pressure (5E2P) model of non-ideal isentropic compressible gas–liquid two-phase flows. This system is more complex due to the extended thermodynamics model for van der Waals gases, that is, typical real gases for gas phase and Tait׳s equation of state for liquid phase. The overall model is strictly hyperbolic and non-conservative form. We investigate the structure of Riemann problem and construct the solution for it. To construct solution of Riemann problem approximately assuming that all waves corresponding to the genuinely non-linear characteristic fields are rarefaction and then we discuss their properties. Lastly, we discuss numerical examples and study the solution influenced by the van der Waals excluded volume.     

The Riemann problem for a hyperbolic model of two-phase flow in conservative form

This article is to continue the present author's work (International Journal of Computational Fluid Dynamics (2009) 23 (9), 623–641) on studying the structure of solutions of the Riemann problem for a system of three conservation laws governing two-phase flows. While existing solutions are limited and found quite recently for the Baer and Nunziato equations, this article presents the first instance of an exact solution of the Riemann problem for two-phase flow in gas–liquid mixture. To demonstrate the structure of the solution, we use a hyperbolic conservative model with mechanical equilibrium and without velocity equilibrium. The Riemann problem solution for the model equations comprises a set of elementary waves, rarefaction and discontinuous waves of various types. In particular, such a solution treats both the wave structure and the intermediate states of the two-phase gas–liquid mixture. The resulting exact Riemann solver is fully non-linear, direct and complete. On this basis then, we use locally the exact Riemann solver for the two-phase flow in gas–liquid mixture within the framework of finite volume upwind Godunov methods. In order to demonstrate the effectiveness and accuracy of the proposed solver, we consider a series of test problems selected from the open literature and compare the exact and numerical results by using upwind Godunov methods, showing excellent oscillation-free results in two-phase fluid flow problems.     

Analytic Solution of an Inhomogeneous Kinetic Equation in the Problem of Second-Order Thermal Creep

An analytic solution of the problem of second-order thermal creep is obtained. A method for solving the half-space boundary value problem for an inhomogeneous linearized kinetic BGK equation forms the basis of the solution. The general solution of the input equation is constructed in the form of an expansion of the corresponding characteristic equation in terms of the eigenfunctions. Substitution of the solution in the boundary conditions leads to a Riemann boundary value problem. The unknown thermal creep velocity is found from the condition of solvability of the boundary value problem. The numerical analysis performed confirms the existence of negative thermophoresis (in the direction of the temperature gradient) for high-conductivity aerosol particles at low Knudsen numbers.     

Exact solutions for some non-conservative hyperbolic systems

In this paper, we study initial value problem for some non-conservative hyperbolic systems of partial differential equations of first order. The first one is the Riemann problem for a model in elastodynamics and the second one the initial value problem for a system which is a generalization of the Hopf equation. The non-conservative products which appear in the equations do not make sense in the classical theory of distributions and are understood in the sense of Volpert (Math. USSR Sb. 2 (1967) 225). Following Lax (Comm. Pure Appl. Math. 10 (1957) 537) and Dal Maso et al. (J. Math. Pures Appl. 74 (1995) 483), we give an explicit solution for the Riemann problem for the elastodynamics equation. The coupled Hopf equation is studied using a generalization of the method of Hopf (Comm. Pure Appl. Math. 3 (1950) 201).     

On Goursat and Dirichlet problems for one equation of the third order

We study the problem of the unique solvability of Goursat and Dirichlet problems for one partial differential equation of the third order. We construct a Riemann function for a linear third-order equation with a hyperbolic operator in the principal part, study some properties of the Riemann function, and then use them to prove theorems on the existence and uniqueness of a solution of the problems indicated. Translated from Neliniini Kolyvannya, Vol. 11, No. 3, pp. 305–315, July–September, 2008.     

The reactive Riemann problem for thermally perfect gases at all combustion regimes

In this work we analyze the reactive Riemann problem for thermally perfect gases in the deflagration or detonation regimes. We restrict our attention to the case of one irreversible infinitely fast chemical reaction; we also suppose that, in the initial condition, one state (for instance the left one) is burnt and the other one is unburnt. The indeterminacy of the deflagration regime is removed by imposing a (constant) value for the fundamental flame speed of the reactive shock. An iterative algorithm is proposed for the solution of the reactive Riemann problem. Then the reactive Riemann problem and the proposed algorithm are investigated from a numerical point of view in the case in which the unburnt state consists of a stoichiometric mixture of hydrogen and air at almost atmospheric condition. In particular, we revisit the problem of 1D plane‐symmetric steady flames in a semi‐infinite domain and we verify that the transition from one combustion regime to another occurs continuously with respect to the fundamental flame speed and the so‐called piston velocity. Finally, we use the ‘all shock’ solution of the reactive Riemann problem to design an approximate (‘all shock’) Riemann solver. 1D and 2D flows at different combustion regimes are computed, which shows that the approximate Riemann solver, and thus the algorithm we use for the solution of the reactive Riemann problem, is robust in both the deflagration and detonation regimes. Copyright © 2009 John Wiley & Sons, Ltd.     

一种基于黎曼解处理大密度比多相流SPH的改进算法

多相流界面存在密度、黏性等物理场间断,直接采用传统光滑粒子水动力学(smoothedparticle hydrodynamics,SPH)方法进行数值模拟,界面附近的压力和速度存在震荡.一套基于黎曼解能够处理大密度比的多相流SPH计算模型被提出,该模型利用黎曼解在处理接触间断问题方面的优势,将黎曼解引入到SPH多相流计算模型中,为了能够准确求解多相流体物理黏性、减小黎曼耗散,对黎曼形式的SPH动量方程进行了改进,又将Adami固壁边界与黎曼单侧问题相结合来施加多相流SPH固壁边界,同时模型中考虑了表面张力对小尺度异相界面的影响,该模型没有添加任何人工黏性、人工耗散和非物理人工处理技术,能够反应多相流真实物理黏性和物理演变状态.采用该模型首先对三种不同粒子间距离散下方形液滴震荡问题进行了数值模拟,验证了该模型在处理异相界面的正确性和模型本身的收敛性;后又通过对Rayleigh--Taylor不稳定、单气泡上浮、双气泡上浮问题进行了模拟计算,结果与文献对比吻合度高,异相界面捕捉清晰,结果表明,本文改进的多相流SPH模型能够稳定、有效的模拟大密度比和黏性比的多相流问题.      

MUSTA schemes for multi‐dimensional hyperbolic systems: analysis and improvements

We develop and analyse an improved version of the multi‐stage (MUSTA) approach to the construction of upwind Godunov‐type fluxes whereby the solution of the Riemann problem, approximate or exact, is not required. The new MUSTA schemes improve upon the original schemes in terms of monotonicity properties, accuracy and stability in multiple space dimensions. We incorporate the MUSTA technology into the framework of finite‐volume weighted essentially nonoscillatory schemes as applied to the Euler equations of compressible gas dynamics. The results demonstrate that our new schemes are good alternatives to current centred methods and to conventional upwind methods as applied to complicated hyperbolic systems for which the solution of the Riemann problem is costly or unknown. Copyright © 2005 John Wiley & Sons, Ltd.     

A fast riemann solver with constant covolume applied to the random choice method

The Riemann problem for the unsteady one-dimensional Euler equations together with the constant-covolume equation of state is solved exactly. The solution is then applied to the random choice method to solve the general initial-boundary value problem for the Euler equations. The iterative procedure to find p*, the pressure between the acoustic waves, involves a single algebraic (non-linear) equation, all other quantities follow directly throughout the x–t plane, except within rarefaction fans where an extra iterative procedure is required. The solution is validated against existing exact results both directly and in conjunction with the random choice method.     

Riemann solution for ideal isentropic magnetogasdynamics

In the present paper, we study the Riemann problem for quasilinear hyperbolic system of partial differential equations governing the one dimensional ideal isentropic magnetogasdynamics with transverse magnetic field. We discuss the properties of rarefaction waves, shocks and contact discontinuities. Differently from single equation methods rooted in the ideal gasdynamics, the new approach is based on the system of two nonlinear equations imposing the equality of total pressure and velocity, assuming as unknowns the two values of densities, on both sides of the contact discontinuity. Newton iterative method is used to obtain densities. The resulting exact solver is implemented with the examples of general applicability of the proposed approach. For comparisons with exact solution we also shown numerical results obtained by the total variation diminishing slope limiter centre scheme. It is shown that both analytical and numerical results demonstrate the broad applicability and robustness of the new Riemann solver.     

THE UNIQUENESS AND EXISTENCE OF SOLUTION OF THE CHABACTERISTIC PROBLEM ON THE GENERALIZED KdV EQUATION

ＴＨＥＵＮＩＱＵＥＮＥＳＳＡＮＤＥＸＩＳＴＥＮＣＥＯＦＳＯＬＵＴＩＯＮＯＦＴＨＥＣＨＡＢＡＣＴＥＲＩＳＴＩＣＰＲＯＢＬＥＭＯＮＴＨＥＧＥＮＥＲＡＬＩＺＥＤＫｄＶＥＱＵＡＴＩＯＮＬｉＷｅｎ－ｓｈｅｎ（李文深）（ＮｏｒｔｈｅａｓｔＦｏｒｅｓｔｒｙＵｎｉｖ...     

含有非圆形双孔的无限平板中应力的解析解研究

无限平板中含有任意形状单个孔的问题可以使用复变函数方法获得其应力解析解.对于无限平板中含有两个圆孔或两个椭圆孔的双连通域问题,也可以利用多种方法进行求解,比如双极坐标法、应力函数法、复变函数法以及施瓦茨交替法等.其中复变函数中的保角变换方法是获得应力解析解的一个重要方法.但目前尚未见到用此方法求解无限板中含有一个正方形孔和一个椭圆孔的问题.当板在无穷远处受有均布载荷和孔边作用垂直均布压力时,利用保角变换方法可以求解板中含有两个特定形状孔的问题.该方法将所讨论的区域映射成象平面里的一个圆环,其中最关键的一步是找出相应的映射函数.基于黎曼映射定理,提出了该映射函数一般形式,并利用最优化方法,找到了该问题的具体映射函数,然后通过孔边应力边界条件建立了求解两个解析函数的基本方程,获得了该问题的应力解析解.运用ANSYS有限单元法与结果进行了对比.研究了孔距、椭圆形孔大小和两孔布置方位对边界切向应力的影响,以及不同载荷下两孔中心线上应力分布规律.      

固定网格下的特征线法求解溃坝问题

溃坝问题是典型的非线性双曲方程的Riemann问题，其数值求解的难点在于对间断面的捕捉以及避免间断面处在数值计算过程中产生数值色散，因而为求解此问题所产生的各种数值计算方法的优劣也体现在这两个方面。本文针对溃坝问题提出一种新的计算方法。该方法基于对偶变量推导的浅水波方程，根据方程的特点，从方程的特征值和黎曼不变量出发，采用高精度的激波捕捉方法计算黎曼不变量的位置随时间的变化，然后映射至不随时间变化的固定网格。根据黎曼不变量的位置，采用保形分段三次Hermite插值将物理量映射至网格节点。计算结果显示，该方法不仅操作简单，计算量小，而且结果准确。     

Riemann–Hilbert method and multi-soliton solutions of the Kundu-nonlinear Schrödinger equation

In this work, we study the Kundu-nonlinear Schrödinger (Kundu-NLS) equation (so-called the extended NLS equation), which can describe the propagation of the waves in dispersive media. A Lax spectral problem is used to construct the Riemann–Hilbert problem, via a matrix transformation. Based on the inverse scattering transformation, the general solutions of the Kundu-NLS equation are calculated. In the reflection-less case, the special matrix Riemann–Hilbert problem is carefully proposed to derive the multi-soliton solutions. Finally, some novel dynamics behaviors of the nonlinear system are theoretically and graphically discussed.

     

Dispersive shock waves in three dimensional Benjamin–Ono equation

Dispersive shock waves (DSWs) in the three dimensional Benjamin–Ono (3DBO) equation are studied with step-like initial condition along a paraboloid front. By using a similarity reduction, the problem of studying DSWs in three space one time (3+1) dimensions reduces to finding DSW solution of a (1+1) dimensional equation. By using a special ansatz, the 3DBO equation exactly reduces to the spherical Benjamin–Ono (sBO) equation. Whitham modulation equations are derived which describes DSW evolution in the sBO equation by using a perturbation method. These equations are written in terms of appropriate Riemann type variables to obtain the sBO-Whitham system. DSW solution which is obtained from the numerical solutions of the Whitham system and the direct numerical solution of the sBO equation are compared. In this comparison, a good agreement is found between these solutions. Also, some physical qualitative results about DSWs in sBO equation are presented. It is concluded that DSW solutions in the reduced sBO equation provide some information about DSW behavior along the paraboloid fronts in the 3DBO equation.