Immersed Interface Finite Element Methods for Elasticity Interface Problems with Non-Homogeneous Jump Conditions

<正>In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.     

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.     

The Immersed Interface Method for a Nonlinear Chemical Diffusion Equation with Local Sites of Reactions

A diffusion equation with nonlinear localized chemical reactions is considered in this paper. As a result of the reactions, although the equation is parabolic, the derivatives of the solution are discontinuous across the interfaces (local sites of reactions). A second-order accurate immersed interface method is constructed for the diffusion equation involving interfaces. The new method is more accurate than the standard approach and it does not require the interfaces to be grid points. Several experiments that confirm second-order accuracy are presented. The efficiency of the proposed algorithm is also demonstrated for solving blow up problems. The proposed technique could be extended for construction of efficient numerical algorithms on uniform grids for the present equations with moving interfaces [9] but more analysis is required.     

Augmented Lagrangian Methods for $p$-Harmonic Flows with the Generalized Penalization Terms and Application to Image Processing

In this paper, we propose a generalized penalization technique and a convex constraint minimization approach for the $p$-harmonic flow problem following the ideas in [Kang & March, IEEE T. Image Process., 16 (2007), 2251-2261]. We use fast algorithms to solve the subproblems, such as the dual projection methods, primal-dual methods and augmented Lagrangian methods. With a special penalization term, some special algorithms are presented. Numerical experiments are given to demonstrate the performance of the proposed methods. We successfully show that our algorithms are effective and efficient due to two reasons: the solver for subproblem is fast in essence and there is no need to solve the subproblem accurately (even 2 inner iterations of the subproblem are enough). It is also observed that better PSNR values are produced using the new algorithms.     

Second Order Multigrid Methods for Elliptic Problems with Discontinuous Coefficients on an Arbitrary Interface,I: One Dimensional Problems

In this paper we present a one dimensional second order accurate method to solve Elliptic equations with discontinuous coefficients on an arbitrary interface. Second order accuracy for the first derivative is obtained as well. The method is based on the Ghost Fluid Method, making use of ghost points on which the value is defined by suitable interface conditions. The multi-domain formulation is adopted, where the problem is split in two sub-problems and interface conditions will be enforced to close the problem. Interface conditions are relaxed together with the internal equations (following the approach proposed in [10] in the case of smooth coefficients), leading to an iterative method on all the set of grid values (inside points and ghost points). A multigrid approach with a suitable definition of the restriction operator is provided. The restriction of the defect is performed separately for both sub-problems, providing a convergence factor close to the one measured in the case of smooth coefficient and independent on the magnitude of the jump in the coefficient. Numerical tests will confirm the second order accuracy. Although the method is proposed in one dimension, the extension in higher dimension is currently underway [12] and it will be carried out by combining the discretization of [10] with the multigrid approach of [11] for Elliptic problems with non-eliminated boundary conditions in arbitrary domain.     

An $L^{\infty}$ Second Order Cartesian Method for 3D Anisotropic Interface Problems

A second order accurate method in the infinity norm is proposed for general three dimensional anisotropic elliptic interface problems in which the solution and its derivatives, the coefficients, and source terms all can have finite jumps across one or several arbitrary smooth interfaces. The method is based on the 2D finite element-finite difference (FE-FD) method but with substantial differences in method derivation, implementation, and convergence analysis. One of challenges is to derive 3D interface relations since there is no invariance anymore under coordinate system transforms for the partial differential equations and the jump conditions. A finite element discretization whose coefficient matrix is a symmetric semi-positive definite is used away from the interface; and the maximum preserving finite difference discretization whose coefficient matrix part is an M-matrix is constructed at irregular elements where the interface cuts through. We aim to get a sharp interface method that can have second order accuracy in the point-wise norm. We show the convergence analysis by splitting errors into several parts. Nontrivial numerical examples are presented to confirm the convergence analysis.     

A Finite Volume Method Preserving Maximum Principle for the Conjugate Heat Transfer Problems with General Interface Conditions

In this paper, we present a unified finite volume method preserving discrete maximum principle (DMP) for the conjugate heat transfer problems with general interface conditions. We prove the existence of the numerical solution and the DMP-preserving property. Numerical experiments show that the nonlinear iteration numbers of the scheme in [24] increase rapidly when the interfacial coefficients decrease to zero. In contrast, the nonlinear iteration numbers of the unified scheme do not increase when the interfacial coefficients decrease to zero, which reveals that the unified scheme is more robust than the scheme in [24]. The accuracy and DMP-preserving property of the scheme are also veried in the numerical experiments.     

The Improved Modified Decomposition Method for the Treatment of Systems of Singular Nonlinear Partial Differential Equations with Initial Conditions

New developed inverse differential operators incorporated into the semi- analytical treatment of the modified decomposition method （MDM） are used to solve the systems of first and second-order singular nonlinear partial differential equations （PDEs） with initial conditions arising in physics. The new proposed method is called the improved modified decomposition method （IMDM）, and is used to the treatment of a few case study initial-value problems. The results obtained by the IMDM are in full agreement with the existing exact analytical solutions.     

Analysis of a Fully Discrete Finite Element Method for Parabolic Interface Problems with Nonsmooth Initial Data

This article concerns numerical approximation of a parabolic interface problem with general $L^2$ initial value. The problem is discretized by a finite element method with a quasi-uniform triangulation of the domain fitting the interface, with piecewise linear approximation to the interface. The semi-discrete finite element problem is furthermore discretized in time by the $k$-step backward difference formula with $ k=1,\ldots,6 $. To maintain high-order convergence in time for possibly nonsmooth $L^2$ initial value, we modify the standard backward difference formula at the first $k-1$ time levels by using a method recently developed for fractional evolution equations. An error bound of $\mathcal{O}(t_n^{-k}\tau^k+t_n^{-1}h^2|\log h|)$ is established for the fully discrete finite element method for general $L^2$ initial data.     

解非线性有限元方程的逐层校正迭代法

本文提出一种求解非线性有限元方程的逐层校正迭代法.有关数值分析表明,当网格分划较细,网格分划参数h_j较小时,仅需一次简单的迭代和校正步骤就可满足数值计算的要求,使用该方法的计算复杂性是最佳阶的,即为O(N_j),其中N_j为最细网格层上离散结点变量的数目.