低速二维流动的粒子仿真研究

本文运用信息保存法对低速二维的流动现象进行模拟，考察了低速条件下的有限平板绕流以及微槽道气体流动问题。研究表明：在对低速流动的模拟过程中，运用IP法在能够获得较好的结果的同时，具有比DSMC方法更高的计算效率。     

Efficient discretization of pressure‐correction equations on non‐orthogonal grids

An alternative discretization of pressure‐correction equations within pressure‐correction schemes for the solution of the incompressible Navier–Stokes equations is introduced, which improves the convergence and robustness properties of such schemes for non‐orthogonal grids. As against standard approaches, where the non‐orthogonal terms usually are just neglected, the approach allows for a simplification of the pressure‐correction equation to correspond to 5‐point or 7‐point computational molecules in two or three dimensions, respectively, but still incorporates the effects of non‐orthogonality. As a result a wide range (including rather high values) of underrelaxation factors can be used, resulting in an increased overall performance of the underlying pressure‐correction schemes. Within this context, a second issue of the paper is the investigation of the accuracy to which the pressure‐correction equation should be solved in each pressure‐correction iteration. The scheme is investigated for standard test cases and, in order to show its applicability to practical flow problems, for a more complex configuration of a micro heat exchanger. Copyright © 2003 John Wiley & Sons, Ltd.     

Axisymmetric frictionless contact of functionally graded materials

The main interest of this study is a new method to solve the axisymmetric frictionless contact problem of functionally graded materials (FGMs). Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into a series of sub-layers with shear modulus varying linearly in each sub-layer and continuous at the sub-interfaces. With this model, the axisymmetric frictionless contact problem of a functionally graded coated half-space is investigated. By using the transfer matrix method and Hankel integral transform technique, the problem is reduced to a Cauchy singular integral equation. The contact pressure, contact region and indentation are calculated for various indenters by solving the equations numerically. An erratum to this article can be found at     

Numerical method and models for anti-plane strain of a system with a thin elastic wedge

The paper presents a general method to find asymptotics for a (multi-)wedge system containing a thin wedge. It employs separation of the symmetric and anti-symmetric parts of the boundary displacements and tractions of the wedge. The method is applicable when the angle of the thin wedge turns to zero. A physical interpretation of the derived equations is obtained by using power expansions of non-polynomial functions, which appear after the Mellin transform. We establish that the first term in the expansion of the symmetric part corresponds to shear, while the first term of the anti-symmetric part describes deflection of the wedge axis. Numerical experiments, performed by using a code developed on the basis of the theory, show that using only the first terms of the expansions insignificantly influence accuracy: the approximate results coincide with the exact values of roots to the third significant digit even for the wedge angle of 30°.     

有限元节点编号优化

从分析最简单最优的矩形单元节点编号着手，分析其特点，得出各个节点的相邻单元节点编号总和除以该节点相邻单元数所得的节点商以及相邻单元最小最大节点编号和都为升序排列。总结出一种新型的有限元节点重编号方法，并举例说明该算法能减少有限元的带宽。     

A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.     

Diffusion via splitting and remeshing via merging in vortex methods

The technique of splitting a fat vortex element (with a core width larger than some threshold) into some thin ones in order to fix the convergence problem of the core‐spreading vortex methods is convenient and efficient. In particular, it keeps the method purely Lagrangian. In the present investigation, the splitting process is further viewed as part of the physical diffusion process. A new splitting method in which several weaker child vortices surround a thinned but still strong parent vortex is proposed. It is found that because of the survival of the parent vortex, the error arising from the splitting events can be largely reduced. The computational amount on the other hand is kept reasonably large by merging similar and close‐by vortices. The merging scheme designed herein not only involves fewer restrictions but also allows merging vortices of opposite rotations through the viewpoint of remeshing. The validity and accuracy of these techniques, proposed particularly for simulations undergoing lots of splitting and merging events, are verified by successfully simulating the interactions between two Burgers vortices under an external straining field. Copyright © 2005 John Wiley & Sons, Ltd.     

On the use of the Galerkin method for 3D numerical modelling of the general circulation: the South Atlantic eperiment

A linear three‐dimensional hydrodynamical numerical model, with the application of the Galerkin Method for the vertical dependence, is here presented. The spherical coordinate system is used, in order to allow large‐scale simulations. The equations and mathematical development of the model are shown in detail, together with the boundary and initial conditions, and the sequence of equations' solution. The model is applied to the South Atlantic Ocean, for estimating typical seasonal circulations, and the results are summarized in maps of currents at surface and 1000 m depth, and in transport values of the Brazil Current between 30°S and 40°S. Copyright © 2002 John Wiley & Sons, Ltd.     

Least‐squares finite element method on adaptive grid for PDEs with shocks

We apply the least‐squares finite element method with adaptive grid to nonlinear time‐dependent PDEs with shocks. The least‐squares finite element method is also used in applying the deformation method to generate the adaptive moving grids. The effectiveness of this method is demonstrated by solving a Burgers' equation with shocks. Computational results on uniform grids and adaptive grids are compared for the purpose of evaluation. The results show that the adaptive grids can capture the shock more sharply with significantly less computational time. For moving shock, the adaptive grid moves correctly with the shock. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Meshless simulation of equilibrium swelling/deswelling of pH‐sensitive hydrogels

Hydrogels have been widely used in microelectromechanical systems (MEMS) and Bio‐MEMS devices. In this article, the equilibrium swelling/deswelling of the pH‐stimulus cylindrical hydrogel in the microchannel is studied and simulated by the meshless method. The multi‐field coupling model, called multi‐effect‐coupling pH‐stimulus (MECpH) model, is presented and used to describe the chemical field, electric field, and the mechanical field involved in the problem. The partial differential equations (PDEs) describing these three fields are either nonlinear or coupled together. This multi‐field coupling and high nonlinear characteristics produce difficulties for the conventional numerical methods (e.g., the finite element method or the finite difference method), so an alternative—meshless method is developed to discretize the PDEs, and the efficient iteration technique is adopted to solve the nonlinear problem. The computational results for the swelling/deswelling diameter of the hydrogel under the different pH values are firstly compared with experimental results, and they have a good agreement. The influences of other parameters on the mechanical properties of the hydrogel are also investigated in detail. It is shown that the multi‐field coupling model and the developed meshless method are efficient, stable, and accurate for simulation of the properties of the stimuli‐sensitive hydrogel. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 326–337, 2006