An approach of dynamic mesh adaptation for simulating 3‐dimensional unsteady moving‐immersed‐boundary flows

In this paper, we present an approach of dynamic mesh adaptation for simulating complex 3‐dimensional incompressible moving‐boundary flows by immersed boundary methods. Tetrahedral meshes are adapted by a hierarchical refining/coarsening algorithm. Regular refinement is accomplished by dividing 1 tetrahedron into 8 subcells, and irregular refinement is only for eliminating the hanging points. Merging the 8 subcells obtained by regular refinement, the mesh is coarsened. With hierarchical refining/coarsening, mesh adaptivity can be achieved by adjusting the mesh only 1 time for each adaptation period. The level difference between 2 neighboring cells never exceeds 1, and the geometrical quality of mesh does not degrade as the level of adaptive mesh increases. A predictor‐corrector scheme is introduced to eliminate the phase lag between adapted mesh and unsteady solution. The error caused by each solution transferring from the old mesh to the new adapted one is small because most of the nodes on the 2 meshes are coincident. An immersed boundary method named local domain‐free discretization is employed to solve the flow equations. Several numerical experiments have been conducted for 3‐dimensional incompressible moving‐boundary flows. By using the present approach, the number of mesh nodes is reduced greatly while the accuracy of solution can be preserved.     

为什么要引入溶剂的渗透因子

阐明了引入溶剂渗透因子的必要性,较详细地介绍了两种渗透因子及它们间的关系。     

A short note on the entrainment and exit boundary conditions

An attempt is made to find out the suitable entrainment and exit boundary conditions in laminar flow situations. Streamfunction vorticity formulation of the Navier–Stokes equations are solved by ADI method. Two‐dimensional laminar plane wall jet flow is used to test different forms of the boundary conditions. Results are compared with the experimental and similarity solution and the proper boundary condition is suggested. The Kind 1 boundary condition is recommended. It consists of zero first derivative condition for velocity variable and for streamfunction equation, mixed derivative at the entrainment and exit boundaries. Copyright © 2005 John Wiley & Sons, Ltd.     

Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

The destabilizing effect of boundary slip on Bénard convection

We investigate the influence of slip boundary conditions on the onset of Bénard convection in an infinite fluid layer. It is shown that the critical Rayleigh number is a decreasing function of the slip length, and therefore boundary slip is seen to have a destabilizing effect. Chebyshev‐tau and compound matrix formulations for solving the eigenvalue problem are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Banach空间二阶非线性微分方程的弱Carathéodory解

本文定义了二阶微分方程的弱 Carathéodory解 ,在不涉及紧型条件的情形下 ,直接用迭代法证明了 Banach空间二阶非线性常微分方程两点边值问题存在唯一解 ,并给出逼近解迭代序列的误差估计 ,对周期边值问题得到类似的结果     

Method of lines with boundary elements for 1‐D transient diffusion‐reaction problems

Time‐dependent differential equations can be solved using the concept of method of lines (MOL) together with the boundary element (BE) representation for the spatial linear part of the equation. The BE method alleviates the need for spatial discretization and casts the problem in an integral format. Hence errors associated with the numerical approximation of the spatial derivatives are totally eliminated. An element level local cubic approximation is used for the variable at each time step to facilitate the time marching and the nonlinear terms are represented in a semi‐implicit manner by a local linearization at each time step. The accuracy of the method has been illustrated on a number of test problems of engineering significance. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Roles of SiH4 and SiF4 in growth and structural changes of poly-Si films

The structural properties of polycrystalline silicon films, prepared by plasma enhanced chemical vapor deposition system, with different flow rates of SiH₄/SiF₄ mixtures at 300 °C were investigated. This study indicates that the low hydrogen coverage on the growing surface, under optimum fluorine radicals, will be leaded to an improvement of crystallized area as compared with case of high hydrogen coverage surface. Moreover, the studies of the role of SiH₄ and SiF₄ radicals show that the SiH₄ radicals are important in the nucleation and growth of grains. However, SiF₄ radicals are effective in the structural change of grain boundaries regions and by this way, in the present system, establish the growth of grains under the dominant 〈1 1 0〉 direction. The stress investigation indicates that addition of high flow rate of SiF₄ in amorphous film, results in the nearly stress free films. Finally, we found that the changes in g-value reflect the changes in the intrinsic compressive and tensile stress in the both polycrystalline and amorphous silicon films.     

超临界流体触媒作用下石墨-金刚石转变

 研究了超临界流体CO₂在石墨-金刚石转变中的触媒作用。实验中，采用Ag₂O作为流体触媒的先驱材料，在7.7 GPa压力下，Ag₂O在1 200 ℃分解成Ag和O₂，O₂与石墨套管在高温高压下反应形成CO₂超临界流体。研究结果表明，在7.7 GPa和1 500 ℃以上温度条件下，石墨在CO₂流体触媒的作用下可转变为金刚石晶体，在1 500～1 700 ℃温度范围内合成出的金刚石具有完好的八面体形貌，与天然金刚石的生长特征非常相似。