Equilibrium schemes for scalar conservation laws with stiff sources

We consider a simple model case of stiff source terms in hyperbolic conservation laws, namely, the case of scalar conservation laws with a zeroth order source with low regularity. It is well known that a direct treatment of the source term by finite volume schemes gives unsatisfactory results for both the reduced CFL condition and refined meshes required because of the lack of accuracy on equilibrium states. The source term should be taken into account in the upwinding and discretized at the nodes of the grid. In order to solve numerically the problem, we introduce a so-called equilibrium schemes with the properties that (i) the maximum principle holds true; (ii) discrete entropy inequalities are satisfied; (iii) steady state solutions of the problem are maintained. One of the difficulties in studying the convergence is that there are no estimates for this problem. We therefore introduce a kinetic interpretation of upwinding taking into account the source terms. Based on the kinetic formulation we give a new convergence proof that only uses property (ii) in order to ensure desired compactness framework for a family of approximate solutions and that relies on minimal assumptions. The computational efficiency of our equilibrium schemes is demonstrated by numerical tests that show that, in comparison with an usual upwind scheme, the corresponding equilibrium version is far more accurate. Furthermore, numerical computations show that equilibrium schemes enable us to treat efficiently the sources with singularities and oscillating coefficients.

     

Second class particles as microscopic characteristics in totally asymmetric nearest-neighbor -exclusion processes

We prove laws of large numbers for a second class particle in one-dimensional totally asymmetric -exclusion processes, under hydrodynamic Euler scaling. The assumption required is that initially the ambient particle configuration converges to a limiting profile. The macroscopic trajectories of second class particles are characteristics and shocks of the conservation law of the particle density. The proof uses a variational representation of a second class particle, to overcome the problem of lack of information about invariant distributions. But we cannot rule out the possibility that the flux function of the conservation law may be neither differentiable nor strictly concave. To give a complete picture we discuss the construction, uniqueness, and other properties of the weak solution that the particle density obeys.

     

非齐次守恒律方程分片光滑解的粘性方法

1　引　　言考虑非齐次守恒律方程ut+f(u) x =g(u) ,　　 -∞ 0 ,(1 .1 )u(x,0 ) =u0 (x) ,　　 -∞ 0 , (1 .5)g∈ C3且 g是 Lipschitz连续的 ,Lipschitz系数为 L . (1 .6 )对于一般守恒律齐次方程 ,粘性解逼近熵解的收敛阶为 O(ε ) [1 ] .在 f严格凸的条件下 ,其收敛速度可以提高到 O(ε|lnε|+ε) [2 ] ,[3] .本文考虑具有条件 (1 .5) (1 .6 )的非齐次方程(1 .1 ) ,在较广泛的一类初值条件下…     

声子气的状态方程和声子气运动的守恒方程

根据爱因斯坦的狭义相对论中质能的等效关系，把固体(本文指非导体)晶格(原子)的质量分为晶格(原子)的静质量和晶格热振动能量的等效质量两个部分，后者就是固体中声子气的等效质量.晶格(原子)热振动的能量则分为晶格(原子)静质量具有的热能以及声子气质量具有的热能.基于固体的状态方程，导得了晶格静质量热振动的状态方程和声子气的状态方程.声子气在固体介质中的宏观运动就是热量在固体中的传递过程.建立了声子气运动的守恒方程组，分析表明，忽略惯性力时声子气的动量守恒方程就退化为傅里叶导热定律，阐明了傅里叶导热定律的物理本质是声子气驱动力与阻力的平衡方程.当热流密度很大惯性力不能忽略时，傅里叶导热定律不再适用. 关键词： 非傅里叶导热 声子气 声子气质量 状态方程 守恒方程     

光的折射定律(Snell定律)演示实验装置

基于光子动量在不同媒质分界面上沿分界面方向守恒原理,制作了光的折射定律演示装置.使用该装置可演示入射角连续变化,折射角也自动遵循折射定律相应变化.本文介绍了该演示装置的原理、构造及使用方法.     

The Tunneling Radiation Characteristics of Kerr-Newman de Sitter Black Hole

The strictly thermal spectrum in dragging coordinate system and the tunneling radiation characteristics of stationary axisymmetry Kerr-Newman de Sitter black hole is studied. The result shows that the tunneling rates at the event and cosmological horizon are related to the change of Bekenstein-Hawking entropy and that the factual radiation spectrum is not strictly pure thermal. Thus an exact correction to the Hawking thermal spectrum is present.     

NONLOCAL CROWD DYNAMICS MODELS FOR SEVERAL POPULATIONS

This paper develops the basic analytical theory related to some recently introduced crowd dynamics models.Where well posedness was known only locally in time,it is here extended to all of R+.The result...     

Local time-stepping procedures for the space-time conservation element and solution element method

A local time-stepping procedure for the space-time conservation element and solution element (CESE) method has been developed. This new procedure allows for variation of time-step size in both space and time, and can also be extended to become multi-dimensional solvers with structured/unstructured spatial grids. Moreover, it differs substantially in concept and methodology from the existing approaches. By taking full advantage of key concepts of the CESE method, in a simple and efficient manner it can enforce flux conservation across an interface separating grid zones of different time-step sizes. In particular, no correction pass is needed. Numerical experiments show that, for a variety of flow problems involving moving shock and flame discontinuities, accurate and robust numerical simulations can be achieved even with a reduction in time-step size on the order of 10 or higher for grids across a single interface.     

A conservative unstructured scheme for rapidly varied flows

This article introduces a new semi‐implicit, staggered finite volume scheme on unstructured meshes for the modelling of rapidly varied shallow water flows. Rapidly varied flows occur in the inundation of dry land during flooding situations. They typically involve bores and hydraulic jumps after obstacles such as road banks. Near such sudden flow transitions, the grid resolution is often low compared with the gradients of the bathymetry. Locally the hydrostatic pressure assumption may become invalid. In these situations, it is crucial to apply the correct conservation properties to obtain accurate results. An important feature of this scheme is therefore its ability to conserve momentum locally or, by choice, preserve constant energy head along a streamline. This is achieved using a special interpolation method and control volumes for momentum. The efficiency of inundation calculations with locally very high velocities, and in the case of unstructured meshes locally very small grid distances, is severely hampered by the Courant condition. This article provides a solution in the form of a locally implicit time integration for the advective terms that allows for an explicit calculation in most of the domain, while maintaining unconditional stability by implicit calculations only where necessary. The complex geometry of flooded urban areas asks for the flexibility of unstructured meshes. The efficient calculation of the pressure gradient in this, and other semi‐implicit staggered schemes, requires, however, an orthogonality condition to be put on the grid. In this article a simple method is introduced to generate unstructured hybrid meshes that fulfil this requirement. Copyright © 2008 John Wiley & Sons, Ltd.     

A volume-amending method to improve mass conservation of level set approach for incompressible two-phase flows

A volume-amending method is developed both to keep the level set function as an algebraic distance function and to preserve the bubble mass in a level set approach for incompressible two-phase flows with the significantly deformed free interface. After the traditional reinitialization procedure, a vol-ume-amending method is added for correcting the position of the interface according to mass loss/gain error until the mass error falls in the allowable range designated in advance. The level set approach with this volume-amending method incorporated has been validated by three test cases: the motion of a single axisymmetrical bubble or drop in liquid, the motion of a two-dimensional water drop falling through the air into a water pool, and the interactional motion of two buoyancy-driven three- dimensional deformable bubbles. The computational results with this volume-amending method in-corporated are in good agreement with the reported experimental data and the mass is well preserved in all cases.