Coupling of free surface and groundwater flows

In this paper we present some preliminary results about the coupling of shallow water equations for free surface flows and Darcy equation for groundwater flows. A suitable set of interface conditions is discussed: the Beavers and Joseph formula for the bottom stress is used. An iterative algorithm to solve the coupled problem is proposed and some numerical results are presented.     

Stabilization Methods for Spectral Element Computations of Incompressible Flows

A stabilization method for the spectral element computation of incompressible flow problems is investigated. It is based on a filtering procedure which consists in filtering the velocity field by a spectral vanishing Helmholtz-type operator at each time step. Relationship between this filtering procedure and SVV-stabilization method, introduced recently in [JCP, 2004, 196(2), p680], is established. A number of numerical examples are presented to show the accuracy and stabilization capability of the method.     

Implicit–Explicit Schemes for BGK Kinetic Equations

In this work a new class of numerical methods for the BGK model of kinetic equations is presented. In principle, schemes of any order of accuracy in both space and time can be constructed with this technique. The methods proposed are based on an explicit–implicit time discretization. In particular the convective terms are treated explicitly, while the source terms are implicit. In this fashion even problems with infinite stiffness can be integrated with relatively large time steps. The conservation properties of the schemes are investigated. Numerical results are shown for schemes of order 1, 2 and 5 in space, and up to third-order accurate in time.     

Viscoelastic State for the Solutions of Spherically Symmetric Viscoelasticity Problems

The solutions of the spherically symmetric, linear, isothermal, and transient viscoelasticity problems via reciprocity theorem have been investigated for a specific material. The integral form of stress–strain relations has been used. The Laplace transform of a viscoelastic state, which is necessary for the integral equation arising as a result of reciprocity theorem, has been derived. This integral equation has been solved by Laplace transform. A sample problem has been solved to test the presented formulation. A numerical application of the analytic solution of this problem has been given.     

一种复合型数值方法的改进与算法实现

文中对原有的复合方法进行了改进,提出了一种新的复合型数值方法,这种方法解决了原方法求解含源项的问题时存在数值震荡的缺点。为了验证算法的有效性,对计算流体力学中一阶浅水波问题的两个典型的算例进行了数值模拟,在解的光滑性、锐利性等几个方面做了分析,结果表明改进后的方法既消除了数值震荡现象,又能准确地逼近真解,求出锐利的间断,是一种行之有效的数值方法。     

泥沙输移矢量式及其应用

本文基于Parker等人泥沙输移矢量式的分析和研究，增加了上举力的作用，并考虑了泥沙颗粒受力非平衡条件，采用牛顿力学的矢量分析法，推导了泥沙输移的通用矢量式。经对公式简化，且用Ikeda和笔者试验资料的验证，证明计算值与试验值吻合较好。     

Affine Dependence Classification for Communications Minimization

This paper introduces results on placement and communications minimization for systems of affine recurrence equations. We show how to classify the dependences according to the number and nature of communications they may result in. We give both communication-free conditions and conditions for an efficient use of broadcast or neighbor-to-neighbor communication primitives. Since the dependences of a problem can generally not be all communication-free, we finally introduce a heuristic to globally minimize the communications based on the classification of dependences.     

ON THE MODELLING OF PARTICLE-BODY INTERACTIONS IN STOKES FLOWS INVOLVING A SPHERE AND CIRCULAR DISC OR A TORUS AND CIRCULAR CYLINDER USING POINT SINGULARITIES

Few exact solutions of the Stokes equations are known, even for steady or quasi-steady flows, involving finite sized bodies, and numerical techniques generally have to be resorted to for finding solutions. However, quite effective modelling of flows involving complicated boundary geometries is possible using the three-dimensional Stokeslet and rotlet point singularities. Two problems are studied in detail. In the first example, exact solutions for the three-dimensional Stokeslet and rotlet placed axisymmetrically along the axis of a circular disc are found and combined with Brenner's first order interaction formulae to determine the effect of the presence of the disc on the force and torque acting on a particle whose dimensions are small compared with its distance from the disc. The results are compared with those of a full numerical integration of the Stokes equations for a sphere translating towards a disc. In the second example, Brenner's first order wall correction theory is applied to the motion of a particle in a circular cylinder using the exact solutions for a torus translating or rotating in isolation. The theoretical predictions for the drag on a torus settling symmetrically in a circular cylinder are compared with those determined experimentally.     

A semi-analytic boundary element method for parabolic problems

A new semi-analytic solution method is proposed for solving linear parabolic problems using the boundary element method. This method constructs a solution as an eigenfunction expansion using separation of variables. The eigenfunctions are determined using the dual reciprocity boundary element method. This separation of variables-dual reciprocity method (SOV-DRM) allows a solution to be determined without requiring either time-stepping or domain discretisation. The accuracy and computational efficiency of the SOV-DRM is found to improve as time increases. These properties make the SOV-DRM an attractive technique for solving parabolic problems.     

Flows in rectangular microchannels driven by capillary force and gravity

Capillary phenomena was studied and discussed by the scholars about 200 years ago, but the progress was slow due to the limited equipment and manufacture precision of the microchannel. In recent years, because of the rapid development of MEMS and micromachining, many applications of the capillary flow is widely developing in some modern processes, such as underfilling of flip chip, flow in microfluidic chip or biochip, and a variety of other fields.