A Pseudo-Stokes Mesh Motion Algorithm

This work presents a moving mesh methodology based on the solution of a pseudo flow problem. The mesh motion is modeled as a pseudo Stokes problem solved by an explicit finite element projection method. The mesh quality requirements are satisfied by employing a null divergent velocity condition. This methodology is applied to triangular unstructured meshes and compared to well known approaches such as the ones based on diffusion and pseudo structural problems. One of the test cases is an airfoil with a fully meshed domain. A specific rotation velocity is imposed as the airfoil boundary condition. The other test is a set of two cylinders that move toward each other. A mesh quality criterion is employed to identify critically distorted elements and to evaluate the performance of each mesh motion approach. The results obtained for each test case show that the pseudo-flow methodology produces satisfactory meshes during the moving process.     

NESTED GRID METHODS FOR AN OCEAN MODEL: A COMPARATIVE STUDY

In this paper a comparison is carried out between three correction methods for multigrid local mesh refinement in oceanic applications: FIC, LDC and the direct method (DM) proposed by Spall and Holland. This study is based on a nested primitive equation model developed by Laugier on the basis of the code OPA (LODYC). The external barotropic problem is solved using any of the three local grid correction algorithms yielding an interactive nested grid model. The non-linear elliptic equation for the barotropic streamfunction tendency is solved on two nested grids, called the global and the zoom grid, that interact between themselves. The zoom grid is entirely embedded within the global domain with a horizontal grid step ratio of 3:1. The computation on the global grid supplies the boundary conditions for the zoom grid region and the fine grid fields are used to correct the global coarse solution. The three local correction methods are tested on two problems relevant to oceanic circulation phenomena proposed by Spall and Holland: a barotropic modon and an anticyclonic vortex. The results show that the nesting technique is a very efficient way to solve these problems in terms of a gain in precision compared with the required CPU time. The two-domain model with local mesh refinement allows one both to manage effectively the open boundary conditions for the local grid and to correct the global solution thanks to the zoom solution. In the case of the modon propagation the three local correction methods provide approximately the same results. For the baroclinic vortex it appears that the two iterative methods are more efficient than the direct one.     

On the construction of interpolation mesh surfaces

A method for constructing two-dimensional interpolation mesh functions is proposed that is more flexible than the classical cubic spline method because it makes it possible to construct interpolation surfaces that fit the given function at specified points by varying certain parameters. The method is relatively simple and is well suited for practical implementation.     

NUMERICAL SOLUTION OF THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON THE CURVILINEAR HALF-STAGGERED MESH NUMERICAL SOLUTION OF THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON THE CURVILINEAR HALF-STAGGERED MESH

1. IntroductionLet us consider the unsteady incompressible Navier--Stokes equatinns (qSE)on a twedimensional reginn n with boundary On. Here w is the velocity vector; interms of its Caxtesian,components w = (u, v)"; p is the Pressure. The initial conditionis given assatisheng (1.2). We are concerned mainly with the solid wall boundary conditiollbut we will also briefly discuss the outflow boundary conditinn. Note the convectionterm can ajBo be written in conservative form (w' glad)w = div(…     

Uniform Convergence Analysis of Finite Difference Scheme for Singularly Perturbed Delay Differential Equation on an Adaptively Generated Grid

<正>Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind scheme for singularly perturbed differential-difference equation on a grid which is formed by equidistributing arc-length monitor function.It is shown that the discrete solution obtained converges uniformly with respect to the perturbation parameter.Numerical experiments illustrate in practice the result of convergence proved theoretically.     

A fast and efficient mesh segmentation method based on improved region growing

Mesh segmentation is one of the important issues in digital geometry processing. Region growing method has been proven to be a efficient method for 3D mesh segmentation. However, in mesh segmentation, feature line extraction algorithm is computationally costly, and the over-segmentation problem still exists during region merging processing. In order to tackle these problems, a fast and efficient mesh segmentation method based on improved region growing is proposed in this paper. Firstly, the dihedral angle of each non-boundary edge is defined and computed simply, then the sharp edges are detected and feature lines are extracted. After region growing process is finished, an improved region merging method will be performed in two steps by considering some geometric criteria. The experiment results show the feature line extraction algorithm can obtain the same geometric information fast with less computational costs and the improved region merging method can solve over-segmentation well.     

ADAPTIVE QUADRILATERAL AND HEXAHEDRAL FINITE ELEMENT METHODS WITH HANGING NODES AND CONVERGENCE ANALYSIS

In this paper we study the convergence of adaptive finite element methods for the gen- eral non-attine equivalent quadrilateral and hexahedral elements on 1-irregular meshes with hanging nodes. Based on several basic ingredients, such as quasi-orthogonality, estimator reduction and D6fler marking strategy, convergence of the adaptive finite element methods for the general second-order elliptic partial equations is proved. Our analysis is effective for all conforming Qm elements which covers both the two- and three-dimensional cases in a unified fashion.