In this paper, the Riemann solutions for scalar conservation laws with discontinuous flux function were constructed. The interactions of elementary waves of the conservation laws were concerned, and the numerical simu...In this paper, the Riemann solutions for scalar conservation laws with discontinuous flux function were constructed. The interactions of elementary waves of the conservation laws were concerned, and the numerical simulations were given.展开更多
The Riemann problem for a two-dimensional 2 x 2 nonstrictly hyperbolic system of nonlinear conservation laws has been considered for constant initial data having discontinuities on three rays with vertex at the origin...The Riemann problem for a two-dimensional 2 x 2 nonstrictly hyperbolic system of nonlinear conservation laws has been considered for constant initial data having discontinuities on three rays with vertex at the origin. The solutions are constructed for some one-J and non-R initial data. One kind of new discontinuity, which is labelled as the delta-shock wave, appears in some solutions. The delta-shock wave is a discontinuity plane that is the suport of a generalized function.展开更多
In this article, we get non-selfsimilar elementary waves of the conservation laws in another kind of view, which is different from the usual self-similar transformation. The solution has different global structure. Th...In this article, we get non-selfsimilar elementary waves of the conservation laws in another kind of view, which is different from the usual self-similar transformation. The solution has different global structure. This article is divided into three parts. The first part is introduction. In the second part, we discuss non-selfsimilar elementary waves and their interactions of a class of twodimensional conservation laws. In this case, we consider the case that the initial discontinuity is parabola with u+ 〉 0, while explicit non-selfsirnilar rarefaction wave can be obtained. In the second part, we consider the solution structure of case u+ 〈 0. The new solution structures are obtained by the interactions between different elementary waves, and will continue to interact with other states. Global solutions would be very different from the situation of one dimension.展开更多
We are concerned with the derivation and analysis of one-dimensional hyperbolic systems of conservation laws modelling fluid flows such as the blood flow through compliant axisyminetric vessels. Early models derived a...We are concerned with the derivation and analysis of one-dimensional hyperbolic systems of conservation laws modelling fluid flows such as the blood flow through compliant axisyminetric vessels. Early models derived are nonconservative and/or nonho- mogeneous with measure source terms, which are endowed with infinitely many Riemann solutions for some Riemann data. In this paper, we derive a one-dimensional hyperbolic system that is conservative and homogeneous. Moreover, there exists a unique global Riemann solution for the Riemann problem for two vessels with arbitrarily large Riemann data, under a natural stability entropy criterion. The Riemann solutions may consist of four waves for some cases. The system can also be written as a 3 × 3 system for which strict hyperbolicity fails and the standing waves can be regarded as the contact discontinuities corresponding to the second family with zero eigenvalue.展开更多
For the two-dimensional(2D)scalar conservation law,when the initial data contain two different constant states and the initial discontinuous curve is a general curve,then complex structures of wave interactions will b...For the two-dimensional(2D)scalar conservation law,when the initial data contain two different constant states and the initial discontinuous curve is a general curve,then complex structures of wave interactions will be generated.In this paper,by proposing and investigating the plus envelope,the minus envelope,and the mixed envelope of 2D non-selfsimilar rarefaction wave surfaces,we obtain and the prove the new structures and classifications of interactions between the 2D non-selfsimilar shock wave and the rarefaction wave.For the cases of the plus envelope and the minus envelope,we get and prove the necessary and sufficient criterion to judge these two envelopes and correspondingly get more general new structures of 2D solutions.展开更多
This paper deals with the Burgers equation which is the most common model used in the nonlinear conservation laws. Here the theoretical aspect of conservation law is discussed by using inviscid Burgers equation. At fi...This paper deals with the Burgers equation which is the most common model used in the nonlinear conservation laws. Here the theoretical aspect of conservation law is discussed by using inviscid Burgers equation. At first, we introduce the general non-linear conservation law as a partial differential equation and its solution procedure by the method of characteristic. Next, we present the weak solution of the problem with entropy condition. Taking into account shock wave and rarefaction wave, the Riemann problem has also been discussed. Finally, the finite volume method is considered to approximate the numerical solution of the inviscid Burgers equation with continuous and discontinuous initial data. An illustration of the problem is provided by some examples. Moreover, the Godunov method provides a good approximation for the problem.展开更多
This paper is concerned with the initial-boundary value problem of a nonlinear conservation law in the half space R+= {x |x > 0} where a>0 , u(x,t) is an unknown function of x ∈ R+ and t>0 , u ± , um ar...This paper is concerned with the initial-boundary value problem of a nonlinear conservation law in the half space R+= {x |x > 0} where a>0 , u(x,t) is an unknown function of x ∈ R+ and t>0 , u ± , um are three given constants satisfying um=u+≠u- or um=u-≠u+ , and the flux function f is a given continuous function with a weak discontinuous point ud. The main purpose of our present manuscript is devoted to studying the structure of the global weak entropy solution for the above initial-boundary value problem under the condition of f '-(ud) > f '+(ud). By the characteristic method and the truncation method, we construct the global weak entropy solution of this initial-boundary value problem, and investigate the interaction of elementary waves with the boundary and the boundary behavior of the weak entropy solution.展开更多
In this paper,it is proved that the weak solution to the Cauchy problem for the scalar viscous conservation law,with nonlinear viscosity,different far field states and periodic perturbations,not only exists globally i...In this paper,it is proved that the weak solution to the Cauchy problem for the scalar viscous conservation law,with nonlinear viscosity,different far field states and periodic perturbations,not only exists globally in time,but also converges towards the viscous shock wave of the corresponding Riemann problem as time goes to infinity.Furthermore,the decay rate is shown.The proof is given by a technical energy method.展开更多
In this paper,we study the shock waves for a mixed-type system from chemotaxis.We are concerned with the jump conditions for the left state which is located in the elliptical region and the right state in the hyperbol...In this paper,we study the shock waves for a mixed-type system from chemotaxis.We are concerned with the jump conditions for the left state which is located in the elliptical region and the right state in the hyperbolic region.Under the generalized entropy conditions,we find that there are different shock wave structures for different parameters.To guarantee the uniqueness of the solutions,we obtain the admissible shock waves which satisfy the generalized entropy condition in both parameters.Finally,we construct the Riemann solutions in some solvable regions.展开更多
We develop an embedded boundary finite difference technique for solving the compressible two-or three-dimensional Euler equations in complex geometries on a Cartesian grid.The method is second order accurate with an e...We develop an embedded boundary finite difference technique for solving the compressible two-or three-dimensional Euler equations in complex geometries on a Cartesian grid.The method is second order accurate with an explicit time step determined by the grid size away from the boundary.Slope limiters are used on the embedded boundary to avoid non-physical oscillations near shock waves.We show computed examples of supersonic flow past a cylinder and compare with results computed on a body fitted grid.Furthermore,we discuss the implementation of the method for thin geometries,and show computed examples of transonic flow past an airfoil.展开更多
A new approach to obtaining space-time conservation schemes is proposed.This approach is very simple and easy to apply in high dimension problems. Numerica1 re-su1ts show that the new schemes obtained by this method n...A new approach to obtaining space-time conservation schemes is proposed.This approach is very simple and easy to apply in high dimension problems. Numerica1 re-su1ts show that the new schemes obtained by this method not only are very efficient andaccurate, but also have very high shock resolution.展开更多
This paper addresses tensile shock physics in thermoviscoelastic (TVE) solids without memory. The mathematical model is derived using conservation and balance laws (CBL) of classical continuum mechanics (CCM), incorpo...This paper addresses tensile shock physics in thermoviscoelastic (TVE) solids without memory. The mathematical model is derived using conservation and balance laws (CBL) of classical continuum mechanics (CCM), incorporating the contravariant second Piola-Kirchhoff stress tensor, the covariant Green’s strain tensor, and its rates up to order n. This mathematical model permits the study of finite deformation and finite strain compressible deformation physics with an ordered rate dissipation mechanism. Constitutive theories are derived using conjugate pairs in entropy inequality and the representation theorem. The resulting mathematical model is both thermodynamically and mathematically consistent and has closure. The solution of the initial value problems (IVPs) describing evolutions is obtained using a variationally consistent space-time coupled finite element method, derived using space-time residual functional in which the local approximations are in hpk higher-order scalar product spaces. This permits accurate description problem physics over the discretization and also permits precise a posteriori computation of the space-time residual functional, an accurate measure of the accuracy of the computed solution. Model problem studies are presented to demonstrate tensile shock formation, propagation, reflection, and interaction. A unique feature of this research is that tensile shocks can only exist in solid matter, as their existence requires a medium to be elastic (presence of strain), which is only possible in a solid medium. In tensile shock physics, a decrease in the density of the medium caused by tensile waves leads to shock formation ahead of the wave. In contrast, in compressive shocks, an increase in density and the corresponding compressive waves result in the formation of compression shocks behind of the wave. Although these are two similar phenomena, they are inherently different in nature. To our knowledge, this work has not been reported in the published literature.展开更多
In this paper, we investigate the elementary wave interactions of the Aw-Rascle model for the generalized Chaplygin gas. We construct the unique solution by the characteristic analysis method and obtain the stability ...In this paper, we investigate the elementary wave interactions of the Aw-Rascle model for the generalized Chaplygin gas. We construct the unique solution by the characteristic analysis method and obtain the stability of the corresponding Riemann solutions under such small perturbations on the initial values. We find that the elementary wave interactions have a much more simple structure for Temple class than general systems of conservation laws. It is important to study the elementary waves interactions of the traffic flow system for the generalized Chaplygin gas not only because of their significance in practical applications in the traffic flow system, but also because of their basic role for the general mathematical theory.展开更多
In this paper, we investigate the elementary wave interactions for the Suliciu relaxation system and construct uniquely the solution by the characteristic analysis method in the phase plane. We find that the elementar...In this paper, we investigate the elementary wave interactions for the Suliciu relaxation system and construct uniquely the solution by the characteristic analysis method in the phase plane. We find that the elementary wave interactions have a much simpler structure for the Temple class than the general systems of conservation laws. It is observed that the Riemann solutions of the Suliciu relaxation system are stable under the small perturbation on the Riemann initial data.展开更多
For a conservation law with convex condition and initial data in L∞(R), it had been commonly believed that the number of discontinuity lines (or shock waves) of the solution is at most countable since Theorem 1 in Ol...For a conservation law with convex condition and initial data in L∞(R), it had been commonly believed that the number of discontinuity lines (or shock waves) of the solution is at most countable since Theorem 1 in Oleinik's seminal paper published in 1956 asserted this fact. In 1977, the author gave an example to show that there is an initial data in C∞(R) ∩ L∞(R) such that the number of shock waves is uncountable. And in 1980, he gave an example to show that there is an initial data in C(R)∩L∞(R) such that the measure of original points of shock waves on the real axis is positive. In this paper, he proves further that the set consisting of initial data in C(R) ∩ L∞(R) with the property: almost all points on the real axis are original points of shock waves, is dense in C(R) ∩ L∞(R). All these results show that Oleinik's assertion on the countability of discontinuity lines is wrong.展开更多
基金Project supported by National Natural Science Foundation of China(Grant No .10271072)
文摘In this paper, the Riemann solutions for scalar conservation laws with discontinuous flux function were constructed. The interactions of elementary waves of the conservation laws were concerned, and the numerical simulations were given.
文摘The Riemann problem for a two-dimensional 2 x 2 nonstrictly hyperbolic system of nonlinear conservation laws has been considered for constant initial data having discontinuities on three rays with vertex at the origin. The solutions are constructed for some one-J and non-R initial data. One kind of new discontinuity, which is labelled as the delta-shock wave, appears in some solutions. The delta-shock wave is a discontinuity plane that is the suport of a generalized function.
基金Sponsored by the National Natural Science Foundation of China (10671116,10871199, and 10001023)Hou Yingdong Fellowship (81004), The China Scholarship Council, Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry, Natural Science Foundation of Guangdong (06027210 and 000804)Natural Science Foundation of Guangdong Education Bureau (200030)
文摘In this article, we get non-selfsimilar elementary waves of the conservation laws in another kind of view, which is different from the usual self-similar transformation. The solution has different global structure. This article is divided into three parts. The first part is introduction. In the second part, we discuss non-selfsimilar elementary waves and their interactions of a class of twodimensional conservation laws. In this case, we consider the case that the initial discontinuity is parabola with u+ 〉 0, while explicit non-selfsirnilar rarefaction wave can be obtained. In the second part, we consider the solution structure of case u+ 〈 0. The new solution structures are obtained by the interactions between different elementary waves, and will continue to interact with other states. Global solutions would be very different from the situation of one dimension.
基金supported in part by the National Science Foundation under Grants DMS-0935967the National Science Foundation under Grants DMS-0807551+2 种基金the National Science Foundation under Grants DMS-0720925the National Science Foundation under Grants DMS-0505473the Natural Science Foundation of China under Grant NSFC-10728101,and the Royal Society-Wolfson Research Merit Award (UK)
文摘We are concerned with the derivation and analysis of one-dimensional hyperbolic systems of conservation laws modelling fluid flows such as the blood flow through compliant axisyminetric vessels. Early models derived are nonconservative and/or nonho- mogeneous with measure source terms, which are endowed with infinitely many Riemann solutions for some Riemann data. In this paper, we derive a one-dimensional hyperbolic system that is conservative and homogeneous. Moreover, there exists a unique global Riemann solution for the Riemann problem for two vessels with arbitrarily large Riemann data, under a natural stability entropy criterion. The Riemann solutions may consist of four waves for some cases. The system can also be written as a 3 × 3 system for which strict hyperbolicity fails and the standing waves can be regarded as the contact discontinuities corresponding to the second family with zero eigenvalue.
基金supported in part by the NSFC(Grant No.11471332)The research of Gao-wei Cao was supported in part by the NSFC(Grant No.11701551).
文摘For the two-dimensional(2D)scalar conservation law,when the initial data contain two different constant states and the initial discontinuous curve is a general curve,then complex structures of wave interactions will be generated.In this paper,by proposing and investigating the plus envelope,the minus envelope,and the mixed envelope of 2D non-selfsimilar rarefaction wave surfaces,we obtain and the prove the new structures and classifications of interactions between the 2D non-selfsimilar shock wave and the rarefaction wave.For the cases of the plus envelope and the minus envelope,we get and prove the necessary and sufficient criterion to judge these two envelopes and correspondingly get more general new structures of 2D solutions.
文摘This paper deals with the Burgers equation which is the most common model used in the nonlinear conservation laws. Here the theoretical aspect of conservation law is discussed by using inviscid Burgers equation. At first, we introduce the general non-linear conservation law as a partial differential equation and its solution procedure by the method of characteristic. Next, we present the weak solution of the problem with entropy condition. Taking into account shock wave and rarefaction wave, the Riemann problem has also been discussed. Finally, the finite volume method is considered to approximate the numerical solution of the inviscid Burgers equation with continuous and discontinuous initial data. An illustration of the problem is provided by some examples. Moreover, the Godunov method provides a good approximation for the problem.
文摘This paper is concerned with the initial-boundary value problem of a nonlinear conservation law in the half space R+= {x |x > 0} where a>0 , u(x,t) is an unknown function of x ∈ R+ and t>0 , u ± , um are three given constants satisfying um=u+≠u- or um=u-≠u+ , and the flux function f is a given continuous function with a weak discontinuous point ud. The main purpose of our present manuscript is devoted to studying the structure of the global weak entropy solution for the above initial-boundary value problem under the condition of f '-(ud) > f '+(ud). By the characteristic method and the truncation method, we construct the global weak entropy solution of this initial-boundary value problem, and investigate the interaction of elementary waves with the boundary and the boundary behavior of the weak entropy solution.
文摘In this paper,it is proved that the weak solution to the Cauchy problem for the scalar viscous conservation law,with nonlinear viscosity,different far field states and periodic perturbations,not only exists globally in time,but also converges towards the viscous shock wave of the corresponding Riemann problem as time goes to infinity.Furthermore,the decay rate is shown.The proof is given by a technical energy method.
基金the National Natural Science Foundation of China(11771442)。
文摘In this paper,we study the shock waves for a mixed-type system from chemotaxis.We are concerned with the jump conditions for the left state which is located in the elliptical region and the right state in the hyperbolic region.Under the generalized entropy conditions,we find that there are different shock wave structures for different parameters.To guarantee the uniqueness of the solutions,we obtain the admissible shock waves which satisfy the generalized entropy condition in both parameters.Finally,we construct the Riemann solutions in some solvable regions.
基金performed under the auspices of the U.S.Department of Energy by University of California Lawrence Livermore National Laboratory under contract No.W7405-Eng-48.
文摘We develop an embedded boundary finite difference technique for solving the compressible two-or three-dimensional Euler equations in complex geometries on a Cartesian grid.The method is second order accurate with an explicit time step determined by the grid size away from the boundary.Slope limiters are used on the embedded boundary to avoid non-physical oscillations near shock waves.We show computed examples of supersonic flow past a cylinder and compare with results computed on a body fitted grid.Furthermore,we discuss the implementation of the method for thin geometries,and show computed examples of transonic flow past an airfoil.
文摘A new approach to obtaining space-time conservation schemes is proposed.This approach is very simple and easy to apply in high dimension problems. Numerica1 re-su1ts show that the new schemes obtained by this method not only are very efficient andaccurate, but also have very high shock resolution.
文摘This paper addresses tensile shock physics in thermoviscoelastic (TVE) solids without memory. The mathematical model is derived using conservation and balance laws (CBL) of classical continuum mechanics (CCM), incorporating the contravariant second Piola-Kirchhoff stress tensor, the covariant Green’s strain tensor, and its rates up to order n. This mathematical model permits the study of finite deformation and finite strain compressible deformation physics with an ordered rate dissipation mechanism. Constitutive theories are derived using conjugate pairs in entropy inequality and the representation theorem. The resulting mathematical model is both thermodynamically and mathematically consistent and has closure. The solution of the initial value problems (IVPs) describing evolutions is obtained using a variationally consistent space-time coupled finite element method, derived using space-time residual functional in which the local approximations are in hpk higher-order scalar product spaces. This permits accurate description problem physics over the discretization and also permits precise a posteriori computation of the space-time residual functional, an accurate measure of the accuracy of the computed solution. Model problem studies are presented to demonstrate tensile shock formation, propagation, reflection, and interaction. A unique feature of this research is that tensile shocks can only exist in solid matter, as their existence requires a medium to be elastic (presence of strain), which is only possible in a solid medium. In tensile shock physics, a decrease in the density of the medium caused by tensile waves leads to shock formation ahead of the wave. In contrast, in compressive shocks, an increase in density and the corresponding compressive waves result in the formation of compression shocks behind of the wave. Although these are two similar phenomena, they are inherently different in nature. To our knowledge, this work has not been reported in the published literature.
文摘In this paper, we investigate the elementary wave interactions of the Aw-Rascle model for the generalized Chaplygin gas. We construct the unique solution by the characteristic analysis method and obtain the stability of the corresponding Riemann solutions under such small perturbations on the initial values. We find that the elementary wave interactions have a much more simple structure for Temple class than general systems of conservation laws. It is important to study the elementary waves interactions of the traffic flow system for the generalized Chaplygin gas not only because of their significance in practical applications in the traffic flow system, but also because of their basic role for the general mathematical theory.
文摘In this paper, we investigate the elementary wave interactions for the Suliciu relaxation system and construct uniquely the solution by the characteristic analysis method in the phase plane. We find that the elementary wave interactions have a much simpler structure for the Temple class than the general systems of conservation laws. It is observed that the Riemann solutions of the Suliciu relaxation system are stable under the small perturbation on the Riemann initial data.
基金supported by National Natural Science Foundation of China (Grant No.10771206)
文摘For a conservation law with convex condition and initial data in L∞(R), it had been commonly believed that the number of discontinuity lines (or shock waves) of the solution is at most countable since Theorem 1 in Oleinik's seminal paper published in 1956 asserted this fact. In 1977, the author gave an example to show that there is an initial data in C∞(R) ∩ L∞(R) such that the number of shock waves is uncountable. And in 1980, he gave an example to show that there is an initial data in C(R)∩L∞(R) such that the measure of original points of shock waves on the real axis is positive. In this paper, he proves further that the set consisting of initial data in C(R) ∩ L∞(R) with the property: almost all points on the real axis are original points of shock waves, is dense in C(R) ∩ L∞(R). All these results show that Oleinik's assertion on the countability of discontinuity lines is wrong.