Navier-Stokes方程二阶速度滑移边界条件的检验

对微尺度气体流动，Navier-Stokes方程和一阶速度滑移边界条件的结果与实验数据相比，在滑移区相互符合，在过渡领域则显著偏离，为改善Navier-Stokes方程在过渡领域的表现，有些研究者尝试引入二阶速度滑移边界条件，如Cercignani模型，Deissler模型和Beskok-Karniadakis模型．以微槽道气体流动为例，将Navier-Stokes方程在不同的二阶速度滑移模型下的结果与动理论的直接模拟Monte Carlo（DSMC）方法和信息保存（IP）方法以及实验数据进行比较．在所考察的3种具有代表性的二阶速度滑移模型中，Cercignani模型表现最好，其所给出的质量流率在Knudsen数为0．4时仍与DSMC和IP结果相符；然而，细致比较表明，Cercignani模型给出的物面滑移速度及其附近的速度分布在滑流区和过渡领域的分界处（Kn=0．1）已明显偏离DSMC和IP的结果．

一个基于非均匀网格上的二阶基本无振荡差分格式

在基本无振荡格式的构造中，将通常的对流通量f的逼近方式推广到对通量导数的逼近，这一构造方法可以有效地应用到非均匀或非结构网格。直接基于非均匀网格上，构造了一个二阶的基本无振荡(ENO)差分格式。该格式具有形式简单，对网格的划分灵活，与传统格式相比不增加计算量等优点，儿个数值算例证明了格式的有效性。

Simulation of Coal Combustion by AUSM Turbulence-Chemistry Char Combustion Model and a Full Two-Fluid Model

An algebraic unified second-order moment (AUSM) turbulence-chemistry model of char combustion is introduced in this paper, to calculate the effect of particle temperature fluctuation on char combustion. The AUSM model is used to simulate gas-particle flows, in coal combustion in a pulverized coal combustor, together with a full two-fluid model for reacting gas-particle flows and coal combustion, including the sub-models as the k-epsilon-k(p) two-phase turbulence niodel, the EBU-Arrhenius volatile and CO combustion model, and the six-flux radiation model. A new method for calculating particle mass flow rate is also used in this model to correct particle outflow rate and mass flow rate for inside sections, which can obey the principle of mass conservation for the particle phase and can also speed up the iterating convergence of the computation procedure effectively. The simulation results indicate that, the AUSM char combustion model is more preferable to the old char combustion model, since the later totally eliminate the influence of particle temperature fluctuation on char combustion rate.

Studies of the Effect of a Coal Concentrator on NO Formation in Swirling Coal Combustion

To develop low-pollution burners, the effect of a coal concentrator on NO formation in swirling coal combustion is studied using both numerical simulation and experiments. The isothermal gas-particle two-phase velocities and particle concentration in a cold model of swirl burners with and without coal concentrators were measured using the phase Doppler particle anemometer (PDPA). A full two-fluid model of reacting gas-particle flows and coal combustion with an algebraic unified second-order moment (AUSM) turbulence-chemistry model for the turbulent reaction rate of NO formation are used to simulate swirling coal combustion and NO formation with different coal concentrators. The results give the turbulent kinetic energy, particle concentration, temperature and NO concentration in cases of with and without coal concentrators. The predicted results for cold two-phase flows are in good agreement with the PDPA measurement results, showing that the coal concentrator increases the turbulence and particle concentration in the recirculation zone. The combustion modeling results indicate that although the coal concentrator increases the turbulence and combustion temperature, but still can remarkably reduce the NO formation due to creating high coal concentration in the recirculation zone.

An Improved Ce/Se Scheme And Its Application To Detonation Propagation

An improved two-dimensional space-time conservation element and solution element ( CE/ SE) method with second-order accuracy is proposed, examined and extended to simulate the detonation propagations using detailed chemical reaction models. The numerical results of planar and cellular detonation are compared with corresponding results by the Chapman-Jouguet theory and experiments, and prove that the method is a new reliable way for numerical simulations of detonation propagation.

Quasi-Dammann Grating With Proportional Intensity Array Spots

A quasi-Dammann grating is proposed to generate array spots with proportional-intensity orders in the far field. To describe the performance of the grating, the uniformities of the array spots are redefined. A two-dimensional even-sampling encode scheme is adopted to design the quasi-Dammann grating. Numerical solutions of the binary-phase quasi-Dammann grating with proportional-intensity orders are given. The experimental results with a third-order quasi-Dammann grating, which has an intensity proportion of 3:2:1 from zero order to second order, are presented. (C) 2008 Optical Society of America

Coupling Of Evaporation And Thermocapillary Convection In A Liquid Layer With Mass And Heat Exchanging Interface

We propose and analyse a new model of thermocapillary convection with evaporation in a cavity subjected to horizontal temperature gradient, rather than the previously studied model without evaporation. The pure liquid layer with a top free surface in contact with its own vapour is considered in microgravity condition. The computing programme developed for simulating this model integrates the two-dimensional, time-dependent Navier-Stokes equations and energy equation by a second-order accurate projection method. We focus on the coupling of evaporation and thermocapillary convection by investigating the influence of evaporation Biot number and Marangoni number on the interfacial mass and heat transfer. Three different regimes of the coupling mechanisms are found and explained from our numerical results.

Hydrodynamic interaction between two vertical cylinders in water waves

The hydrodynamic interaction between two vertical cylinders in water waves is investigated based on the linearized potential flow theory. One of the two cylinders is fixed at the bottom while the other is articulated at the bottom and oscillates with small amplitudes in the direction of the incident wave. Both the diffracted wave and the radiation wave are studied in the present paper. A simple analytical expression for the velocity potential on the surface of each cylinder is obtained by means of Graf's addition theorem. The wave-excited forces and moments on the cylinders, the added masses and the radiation damping coefficients of the oscillating cylinder are all expressed explicitly in series form. The coefficients of the series are determined by solving algebraic equations. Several numerical examples are given to illustrate the effects of various parameters, such as the separation distance, the relative size of the cylinders, and the incident angle, on the first-order and steady second-order forces, the added masses and radiation-damping coefficients as well as the response of the oscillating cylinder.

Influence of thermal residual stresses on the composite macroscopic behavior

The influence of the thermal residual stress on the deformation behavior of a composite has been analyzed with a new micromechanical method. The method is based on secant moduli approximation and a new homogenized effective stress to characterize the plastic state of the matrix. It is found that the generated thermal residual stresses after cooling and their influence on the subsequent deformation behavior depends significantly on the aspect ratio of the inclusions. With prolate inclusions, the presence of thermal residual stresses generate a higher compressive hardening curves of the composite, but it is reversed with oblate inclusions. For particle reinforced composite, thermal residual stresses induce a tensile hardening curve higher than the compressive one and this is in agreement with experimental observations. (C) 1998 Elsevier Science Ltd.

A micromechanical model of influence of particle fracture and particle cluster on mechanical properties of metal matrix composites

A general analytical model for a composite with an isotropic matrix and two populations of spherical inclusions is proposed. The method is based on the second order moment of stress for evaluating the homogenised effective stress in the matrix and on the secant moduli concept for the plastic deformation. With Webull's statistical law for the strength of SiCp particles, the model can quantitatively predict the influence of particle fracture on the mechanical properties of PMMCs. Application of the proposed model to the particle cluster shows that the particle cluster has neglected influence on the strain and stress curves of the composite. (C) 1998 Elsevier Science B.V.

Rheological effect on thermocapillary flow of a liquid film jet painted on a moving boundary

In the present paper, a liquid (or melt) film of relatively high temperature ejected from a vessel and painted on the-moving solid film is analyzed by using the second-order fluid model of the non-Newtonian fluid. The thermocapillary flow driven by the temperature gradient on the free surface of a Newtonian liquid film was discussed before. The effect of rheological fluid on thermocapillary flow is considered in the present paper. The analysis is based on the approximations of lubrication theory and perturbation theory. The equation of liquid height and the process of thermal hydrodynamics of the non-Newtonian liquid film are obtained, and the case of weak effect of the rheological fluid is solved in detail.

Numerical investigation of dynamic stall of an oscillating aerofoil

The flow structure around an NACA 0012 aerofoil oscillating in pitch around the quarter-chord is numerically investigated by solving the two-dimensional compressible N-S equations using a special matrix-splitting scheme. This scheme is of second-order accuracy in time and space and is computationally more efficient than the conventional flux-splitting scheme. A 'rigid' C-grid with 149 x 51 points is used for the computation of unsteady flow. The freestream Mach number varies from 0.2 to 0.6 and the Reynolds number from 5000 to 20,000. The reduced frequency equals 0.25-0.5. The basic flow structure of dynamic stall is described and the Reynolds number effect on dynamic stall is briefly discussed. The influence of the compressibility on dynamic stall is analysed in detail. Numerical results show that there is a significant influence of the compressibility on the formation and convection of the dynamic stall vortex. There is a certain influence of the Reynolds number on the flow structure. The average convection velocity of the dynamic stall vortex is approximately 0.348 times the freestream velocity.

Computation of turbulent-flow in a square duct - aspects of the secondary flow and its origin

A numerical study of turbulent flow in a straight duct of square cross-section is made. An order-of-magnitude analysis of the 3-D, time-averaged Navier-Stokes equations resulted in a parabolic form of the Navier-Stokes equations. The governing equations, expressed in terms of a new vector-potential formulation, are expanded as a multi-deck structure with each deck characterized by its dominant physical forces. The resulting equations are solved using a finite-element approach with a bicubic element representation on each cross-sectional plane. The numerical integration along the streamwise direction is carried out with finite-difference approximations until a fully-developed state is reached. The computed results agree well with other numerical studies and compare very favorably with the available experimental data. One important outcome of the current investigation is the interpretation analytically that the driving force of the secondary flow in a square duct comes mainly from the second-order terms of the difference in the gradients of the normal and transverse Reynolds stresses in the axial vorticity equation.

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.

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.