A perturbational h4 exponential finite-difference scheme for the convective diffusion equation

A perturbational h4 compact exponential finite difference scheme with diagonally dominant coefficient matrix and upwind effect is developed for the convective diffusion equation. Perturbations of second order are exerted on the convective coefficients and source term of an h2 exponential finite difference scheme proposed in this paper based on a transformation to eliminate the upwind effect of the convective diffusion equation. Four numerical examples including one- to three-dimensional model equations of fluid flow and a problem of natural convective heat transfer are given to illustrate the excellent behavior of the present exponential schemes, the h4 accuracy of the perturbational scheme is verified using double precision arithmetic.

Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

A kind of extended Korteweg-de Vries equation and solitary wave solutions for interfacial waves in a two-fluid system

This paper considers interfacial waves propagating along the interface between a two-dimensional two-fluid with a flat bottom and a rigid upper boundary. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. It just focuses on the weakly non-linear small amplitude waves by introducing two small independent parameters: the nonlinearity ratio epsilon, represented by the ratio of amplitude to depth, and the dispersion ratio mu, represented by the square of the ratio of depth to wave length, which quantify the relative importance of nonlinearity and dispersion. It derives an extended KdV equation of the interfacial waves using the method adopted by Dullin et al in the study of the surface waves when considering the order up to O(mu(2)). As expected, the equation derived from the present work includes, as special cases, those obtained by Dullin et al for surface waves when the surface tension is neglected. The equation derived using an alternative method here is the same as the equation presented by Choi and Camassa. Also it solves the equation by borrowing the method presented by Marchant used for surface waves, and obtains its asymptotic solitary wave solutions when the weakly nonlinear and weakly dispersive terms are balanced in the extended KdV equation.

Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.

对流扩散方程的变步长摄动有限差分格式

摄动有限差分(PFD)方法是构造高精度差分格式的一种新方法。变步长摄动有限差分方法是等步长摄动有限差分方法的发展和推广。对需要局部加密网格的计算问题，变步长PFD格式不需要对自变量进行数学变换，且和等步长PFD格式一样，具有如下的共同特点：从变步长一阶迎风格式出发，通过把非微商项(对流系数和源项)作变步长摄动展开，展开幂级数系数通过消去摄动格式修正微分方程的截断误差项求出，由此获得高精度变步长PFD格式。该格式在一、二和三维情况下分别仅使用三、五和七个基点，且具有迎风性。文中利用变步长PFD格式对对流扩散反应模型方程，变系数方程及Burgers方程等进行了数值模拟，并与一阶迎风和二阶中心格式及其问题的精确解作了比较。数值试验表明，与一阶迎风和二阶中心格式相比，变步长PFD格式具有精度高，稳定性与收敛性好的特点。变步长PFD格式与等步长PFD格式相比，变步长PFD解在薄边界层型区域的分辨率得到了明显的提高。

Direct numerical test of the B-G-K model equation by the DSMC method

In the present paper the rarefied gas how caused by the sudden change of the wall temperature and the Rayleigh problem are simulated by the DSMC method which has been validated by experiments both in global flour field and velocity distribution function level. The comparison of the simulated results with the accurate numerical solutions of the B-G-K model equation shows that near equilibrium the BG-K equation with corrected collision frequency can give accurate result but as farther away from equilibrium the B-G-K equation is not accurate. This is for the first time that the error caused by the B-G-K model equation has been revealed.

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

A lattice Boltzmann method for KDV equation

We prepose a 5-bit lattice Boltzmann model for KdV equation. Using Chapman-Enskog expansion and multiscale technique, we obtained high order moments of equilibrium distribution function, and the 3rd dispersion coefficient and 4th order viscosity. The parameters of this scheme can be determined by analysing the energy dissipation.

Modified Simplified Rate Equation Model of Flowing Chemical Oxygen-Iodine Laser and its Application

A modified simplified rate equation (RE) model of flowing chemical oxygen-iodine laser (COIL), which is adapted to both the condition of homogeneous broadening and inhomogeneous broadening being of importance and the condition of inhomogeneous broadening being predominant, is presented for performance analyses of a COIL. By using the Voigt profile function and the gain-equal-loss approximation, a gain expression has been deduced from the rate equations of upper and lower level laser species. This gain expression is adapted to the conditions of very low gas pressure up to quite high pressure and can deal with the condition of lasing frequency being not equal to the central one of spectral profile. The expressions of output power and extraction efficiency in a flowing COIL can be obtained by solving the coupling equations of the deduced gain expression and the energy equation which expresses the complete transformation of the energy stored in singlet delta state oxygen into laser energy. By using these expressions, the RotoCOIL experiment is simulated, and obtained results agree well with experiment data. Effects of various adjustable parameters on the performances of COIL are also presented.

High-Resolution Finite Compact Difference Schemes for Hyperbolic Conservation Laws

A finite compact (FC) difference scheme requiring only bi-diagonal matrix inversion is proposed by using the known high-resolution flux. Introducing TVD or ENO limiters in the numerical flux, several high-resolution FC-schemes of hyperbolic conservation law are developed, including the FC-TVD, third-order FC-ENO and fifth-order FC-ENO schemes. Boundary conditions formulated need only one unknown variable for third-order FC-ENO scheme and two unknown variables for fifth-order FC-ENO scheme. Numerical test results of the proposed FC-scheme were compared with traditional TVD, ENO and WENO schemes to demonstrate its high-order accuracy and high-resolution.