High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

改进GAs算法在大规模资源分配问题中的应用

采用改进ＧＡｓ算法建立了求解大规模规划的资源分配模型．针对大规模资源分配问题的具体特点，设计了合适的ＧＡｓ算子，并以实例验证了算法的合理性及有效性。

字符方程的拟定

本文在文[1]所提出的α指数技巧的基础上,进一步提出了线条叠切法、切叠法、叠切叠法,应用这些方法建立了字符状图样的单一方程,并说明所有常用字符均可列出其近似的单一方程式。最后指出了叠切叠法在一般工程问题中的应用。

Boundary knot method for 2D and 3D Helmholtz and convection-diffusion problems under complicated geometry

The boundary knot method (BKM) of very recent origin is an inherently meshless, integration-free, boundary-type, radial basis function collocation technique for the numerical discretization of general partial differential equation systems. Unlike the method of fundamental solutions, the use of non-singular general solution in the BKM avoids the unnecessary requirement of constructing a controversial artificial boundary outside the physical domain. The purpose of this paper is to extend the BKM to solve 2D Helmholtz and convection-diffusion problems under rather complicated irregular geometry. The method is also first applied to 3D problems. Numerical experiments validate that the BKM can produce highly accurate solutions using a relatively small number of knots. For inhomogeneous cases, some inner knots are found necessary to guarantee accuracy and stability. The stability and convergence of the BKM are numerically illustrated and the completeness issue is also discussed.

The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress,intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

DISTRIBUTION FUNCTION OF A DRIFTING ELECTRON-GAS FOR GENERAL ENERGY-BAND STRUCTURES

The usual application of the Lei-Ting balance equation method for treating electron transport problems makes use of a Fermi distribution function for the electron motion relative to the center of mass. It is pointed out that this presumes the existence of a moving frame of reference that is dynamically equivalent to the rest frame of reference, and this is only true for electrons with a constant effective mass. The method is thus inapplicable to problems where electrons governed by a general energy-band dispersion E(k) are important (such as in miniband conduction). It is demonstrated that this difficulty can be overcome by introducing a distribution function for a drifting electron gas by maximizing the entropy subject to a prescribed average drift velocity. The distribution function reduces directly to the usual Fermi distribution for electron motion relative to the center of mass in the special case of E(k)=($) over bar h(2)\k\(2)/2m*. This maximum entropy treatment of a drifting electron gas provides a physically more direct as well as a more general basis for the application of the balance equation method.