In this paper,we consider the plasma-vacuum interface problem in a cylindrical tube region impressed by a special background magnetic field.The interior region is occupied with plasma,which is governed by the incompre...In this paper,we consider the plasma-vacuum interface problem in a cylindrical tube region impressed by a special background magnetic field.The interior region is occupied with plasma,which is governed by the incompressible inviscid and resistive MHD system without damping term.The exterior vacuum region is governed by the so-called the“pre-Maxwell equations”.And on the free interface,additionally,the effect of surface tension is taken into account.The original region can be transformed into a horizontally periodic slab through the cylindrical coordinate transformation,which will be impressed by a uniform nonhorizontal magnetic field.Appending with the appropriate physical boundary conditions,the global well-posedness of the problem is established by the energy method.展开更多
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...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.展开更多
The lowest order Pl-nonconforming triangular finite element method （FEM） for elliptic and parabolic interface problems is investigated. Under some reasonable regularity assumptions on the exact solutions, the optima...The lowest order Pl-nonconforming triangular finite element method （FEM） for elliptic and parabolic interface problems is investigated. Under some reasonable regularity assumptions on the exact solutions, the optimal order error estimates are obtained in the broken energy norm. Finally, some numerical results are provided to verify the theoretical analysis.展开更多
A second order isoparametric finite element method （IPFEM） is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the...A second order isoparametric finite element method （IPFEM） is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the immersed curved interface is discontinuous. Based on an initial Cartesian mesh, a mesh optimization strategy is presented by employing curved boundary elements at the interface, and an incomplete quadratic finite element space is constructed on the optimized mesh. It turns out that the number of curved boundary elements is far less than that of the straight one, and the total degree of freedom is almost the same as the uniform Cartesian mesh. Numerical examples with simple and complicated geometrical interfaces demonstrate the efficiency of the proposed method.展开更多
In this paper,the minimal dissipation local discontinuous Galerkin method is studied to solve the elliptic interface problems in two-dimensional domains.The interface may be arbitrary smooth curves.It is shown that th...In this paper,the minimal dissipation local discontinuous Galerkin method is studied to solve the elliptic interface problems in two-dimensional domains.The interface may be arbitrary smooth curves.It is shown that the error estimates in L;-norm for the solution and the flux are O（h;|log h|）and O(h|log h|;),respectively.In numerical experiments,the successive substitution iterative methods are used to solve the LDG schemes.Numerical results verify the efficiency and accuracy of the method.展开更多
The aim of this paper is to develop a fast multigrid solver for interpolation-free finite volume (FV) discretization of anisotropic elliptic interface problems on general bounded domains that can be described as a uni...The aim of this paper is to develop a fast multigrid solver for interpolation-free finite volume (FV) discretization of anisotropic elliptic interface problems on general bounded domains that can be described as a union of blocks. We assume that the curved interface falls exactly on the boundaries of blocks. The transfinite interpolation technique is applied to generate block-wise distorted quadrilateral meshes, which can resolve the interface with fine geometric details. By an extensive study of the harmonic average point method, an interpolation-free nine-point FV scheme is then derived on such multi-block grids for anisotropic elliptic interface problems with non-homogeneous jump conditions. Moreover, for the resulting linear algebraic systems from cell-centered FV discretization, a high-order prolongation operator based fast cascadic multigrid solver is developed and shown to be robust with respect to both the problem size and the jump of the diffusion coefficients. Various non-trivial examples including four interface problems and an elliptic problem in complex domain without interface, all with tens of millions of unknowns, are provided to show that the proposed multigrid solver is dozens of times faster than the classical algebraic multigrid method as implemented in the code AMG1R5 by Stüben.展开更多
This paper proposes a deep unfitted Nitsche method for solving elliptic interface problems with high contrasts in high dimensions.To capture discontinuities of the solution caused by interfaces,we reformulate the prob...This paper proposes a deep unfitted Nitsche method for solving elliptic interface problems with high contrasts in high dimensions.To capture discontinuities of the solution caused by interfaces,we reformulate the problem as an energy minimization problem involving two weakly coupled components.This enables us to train two deep neural networks to represent two components of the solution in highdimensional space.The curse of dimensionality is alleviated by using theMonte-Carlo method to discretize the unfittedNitsche energy functional.We present several numerical examples to show the performance of the proposed method.展开更多
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 lo...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.展开更多
With the remarkable empirical success of neural networks across diverse scientific disciplines,rigorous error and convergence analysis are also being developed and enriched.However,there has been little theoretical wo...With the remarkable empirical success of neural networks across diverse scientific disciplines,rigorous error and convergence analysis are also being developed and enriched.However,there has been little theoretical work focusing on neural networks in solving interface problems.In this paper,we perform a convergence analysis of physics-informed neural networks(PINNs)for solving second-order elliptic interface problems.Specifically,we consider PINNs with domain decomposition technologies and introduce gradient-enhanced strategies on the interfaces to deal with boundary and interface jump conditions.It is shown that the neural network sequence obtained by minimizing a Lipschitz regularized loss function converges to the unique solution to the interface problem in H2 as the number of samples increases.Numerical experiments are provided to demonstrate our theoretical analysis.展开更多
A new approach is given to analyse the regularity of solutions near singular points for the interface problems of second order elliptic partial differential equations. For general equations with nonsymmetric dominant ...A new approach is given to analyse the regularity of solutions near singular points for the interface problems of second order elliptic partial differential equations. For general equations with nonsymmetric dominant terms and discontinuous piecewise smooth coefficients, it is proved that solutions in H 1 can be docomposed into two parts, one of which is a finite sum of particular solutions to the corresponding homogeneous equations with piecewise constant coefficients, and the other one of which is the regular part. Moreover a priori estimations are proven.展开更多
In this paper, we study Nitsche extended finite element method （XFEM） for the inter- face problem of a two dimensional diffusion equation. Specifically, we study the quadratic XFEM scheme on some shape-regular famil...In this paper, we study Nitsche extended finite element method （XFEM） for the inter- face problem of a two dimensional diffusion equation. Specifically, we study the quadratic XFEM scheme on some shape-regular family of grids and prove the optimal convergence rate of the scheme with respect to the mesh size. Main efforts are devoted onto classifying the cases of intersection between the elements and the interface and prove a weighted trace inequality for the extended finite element functions needed, and the general framework of analysing XFEM c^n be implemented then.展开更多
The embedded boundary method for solving elliptic and parabolic problems in geometrically complex domains using Cartesian meshes by Johansen and Colella （1998, J. Comput. Phys. 147, 60） has been extended for ellipti...The embedded boundary method for solving elliptic and parabolic problems in geometrically complex domains using Cartesian meshes by Johansen and Colella （1998, J. Comput. Phys. 147, 60） has been extended for elliptic and parabolic problems with interior boundaries or interfaces of discontinuities of material properties or solutions. Second order accuracy is achieved in space and time for both stationary and moving interface problems. The method is conservative for elliptic and parabolic problems with fixed interfaces. Based on this method, a front tracking algorithm for the Stefan problem has been developed. The accuracy of the method is measured through comparison with exact solution to a two-dimensional Stefan problem. The algorithm has been used for the study of melting and solidification problems.展开更多
.In this paper,an augmented Lagrangian Uzawa iterative method is developed and analyzed for solving a class of double saddle-point systems with semidefinite(2,2)block.Convergence of the iterativemethod is proved under....In this paper,an augmented Lagrangian Uzawa iterative method is developed and analyzed for solving a class of double saddle-point systems with semidefinite(2,2)block.Convergence of the iterativemethod is proved under the assumption that the double saddle-point problem exists a unique solution.An application of the iterative method to the double saddle-point systems arising from the distributed Lagrange multiplier/fictitious domain(DLM/FD)finite element method for solving elliptic interface problems is also presented,in which the existence and uniqueness of the double saddle-point system is guaranteed by the analysis of the DLM/FD finite element method.Numerical experiments are conducted to validate the theoretical results and to study the performance of the proposed iterative method.展开更多
In this paper, we propose adaptive finite element methods with error control for solving elasticity problems with discontinuous coefficients. The meshes in the methods do not need to fit the interfaces. We establish a...In this paper, we propose adaptive finite element methods with error control for solving elasticity problems with discontinuous coefficients. The meshes in the methods do not need to fit the interfaces. We establish a residual-based a posteriori error estimate which is λ- independent multiplicative constants; the Lame constant λ steers the incompressibility. The error estimators are then implemented and tested with promising numerical results which will show the competitive behavior of the adaptive algorithm.展开更多
This article is to discuss the bilinear and linear immersed finite element(IFE)solutions generated from the algebraic multigrid solver for both stationary and moving interface problems.For the numerical methods based ...This article is to discuss the bilinear and linear immersed finite element(IFE)solutions generated from the algebraic multigrid solver for both stationary and moving interface problems.For the numerical methods based on finite difference formulation and a structured mesh independent of the interface,the stiffness matrix of the linear system is usually not symmetric positive-definite,which demands extra efforts to design efficient multigrid methods.On the other hand,the stiffness matrix arising from the IFE methods are naturally symmetric positive-definite.Hence the IFE-AMG algorithm is proposed to solve the linear systems of the bilinear and linear IFE methods for both stationary and moving interface problems.The numerical examples demonstrate the features of the proposed algorithms,including the optimal convergence in both L 2 and semi-H1 norms of the IFE-AMG solutions,the high efficiency with proper choice of the components and parameters of AMG,the influence of the tolerance and the smoother type of AMG on the convergence of the IFE solutions for the interface problems,and the relationship between the cost and the moving interface location.展开更多
The author studies the structure of solutions to the interface problems for second order linear elliptic partial differential equations in three space dimension. The set of singular points consists of some singular li...The author studies the structure of solutions to the interface problems for second order linear elliptic partial differential equations in three space dimension. The set of singular points consists of some singular lines and some isolated singular points. It is proved that near a singular line or a singular point, each weak solution can be decomposed into two parts, a singular part and a regular part. The singular parts are some finite sum of particular solutions to some simpler equations, and the regular parts are bounded in some norms, which are slightly weaker than that in the Sobolev space H^2.展开更多
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 fitt...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,...,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 O(t−k nτk+t−1 n h 2|log h|)is established for the fully discrete finite element method for general L 2 initial data.展开更多
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 meth...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.展开更多
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 explic...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.展开更多
Interface problems for second order quasi-linear elliptic partial differential equations in a two-dimensional space are studied.We prove that each weak solution can be decomposed into two parts near singular points,on...Interface problems for second order quasi-linear elliptic partial differential equations in a two-dimensional space are studied.We prove that each weak solution can be decomposed into two parts near singular points,one of which is a finite sum of functions of the form cr~α log^m r(?)(θ),where the coefficients c depend on the H^1-norm of the solution,the C^(0,δ)-norm of the solution,and the equation only;and the other one of which is a regular one,the norm of which is also estimated.展开更多
基金supported by the NSFC(11571177)the National Key Research and Development Program of China(2020YFA0713803).
文摘In this paper,we consider the plasma-vacuum interface problem in a cylindrical tube region impressed by a special background magnetic field.The interior region is occupied with plasma,which is governed by the incompressible inviscid and resistive MHD system without damping term.The exterior vacuum region is governed by the so-called the“pre-Maxwell equations”.And on the free interface,additionally,the effect of surface tension is taken into account.The original region can be transformed into a horizontally periodic slab through the cylindrical coordinate transformation,which will be impressed by a uniform nonhorizontal magnetic field.Appending with the appropriate physical boundary conditions,the global well-posedness of the problem is established by the energy method.
基金supported by the US ARO grants 49308-MA and 56349-MAthe US AFSOR grant FA9550-06-1-024+1 种基金he US NSF grant DMS-0911434the State Key Laboratory of Scientific and Engineering Computing of Chinese Academy of Sciences during a visit by Z.Li between July-August,2008.
文摘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.
基金Project supported by the National Natural Science Foundation of China(No.11271340)
文摘The lowest order Pl-nonconforming triangular finite element method （FEM） for elliptic and parabolic interface problems is investigated. Under some reasonable regularity assumptions on the exact solutions, the optimal order error estimates are obtained in the broken energy norm. Finally, some numerical results are provided to verify the theoretical analysis.
基金Supported by the National Natural Science Foundation of China (11071216 and 11101361)
文摘A second order isoparametric finite element method （IPFEM） is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the immersed curved interface is discontinuous. Based on an initial Cartesian mesh, a mesh optimization strategy is presented by employing curved boundary elements at the interface, and an incomplete quadratic finite element space is constructed on the optimized mesh. It turns out that the number of curved boundary elements is far less than that of the straight one, and the total degree of freedom is almost the same as the uniform Cartesian mesh. Numerical examples with simple and complicated geometrical interfaces demonstrate the efficiency of the proposed method.
基金Supported by National Natural Science Foundation of China(11571002,11461046)Natural Science Foundation of Jiangxi Province,China(20151BAB211013,20161ACB21005)+2 种基金Science and Technology Project of Jiangxi Provincial Department of Education,China(150172)Science Foundation of China Academy of Engineering Physics(2015B0101021)Defense Industrial Technology Development Program(B1520133015)
文摘In this paper,the minimal dissipation local discontinuous Galerkin method is studied to solve the elliptic interface problems in two-dimensional domains.The interface may be arbitrary smooth curves.It is shown that the error estimates in L;-norm for the solution and the flux are O（h;|log h|）and O(h|log h|;),respectively.In numerical experiments,the successive substitution iterative methods are used to solve the LDG schemes.Numerical results verify the efficiency and accuracy of the method.
基金supported by the National Natural Science Foundation of China(Grant No.42274101)X.X.Wu was supported by the Fundamental Research Funds for the Central Universities of Central South University(Grant No.2020zzts354)+2 种基金H.L.Hu was supported by the National Natural Science Foundation of China(Grant No.12071128)by the Natural Science Foundation of Hunan Province(Grant No.2021JJ30434)Z.L.Li was supported by a Simons Grant No.633724.
文摘The aim of this paper is to develop a fast multigrid solver for interpolation-free finite volume (FV) discretization of anisotropic elliptic interface problems on general bounded domains that can be described as a union of blocks. We assume that the curved interface falls exactly on the boundaries of blocks. The transfinite interpolation technique is applied to generate block-wise distorted quadrilateral meshes, which can resolve the interface with fine geometric details. By an extensive study of the harmonic average point method, an interpolation-free nine-point FV scheme is then derived on such multi-block grids for anisotropic elliptic interface problems with non-homogeneous jump conditions. Moreover, for the resulting linear algebraic systems from cell-centered FV discretization, a high-order prolongation operator based fast cascadic multigrid solver is developed and shown to be robust with respect to both the problem size and the jump of the diffusion coefficients. Various non-trivial examples including four interface problems and an elliptic problem in complex domain without interface, all with tens of millions of unknowns, are provided to show that the proposed multigrid solver is dozens of times faster than the classical algebraic multigrid method as implemented in the code AMG1R5 by Stüben.
基金supported by Andrew Sisson Fund of the University of MelbourneX.Y.was partially supported by the NSF grants DMS-1818592 and DMS-2109116.
文摘This paper proposes a deep unfitted Nitsche method for solving elliptic interface problems with high contrasts in high dimensions.To capture discontinuities of the solution caused by interfaces,we reformulate the problem as an energy minimization problem involving two weakly coupled components.This enables us to train two deep neural networks to represent two components of the solution in highdimensional space.The curse of dimensionality is alleviated by using theMonte-Carlo method to discretize the unfittedNitsche energy functional.We present several numerical examples to show the performance of the proposed method.
基金the supports by National Science and Technology Council,Taiwan,under the research grants 111-2115-M-008-009-MY3,111-2628-M-A49-008-MY4,111-2115-M-390-002,and 110-2115-M-A49-011-MY3,respectivelythe supports by National Center for Theoretical Sciences,Taiwan.
文摘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.
基金the National Natural Science Foundation of China(Grant Nos.11771435,22073110 and 12171466).
文摘With the remarkable empirical success of neural networks across diverse scientific disciplines,rigorous error and convergence analysis are also being developed and enriched.However,there has been little theoretical work focusing on neural networks in solving interface problems.In this paper,we perform a convergence analysis of physics-informed neural networks(PINNs)for solving second-order elliptic interface problems.Specifically,we consider PINNs with domain decomposition technologies and introduce gradient-enhanced strategies on the interfaces to deal with boundary and interface jump conditions.It is shown that the neural network sequence obtained by minimizing a Lipschitz regularized loss function converges to the unique solution to the interface problem in H2 as the number of samples increases.Numerical experiments are provided to demonstrate our theoretical analysis.
文摘A new approach is given to analyse the regularity of solutions near singular points for the interface problems of second order elliptic partial differential equations. For general equations with nonsymmetric dominant terms and discontinuous piecewise smooth coefficients, it is proved that solutions in H 1 can be docomposed into two parts, one of which is a finite sum of particular solutions to the corresponding homogeneous equations with piecewise constant coefficients, and the other one of which is the regular part. Moreover a priori estimations are proven.
基金The first author is partially supported by the U.S. Department of Energy, Office of Science, Office of Advanced Scientific Computing Research as part of the Collaboratory on Mathematics for Mesoscopic Modeling of Materials under Award Number DE-SC-0009249, and the Key Program of National Natural Science Foundation of China with Grant No. 91430215. The second author is supported by State Key Laboratory of Scientific and Engineering Computing （LSEC）, National Center for Mathematics and Interdisciplinary Sciences of Chinese Academy of Sciences （NCMIS）, and National Natural Science Foundation of China with Grant No. 11471026 he is thankful to the Center for Computational Mathematics and Applications, the Pennsylvania State University, where he worked on this manuscript as a visiting scholar. The authors are grateful to Professor Jinchao Xu, Dr. Yuanming Xiao and Dr. Maximilian Metti for their valuable suggestions and discussions, to Professor Haijun Wu for his valuable help on preparing the numerical example, and to the anonymous referee for the valuable comments and suggestion which lead to improvements of the paper.
文摘In this paper, we study Nitsche extended finite element method （XFEM） for the inter- face problem of a two dimensional diffusion equation. Specifically, we study the quadratic XFEM scheme on some shape-regular family of grids and prove the optimal convergence rate of the scheme with respect to the mesh size. Main efforts are devoted onto classifying the cases of intersection between the elements and the interface and prove a weighted trace inequality for the extended finite element functions needed, and the general framework of analysing XFEM c^n be implemented then.
基金supported by the U.S.Department of Energy under Contract No.DE-AC02-98CH10886 and by the State of New York
文摘The embedded boundary method for solving elliptic and parabolic problems in geometrically complex domains using Cartesian meshes by Johansen and Colella （1998, J. Comput. Phys. 147, 60） has been extended for elliptic and parabolic problems with interior boundaries or interfaces of discontinuities of material properties or solutions. Second order accuracy is achieved in space and time for both stationary and moving interface problems. The method is conservative for elliptic and parabolic problems with fixed interfaces. Based on this method, a front tracking algorithm for the Stefan problem has been developed. The accuracy of the method is measured through comparison with exact solution to a two-dimensional Stefan problem. The algorithm has been used for the study of melting and solidification problems.
基金supported by the 10 plus 10 project of Tongji University(No.4260141304/004/010).
文摘.In this paper,an augmented Lagrangian Uzawa iterative method is developed and analyzed for solving a class of double saddle-point systems with semidefinite(2,2)block.Convergence of the iterativemethod is proved under the assumption that the double saddle-point problem exists a unique solution.An application of the iterative method to the double saddle-point systems arising from the distributed Lagrange multiplier/fictitious domain(DLM/FD)finite element method for solving elliptic interface problems is also presented,in which the existence and uniqueness of the double saddle-point system is guaranteed by the analysis of the DLM/FD finite element method.Numerical experiments are conducted to validate the theoretical results and to study the performance of the proposed iterative method.
文摘In this paper, we propose adaptive finite element methods with error control for solving elasticity problems with discontinuous coefficients. The meshes in the methods do not need to fit the interfaces. We establish a residual-based a posteriori error estimate which is λ- independent multiplicative constants; the Lame constant λ steers the incompressibility. The error estimators are then implemented and tested with promising numerical results which will show the competitive behavior of the adaptive algorithm.
基金supported by DOE grant DE-FE0009843National Natural Science Foundation of China(11175052)GRF of HKSAR#501012 and NSERC(Canada).
文摘This article is to discuss the bilinear and linear immersed finite element(IFE)solutions generated from the algebraic multigrid solver for both stationary and moving interface problems.For the numerical methods based on finite difference formulation and a structured mesh independent of the interface,the stiffness matrix of the linear system is usually not symmetric positive-definite,which demands extra efforts to design efficient multigrid methods.On the other hand,the stiffness matrix arising from the IFE methods are naturally symmetric positive-definite.Hence the IFE-AMG algorithm is proposed to solve the linear systems of the bilinear and linear IFE methods for both stationary and moving interface problems.The numerical examples demonstrate the features of the proposed algorithms,including the optimal convergence in both L 2 and semi-H1 norms of the IFE-AMG solutions,the high efficiency with proper choice of the components and parameters of AMG,the influence of the tolerance and the smoother type of AMG on the convergence of the IFE solutions for the interface problems,and the relationship between the cost and the moving interface location.
文摘The author studies the structure of solutions to the interface problems for second order linear elliptic partial differential equations in three space dimension. The set of singular points consists of some singular lines and some isolated singular points. It is proved that near a singular line or a singular point, each weak solution can be decomposed into two parts, a singular part and a regular part. The singular parts are some finite sum of particular solutions to some simpler equations, and the regular parts are bounded in some norms, which are slightly weaker than that in the Sobolev space H^2.
文摘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,...,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 O(t−k nτk+t−1 n h 2|log h|)is established for the fully discrete finite element method for general L 2 initial data.
基金supported by the National Key R&D Program of China(Project No.2020YFA0712000)supported by the Shanghai Science and Technology Innovation Action Plan in Basic Research Area(Project No.22JC1401700)+1 种基金the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDA25010405)the National Natural Science Foundation of China(Grant No.DMS-11771290).
文摘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.
基金the National Natural Science Foundation of China(Grant No.DMS-12101553,Grant No.DMS-11771290)the Natural Science Foundation of Zhejiang Province(Grant No.LQ22A010017)+4 种基金the National Key Research and Development Program of China(Project No.2020YFA0712000)the Science Challenge Project of China(Grant No.TZ2016002)the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDA25000400)the National Science Foundation of America(Grant No.ECCS-1927432)also partially supported by the National Science Foundation of America(Grant No.DMS-1720420).
文摘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.
基金supported by the China State Major Key Project for Basic Researchesthe Science Fund of the Ministry of Education of China
文摘Interface problems for second order quasi-linear elliptic partial differential equations in a two-dimensional space are studied.We prove that each weak solution can be decomposed into two parts near singular points,one of which is a finite sum of functions of the form cr~α log^m r(?)(θ),where the coefficients c depend on the H^1-norm of the solution,the C^(0,δ)-norm of the solution,and the equation only;and the other one of which is a regular one,the norm of which is also estimated.
