动态混合网格生成及隐式非定常计算方法

建立了一种基于动态混合网格的非定常数值计算方法. 混合网格由贴体的四边形网格、外场 的多层次矩形网格和中间的三角形网格构成. 当物体运动时，贴体四边形网格随物体运动而 运动，而外场的矩形网格保持静止，中间的三角形网格随之变形；当物体运动位移较大，导 致三角形网格的质量降低，甚至导致网格相交时，在局部重新生成网格. 新网格上的物理量 由旧网格上的物理量插值而得. 为了提高计算效率，采用了双时间步和子迭代相结合的隐式 有限体积格式计算非定常Navier-Stokes方程. 子迭代采用高效的块LU-SGS方法. 利用该 方法数值模拟了NACA0012振荡翼型的无黏和黏性绕流，得到了与实验和他人计算相当一致 的结果.     

A hybrid building‐block and gridless method for compressible flows

A hybrid building‐block Cartesian grid and gridless method is presented to compute unsteady compressible flows for complex geometries. In this method, a Cartesian mesh based on a building‐block grid is used as a baseline mesh to cover the computational domain, while the boundary surfaces are represented using a set of gridless points. This hybrid method combines the efficiency of a Cartesian grid method and the flexibility of a gridless method for the complex geometries. The developed method is used to compute a number of test cases to validate the accuracy and efficiency of the method. The numerical results obtained indicate that the use of this hybrid method leads to a significant improvement in performance over its unstructured grid counterpart for the time‐accurate solution of the compressible Euler equations. An overall speed‐up factor from six to more than one order of magnitude and a saving in storage requirements up to one order of magnitude for all test cases in comparison with the unstructured grid method are demonstrated. Copyright © 2008 John Wiley & Sons, Ltd.     

Orthogonal grid generation for Navier–Stokes computations

Body conforming orthogonal grids were generated using a fast hyperbolic method for aerofoils, and were used to solve the Navier–Stokes equation in the generalized orthogonal system for the first time for time accurate simulation of incompressible flow. For grid generation, the Beltrami equation and the definition equation for the orthogonality are solved using a finite difference method. The grids generated around aerofoils by this method have better orthogonality than the results published by earlier investigators. The Navier–Stokes equation at Reynolds numbers of 3000 and 35 000 for NACA 0012 and NACA 0015 respectively, have been solved as an application. The obtained results match quite well with the corresponding experimental results. © 1998 John Wiley & Sons, Ltd.     

Convergence of steady and unsteady formulations for inviscid hovering rotor solutions

An upwind Euler solver is presented, and applied to multibladed lifting hovering rotor flow. These flows can be simulated as a steady case, in a blade‐fixed rotating co‐ordinate system. However, forward flight simulation will always require an unsteady solution. Hence, as a stepping stone in the development of a forward flight simulation tool, both explicit steady and implicit unsteady simulations of the same hovering case are presented. Convergence of the two approaches is examined and compared, in terms of residual history, cost, and solution evolution, as a means of both validating the unsteady formulation and considering implications for forward flight simulation. Consideration of the solution evolution and wake capturing shows that for hovering rotor cases, the unsteady and steady solutions are the same, but the unsteady solution is more expensive in terms of CPU time. It is also shown that for hover, the fewer real time‐steps taken per revolution the more efficient the implicit scheme is. However, this is a characteristic of the case, which results in smooth solution variation between time steps. It is also demonstrated that for rotary flow simulation, the global residual is not a useful quantity to assess convergence. The residual reaches a very low (constant in the implicit case) value while the solution is still evolving. Copyright © 2003 John Wiley & Sons, Ltd.     

A sliding interface method for unsteady unstructured flow simulations

The objective of this work is to develop a sliding interface method for simulations involving relative grid motion that is fast and efficient and involves no grid deformation, remeshing, or hole cutting. The method is implemented into a parallel, node‐centred finite volume, unstructured viscous flow solver. The rotational motion is accomplished by rigidly rotating the subdomain representing the moving component. At the subdomain interface boundary, the faces along the interface are extruded into the adjacent subdomain to create new volume elements forming a one‐cell overlap. These new volume elements are used to compute a flux across the subdomain interface. An interface flux is computed independently for each subdomain. The values of the solution variables and other quantities for the nodes created by the extrusion process are determined by linear interpolation. The extrusion is done so that the interpolation will maintain information as localized as possible. The grid on the interface surface is arbitrary. The boundary between the two subdomains is completely independent from one another; meaning that they do not have to connect in a one‐to‐one manner and no symmetry or pattern restrictions are placed on the grid. A variety of numerical simulations were performed on model problems and large‐scale applications to examine conservation of the interface flux. Overall solution errors were found to be comparable to that for fully connected and fully conservative simulations. Excellent agreement is obtained with theoretical results and results from other solution methodologies. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical solution of the unsteady Euler equations for airfoils using approximate boundary conditions

This paper presents an efficient numerical method for solving the unsteady Euler equations on stationary rectilinear grids. Boundary conditions on the surface of an airfoil are implemented by using their first-order expansions on the mean chord line. The method is not restricted to flows with small disturbances since there are no restrictions on the mean angle of attack of the airfoil. The mathematical formulation and the numerical implementation of the wall boundary conditions in a fully implicit time-accurate finite-volume Euler scheme are described. Unsteady transonic flows about an oscillating NACA 0012 airfoil are calculated. Computational results compare well with Euler solutions by the full boundary conditions on a body-fitted curvilinear grid and published experimental data. This study establishes the feasibility for computing unsteady fluid-structure interaction problems, where the use of a stationary rectilinear grid offers substantial advantages in saving computer time and program design since it does not require the generation and implementation of time-dependent body-fitted grids.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid

This paper presents a finite element method for incompressible multiphase flows with capillary interfaces based on a (formally) second‐order projection scheme. The discretization is on a fixed Eulerian grid. The fluid phases are identified and advected using a level set function. The grid is temporarily adapted around the interfaces in order to maintain optimal interpolations accounting for the pressure jump and the discontinuity of the normal velocity derivatives. The least‐squares method for computing the curvature is used, combined with piecewise linear approximation to the interface. The time integration is based on a formally second order splitting scheme. The convection substep is integrated over an Eulerian grid using an explicit scheme. The remaining generalized Stokes problem is solved by means of a formally second order pressure‐stabilized projection scheme. The pressure boundary condition on the free interface is imposed in a strong form (pointwise) at the pressure‐computation substep. This allows capturing significant pressure jumps across the interface without creating spurious instabilities. This method is simple and efficient, as demonstrated by the numerical experiments on a wide range of free‐surface problems. Copyright © 2004 John Wiley & Sons, Ltd.     

A fast nested multi‐grid viscous flow solver for adaptive Cartesian/Quad grids

A nested multi‐grid solution algorithm has been developed for an adaptive Cartesian/Quad grid viscous flow solver. Body‐fitted adaptive Quad (quadrilateral) grids are generated around solid bodies through ‘surface extrusion’. The Quad grids are then overlapped with an adaptive Cartesian grid. Quadtree data structures are employed to record both the Quad and Cartesian grids. The Cartesian grid is generated through recursive sub‐division of a single root, whereas the Quad grids start from multiple roots—a forest of Quadtrees, representing the coarsest possible Quad grids. Cell‐cutting is performed at the Cartesian/Quad grid interface to merge the Cartesian and Quad grids into a single unstructured grid with arbitrary cell topologies (i.e., arbitrary polygons). Because of the hierarchical nature of the data structure, many levels of coarse grids have already been built in. The coarsening of the unstructured grid is based on the Quadtree data structure through reverse tree traversal. Issues arising from grid coarsening are discussed and solutions are developed. The flow solver is based on a cell‐centered finite volume discretization, Roe's flux splitting, a least‐squares linear reconstruction, and a differentiable limiter developed by Venkatakrishnan in a modified form. A local time stepping scheme is used to handle very small cut cells produced in cell‐cutting. Several cycling strategies, such as the saw‐tooth, W‐ and V‐cycles, have been studies. The V‐cycle has been found to be the most efficient. In general, the multi‐grid solution algorithm has been shown to greatly speed up convergence to steady state—by one to two orders. Copyright © 2000 John Wiley & Sons, Ltd.     

An adaptive wavelet-collocation method for shock computations

A simple and robust method for solving hyperbolic conservation equations based on the adaptive wavelet-collocation method, which uses a dynamically adaptive grid, are presented. The method utilises natural ability of wavelet analysis to sense localised structures and is based on analysis of wavelet coefficients on the finest level of resolution to create a discontinuity locator function Φ. Using this function, an artificial viscous term is explicitly added in the needed regions using a localised numerical viscosity that ensures the positivity and TVD non-linear stability conditions. Once the wavelet coefficients on the finest level of resolution are below the error threshold parameter ?, the artificial viscosity is shut off and any remaining physical waves are free to propagate undamped. Multiple examples in one and two dimensions are presented to demonstrate the method's robustness, simplicity and ease of extending to more complex problems.     

A preconditioning technique for steady Euler solutions

This paper presents a preconditioning technique for solving a two‐dimensional system of hyperbolic equations. The main attractive feature of this approach is that, unlike a technique based on simply extending the solver for a one‐dimensional hyperbolic system, convergence and stability analysis can be investigated. This method represents a genuine numerical algorithm for multi‐dimensional hyperbolic systems. In order to demonstrate the effectiveness of this approach, applications to solving a two‐dimensional system of Euler equations in supersonic flows are reported. It is shown that the Lax–Friedrichs scheme diverges when applied to the original Euler equations. However, convergence is achieved when the same numerical scheme is employed using the same CFL number to solve the equivalent preconditioned Euler system. Copyright © 1999 John Wiley & Sons, Ltd.     

Study of a staggered fourth‐order compact scheme for unsteady incompressible viscous flows

The objective of this paper is the development and assessment of a fourth‐order compact scheme for unsteady incompressible viscous flows. A brief review of the main developments of compact and high‐order schemes for incompressible flows is given. A numerical method is then presented for the simulation of unsteady incompressible flows based on fourth‐order compact discretization with physical boundary conditions implemented directly into the scheme. The equations are discretized on a staggered Cartesian non‐uniform grid and preserve a form of kinetic energy in the inviscid limit when a skew‐symmetric form of the convective terms is used. The accuracy and efficiency of the method are demonstrated in several inviscid and viscous flow problems. Results obtained with different combinations of second‐ and fourth‐order spatial discretizations and together with either the skew‐symmetric or divergence form of the convective term are compared. The performance of these schemes is further demonstrated by two challenging flow problems, linear instability in plane channel flow and a two‐dimensional dipole–wall interaction. Results show that the compact scheme is efficient and that the divergence and skew‐symmetric forms of the convective terms produce very similar results. In some but not all cases, a gain in accuracy and computational time is obtained with a high‐order discretization of only the convective and diffusive terms. Finally, the benefits of compact schemes with respect to second‐order schemes is discussed in the case of the fully developed turbulent channel flow. Copyright © 2008 John Wiley & Sons, Ltd.     

A large-scale parallel hybrid grid generation technique for realistic complex geometry

High-Performance Computing (HPC) systems and Computational Fluid Dynamics (CFD) have made significant progress in recent years; however, as the basis of the large-scale parallel computing, the massive grid generation of billions of cells has become a bottleneck problem. In this study, a parallel grid generation technique is proposed to generate large-scale mixed grids with arbitrary cell types and scales. The basic idea of our method is analogous to the global mesh refinement technique. An initial coarse grid with arbitrary cell types is regarded as a background mesh which is partitioned into subzones, and subzones are assigned onto different CPU cores. After the cells and faces in each subzone are split, the inserted new points of the solid wall are projected onto the original CAD entities to preserve the geometry accurately. Finally, the tangled cells caused by the projection in the boundary layer are untangled by a local Radial Basis Function mesh deformation technique. Furthermore, a parallel partition approach and an efficient wall distance computing technique for massive grids are developed also to shorten the preprocessing time. The tests show that the preprocessing efficiency has been increased by two or three orders compared with traditional methods. Billions of grids are generated for the AIAA JSM high-lift model and the Chinese CHN-T1 transport model to test the ability of the parallel grid generation technique. The maximum scale up to 19 billion mixed elements is generated using 16 384 CPU cores in parallel, and the mesh quality is acceptable for CFD simulations.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid: 3D validation

A three‐dimensional finite element method for incompressible multiphase flows with capillary interfaces is developed based on a (formally) second‐order projection scheme. The discretization is on a fixed (Eulerian) reference grid with an edge‐based local h‐refinement in the neighbourhood of the interfaces. The fluid phases are identified and advected using the level‐set function. The reference grid is then temporarily reconnected around the interface to maintain optimal interpolations accounting for the singularities of the primary variables. Using a time splitting procedure, the convection substep is integrated with an explicit scheme. The remaining generalized Stokes problem is solved by means of a pressure‐stabilized projection. This method is simple and efficient, as demonstrated by a wide range of difficult free‐surface validation problems, considered in the paper. Copyright © 2005 John Wiley & Sons, Ltd.     

An improved unstructured WENO method for compressible multi‐fluids flow

An improved high‐order accurate WENO finite volume method based on unstructured grids for compressible multi‐fluids flow is proposed in this paper. The third‐order accuracy WENO finite volume method based on triangle cell is used to discretize the governing equations. To have higher order of accuracy, the P1 polynomial is reconstructed firstly. After that, the P2 polynomial is reconstructed from the combination of the P1. The reconstructed coefficients are calculated by analytical form of inverse matrix rather than the numerical inversion. This greatly improved the efficiency and the robustness. Four examples are presented to examine this algorithm. Numerical results show that there is no spurious oscillation of velocity and pressure across the interface and high‐order accurate result can be achieved. Copyright © 2016 John Wiley & Sons, Ltd.     

Development of circular function‐based gas‐kinetic scheme (CGKS) on moving grids for unsteady flows through oscillating cascades

In this paper, the circular function‐based gas‐kinetic scheme (CGKS), which was originally developed for simulation of flows on stationary grids, is extended to solve moving boundary problems on moving grids. Particularly, the unsteady flows through oscillating cascades are our major interests. The main idea of the CGKS is to discretize the macroscopic equations by the finite volume method while the fluxes at the cell interface are evaluated by locally reconstructing the solution of the continuous Boltzmann Bhatnagar–Gross–Krook equation. The present solver is based on the fact that the modified Boltzmann equation, which is expressed in a moving frame of reference, can recover the corresponding macroscopic equations with Chapman–Enskog expansion analysis. Different from the original Maxwellian function‐based gas‐kinetic scheme, in improving the computational efficiency, a simple circular function is used to describe the equilibrium state of distribution function. Considering that the concerned cascade oscillating problems belong to cases that the motion of surface boundary is known a priori, the dynamic mesh method is suitable and is adopted in the present work. In achieving the mesh deformation with high quality and efficiency, a hybrid dynamic mesh method named radial basic functions‐transfinite interpolation is presented and applied for cascade geometries. For validation, several numerical test cases involving a wide range are investigated. Numerical results show that the developed CGKS on moving grids is well applied for cascade oscillating flows. And for some cases where nonlinear effects are strong, the solution accuracy could be effectively improved by using the present method. Copyright © 2017 John Wiley & Sons, Ltd.     

A robust multi‐grid pressure‐based algorithm for multi‐fluid flow at all speeds

This paper reports on the implementation and testing, within a full non‐linear multi‐grid environment, of a new pressure‐based algorithm for the prediction of multi‐fluid flow at all speeds. The algorithm is part of the mass conservation‐based algorithms (MCBA) group in which the pressure correction equation is derived from overall mass conservation. The performance of the new method is assessed by solving a series of two‐dimensional two‐fluid flow test problems varying from turbulent low Mach number to supersonic flows, and from very low to high fluid density ratios. Solutions are generated for several grid sizes using the single grid (SG), the prolongation grid (PG), and the full non‐linear multi‐grid (FMG) methods. The main outcomes of this study are: (i) a clear demonstration of the ability of the FMG method to tackle the added non‐linearity of multi‐fluid flows, which is manifested through the performance jump observed when using the non‐linear multi‐grid approach as compared to the SG and PG methods; (ii) the extension of the FMG method to predict turbulent multi‐fluid flows at all speeds. The convergence history plots and CPU‐times presented indicate that the FMG method is far more efficient than the PG method and accelerates the convergence rate over the SG method, for the problems solved and the grids used, by a factor reaching a value as high as 15. Copyright © 2003 John Wiley & Sons, Ltd.     

A note on compressive limiting for two‐material flows

In this short note we describe a simple extension to the multi‐material shock‐capturing algorithm presented in (J. Comput. Phys. 2007; 223 :262–297) that can be used to maintain sharp material interfaces. The method takes the form of an artificial compression which is designed so that the material indicator jumps across only a few cells but which does not excite physical instabilities in the flow. The advantages of the approach include its simplicity and flexibility in that it provides a parameter that effectively determines the captured interface thickness. Copyright © 2009 John Wiley & Sons, Ltd.     

A boundary conforming structured grid for global ocean circulation studies

A boundary conforming two-dimensional structured grid for the irregular domain of the world's ocean is generated numerically using differential equation techniques. It is calculated using block structured methods which allow the inclusion of all major bodies of water including seas and basins, and which preserve slope continuity of the co-ordinate lines across the global domain. The block structure is coupled with an innovative blown-up cube model of the Earth which permits all areas of the global ocean to be modeled with the same resolution, eliminating problems associated with polar singularities. The grid is generated on the curved surface of the Earth (rather than the longitude–latitude plane) by employing the Beltrami operator instead of the standard Laplacian operator. Application of the grid to a steady state heat conduction problem shows the relative computational accuracy and the potential to resolve the complex, smaller scale oceanographic phenomena of great importance to global circulation studies. © 1998 John Wiley & Sons, Ltd.     

A novel non‐upwind,interconnected, multi‐grid,overlapping numerical procedure for problems involving fluid flow

A novel numerical procedure for heat, mass and momentum transfer in fluid flow is presented. The new scheme is passed on a non‐upwind, interconnected, multi‐grid, overlapping (NIMO) finite‐difference algorithm. In 2D flows, the NIMO algorithm solves finite‐difference equations for each dependent variable on four overlapping grids. The finite‐difference equations are formulated using the control‐volume approach, such that no interpolations are needed for computing the convective fluxes. For a particular dependent variable, four fields of values are produced. The NIMO numerical procedure is tested against the exact solution of two test problems. The first test problem is an oblique laminar 2D flow with a double step abrupt change in a passive scalar variable for infinite Peclet number. The second test problem is a rotating radial flow in an annular sector with a single step abrupt change in a passive scalar variable for infinite Peclet number. The NIMO scheme produced essentially the exact solution using different uniform and non‐uniform square and rectangular grids for 45 and 30° angle of inclination. All other schemes were unable to capture the exact solution, especially for the rectangular and non‐uniform grids. The NIMO scheme was also successful in predicting the exact solution for the rotating radial flow, using a uniform cylindrical‐polar coordinate grid. Copyright © 2008 John Wiley & Sons, Ltd.