A Fast Cartesian Grid-Based Integral Equation Method for Unbounded Interface Problems with Non-Homogeneous Source Terms

This work presents a fast Cartesian grid-based integral equation method for unbounded interface problems with non-homogeneous source terms.The unbounded interface problem is solved with boundary integral equation methods such that infinite boundary conditions are satisfied naturally.This work overcomes two difficulties.The first difficulty is the evaluation of singular integrals.Boundary and volume integrals are transformed into equivalent but much simpler bounded interface problems on rectangular domains,which are solved with FFT-based finite difference solvers.The second one is the expensive computational cost for volume integrals.Despite the use of efficient interface problem solvers,the evaluation for volume integrals is still expensive due to the evaluation of boundary conditions for the simple interface problem.The problem is alleviated by introducing an auxiliary circle as a bridge to indirectly evaluate boundary conditions.Since solving boundary integral equations on a circular boundary is so accurate,one only needs to select a fixed number of points for the discretization of the circle to reduce the computational cost.Numerical examples are presented to demonstrate the efficiency and the second-order accuracy of the proposed numerical method.

Development of Cartesian grid method for simulation of violent ship-wave interactions

We present a Cartesian grid method for numerical simulation of strongly nonlinear phenomena of ship-wave interactions. The Constraint Interpolation Profile （CIP） method is applied to the flow solver, which can efficiently increase the discretization accuracy on the moving boundaries for the Cartesian grid method. Tangent of Hyperbola for Interface Capturing （THINC） is imple- mented as an interface capturing scheme for free surface calculation. An improved immersed boundary method is developed to treat moving bodies with complex-shaped geometries. In this paper, the main features and some recent improvements of the Cartesian grid method are described and several numerical simulation results are presented to discuss its performance.

A fast and automatic full-potential finite volume solver on Cartesian grids for unconventional configurations

To meet the requirements of fast and automatic computation of subsonic and transonic aerodynamics in aircraft conceptual design,a novel finite volume solver for full potential flows on adaptive Cartesian grids is developed in this paper.Cartesian grids with geometric adaptation are firstly generated automatically with boundary cells processed by cell-cutting and cell-merging algorithms.The nonlinear full potential equation is discretized by a finite volume scheme on these Cartesian grids and iteratively solved in an implicit fashion with a generalized minimum residual（GMRES） algorithm.During computation,solution-based mesh adaptation is also applied so as to capture flow features more accurately.An improved ghost-cell method is proposed to implement the non-penetration wall boundary condition where the velocity-potential of a ghost cell is modified by an analytic method instead.According to the characteristics of the Cartesian grids,the Kutta condition is applied by specially computing the gradients on Kutta-faces without directly assigning the potential jump to cells adjacent wake faces,which can significantly improve the solution converging speed.The feasibility and accuracy of the proposed method are validated by several typical cases of sub/transonic flows around an ONERA M6 wing,a DLR-F4 wing-body,and an unconventional figuration of a blended wing body（BWB）.The validation cases demonstrate a fast convergence with fully automatic grid treatment and computation,and the results suggest its capacity in application for aircraft conceptual design.

Sharp interface direct forcing immersed boundary methods： A summary of some algorithms and applications

Body-fitted mesh generation has long been the bottleneck of simulating fluid flows involving complex geometries. Immersed boundary methods are non-boundary-conforming methods that have gained great popularity in the last two decades for their simplicity and flexibility, as well as their non-compromised accuracy. This paper presents a summary of some numerical algori- thms along the line of sharp interface direct forcing approaches and their applications in some practical problems. The algorithms include basic Navier-Stokes solvers, immersed boundary setup procedures, treatments of stationary and moving immersed bounda- ries, and fluid-structure coupling schemes. Applications of these algorithms in particulate flows, flow-induced vibrations, biofluid dynamics, and free-surface hydrodynamics are demonstrated. Some concluding remarks are made, including several future research directions that can further expand the application regime of immersed boundary methods.

An Adaptive Mesh Refinement Strategy for Immersed Boundary/Interface Methods

An adaptive mesh refinement strategy is proposed in this paper for the Immersed Boundary and Immersed Interface methods for two-dimensional elliptic interface problems involving singular sources.The interface is represented by the zero level set of a Lipschitz functionϕ(x,y).Our adaptive mesh refinement is done within a small tube of|ϕ(x,y)|δwith finer Cartesian meshes.The discrete linear system of equations is solved by a multigrid solver.The AMR methods could obtain solutions with accuracy that is similar to those on a uniform fine grid by distributing the mesh more economically,therefore,reduce the size of the linear system of the equations.Numerical examples presented show the efficiency of the grid refinement strategy.

A Fourth-Order Kernel-Free Boundary Integral Method for Interface Problems

This paper presents a fourth-order Cartesian grid based boundary integral method(BIM)for heterogeneous interface problems in two and three dimensional space,where the problem interfaces are irregular and can be explicitly given by parametric curves or implicitly defined by level set functions.The method reformulates the governing equation with interface conditions into boundary integral equations(BIEs)and reinterprets the involved integrals as solutions to some simple interface problems in an extended regular region.Solution of the simple equivalent interface problems for integral evaluation relies on a fourth-order finite difference method with an FFT-based fast elliptic solver.The structure of the coefficient matrix is preserved even with the existence of the interface.In the whole calculation process,analytical expressions of Green’s functions are never determined,formulated or computed.This is the novelty of the proposed kernel-free boundary integral(KFBI)method.Numerical experiments in both two and three dimensions are shown to demonstrate the algorithm efficiency and solution accuracy even for problems with a large diffusion coefficient ratio.

The Immersed Interface Method for Simulating Two-Fluid Flows

We develop the immersed interface method(IIM)to simulate a two-fluid flow of two immiscible fluids with different density and viscosity.Due to the surface tension and the discontinuous fluid properties,the two-fluid flow has nonsmooth velocity and discontinuous pressure across the moving sharp interface separating the two fluids.The IIM computes the flow on a fixed Cartesian grid by incorporating into numerical schemes the necessary jump conditions induced by the interface.We present how to compute these necessary jump conditions from the analytical principal jump conditions derived in[Xu,DCDS,Supplement 2009,pp.838-845].We test our method on some canonical two-fluid flows.The results demonstrate that the method can handle large density and viscosity ratios,is second-order accurate in the infinity norm,and conserves mass inside a closed interface.