A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting

An explicit moving boundary method for the numerical solution of time-dependent hyperbolic conservation laws on grids produced by the intersection of complex geometries with a regular Cartesian grid is presented. As it employs directional operator splitting, implementation of the scheme is rather straightforward. Extending the method for static walls from Klein et al., Phil. Trans. Roy. Soc., A367, no. 1907, 4559-4575 (2009), the scheme calculates fluxes needed for a conservative update of the near-wall cut-cells as linear combinations of standard fluxes from a one-dimensional extended stencil. Here the standard fluxes are those obtained without regard to the small sub-cell problem, and the linear combination weights involve detailed information regarding the cut-cell geometry. This linear combination of standard fluxes stabilizes the updates such that the time-step yielding marginal stability for arbitrarily small cut-cells is of the same order as that for regular cells. Moreover, it renders the approach compatible with a wide range of existing numerical flux-approximation methods. The scheme is extended here to time dependent rigid boundaries by reformulating the linear combination weights of the stabilizing flux stencil to account for the time dependence of cut-cell volume and interface area fractions. The two-dimensional tests discussed include advection in a channel oriented at an oblique angle to the Cartesian computational mesh, cylinders with circular and triangular cross-section passing through a stationary shock wave, a piston moving through an open-ended shock tube, and the flow around an oscillating NACA 0012 aerofoil profile.Comment: 30 pages, 27 figures, 3 table

A comparative study of immersed-boundary interpolation methods for a flow around a stationary cylinder at low Reynolds number

The accuracy and computational efficiency of various interpolation methods for the implementation of non grid-confirming boundaries is assessed. The aim of the research is to select an interpolation method that is both efficient and sufficiently accurate to be used in the simulation of vortex induced vibration of the flow around a deformable cylinder. Results are presented of an immersed boundary implementation in which the velocities near nonconfirming boundaries were interpolated in the normal direction to the walls. The flow field is solved on a Cartesian grid using a finite volume method with a staggered variable arrangement. The Strouhal number and Drag coefficient for various cases are reported. The results show a good agreement with the literature. Also, the drag coefficient and Strouhal number results for five different interpolation methods were compared it was shown that for a stationary cylinder at low Reynolds number, the interpolation method could affect the drag coefficient by a maximum 2% and the Strouhal number by maximum of 3%. In addition, the bi-liner interpolation method took about 2% more computational time per vortex shedding cycle in companion to the other methods

Two fluid space-time discontinuous Galerkin finite element method. Part I: numerical algorithm

A novel numerical method for two fluid flow computations is presented, which combines the space-time discontinuous Galerkin finite element discretization with the level set method and cut-cell based interface tracking. The space-time discontinuous Galerkin (STDG) finite element method offers high accuracy, an inherent ability to handle discontinuities and a very local stencil, making it relatively easy to combine with local {\it hp}-refinement. The front tracking is incorporated via cut-cell mesh refinement to ensure a sharp interface between the fluids. To compute the interface dynamics the level set method (LSM) is used because of its ability to deal with merging and breakup. Also, the LSM is easy to extend to higher dimensions. Small cells arising from the cut-cell refinement are merged to improve the stability and performance. The interface conditions are incorporated in the numerical flux at the interface and the STDG discretization ensures that the scheme is conservative as long as the numerical fluxes are conservative

High-order numerical methods for 2D parabolic problems in single and composite domains

In this work, we discuss and compare three methods for the numerical approximation of constant- and variable-coefficient diffusion equations in both single and composite domains with possible discontinuity in the solution/flux at interfaces, considering (i) the Cut Finite Element Method; (ii) the Difference Potentials Method; and (iii) the summation-by-parts Finite Difference Method. First we give a brief introduction for each of the three methods. Next, we propose benchmark problems, and consider numerical tests-with respect to accuracy and convergence-for linear parabolic problems on a single domain, and continue with similar tests for linear parabolic problems on a composite domain (with the interface defined either explicitly or implicitly). Lastly, a comparative discussion of the methods and numerical results will be given.Comment: 45 pages, 12 figures, in revision for Journal of Scientific Computin

A space-time discontinuous Galerkin finite element method for two-fluid problems

A space-time discontinuous Galerkin finite element method for two fluid flow problems is presented. By using a combination of level set and cut-cell methods the interface between two fluids is tracked in space-time. The movement of the interface in space-time is calculated by solving the level set equation, where the interface geometry is identified with the 0-level set. To enhance the accuracy of the interface approximation the level set function is advected with the interface velocity, which for this purpose is extended into the domain. Close to the interface the mesh is locally refined in such a way that the 0-level set coincides with a set of faces in the mesh. The two fluid flow equations are solved on this refined mesh. The procedure is repeated until both the mesh and the flow solution have converged to a reasonable accuracy.\ud The method is tested on linear advection and Euler shock tube problems involving ideal gas and compressible bubbly magma. Oscillations around the interface are eliminated by choosing a suitable interface flux

Grid generation for the solution of partial differential equations

A general survey of grid generators is presented with a concern for understanding why grids are necessary, how they are applied, and how they are generated. After an examination of the need for meshes, the overall applications setting is established with a categorization of the various connectivity patterns. This is split between structured grids and unstructured meshes. Altogether, the categorization establishes the foundation upon which grid generation techniques are developed. The two primary categories are algebraic techniques and partial differential equation techniques. These are each split into basic parts, and accordingly are individually examined in some detail. In the process, the interrelations between the various parts are accented. From the established background in the primary techniques, consideration is shifted to the topic of interactive grid generation and then to adaptive meshes. The setting for adaptivity is established with a suitable means to monitor severe solution behavior. Adaptive grids are considered first and are followed by adaptive triangular meshes. Then the consideration shifts to the temporal coupling between grid generators and PDE-solvers. To conclude, a reflection upon the discussion, herein, is given