A parallel cell‐based DSMC method on unstructured adaptive meshes

A parallel DSMC method based on a cell‐based data structure is developed for the efficient simulation of rarefied gas flows on PC‐clusters. Parallel computation is made by decomposing the computational domain into several subdomains. Dynamic load balancing between processors is achieved based on the number of simulation particles and the number of cells allocated in each subdomain. Adjustment of cell size is also made through mesh adaptation for the improvement of solution accuracy and the efficient usage of meshes. Applications were made for a two‐dimensional supersonic leading‐edge flow, the axi‐symmetric Rothe's nozzle, and the open hollow cylinder flare flow for validation. It was found that the present method is an efficient tool for the simulation of rarefied gas flows on PC‐based parallel machines. Copyright © 2004 John Wiley & Sons, Ltd.     

Compact high‐resolution algorithms for time‐dependent advection on unstructured grids

A technique for constructing monotone, high resolution, multi‐dimensional upwind fluctuation distribution schemes for the scalar advection equation is presented. The method combines the second‐order Lax–Wendroff scheme with the upwind positive streamwise invariant (PSI) scheme via a fluctuation redistribution step, which ensures monotonicity (and which is a generalization of the flux‐corrected transport approach for fluctuation distribution schemes). Furthermore, the concept of a distribution point is introduced, which, when related to the equivalent equation for the scheme, leads to a ‘preferred direction’ for the limiting procedure, and hence to a new distribution of the fluctuation, which retains second‐order accuracy from the Lax–Wendroff scheme, even when the solution contains turning points. Experimental comparisons show that the new method compares favourably in terms of speed, accuracy and robustness with other, similar, techniques. Copyright © 2000 John Wiley & Sons, Ltd.     

An hp‐adaptive unstructured finite volume solver for compressible flows

The idea of hp‐adaptation, which has originally been developed for compact schemes (such as finite element methods), suggests an adaptation scheme using a mixture of mesh refinement and order enrichment based on the smoothness of the solution to obtain an accurate solution efficiently. In this paper, we develop an hp‐adaptation framework for unstructured finite volume methods using residual‐based and adjoint‐based error indicators. For the residual‐based error indicator, we use a higher‐order discrete operator to estimate the truncation error, whereas this estimate is weighted by the solution of the discrete adjoint problem for an output of interest to form the adaptation indicator for adjoint‐based adaptations. We perform our adaptation by local subdivision of cells with nonconforming interfaces allowed and local reconstruction of higher‐order polynomials for solution approximations. We present our results for two‐dimensional compressible flow problems including subsonic inviscid, transonic inviscid, and subsonic laminar flow around the NACA 0012 airfoil and also turbulent flow over a flat plate. Our numerical results suggest the efficiency and accuracy advantages of adjoint‐based hp‐adaptations over uniform refinement and also over residual‐based adaptation for flows with and without singularities.     

A semi‐implicit scheme for 3D free surface flows with high‐order velocity reconstruction on unstructured Voronoi meshes

In this paper, we present a computationally efficient semi‐implicit scheme for the simulation of three‐dimensional hydrostatic free surface flow problems on staggered unstructured Voronoi meshes. For each polygonal control volume, the pressure is defined in the cell center, whereas the discrete velocity field is given by the normal velocity component at the cell faces. A piecewise high‐order polynomial vector velocity field is then reconstructed from the scalar normal velocities at the cell faces by using a new high‐order constrained least‐squares reconstruction operator. The reconstructed high‐order piecewise polynomial velocity field is used for trajectory integration in a semi‐Lagrangian approach to discretize the nonlinear convective terms in the governing PDE. For that purpose, a high‐order Taylor method is used as ODE integrator. The resulting semi‐implicit algorithm is extensively validated on a large set of different academic test problems with exact analytical solution and is finally applied to a real‐world engineering problem consisting of a curved channel upstream of two micro‐turbines of a hydroelectric power plant. For this realistic case, some experimental reference data are available from field measurements. Copyright © 2012 John Wiley & Sons, Ltd.     

A high‐order element based adaptive mesh refinement strategy for three‐dimensional unstructured grid

Adaptive mesh refinement (AMR) shows attractive properties in automatically refining the flow region of interest, and with AMR, better prediction can be obtained with much less labor work and cost compared to manually remeshing or the global mesh refinement. Cartesian AMR is well established; however, AMR on hybrid unstructured mesh, which is heavily used in the high‐Reynolds number flow simulation, is less matured and existing methods may result in degraded mesh quality, which mostly happens in the boundary layer or near the sharp geometric features. User intervention or additional constraints, such as freezing all boundary layer elements or refining the whole boundary layer, are required to assist the refinement process. In this work, a novel AMR strategy is developed to handle existing difficulties. In the new method, high‐order unstructured elements are first generated based on the baseline mesh; then the refinement is conducted in the parametric space; at last, the mesh suitable for the solver is output. Generating refined elements in the parametric space with high‐order elements is the key of this method and this helps to guarantee both the accuracy and robustness. With the current method, 3‐dimensional hybrid unstructured mesh of huge size and complex geometry can be automatically refined, without user intervention nor additional constraints. With test cases including the 2‐dimensional airfoil and 3‐dimensional full aircraft, the current AMR method proves to be accurate, simple, and robust.     

On the implicit time integration of semi‐discrete viscous fluxes on unstructured dynamic meshes

Numerical simulations of viscous flow problems with complex moving and/or deforming boundaries commonly require the solution of the corresponding fluid equations of motion on unstructured dynamic meshes. In this paper, a systematic investigation of the importance of the choice of the mesh configuration for evaluating the viscous fluxes is performed when the semi‐discrete Navier–Stokes equations are time‐integrated using the popular second‐order implicit backward difference algorithm. The findings are illustrated with the simulation of a laminar viscous flow problem around an oscillating airfoil. Copyright © 1999 John Wiley & Sons, Ltd.     

A gas‐kinetic scheme for the simulation of turbulent flows on unstructured meshes

This paper presents a numerical method for simulating turbulent flows via coupling the Boltzmann BGK equation with Spalart–Allmaras one equation turbulence model. Both the Boltzmann BGK equation and the turbulence model equation are carried out using the finite volume method on unstructured meshes, which is different from previous works on structured grid. The application of the gas‐kinetic scheme is extended to the simulation of turbulent flows with arbitrary geometries. The adaptive mesh refinement technique is also adopted to reduce the computational cost and improve the efficiency of meshes. To organize the unstructured mesh data structure efficiently, a non‐manifold hybrid mesh data structure is extended for polygonal cells. Numerical experiments are performed on incompressible flow over a smooth flat plate and compressible turbulent flows around a NACA 0012 airfoil using unstructured hybrid meshes. These numerical results are found to be in good agreement with experimental data and/or other numerical solutions, demonstrating the applicability of the proposed method to simulate both subsonic and transonic turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

A Voronoi‐based ALE solver for the calculation of incompressible flow on deforming unstructured meshes

An Arbitrary Lagrangian–Eulerian method for the calculation of incompressible Navier–Stokes equations in deforming geometries is described. The mesh node connectivity is defined by a Delaunay triangulation of the nodes, whereas the discretized equations are solved using finite volumes defined by the Voronoi dual of the triangulation. For prescribed boundary motion, an automatic node motion algorithm provides smooth motion of the interior nodes. Changes in the connectivity of the nodes are made through the use of local transformations to maintain the mesh as Delaunay. This allows the nodes and their associated Voronoi finite volumes to migrate through the domain in a free manner, without compromising the quality of the mesh. An MAC finite volume solver is applied on the Voronoi dual using a cell‐centred non‐staggered formulation, with cell‐face velocities being calculated by the Rhie–Chow momentum interpolation. Advective fluxes are approximated with the third‐order QUICK differencing scheme. The solver is demonstrated via its application to a driven cavity flow, and the flow about flapping aerofoil geometries. Copyright © 2010 John Wiley & Sons, Ltd.     

A cost‐effective curvature calculation approach for interfacial flows on unstructured meshes

We present a simple and cost‐effective curvature calculation approach for simulations of interfacial flows on structured and unstructured grids. The interface is defined using volume fractions, and the interface curvature is obtained as a function of the gradients of volume fractions. The gradient computation is based on a recently proposed gradient recovery method that mimicks the least squares approach without the need to solve a system of equations and is quite easy to implement on arbitrary polygonal meshes. The resulting interface curvature is used in a continuum surface force formulation within the framework of a well‐balanced finite‐volume algorithm to simulate multiphase flows dominated by surface tension. We show that the proposed curvature calculation is at least as accurate as some of the existing approaches on unstructured meshes while being straightforward to implement on any mesh topology. Numerical investigations also show that spurious currents in stationary problems that are dependent on the curvature calculation methodology are also acceptably low using the proposed approach. Studies on capillary waves and rising bubbles in viscous flows lend credence to the ability of the proposed method as an inexpensive, robust, and reasonably accurate approach for curvature calculation and numerical simulation of multiphase flows.     

A robust and well‐balanced scheme for the 2D Saint‐Venant system on unstructured meshes with friction source term

In the following lines, we propose a numerical scheme for the shallow‐water system supplemented by topography and friction source terms, in a 2D unstructured context. This work proposes an improved version of the well‐balanced and robust numerical model recently introduced by Duran et al. (J. Comp. Phys., 235 , 565–586, 2013) for the pre‐balanced shallow‐water equations, accounting for varying topography. The present work aims at relaxing the robustness condition and includes a friction term. To this purpose, the scheme is modified using a recent method, entirely based on a modified Riemann solver. This approach preserves the robustness and well‐balanced properties of the original scheme and prevents unstable computations in the presence of low water depths. A series of numerical experiments are devoted to highlighting the performances of the resulting scheme. Simulations involving dry areas, complex geometry and topography are proposed to validate the stability of the numerical model in the neighbourhood of wet/dry transitions. Copyright © 2015 John Wiley & Sons, Ltd.     

An eigenvector‐based linear reconstruction scheme for the shallow‐water equations on two‐dimensional unstructured meshes

This paper presents a new approach to MUSCL reconstruction for solving the shallow‐water equations on two‐dimensional unstructured meshes. The approach takes advantage of the particular structure of the shallow‐water equations. Indeed, their hyperbolic nature allows the flow variables to be expressed as a linear combination of the eigenvectors of the system. The particularity of the shallow‐water equations is that the coefficients of this combination only depend upon the water depth. Reconstructing only the water depth with second‐order accuracy and using only a first‐order reconstruction for the flow velocity proves to be as accurate as the classical MUSCL approach. The method also appears to be more robust in cases with very strong depth gradients such as the propagation of a wave on a dry bed. Since only one reconstruction is needed (against three reconstructions in the MUSCL approach) the EVR method is shown to be 1.4–5 times as fast as the classical MUSCL scheme, depending on the computational application. Copyright © 2006 John Wiley & Sons, Ltd.     

Anisotropic mesh adaption for time‐dependent problems

We propose a space–time adaptive procedure for a model parabolic problem based on a theoretically sound anisotropic a posteriori error analysis. A space–time finite element scheme (continuous in space but discontinuous in time) is employed to discretize this problem, thus allowing for non‐matching meshes at different time levels. Copyright © 2008 John Wiley & Sons, Ltd.     

Weighted essential non‐oscillatory schemes for tidal bore on unstructured meshes

In this paper, the third‐order weighted essential non‐oscillatory (WENO) schemes are used to simulate the two‐dimensional shallow water equations with the source terms on unstructured meshes. The balance of the flux and the source terms makes the shallow water equations fit to non‐flat bottom questions. The simulation of a tidal bore on an estuary with trumpet shape and Qiantang river is performed; the results show that the schemes can be used to simulate the current flow accurately and catch the stronger discontinuous in water wave, such as dam break and tidal bore effectively. Copyright © 2008 John Wiley & Sons, Ltd.     

Mesh adaptation for simulation of unsteady flow with moving immersed boundaries

In this work, an approach for performing mesh adaptation in the numerical simulation of two‐dimensional unsteady flow with moving immersed boundaries is presented. In each adaptation period, the mesh is refined in the regions where the solution evolves or the moving bodies pass and is unrefined in the regions where the phenomena or the bodies deviate. The flow field and the fluid–solid interface are recomputed on the adapted mesh. The adaptation indicator is defined according to the magnitude of the vorticity in the flow field. There is no lag between the adapted mesh and the computed solution, and the adaptation frequency can be controlled to reduce the errors due to the solution transferring between the old mesh and the new one. The preservation of conservation property is mandatory in long‐time scale simulations, so a P1‐conservative interpolation is used in the solution transferring. A nonboundary‐conforming method is employed to solve the flow equations. Therefore, the moving‐boundary flows can be simulated on a fixed mesh, and there is no need to update the mesh at each time step to follow the motion or the deformation of the solid boundary. To validate the present mesh adaptation method, we have simulated several unsteady flows over a circular cylinder stationary or with forced oscillation, a single self‐propelled swimming fish, and two fish swimming in the same or different directions. Copyright © 2012 John Wiley & Sons, Ltd.     

Reference values for drag and lift of a two‐dimensional time‐dependent flow around a cylinder

This paper presents a numerical study of a two‐dimensional time‐dependent flow around a cylinder. Its main objective is to provide accurate reference values for the maximal drag and lift coefficient at the cylinder and for the pressure difference between the front and the back of the cylinder at the final time. In addition, the accuracy of these values obtained with different time stepping schemes and different finite element methods is studied. Copyright © 2004 John Wiley & Sons, Ltd.     

A two‐dimensional finite volume morphodynamic model on unstructured triangular grids

We discuss the application of a finite volume method to morphodynamic models on unstructured triangular meshes. The model is based on coupling the shallow water equations for the hydrodynamics with a sediment transport equation for the morphodynamics. The finite volume method is formulated for the quasi‐steady approach and the coupled approach. In the first approach, the steady hydrodynamic state is calculated first and the corresponding water velocity is used in the sediment transport equation to be solved subsequently. The second approach solves the coupled hydrodynamics and sediment transport system within the same time step. The gradient fluxes are discretized using a modified Roe's scheme incorporating the sign of the Jacobian matrix in the morphodynamic system. A well‐balanced discretization is used for the treatment of source terms. We also describe an adaptive procedure in the finite volume method by monitoring the bed–load in the computational domain during its transport process. The method uses unstructured meshes, incorporates upwinded numerical fluxes and slope limiters to provide sharp resolution of steep bed gradients that may form in the approximate solution. Numerical results are shown for a test problem in the evolution of an initially hump‐shaped bed in a squared channel. For the considered morphodynamical regimes, the obtained results point out that the coupled approach performs better than the quasi‐steady approach only when the bed–load rapidly interacts with the hydrodynamics. Copyright © 2009 John Wiley & Sons, Ltd.     

Low‐cost implicit schemes for all‐speed flows on unstructured meshes

Matrix‐free implicit treatments are now commonly used for computing compressible flow problems: a reduced cost per iteration and low‐memory requirements are their most attractive features. This paper explains how it is possible to preserve these features for all‐speed flows, in spite of the use of a low‐Mach preconditioning matrix. The proposed approach exploits a particular property of a widely used low‐Mach preconditioner proposed by Turkel. Its efficiency is demonstrated on some steady and unsteady applications. Copyright © 2008 John Wiley & Sons, Ltd.     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

Anisotropic mesh adaptation: towards user‐independent,mesh‐independent and solver‐independent CFD. Part III. Unstructured meshes

The present paper is the third article in a three‐part series on anisotropic mesh adaptation and its application to two‐ and three‐dimensional, structured and unstructured meshes. This third paper concerns the application of the full adaptation methodology to 2‐D unstructured meshes, including all four mesh modification strategies presented in Part I, i.e. refinement/coarsening, edge swapping and node movement. The mesh adaptation procedure is validated through a careful monitoring of a single adaptation step and of the solution–adaptation loop. Independence from the initial mesh and from the flow solver is illustrated. The efficiency of the overall methodology is investigated on relevant laminar and turbulent flow benchmarks. Copyright © 2002 John Wiley & Sons, Ltd.     

A high‐precision unstructured adaptive mesh technique for gas–liquid two‐phase flows

Adaptive mesh techniques are used widely in the numerical simulations of fluid flows, and the simulation results with high accuracies are obtained by appropriate mesh adaptations. However, gas–liquid two‐phase flows are still difficult to be simulated on adaptive meshes, especially on unstructured adaptive meshes, because the physical phenomena near gas–liquid interfaces are highly complicated and in general, not modeled appropriately on adaptive meshes. In this paper, a high‐precision unstructured adaptive mesh technique for gas–liquid two‐phase flows is developed and verified/validated. In the unstructured adaptive mesh technique, the PLIC algorithm is employed to simulate interfacial dynamic behaviors and, therefore, the reconstruction method for the interfaces in refined cells is developed, which satisfies the gas and liquid volume conservations and geometrical conservations of interfaces. In addition, the physics‐based consideration is performed on the momentum calculations near interfaces, and the calculation method with gas and liquid momentum conservations is developed. For verification, the slotted‐disk revolution problem is solved. As a result, the unstructured adaptive mesh technique succeeds in reproducing the slotted‐disk shape accurately and well maintaining the shape after one full‐revolution. The dam‐break problem is also simulated and the momentum conservative calculation method succeeds in providing physically appropriate results, which show good agreements with experimental data. Therefore, it is confirmed that the developed unstructured adaptive mesh technique is very efficient to simulate gas–liquid two‐phase flows accurately. Copyright © 2010 John Wiley & Sons, Ltd.