A three‐dimensional Cartesian cut cell method for incompressible viscous flow with irregular domains

A three‐dimensional Cartesion cut cell method is presented for the simulations of incompressible viscous flows with irregular domains. A new model (referred to as ‘6+N’ model) is proposed to describe arbitrarily shaped cut cells and treat all the cells as polyhedrons with 6+N faces. The finite volume discretization of the Navier–Stokes equation is then implemented by using the ‘6+N’ model to separate the surface flux integrals into two parts, that is, the fluxes through the basic face of the hexahedron and those through the cutting surfaces. The previously proposed Kitta Cube algorithm and volume computer‐aided design platform (J. Comput. Aided. Des. 2005; 37(4): 1509–1520. Doi:10.1016/j.cad.2005.03.006) are adopted to generate cut cells and provide shape data and physical attributes for the numerical analysis. A modified SIMPLE‐based smoothing pressure correction scheme is applied to suppress checkerboard pressure oscillations caused by the collocated arrangement of velocities and pressure. The calculation accuracy of the numerical method expressed by L₁ and L _∞ norm errors is first demonstrated by the simulation of a pipe flow. Then its feasibility, efficiency, and potential in engineering applications are verified by applying it to solve natural convections between concentric spheres and between eccentric spheres. The heat transfer patterns in eccentric spheres are also obtained by using the numerical method. Copyright © 2011 John Wiley & Sons, Ltd.     

The development of a Cartesian cut cell method for incompressible viscous flows

This paper describes the extension of the Cartesian cut cell method to applications involving unsteady incompressible viscous fluid flow. The underlying scheme is based on the solution of the full Navier–Stokes equations for a variable density fluid system using the artificial compressibility technique together with a Jameson‐type dual time iteration. The computational domain encompasses two fluid regions and the interface between them is treated as a contact discontinuity in the density field, thereby eliminating the need for special free surface tracking procedures. The Cartesian cut cell technique is used for fitting the complex geometry of solid boundaries across a stationary background Cartesian grid which is located inside the computational domain. A time accurate solution is achieved by using an implicit dual‐time iteration technique based on a slope‐limited, high‐order, Godunov‐type scheme for the inviscid fluxes, while the viscous fluxes are estimated using central differencing. Validation of the new technique is by modelling the unsteady Couette flow and the Rayleigh–Taylor instability problems. Finally, a test case for wave run‐up and overtopping over an impermeable sea dike is performed. Copyright © 2006 John Wiley & Sons, Ltd.     

A Cartesian cut cell free surface capturing method for 3D water impact problems

A Cartesian cut cell mesh generation procedure is developed together with a finite volume Euler solver for a two‐fluid system with a free surface. A fast and robust triangle to triangle overlap scheme is used to determine the intersection of a body‐surface with the background Cartesian mesh. Improvements to the cut cell routines include a new treatment for multiple cuts within a single cell and a surface trimming procedure to ensure a good quality mesh around solid boundaries. The formulae for calculating all necessary information about a cut cell are also presented. These are generic and can be used for arbitrarily irregular boundary elements. A collocated finite volume method with a high resolution Godunov‐type scheme in space is used for discretization of the governing flow equations. By computing in both the air and water regions simultaneously in a consistent manner, the free surface is automatically captured as a contact discontinuity in the density field without the need for any special free surface tracking method. The algorithm incorporates the artificial compressibility method with a dual time stepping strategy to maintain a divergence free velocity field. The mathematical formulation including its numerical implementation of the method is reviewed and results for a number of test cases are also presented. Copyright © 2012 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.     

A simplified adaptive Cartesian grid system for solving the 2D shallow water equations

This paper presents a new simplified grid system that provides local refinement and dynamic adaptation for solving the 2D shallow water equations (SWEs). Local refinement is realized by simply specifying different subdivision levels to the cells on a background uniform coarse grid that covers the computational domain. On such a non‐uniform grid, the structured property of a regular Cartesian mesh is maintained and neighbor information is determined by simple algebraic relationships, i.e. data structure becomes unnecessary. Dynamic grid adaptation is achieved by changing the subdivision level of a background cell. Therefore, grid generation and adaptation is greatly simplified and straightforward to implement. The new adaptive grid‐based SWE solver is tested by applying it to simulate three idealized test cases and promising results are obtained. The new grid system offers a simplified alternative to the existing approaches for providing adaptive mesh refinement in computational fluid dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

An immersed boundary solver for inviscid compressible flows

In this paper, a simple and efficient immersed boundary (IB) method is developed for the numerical simulation of inviscid compressible Euler equations. We propose a method based on coordinate transformation to calculate the unknowns of ghost points. In the present study, the body‐grid intercept points are used to build a complete bilinear (2‐D)/trilinear (3‐D) interpolation. A third‐order weighted essentially nonoscillation scheme with a new reference smoothness indicator is proposed to improve the accuracy at the extrema and discontinuity region. The dynamic blocked structured adaptive mesh is used to enhance the computational efficiency. The parallel computation with loading balance is applied to save the computational cost for 3‐D problems. Numerical tests show that the present method has second‐order overall spatial accuracy. The double Mach reflection test indicates that the present IB method gives almost identical solution as that of the boundary‐fitted method. The accuracy of the solver is further validated by subsonic and transonic flow past NACA2012 airfoil. Finally, the present IB method with adaptive mesh is validated by simulation of transonic flow past 3‐D ONERA M6 Wing. Global agreement with experimental and other numerical results are obtained.     

A symmetry-preserving Cartesian grid method for computing a viscous flow past a circular cylinder

This article deals with a numerical method for solving the unsteady, incompressible Navier–Stokes equations in domains with arbitrarily-shaped boundaries, where the boundary is represented using the Cartesian grid approach. We introduce a novel cut-cell discretization which preserves the spectral properties of convection and diffusion. Here, convection is discretized by a skew-symmetric operator and diffusion is approximated by a symmetric, positive-definite coefficient matrix. Such a symmetry-preserving discretization conserves the kinetic energy (if the dissipation is turned off) and is stable on any grid. The method is successfully tested for an incompressible, unsteady flow around a circular cylinder at $\mathit{Re}=100$. To cite this article: R. Verstappen, M. Dröge, C. R. Mecanique 333 (2005).     

A three‐dimensional hydrodynamic model for shallow waters using unstructured Cartesian grids

This paper presents a three‐dimensional unstructured Cartesian grid model for simulating shallow water hydrodynamics in lakes, rivers, estuaries, and coastal waters. It is a flux‐based finite difference model that uses a cut‐cell approach to fit the bottom topography and shorelines and, at the same time, has the flexibility of discretizing complex geometries with Cartesian grids that can be arbitrarily downsized in the two horizontal directions simultaneously. Because of the use of Cartesian grids, the grid generation is very simple and does not suffer the grid generation headache often seen in many other unstructured models, as the unstructured Cartesian grid model does not have any requirements on the orthogonality of the grids. The newly developed unstructured Cartesian grid model was validated against analytical solutions for a three‐dimensional seiching case in a rectangular basin, before it was compared with another three‐dimensional model named LESS3D for circulations and salinity transport processes in an idealized embayment that is driven by tides and freshwater inflows. Model tests show that the numerical procedure used in the unstructured Cartesian grid model is robust. Similar to other unstructured models, a variable grid size has resulted in a smaller number of grids required for a reasonable model simulation, which in turn reduces the CPU time used in the model run. Copyright © 2010 John Wiley & Sons, Ltd.     

An efficient second‐order accurate cut‐cell method for solving the variable coefficient Poisson equation with jump conditions on irregular domains

A numerical method is presented for solving the variable coefficient Poisson equation on a two‐dimensional domain in the presence of irregular interfaces across which both the variable coefficients and the solution itself may be discontinuous. The approach involves using piecewise cubic splines to represent the irregular interface, and applying this representation to calculate the volume and area of each cut cell. The fluxes across the cut‐cell faces and the interface faces are evaluated using a second‐order accurate scheme. The deferred correction approach is used, resulting in a computational stencil for the discretized Poisson equation on an irregular (complex) domain that is identical to that obtained on a regular (simple) domain. In consequence, a highly efficient multigrid solver based on the additive correction multigrid (ACM) method can be applied to solve the current discretized equation system. Several test cases (for which exact solutions to the variable coefficient Poisson equation with and without jump conditions are known) have been used to evaluate the new methodology for discretization on an irregular domain. The numerical solutions show that the new algorithm is second‐order accurate as claimed, even in the presence of jump conditions across an interface. Copyright © 2006 John Wiley & Sons, Ltd.     

An improved ghost‐cell immersed boundary method for compressible flow simulations

This study presents an improved ghost‐cell immersed boundary approach to represent a solid body in compressible flow simulations. In contrast to the commonly used approaches, in the present work, ghost cells are mirrored through the boundary described using a level‐set method to farther image points, incorporating a higher‐order extra/interpolation scheme for the ghost‐cell values. A sensor is introduced to deal with image points near the discontinuities in the flow field. Adaptive mesh refinement is used to improve the representation of the geometry efficiently in the Cartesian grid system. The improved ghost‐cell method is validated against four test cases: (a) double Mach reflections on a ramp, (b) smooth Prandtl–Meyer expansion flows, (c) supersonic flows in a wind tunnel with a forward‐facing step, and (d) supersonic flows over a circular cylinder. It is demonstrated that the improved ghost‐cell method can reach the accuracy of second order in L¹ norm and higher than first order in L^∞ norm. Direct comparisons against the cut‐cell method demonstrate that the improved ghost‐cell method is almost equally accurate with better efficiency for boundary representation in high‐fidelity compressible flow simulations. Copyright © 2016 John Wiley & Sons, Ltd.     

A new characteristic approach for incompressible thermo‐flow in Cartesian and non‐Cartesian grids

A virtual‐characteristic approach is developed for thermo‐flow with finite‐volume methodology in which a multidimensional characteristic (MC) scheme is applied along with artificial compressibility. To obtain compatibility equations and pseudo‐characteristics, energy equation is taken into account in the MC scheme. With this inherent upwinding of convective fluxes, no artificial viscosity is required even at high Reynolds numbers. Another remarkable advantage of the MC scheme lies in its faster convergence rate with respect to the averaging scheme that is found to exhibit substantial delays in convergence. As benchmarks, forced and mixed convections in a cavity and in flow over cylinder and between parallel plates are examined for a wide range of Reynolds, Grashof, and Prandtl numbers. The MC and averaging schemes are applied for simulation purposes. Results show the better performance of the MC scheme in forced and mixed convections. Results confirm the robustness of the MC scheme in terms of accuracy and convergence. Copyright © 2015 John Wiley & Sons, Ltd.     

An unfitted interior penalty discontinuous Galerkin method for incompressible Navier–Stokes two‐phase flow

A discontinuous Galerkin method for the solution of the immiscible and incompressible two‐phase flow problem based on the nonsymmetric interior penalty method is presented. Therefore, the incompressible Navier–Stokes equation is solved for a domain decomposed into two subdomains with different values of viscosity and density as well as a singular surface tension force. On the basis of a piecewise linear approximation of the interface, meshes for both phases are cut out of a structured mesh. The discontinuous finite elements are defined on the resulting Cartesian cut‐cell mesh and may therefore approximate the discontinuities of the pressure and the velocity derivatives across the interface with high accuracy. As the mesh resolves the interface, regularization of the density and viscosity jumps across the interface is not required. This preserves the local conservation property of the velocity field even in the vicinity of the interface and constitutes a significant advantage compared with standard methods that require regularization of these discontinuities and cannot represent the jumps and kinks in pressure and velocity. A powerful subtessellation algorithm is incorporated to allow the usage of standard time integrators (such as Crank–Nicholson) on the time‐dependent mesh. The presented discretization is applicable to both the two‐dimensional and three‐dimensional cases. The performance of our approach is demonstrated by application to a two‐dimensional benchmark problem, allowing for a thorough comparison with other numerical methods. Copyright © 2012 John Wiley & Sons, Ltd.     

A cell‐based smoothed finite element method with semi‐implicit CBS procedures for incompressible laminar viscous flows

In this paper, the cell‐based smoothed finite element method (CS‐FEM) with the semi‐implicit characteristic‐based split (CBS) scheme (CBS/CS‐FEM) is proposed for computational fluid dynamics. The 3‐node triangular (T3) element and 4‐node quadrilateral (Q4) element are used for present CBS/CS‐FEM for two‐dimensional flows. The 8‐node hexahedral element (H8) is used for three‐dimensional flows. Two types of CS‐FEM are implemented in this paper. One is standard CS‐FEM with quadrilateral gradient smoothing cells for Q4 element and hexahedron cells for H8 element. Another is called as n‐sided CS‐FEM (nCS‐FEM) whose gradient smoothing cells are triangles for Q4 element and pyramids for H8 element. To verify the proposed methods, benchmarking problems are tested for two‐dimensional and three‐dimensional flows. The benchmarks show that CBS/CS‐FEM and CBS/nCS‐FEM are capable to solve incompressible laminar flow and can produce reliable results for both steady and unsteady flows. The proposed CBS/CS‐FEM method has merits on better robustness against distorted mesh with only slight more computation time and without losing accuracy, which is important for problems with heavy mesh distortion. The blood flow in carotid bifurcation is also simulated to show capabilities of proposed methods for realistic and complicated flow problems.     

A finite volume‐based high‐order,Cartesian cut‐cell method for wave propagation

Computational aeroacoustics requires numerical techniques capable of yielding low artificial dispersion and dissipation to preserve the amplitude and the frequency characteristics of the physical processes. Furthermore, for engineering applications, the techniques need to handle irregular geometries associated with realistic configurations. We address these issues by developing an optimized prefactored compact finite volume (OPC‐fv) scheme along with a Cartesian cut‐cell technique. The OPC‐fv scheme seeks to minimize numerical dispersion and dissipation while satisfying the conservation laws. The cut‐cell approach treats irregularly shaped boundaries using divide‐and‐merge procedures for the Cartesian cells while maintaining a desirable level of accuracy. We assess these techniques using several canonical test problems, involving different levels of physical and geometric complexities. Richardson extrapolation is an effective tool for evaluating solutions of no high gradients or discontinuities, and is used to evaluate the performance of the solution technique. It is demonstrated that while the cut‐cell method has a modest effect on the order of accuracy, it is a robust method. The combined OPC‐fv scheme and the Cartesian cut‐cell technique offer good accuracy as well as geometric flexibility. Copyright © 2008 John Wiley & Sons, Ltd.     

Vortex‐in‐cell method combined with a boundary element method for incompressible viscous flow analysis

In this study, an immersed boundary vortex‐in‐cell (VIC) method for simulating the incompressible flow external to two‐dimensional and three‐dimensional bodies is presented. The vorticity transport equation, which is the governing equation of the VIC method, is represented in a Lagrangian form and solved by the vortex blob representation of the flow field. In the present scheme, the treatment of convection and diffusion is based on the classical fractional step algorithm. The rotational component of the velocity is obtained by solving Poisson's equation using an FFT method on a regular Cartesian grid, and the solenoidal component is determined from solving an integral equation using the panel method for the convection term, and the diffusion term is implemented by a particle strength exchange scheme. Both the no‐slip and no‐through flow conditions associated with the surface boundary condition are satisfied by diffusing vortex sheet and distributing singularities on the body, respectively. The present method is distinguished from other methods by the use of the panel method for the enforcement of the no‐through flow condition. The panel method completes making use of the immersed boundary nature inherent in the VIC method and can be also adopted for the calculation of the pressure field. The overall process is parallelized using message passing interface to manage the extensive computational load in the three‐dimensional flow simulations. Copyright © 2011 John Wiley & Sons, Ltd.     

A high‐order upwind control‐volume method based on integrated RBFs for fluid‐flow problems

This paper is concerned with the development of a high‐order upwind conservative discretization method for the simulation of flows of a Newtonian fluid in two dimensions. The fluid‐flow domain is discretized using a Cartesian grid from which non‐overlapping rectangular control volumes are formed. Line integrals arising from the integration of the diffusion and convection terms over control volumes are evaluated using the middle‐point rule. One‐dimensional integrated radial basis function schemes using the multiquadric basis function are employed to represent the variations of the field variables along the grid lines. The convection term is effectively treated using an upwind scheme with the deferred‐correction strategy. Several highly non‐linear test problems governed by the Burgers and the Navier–Stokes equations are simulated, which show that the proposed technique is stable, accurate and converges well. Copyright © 2010 John Wiley & Sons, Ltd.     

A Navier–Stokes solver for complex three-dimensional turbulent flows adopting non-linear modelling of the Reynolds stresses

A non-linear modelling of the Reynolds stresses has been incorporated into a Navier–Stokes solver for complex three-dimensional geometries. A k–ε model, adopting a modelling of the turbulent transport which is not based on the eddy viscosity, has been written in generalised co-ordinates and solved with a finite volume approach, using both a GMRES solver and a direct solver for the solution of the linear systems of equations. An additional term, quadratic in the main strain rate, has been introduced into the modelling of the Reynolds stresses to the basic Boussinesq's form; the corresponding constant has been evaluated through comparison with the experimental data. The computational procedure is implemented for the flow analysis in a 90° square section bend and the obtained results show that with the non-linear modelling a much better agreement with the measured data is obtained, both for the velocity and the pressure. The importance of the convection scheme is also discussed, showing how the effect of the non-linear correction added to the Reynolds stresses is effectively hidden by the additional numerical diffusion introduced by a low-order convection scheme as the first-order upwind scheme, thus making the use of higher order schemes necessary. © 1998 John Wiley & Sons, Ltd.