Efficient high-order immersed interface methods for heat equations with interfaces

An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.

Superconvergence of Direct Discontinuous Galerkin Methods:Eigen-structure Analysis Based on Fourier Approach

This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.

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.

Formation of Uniform Oil-Soluble Fe_3O_4 Nanoparticles via Oil-Water Interface System

Nowadays,the novel oil water interface method has attracted a considerable attention owing to the advantages of mild reaction conditions,simple operation,low cost,and high efficiency.In this paper,uniform oil-soluble Fe_3O_4 nanoparticles(NPs) were synthesized by oil-water interface method from mixing iron tristearate of 0.067mol/L in cyclohexane with ferrous sulfate in water.The as-prepared products were characterized by X-ray diffraction(XRD),transmission electron microscopy(TEM),vibrating sample magnetometer(VSM),Fourier transform infrared spectroscopy(FT-IR) and thermogravimetric analyzer(TGA).TEM images and XRD profiles showed that the size of the oil-soluble products ranged in 1.7-6.9 nm.VSM indicated that the Fe_3O_4 NPs were superparamagetic.FT-IR and TGA proved that oleic acid was combined to the surface of Fe_3O_4 NPs closely.TEM images and XRD profiles revealed that the most suitable reaction concentration of NH_3·H_2O,oleic acid/water in volume,reaction temperature and reaction time were 4.5 mol/L,50:1 000,80℃ and 6 h,respectively.The formation mechanism of the nearly monodispersed Fe_3O_4 NPs was that the preformed Fe_3O_4 nuclei were capped by oleic acid as early as the nucleation occurred in oil-water interface and subsequently entered into oil phase to stop growing.

The Algebraic Immersed Interface and Boundary Method for Elliptic Equations with Jump Conditions

A new simple fictitious domain method, the algebraic immersed interface and boundary (AIIB) method, is presented for elliptic equations with immersed interface conditions. This method allows jump conditions on immersed interfaces to be discretized with a good accuracy on a compact stencil. Auxiliary unknowns are created at existing grid locations to increase the degrees of freedom of the initial problem. These auxiliary unknowns allow imposing various constraints to the system on interfaces of complex shapes. For instance, the method is able to deal with immersed interfaces for elliptic equations with jump conditions on the solution or discontinuous coefficients with a second order of spatial accuracy. As the AIIB method acts on an algebraic level and only changes the problem matrix, no particular attention to the initial discretization is required. The method can be easily implemented in any structured grid code and can deal with immersed boundary problems too. Several validation problems are presented to demonstrate the interest and accuracy of the method.

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.

A sharp interface approach for cavitation modeling using volume-of-fluid and ghost-fluid methods

This paper describes a novel sharp interface approach for modeling the cavitation phenomena in incompressible viscous flows. A one-field formulation is adopted for the vapor-liquid two-phase flow and the interface is tracked using a volume of fluid（VOF） method. Phase change at the interface is modeled using a simplification of the Rayleigh-Plesset equation. Interface jump conditions in velocity and pressure field are treated using a level set based ghost fluid method. The level set function is constructed from the volume fraction function. A marching cubes method is used to compute the interface area at the interface grid cells. A parallel fast marching method is employed to propagate interface information into the field. A description of the equations and numerical methods is presented. Results for a cavitating hydrofoil are compared with experimental data.

An Efficient Neural-Network and Finite-Difference Hybrid Method for Elliptic Interface Problems with Applications

A new and efficient neural-network and finite-difference hybrid method is developed for solving Poisson equation in a regular domain with jump discontinuities on embedded irregular interfaces.Since the solution has low regularity across the interface,when applying finite difference discretization to this problem,an additional treatment accounting for the jump discontinuities must be employed.Here,we aim to elevate such an extra effort to ease our implementation by machine learning methodology.The key idea is to decompose the solution into singular and regular parts.The neural network learning machinery incorporating the given jump conditions finds the singular solution,while the standard five-point Laplacian discretization is used to obtain the regular solution with associated boundary conditions.Regardless of the interface geometry,these two tasks only require supervised learning for function approximation and a fast direct solver for Poisson equation,making the hybrid method easy to implement and efficient.The two-and three-dimensional numerical results show that the present hybrid method preserves second-order accuracy for the solution and its derivatives,and it is comparable with the traditional immersed interface method in the literature.As an application,we solve the Stokes equations with singular forces to demonstrate the robustness of the present method.

Interface Finite Volume Method with Low Calculating Amount for Step Size Varying Electro-Quasistatic Field Problem

In this paper,a novel numerical method called interface finite volume method(I-FVM)for calculation of step-varying Electro-quasistatic(EQS)field is proposed.First,the principle of I-FVM is derived.Then,with numerical example of double layers parallel plate structure under step voltage which has an analytical solution,effectiveness and correctness of the I-FVM are verified.It can be found the calculating time of I-FVM is only 30%of normal FVM without decreasing accuracy during the whole calculating process.Furthermore,an engineering example about the electric field of DBC(Direct Bonding Copper)structure in a high voltage IGBT device is given.It can be found that accuracy of the I-FVM is the same as normal FVM,while time cost of I-FVM is only 20.8%of normal FVM.At last,the I-FVM is extended to one dimension based on the two-direction tri-diagonal matrix algorithm(TDMA)method given in this paper which can save processing of LU decomposition compared to one-dimensional traditional TDMA.In conclusion,the novel method called I-FVM proposed in this paper can decrease calculating amount for a step size varying electro-quasistatic field calculation problem.It may be a good method for large-scale EQS field calculation.

A Discontinuity and Cusp Capturing PINN for Stokes Interface Problems with Discontinuous Viscosity and Singular Forces

In this paper,we present a discontinuity and cusp capturing physicsinformed neural network(PINN)to solve Stokes equations with a piecewiseconstant viscosity and singular force along an interface.We first reformulate the governing equations in each fluid domain separately and replace the singular force effect with the traction balance equation between solutions in two sides along the interface.Since the pressure is discontinuous and the velocity has discontinuous derivatives across the interface,we hereby use a network consisting of two fully-connected sub-networks that approximate the pressure and velocity,respectively.The two sub-networks share the same primary coordinate input arguments but with different augmented feature inputs.These two augmented inputs provide the interface information,so we assume that a level set function is given and its zero level set indicates the position of the interface.The pressure sub-network uses an indicator function as an augmented input to capture the function discontinuity,while the velocity sub-network uses a cusp-enforced level set function to capture the derivative discontinuities via the traction balance equation.We perform a series of numerical experiments to solve two-and three-dimensional Stokes interface problems and perform an accuracy comparison with the augmented immersed interface methods in literature.Our results indicate that even a shallow network with a moderate number of neurons and sufficient training data points can achieve prediction accuracy comparable to that of immersed interface methods.

A Hybrid Immersed Interface Method for Driven Stokes Flow in an Elastic Tube

We present a hybrid numerical method for simulating fluid flow through a compliant,closed tube,driven by an internal source and sink.Fluid is assumed to be highly viscous with its motion described by Stokes flow.Model geometry is assumed to be axisymmetric,and the governing equations are implemented in axisymmetric cylindrical coordinates,which capture 3D flow dynamics with only 2D computations.We solve the model equations using a hybrid approach:we decompose the pressure and velocity fields into parts due to the surface forcings and due to the source and sink,with each part handled separately by means of an appropriate method.Because the singularly-supported surface forcings yield an unsmooth solution,that part of the solution is computed using the immersed interface method.Jump conditions are derived for the axisymmetric cylindrical coordinates.The velocity due to the source and sink is calculated along the tubular surface using boundary integrals.Numerical results are presented that indicate second-order accuracy of the method.

Phase Field Models Versus Parametric Front Tracking Methods: Are They Accurate and Computationally Efficient?

We critically compare the practicality and accuracy of numerical approximations of phase field models and sharp interface models of solidification.Here we focus on Stefan problems,and their quasi-static variants,with applications to crystal growth.New approaches with a high mesh quality for the parametric approximations of the resulting free boundary problems and new stable discretizations of the anisotropic phase field system are taken into account in a comparison involving benchmark problems based on exact solutions of the free boundary problem.

Multi-domain Spectral Immersed Interface Method for Solving Elliptic Equation with a Global Description of Discontinuous Functions

This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed interlace method（IIM） divides the whole computation domain into the smooth and discontinuous parts. Fewer points on the smooth domains are used via taking advantage of the high accuracy property of the spectral method, but more points on the discontinuous domains are employed to enhance the resolution of the calculation. Two that the domain decomposition technique can placed around the discontinuity. The present reached, in spite of the enlarged computational discontinuous problems are tested to verify the present method. The results show reduce the error of the spectral IIM, especially when more collocation points are method is t：avorable for the reason that the same level of the accuracy can be domain.

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.

A Hybrid Fluid-Solid Interaction Scheme Combining the Multi-Component Diffuse Interface Method and the Material Point Method

We propose a hybrid scheme combing the diffuse interface method and the material point method to simulate the complex interactions between the multiphase compressible flow and elastoplastic solid.The multiphase flow is modelled by the multi-component model and solved using a generalized Godunov method in the Eulerian grids,while the elastoplastic solid is solved by the classical material point method in a combination of Lagrangian particles and Eulerian background grids.In order to facilitate the simulation of fluid-solid interactions,the solid variables are further interpolated to the cell center and coexist with the fluid in the same cell.An instantaneous relaxation procedure of velocity and pressure is adopted to simulate the momentum and energy transfers between various materials,and to keep the system within a tightly coupled interaction.Several numerical examples,including shock tube problem,gasbubble problem,air blast,underwater explosion and high speed impact applications are presented to validate the numerical scheme.

High Order Scheme for Schrödinger Equation with Discontinuous Potential I: Immersed Interface Method

The immersed interface method is modified to compute Schrödinger equation with discontinuous potential.By building the jump conditions of the solution into the finite difference approximation near the interface,this method can give at least second order convergence rate for the numerical solution on uniform cartesian grids.The accuracy of this algorithm is tested via several numerical examples.

An Immersed Interface Method for the Simulation of Inextensible Interfaces in Viscous Fluids

In this paper,an immersed interface method is presented to simulate the dynamics of inextensible interfaces in an incompressible flow.The tension is introduced as an augmented variable to satisfy the constraint of interface inextensibility,and the resulting augmented system is solved by the GMRES method.In this work,the arclength of the interface is locally and globally conserved as the enclosed region undergoes deformation.The forces at the interface are calculated from the configuration of the interface and the computed augmented variable,and then applied to the fluid through the related jump conditions.The governing equations are discretized on a MAC grid via a second-order finite difference scheme which incorporates jump contributions and solved by the conjugate gradient Uzawa-type method.The proposed method is applied to several examples including the deformation of a liquid capsule with inextensible interfaces in a shear flow.Numerical results reveal that both the area enclosed by interface and arclength of interface are conserved well simultaneously.These provide further evidence on the capability of the present method to simulate incompressible flows involving inextensible interfaces.

A Simple 3D Immersed InterfaceMethod for Stokes Flow with Singular Forces on Staggered Grids

In this paper,a fairly simple 3D immersed interface method based on the CG-Uzawa type method and the level set representation of the interface is employed for solving three-dimensional Stokes flow with singular forces along the interface.The method is to apply the Taylor’s expansions only along the normal direction and incorporate the jump conditions up to the second normal derivatives into the finite difference schemes.A second order geometric iteration algorithm is employed for computing orthogonal projections on the surface with third-order accuracy.The Stokes equations are discretized involving the correction terms on staggered grids and then solved by the conjugate gradient Uzawa type method.The major advantages of the present method are the special simplicity,the ability in handling the Dirichlet boundary conditions,and no need of the pressure boundary condition.The method can also preserve the volume conservation and the discrete divergence free condition very well.The numerical results show that the proposed method is second order accurate and efficient.

An Iterative Two-Fluid Pressure Solver Based on the Immersed Interface Method

An iterative solver based on the immersed interface method is proposed to solve the pressure in a two-fluid flow on a Cartesian grid with second-order accuracy in the infinity norm.The iteration is constructed by introducing an unsteady term in the pressure Poisson equation.In each iteration step,a Helmholtz equation is solved on the Cartesian grid using FFT.The combination of the iteration and the immersed interface method enables the solver to handle various jump conditions across twofluid interfaces.This solver can also be used to solve Poisson equations on irregular domains.

Immersed Interface CIP for One Dimensional Hyperbolic Equations

The immersed interface technique is incorporated into CIP method to solve one-dimensional hyperbolic equations with piecewise constant coefficients.The proposed method achieves the third order of accuracy in time and space in the vicinity of the interface where the coefficients have jump discontinuities,which is the same order of accuracy of the standard CIP scheme.Some numerical tests are given to verify the accuracy of the proposed method.