Some properties of the equations governing a two-dimensional quasi-gasdynamic model of traffic flows

A recently proposed two-dimensional quasi-gasdynamic model of traffic flows is considered. Its Petrovskii parabolicity is analyzed, and the stability of small perturbations against a constant background is investigated. In a nonlinear setting, an energy equality is derived and an energy estimate of the solution is obtained.     

Numerical analysis of a quasi-gasdynamic algorithm as applied to the Euler equations

The computational properties of a quasi-gasdynamic algorithm are analyzed as applied to strong-shock Riemann problems and acoustic disturbance propagation. It is shown that a unified treatment can be given to these problems in the framework of the algorithm. Stability conditions are numerically obtained, and the error and computational complexity of the difference scheme are estimated.     

Finite element simulation of compressible particle-laden gas flows

A macroscopic two-fluid model of compressible particle-laden gas flows is considered. The governing equations are discretized by a high-resolution finite element method based on algebraic flux correction. A multidimensional limiter of TVD type is employed to constrain the local characteristic variables for the continuous gas phase and conservative fluxes for a suspension of solid particles. Special emphasis is laid on the efficient computation of steady state solutions at arbitrary Mach numbers. To avoid stability restrictions and convergence problems, the characteristic boundary conditions are imposed weakly and treated in a fully implicit manner. A two-way coupling via the interphase drag force is implemented using operator splitting. The Douglas-Rachford scheme is found to provide a robust treatment of the interphase exchange terms within the framework of a fractional-step solution strategy. Two-dimensional simulation results are presented for a moving shock wave and for a steady nozzle flow.     

Stability criterion for small perturbations for a quasi-gasdynamic system of equations

The stability of small perturbations against a constant background is studied for a system of quasi-gasdynamic equations in an arbitrary number of space variables. It is established that, for a fixed adiabatic exponent γ, the stability is determined only by the background Mach number, and a necessary and sufficient condition for stability at any Mach number is $\gamma \leqslant \bar \gamma $ , where $\bar \gamma \approx 6.2479$ . The proof is based on a direct analysis of the corresponding complex characteristic numbers depending on several parameters. The multidimensional case is successfully reduced to the one-dimensional one. Then, the generalized Routh-Hurwitz criterion is applied in conjunction with analytical calculations based on Mathematica.     

From gas dynamics to pressureless gas dynamics

This paper is devoted to the convergence of solutions of the compressible Euler equations towards solutions of the pressureless gas dynamics system, when the pressure tends to 0. The goal is to prove accurate uniform bounds for particular solutions of the Euler equations.

     

A fast BIE iteration method for an arbitrary body in a flow of incompressible inviscid fluid

The paper is concerned with the new iteration algorithm to solve boundary integral equations arising in boundary value problems of mathematical physics. The stability of the algorithm is demonstrated on the problem of a flow around bodies placed in the incompressible inviscid fluid. With a discrete numerical treatment, we approximate the exact matrix by a certain Töeplitz one and then apply a fast algorithm for this matrix, on each iteration step. We illustrate the convergence of this iteration scheme by a number of numerical examples, both for hard and soft boundary conditions. It appears that the method is highly efficient for hard boundaries, being much less efficient for soft boundaries.     

Numerical study of the interaction between shocks and rarefaction waves in an ideal gas

The interaction between shock waves and rarefaction waves is numerically studied using the one-dimensional Euler equations for an ideal gas. A specific form of solutions, which are called contact regions, is detected. They represent extended zones with continuously varying density and temperature at constant pressure and velocity. It is shown that, at long times, the solutions to the interaction problem tend to those to the Riemann problems with the contact discontinuity replaced by a contact region.     

The inviscid limit for an inflow problem of compressible viscous gas in presence of both shocks and boundary layers

In this paper, we study the inviscid limit problem for the Navier-Stokes equations of one-dimensional compressible viscous gas on half plane. We prove that if the solution of the inviscid Euler system on half plane is piecewise smooth with a single shock satisfying the entropy condition, then there exist solutions to Navier-Stokes equations which converge to the inviscid solution away from the shock discontinuity and the boundary at an optimal rate of ε¹ as the viscosity ε tends to zero.     

Error analysis of the L2 least‐squares finite element method for incompressible inviscid rotational flows

In this article we analyze the L² least‐squares finite element approximations to the incompressible inviscid rotational flow problem, which is recast into the velocity‐vorticity‐pressure formulation. The least‐squares functional is defined in terms of the sum of the squared L² norms of the residual equations over a suitable product function space. We first derive a coercivity type a priori estimate for the first‐order system problem that will play the crucial role in the error analysis. We then show that the method exhibits an optimal rate of convergence in the H¹ norm for velocity and pressure and a suboptimal rate of convergence in the L² norm for vorticity. A numerical example in two dimensions is presented, which confirms the theoretical error estimates. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

Network simplex algorithm for the general equal flow problem

In this paper the general equal flow problem is considered. This is a minimum cost network flow problem with additional side constraints requiring the flow of arcs in some given sets of arcs to take on the same value. This model can be applied to approach water resource system management problems or multiperiod logistic problems in general involving policy restrictions which require some arcs to carry the same amount of flow through the given study period. Although the bases of the general equal flow problem are no longer spanning trees, it is possible to recognize a similar structure that allows us to take advantage of the practical computational capabilities of network models. After characterizing the bases of the problem as good (r+1)-forests, a simplex primal algorithm is developed that exploits the network structure of the problem and requires only slight modifications of the well-known network simplex algorithm.     

Direct numerical simulation of the turbulent flow in an elliptical pipe

The turbulent flow in a pipe with an elliptical cross section is directly simulated at Re = 4000 (where the Reynolds number Re is calculated in terms of the mean velocity and the hydraulic diameter). The incompressible Navier-Stokes equations are solved in curvilinear orthogonal coordinates by using a central-difference approximation in space and a third-order accurate semi-implicit Runge-Kutta method for time integration. The discrete equations inherit some properties of the original differential equations, in particular, the neutrality of the convective terms and of the pressure gradient in the kinetic energy production. The distributions of the mean and fluctuation characteristics of the turbulent motion over the pipe’s cross section are computed.     

A large deviations principle for stochastic flows of viscous fluids

We study the well-posedness of a stochastic differential equation on the two dimensional torus ${\mathbb{T}}^2$, driven by an infinite dimensional Wiener process with drift in the Sobolev space $L^2(0,T;H^1\left({\mathbb{T}}^2\right))$. The solution corresponds to a stochastic Lagrangian flow in the sense of DiPerna Lions. By taking into account that the motion of a viscous incompressible fluid on the torus can be described through a suitable stochastic differential equation of the previous type, we study the inviscid limit. By establishing a large deviations principle, we show that, as the viscosity goes to zero, the Lagrangian stochastic Navier–Stokes flow approaches the Euler deterministic Lagrangian flow with an exponential rate function.     

The linear stability of inviscid shear flow of a vibrationally excited diatomic gas

The problem of the linear stability of plane-parallel shear flows of a vibrationally excited compressible diatomic gas is investigated using a two-temperature gas dynamics model. The necessary and sufficient conditions for stability of the flows considered are obtained using the energy integrals of the corresponding linearized system for the perturbations. It is proved that thermal relaxation produces an additional dissipation factor, which enhances the flow stability. A region of eigenvalues of unstable perturbations is distinguished in the upper complex half-plane. Numerical calculations of the eigenvalues and eigenfunctions of the unstable inviscid modes are carried out. The dependence on the Mach number of the carrier stream, the vibrational relaxation time τ and the degree of non-equilibrium of the vibrational mode is analysed. The most unstable modes with maximum growth rate are obtained. It is shown that in the limit there is a continuous transition to well-known results for an ideal fluid as the Mach number and τ approach zero and for an ideal gas when τ → 0.     

A capacity-rounding algorithm for the minimum-cost circulation problem: A dual framework of the Tardos algorithm

Recently, É. Tardos gave a strongly polynomial algorithm for the minimum-cost circulation problem and solved the open problem posed in 1972 by J. Edmonds and R.M. Karp. Her algorithm runs in O(m ² T(m, n) logm) time, wherem is the number of arcs,n is the number of vertices, andT(m, n) is the time required for solving a maximum flow problem in a network withm arcs andn vertices. In the present paper, taking an approach that is a dual of Tardos's, we also give a strongly polynomial algorithm for the minimum-cost circulation problem. Our algorithm runs in O(m ² S(m, n) logm) time and reduces the computational complexity, whereS(m, n) is the time required for solving a shortest path problem with a fixed origin in a network withm arcs,n vertices, and a nonnegative arc length function. The complexity is the same as that of Orlin's algorithm, recently developed by efficiently implementing the Edmonds-Karp scaling algorithm.     

Note on the uniqueness of subsonic Euler flows in an axisymmetric nozzle

In this paper, we consider the uniqueness of globally subsonic compressible flows through an infinitely long axisymmetric nozzle. The flow is governed by the steady Euler equations and satisfies no-flow boundary conditions on the nozzle walls. We will show that for given mass flux and Bernoulli’s function in the upstream, the subsonic flow is unique in the class of all axisymmetric solutions, which possess the asymptotic behaviors at the far fields. This result extends the uniqueness of solutions in the previous paper Du and Duan (2011) [1].