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.展开更多
In this paper two dimensional elliptic interface problem with imperfect contact is considered,which is featured by the implicit jump condition imposed on the imperfect contact interface,and the jumping quantity of the...In this paper two dimensional elliptic interface problem with imperfect contact is considered,which is featured by the implicit jump condition imposed on the imperfect contact interface,and the jumping quantity of the unknown is related to the flux across the interface.A finite difference method is constructed for the 2D elliptic interface problems with straight and curve interface shapes.Then,the stability and convergence analysis are given for the constructed scheme.Further,in particular case,it is proved to be monotone.Numerical examples for elliptic interface problems with straight and curve interface shapes are tested to verify the performance of the scheme.The numerical results demonstrate that it obtains approximately second-order accuracy for elliptic interface equations with implicit jump condition.展开更多
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.展开更多
.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.展开更多
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,we introduce a nonconforming Nitsche’s extended finite element method(NXFEM)for elliptic interface problems on unfitted triangulation elements.The solution on each side of the interface is separately ex...In this paper,we introduce a nonconforming Nitsche’s extended finite element method(NXFEM)for elliptic interface problems on unfitted triangulation elements.The solution on each side of the interface is separately expanded in the standard nonconforming piecewise linear polynomials with the edge averages as degrees of freedom.The jump conditions on the interface and the discontinuities on the cut edges(the segment of edges cut by the interface)are weakly enforced by the Nitsche’s approach.In the method,the harmonic weighted fluxes are used and the extra stabilization terms on the interface edges and cut edges are added to guarantee the stability and the well conditioning.We prove that the convergence order of the errors in energy and L 2 norms are optimal.Moreover,the errors are independent of the position of the interface relative to the mesh and the ratio of the discontinuous coefficients.Furthermore,we prove that the condition number of the system matrix is independent of the interface position.Numerical examples are given to confirm the theoretical results.展开更多
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.展开更多
Deep learning techniques for solving elliptic interface problems have gained significant attentions.In this paper,we introduce a hybrid residual and weak form(HRW)loss aimed at mitigating the challenge of model traini...Deep learning techniques for solving elliptic interface problems have gained significant attentions.In this paper,we introduce a hybrid residual and weak form(HRW)loss aimed at mitigating the challenge of model training.HRW utilizes the functions residual loss and Ritz method in an adversary-system,which enhances the probability of jumping out of the local optimum even when the loss landscape comprises multiple soft constraints(regularization terms),thus improving model’s capability and robustness.For the problem with interface conditions,unlike existing methods that use the domain decomposition,we design a Pre-activated ResNet of ResNet(PRoR)network structure employing a single network to feed both coordinates and corresponding subdomain indicators,thus reduces the number of parameters.The effectiveness and improvements of the PRoR with HRW are verified on two-dimensional interface problems with regular or irregular interfaces.We then apply the PRoR with HRW to solve the size-modified Poisson-Boltzmann equation,an improved dielectric continuum model for predicting the electrostatic potentials in an ionic solvent by considering the steric effects.Our findings demonstrate that the PRoR with HRW accurately approximates solvation free-energies of three proteins with irregular interfaces,showing the competitive results compared to the ones obtained using the finite element method.展开更多
We present a parallel Cartesianmethod to solve elliptic problems with complex immersed interfaces.This method is based on a finite-difference scheme and is second-order accurate in the whole domain.The originality of ...We present a parallel Cartesianmethod to solve elliptic problems with complex immersed interfaces.This method is based on a finite-difference scheme and is second-order accurate in the whole domain.The originality of the method lies in the use of additional unknowns located on the interface,allowing to express straightforwardly the interface transmission conditions.We describe the method and the details of its parallelization performed with the PETSc library.Then we present numerical validations in two dimensions,assorted with comparisons to other related methods,and a numerical study of the parallelized method.展开更多
This paper introduces a new spontaneous potential log model for the case in which formation resistivity is not piecewise constant. The spontaneous potential satisfies an elliptic boundary value problem with jump condi...This paper introduces a new spontaneous potential log model for the case in which formation resistivity is not piecewise constant. The spontaneous potential satisfies an elliptic boundary value problem with jump conditions on the interfaces. It has beer/ shown that the elliptic interface problem has a unique weak solution. Furthermore, a jump condition capturing finite difference scheme is proposed and applied to solve such elliptic problems. Numerical results show validity and effectiveness of the proposed 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.展开更多
Based on an error estimate in terms of element edge vectors on arbitrary unstructured simplex meshes,we propose a new edge-based anisotropic mesh refinement algorithm.As the mesh adaptation indicator,the error estimat...Based on an error estimate in terms of element edge vectors on arbitrary unstructured simplex meshes,we propose a new edge-based anisotropic mesh refinement algorithm.As the mesh adaptation indicator,the error estimate involves only the gradient of error rather than higher order derivatives.The preferred refinement edge is chosen to reduce the maximal term in the error estimate.The algorithm is implemented in both two-and three-dimensional cases,and applied to the singular function interpolation and the elliptic interface problem.The numerical results demonstrate that the convergence order obtained by using the proposed anisotropic mesh refinement algorithm can be higher than that given by the isotropic one.展开更多
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 repre...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.展开更多
基金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(Grants 12261067,12161067,12361088,62201298,12001015,51961031)the Inner Mongolia Autonomous Region"Youth Science and Technology Talents"support program(Grant NJYT20B15)+1 种基金the Inner Mongolia Scientific Fund Project(Grants 2020MS06010,2021LHMS01006,2022MS01008)by the Innovation fund project of Inner Mongolia University of science and technology-Excellent Youth Science Fund Project(Grant 2019YQL02).
文摘In this paper two dimensional elliptic interface problem with imperfect contact is considered,which is featured by the implicit jump condition imposed on the imperfect contact interface,and the jumping quantity of the unknown is related to the flux across the interface.A finite difference method is constructed for the 2D elliptic interface problems with straight and curve interface shapes.Then,the stability and convergence analysis are given for the constructed scheme.Further,in particular case,it is proved to be monotone.Numerical examples for elliptic interface problems with straight and curve interface shapes are tested to verify the performance of the scheme.The numerical results demonstrate that it obtains approximately second-order accuracy for elliptic interface equations with implicit jump condition.
基金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.
基金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.
基金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.
基金The work of the second author was partially supported by the Natural Science Foundation of the Jiangsu Higher Institutions of China(No.18KJB110015)by No.GXL2018024+1 种基金The work of the third author was partially supported by the the NSF of China grant No.10971096by the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions.
文摘In this paper,we introduce a nonconforming Nitsche’s extended finite element method(NXFEM)for elliptic interface problems on unfitted triangulation elements.The solution on each side of the interface is separately expanded in the standard nonconforming piecewise linear polynomials with the edge averages as degrees of freedom.The jump conditions on the interface and the discontinuities on the cut edges(the segment of edges cut by the interface)are weakly enforced by the Nitsche’s approach.In the method,the harmonic weighted fluxes are used and the extra stabilization terms on the interface edges and cut edges are added to guarantee the stability and the well conditioning.We prove that the convergence order of the errors in energy and L 2 norms are optimal.Moreover,the errors are independent of the position of the interface relative to the mesh and the ratio of the discontinuous coefficients.Furthermore,we prove that the condition number of the system matrix is independent of the interface position.Numerical examples are given to confirm the theoretical results.
基金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.
文摘Deep learning techniques for solving elliptic interface problems have gained significant attentions.In this paper,we introduce a hybrid residual and weak form(HRW)loss aimed at mitigating the challenge of model training.HRW utilizes the functions residual loss and Ritz method in an adversary-system,which enhances the probability of jumping out of the local optimum even when the loss landscape comprises multiple soft constraints(regularization terms),thus improving model’s capability and robustness.For the problem with interface conditions,unlike existing methods that use the domain decomposition,we design a Pre-activated ResNet of ResNet(PRoR)network structure employing a single network to feed both coordinates and corresponding subdomain indicators,thus reduces the number of parameters.The effectiveness and improvements of the PRoR with HRW are verified on two-dimensional interface problems with regular or irregular interfaces.We then apply the PRoR with HRW to solve the size-modified Poisson-Boltzmann equation,an improved dielectric continuum model for predicting the electrostatic potentials in an ionic solvent by considering the steric effects.Our findings demonstrate that the PRoR with HRW accurately approximates solvation free-energies of three proteins with irregular interfaces,showing the competitive results compared to the ones obtained using the finite element method.
文摘We present a parallel Cartesianmethod to solve elliptic problems with complex immersed interfaces.This method is based on a finite-difference scheme and is second-order accurate in the whole domain.The originality of the method lies in the use of additional unknowns located on the interface,allowing to express straightforwardly the interface transmission conditions.We describe the method and the details of its parallelization performed with the PETSc library.Then we present numerical validations in two dimensions,assorted with comparisons to other related methods,and a numerical study of the parallelized method.
基金supported by the National Natural Science Foundation of China (No. 10431030)the Shanghai Natural Science Foundation (No. 08ZR1401100)
文摘This paper introduces a new spontaneous potential log model for the case in which formation resistivity is not piecewise constant. The spontaneous potential satisfies an elliptic boundary value problem with jump conditions on the interfaces. It has beer/ shown that the elliptic interface problem has a unique weak solution. Furthermore, a jump condition capturing finite difference scheme is proposed and applied to solve such elliptic problems. Numerical results show validity and effectiveness of the proposed 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 National Basic Research Program under the Grant 2005CB321701the National Science Foundation of China under the grant 10771008 and 10771211partial supported by A Foundation for the Author of National Excellent Doctoral Dissertation of PRC.
文摘Based on an error estimate in terms of element edge vectors on arbitrary unstructured simplex meshes,we propose a new edge-based anisotropic mesh refinement algorithm.As the mesh adaptation indicator,the error estimate involves only the gradient of error rather than higher order derivatives.The preferred refinement edge is chosen to reduce the maximal term in the error estimate.The algorithm is implemented in both two-and three-dimensional cases,and applied to the singular function interpolation and the elliptic interface problem.The numerical results demonstrate that the convergence order obtained by using the proposed anisotropic mesh refinement algorithm can be higher than that given by the isotropic one.
文摘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.