Direct Forcing Immersed Boundary Methods: Finite Element Versus Finite Volume Approach

Two immersed boundary methods (IBM) for the simulation of conjugate heat transfer problems with complex geometries are introduced: a finite element (IFEM) and a finite volume (IFVM) immersed boundary methods are discussed. In the IFEM a projection approach is presented for the coupled system of time-dependent incompressible Navier-Stokes equations (NSEs) and energy equation in conjunction with the immersed boundary method for solving fluid flow and heat transfer problems in the presence of rigid objects not represented by the underlying mesh. The IBM allows solving the flow for geometries with complex objects without the need of generating a body-fitted mesh. Dirichlet boundary constraints are satisfied applying a boundary force at the immersed body surface. Using projection and interpolation operators from the fluid volume mesh to the solid surface mesh (i.e., the “immersed” boundary) and vice versa, it is possible to impose the extra constraint to the NSEs as a Lagrange multiplier in a fashion very similar to the effect pressure has on the momentum equations to satisfy the divergence-free constraint. The IFEM approach presented shows third order accuracy in space and second order accuracy in time when the simulation results for the Taylor-Green decaying vortex are compared to the analytical solution. For the IFVM a ghost-cell approach with sharp interface scheme is used to enforce the boundary condition at the fluid/solid interface. The interpolation procedure at the immersed boundary preserves the overall second order accuracy of the base solver. The developed ghost-cell method is applied on a staggered configuration with the Semi-Implicit Method for Pressure-Linked Equations Revised algorithm. Second order accuracy in space and first order accuracy in time are obtained when the Taylor-Green decaying vortex test case is compared to the IFVM analytical solution. Computations were performed using the IFEM and IFVM approaches for the two-dimensional flow over a backward-facing step, two-dimensional flow past a stationary circular cylinder, three-dimensional flow past a sphere and two and three-dimensional natural convection in an enclosure with/without immersed body. The numerical results obtained with the discussed IFEM and IFVM were compared against other IBMs available in literature and simulations performed with the commercial computational fluid dynamics code STAR-CCM+/V7.04.006. The benchmark test cases showed that the numerical results obtained with the implemented immersed boundary methods are in good agreement with the predictions from STAR-CCM+ and the numerical data from the other IBMs. The immersed boundary method based of finite element approach is numerically more accurate than the IBM based on finite volume discretization. In contrast, the latter is computationally more efficient than the former

Numerical simulation of the tip vortex off a low-aspect-ratio wing at transonic speed

The viscous transonic flow around a low aspect ratio wing was computed by an implicit, three dimensional, thin-layer Navier-Stokes solver. The grid around the geometry of interest is obtained numerically as a solution to a Dirichlet problem for the cube. A low aspect ratio wing with large sweep, twist, taper, and camber is the chosen geometry. The topology chosen to wrap the mesh around the wing with good tip resolution is a C-O type mesh. The flow around the wing was computed for a free stream Mach number of 0.82 at an angle of attack of 5 deg. At this Mach number, an oblique shock forms on the upper surface of the wing, and a tip vortex and three dimensional flow separation off the wind surface are observed. Particle path lines indicate that the three dimensional flow separation on the wing surface is part of the roots of the tip vortex formation. The lifting of the tip vortex before the wing trailing edge is observed by following the trajectory of particles release around the wing tip

Use of a hyperbolic grid generation scheme in simulating supersonic viscous flow about three-dimensional winged configuration

The present paper describes a numerical mesh generation technique to be used with an implicit finite difference method for simulating visous supersonic flow about low-aspect-ratio wing body configurations using a single grid strategy. The computational domain is segmented into multiple regions, with borders located in supersonic areas to avoid the otherwise costly interfacing procedure between adjacent segments. The numerical procedure is applied to calculate the turbulent flow around the shuttle orbiter and a canard projectile at supersonic free stream Mach number

A fast immersed boundary method for external incompressible viscous flows using lattice Green's functions

A new parallel, computationally efficient immersed boundary method for solving three-dimensional, viscous, incompressible flows on unbounded domains is presented. Immersed surfaces with prescribed motions are generated using the interpolation and regularization operators obtained from the discrete delta function approach of the original (Peskin's) immersed boundary method. Unlike Peskin's method, boundary forces are regarded as Lagrange multipliers that are used to satisfy the no-slip condition. The incompressible Navier-Stokes equations are discretized on an unbounded staggered Cartesian grid and are solved in a finite number of operations using lattice Green's function techniques. These techniques are used to automatically enforce the natural free-space boundary conditions and to implement a novel block-wise adaptive grid that significantly reduces the run-time cost of solutions by limiting operations to grid cells in the immediate vicinity and near-wake region of the immersed surface. These techniques also enable the construction of practical discrete viscous integrating factors that are used in combination with specialized half-explicit Runge-Kutta schemes to accurately and efficiently solve the differential algebraic equations describing the discrete momentum equation, incompressibility constraint, and no-slip constraint. Linear systems of equations resulting from the time integration scheme are efficiently solved using an approximation-free nested projection technique. The algebraic properties of the discrete operators are used to reduce projection steps to simple discrete elliptic problems, e.g. discrete Poisson problems, that are compatible with recent parallel fast multipole methods for difference equations. Numerical experiments on low-aspect-ratio flat plates and spheres at Reynolds numbers up to 3,700 are used to verify the accuracy and physical fidelity of the formulation.Comment: 32 pages, 9 figures; preprint submitted to Journal of Computational Physic