HIGH RESOLUTION SCHEMES FOR CONSERVATION LAWS AND CONVECTION-DIFFUSION EQUATIONS WITH VARYING TIME AND SPACE GRIDS

This paper presents a class of high resolution local time step schemes for nonlinear hyperbolic conservation laws and the closely related convection-diffusion equations, by projecting the solution increments of the underlying partial differential equations （PDE） at each local time step. The main advantages are that they are of good consistency, and it is convenient to implement them. The schemes are L^∞ stable, satisfy a cell entropy inequality, and may be extended to the initial boundary value problem of general unsteady PDEs with higher-order spatial derivatives. The high resolution schemes are given by combining the reconstruction technique with a second order TVD Runge-Kutta scheme or a Lax-Wendroff type method, respectively. The schemes are used to solve a linear convection-diffusion equation, the nonlinear inviscid Burgers＇ equation, the one- and two-dimensional compressible Euler equations, and the two-dimensional incompressible Navier-Stokes equations. The numerical results show that the schemes are of higher-order accuracy, and efficient in saving computational cost, especially, for the case of combining the present schemes with the adaptive mesh method [15]. The correct locations of the slow moving or stronger discontinuities are also obtained, although the schemes are slightly nonconservative.

Study on a Cross Diffusion Parabolic System

This paper considers a kind of strongly coupled cross diffusion parabolic system, which can be used as the multi-dimensional Lyurnkis energy transport model in semiconductor science. The global existence and large time behavior are obtained for smooth solution to the initial boundary value problem. When the initial data are a small perturbation of an isothermal stationary solution, the smooth solution of the problem under the insulating boundary condition, converges to that stationary solution exponentially fast as time goes to infinity.

A Space-Time Conservative Method for Hyperbolic Systems with Stiff and Non Stiff Source Terms

In this article we propose a higher-order space-time conservative method for hyperbolic systems with stiff and non stiff source terms as well as relaxation systems.We call the scheme a slope propagation(SP)method.It is an extension of our scheme derived for homogeneous hyperbolic systems[1].In the present inhomogeneous systems the relaxation time may vary from order of one to a very small value.These small values make the relaxation term stronger and highly stiff.In such situations underresolved numerical schemes may produce spurious numerical results.However,our present scheme has the capability to correctly capture the behavior of the physical phenomena with high order accuracy even if the initial layer and the small relaxation time are not numerically resolved.The scheme treats the space and time in a unified manner.The flow variables and their slopes are the basic unknowns in the scheme.The source term is treated by its volumetric integration over the space-time control volume and is a direct part of the overall space-time flux balance.We use two approaches for the slope calculations of the flow variables,the first one results directly from the flux balance over the control volumes,while in the second one we use a finite difference approach.The main features of the scheme are its simplicity,its Jacobian-free and Riemann solver-free recipe,as well as its efficiency and high of order accuracy.In particular we show that the scheme has a discrete analog of the continuous asymptotic limit.We have implemented our scheme for various test models available in the literature such as the Broadwell model,the extended thermodynamics equations,the shallow water equations,traffic flow and the Euler equations with heat transfer.The numerical results validate the accuracy,versatility and robustness of the present scheme.

Modeling of aggregation kernels for fluidized beds using discrete particle model simulations

Aggregation is one of the many important processes in chemical and process engineering. Several researchers have attempted to understand this complex process in fluidized beds using the macro-model of population balance equations （PBEs）. The aggregation kernel is an effective parameter in PBEs, and is defined as the product of the aggregation efficiency and collision frequency functions. Attempts to derive this kernel have taken different approaches, including theoretical, experimental, and empirical techniques. The present paper calculates the aggregation kernel using micro-model computer simulations, i.e., a discrete particle model. We simulate the micro-model without aggregation for various initial conditions, and observe that the collision frequency function is in good agreement with the shear kernel. We then simulate the micro-model with aggregation and calculate the aggregation efficiency rate.