An algebraic multigrid solver for Navier-Stokes problems

An efficient numerical method is presented for solving the equations of motion for viscous fluids. The equations are discretized on the basis of unstructured finite element meshes and then solved by direct iteration. Advective fluxes are temporarily fixed at each iteration to provide a linearized set of coupled equations which are then also solved by iteration using a fully implicit algebraic multigrid (AMG) scheme. A rapid convergence to machine accuracy is achieved that is almost mesh-independent. The scaling of computing time with mesh size is therefore close to the optimum.     

A multigrid method with unstructured adaptive grids for steady Euler equations

The multigrid method based on multi-stage Jacobi relaxation, earlier developed by the authors for structured grid calculations with Euler equations, is extended to unstructured grid applications. The meshes are generated with Delaunay triangulation algorithms and are adapted to the flow solution.     

Resolution of computational aeroacoustics problems on unstructured grids with a higher-order finite volume scheme

Computational fluid dynamics (CFD) has become increasingly used in the industry for the simulation of flows. Nevertheless, the complex configurations of real engineering problems make the application of very accurate methods that only work on structured grids difficult. From this point of view, the development of higher-order methods for unstructured grids is desirable. The finite volume method can be used with unstructured grids, but unfortunately it is difficult to achieve an order of accuracy higher than two, and the common approach is a simple extension of the one-dimensional case. The increase of the order of accuracy in finite volume methods on general unstructured grids has been limited due to the difficulty in the evaluation of field derivatives. This problem is overcome with the application of the Moving Least Squares (MLS) technique on a finite volume framework. In this work we present the application of this method (FV-MLS) to the solution of aeroacoustic problems.     

A nonlinear PDE model for reconstructing a regular surface from sampled data using a level set formulation on triangular meshes

In this paper, we propose a nonlinear PDE model for reconstructing a regular surface from sampled data. At first, we show the existence and the uniqueness of a viscosity solution to this problem. Then we propose a numerical scheme for solving the nonlinear level set equation on unstructured triangulations adapted to the data sample. We show the consistency of this scheme. In addition, we show how to compute nodewise first and second order derivatives. Some application examples of curve or surface construction are provided to illustrate the potential and to demonstrate the accuracy of this method.     

Anti-diffusion method for interface steepening in two-phase incompressible flow

In this paper, we present a method for obtaining sharp interfaces in two-phase incompressible flows by an anti-diffusion correction, that is applicable in a straight-forward fashion for the improvement of two-phase flow solution schemes typically employed in practical applications. The underlying discretization is based on the volume-of-fluid (VOF) interface-capturing method on unstructured meshes. The key idea is to steepen the interface, independently of the underlying volume-fraction transport equation, by solving a diffusion equation with reverse time, i.e. an anti-diffusion equation, after each advection time step of the volume fraction. As the solution of the anti-diffusion equation requires regularization, a limiter based on the directional derivative is developed for calculating the gradient of the volume fraction. This limiter ensures the boundedness of the volume fraction. In order to control the amount of anti-diffusion introduced by the correction algorithm we propose a suitable stopping criterion for interface steepening. The formulation of the limiter and the algorithm for solving the anti-diffusion equation are applicable to 3-dimensional unstructured meshes. Validation computations are performed for passive advection of an interface, for 2-dimensional and 3-dimensional rising-bubbles, and for a rising drop in a periodically constricted channel. The results demonstrate that sharp interfaces can be recovered reliably. They show that the accuracy is similar to or even better than that of level-set methods using comparable discretizations for the flow and the level-set evolution. Also, we observe a good agreement with experimental results for the rising drop where proper interface evolution requires accurate mass conservation.     

The multi-dimensional limiters for solving hyperbolic conservation laws on unstructured grids

Novel limiters based on the weighted average procedure are developed for finite volume methods solving multi-dimensional hyperbolic conservation laws on unstructured grids. The development of these limiters is inspired by the biased averaging procedure of Choi and Liu [10]. The remarkable features of the present limiters are the new biased functions and the weighted average procedure, which enable the present limiter to capture strong shock waves and achieve excellent convergence for steady state computations. The mechanism of the developed limiters for eliminating spurious oscillations in the vicinity of discontinuities is revealed by studying the asymptotic behavior of the limiters. Numerical experiments for a variety of test cases are presented to demonstrate the superior performance of the proposed limiters.     

A parallel adaptive mesh approach for flowfields with shock waves

A simple and efficient implementation of Adaptive Mesh Refinement (AMR) with distributed-memory approach is presented. Introducing a lookup table including grid connectivity information and simplified algorithms for AMR, the procedures for reconstructing adaptive grids are carried out in parallel, with local data to a large extent. A simple static load-balancing scheme is adopted, and the grids are not repartitioned and no data redistribution is performed. A numerical example on two different parallel computers shows that the proposed implementation of AMR is effective to reduce the computational time for unsteady flows with shock waves. Received 23 October 2000 / Accepted 30 March 2001 Published online 11 June 2002     

Shock waves from an open-ended shock tube with different shapes

A new method for decreasing the attenuation of a shock wave emerging from an open-ended shock tube exit into a large free space has been developed to improve the shock wave technique for cleaning deposits on the surfaces in industrial equipments by changing the tube exit geometry. Three tube exits (the simple tube exit, a tube exit with ring and a coaxial tube exit) were used to study the propagation processes of the shock waves. The detailed flow features were experimentally investigated by use of a two-dimensional color schlieren method and by pressure measurements. By comparing the results for different tube exits, it is shown that the expansion of the shock waves near the mouth can be restricted by using the tube exit with ring or the coaxial tube exit. Thus, the attenuation of the shock waves is reduced. The time histories of overpressure have illustrated that the best results are obtained for the coaxial tube exit. But the pressure signals for the tube exit with ring showed comparable results with the advantage of a relatively simple geometry. The flow structures of diffracting shock waves have also been simulated by using an upwind finite volume scheme based on a high order extension of Godunov's method as well as an adaptive unstructured triangular mesh refinement/unrefinement algorithm. The numberical results agree remarkably with the experimental ones.     

A reduced domain strategy for local mesh movement application in unstructured grids

Automatic control of mesh movement is mandatory in many fluid flow and fluid-solid interaction problems. This paper presents a new strategy, called reduced domain strategy (RDS), which enhances the efficiency of node connectivity-based mesh movement methods and moves the unstructured grid locally and effectively. The strategy dramatically reduces the grid computations by dividing the unstructured grid into two active and inactive zones. After any local boundary movement, the grid movement is performed only within the active zone. To enhance the efficiency of our strategy, we also develop an automatic mesh partitioning scheme. This scheme benefits from a new quasi-structured mesh data ordering, which determines the boundary of active zone in the original unstructured grid very easily. Indeed, the new partitioning scheme eliminates the need for sequential reordering of the original unstructured grid data in different mesh movement applications. We choose the spring analogy method and apply our new strategy to perform local mesh movements in two boundary movement problems including a multi-element airfoil with moving slat or deforming main body section. We show that the RDS is robust and cost effective. It can be readily employed in different node connectivity-based mesh movement methods. Indeed, the RDS provides a flexible local grid deformation tool for moving grid applications.