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.     

Inequalities for the Moments and Distribution of the Ladder Height of a Random Walk

We obtain upper bounds for the tail distribution of the first nonnegative sum of a random walk and for the moments of the overshoot over an arbitrary nonnegative level if the expectation of jumps is positive and close to zero. In addition, we find an estimate for the expectation of the first ladder epoch.     

Coupled intervals in the discrete calculus of variations: necessity and sufficiency

In this work we study nonnegativity and positivity of a discrete quadratic functional with separately varying endpoints. We introduce a notion of an interval coupled with 0, and hence, extend the notion of conjugate interval to 0 from the case of fixed to variable endpoint(s). We show that the nonnegativity of the discrete quadratic functional is equivalent to each of the following conditions: The nonexistence of intervals coupled with 0, the existence of a solution to Riccati matrix equation and its boundary conditions. Natural strengthening of each of these conditions yields a characterization of the positivity of the discrete quadratic functional. Since the quadratic functional under consideration could be a second variation of a discrete calculus of variations problem with varying endpoints, we apply our results to obtain necessary and sufficient optimality conditions for such problems. This paper generalizes our recent work in [R. Hilscher, V. Zeidan, Comput. Math. Appl., to appear], where the right endpoint is fixed.