Spectral element-FCT method for scalar hyperbolic conservation laws

A new algorithm based on spectral element discretizations and flux-corrected transport (FCT) ideas is developed for the solution of discontinuous hyperbolic problems. A conservative formulation is proposed, based on cell averaging and reconstruction procedures, that employs a staggered grid of Gauss–Chebyshev and Gauss–Lobatto–Chebyshev discretizations. In addition, high-order time-differencing schemes, a flux limiter and a general spectral filter are employed to improve the quality of the solution. It is demonstrated through model problems of linear advection and examples of one-dimensional shock formation that the proposed algorithm leads to stable, non-oscillatory solutions of high accuracy away from discontinuities. Typically, spectral or spectral element methods perform very poorly in the presence of even weak discontinuities, although they produce only exponentialy small errors for smooth solutions. Spectral element–FCT methods can provide spectral properties (i.e. minimum dispersion and diffusion errors) as well as great flexibility in the discretization, since a variable number of macroelements or collocation points per element can be employed to accommodate both accuracy and geometric requirements.     

The regularized Chapman-Enskog expansion for scalar conservation laws

Rosenau [R] has recently proposed a regularized version of the Chapman-Enskog expansion of hydrodynamics. This regularized expansion resembles the usual Navier-Stokes viscosity terms at low wave numbers, but unlike the latter, it has the advantage of being a bounded macroscopic approximation to the linearized collision operator.In this paper we study the behavior of the Rosenau regularization of the Chapman-Enskog expansion (R-C-E) in the context of scalar conservation laws. We show that this R-C-E model retains the essential properties of the usual viscosity approximation, e.g., existence of travelling waves, monotonicity, upper-Lipschitz continuity, etc., and at the same time, it sharpens the standard viscous shock layers. We prove that the regularized R-C-E approximation converges to the underlying inviscid entropy solution as its mean-free-path 0, and we estimate the convergence rate.     

Higher-order jump conditions for conservation laws

The hyperbolic conservation laws admit discontinuous solutions where the solution variables can have finite jumps in space and time. The jump conditions for conservation laws are expressed in terms of the speed of the discontinuity and the state variables on both sides. An example from the Gas Dynamics is the Rankine–Hugoniot conditions for the shock speed. Here, we provide an expression for the acceleration of the discontinuity in terms of the state variables and their spatial derivatives on both sides. We derive a jump condition for the shock acceleration. Using this general expression, we show how to obtain explicit shock acceleration formulas for nonlinear hyperbolic conservation laws. We start with the Burgers’ equation and check the derived formula with an analytical solution. We next derive formulas for the Shallow Water Equations and the Euler Equations of Gas Dynamics. We will verify our formulas for the Euler Equations using an exact solution for the spherically symmetric blast wave problem. In addition, we discuss the potential use of these formulas for the implementation of shock fitting methods.     

Decay of solutions of hyperbolic systems of conservation laws with a convex extension

Archive for Rational Mechanics and Analysis -     

The method of quasidecoupling for discontinuous solutions to conservation laws

Communicated by C. M. Dafermos     

Viscous limits for piecewise smooth solutions to systems of conservation laws

In this paper we study the zero dissipation problem for a general system of conservation laws with positive viscosity. It is shown that if the solution of the problem with zero viscosity is piecewise smooth with a finite number of noninteracting shocks satisfying the entropy condition, then there exist solutions to the corresponding system with viscosity that converge to the solutions of the system without viscosity away from shock discontinuities at a rate of order as the viscosity coefficient goes to zero. The proof uses a matched asymptotic analysis and an energy estimate related to the stability theory for viscous shock profiles.     

Method of homogeneous solutions in problems with mixed boundary conditions

Kramatorsk Industrial Institute. Donetsk University. Translated from Prikladnaya Mekhanika, Vol. 25, No. 9, pp. 57–61, September, 1989.     

Existence and uniqueness of solutions for some hyperbolic systems of conservation laws

We study the Cauchy problem for systems of conservation laws which belong to the Temple class. The compensated-compactness theory is used to prove existence of solutions. Some uniqueness results are established by means of Holmgren's principle.     

On symmetries, conservation laws and invariant solutions of the foam-drainage equation

This study deals with symmetry group properties and conservation laws of the foam-drainage equation. Firstly, we study the classical Lie symmetries, optimal systems, similarity reductions and similarity solutions of the foam-drainage equation which are obtained through the Lie group method of infinitesimal transformations. Secondly, using the new general theorem on non-local conservation laws and partial Lagrangian approach, local and non-local conservation laws are also studied and, finally, non-classical symmetries are derived.     

Analytical solutions for thermal vibration of nanobeams with elastic boundary conditions

A nonlocal Euler beam model with second-order gradient of stress taken into consideration is used to study the thermal vibration of nanobeams with elastic boundary.An analytical solution is proposed to investigate the free vibration of nonlocal Euler beams subjected to axial thermal stress.The effects of the nonlocal parameter,thermal stress and stiffness of boundary constraint on the vibration behaviors of nanobeams are revealed.The results show that natural frequencies including the thermal stress are lower than those without the thermal stress when temperature rises.The boundary-constrained springs have significant effects on the vibration of nanobeams.In addition,numerical simulations also indicate the importance of small-scale effect on the vibration of nanobeams.     

Supercritical equilibrium solutions of axially moving beams with hybrid boundary conditions

In this paper supercritical equilibria and critical speeds of axially moving beams constrained by sleeves with torsion springs are deduced. Transverse vibration of the beams is governed by a nonlinear integro-partial-differential equation. In the supercritical regime, the corresponding static equilibrium equation for the hybrid boundary conditions is analytically solved for the equilibria and the critical speeds. In the view of the non-trivial equilibrium, comparisons are made among the integro-partial-differential equation, a nonlinear partial-differential equation for transverse vibration, and coupled equations for planar motion under the hybrid boundary conditions.     

Evaluation of discontinuous Galerkin and spectral volume methods for scalar and system conservation laws on unstructured grids

The discontinuous Galerkin (DG) and spectral volume (SV) methods are two recently developed high‐order methods for hyperbolic conservation laws capable of handling unstructured grids. In this paper, their overall performance in terms of efficiency, accuracy and memory requirement is evaluated using a 2D scalar conservation laws and the 2D Euler equations. To measure their accuracy, problems with analytical solutions are used. Both methods are also used to solve problems with strong discontinuities to test their ability in discontinuity capturing. Both the DG and SV methods are capable of achieving their formal order of accuracy while the DG method has a lower error magnitude and takes more memory. They are also similar in efficiency. The SV method appears to have a higher resolution for discontinuities because the data limiting can be done at the sub‐element level. Copyright © 2004 John Wiley & Sons, Ltd.     

Optimal system,group invariant solutions and conservation laws of the CGKP equation

In this paper, the (2 + 1)-dimensional cubic generalized Kadomtsev–Petviashvili (CGKP) equation that is derived from the Maxwell–Bloch equations is investigated. By means of Lie symmetry analysis method, we obtain the Lie point symmetries for the equation and the optimal system of the symmetry algebra. Based on the optimal system, a lot of group invariant solutions are obtained. In addition, explicit conservation laws of the equation are studied.     

Exact solutions of multi-term fractional difusion-wave equations with Robin type boundary conditions

General exact solutions in terms of wavelet expansion are obtained for multiterm time-fractional difusion-wave equations with Robin type boundary conditions. By proposing a new method of integral transform for solving boundary value problems, such fractional partial diferential equations are converted into time-fractional ordinary diferential equations, which are further reduced to algebraic equations by using the Laplace transform. Then, with a wavelet-based exact formula of Laplace inversion, the resulting exact solutions in the Laplace transform domain are reversed to the time-space domain. Three examples of wave-difusion problems are given to validate the proposed analytical method.